Pharmacogenetics of pemetrexed combination therapy in lung cancer: pathway analysis reveals novel toxicity associations

Identification of molecular biomarkers associated with non-small-cell lung carcinoma (NSCLC) using whole-exome sequencing

ObjectivesSignificant progress has been made in the treatment of patients with pulmonary adenocarcinoma (ADCA) based on molecular profiling. However, no such molecular target exists for squamous cell carcinoma (SQCC). An exome sequence may provide new markers for personalized medicine for lung cancer patients of all subtypes. The current study aims to discover new genetic markers that can be used as universal biomarkers for non-small cell lung cancer (NSCLC).MethodsWES of 19 advanced NSCLC patients (10 ADCA and 9 SQCC) was performed using Illumina HiSeq 2000. Variant calling was performed using GATK HaplotypeCaller and then the impacts of variants on protein structure or function were predicted using SnpEff and ANNOVAR. The clinical impact of somatic variants in cancer was assessed using cancer archives. Somatic variants were further prioritized using a knowledge-driven variant interpretation approach. Sanger sequencing was used to validate functionally important variants.ResultsWe identified 24 rare single-nucleotide variants (SNVs) including 17 non-synonymous SNVs, and 7 INDELs in 18 genes possibly linked to lung carcinoma. Variants were classified as known somatic (n = 10), deleterious (n = 8), and variant of uncertain significance (n = 6). We found TBP and MPRIP genes exclusively associated with ADCA subtypes, FBOX6 with SQCC subtypes and GPRIN2, KCNJ18 and TEKT4 genes mutated in all the patients. The Sanger sequencing of 10 high-confidence somatic SNVs showed 100% concordance in 7 genes, and 80% concordance in the remaining 3 genes.ConclusionsOur bioinformatics analysis identified KCNJ18, GPRIN2, TEKT4, HRNR, FOLR3, ESSRA, CTBP2, MPRIP, TBP, and FBXO6 may contribute to progression in NSCLC and could be used as new biomarkers for the treatment. The mechanism by which GPRIN2, KCNJ12, and TEKT4 contribute to tumorigenesis is unclear, but our results suggest they may play an important role in NSCLC and it is worth investigating in future.

Genetic polymorphisms and platinum-induced hematological toxicity: a systematic review

Background

Platinum-based chemotherapy bring severe hematological toxicity that can lead to dose reduction or discontinuation of therapy. Genetic variations have been reported to influence the risk and extent of hematological toxicity; however, the results are controversial and a comprehensive overview is lacking. This systematic review aimed to identify genetic biomarkers of platinum-induced hematological toxicity.

Method

Pubmed, Embase and Web of science database were systematically reviewed for studies that evaluated the association of genetic variants and platinum-related hematological toxicity in tumor patients with no prior history of chemotherapy or radiation, published from inception to the 28th of January 2022. The studies should have specific toxicity scoring system as well as defined toxicity end-point. The quality of reporting was assessed using the Strengthening the Reporting of Genetic Association Studies (STREGA) checklist. Results were summarized using narrative synthesis.

Results

83 studies were eligible with over 682 single-nucleotide polymorphisms across 110 genes. The results are inconsistent and diverse with methodological issues including insufficient sample size, population stratification, various treatment schedule and toxicity end-point, and inappropriate statistics. 11 SNPs from 10 genes (ABCB1 rs1128503, GSTP1 rs1695, GSTM1 gene deletion, ERCC1 rs11615, ERCC1 rs3212986, ERCC2 rs238406, XPC rs2228001, XPCC1 rs25487, MTHFR rs1801133, MDM2 rs2279744, TP53 rs1042522) had consistent results in more than two independent populations. Among them, GSTP1 rs1695, ERCC1 rs11615, ERCC1 rs3212986, and XRCC1 rs25487 present the most promising results.

Conclusion

Even though the results are inconsistent and several methodological concerns exist, this systematic review identified several genetic variations that deserve validation in well-defined studies with larger sample size and robust methodology.

