Abstract
Objective
The TP73 G4C14-A4T14 variant has been associated with elevated cancer risk, but the evidence is inconclusive. We performed a meta-analysis to clarify the role of this variant in cancer development.
Methods
Eligible literature was selected by searching PubMed, Google Scholar, Cochrane Library, and Embase. The meta-analysis was performed using Review Manager 5.4.
Results
A meta-analysis of 55 case–control studies showed that the G4C14-A4T14 variant was significantly associated with overall cancer development in five genetic models, including the allele model (AM), codominant model 1 (COD1), COD2, dominant model (DM), and over-dominant model (OD). Sub-group analysis based on ethnicity showed significantly higher risks in Africans in COD2 and RM and in Whites in AM, COD2, DM, and recessive model (RM). Cancer-specific subgroup analysis identified significant risks of gynecological (ovarian, cervical, and endometrial cancer), colorectal, oral, head and neck, and other cancers. Moreover, hospital-based controls revealed significant cancer risks in the AM, COD1, COD2, DM, and RM genetic models. Our findings were confirmed by trial sequential analysis.
Conclusion
This meta-analysis confirmed that TP73 G4C14-A4T14 significantly elevates the overall cancer risk, especially in White, African, and hospital-based populations, and specifically predisposes individuals to gynecological, colorectal, oral, and head and neck cancers.
This meta-analysis was registered at INPLASY (registration number: INPLASY202210070).
Introduction
Cancer is an evolving health problem and a major cause of death worldwide, with 19.3 million new cancer cases in 2020, and 10 million deaths due to various cancers. 1 Malignancies involve the accumulation of multiple genetic mutations, and scientists have discovered more than 10,000 genetic risk variants associated with susceptibility to cancer development. Mutations in tumor suppressor genes such as the TP53 family, especially loss of function mutations that suppress the actions of the genes, are among the most important factors associated with carcinogenesis. TP53 is the most widely investigated and common tumor suppressor gene, and has been found to be associated with almost all types of cancers. Researchers are now focusing on rare genetic variants to provide more specific information on cancer genetics.2,3
TP73 is a vital gene that encodes p73, an essential member of the p53 family that is structurally and functionally homologous to p53 (63% homologous amino acid sequence). This protein, also known as p53-like transcription factor, is involved in cellular proliferation, programmed cell death (apoptosis), cell cycle regulation or arrest, and transactivation of overlapping target genes such as the p21 gene.4–8 However, unlike TP53, mutations in TP73 are rare. During DNA damage, p73 is overexpressed in malignancies resulting from p53 mutation. It mimics the tumor suppression function of p53 by initiating the transcription of genes involved in cell cycle regulation, which are usually responsive to p53, repairing damaged DNA, promoting apoptosis, and preventing uncontrolled cellular growth and proliferation via blocking the G1 cell cycle checkpoint.9–14 p73 thus helps to maintain cellular homeostasis through compensating for the TP53 loss of function polymorphism.7,14,15 Although mutations in TP73 have been detected in less than 2% of all cancers, the gene is highly polymorphic and loss of heterozygosity polymorphisms have been reported in different types of tumors. TP73 is located at chromosomal region 1p36-33, which is deleted in many human cancers. This suggests that p73 might be strongly related to cancer susceptibility.16–18
Nineteen exonic and intronic single nucleotide polymorphisms (SNPs) have been identified in TP73, but none of these result in miscoded amino acids.19,20 Two common SNPs, rs2273953 and rs1801173, are located at positions 4 (G>A) and 14 (C>T), respectively, within a noncoding 5ʹ-untranslated region upstream of the TP73 promoter in exon 2. The distance between the two polymorphisms is short, with a tendency for non-random associations between them. The two polymorphisms are in complete disequilibrium with each other and are jointly referred to as G4C14-A4T14. This set of polymorphisms is located just above the translation initiation site and has been shown to affect TP73 gene expression levels by forming a stem-loop-like structure.19,21–23
Given its ability to modify the tumor suppression activity of TP73, the association between G4C14-A4T14 and carcinogenesis has recently been investigated in genome-wide association studies in multiple cancer types, including lung, colorectal, breast, cervical, gastric, esophageal, endometrial, oral, and ovarian cancer, in addition to head and neck squamous cell carcinoma, lymphoma, and cutaneous melanoma.2,24–74 However, the findings of these studies were inconsistent. Although previous meta-analyses have summarized the evidence regarding the roles of the G4C14-A4T14 polymorphism in different cancers, the numbers of studies included in those meta-analyses were limited,75–78 While a larger sample size provides firmer evidence in population-based genetic association studies.
In this study, we performed a comprehensive meta-analysis of 55 case–control studies to resolve previous controversies and provide systematic evidence for the association between the TP73 G4C14-A4T14 polymorphism and cancer development.
Materials and methods
This meta-analysis was performed following the updated PRISMA 2020 guidelines (available at https://www.bmj.com/content/372/bmj.n160). The need for obtaining informed consent from patients or controls was not applicable as no participants were directly involved in this study.
Literature search strategy
We carried out a comprehensive literature search of the PubMed, Google Scholar, Embase, Cochrane Library, Web of Science, and China National Knowledge Infrastructure electronic databases up to 20 July 2021, using the following key terms: ‘TP73 or p73’, ‘Cancer or tumor’, ‘G4C14-A4T14 polymorphism’, ‘rs2273953 and rs1801173’, ‘TP73 polymorphism and cancer’, and ‘association between TP73 G4C14-A4T14 polymorphism and cancer’. Additional studies were extracted from the reference lists of the selected literature. We also screened the ‘similar studies’ options in the above databases. Finally, published studies were included avoiding any language barriers.
Publication selection and eligibility criteria
The overall selection process was completed according to the authors’ predesigned protocol. Eligible studies containing the required data were selected and the data were organized for further analysis by comprehensive screening. The overall study selection method is outlined in a PRISMA flow diagram (Figure 1). Two authors (SJ and MAA) carefully revised the whole procedure, and the other author (MSI) conducted a final screening to reduce the chances of disagreement. This meta-analysis was retrospectively registered at INPLASY (https://inplasy.com/, registration number: INPLASY202210070).

Systematic flow diagram of study selection process.
The inclusion criteria of the selected studies were case–control studies examining the association between TP73 G4C14-A4T14 polymorphism and cancer susceptibility, studies with detailed comparative genotypic information for both controls and patients, and study population in agreement with the Hardy–Weinberg equilibrium (HWE) after adjustments. If the selected studies contained genotypic data on other SNPs, as well as the selected SNP, we only extracted data on the selected SNP for inclusion in this meta-analysis. We excluded studies without G4C14-A4T14 genotypic data for cancer patients and controls, studies lacking a control population data or with incomplete genotypic information, systematic reviews and meta-analyses, and studies conducted on cell lines or animal models.
Data extraction and quality assessment
We extracted the following information from the selected studies: study ID, date of publication, country, ethnicity or region of the recruited population, type of cancer, category of control population, type of genotyping method used, sample and control sizes, and genotypic data for the selected SNP. In addition, the HWE p-value was collected and adjusted (corrected) by Benjamini and Hochberg’s (1995) false discovery rate, 79 and the Newcastle–Ottawa Scale (NOS) score 80 was calculated from each selected study by the authors to maintain the quality of the selected studies. Two authors (SJ and MAA) extracted the above data from each study, and the other author (MSI) carried out a final screening of the organized data to avoid mistakes and misinterpretation.
Statistical analysis
The overall statistical analysis was carried out using Review Manager (RevMan) software version 5.4 (Cochrane Collaboration, 2020) to elucidate the impact of the TP73 G4C14-A4T14 variant on susceptibility to different cancers. We applied seven genetic association models to evaluate the association: the allele model (AM) (AT vs. GC), codominant model 1 (COD1) (GC/AT vs. GC/GC), codominant model 2 (COD2) (AT/AT vs. GC/GC), codominant model 3 (COD3) (AT/AT vs. GC/AT), dominant model (DM) (AT/AT+ GC/AT vs. GC/GC), recessive model (RM) (AT/AT vs. GC/AT+GC/GC), and over-dominant model (OD) (GC/AT vs. AT/AT+ GC/GC). We also conducted a subgroup analysis in which the controls were divided into hospital-based (HB) and population-based (PB) control populations, while the case or experimental arm included patients with different cancers carrying the TP73 G4C14-A4T14 variant. We also conducted subgroup analysis according to ethnicity in Asian, White, and African populations. The degree of cancer risk was estimated as an odds ratio (OR) with 95% confidence intervals (CIs), and the significance level (Pz) was set to Pz<0.05. A fixed-effects or random-effects model was applied based on the results of the heterogeneity test (Q-test): when heterogeneity was significant (PH<0.10), the random-effects model (DerSimonian–Laird) was applied, and when heterogeneity was non-significant, the fixed-effects model (Mantel–Haenszel) was applied. Visual inspection of funnel plots as well as the results of Egger’s regression and Begg–Mazumdar tests were used to estimate publication bias. Sensitivity analysis was performed to assess the reliability of the results by subtracting the studies one by one. Trial sequential analysis (TSA) was performed using TSA software (version 0.9.5.10 Beta), maintaining an overall 5% risk of a type I error, a relative risk reduction of 20%, and a power of 80%.
Results
Study characteristics
Fifty-five case–control studies2,24–74 including 15,648 cancer cases and 19,159 controls met the eligibility criteria and were finally included in this meta-analysis (Figure 1). A total of 194 studies were excluded after screening the title, abstract and full-text, because of irrelevant information, incomplete genetic data, or duplicate contents. Among the 55 included studies, 11 focused on lung cancer (LC), 10 on gynecological cancers [cervical cancer (CC), endometrial cancer (EM) and ovarian cancer (OVC)], six on colorectal cancer (CRC), five on gastric cancer (GC), four each on esophageal cancer (EC), breast cancer (BC), and oral cancer (OC), three on prostate cancer (PC), one on bladder cancer (UBC), and the others on hepatocellular carcinoma (HCC), non-Hodgkin’s lymphoma (NHL), and neuroblastoma (NB). The included studies were grouped according to ethnicity, including 38 studies of Asian populations, 13 in White populations, three in African populations, and one in a mixed population. In addition, 31 studies recruited controls from HB sources and 24 recruited controls from PB sources. Regarding quality assessment, we determined the NOS score and excluded studies that scored less than 6 points. Detailed demographic information on the included studies is presented in Table 1.
Baseline demographic information of the included studies.
DD, GC/GC; DE, GC/AT; EE, AT/AT; NB, neuroblastoma; BC, breast cancer; EC, esophageal cancer; GC, gastric cancer; SC, skin cancer; CRC, colorectal cancer; LC, lung cancer; NHC, non-Hodgkin’s lymphoma; HNC, head and neck cancer; EMC, endometrial cancer; OC, oral cancer; PC, prostate cancer; CC, cervical cancer; OVC, ovarian cancer; HCC, hepatocellular carcinoma; UBC, bladder cancer; HWE, Hardy–Weinberg equilibrium; PCT, polymerase chain reaction; RFLP, restriction fragment length polymorphism; HRM, high-resolution melting; CTPP, confronting two-pair primers.
Association of TP73 G4C14-A4T14 variant with cancer
We evaluated the overall impact of the TP73 G4C14-A4T14 variant on cancer in a meta-analysis of 55 studies, using seven common genetic models. Five of the genetic models showed significant risk associations with overall cancer, including AM, COD1, COD2, DM, and OD. COD3 and RM did not confirm a significant association between TP73 G4C14-A4T14 and cancer susceptibility (Table 2, Figure 2).
Associations of TP73 G4C14-A4T14 polymorphism with cancer risk in different ethnicities.
AM, allele model; COD1, codominant model 1; COD2, codominant model 2; COD3, codominant model 3; DM, dominant model; RM, recessive model; OD, overdominant model; OR, odds ratio; CI, confidence interval.

Forest plots of results of different genetic models for the association between TP73 G4C14-A4T14 polymorphism and cancer development.
Subgroup analysis based on ethnicity
We compared the results of the seven genetic models among the three ethnic populations: Asian, African, and White (Table 2). There was no significant association between TP73 G4C14-A4T14 and cancer susceptibility in the Asian population. Only the COD2 and RM models showed significant high-risk associations in African populations, while the AM, COD2, DM, and RM models showed significantly increased cancer risks in carriers of the TP73 G4C14-A4T14 in White populations. Forest plots of the results of the AM model for the association of TP73 G4C14-A4T14 with cancer development in different ethnic populations are shown in Figure 3.

Forest plots of results of allele model (AM) on association between TP73 G4C14-A4T14 polymorphism and cancer development in relation to ethnicity.
Subgroup analysis based on cancer types
All the genetic models were applied to analyze the correlation between the TP73 G4C14-A4T14 variant and each cancer type (Table 3). The AM and DM models demonstrated significantly increased susceptibility to gynecological cancers (OVC, CC and EM) in carriers of the TP73 G4C14-A4T14 variant. Five of the genetic models, including AM, COD2, COD3, DM, and RM indicated significant a significant association of the variant with susceptibility to CRC. The G4C14-A4T14 variant was only associated with oral cancer (OC) risk in the COD2 model. Four genetic models implied significant risk susceptibility for HNC, including AM, COD1, DM, and OD model. Cancers in ‘others’ category (HCC+NHL+NB) also showed significant risk association with TP73 G4C14-A4T14 variant in four genetic models- AM, COD1, DM, and OD model. No connection of this polymorphism was found with the risk of LC, EC, GC, BC, UBC+PC, and SC development. Forest plots presenting AM on the cancer type-based association of TP73 G4C14-A4T14 variant with cancer development are presented in Figure 4.
Associations of TP73 G4C14-A4T14 polymorphism with risks of different cancer types.
AM, allele model; COD1, codominant model 1; COD2, codominant model 2; COD3, codominant model 3; DM, dominant model; RM, recessive model; OD, overdominant model; BC, breast cancer; EC, esophageal cancer; GC, gastric cancer; SC, skin cancer; CRC, colorectal cancer; LC, lung cancer; HNC, head and neck cancer; EM, endometrial cancer; OC, oral cancer; PC, prostate cancer; CC, cervical cancer; OVC, ovarian cancer; UBC, bladder cancer; OR, odds ratio; CI, confidence interval.

Forest plots of results of allele model (AM) on association between TP73 G4C14-A4T14 polymorphism and cancer development in relation to cancer type.
Subgroup analysis based on control sources
Among the two types of controls, only studies with HB controls revealed a significant risk susceptibility of the TP73 G4C14-A4T14 variant for cancer development. Five genetic models supported this association namely, the AM, COD1, COD2, DM, and RM models. Studies with PB controls did not reveal any significant risk susceptibility for cancer in relation to the TP73 G4C14-A4T14 variant (Table 4).
Associations of TP73 G4C14-A4T14 polymorphism with cancer risk based on control source.
HB, hospital-based; PB, population-based; OR, odds ratio; CI, confidence interval.
Test of heterogeneity
We determined the level of heterogeneity in this meta-analysis by Q-test. The level of significance was determined by PH and the level of heterogeneity was estimated by I2 statistics. Heterogeneity was significant in the maximum subgroup analysis models (PH<0.1) and random-effects models were applied, while fixed-effects models were used for analyses with PH>0.10. There was significant heterogeneity in all the genetic models for overall cancer. The results for the heterogeneity test of heterogeneity are displayed in Tables 2–4.
Publication bias and sensitivity analysis
Publication bias was determined using Egger’s and Begg–Mazumdar’s tests (Table 5). The funnel plots are shown in Figure 5. We conducted the bias study on the overall analysis with 55 studies using seven genetic models. There was no noticeable visual asymmetry signifying the presence of publication bias. Moreover, the pooled outcomes of this study were considered to be free from publication bias because the p-values were not significant in any of the seven genetic models.
Publication bias analysis.
AM, allele model; COD1, codominant model 1; COD2, codominant model 2; COD3, codominant model 3; DM, dominant model; RM, recessive model; OD, overdominant model.

Funnel plots indicating publication bias of included studies for different models.
To confirm the authenticity of the final findings, we conducted a sensitivity analysis of the studies by sequential elimination of the studies. The impact of each study on the final pooled ORs was checked, and none of the studies affected the pooled ORs. The sensitivity analysis thus confirmed the credibility and robustness of this meta-analysis (Table 6).
Sensitivity analysis of the included studies.
All values presented represent odds ratios with 95% confidence intervals; AM, allele model; COD1, codominant model 1; COD2, codominant model 2; COD3, codominant model 3; DM, dominant model; RM, recessive model; OD, overdominant model.
TSA outcomes
The TSA plots (Figure 6) indicated that the Z-curves exceeded the required information size in the overall population and in Whites and Asians, indicating that the total cases and controls were sufficient to confirm the outcomes, and no further studies were required. However, the Z-curve did not exceed the required information size in Africans, and further studies are therefore required to confirm the outcome.

Trial sequential analysis (TSA) of association between TP73 G4C14-A4T14 polymorphism and cancer risk in allele model. (a) Overall population; (b) Whites; (c) Asians; and (d) Africans.
Discussion
The TP73 gene encodes multiple protein isoforms with similar or opposite functions. The protein shows almost 63% homology with the tumor suppressor protein p53 in terms of its DNA-binding capability, oligomerization of the domains, and gene transactivation. 21 The protein isoforms of p73 arise from the utilization of different promoter sites and alternative mRNA splicing. Two common isoforms of p73 are TAp73 (TA domain present) and ΔNp73 (TA domain absent). Of these, TAp73 mimics the tumor suppression activities of p53 by inducing apoptosis, arresting G1 cell cycle checkpoint, and regulating the transcription of p53-related genes, while ΔNp73 exerts opposing functions by promoting oncogenic activities due to the lack of TA domain. ΔNp73 acts as an inhibitor of both p53 and p73 proteins.10,22,81,82 The TP73 G4C14-A4T14 variant of exon 2 potentially influences the translation of p73 by forming a stem-loop structure. 21 A recent study identified a significant association between the TP73 G4C14-A4T14 variant and ΔNp73 tumoral immunostaining in 91% of cancer patients. 83 High expression of ΔNp73 in carriers of the G4C14-A4T14 variant demonstrates the potential role of this polymorphism in carcinogenesis. Numerous studies have evaluated this association, but most of the findings have been inconclusive. We carried out the current meta-analysis to address these inconsistencies, and showed that the TP73 G4C14-A4T14 variant could significantly elevate the risk of cancer.
Recent studies have examined the association between the G4C14-A4T14 variant and lung cancer risk in patients with non-small cell lung cancer (NSCLC). Most of these studies reported no significant risk association, and no significant difference in the frequency of the variant between patients and controls. However, some studies found that this variant was associated with a reduced risk of NSCLC in AT/AT carriers compared with GC/GC carriers.43,45,47,59,60,64,71 In contrast, other studies showed that this polymorphism could significantly increase the risk of NSCLC among the variant GC/AT and AT/AT genotype carriers, and a high GC content increased the risk. TP73 and MDM2 variants jointly increase the risk of lung cancer, depending on the number of variant alleles.32,54,55,68 The G4C14-A4T14 variant also increased the risk of gynecological cancers, with a two-fold increase in susceptibility to high-grade squamous intraepithelial lesion in women carrying the TP73 AT allele. 26 The risks of CC and EM cancers were also increased among carriers of the TP73 polymorphism who were passive smokers.27,28,41,58 However, some association studies failed to detect any significant association between TP73 genotype and tumor stage, histological type, or lymph node metastasis in patients with gynecological cancers.29,48,52 Regarding CRC, AT/AT homozygous genotype of TP73 was associated with an increased risk of CRC and a poor prognosis, whereas AT allele carriers had a better prognosis. However, another study failed to observe any significant association between the TP73 GC/AT variant genotype and allele distribution and clinical parameters of CRC.30,35,42,44,57 Some previous studies identified the AC/GT genotype of G4C14-A4T14 as a significant risk factor for GCs, although other studies found no such association.31,36,40,67,70
Decreased expression of p73 mRNA was identified in both inflammatory and non-inflammatory BC cells compared with normal breast epithelial cells, indicating that this variant might increase the risk of BC by reducing the expression of p73. A recent study postulated that the TP73 GC/AT and AT/AT genotypes could increase the susceptibility to BC, while another study found that the GC/GC genotype was associated with an increased risk of triple-negative BC,2,37,66 and yet another study found no significant association between this polymorphism and BC. 62 Similar findings were observed in EC studies with contradictory conclusions.31,46,67,72 TP73 G4C14-A4T14 was recently identified as a risk factor for OC development.33,51,74 Although the risk variant was associated with an increased risk of UBC, it showed a significant inverse relationship with PC.39,50,56 Among other studies of the association between this variant and OC, SC, HNC, and other cancers (HCC+NHL+ NB), most identified G4C14-A4T14 polymorphism as a risk variant for cancer.24,25,49,53,61,65,69,73
The current meta-analysis of 55 case–control studies found that the TP73 G4C14-A4T14 variant was linked to an increased risk of overall cancer development. Five of the tested genetic models (AM, COD1, COD2, DM, and OD) showed a significantly increased risk of overall cancer (1.10, 1.09, 1.18, 1.11, and 1.08-fold, respectively). Subgroup analysis based on ethnicity also showed a significant association between the variant and cancer risk in Africans in two genetic models (COD2, 2.12-fold; RM, 2.00-fold), while four genetic models reported significantly elevated cancer risks among TP73 G4C14-A4T14 variant carriers in White populations (AM, 1.14-fold; COD2, 1.30-fold; DM, 1.15-fold; RM, 1.24-fold). In terms of specific cancers, sub-group analysis identified significant associations between the TP73 G4C14-A4T14 variant and the risks of gynecological cancer (OVC, CC and EM), CRC, OC, HNC, and other cancers (HCC+NHL+NB). An increased susceptibility to gynecological cancers was reported in two genetic models (AM, OR = 1.16; DM, OR = 1.18), an increased risk of CRC in five genetic models (AM, OR = 1.26; COD2, OR = 1.97; COD3, OR = 1.81; DM, OR = 1.20; RM, OR = 1.89). The G4C14-A4T14 variant was only associated with OC according to the COD2 model (OR = 1.51) and with HNC according to the AM (OR = 1.25), COD1 (OR = 1.39), DM (OR = 1.36), and OD models (OR = 1.38). The variant was significantly associated with ‘other cancers’ according to the AM (1.57-fold), COD1 (1.80-fold), DM (1.82-fold), and OD (1.67-fold) genetic association models. Moreover, studies with HB controls revealed significant susceptibility of G4C14-A4T14 variant carriers to cancer according to the AM (OR = 1.13), COD1 (OR = 1.11), COD2 (OR = 1.27), DM (OR = 1.14), and RM models (OR = 1.22).
Some previous systematic meta-analyses examined the relationship between various cancer types and the TP73 G4C14-A4T14 variant. Yu and colleagues performed a meta-analysis of 23 case–control studies and reported that this polymorphism might be significantly associated with cancer risk in Asian and White populations. 75 Another meta-analysis of 27 case–control studies concluded that carriers of the AT/AT genotype might be at high-risk of developing cancer among Asians and Whites. 76 A further meta-analysis of five case–control studies in 2017 confirmed that the polymorphism was associated with CC risk, but the number of included studies was small. 77 Meng et al. performed a recent meta-analysis of 36 case–control studies and found that the TP73 G4C14-A4T14 variant was associated with an increased cancer risk, especially among Whites. 78 In contrast to these previous meta-analyses, the current meta-analysis included a large number of studies (55 case–control studies) that provided more consistent outcomes than previous studies. Moreover, we validated the stability and consistency of our findings by carrying out heterogeneity, publication bias, and sensitivity analyses, as well as TSA. The results of this study provide strong evidence for an association between the TP73 G4C14-A4T14 variant and cancer development, by successfully avoiding publication bias. The quality of the included studies was also evaluated by NOS scoring, and low-quality studies were excluded to maintain the robustness of the final findings.
Although the present meta-analysis was conducted carefully, some limitations could not be avoided. The number of studies included in some of the subgroups was small, due to the lack of available information. In addition, some basic information on both the patients and controls was lacking, such as age, sex, medication, and body mass index, which could have further enriched the analysis. Further analyses should therefore be conducted, including more studies, to confirm the relationship between TP73 G4C14-A4T14 and cancer risk.
Conclusion
This updated meta-analysis provides strong evidence indicating that the TP73 G4C14-A4T14 variant may elevate the overall cancer risk, especially in White and African populations. Carriers of the G4C14-A4T14 variant have increased risks of developing gynecological cancers, such as cervical, ovarian, and endometrial cancer, as well as colorectal, head and neck, and oral cancers, non-Hodgkin’s lymphoma, and neuroblastoma. Moreover, studies recruiting HB controls revealed a significant association between the G4C14-A4T14 variant and cancer risk.
Footnotes
Acknowledgments
The authors are thankful to the Department of Pharmacy, Noakhali Science and Technology University, for supporting the authors during the study.
Availability of data and material
All relevant data that support the study's results are accessible upon request from the corresponding author.
Author contributions
Mohammad Safiqul Islam: conceptualization, supervision, data analysis, software; Sarah Jafrin and Md. Abdul Aziz: literature search; Sarah Jafrin: writing- original draft preparation, methodology; Md. Abdul Aziz: writing – original draft preparation, methodology; writing – reviewing and editing; Mohammad Safiqul Islam: writing – reviewing and editing.
Declaration of conflicting interests
The authors declared no potential conflicts of interest
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
