Time trends in social contacts before and during the COVID-19 pandemic: the CONNECT study

Study design

CONNECT (CONtact and Network Estimation to Control Transmission) is a population-based survey of epidemiologically relevant social contacts and mixing patterns conducted in the province of Quebec, Canada. The first phase of CONNECT was conducted in 2018-2019 (February 2018 to March 2019), one year before the COVID-19 pandemic. Four additional phases of CONNECT were undertaken to document changes in social contacts during the COVID-19 pandemic period (CONNECT2: April 21^st-May 25^th 2020 and CONNECT3,4,5: July 3^rd 2020-February 26^th 2021) (Additional file 1: Table S1). All CONNECT phases were conducted with the same methodology.

Recruitment of participants

The target population of CONNECT consisted of all non-institutionalized Quebecers without any age restriction (e.g., elderly living in retirement homes who generally have personal phone lines were eligible but those living in long-term care homes (nursing homes, Quebec CHSLD) were not eligible). We used random digit dialling to recruit participants. The randomly generated landline and mobile phone number sample was provided by ASDE, a Canadian firm specialized in survey sampling [4]. After having explained the study, verified eligibility of the household and documented the age and sex of all household members, we randomly selected one person per household to participate in the study, using a probability sample stratified by age. This recruitment procedure was sequentially repeated for every new phase of CONNECT (i.e, new participants were recruited for every CONNECT phase).

Data collection

We collected data using a self-administered web-based questionnaire. A secured individualized web link to the questionnaire and information about the study were sent by email to each selected participant who consented to participate in the study. Parents of children aged less than 12 years were asked to complete the questionnaire on behalf of their child, whereas teenagers aged 12-17 years completed their own questionnaire, after parental consent.

The same questionnaire was used for all CONNECT phases. The first section of the questionnaire documented key socio-demographic characteristics. The second section was a social contact diary, based on instruments previously used in Polymod and other similar studies[5-7] (an example of the diary is provided in the Additional file 1: Figure S1). Briefly, participants were assigned two random days of the week (one week day and one weekend day) to record every different person they had contact with between 5 am and 5 am the following morning. A contact was defined as either physical (handshake, hug, kiss) or nonphysical (two-way conversation in the physical presence of the person, at a distance equal or less than 2 meters, irrespective of masking). Participants provided the characteristics of the contact persons (age, sex, ethnicity, and relationship to themselves (e.g., household member, friend, colleague)) as well as characteristics of the contacts with this person: location where the contact(s) occurred (home, work, daycare/school, public transport, leisure, other location), duration, usual frequency of contact with that person, and whether the contact was physical or not. Participants reporting more than 20 professional contacts per day were asked not to report all their professional contacts in the diary. Instead they were asked general questions about these professional contacts: age groups of the majority of contact persons, average durations of contacts and whether physical contacts were generally involved or not. Additional questions about teleworking were included from CONNECT2 onwards.

All CONNECT phases were approved by the ethics committee of the CHU de Québec research center (project 2016-2172) and we commissioned the market company Advanis for recruitment and data collection. All participants gave their consent to participate in the study during the recruitment phone call. Informed consent was taken from a parent and/or legal guardian for study participation in the case of minors.

Analyses

We weighted the participants of the CONNECT 1-5 surveys by age, sex, region (Greater Montreal and other Quebec regions), and household composition (households without 0-17-year-olds, households with 0-5-year-olds, if not with 6-14-year-olds, if not with 15-17-year-olds), using the Quebec data of the 2016 Canadian census (Additional file 1: Table S2) and we verified that they were representative of the Quebec population for key socio-demographic characteristics. To obtain daily number of social contacts on a weekly basis, we weighted the number of daily contacts reported during the week (5/7) and the weekend (2/7). We classified the type of employment of workers using the 2016 National occupation classification (NOC)[8].

