The Prevalence of Pharmacogenomics Variants and Their Clinical Relevance Among the Pakistani Population.

Background: Pharmacogenomics (PGx), forming the basis of precision medicine, has revolutionized traditional medical practice. Currently, drug responses such as drug efficacy, drug dosage, and drug adverse reactions can be anticipated based on the genetic makeup of the patients. The pharmacogenomic data of Pakistani populations are limited. This study investigates the frequencies of pharmacogenetic variants and their clinical relevance among ethnic groups in Pakistan. Methods: The Pharmacogenomics Knowledge Base (PharmGKB) database was used to extract pharmacogenetic variants that are involved in medical conditions with high (1A + 1B) to moderate (2A + 2B) clinical evidence. Subsequently, the allele frequencies of these variants were searched among multiethnic groups of Pakistan (Balochi, Brahui, Burusho, Hazara, Kalash, Pashtun, Punjabi, and Sindhi) using the 1000 Genomes Project (1KGP) and ALlele FREquency Database (ALFRED). Furthermore, the published Pharmacogenomics literature on the Pakistani population was reviewed in PubMed and Google Scholar. Results: Our search retrieved (n = 29) pharmacogenetic genes and their (n = 44) variants with high to moderate evidence of clinical association. These pharmacogenetic variants correspond to drug-metabolizing enzymes (n = 22), drug-metabolizing transporters (n = 8), and PGx gene regulators, etc. (n = 14). We found 5 pharmacogenetic variants present at >50% among 8 ethnic groups of Pakistan. These pharmacogenetic variants include CYP2B6 (rs2279345, C; 70%-86%), CYP3A5 (rs776746, C; 64%-88%), FLT3 (rs1933437, T; 54%-74%), CETP (rs1532624, A; 50%-70%), and DPP6 (rs6977820, C; 61%-86%) genes that are involved in drug response for acquired immune deficiency syndrome, transplantation, cancer, heart disease, and mental health therapy, respectively. Conclusions: This study highlights the frequency of important clinical pharmacogenetic variants (1A, 1B, 2A, and 2B) among multi-ethnic Pakistani populations. The high prevalence (>50%) of single nucleotide pharmacogenetic variants may contribute to the drug response/diseases outcome. These PGx data could be used as pharmacogenetic markers in the selection of appropriate therapeutic regimens for specific ethnic groups of Pakistan.

Germline Mutation Analysis in Sporadic Breast Cancer Cases With Clinical Correlations.

Demographics for breast cancers vary widely among nations. The frequency of germline mutations in breast cancers, which reflects the hereditary cases, has not been investigated adequately and accurately in highly-consanguineous Pakistani population. In the present discovery case series, germ-line mutations in twenty-seven breast cancer candidate genes were investigated in eighty-four sporadic breast cancer patients along with the clinical correlations. The germ-line variants were also assessed in two healthy gender-matched controls. The clinico-pathological features were evaluated by descriptive analysis and Pearson χ2 test (with significant p-value <0.05). The most frequent parameters associated with hereditary cancer cases are age and ethnicity. Therefore, the analyses were stratified on the basis of age (≤40 years vs. >40 years) and ethnicity. The breast cancer gene panel assay was carried out by BROCA, which is a genomic capture, massively parallel next generation sequencing assay on Illumina Hiseq2000 with 100bp read lengths. Copy number variations were determined by partially-mapped read algorithm. Once the mutation was identified, it was validated by Sanger sequencing. The ethnic analysis stratified on the basis of age showed that the frequency of breast cancer at young age (≤40 years) was higher in Sindhis (n = 12/19; 64%) in contrast to patients in other ethnic groups. Majority of the patients had stage III (38.1%), grade III (50%), tumor size 2-5 cm (54.8%), and invasive ductal carcinoma (81%). Overall, the analysis revealed germ-line mutations in 11.9% of the patients, which was not significantly associated with younger age or any particular ethnicity. The mutational spectrum was restricted to three genes: BRCA1, BRCA2, and TP53. The identified mutations consist of seven novel germ-line mutations, while three mutations have been reported previously. All the mutations are predicted to result in protein truncation. No mutations were identified in the remaining twenty-four candidate breast cancer genes. The present study provides the framework for the development of hereditary-based preventive and treatment strategies against breast cancers in Pakistani population.