Systematic Review Registration

https://www.crd.york.ac.uk/, identifier CRD42021234164.

MTHFR polymorphism as a predictive biomarker for gastrointestinal and hematological toxicity in North Indian adenocarcinoma patients

In the present study, we investigated the relationship between the methylenetetrahydrofolate reductase (MTHFR) C677T and A1298C polymorphisms and overall survival, toxicity and treatment response for North Indian adenocarcinoma patients. The polymorphisms of MTHFR gene in north Indian adenocarcinoma patients were assessed by PCR-RFLP. Our data observed that patients with mutant genotype (C/C) for 1298 A>C) polymorphism showed higher trend of median survival time compared to patients bearing the wild type genotype (A/A) (MST= 13.93 vs. 7.97, p=0.12). Further, we observed patients with the heterozygous genotype for A1298C polymorphism had 12-fold risk of diarrhea (AOR =12.54, 95% CI = 1.54-101.86, p=0.018). The patients with heterozygous genotype (CT) of the C677T polymorphism had 5.34-fold increased risk of developing neutropenia (AOR=5.34, 95% CI=1.49-19.06, p=0.009). Our results suggest that MTHFR polymorphisms are associated with hematological toxicity. MTHFR polymorphism might impact the development of pemetrexed and platinum-related toxicities but not as a clinical predictor of efficiency.

Genetic Polymorphisms of Very Important Pharmacogene Variants in the Blang Population from Yunnan Province in China

Background

We aimed to enrich the pharmacogenomic information of a Blang population (BP) from Yunnan Province in China.

Methods

We genotyped 55 very important pharmacogene (VIP) variants from the PharmGKB database and compared their genotype distribution (GD) in a BP with that of 26 populations by the χ² test. The minor allele frequency (MAF) distribution of seven significantly different single-nucleotide polymorphisms (SNPs) was conducted to compare the difference between the BP and 26 other populations.

Results

Compared with the GD of 55 loci in the BP, among 26 studied populations, GWD, YRI, GIH, ESN, MSL, TSI, PJL, ACB, FIN and IBS were the top-10 populations, which showed a significantly different GD >35 loci. CHB, JPT, CDX, CHS, and KHV populations had a significantly different GD <20 loci. A GD difference of 27–34 loci was found between the BP and 11 populations (LWK, CEU, ITU, STU, PUR, CLM, GBR, ASW, BEB, MXL and PEL). The GD of five loci (rs750155 (SULT1A1), rs4291 (ACE), rs1051298 (SLC19A1), rs1131596 (SLC19A1) and rs1051296 (SLC19A1)) were the most significantly different in the BP as compared with that of the other 26 populations. The genotype frequency of rs1800764 (ACE) and rs1065852 (CYP2D6) was different in all populations except for PEL and LWK, respectively. MAFs of rs1065852 (CYP2D6) and rs750155 (SULT1A1) showed the largest fluctuation between the BP and SAS, EUR, AFR and AMR populations.

Conclusion

Our data can provide theoretical guidance for safe and efficacious personalized drug use in the Blang population.

Systematic review and meta-analysis of the predictive power of MTHFR polymorphisms for pemetrexed drug efficacy and toxicity in non-small cell lung cancer patients

We conducted a meta-analysis to determine if MTHFR polymorphisms are effective biomarkers for non-small cell lung cancer (NSCLC) patient survival and pemetrexed (PEM) treatment toxicity. Because of data heterogeneity, fixed or random effects models were chosen, and pooled HRs and 95% confidence intervals (CIs) were calculated. No correlation between MTHFR 677 C > T polymorphism and progression-free survival (PFS) or overall survival (OS) was detected in NSCLC patients; however, patients with the T allele benefited more than those with the wild-type allele. Two papers reported hematologic toxicity of single-agent PEM treatment in patients with the MTHFR 677 C > T polymorphism. However, data on MTHFR polymorphisms and toxicity could not be combined, even though publication bias and sensitivity analysis results were stable and reliable. We conclude that the MTHFR 677 C > T polymorphism could not predict PEM efficacy in NSCLC patients; however, the T allele may increase the risk of haematological toxicity. A large-scale clinical trial is recommended.

