Abstract
Background:
Brain tumors are a common cause of morbidity, disability, cognitive deterioration and mortality in children, even after treatment. Little is know about the specific causes. The study aimed to assess potential socio-demographic and antenatal factors in primary brain tumor (PBTs) in children and young people (CYP) in Karachi, Pakistan.
Designs and methods:
A single center hospital based matched case control study in Karachi, Pakistan. Cases were defined as CYP aged between 5 and 21 years with any histological type and grade of primary brain tumor of any histology, stage or grade. Data were collected from parents of 244 patients at the selected center between 2017 and 2021 via telephonic interview. Controls were 5–21 years old CYP admitted with non-oncological diagnoses matched on age and sex. Matched Odds Ratios for predictors of brain tumor in children were derived. Those of statistical significance were included in a multivariable logistic regression model.
Results:
In the adjusted model, lower paternal education (matched adjusted odds ratio (maOR) 2.46; 95% CI 1.09–5.55), higher household monthly income (maOR 3.4; 95% CI 1.1–10.2), antenatal paternal use of addictive substances (maOR 19.5; 95% CI 2.1–179.8), and antenatal maternal use of analgesics during pregnancy (maOR 3.0; 95% CI 1.2–7.9) were all independently predictive of brain tumors.
Conclusion:
This matched case-control study found novel associations between maternal use of analgesics, paternal use of addictive substances, higher household income, and lower paternal education and Primary Brain Tumors in Children and Young People. Longitudinal multicenter studies will be required to test these associations prospectively.
Keywords
Introduction
Central nervous system (CNS) tumors account for approximately 20% of all childhood malignancies 1 and the second most common group of cancers in this age group. In Pakistan, data from the Karachi Cancer Registry shows the age-standardized incidence rate of CNS tumors among children of 3.14 per 100,000 while in adolescent 0.58 per 100,000. 2 Primary brain tumors (PBTs) in children and young people (CYP) age 0–24 years old 1 are a diverse group of neoplasms with unique histopathology, molecular features, and etiology. Tumor histology, location, age at diagnosis, sex, race, and ethnicity all correlate prognosis. 1
Little is known regarding the pathogenesis of most PBTs in CYP. However, certain genetic syndromes, such as Neurofibromatosis type 1 and 2, Fanconi anemia and tuberous sclerosis have been associated with increased risk to develop brain tumors. 3 Moreover, family history of brain tumor has also found to be associated with PBTs in CYP. 4 Other reported or suspected nongenetic risk factors include, younger age of the child, 5 male gender, 6 higher parental age at the time of conception, 3 non-ionizing radiation, exposure to allergens, high birth weight of the child, exposure to infections early in life, greater parental age, high socioeconomic status, and exposure to pesticides.1,4,7–9 N-nitroso compounds, parental smoking, and maternal antenatal use of hair dye are hypothesized to be associated with PBTs as well, though conclusive evidence is lacking. 10
Most of the literature comes from high-income countries and data from lower-middle-income countries (LMICs) like Pakistan is scarce. The poor prognosis compels the continued search for potentially reversible factors that might be addressed to prevent disease. This study aimed to identify socio-demographic and antenatal factors by using a matched control study design.
Methodology
Study design and setting
A matched case-control study was conducted. A matched cases and controls based on specific characteristics allows for better control of potential confounding variables. By matching cases and controls on age and gender that are the potential confounders we tried to minimize the potential influence of these variables on the outcome leading to a more accurate assessment of the exposure-disease relationship. Individual matching was done that is, for each case a control was selected who closely resembled the case in terms of the matching criteria (age and gender). The patients were recruited from the Aga Khan University Hospital (AKUH), which is a Joint Commission International Accreditation (JCIA-accredited) private tertiary care hospital in Karachi. AKU is the largest children′s oncology referral center in Pakistan. Children are referred from the city, province and whole country.
Study participants
Cases were defined as children and young people (CYP) aged between 5 and 21 years 1 with PBTs of any histology, stage or grade presenting to a tertiary care hospital of Karachi, Pakistan between 2017 and 2021. Controls were 5 and 21 years old CYP admitted with non-oncological diagnoses matched on age and sex presenting to AKUH contemporaneously. Only CYP with parents who spoke and understood either Urdu or English were considered. CYP whose parents refused to consent were excluded.
Sampling strategy and data collection
Purposive sampling technique was employed to select participants. We identified cases diagnosed between 1st January 2017 and 31st December 2021 from Health Information Management Services (HIMS) in AKUH. Data were collected from the parents of the cases and controls via telephonic interview using a predesigned questionnaire developed by the investigators based on previously reported antenatal factors of brain tumor (Figure 1). The questionnaire comprised of child related factors such as; demographics (age, gender, province of residence, mother tongue), family history of brain tumor or other cancer, gestational age and birth order. The questionnaire also included information on parental factors such as; demographics (age, educational status, and household income) and antenatal factors such as; parental age at child’s birth, parental use of addictives, parental smoking, maternal use of medications, maternal consumption of cured meat, maternal exposure to pesticides, ultrasound & X-ray, maternal use of hair spray or hair color and maternal comorbidities. 4