We estimated the number of social contacts per person and per day, for all locations combined and for 6 different locations: home, work, school, public transportation, leisure, and other locations. To do so, several steps were necessary. First, for a contact person met in more than a single location during a single day, the location of the contact was assigned in the following hierarchical order, according to risk of transmission: home, work, school, public transport, leisure and other locations [9]. For example, if a parent reported contacts with his child at home, in public transportation and in a leisure activity, we only considered the home contact to avoid counting contacts with the same person multiple times. Second, for workers reporting more than 20 professional contacts per working day, we added their reported number of professional contacts to the work location for their working day(s). Similar to other studies which allowed a maximal number of contacts per day [5, 6, 10], we truncated professional contacts at a maximum of 40 per day to eliminate extreme values and contacts at low risk of transmission of infectious diseases. Third, we identified all workers in schools through their NOC code and job descriptions and attributed their professional contacts to the school location. We did so to describe social contacts in schools, not only between students, but also between students and their teachers, educators, and other school’s workers. Unless specified, we estimated the mean number of contacts in the different locations using a population-based denominator. With this method, all individuals were considered in the denominator of each location had they reported contacts or not for that location. The sum of contacts in the different locations gives the total number of contacts.

Using data available from CONNECT1-5, we determined different periods to reflect the Quebec COVID-19 epidemiology, their related physical distancing measures, and expected seasonality in social contacts (Additional file 1: Figure S2). We used data collected from February 1^st 2018 to March 17^th 2019 as our pre-COVID period. We used data collected from April 21^st to May 25^th 2020 to represent the first wave, data collected from July 3^rd to August 31^st 2020 to represent the summer, and data collected from September 1^st 2020 to February 26^th 2021 to represent the second wave. We further stratified the second wave to represent periods of expected seasonality in social contacts: September with the return to school and at work, fall with gathering restrictions, the Christmas holidays with school and work vacations and closure of non-essential business, January and February 2021 with the gradual return to work and school after Christmas vacations and school/non-essential business closures, and the introduction of a curfew. We used a Canadian stringency index, adapted from the Oxford COVID-19 Government Response Tracker (OxCGRT) [11], to quantify the intensity of public health measures in Quebec over time [12]. This index is obtained by averaging the intensity score of 12 policy indicators (e.g., school closures, workplace closures, gathering restrictions, stay-at-home requirements, etc.) and higher values indicate stricter measures. We estimated the mean stringency index for each of the 8 periods described previously by averaging the daily values of the index.

We compared the mean number of social contacts over time (total or by location) using weighted generalized linear models with a Poisson distribution and an identity link. Generalized estimating equations with robust variance [13] were used to account for the correlation between the two days of diary data collection and overdispersion. A categorical period effect was included in the model and is presented as the absolute difference in the mean number of contacts compared to the previous period. We also compared the mean number of social contacts according to different key socio-demographic characteristics using the same model with a period-by-covariate interaction and adjusting for age (in 8 categories). In this model, period and characteristic effects were tested using contrasts: each period was compared to the previous period within each level of the covariate, and the global effect of the characteristic was tested within each period.. We also examined the association between the mean number of social contacts and the stringency index (in 5 categories), irrespective of periods, using a model similar to the one comparing periods.

Finally, we estimated mixing matrices. The entries of the mixing matrix represent the mean number of social contacts per person per day according to the age of the respondent (column) and the age of his contacts (row). Mixing matrices were estimated separately for the 8 periods described previously and for 3 categories of contact locations: all locations, home (contacts with household members and visitors), any location outside home. The matrices were obtained by maximizing a constrained log likelihood of the number of reported contacts per day among CONNECT participants weighted by age, sex, household composition and region. The number of contacts was assumed to follow a negative binomial distribution. The likelihood constraint ensured that the total number of contacts between individuals of age i and age j is the same whether it is estimated from entry (i,j) or entry (j,i) of the total mixing matrix including contacts in all locations (i.e., reciprocity of the mixing matrix).

All statistical analyses were performed with SAS version 9.4. Maximization of the log likelihood for the mixing matrices was performed using a nonlinear programming algorithm (nlminb2 function from the ROI package in R).

Participants

A total of 1291, 546, and 4970 Quebecers completed the social contact questionnaires during the pre-COVID period (CONNECT1), the first wave (CONNECT2), and summer 2020 and second wave (CONNECT3-5), respectively. Participation rates (number of questionnaires completed among consenting participant) were 30%, 38%, and 34% for CONNECT 1,2, and 3-5 respectively (Additional file 1: Figure S3). These participants were generally representative of the Quebec general population, and they were comparable across the different phases of CONNECT (Table 1).

Time trends in the number of social contacts