Absence of Glutathione S-Transferase Theta 1 Gene Is Significantly Associated With Breast Cancer Susceptibility in Pakistani Population and Poor Overall Survival in Breast Cancer Patients: A Case-Control and Case Series Analysis.

PURPOSE: Deletion of Glutathione S-Transferase Theta 1 (GSTT1) encoding gene is implicated in breast cancer susceptibility, clinical outcomes, and survival. Contradictory results have been reported in different studies. The present investigation based on a representative Pakistani population evaluated the GSTT1-absent genotype in breast cancer risk and prognosis. METHODS: A prospective study comprising case-control analysis and case series analysis components was designed. Peripheral blood samples were collected from enrolled participants. After DNA extraction, GSTT1 genotyping was carried out by a multiplex PCR with ß-globin as an amplification control. Association evaluation of GSTT1 genotypes with breast cancer risk, specific tumor characteristics, and survival were the primary endpoints. RESULTS: A total of 264 participants were enrolled in the molecular investigation (3 institutions). The study included 121 primary breast cancer patients as cases and 143 age-matched female subjects, with no history of any cancer, as controls. A significant genetic association between GSTT1-absent genotype and breast cancer susceptibility (p-value: 0.03; OR: 2.13; 95% CI: 1.08-4.29) was reported. The case-series analysis showed lack of association of GSTT1 genotypes with menopause (p-value: 0.86), tumor stage (p-value: 0.12), grade (p-value: 0.32), and size (p-value: 0.07). The survival analysis revealed that GSTT1-absent genotype cases had a statistically significant shorter overall survival (OS) than those with the GSTT1-present genotype cases (mean OS: 23 months vs 33 months). The HR (95% CI) for OS in patients carrying GSTT1-absent genotype was 8.13 (2.91-22.96) when compared with the GSTT1-present genotype. CONCLUSIONS: The present study is the first report of an independent significant genetic association between GSTT1-absent genotype and breast cancer susceptibility in a Pakistani population. It is also the foremost report of the association of this genotype with OS in breast cancer cases. Upon further validation, GSTT1 variation may serve as a marker for devising better population-specific strategies. The information may have translational implications in the screening and treatment of breast cancers.

Association of specific single nucleotide variants (SNVs) in the promoter and 3'-Untranslated region of Vascular Endothelial growth factor (VEGF) gene with risk and higher tumour grade of head and neck cancers.

BACKGROUND: Head and Neck Cancers (HNCs)comprise one of the most frequent cancers in South-Asian region. Vascular Endothelial Growth Factor (VEGF) has a potent role in tumorigenesis and metastasis. Certain common single nucleotide variants (SNVs) in the highly polymorphic VEGF gene are correlated with variations in VEGF functions. The data for these SNVs in HNCs is scarce for South Asian populations. The present study addresses this shortfall. It investigates the association of two VEGF SNVs, -2578C/A (rs699947) in the promoter region and + 936C/T (rs3025039) in 3'-UTR, with the risk of HNCs and tumour characteristics. METHODS: The study comprised 323 participants with 121 HNC patients and 202 controls. Germline DNA was isolated from peripheral blood samples. PCR-RFLP methods were optimized and validated by Sanger sequencing. After Hardy-Weinberg evaluation, the independent associations were analyzed under the assumptions of different genetic models. The χ2 test of independence or Fisher's Exact test (significant p-values at < 0.05) were performed and ORs (odds ratios) with 95% confidence interval were tabulated. RESULTS: VEGF -2578 A-allele, CA + AA, and AA genotypes had significant protective association against HNCs. The respective ORs were: 0.651 (0.469-0.904), 0.613 (0.381 - 0.985), and 0.393 (0.193-0.804). VEGF + 936 T-allele, CT, and CT + TT genotypes had significantly increased susceptibility for HNCs. The respective ORs were 1.882 (1.001 - 3.536), 2.060 (1.035 - 4.102), and 2.023 (1.032 - 3.966). Additionally, VEGF + 936 CT and CT + TT genotypes showed significant associations with higher tumour grade (p-values < 0.029, and < 0.037, respectively). CONCLUSION: The present study is the foremost report of independent and unique associations of the investigated VEGF SNVs with HNCs.