Flow of the study.
Sample size
The sample size was calculated on PASS (Power Analysis and Sample Size) 11 software licensed by NCSS. In previous studies, the probability of exposure (higher parental age, parental education, socioeconomic status and antenatal factors) among sampled control patients ranged between 16% and 70%4,7,10,11 with 26% of the controls having a high socioeconomic status. 8 A minimum sample of 122 CYP with PBTs (cases) was calculated. For each CYP with PBT, a matching sample of 1 control CYP, matched on age and sex was also obtained. This sample of 244 patients achieved an 80% power to detect an odds ratio of 2.5 versus the alternative of equal odds using a chi-square test with a 0.05 significance level.4,7,10,11
Ethical considerations
Ethical/institutional review committee AKU-ERC approval was obtained (ERC# 2020-4859-11855). Verbal consent was taken from the parents of the cases and controls after explaining the study procedure and its potential risks and benefits to them in Urdu or English as appropriate. All study materials containing personal identifiers were kept in a locked file cabinet. A unique study identification number was assigned to each participant. Data were entered in a password-protected electronic database that was only accessible to the research team.
Statistical analysis
Data were analyzed using STATA version 12. Results were presented as mean and standard deviation/ median (IQR) for normally and non-normally numeric variables. Comparison of numeric variables were made using independent sample t-tests/Mann-Whitney U-Test, as appropriate. Categorical variables were reported as frequency and percentages and compared using Chi-squared/Fisher exact tests as appropriate. Correlations to identify collinearity were performed using Pearson’s correlation. Unadjusted and adjusted matched odds ratio (maOR) with their 95% CI were reported by using conditional logistic regression to determine the association of socio-economic and antenatal factors with brain tumor in CYP. A p-value < 0.2 was used as the screening cut-off on univariate analysis. A p < 0.05 for multivariable analysis was considered statistically significant throughout the study.
Results
A total of 122 cases and 122 controls were recruited. Description of cases and controls are shown in Table 1.
Diagnosis of children and young people with and without Brain tumor.
Socio-demographic factors
A higher proportion of cases versus controls were from the province of Punjab. Both maternal and paternal years of education were higher among parents of controls than cases. Fathers of controls were more likely to have received education of secondary and above as compared to the cases (Table 2).
Socio-demographic characteristics of the study participants.
Others for the variable mother tongue includes: Brushia (1), Burushaski (3), Chitrali (2), Haryanvi (1), Hazargi (1), Hindko (2), Kachi (1), Kashmiri (1), Kostan (1), Marvari (1), Memoni (7), Satri (1), and Shina (7).
Exchange rate USD to PKR 202; Euro to PKR 215.
Significant at p value <0.05 by using Chi square or Fisher exact test.
Significant at p value <0.05 by using Mann Whitney U test.
Antenatal factors
The paternal mean age at child’s birth was significantly lower in cases than controls. Paternal use of addictive substances was higher in cases (9%) than controls (0.8%). A higher proportion of mothers of cases used pain relievers and anti-emetics during pregnancy (16.4% vs 7.4% and 23% vs 11.5%) respectively as compared to the controls. These findings are shown in Table 3.
Antenatal characteristics of study participants.
Significant at p value <0.05 by using Chi square or Fisher exact test.
Significant at p value <0.05 by using independent t test or Mann Whitney U test.
Logistic regression
Multivariable logistic regression was performed, matched for age and sex, while adjusting for covariates identified via univariate screening (Table 4). The factors that were associated with the outcome included: high family income, paternal use of addictive substances, maternal use of pain killers and lower paternal education. It was observed that the odds of a higher monthly income of >Rs 80,000 (>$397.6)/(>Є373.57) was 3.4 times significantly higher in cases as compared to controls. Moreover, the odds of paternal use of addictive was 19.5 times significantly higher in cases as compared to controls. Maternal use of pain killers was also 3.7 times significantly higher in cases versus controls. However, the odds of fathers having no education or primary education was 2.46 times significantly higher in cases as compared to controls.
Univariate and multivariable to determine independent socio-demographic and antenatal risk factors of PBTs in children and young people (CYP).
p < 0.2 on univariate regression.
p < 0.05 on multivariable regression.
Discussion
Our results found maternal use of analgesics, paternal use of addictive substances, higher household income, and lower paternal education, to be statistically significant, independent factors of PBTs in CYP. Our study also found that analgesics use during pregnancy were reported by more than 60% women, whether by prescription or self-medication. This is a new finding: previous literature has shown associations with perinatal opioid use, but not other non-opioid analgesia and pain-killers in general are not associated with the development of pediatric cancers.8,12–14
Our study also found that paternal use of addictive substances to be associated with a 20-times increased risk of offspring developing PBTs. Though a 2014 meta-analysis found no relationship between paternal smoking before or during pregnancy, 12 the strength of the relationship in our study warrants a closer look. A meta-analysis in 2016 concluded that maternal smoking during pregnancy increases the risk of PBTs in offspring, 13 suggesting that the exposure to carcinogenic substances in cigarette smoke may be a risk factor. However, in our study, there were too few maternal smokers to explore this. LMICs, like Pakistan, have significantly higher rates of indoor. smokers, which increases second-hand exposure to cigarette smoke amongst household members. 14 In Pakistan, more than 40% of pregnant women are exposed to second-hand smoke in quantities sufficient to cause adverse health outcomes in the fetus and offspring. 15 Thus, it is possible that maternal exposure to second-hand smoke, due to paternal use of addictive substance such as smoking during pregnancy, may account for the statistically significant association with PBTs in CYP observed in our study.
We also identified higher household income to be statistically significant, independent risk factor for the development of PBTs in CYP. The trend of better economic status being associated with development of PBTs has been observed in different settings globally.8,9 While the reasons for this association are yet poorly understood and likely to vary across different countries and settings, a few possible links have been identified. These include exposure to causal risk factors, behavioral differences, underlying racial differences, and even greater case ascertainment amongst higher economic groups due to better access to healthcare. 9 Given that Pakistan is a LMIC where access to cancer care is severely restricted amongst the lower socioeconomic strata of the population, 16 the latter explanation likely contributes to our findings. This case ascertainment bias is further amplified in Pakistan, where the lack of insurance or government coverage means that the majority of healthcare costs are borne out-of-pocket by patients. 17 Our study may additionally have been prone to a selection bias, since both cases and controls were recruited from one of the best equipped private hospitals in the country, where pediatric neuro-oncologic care can be afforded out-of-pocket only by patients from more privileged economic backgrounds. The distribution of provinces of residency across our study participants corroborates this point, as more than 40% of CYP with PBTs had traveled from different provinces across Pakistan to seek care at AKUH, indicating the financial means to seek healthcare at one of the country’s top hospitals.
Lastly, our findings also demonstrated a higher risk of PBTs in CYP whose fathers were less educated. While data from high-income countries (HICs) reveals either no 18 or the opposite association, 10 the educational landscape in Pakistan differs significantly from HICs. Just over half the population in Pakistan attends school, with an average of less than 6 years of formal schooling. 19 Similar to economic status, educational background likely affects disease risk through a variety of factors, such as exposures, behavioral practices, ethnicity and race, disease-related awareness, and health-seeking behavior. These factors remain to be studied in the local context of Pakistan, and are an avenue inviting future investigation.
Limitations
Our study has several limitations. Firstly, being a single-center study, its data may not be generalizable to other centers across Pakistan. Secondly, as data were collected using a questionnaire, responses were prone to recall bias. Recall bias was mitigated by assessing exposure information from both cases and controls in a consistent manner. Thirdly, we anticipate reporting bias regarding maternal addictive use because addictive use and substance abuse are often stigmatized in many societies, including Pakistan. Women may face additional social and cultural stigma associated with substance use, leading to underreporting due to fear of judgment, shame, or negative consequences. Thirdly, while our study identified a three times higher risk of PBT with analgesic use, we were not able to explore the relationship for opioid and non-opioid analgesics separately and did not have sufficiently granular data to test a dose-response effect to strengthen causal inference. Nevertheless, our findings lay the ground for future detailed investigation into the risk conferred by different types of analgesics as well as the mechanisms involved in increasing risk.
Strengths
Our study has several strengths. Firstly, it is a matched case control study and matching on age and gender allowed the comparability of cases and controls in terms of these variables. This reduced the potential for bias due to these factors and strengthens the ability to attribute observed differences to the exposure of interest. Secondly, this is the first of its kind of study from our part of the world.
Conclusion
This matched case-control study, the first of its kind from Pakistan, found an association between maternal use of analgesics, paternal use of addictive substances, higher household income, and lower paternal education and primary brain tumors in children and young people which were robust to known potential confounders. Longitudinal multicenter studies will be required to test these associations prospectively
Clinical implications
Identifying these risk factors can help healthcare providers recognize children at higher risk, leading to early detection and timely treatment initiation. Moreover, knowledge of specific risk factors would enable healthcare providers to counsel expectant parents about the potential risks and empower them to make informed decisions regarding prenatal care. Understanding these risk factors can also contribute to development of strategies to minimize or avoid them during pregnancy, potentially reducing the incidence of brain tumors. Findings from the study can lead to the formulation of specific recommendations in prenatal care guidelines for brain tumor, enabling tailored care and early management of risk factors.
Patient and public involvement
This was a matched case control study design and the parents of cases and controls were interviewed via telephone by trained research assistants. The study findings will be disseminated to different stakeholders, such as healthcare professionals and cancer patients through: publications at local, national and international journals, presentations at conferences and workshops and through research briefs.
Footnotes
Acknowledgements
The authors would like to acknowledge Ms Aneesa Hassan for facilitating in data collection. The authors thank all families of the patients for their generous participation.
Author contributions
NZ conceived the study, analyzed data, and critically reviewed the manuscript. SAE, TM and AM overlooked the study and critically reviewed the manuscript. FU and RSM revised, and critically reviewed the manuscript. FK collected the data. NM, MM, NM, SA and SK critically reviewed the manuscript. NB overlooked the study and critically reviewed the manuscript.
Availability of data
Data will be available upon request from the corresponding author.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The study was funded by faculty development award (FDA), Aga Khan University Karachi, Pakistan. The funders have no role in the preparation of data or the manuscript.
Ethical approval
The research was conducted ethically in accordance with the World Medical Association Declaration of Helsinki. Ethical approval was obtained from the Aga Khan University ethics review committee. Ethical/institutional review committee AKU-ERC approval was obtained (ERC# 2020-4859-11855). Verbal consent was taken from the parents of the cases and controls after explaining the study procedure and its potential risks and benefits to them in Urdu or English as appropriate. All study materials containing personal identifiers were kept in a locked file cabinet. A unique study identification number was assigned to each participant. Data were entered in a password-protected electronic database that was only accessible to the research team.