During the pre-pandemic period, the mean number of social contacts per person per day was 7.8 (95% confidence interval (CI):7.2-8.5) (Figure 1 and Additional file 1: Table S3). This number decreased significantly by 60% during the spring 2020 lockdown to 3.1 (95% CI:2.6-3.5). It then increased gradually during summer 2020 and peaked at 5.0 (95% CI:4.3-5.8) contacts/day in September 2020; this peak coincided with the return to school and work. The mean number of contacts decreased significantly again during fall 2020 to 4.1 (95% CI:3.7-4.5) when physical distancing measures were intensified in Quebec to control the second wave. The mean number of social contacts also decreased significantly during the Christmas holidays at 2.9 (95% CI:2.7-3.1) because of school and work vacations and closure of non-essential businesses. There was a trend towards increasing numbers of contacts in January (3.5, 95% CI (3.0-3.9)) with the gradual return to school and in February 2021 (4.0, 95% CI (3.3-4.6)) with the re-opening of non-essential businesses. These time trends in social contacts closely followed the intensity of public health measures as quantified by the stringency index (Figure 1). The mean number of contacts was also significantly associated with the stringency index, irrespective of periods (Additional file 1: Table S4 and Figure S4).

• Download figure
• Open in new tab

Figure 1. Time trends in the number of social contacts in the province of Quebec.

Pre-COVID: February 1^st 2018 to March 17^th 2019; Spring 2020: April 21^st to May 25^th 2020; Summer 2020: July 3^rd to August 31^st 2020; Fall 2020: October 1^st to December 16^th 2020; Holidays 2020-2021: December 17^th 2020 to January 8^th 2021.

1^st wave: Spring 2020; 2^nd wave: September 2020 to February 2021.

* Contacts at work were truncated to a maximum of 40 contacts per day.

^† Contacts for workers in schools were included in the school location.

Stringency index: Higher values indicate stricter measures; the mean stringency index for each of the 8 periods was obtained by averaging the daily values of the index.

Error bars represent the 95% confidence interval of the total number of social contacts

During the pre-pandemic period, the great majority of contacts occurred at home (2.3 contacts: 1.2 with household members and 1.1 with visitors), at work (2.7 contacts) and at school (1.6 contacts) (Figure 1 and Additional file 1: Table S3). The mean number of contacts at home with household members remained relatively constant over time (1.2 to 1.4 contacts), whereas the number of contacts at home with visitors varied significantly through the study period with lower numbers observed during the spring 2020 lockdown, in January and February 2021 (0.2-0.3 contacts). Compared to the pre-pandemic period, contacts at work and school decreased significantly during the spring lockdown (1.2 and 0.0, respectively), summer (1.2 and 0.2) and the holidays (0.8 and 0.0) and peaked in September 2020 (1.5 and 0.9) with the return at work and school. Contacts in the other locations (transport, leisure and other locations) represented a small proportion of overall contacts during the pre-pandemic period (1.3 contacts). They also decreased significantly since the beginning of the pandemic and stayed low through the study period.

Time trends in the number of social contacts by age

The location of social contacts varied substantially by age (Figure 2, Additional file 1: Table S5). Contacts in households represented an important part of contacts for all age groups. Contacts in other locations were highly dependent on age. The main locations of contacts away from home for individuals aged 0-17, 18-65, and >65 years were, respectively, school, work, and other locations.

• Download figure
• Open in new tab

Figure 2. Time trends in the number of social contacts in the province of Quebec, according to age of participants and location of contacts.

A) 0-17 years old *

B) 18-65 years old ^†

C) C) > 65 years old

Pre-COVID: February 1^st 2018 to March 17^th 2019; Spring 2020: April 21^st to May 25^th 2020; Summer 2020: July 3^rd to August 31^st 2020; Fall 2020: October 1^st to December 16^th 2020; Holidays 2020-2021: December 17^th 2020 to January 8^th 2021.

1^st wave: Spring 2020; 2^nd wave: September 2020 to February 2021.

* Only the main locations of contacts were included in this figure for 0-17-year-olds (contacts at work were excluded)

^† For adults: contacts at work are truncated to a maximum of 40 contacts per day and they include contacts in schools for school workers and students.

Contacts in other locations include transports, leisure and other locations.

Error bars represent the 95% confidence interval of the total number of social contacts

During the pre-pandemic period, the mean number of social contacts at school/daycare for youth aged 0-17 years was 3.3 contacts (Figure 2A). These contacts significantly decreased to nearly 0 during the spring 2020 lockdown and the Christmas holidays. They reached the pre-pandemic level with the return to school in September (4.0 contacts), during fall 2020 (3.0) and in February 2021 (3.7). Except for post-secondary, similar time trends in contacts at school/daycare were observed by education level (daycare, elementary, high school) (Additional file 1: Table S6).