Genetic analysis of pharmacogenomic VIP variants in the Wa population from Yunnan Province of China

Background

The variation of drug responses and target does among individuals is mostly determined by genes. With the development of pharmacogenetics and pharmacogenomics, the differences in drug response between different races seem to be mainly caused by the genetic diversity of pharmacodynamics and pharmacokinetics genes. Very important pharmacogenetic (VIP) variants mean that genes or variants play important and vital roles in drug response, which have been listed in pharmacogenomics databases, such as Pharmacogenomics Knowledge Base (PharmGKB). The information of Chinese ethnic minorities such as the Wa ethnic group is scarce. This study aimed to uncover the significantly different loci in the Wa population in Yunnan Province of China from the perspective of pharmacogenomics, to provide a theoretical basis for the future medication guidance, and to ultimately achieve the best treatment in the future.

Results

In this study, we recruited 200 unrelated healthy Wa adults from the Yunnan province of China, selected 52 VIP variants from the PharmGKB for genotyping. We also compared the genotype frequency and allele distribution of VIP variants between Wa population and the other 26 populations from the 1000 Genomes Project (http://www.1000Genomes.org/). Next, χ² test was used to determine the significant points between these populations. The study results showed that compared with the other 26 population groups, five variants rs776746 (CYP3A5), rs4291 (ACE), rs3093105 (CYP4F2), rs1051298 (SLC19A1), and rs1065852 (CYP2D6) had higher frequencies in the Wa population. The genotype frequencies rs4291-TA, rs3093105-CA, rs1051298-AG and rs1065852-GA were higher than those of the other populations, and the allele distributions of rs4291-T and rs3093105-C were significantly different. Additionally, the difference between the Wa ethnic group and East Asian populations, such as CDX, CHB, and CHS, was the smallest.

Conclusions

Our research results show that there is a significant difference in the distribution of VIP variants between the Wa ethnic group and the other 26 populations. The study results will have an effect on supplementing the pharmacogenomics information for the Wa population and providing a theoretical basis for individualised medication for the Wa population.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-00999-8.

Identification and Prognostic Value Exploration of Radiotherapy Sensitivity-Associated Genes in Non-Small-Cell Lung Cancer

Background

Non-small-cell lung cancer (NSCLC) is a prevalent malignancy with high mortality and poor prognosis. The radiotherapy is one of the most common treatments of NSCLC, and the radiotherapy sensitivity of patients could affect the individual prognosis of NSCLC. However, the prognostic signatures related to radiotherapy response still remain limited. Here, we explored the radiosensitivity-associated genes and constructed the prognostically predictive model of NSCLC cases.

Methods

The NSCLC samples with radiotherapy records were obtained from The Cancer Genome Atlas database, and the mRNA expression profiles of NSCLC patients from the GSE30219 and GSE31210 datasets were obtained from the Gene Expression Omnibus database. The Weighted Gene Coexpression Network Analysis (WGCNA), univariate, least absolute shrinkage and selection operator (LASSO), multivariate Cox regression analysis, and nomogram were conducted to identify and validate the radiotherapy sensitivity-related signature.

Results

WGCNA revealed that 365 genes were significantly correlated with radiotherapy response. LASSO Cox regression analysis identified 8 genes, including FOLR3, SLC6A11, ALPP, IGFN1, KCNJ12, RPS4XP22, HIST1H2BH, and BLACAT1. The overall survival (OS) of the low-risk group was better than that of the high-risk group separated by the Risk Score based on these 8 genes for the NSCLC patients. Furthermore, the immune infiltration analysis showed that monocytes and activated memory CD4 T cells had different relative proportions in the low-risk group compared with the high-risk group. The Risk Score was correlated with immune checkpoints, including CTLA4, PDL1, LAG3, and TIGIT.

Conclusion

We identified 365 genes potentially correlated with the radiotherapy response of NSCLC patients. The Risk Score model based on the identified 8 genes can predict the prognosis of NSCLC patients.

Impact of genetic factors on platinum-induced gastrointestinal toxicity

Severe gastrointestinal (GI) toxicity is a common side effect after platinum-based chemotherapy. The incidence and severity of GI toxicity vary among patients with the same chemotherapy. Genetic factors involved in platinum transport, metabolism, detoxification, DNA repair, cell cycle control, and apoptosis pathways may account for the interindividual difference in GI toxicity. The influence of gene polymorphisms in the platinum pathway on GI toxicity has been extensively analyzed. Variations in study sample size, ethnicity, design, treatment schedule, dosing, endpoint definition, and assessment of toxicity make it difficult to precisely interpret the results. Hence, we conducted a review to summarize the most recent pharmacogenomics studies of GI toxicity in platinum-based chemotherapy and identify the most promising avenues for further research.

Mutational Profile of Malignant Pleural Mesothelioma (MPM) in the Phase II RAMES Study

Simple Summary

Malignant pleura Mesothelioma (MPM) is an aggressive cancer arising from the mesothelial cells of the pleura. About 80% of mesothelioma cases are linked to asbestos exposure; the remainder may be related to prior chest radiation, genetic predisposition or spontaneous occurrence. Understanding the genetic alterations that drive MPM is critical for successful development of diagnostics, prognostics and personalized therapeutic modalities. Because MPM is rare, genomic studies are limited and have typically involved a small number of samples. The aim of our study was to understand the mutational landscape of MPM. Our analysis identified significantly mutated genes. This study on a rare tumor type will be important for patients “in real life”.

Abstract

Purpose: Malignant pleural mesothelioma (MPM) is an aggressive cancer. Data are not available in prospective trials on correlations between genetic alterations and outcomes of therapies. In this study, we assessed the genetic profile of MPM patients (pts) in tissue samples. Patients and Methods: From December 2016 to July 2018 (end of enrolment), 164 pts were enrolled. We evaluated by targeted sequencing the mutational profile of a panel of 34 genes: ACTB, ACTG1, ACTG2, ACTR1A, BAP1, CDH8, CDK4, CDKN2A, CDKN2B, COL3A1, COL5A2, CUL1, DHFR, GOT1, KDR, KIT, MXRA5, NF2, NFRKB, NKX6-2, NOD2, PCBD2, PDZK1IP1, PIK3CA, PIK3CB, PSMD13, RAPGEF6, RDX, SETDB1, TAOK1, TP53, TXNRD1, UQCRC1, XRCC6. Genetic profiling was correlated with clinical and pathological variables. Results: Overall, 110 pts (67%) from both treatment arms had samples available for molecular analysis. Median age was 63 years (45–81), 25.5% (n = 28) were females, and 74.5% (n = 82) were males. Tumor histotype was 81.8% (n = 90) epithelioid and 18.2% (n = 20) non-epithelioid; 28.5% of the tumors (n = 42) were stage IV, 71.5% (n = 68) were stage III. Targeted sequencing of tissue specimens identified 275 functional somatic mutations in the 34 genes analyzed. The number of mutated genes was positively associated with higher stage and metastatic disease (p = 0.025). RDX (42%), MXRA5 (23%), BAP1 (14%), and NF2 (11%) were the most frequently mutated genes. Mutations in RAPGEF6 (p = 0.03) and ACTG1 (p = 0.02) were associated with the non-epithelioid subtype, and mutations in BAP1 (p = 0.04) were related to progression-free survival (PFS) > 6 months. Conclusions: In the Ramucirumab Mesothelioma clinical trial (RAMES), mutation of the gene BAP1 is related to a prolonged PFS for patients treated with platinum/pemetrexed regimens (p = 0.04).

The genetic landscape of the human solute carrier (SLC) transporter superfamily

The human solute carrier (SLC) superfamily of transporters is comprised of over 400 membrane-bound proteins, and plays essential roles in a multitude of physiological and pharmacological processes. In addition, perturbation of SLC transporter function underlies numerous human diseases, which renders SLC transporters attractive drug targets. Common genetic polymorphisms in SLC genes have been associated with inter-individual differences in drug efficacy and toxicity. However, despite their tremendous clinical relevance, epidemiological data of these variants are mostly derived from heterogeneous cohorts of small sample size and the genetic SLC landscape beyond these common variants has not been comprehensively assessed. In this study, we analyzed Next-Generation Sequencing data from 141,456 individuals from seven major human populations to evaluate genetic variability, its functional consequences, and ethnogeographic patterns across the entire SLC superfamily of transporters. Importantly, of the 204,287 exonic single-nucleotide variants (SNVs) which we identified, 99.8% were present in less than 1% of analyzed alleles. Comprehensive computational analyses using 13 partially orthogonal algorithms that predict the functional impact of genetic variations based on sequence information, evolutionary conservation, structural considerations, and functional genomics data revealed that each individual genome harbors 29.7 variants with putative functional effects, of which rare variants account for 18%. Inter-ethnic variability was found to be extensive, and 83% of deleterious SLC variants were only identified in a single population. Interestingly, population-specific carrier frequencies of loss-of-function variants in SLC genes associated with recessive Mendelian disease recapitulated the ethnogeographic variation of the corresponding disorders, including cystinuria in Jewish individuals, type II citrullinemia in East Asians, and lysinuric protein intolerance in Finns, thus providing a powerful resource for clinical geneticists to inform about population-specific prevalence and allelic composition of Mendelian SLC diseases. In summary, we present the most comprehensive data set of SLC variability published to date, which can provide insights into inter-individual differences in SLC transporter function and guide the optimization of population-specific genotyping strategies in the bourgeoning fields of personalized medicine and precision public health.

Discovery of Novel Biomarkers of Therapeutic Responses in Han Chinese Pemetrexed-Based Treated Advanced NSCLC Patients

Pemetrexed, one of the most commonly used drugs in advanced non–small cell lung cancer (NSCLC) therapies, often leads to various therapeutic responses in patients. These therapeutic responses to pemetrexed, including adverse drug reactions (ADRs) and its intended therapeutic effects, have been demonstrated to be highly individual-specific. Such difference in therapeutic responses across individuals may be caused by the unique genetic variations in each patient. However, only a few pemetrexed-based studies have been performed using Han Chinese patients. In this study, we aimed to identify genetic signatures of therapeutic responses of pemetrexed-based treatment using 203 Han Chinese patients with advanced NSCLC. All the participants received two different types of therapies: 1) treatment with only pemetrexed and 2) treatment with both pemetrexed and platinum (mainly cisplatin and carboplatin). We then performed a genetic association analysis on 16 selected single-nucleotide polymorphisms (SNPs) in 7 genes using these 2 groups. The analysis of patients receiving only pemetrexed suggests that the SNP rs1051298 on the SLC19A1 gene (c.*746C > T) increased the risk of all ADRs (collected all types of ADRs) in different cycles of pemetrexed therapy [1-2 cycles: P = 0.0059, odds ratio (OR) = 3.143; 1-4 cycles: P = 0.0072, OR = 2.340; 1-6 cycles: P = 0.0071, OR = 2.243]. This influence of rs1051298 is particularly significant in terms of liver injury (1-4 cycles: P = 0.0056, OR = 3.863; 1-6 cycles: P = 0.0071, OR = 3.466). In all the patients, including patients who received both pemetrexed and platinum, SNP rs1801133 on the MTHFR gene (665C > T) was found to be significantly associated with hematological ADRs in 1 to 2 cycles (P = 0.0079, OR = 3.566). Additionally, we discovered that SNP rs12995526 (c.815-102T > C) in the ATIC gene and SNP rs11545077 (c.91G > T) in the GGH gene were associated with both ADRs and therapeutic effects. In summary, our study identified several potential biomarkers that were significantly associated with ADRs and therapeutic effects of pemetrexed-related treatments using Han Chinese patients. Our discoveries will provide important clues for personalized pemetrexed-based treatment design for Han Chinese NSCLC patients in the future.

