Study design
This study aims to evaluate the relationship between rotavirus vaccination and the risk of intussusception via the design of TTE, which incorporates eligibility criteria, treatment strategies, treatment assignment, outcomes, start and end of follow-up, causal contrasts, and a data-analysis plan(20). The detailed TTE design is shown in Supplementary Table 1.
Eligibility criteria: This study selects infants aged 6 to 12 weeks who were born locally in the study area and were vaccinated with BCG and the first dose of hepatitis B vaccine. The exclusion criteria includes a history of previous rotavirus vaccination, intussusception, other serious diseases, such as malignant tumors, immunodeficiency or autoimmune diseases, and medication within 14 days (such as high fever, infection, and the use of immune preparations, etc.).
Treatment strategy: The intervention group of the study is infants vaccinated with rotavirus vaccine. Exposure is defined as at least one dose of rotavirus vaccine (RotaTeq or LLR) in the baseline period. The non-exposed group is infants who have not received any rotavirus vaccine at baseline or during the follow-up period.
Treatment assignment: The groups are divided according to the vaccination status at baseline. Furthermore, the randomization is simulated by adjusting the baseline confounding factors to ensure the comparability between the groups.
Outcome definition: The main outcome event of the study is the first confirmed case of intussusception in children aged 0-3 years, and the diagnosis is based on medical records such as outpatient, emergency, and hospitalization records, as well as intussusception confirmed by imaging (such as X-ray, B-scan ultrasonography, and CT).
Start and end of follow-up: Follow-up starts from the baseline period and ends at the first diagnosis of intussusception, death, loss to follow-up, the expiration of the 90-day follow-up period, or the end of the study, whichever occurs first.
Causal contrasts: This study adopts intention-to-treat analysis and per-protocol analysis. The intention-to-treat analysis is grouped by baseline period. The per-protocol analysis addresses protocol violation cases. Protocol violation cases are defined as incomplete vaccination or rotavirus vaccination in the control group.
Statistical analysis: Pooled logistic regression models will be used to estimate risk ratios and generate survival curves, with follow-up time (months) set as linear or quadratic terms. In addition, subgroup analysis (such as analysis of vaccine type, gender, and number of doses) and sensitivity analysis will be performed.
The study plans to include a birth cohort from 2015 to 2023. Given the 99% coverage rate of the national immunization program vaccines for children of appropriate age by 2020, it is assumed that the immunization planning information system can cover all locally born resident children. "Locally born" is defined as being born in a local hospital or having a household registration address in the region (province/city). "Resident" is defined as having a permanent or residential address in the region (province/city). The rotavirus vaccination schedules are as follows: LLR is recommended for infants aged 2 months to 3 years, with one dose annually. Rotateq is recommended at 2, 4, and 6 months of age, with the first dose given at 6-12 weeks, subsequent doses at 4-10 week intervals, and all three doses completed by 32 weeks of age. Considering the above immunization schedules, a 90-day postvaccination risk period, and potential delays in diagnosis, the study observation period for the rotavirus vaccine is set from birth to 39 months of age, and that for Rotateq is set from birth to 11 months of age. The follow-up will utilize the vaccination schedules for children aged 0-3 years (including both immunization program vaccines and nonimmunization program vaccines), migration, and death data.
The cohort start time is defined as the baseline period (6-12 weeks of age for Rotateq, 8-52 weeks for LLR), with the first dose administration time for the exposed group. The cohort end time is defined as the earliest occurrence of any of the following: 1) first occurrence of intussusception, 2) death, 3) loss to follow-up, 4) 42 days after the last dose, or 5) end of the study. Loss to follow-up is defined as any of the following events during the observation period: death, migration, or failure to receive any vaccine doses (including both immunization program vaccines and nonimmunization program vaccines) at 8-12 months, 12-24 months, or 25-36 months of age (for RotaTeq recipients, loss to follow-up includes failure to receive any vaccine doses at 8-12 months of age).
During the observation period, the birth cohort will be divided into different subgroups on the basis of rotavirus vaccination status. The exposed group includes children who received at least one dose of rotavirus vaccine and were subdivided into rotavirus vaccine and Rotateq groups. The control group includes children who did not receive any rotavirus vaccine during the baseline and follow-up periods.
Data sources
The study will be conducted on regional health information platforms in Ningbo, Chongqing, and Jiangsu, which have been validated for data quality and sample size sufficiency. Ningbo, a subprovincial city in Zhejiang Province located on the southeast coast of China, had a permanent population of 9.618 million at the end of 2022, with an urbanization rate of 78.9%, a birth rate of 5.58‰, and a death rate of 6.83‰. Chongqing, a directly administered municipality and national central city located in Southwest China, had a permanent population of 32.133 million at the end of 2022, with an urbanization rate of 70.96%, a birth rate of 5.98‰, and a death rate of 8.09‰. Jiangsu Province, located on the east coast of China, had a permanent population of 85.15 million at the end of 2022, with an urbanization rate of 75.0%, a birth rate of 5.23‰, and a death rate of 7.04‰. These four regions have stable populations in the millions, which are suitable for establishing dynamic cohorts for the entire lifecycle and conducting active safety monitoring of vaccines for rare clinical outcomes at single or multiple centers. The study uses information from these platforms, including resident records, immunization planning, maternal and child health, electronic medical records, and mortality surveillance, as detailed in the Supplementary Table 2.
Data collection and management
Exposure measurement
The rotavirus vaccination records will be extracted from the immunization planning information, including anonymous vaccination ID, vaccine name or code, vaccination date, dose, manufacturer, and records of other vaccinations administered on the same day and within 14 days before and after each rotavirus vaccination. If the interval between two consecutive doses of rotavirus vaccine is less than the risk period, i.e., the risk period overlaps with the date of the next dose, the 0th day of the next dose is the end of the risk period for the previous dose. Fully vaccinated population is defined as the children with completion of three doses according to the recommended immunization schedule.
Outcome Measurement
The suspected cases of intussusception in children aged 0-3 years will be extracted from electronic medical records, including diagnoses from outpatient, emergency, inpatient, and discharge records with the term "intussusception" or the ICD-10 code K56.1; imaging examinations (X-ray, ultrasound, CT) indicating intussusception; Enema treatments (hydrostatic enema, air enema, ultrasound-guided enema) indicating intussusception; and surgical or autopsy indications of intussusception.
New-onset intussusception is defined as the interval between the current and the most recent episode of intussusception being greater than 72 hours. The onset time is defined as the earliest occurrence of the visit time, or specific clinical symptoms (paroxysmal abdominal pain, bloody stool, or an abdominal mass). Only the first episode of intussusception is included in this study.
In line with the Brighton criteria, the study will establish levels of certainty and criteria for diagnosing intussusception. An adjudication committee, including surgeons, pediatricians, and radiologists related to intussusception, will be formed. Surgeons will serve as adjudicators, with other physicians as auxiliary members. The adjudicators will independently review the data of suspected intussusception cases and assign certainty ratings. If there is a disagreement between two adjudicators, a more senior surgeon will be invited to reach a consensus through a meeting. The data collection and case adjudication process will be blinded to the patients' vaccination history.
Covariate measurement
General demographic information: Anonymous ID (including parents), anonymous vaccination ID, birth date, gender, ethnicity, registration date, migration date, birth hospital, household registration address (city), residential address (city), urban/rural area from the immunization planning information.
Mortality information: Anonymous ID, death date, death diagnosis, and ICD-10 code from mortality registration information.
Previous clinical information: Anonymous ID, diagnosis time, and disease history from 4-12 weeks of age (including high fever, acute respiratory infections, acute gastrointestinal infections, adenovirus infections, and diarrhea); medication history from 4-12 weeks of age (including immunomodulators, antiviral drugs, and antibiotics); malignancies, HIV infection, complement system deficiencies from 0-12 weeks; and any hospitalization records from 0 weeks.
Previous birth information: Anonymous ID, mother's last menstrual period, preterm birth, birth weight, Apgar score from maternal and child health information.
Unmeasured confounding factors: Including breastfeeding, other special health conditions such as congenital birth defects, other serious diseases (including metabolic diseases, neurological diseases, kidney diseases, liver diseases, and hematologic diseases). The negative outcome control is set as hand foot and mouth disease.
Statistical analysis
The statistical analysis will be conducted via SAS 9.4 or R software. Continuous variables will be presented as the means and standard deviations or median and quartiles. Group comparisons will be performed via ANOVA or the Wilcoxon rank-sum test, depending on the distribution of the data. Categorical variables will be presented as frequencies and percentages, and comparisons between groups will be carried out via the chi-square test.
For the primary study, a nested case-control design will be implemented. The analysis will include both an intention-to-treat (ITT) approach and a per-protocol analysis. In the ITT analysis, the incidence rate ratio (IRR) will be estimated via Poisson regression, adjusting for potential confounders such as age, gender, birth year, birth season, and other relevant factors. The per-protocol analysis will account for deviations from the study protocol by applying inverse probability weighting.
Subgroup analyses will further explore the exposure periods postvaccination, considering different time intervals (1-7 days, 8-14 days, 15-28 days, 29-42 days, and 43-90 days). Additionally, sensitivity analyses will be conducted to assess the robustness of the results. These will include day 0 of vaccination analysis, multiple imputations for missing data, self-controlled case series analysis, propensity score matching for comparison groups, and the use of a negative outcome control.
Sample size calculation
The conditional logistic regression module of PASS15 software is used to calculate the sample size required for the case-control study, with case-control ratio = 1:100, exposure rate Pᴇ = 0.3-0.5, and other parameters as above. The logarithm N of the required case-control matching is shown in the following table.
Table 1. Sample size for the nested case-control study
PE
|
OR
|
|
1.1
|
1.3
|
1.5
|
2.0
|
0.3
|
6182
|
816
|
342
|
117
|
0.4
|
5409
|
714
|
299
|
103
|
0.5
|
5193
|
686
|
287
|
99
|
Quality Control and Management
Before the study implementation, senior epidemiology, statistics, and clinical experts will be invited to guide the improvement of the study methodology. On the basis the final protocol, data management, case validation, and statistical analysis plans will be formulated.
During the study implementation, prejob training will be provided to all relevant personnel, including data managers, case adjudicators, and statistical analysts, who will participate in this study upon passing the tests. The study will be conducted in strict accordance with the relevant protocols, with parallel quality control and documentation conducted by the researchers.
At the conclusion of the study, senior epidemiology, statistics, and clinical experts will be invited to review the study report and related materials, evaluate the study's processes, and thoroughly review the study's conclusions to reach a consensus.