Study population:
Residents of Chiang Mai, Thailand, aged 18 years or older, presenting to any community-based testing facilities for a SARS-CoV-2 test between 1 Oct 2021 – 31 Dec 2021 (Delta-predominant) and 1 Feb 2022 – 10 Apr 2022 (Omicron- predominant) were assessed for eligibility to be included in the study. Molecular testing revealed 96.5% delta and 95.6% omicron lineage during 1 Oct 2021 – 31 Dec 2021 and 1 Feb 2022 – 10 Apr 2022 periods respectively. Tests done in Jan 2022 were excluded due to mixed delta-omicron lineage among samples (Omicron 75%, Delta 25%). The data capture ended on 10 April 2022 which was the last date community testing facilities were operational in Chiang Mai, before largely shifting to self-testing.
Subjects were included if they were either close-contacts of COVID-19 cases or attended an event where a COVID-19 outbreak was detected, irrespective of occurrence of symptoms. Those without history of exposure to COVID-19 were excluded. Non-Thai residents (foreigners and migrants) were excluded as the vaccination and other data for this group may be incomplete.
The patient selection flow is presented under Figures 1a and 1b.
Data Sources:
We have previously published the details on creating and implementing the different information systems used in this study.19 Brief descriptions of the data sources are outlined below.
Surveillance System (Epid-CM platform) and COVID-19 case management system in Chiang Mai (CMC-19 platform):
COVID-19 cases were reported in the surveillance system of Chiang Mai Provincial Health Office. According to Communicable Disease Control Act (B.E. 2558), all COVID-19 cases must be reported from every public and private health facilities in the province. Epid-CM, the web-based reporting system of COVID-19 cases in Chiang Mai was created by Chiang Mai Provincial Health Office and Faculty of Public Health of Chiang Mai University since April 2021. After detection of COVID-19 case, the patient details, clinical manifestations, history of exposure with COVID-19 cases or suspected outbreak, laboratory testing method and result were entered into the system under a unique ID for each patient. Epid-CM is synchronized with Chiang Mai hospital management platform for COVID-19 (CMC-19) for bed management in the overall province. Classification of severity was evaluated and entered into CMC-19 by the attending clinical team. Asymptomatic or mild cases were treated as out-patient and isolated in designated place in community, community-isolation, or home-isolation since the beginning of Omicron variant spread throughout Chiang Mai in Jan 2022. Data on progression of the disease and treatments are recorded in each hospital information system. Death cases are reported to Chiang Mai Provincial Health Office and also recorded in Epid-CM.
Community-based testing sites were initiated in the city of Chiang Mai, where most of the cases were located since April 2021, and provided free of charge COVID-19 tests. Those tested included close contacts of COVID-19 cases, those who attended an event with a COVID-19 outbreak, or those who have developed suggestive respiratory symptoms, i.e., fever, cough, rhinorrhea, anosmia, dyspnea. Health personnel perform either RT-PCR testing (mainly November 2021 to January 2022) or antigen testing (February 2022 onwards), and note down test details and results for all cases in Epid-CM.
Ministry of Public Health Immunization Center (MOPH IC) database:
All national vaccination records are deposited in the “Morprom” application, which was launched by the Ministry of Public Health, Thailand. This allows the residents to access vaccination services, including vaccine reservation and tracking as well as updated information on COVID-19. It also includes a feature that enables post-vaccination monitoring to check for side effects.
Ethical considerations:
The study was conducted on routine data collected as part of the national COVID-19 response under the Communicable Disease ACT (B.E. 2558) and was exempted from ethics review. Data were de-identified at source and analysed by Chiang Mai Provincial Health Office and Faculty of Public Health, Chiang Mai University.
Study Design:
A test-negative, case-control analysis was conducted to estimate VE against SARS-CoV-2 infection. Separate analyses were done for delta-predominant and omicron predominant periods. Cases were defined as those with a positive SARS-CoV-2 result, and controls were those with negative SARS-CoV-2 result, either by RT-PCR or medically administered antigen testing.
The type of COVID-19 vaccine, and date of vaccination were extracted from MOPH-IC. For the primary analysis, subjects who received their COVID-19 vaccination within 14 days of the test date were excluded to allow time for adequate immune responses to develop. Additionally, we explored the data incorporating subjects vaccinated less than 14 days but after 7 days prior to inclusion.
Statistical Analysis:
Descriptive statistics are reported separately for the cases and controls, stratified by delta and omicron predominance. Continuous variables are summarized as mean and standard deviation (SD) or median and interquartile range (IQR) depending on the distribution. Categorical variables are summarized as frequency and percentages. Between group comparisons were done using Mann-Whitney-U test or Kruskal Wallis test for continuous variables and Chi-squared test for categorical variables.
Associations between SARS-CoV-2 infection and heterologous vaccination schedule were estimated by comparing the odds of vaccination (exposed) vs no vaccination (unexposed) separately for delta and omicron predominant periods. The OR was used to estimate VE, where VE = (1 – OR) × 100% with 95% Confidence Interval (CI). Separate VEs were calculated for different vaccination schedules and stratified by age group. Forest plots were used to visualise the VEs.
Logistic regression models were used to adjusted for age, gender, calendar day of test (in weekly units), to calculate adjusted OR for SARS-CoV-2 infection separately for delta and omicron predominant periods. Primary vaccination series type, type of last vaccination received, and time since last vaccine dose were examined as separate factors in the model. Interaction terms were explored as appropriate. VE was estimated using adjusted ORs.
In addition to the test-negative case control study described above, we analysed data from the COVID-19 patient management program, which is a hospital patient management system, to examine the association between mortality and vaccination schedules.
All statistical analyses were conducted using stata (version 15.0 SE, College station, TX:StataCorp LP). Significance tests were 2 sided and a p-values <0.05 were considered statistically significant.