Mutations in folate transporter genes and risk for human myelomeningocele

The molecular mechanisms linking folate deficiency and neural tube defect (NTD) risk in offspring remain unclear. Folate transporters (SLC19A1, SLC46A1, SLC25A32, and FOLH1) and folate receptors (FOLR1, FOLR2, and FOLR3) are suggested to play essential roles in transporting folate from maternal intestinal lumen to the developing embryo. Loss of function variants in these genes may affect folate availability and contribute to NTD risk. This study examines whether variants within the folate transporter and receptor genes are associated with an increased risk for myelomeningocele (MM). Exons and their flanking intron sequences of 348 MM subjects were sequenced using the Sanger sequencing method and/or next generation sequencing to identify variants. Frequencies of alleles of single nucleotide polymorphisms (SNPs) in MM subjects were compared to those from ethnically matched reference populations to evaluate alleles' associated risk for MM. We identified eight novel variants in SLC19A1 and twelve novel variants in FOLR1, FOLR2, and FOLR3. Pathogenic variants include c.1265delG in SLC19A1 resulting in an early stop codon, four large insertion deletion variants in FOLR3, and a stop_gain variant in FOLR3. No new variants were identified in SLC46A1, SLC25A32, or FOLH1. In SLC19A1, c.80A>G (rs1051266) was not associated with our MM cohort; we did observe a variant allele G frequency of 61.7%, higher than previously reported in other NTD populations. In conclusion, we discovered novel loss of function variants in genes involved in folate transport in MM subjects. Our results support the growing evidence of associations between genes involved in folate transport and susceptibility to NTDs.