Analysis of the glutathione S-transferase genes polymorphisms in the risk and prognosis of renal cell carcinomas. Case-control and meta-analysis.

BACKGROUND: The Glutathione S-transferases (GSTs) genes deletion polymorphisms have been associated with the progression of several cancers. The association studies between the 2 GSTs (GSTM1 and GSTT1) null polymorphisms with the susceptibility to renal cell carcinoma (RCC) have been inconclusive. Therefore, with the inclusion of our own data, we performed a comprehensive meta-analysis to assess the association between these 2 polymorphisms and the risk of RCC. METHODS: A systematic literature search was carried out for studies published in the PubMed, EMBASE, Cochrane library, and Google Scholar from 1997 to December 2014. Results were stated as pooled odds ratios (ORs) for nonparametric data after heterogeneity analysis with 95% CI using fixed effect or random effect model. RESULTS: We systematically selected 13 relevant studies after thorough searches from the databases. Data showed no association between the GSTM1 and the GSTT1 null genotypes and the risk of RCC (OR = 1.01; CI: 0.92-1.11; P = 0.89 for GSTM1 and OR = 1.14; CI: 0.91-1.42; P = 0.25 for GSTT1). No association was found when the data were stratified according to the geographical/ethnic basis, source of control, and the risk factor evaluation. Subgroup analysis of occupational exposure to pesticides showed an inverse association of the active genotypes of both GSTM1 and GSTT1 polymorphisms with the exposed group of RCC (P<0.00001 and P<0.00001, respectively). The combined null genotype of the GSTM1/GSTT1 significantly increased the susceptibility to RCC by 1.4-fold (P = 0.001). This association remained significant for the Asian populations in subgroup analysis (OR = 1.8; CI: 1.30-2.49; P = 0.0004). CONCLUSION: In conclusion, this meta-analysis suggests that the 2 GSTs deletion polymorphisms independently have no association with the risk of RCC. However, combination of both deletions increases the risk of developing the RCC.

Unique molecular alteration patterns in von Hippel-Lindau (VHL) gene in a cohort of sporadic renal cell carcinoma patients from Pakistan.

BACKGROUND: Renal cell carcinoma (RCC) is the most frequent form of kidney cancer in adults. Somatic mutations that inactivate the von Hippel-Lindau (VHL) gene are the most common cause of RCC. The frequencies of molecular changes in the VHL gene in RCCs vary among different populations. So far, a single chromosomal-based study has been reported from a South Asian population. This report presents, for the first time, the somatic changes and promoter hypermethylation in VHL in a cohort of 300 RCC patients from Pakistan. METHODS: To identify mutations in the VHL gene, direct DNA sequencing was carried out. Epigenetic silencing was investigated by using methylation-specific polymerase chain reaction. RESULTS: Our data showed molecular alterations in the VHL gene in 163 (54%) renal cell carcinoma patients. Somatic mutations were found in 87 (29%) patients and 35 novel mutations were identified. VHL promoter hyper-methylation analysis showed epigenetic changes in 106 (35%) out of 300 patients. Patients who had no evidence of molecular alterations in the VHL gene were significantly younger than patients who carried some molecular change. Molecular alterations in the VHL gene were not restricted to clear-cell RCCs (ccRCCs). CONCLUSIONS: This is the first report that identifies molecular aberrations in the VHL gene from a South Asian population. The frequency of somatic mutation is lower and that of promoter hypermethylation is higher when compared with data from other parts of the world. The data has important implications in the population-specific application of tailored preventive and therapeutic regimens in non-familial RCCs.

Polymorphisms in the methylene tetrahydrofolate reductase gene and their unique combinations are associated with an increased susceptibility to the renal cancers.