During the pre-pandemic period, the mean number of contacts at work for adults aged 18-65 years was 5.6 (Figure 2B). These contacts significantly decreased to 1.7 during the spring 2020 lockdown and thereafter remained significantly lower than the pre-pandemic period (from 1.3 during the Christmas holidays to 2.7 in September). The number of contacts at work varied by the type of occupation and the proportion of workers reporting teleworking, and therefore having no contact at work (Additional file 1: Tables S7,S8). During the pre-pandemic period, the greatest number of contacts at work were reported by workers in the domains of Sales & services (10.8), Management (10.2), and Health (10.1). Contacts at work decreased during the spring 2020 lockdown for the majority of domains and remained lower than the pre-pandemic period thereafter. Except for workers in the domains of Health and Sales & Services, the majority of workers in the other domains (> 50%) reported teleworking since the beginning of the pandemic.

During the pre-pandemic period, the mean number of social contacts in other locations for adults older than 65 years was 1.6 (Figure 2C). These contacts decreased significantly at the beginning of the pandemic and remained low through the study period (between 0.2 and 0.8). Therefore, adults older than 65 years had virtually no contact outside their house during this period.

Time trends in the number of social contacts by key socio-demographic characteristics

During the pre-pandemic period, the mean number of social contacts was significantly higher among individuals living in households with ≥3 individuals (vs households with 1-2 individuals), in households with 0-17-year-olds (vs household without 0-17-year-olds), among native French or English speakers (vs other mother tongues), and among individuals with a university degree (vs no degree) (Figure 3, Additional file 1: Table S5). During the first wave, social contacts significantly decreased for most socio-demographic characteristics. The mean number of social contacts slightly increased after the first wave for all socio-demographic characteristics, although it remained lower than the pre-pandemic period through the study period. During the second wave, the only significant differences between socio-demographic characteristics were a higher number of contacts in households with more individuals and/or households with 0-17-year-olds, mainly explained by the greater number of contacts with household members. Of note, individuals with a university degree had the greatest decrease of their social contacts during the first wave (from 10.7 to 2.5, p<0.0001) and their contacts remained relatively low through the study period (2.5 to 4.4).

• Download figure
• Open in new tab

Figure 3. Time trends in the total number of social contacts in the province of Quebec, according to key socio-demographic characteristics (adjusted for age).

A) Household size

B) Household composition

C) Level of education (among ≥ 25 years old)

D) Mother tongue

Pre-COVID: February 1^st 2018 to March 17^th 2019; 1^st wave: April 21^st to May 25^th 2020; Summer 2020: July 3^rd to August 31^st 2020; 2^nd wave: Sept 1^st 2020 to February 26^th 2021.

Error bars represent the 95% confidence interval of the total number of social contacts

Time trends in social contact matrices

During the pre-pandemic period, the mixing matrices indicated a high assortativity of contacts by age (as illustrated by the central diagonal), and mixing between children and adults, mainly at home (as illustrated by the 2 secondary diagonals) (Figure 4). These general mixing patterns with 3 diagonals remained apparent during the different pandemic periods, even though the number of contacts was substantially reduced. Interestingly, the matrices of contacts in any location outside home clearly illustrate the impact of school closures or holidays (Spring 2020, Summer 2020, Holidays 2020-21) with fewer contacts between children/teenagers. The matrices of contacts at home (with household members and visitors) also illustrate the impact of restrictions on private gatherings (Spring 2020, Fall 2020 to February 2021) with fewer contacts outside the 3 main diagonals and contacts limited to household members (i.e., same-age partners/siblings and children/parents).

• Download figure
• Open in new tab

Figure 4. Time trends in social contacts matrices.

Pre-COVID: February 1^st 2018 to March 17^th 2019; Spring 2020: April 21^st to May 25^th 2020; Summer 2020: July 3^rd to August 31^st 2020; Fall 2020: October 1^st to December 16^th 2020; Holidays 2020-2021: December 17^th 2020 to January 8^th 2021.

1^st wave: Spring 2020; 2^nd wave: September 2020 to February 2021.

The matrices of contacts at home include contacts with household members and visitors. The matrices of contacts in other locations include contacts at work, school, transport, leisure, and other locations.