An angiotensin-converting enzyme (ACE) polymorphism may mitigate the effects of angiotensin-pathway medications on posttraumatic stress symptoms

Angiotensin, which regulates blood pressure may also act within the brain to mediate stress and fear responses. Common antihypertensive medication classes of angiotensin-converting enzyme inhibitors (ACE-Is) and angiotensin receptor blockers (ARBs) have been associated with lower PTSD symptoms. Here we examine the rs4311 SNP in the ACE gene, previously implicated in panic attacks, in the relationship between ACE-I/ARB medications and PTSD symptoms. Participants were recruited from outpatient wait rooms between 2006 and March 2014 (n=  803). We examined the interaction between rs4311 genotype and the presence of blood pressure medication on PTSD symptoms and diagnosis. PTSD symptoms were lower in individuals taking ACE-Is or ARBs (N = 776). The rs4311 was associated with PTSD symptoms and diagnosis (N = 3803), as the T-carriers at the rs4311 SNP had significantly greater likelihood of a PTSD diagnosis. Lastly, the rs4311 genotype modified the effect of ACE-Is or ARBs on PTSD symptoms (N = 443; F1,443 = 4.41, P < 0.05). Individuals with the CC rs4311 genotype showed lower PTSD symptoms in the presence of ACE-Is or ARBs. In contrast, T- carriers showed the opposite, such that the presence of ACE-Is or ARBs was associated with higher PTSD symptoms. These data suggest that the renin-angiotensin system may be important in PTSD, as ACE-I/ARB usage associates with lower symptoms. Furthermore, we provide genetic evidence that some individuals are comparatively more benefitted by ACE-Is/ARBs in PTSD treatment. Future research should examine the mechanisms by which ACE-Is/ARBs affect PTSD symptoms such that pharmaco-genetically informed interventions may be used to treat PTSD.