3.1. Study Setting
This case-control study was conducted at Ali-Asghar Children’s Hospital, a major pediatric referral center in Tehran, Iran, affiliated with Iran University of Medical Sciences. The hospital serves a diverse patient population from across Iran and internationally. The study spanned a 10-year period from 2013 to 2023, encompassing both case and control groups.
Cases comprised all children under 18 years of age diagnosed with leukemia or solid tumors (retinoblastoma, neuroblastoma, hepatoblastoma, Wilms tumor, medulloblastoma, germ cell tumors, or thyroid carcinoma) during the study period at Ali-Asghar Children’s Hospital. The control group consisted of children formally admitted as inpatients to the same hospital for non-oncologic reasons within the same timeframe. This specifically included inpatient children admitted to different departments. Patients who were treated and discharged directly from the emergency ward were not eligible as controls.
The use of hospital-based controls was chosen for practical reasons and because they share a similar referral pattern and geographic catchment area as the cases, making them a comparable population for this setting. Controls were selected from a broad range of inpatient departments, including general pediatrics, infectious diseases, surgery, orthopedics, and otorhinolaryngology, to create a heterogeneous group representative of the general pediatric inpatient population without cancer. To ensure the control group was representative of a "healthy" population in terms of the exposures under investigation, we excluded children admitted for conditions that might share risk factors with childhood cancers. These exclusions included:
-Chronic diseases with potential environmental or genetic links (e.g., congenital heart disease, genetic syndromes, autoimmune disorders).
-Conditions directly related to parental occupational exposures (e.g., specific poisonings).
-Premature infants admitted to the NICU for prolonged stays.
For each case, one control was selected. Controls were randomly chosen using computer-generated random numbers from the hospital’s admission lists of eligible non-oncologic patients. Eligible controls were then approached and recruited for participation in the study.
3.2. Control Selection Process
Controls were selected from children formally admitted as inpatients to the same hospital for non-oncologic reasons during the study period. To ensure methodological rigor and comparable healthcare-seeking behavior, our sampling frame was explicitly restricted to the hospital's inpatient admission lists; children treated solely in the emergency department or outpatient clinics were not eligible. From this inpatient population, we randomly selected one control for each case using computer-generated random numbers. To create a representative sample of the general pediatric population without cancer and to minimize confounding, controls were drawn from a broad range of inpatient departments (e.g., general pediatrics, infectious diseases, surgery) and were excluded if they had chronic diseases, genetic syndromes, or conditions with potential environmental or genetic links to cancer.
3.3. Inclusion Criteria
Participants were included in the study if they met the following criteria:
-Age less than 18 years
-For cases: Confirmed diagnosis of leukemia or solid tumor by histopathology, morphology, and flow cytometry
-For controls: Admission for an acute, non-chronic condition (e.g., acute appendicitis, minor trauma, non-chronic infectious diseases like pneumonia)
Power consideration: The sample size was determined by the total number of eligible cases presenting during the study period. While we captured all available solid tumor cases (n = 30) over the decade, this sample size inherently limits the statistical power for this subgroup. Therefore, the analysis for solid tumors should be interpreted as exploratory and descriptive.
Power analysis: A post-hoc power analysis was conducted using G*Power 3.1 software to quantify the statistical power for the solid tumor subgroup. Given the available 30 solid tumor cases and 114 controls, with an alpha of 0.05 and assuming a small-to-medium effect size (odds ratio = 1.8 - 2.0), the study achieved approximately 35 - 45% power for detecting associations with solid tumors. For the leukemia group (84 cases, 114 controls), under the same conditions, the study achieved 70 - 80% power. These calculations confirm that while the study was adequately powered for leukemia analyses, the solid tumor subgroup was underpowered to detect anything but very large effect sizes, supporting the exploratory interpretation of these results
3.4. Exclusion Criteria
These Participants were excluded:
For cases: The patient had died prior to recruitment. This decision was made for ethical considerations, to avoid causing distress to bereaved families, and due to concerns about the reliability of retrospectively collected exposure data from grieving parents.
For controls: A diagnosis of any chronic illness, genetic syndrome, or condition with a known or suspected environmental etiology that could plausibly be linked to the exposure variables under study (as listed in section 2.1).
3.5. Data Collection
Data were collected by trained healthcare staff under the direct supervision of the principal investigator. A structured, custom-designed checklist was used to systematically gather information during face-to-face or telephone interviews with one parent or guardian of each participant. Interviews were conducted in a quiet, controlled hospital environment to minimize distractions and ensure the quality of the data.
The Data Entry/Data Extraction Form captured the following variables:
-Demographic factors: Age, gender, birth weight, parental age, and parental education level
-Environmental and lifestyle factors: Parental occupational exposures (e.g., chemicals, electromagnetic fields, ionizing radiation), maternal obstetric history (including delivery method), infant feeding practices, parental smoking, and proximity of residence to industrial facilities.
All participants’ parents or guardians received detailed explanations regarding the study’s objectives, protocols, and potential implications before enrollment. The principal investigator provided ongoing oversight to ensure strict adherence to study protocols and to maintain the integrity of data collection.
3.6. Exposure Assessment
Data on environmental and lifestyle exposures were collected via parental report using a structured questionnaire. It is critical to note that variables such as 'house painting during pregnancy', 'proximity to a power pole,' and 'history of severe psychological stress' were operationalized as binary (yes/no) indicators based on parental perception and recall. To improve consistency in reporting, specific thresholds were used as prompts during the interviews: Exposure to house painting was defined as lasting for more than 10 days; proximity to a power pole was defined as living within an estimated 500-meter radius; and severe stress was explicitly defined as stemming from a major traumatic event such as bereavement or marital separation.
3.7. Statistical Analysis
Categorical and continuous variables were described as counts (N) and percentages (%), as well as means with standard deviations (SD). The chi-square test or Fisher's exact test was employed to assess the relationship between categorical variables. The Kolmogorov-Smirnov test was utilized to check the normality assumption for continuous variables. An ANOVA test was conducted to compare continuous variables across groups.
The dependent variable for this analysis was polytomous, comprising three mutually exclusive categories: Control, leukemia, and solid tumor. To examine the relationship between predictor variables and cancer type, we employed multinomial logistic regression, conducting both univariable and multivariable analyses. This technique was selected in preference to performing separate binary logistic regressions, as it allows for the simultaneous estimation of all relative risk ratios using a common referent category (the control group) within a single, coherent statistical framework. A key advantage of this approach is its superior statistical efficiency. Furthermore, it enables direct comparisons of coefficient estimates, thereby permitting an explicit evaluation of how specific risk factors differentially associate with the odds of a leukemia diagnosis relative to a solid tumor diagnosis.
Prior to finalizing the multivariable model, multi-collinearity between the independent variables was assessed using the variance inflation factor (VIF). A VIF value of 10 is often taken as an indicator of severe multi-collinearity, while a value above 5 may suggest moderate correlation. In our final model, all VIF values were below 3, indicating that multi-collinearity was not a significant concern and that the parameter estimates are stable. The VIF values were calculated for all variables in the final model; the maximum VIF was 1.87, indicating no evidence of multicollinearity.
Variables showing a P-value < 0.25 in univariable analysis were included in the initial multivariable multinomial logistic regression model. This liberal threshold is recommended to avoid excluding variables that may be important contributors in the multivariable context (
21). In addition to this statistical screening, variables considered to be theoretically important confounders based on existing literature (e.g., maternal age, child’s sex, and birth weight) were assessed for inclusion in the multivariable model regardless of their statistical significance in univariable analysis (
22).
The assumptions for the multinomial logistic regression were satisfied. Statistical significance was established for adjusted odds ratios (AORs) with 95% confidence intervals (CIs) using a two-tailed P-value < 0.05. Data analysis was conducted with the Statistical Package for Social Sciences (IBM SPSS Statistics for Windows, Version 23.0. Armonk, NY: IBM Corp). Given the limited sample size for the solid tumor subgroup, a post-hoc power analysis was performed using G*Power 3.1. This analysis indicated that with 30 solid tumor cases and 114 controls, the study had approximately 35 - 45% power to detect a small-to-medium effect size (AORs = 1.8 - 2.0) at α = 0.05. This confirms the exploratory nature of the analysis for solid tumors and underscores that the study was only powered to detect very large effect sizes for this outcome.
3.8. Ethical Consideration
This study was approved by the Research Ethical Review of Iran University of Medical Sciences (
IR.IUMS.REC.1399.504). The researcher obtained written informed consent from parents or guardians and assured them of confidentiality. All methods were performed following the relevant guidelines and regulations.