AIMS: Two single nucleotide polymorphisms in the methylene tetrahydrofolate reductase (MTHFR) gene, 677C/T and 1298A/C, encode the thermolabile isoforms of the MTHFR enzyme that adversely affect the folic acid metabolic pathway. In the present study, these polymorphisms were investigated for their associations with the risk and prognosis of the renal cell carcinomas (RCCs) in Pakistani patients. RESULTS: The study included 168 RCC patients and 178 controls. The polymorphisms were analyzed by the polymerase chain reaction-restriction fragment length polymorphism method. Statistical analysis revealed that the C-allele and homozygous C genotype of the MTHFR 1298A/C polymorphism were significantly correlated with the risk of RCCs (odds ratio [OR]=1.60; 95% confidence interval [CI]=1.1-2.34 and OR=3.26; 95% CI=1.27-8.37, respectively). The combined genotype analysis showed that the 677CC+1298CC combination greatly increased the susceptibility to RCCs (OR=8.34; 95% CI=2.7-25.7). The 677CT+1298AA and 677CC+1298CA combinations were also associated with an increased risk of RCC (OR=3.21; 95% CI=1.3-7.8 and OR=2.45; 95% CI=1.3-4.6, respectively). The combined genotype effects were also evident in a semiparametric expectation-maximization-based haplotype analysis. CONCLUSION: The results presented here indicate that the two MTHFR gene polymorphisms are significantly associated with the risk of RCCs in a cohort of Pakistani patients and may be useful as susceptibility markers in other populations of the world as well.

Association of a single-nucleotide polymorphism in the promoter region of the VEGF gene with the risk of renal cell carcinoma.

AIMS: Vascular endothelial growth factor (VEGF) protein plays an important role in tumor development and progression. Polymorphisms in the VEGF gene may lead to over- or underexpression of the protein and may be associated with either risk or progression of malignancy. The aim of this case-control study is to identify and quantify the correlation between VEGF polymorphisms and renal cell carcinoma (RCC). RESULTS: Restriction fragment length polymorphism methods were used for the analysis of VEGF polymorphisms at -2578 and +936 positions in the promoter and 3'-untranslated regions, respectively. The VEGF -2578 A-allele was associated with an increased risk of RCC (odds ratio: 1.6; 95% CI: 1.2-2.3) and A-carrier genotypes were strongly correlated (odds ratio: 2.7; 95% CI: 1.5-4.7) with higher risk. Comparison of VEGF +936 C/T polymorphism between patient and control groups revealed no association with renal carcinoma. Both VEGF -2578 C/A and VEGF +936 C/T polymorphisms showed no significant association with the histopathological parameters of RCC. CONCLUSIONS: This study shows that VEGF -2578 A-allele and A-carrier genotypes are associated with an increased risk of RCC. In groups with higher incidence of RCC, a screening test for this polymorphism may be recommended in conjunction with other established markers.

Altered gene expression profiles define pathways in colorectal cancer cell lines affected by celecoxib.

It is well established that celecoxib, a selective inhibitor of cyclooxygenase-2 (COX-2) and a tested chemopreventive agent, has several COX-2-independent activities. In an attempt to better understand COX-2-independent molecular mechanisms underlying the chemopreventive activity of celecoxib, we did global transcription profiling of celecoxib-treated COX-2-positive and COX-2-deficient colorectal cancer cell lines. Celecoxib treatment resulted in significantly altered expression levels of over 1,000 to 3,000 transcripts in these cell lines, respectively. A pathway/functional analysis of celecoxib-affected transcripts, using Gene Ontology and Biocarta Pathways and exploring biological association networks, revealed that celecoxib modulates expression of numerous genes involved in a variety of cellular processes, including metabolism, cell proliferation, apoptotic signaling, cell cycle check points, lymphocyte activation, and signaling pathways. Among these processes, cell proliferation and apoptotic signaling consistently ranked as the highest-scoring Gene Ontology terms and Biocarta Pathways in both COX-2 expresser and nonexpresser cell lines. Altered expression of many of the genes by celecoxib was confirmed by quantitative PCR and at the protein level by Western blotting. Many novel genes emerged from our analysis of global transcription patterns that were not previously reported to be affected by celecoxib. In the future, in-depth work on selected genes will determine if these genes may serve as potential molecular targets for more effective chemopreventive strategies.