Considerations on diagnosis and surveillance measures of PTEN hamartoma tumor syndrome: clinical and genetic study in a series of Spanish patients

Background

The limited knowledge about the PTEN hamartoma tumor syndrome (PHTS) makes its diagnosis a challenging task. We aimed to define the clinical and genetic characteristics of this syndrome in the Spanish population and to identify new genes potentially associated with the disease.

Results

We reviewed the clinical data collected through a specific questionnaire in a series of 145 Spanish patients with a phenotypic features compatible with PHTS and performed molecular characterization through several approaches including next generation sequencing and whole exome sequencing (WES). Macrocephaly, mucocutaneous lesions, gastrointestinal polyposis and obesity are prevalent phenotypic features in PHTS and help predict the presence of a PTEN germline variant in our population. We also find that PHTS patients are at risk to develop cancer in childhood or adolescence. Furthermore, we observe a high frequency of variants in exon 1 of PTEN, which are associated with renal cancer and overexpression of KLLN and PTEN. Moreover, WES revealed variants in genes like NEDD4 that merit further research.

Conclusions

This study expands previously reported findings in other PHTS population studies and makes new contributions regarding clinical and molecular aspects of PHTS, which are useful for translation to the clinic and for new research lines.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-021-02079-7.

MicroRNAs Targeting HIF-2α, VEGFR1 and/or VEGFR2 as Potential Predictive Biomarkers for VEGFR Tyrosine Kinase and HIF-2α Inhibitors in Metastatic Clear-Cell Renal Cell Carcinoma

Metastatic clear-cell renal cell carcinoma (m-ccRCC) is characterized by increased hypoxia-induced factor (HIF)-2α and vascular endothelial growth factor receptor (VEGFR)-dependent angiogenesis through loss of function of the von Hippel-Lindau protein. VEGFR tyrosine kinase inhibitors (VEGFR-TKIs) are a cornerstone of m-ccRCC treatment, and new treatments targeting HIF-2α are currently under investigation. However, predictive biomarkers for these treatments are lacking. In this retrospective cohort study including 109 patients treated with VEGFR-targeted therapies as first-line treatment, we aimed to study the possible predictive function of microRNAs (miRNAs) targeting HIF-2α, VEGFR1 and VEGFR2. We selected miRNAs inversely correlated with HIF-2α, VEGFR1 and/or VEGFR2 expression and with predicted target sites in the respective genes and subsequently studied their impact on therapeutic outcomes. We identified four miRNAs (miR-34c-5p, miR-221-3p, miR-222-3p and miR-3529-3p) inversely correlated with VEGFR1 and/or VEGFR2 expression and associated with tumor shrinkage and progression-free survival (PFS) upon treatment with VEGFR-TKIs, highlighting the potential predictive value of these miRNAs. Moreover, we identified three miRNAs (miR-185-5p, miR-223-3p and miR-3529-3p) inversely correlated with HIF-2α expression and associated with tumor shrinkage and PFS upon treatment with VEGFR-TKIs. These three miRNAs can have a predictive value not only upon treatment with VEGFR-TKIs but possibly also upon treatment with the upcoming HIF-2α inhibitor belzutifan.

MicroRNAs Possibly Involved in the Development of Bone Metastasis in Clear-Cell Renal Cell Carcinoma

Simple Summary

Bone metastases cause substantial morbidity and implicate worse clinical outcomes for clear-cell renal cell carcinoma patients. MicroRNAs are small RNA molecules that modulate gene translation and are involved in the development of cancer and metastasis. We identified six microRNAs that are potentially specifically involved in metastasis to bone, of which two seem protective and four implicate a higher risk. This aids further understanding of the process of metastasizing to bone. Furthermore, these microRNA hold potential for biomarkers or therapeutic targets.

Abstract

Bone metastasis in clear-cell renal cell carcinoma (ccRCC) leads to substantial morbidity through skeletal related adverse events and implicates worse clinical outcomes. MicroRNAs (miRNA) are small non-protein coding RNA molecules with important regulatory functions in cancer development and metastasis. In this retrospective analysis we present dysregulated miRNA in ccRCC, which are associated with bone metastasis. In particular, miR-23a-3p, miR-27a-3p, miR-20a-5p, and miR-335-3p specifically correlated with the earlier appearance of bone metastasis, compared to metastasis in other organs. In contrast, miR-30b-3p and miR-139-3p were correlated with less occurrence of bone metastasis. These miRNAs are potential biomarkers and attractive targets for miRNA inhibitors or mimics, which could lead to novel therapeutic possibilities for bone targeted treatment in metastatic ccRCC.

Integrative multi-omics analysis identifies a prognostic miRNA signature and a targetable miR-21-3p/TSC2/mTOR axis in metastatic pheochromocytoma/paraganglioma

Rationale: Pheochromocytomas and paragangliomas (PPGLs) are rare neuroendocrine tumors that present variable outcomes. To date, no effective therapies or reliable prognostic markers are available for patients who develop metastatic PPGL (mPPGL). Our aim was to discover robust prognostic markers validated through in vitro models, and define specific therapeutic options according to tumor genomic features.

Methods: We analyzed three PPGL miRNome datasets (n=443), validated candidate markers and assessed them in serum samples (n=36) to find a metastatic miRNA signature. An integrative study of miRNome, transcriptome and proteome was performed to find miRNA targets, which were further characterized in vitro.

Results: A signature of six miRNAs (miR-21-3p, miR-183-5p, miR-182-5p, miR-96-5p, miR-551b-3p, and miR-202-5p) was associated with metastatic risk and time to progression. A higher expression of five of these miRNAs was also detected in PPGL patients' liquid biopsies compared with controls. The combined expression of miR-21-3p/miR-183-5p showed the best power to predict metastasis (AUC=0.804, P=4.67·10^-18), and was found associated in vitro with pro-metastatic features, such as neuroendocrine-mesenchymal transition phenotype, and increased cell migration rate. A pan-cancer multi-omic integrative study correlated miR-21-3p levels with TSC2 expression, mTOR pathway activation, and a predictive signature for mTOR inhibitor-sensitivity in PPGLs and other cancers. Likewise, we demonstrated in vitro a TSC2 repression and an enhanced rapamycin sensitivity upon miR-21-3p expression.

Conclusions: Our findings support the assessment of miR-21-3p/miR-183-5p, in tumors and liquid biopsies, as biomarkers for risk stratification to improve the PPGL patients' management. We propose miR-21-3p to select mPPGL patients who may benefit from mTOR inhibitors.

Concomitant Medications and Risk of Chemotherapy-Induced Peripheral Neuropathy

BACKGROUND

Peripheral neuropathy is the dose-limiting toxicity of many oncology drugs, including paclitaxel. There is large interindividual variability in the neuropathy, and several risk factors have been proposed; however, many have not been replicated. Here we present a comprehensive study aimed at identifying treatment and physiopathology-related paclitaxel-induced neuropathy risk factors in a large cohort of well-characterized patients.

PATIENTS AND METHODS

Analyses included 503 patients with breast or ovarian cancer who received paclitaxel treatment. Paclitaxel dose modifications caused by the neuropathy were extracted from medical records and patients self-reported neuropathy symptoms were collected. Multivariate logistic regression analyses were performed to identify concomitant medications and comorbidities associated with paclitaxel-induced neuropathy.

RESULTS

Older patients had higher neuropathy: for each increase of 1 year of age, the risk of dose modifications and grade 3 neuropathy increased 4% and 5%, respectively. Cardiovascular drugs increased the risk of paclitaxel dose reductions (odds ratio [OR], 2.51; p = .006), with a stronger association for beta-adrenergic antagonists. The total number of concomitant medications also showed an association with dose modifications (OR, 1.25; p = .012 for each concomitant drug increase). A dose modification predictive model that included the new identified factors gave an area under the curve of 0.74 (p = 1.07 × 10^-10). Preexisting nerve compression syndromes seemed to increase neuropathy risk.

CONCLUSION

Baseline characteristics of the patients, including age and concomitant medications, could be used to identify individuals at high risk of neuropathy, personalizing chemotherapy treatment and reducing the risk of severe neuropathy.

IMPLICATIONS FOR PRACTICE

Peripheral neuropathy is a common adverse effect of many cancer drugs, including chemotherapeutics, targeted therapies, and immune checkpoint inhibitors. About 40% of survivors of cancer have functional deficits caused by this toxicity, some of them irreversible. Currently, there are no effective treatments to prevent or treat this neuropathy. This study, performed in a large cohort of well-characterized patients homogenously treated with paclitaxel, identified concomitant medications, comorbidities, and demographic factors associated with peripheral neuropathy. These factors could serve to identify patients at high risk of severe neuropathy for whom alternative non-neurotoxic alternatives may be considered.

PheoSeq: A Targeted Next-Generation Sequencing Assay for Pheochromocytoma and Paraganglioma Diagnostics

Genetic diagnosis is recommended for all pheochromocytoma and paraganglioma (PPGL) cases, as driver mutations are identified in approximately 80% of the cases. As the list of related genes expands, genetic diagnosis becomes more time-consuming, and targeted next-generation sequencing (NGS) has emerged as a cost-effective tool. This study aimed to optimize targeted NGS in PPGL genetic diagnostics. A workflow based on two customized targeted NGS assays was validated to study the 18 main PPGL genes in germline and frozen tumor DNA, with one of them specifically directed toward formalin-fixed paraffin-embedded tissue. The series involved 453 unrelated PPGL patients, of whom 30 had known mutations and were used as controls. Partial screening using Sanger had been performed in 275 patients. NGS results were complemented with the study of gross deletions. NGS assay showed a sensitivity ≥99.4%, regardless of DNA source. We identified 45 variants of unknown significance and 89 pathogenic mutations, the latter being germline in 29 (7.2%) and somatic in 58 (31.7%) of the 183 tumors studied. In 37 patients previously studied by Sanger sequencing, the causal mutation could be identified. We demonstrated that both assays are an efficient and accurate alternative to conventional sequencing. Their application facilitates the study of minor PPGL genes, and enables genetic diagnoses in patients with incongruent or missing clinical data, who would otherwise be missed.

Targeted Sequencing Reveals Low-Frequency Variants in EPHA Genes as Markers of Paclitaxel-Induced Peripheral Neuropathy

Purpose: Neuropathy is the dose-limiting toxicity of paclitaxel and a major cause for decreased quality of life. Genetic factors have been shown to contribute to paclitaxel neuropathy susceptibility; however, the major causes for interindividual differences remain unexplained. In this study, we identified genetic markers associated with paclitaxel-induced neuropathy through massive sequencing of candidate genes.Experimental Design: We sequenced the coding region of 4 EPHA genes, 5 genes involved in paclitaxel pharmacokinetics, and 30 Charcot-Marie-Tooth genes, in 228 cancer patients with no/low neuropathy or high-grade neuropathy during paclitaxel treatment. An independent validation series included 202 paclitaxel-treated patients. Variation-/gene-based analyses were used to compare variant frequencies among neuropathy groups, and Cox regression models were used to analyze neuropathy along treatment.Results: Gene-based analysis identified EPHA6 as the gene most significantly associated with paclitaxel-induced neuropathy. Low-frequency nonsynonymous variants in EPHA6 were present exclusively in patients with high neuropathy, and all affected the ligand-binding domain of the protein. Accumulated dose analysis in the discovery series showed a significantly higher neuropathy risk for EPHA5/6/8 low-frequency nonsynonymous variant carriers [HR, 14.60; 95% confidence interval (CI), 2.33-91.62; P = 0.0042], and an independent cohort confirmed an increased neuropathy risk (HR, 2.07; 95% CI, 1.14-3.77; P = 0.017). Combining the series gave an estimated 2.5-fold higher risk of neuropathy (95% CI, 1.46-4.31; P = 9.1 × 10^-4).Conclusions: This first study sequencing EPHA genes revealed that low-frequency variants in EPHA6, EPHA5, and EPHA8 contribute to the susceptibility to paclitaxel-induced neuropathy. Furthermore, EPHA's neuronal injury repair function suggests that these genes might constitute important neuropathy markers for many neurotoxic drugs. Clin Cancer Res; 23(5); 1227-35. ©2016 AACR.

Deep sequencing reveals microRNAs predictive of antiangiogenic drug response

The majority of metastatic renal cell carcinoma (RCC) patients are treated with tyrosine kinase inhibitors (TKI) in first-line treatment; however, a fraction are refractory to these antiangiogenic drugs. MicroRNAs (miRNAs) are regulatory molecules proven to be accurate biomarkers in cancer. Here, we identified miRNAs predictive of progressive disease under TKI treatment through deep sequencing of 74 metastatic clear cell RCC cases uniformly treated with these drugs. Twenty-nine miRNAs were differentially expressed in the tumors of patients who progressed under TKI therapy (P values from 6 × 10^-9 to 3 × 10^-3). Among 6 miRNAs selected for validation in an independent series, the most relevant associations corresponded to miR-1307-3p, miR-155-5p, and miR-221-3p (P = 4.6 × 10^-3, 6.5 × 10^-3, and 3.4 × 10^-2, respectively). Furthermore, a 2 miRNA-based classifier discriminated individuals with progressive disease upon TKI treatment (AUC = 0.75, 95% CI, 0.64-0.85; P = 1.3 × 10^-4) with better predictive value than clinicopathological risk factors commonly used. We also identified miRNAs significantly associated with progression-free survival and overall survival (P = 6.8 × 10^-8 and 7.8 × 10^-7 for top hits, respectively), and 7 overlapped with early progressive disease. In conclusion, this is the first miRNome comprehensive study, to our knowledge, that demonstrates a predictive value of miRNAs for TKI response and provides a new set of relevant markers that can help rationalize metastatic RCC treatment.

Tumor microRNA expression profiling identifies circulating microRNAs for early breast cancer detection

BACKGROUND

The identification of novel biomarkers for early breast cancer detection would be a great advance. Because of their role in tumorigenesis and stability in body fluids, microRNAs (miRNAs) are emerging as a promising diagnostic tool. Our aim was to identify miRNAs deregulated in breast tumors and evaluate the potential of circulating miRNAs in breast cancer detection.

METHODS

We conducted miRNA expression profiling of 1919 human miRNAs in paraffin-embedded tissue from 122 breast tumors and 11 healthy breast tissue samples. Differential expression analysis was performed, and a microarray classifier was generated. The most relevant miRNAs were analyzed in plasma from 26 healthy individuals and 83 patients with breast cancer (36 before and 47 after treatment) and validated in 116 healthy individuals and 114 patients before treatment.

RESULTS

We identified a large number of miRNAs deregulated in breast cancer and generated a 25-miRNA microarray classifier that discriminated breast tumors with high diagnostic sensitivity and specificity. Ten miRNAs were selected for further investigation, of which 4 (miR-505-5p, miR-125b-5p, miR-21-5p, and miR-96-5p) were significantly overexpressed in pretreated patients with breast cancer compared with healthy individuals in 2 different series of plasma. MiR-505-5p and miR-96-5p were the most valuable biomarkers (area under the curve 0.72). Moreover, the expression levels of miR-3656, miR-505-5p, and miR-21-5p were decreased in a group of treated patients.

CONCLUSIONS

Circulating miRNAs reflect the presence of breast tumors. The identification of deregulated miRNAs in plasma of patients with breast cancer supports the use of circulating miRNAs as a method for early breast cancer detection.

Thyroid cancer GWAS identifies 10q26.12 and 6q14.1 as novel susceptibility loci and reveals genetic heterogeneity among populations

Thyroid cancer is the most heritable cancer of all those not displaying typical Mendelian inheritance. However, most of the genetic factors that would explain the high heritability remain unknown. Our aim was to identify additional common genetic variants associated with susceptibility to this disease. In order to do so, we performed a genome-wide association study in a series of 398 cases and 502 controls from Spain, followed by a replication in four well-defined Southern European case-control collections contributing a total of 1,422 cases and 1,908 controls. The association between the variation at the 9q22 locus near FOXE1 and thyroid cancer risk was consistent across all series, with several SNPs identified (rs7028661: OR = 1.64, p = 1.0 × 10(-22) , rs7037324: OR = 1.54, p = 1.2 × 10(-17) ). Moreover, the rare alleles of three SNPs (rs2997312, rs10788123 and rs1254167) at 10q26.12 showed suggestive evidence of association with higher risk of the disease (OR = 1.35, p = 1.2 × 10(-04) , OR = 1.26, p = 5.2 × 10(-04) and OR = 1.38, p = 5.9 × 10(-05) , respectively). Finally, the rare allele of rs4075570 at 6q14.1 conferred protection in the series studied (OR = 0.82, p = 2.0 × 10(-04) ). This study suggests that heterogeneity in genetic susceptibility between populations is a key feature to take into account when exploring genetic risk factors related to this disease.

Whole-exome sequencing identifies MDH2 as a new familial paraganglioma gene

Disruption of the Krebs cycle is a hallmark of cancer. IDH1 and IDH2 mutations are found in many neoplasms, and germline alterations in SDH genes and FH predispose to pheochromocytoma/paraganglioma and other cancers. We describe a paraganglioma family carrying a germline mutation in MDH2, which encodes a Krebs cycle enzyme. Whole-exome sequencing was applied to tumor DNA obtained from a man age 55 years diagnosed with multiple malignant paragangliomas. Data were analyzed with the two-sided Student's t and Mann-Whitney U tests with Bonferroni correction for multiple comparisons. Between six- and 14-fold lower levels of MDH2 expression were observed in MDH2-mutated tumors compared with control patients. Knockdown (KD) of MDH2 in HeLa cells by shRNA triggered the accumulation of both malate (mean ± SD: wild-type [WT] = 1±0.18; KD = 2.24±0.17, P = .043) and fumarate (WT = 1±0.06; KD = 2.6±0.25, P = .033), which was reversed by transient introduction of WT MDH2 cDNA. Segregation of the mutation with disease and absence of MDH2 in mutated tumors revealed MDH2 as a novel pheochromocytoma/paraganglioma susceptibility gene.

Impact of chemotherapy on telomere length in sporadic and familial breast cancer patients

Recently, we observed that telomeres of BRCA1/2 mutation carriers were shorter than those of controls or sporadic breast cancer patients, suggesting that mutations in these genes might be responsible for this event. Given the contradictory results reported in the literature, we tested whether other parameters, such as chemotherapy, could be modifying telomere length (TL). We performed a cross-sectional study measuring leukocyte TL of 266 sporadic breasts cancer patients treated with first-line chemotherapy, with a median follow-up of 240 days. Additionally, we performed both cross-sectional and longitudinal studies in a series of 236 familial breast cancer patients that included affected and non-affected BRCA1/2 mutation carriers. We have measured in leukocytes from peripheral blood: the TL, percentage of short telomeres (<3 kb), telomerase activity levels and the annual telomere shortening speed. In sporadic cases we found that chemotherapy exerts a transient telomere shortening effect (around 2 years) that varies depending on the drug combination. In familial cases, only patients receiving treatment were associated with telomere shortening but they recovered normal TL after a period of 2 years. Chemotherapy affects TL and should be considered in the studies that correlate TL with disease susceptibility.

VEGF, VEGFR3, and PDGFRB protein expression is influenced by RAS mutations in medullary thyroid carcinoma

BACKGROUND

Tyrosine kinase inhibitors (TKIs) have achieved remarkable clinical results in medullary thyroid carcinoma (MTC) patients. However, the considerable variability in patient response to treatment with TKIs remains largely unexplained. There is evidence that it could be due, at least in part, to alterations in genes associated with the disease via their effect on the expression of TKI targets. The objective of this study was to evaluate the influence of RAS mutations on the expression levels in MTC tumors of eight key TKI target proteins.

METHODS

We assessed by immunohistochemistry the expression of EGFR, KIT, MET, PDGFRB, VEGF, VEGFR1, VEGFR2, and VEGFR3 in a series of 84 primary MTC tumors that had previously been molecularly characterized, including 14 RAS-positive, 18 RET(M918T)-positive, and 24 RET(C634)-positive tumors, as well as 15 wild-type tumors with no mutations in the RET or RAS genes.

RESULTS

In contrast to RET-positive tumors, RAS-positive tumors expressed neither PDGFRB nor MET (p=0.0060 and 0.047, respectively). Similarly, fewer RAS-positive than RET-related tumors expressed VEGFR3 (p=0.00062). Finally, wild-type tumors expressed VEGF more often than both RAS- and RET-positive tumors (p=0.0082 and 0.011, respectively).

CONCLUSIONS

This is the first study identifying that the expression of TKI targets differs according to the presence of RAS mutations in MTC. This information could potentially be used to select the most beneficial TKI treatment for these patients.

DNA methylation profiling of well-differentiated thyroid cancer uncovers markers of recurrence free survival

Thyroid cancer is a heterogeneous disease with several subtypes characterized by cytological, histological and genetic alterations, but the involvement of epigenetics is not well understood. Here, we investigated the role of aberrant DNA methylation in the development of well-differentiated thyroid tumors. We performed genome-wide DNA methylation profiling in the largest well-differentiated thyroid tumor series reported to date, comprising 83 primary tumors as well as 8 samples of adjacent normal tissue. The epigenetic profiles were closely related to not only tumor histology but also the underlying driver mutation; we found that follicular tumors had higher levels of methylation, which seemed to accumulate in a progressive manner along the tumorigenic process from adenomas to carcinomas. Furthermore, tumors harboring a BRAF or RAS mutation had a larger number of hypo- or hypermethylation events, respectively. The aberrant methylation of several candidate genes potentially related to thyroid carcinogenesis was validated in an independent series of 52 samples. Furthermore, through the integration of methylation and transcriptional expression data, we identified genes whose expression is associated with the methylation status of their promoters. Finally, by integrating clinical follow-up information with methylation levels we propose etoposide-induced 2.4 and Wilms tumor 1 as novel prognostic markers related to recurrence-free survival. This comprehensive study provides insights into the role of DNA methylation in well-differentiated thyroid cancer development and identifies novel markers associated with recurrence-free survival.

Integrative analysis of miRNA and mRNA expression profiles in pheochromocytoma and paraganglioma identifies genotype-specific markers and potentially regulated pathways

Pheochromocytomas (PCCs) and paragangliomas (PGLs) are rare neuroendocrine neoplasias of neural crest origin that can be part of several inherited syndromes. Although their mRNA profiles are known to depend on genetic background, a number of questions related to tumor biology and clinical behavior remain unanswered. As microRNAs (miRNAs) are key players in the modulation of gene expression, their comprehensive analysis could resolve some of these issues. Through characterization of miRNA profiles in 69 frozen tumors with germline mutations in the genes SDHD, SDHB, VHL, RET, NF1, TMEM127, and MAX, we identified miRNA signatures specific to, as well as common among, the genetic groups of PCCs/PGLs. miRNA expression profiles were validated in an independent series of 30 composed of VHL-, SDHB-, SDHD-, and RET-related formalin-fixed paraffin-embedded PCC/PGL samples using quantitative real-time PCR. Upregulation of miR-210 in VHL- and SDHB-related PCCs/PGLs was verified, while miR-137 and miR-382 were confirmed as generally upregulated in PCCs/PGLs (except in MAX-related tumors). Also, we confirmed overexpression of miR-133b as VHL-specific miRNAs, miR-488 and miR-885-5p as RET-specific miRNAs, and miR-183 and miR-96 as SDHB-specific miRNAs. To determine the potential roles miRNAs play in PCC/PGL pathogenesis, we performed bioinformatic integration and pathway analysis using matched mRNA profiling data that indicated a common enrichment of pathways associated with neuronal and neuroendocrine-like differentiation. We demonstrated that miR-183 and/or miR-96 impede NGF-induced differentiation in PC12 cells. Finally, global proteomic analysis in SDHB and MAX tumors allowed us to determine that miRNA regulation occurs primarily through mRNA degradation in PCCs/PGLs, which partially confirmed our miRNA-mRNA integration results.

Genome-wide association study identifies ephrin type A receptors implicated in paclitaxel induced peripheral sensory neuropathy

BACKGROUND

Peripheral neuropathy is the dose limiting toxicity of paclitaxel, a chemotherapeutic drug widely used to treat solid tumours. This toxicity exhibits great inter-individual variability of unknown origin. The present study aimed to identify genetic variants associated with paclitaxel induced neuropathy via a whole genome approach.

METHODS

A genome-wide association study (GWAS) was performed in 144 white European patients uniformly treated with paclitaxel/carboplatin and for whom detailed data on neuropathy was available. Per allele single nucleotide polymorphism (SNP) associations were assessed by Cox regression, modelling the cumulative dose of paclitaxel up to the development of grade 2 sensory neuropathy.

RESULTS

The strongest evidence of association was observed for the ephrin type A receptor 4 (EPHA4) locus (rs17348202, p=1.0×10(-6)), and EPHA6 and EPHA5 were among the top 25 and 50 hits (rs301927, p=3.4×10(-5) and rs1159057, p=6.8×10(-5)), respectively. A meta-analysis of EPHA5-rs7349683, the top marker for paclitaxel induced neuropathy in a previous GWAS (r(2)=0.79 with rs1159057), gave a hazard ratio (HR) estimate of 1.68 (p=1.4×10(-9)). Meta-analysis of the second hit of this GWAS, XKR4-rs4737264, gave a HR of 1.71 (p=3.1×10(-8)). Imputed SNPs at LIMK2 locus were also strongly associated with this toxicity (HR=2.78, p=2.0×10(-7)).

CONCLUSIONS

This study provides independent support of EPHA5-rs7349683 and XKR4-rs4737264 as the first markers of risk of paclitaxel induced neuropathy. In addition, it suggests that other EPHA genes also involved in axonal guidance and repair following neural injury, as well as LIMK2 locus, may play an important role in the development of this toxicity. The identified SNPs could form the basis for individualised paclitaxel chemotherapy.

Genetic variation in the SLC19A1 gene and methotrexate toxicity in rheumatoid arthritis patients

AIM

We investigated the clinical relevance of SLC19A1 genetic variability for methotrexate (MTX) toxicity in rheumatoid arthritis patients using a haplotype-based approach.

PATIENTS & METHODS

Two hundred and twelve unrelated rheumatoid arthritis patients and 89 lymphoblastoid cell lines were used to investigate the effect of SLC19A1 SNPs and haplotypes on MTX adverse events and treatment discontinuation.

RESULTS

Two putatively functional SNPs in high linkage disequilibrium, rs1051266 and rs1131596, were associated with protection (hazard ratio: 0.33; 95% CI: 0.16-0.69; adjusted p = 0.021 and hazard ratio: 0.38; 95% CI: 0.17-0.27; adjusted p = 0.021, respectively) of discontinuation of MTX treatment owing to toxicity. These SNPs were also associated with protection from infections. SLC19A1 haplotype analysis found significant associations with MTX discontinuation owing to toxicity (p = 0.025). Quantification of SLC19A1 mRNA in cell lines suggested that rs1131596 was not a major causal variant.

CONCLUSION

Individual SNP and haplotype analyses suggest that rs1051266 could be a functional variant altering MTX toxicity; however, validation in independent studies is needed.

Hematologic β-tubulin VI isoform exhibits genetic variability that influences paclitaxel toxicity

Cellular microtubules composed of α-β-tubulin heterodimers that are essential for cell shape, division, and intracellular transport are valid targets for anticancer therapy. However, not all the conserved but differentially expressed members of the β-tubulin gene superfamily have been investigated for their role in these settings. In this study, we examined roles for the hematologic isoform β-tubulin VI and functional genetic variants in the gene. β-tubulin VI was highly expressed in blood cells with a substantial interindividual variability (seven-fold variation in mRNA). We characterized DNA missense variations leading to Q43P, T274M, and R307H, and a rare nonsense variant, Y55X. Because variations in the hematologic target of microtubule-binding drugs might alter their myelosuppressive action, we tested their effect in cell lines stably expressing the different β-tubulin VI full-length variants, finding that the T274M change significantly decreased sensitivity to paclitaxel-induced tubulin polymerization. Furthermore, patients treated with paclitaxel and carrying β-tubulin VI T274M exhibited a significantly lower thrombocytopenia than wild-type homozygous patients (P = 0.031). Together, our findings define β-tubulin VI as a hematologic isotype with significant genetic variation in humans that may affect the myelosuppresive action of microtubule-binding drugs. A polymorphism found in a tubulin isoform expressed only in hemapoietic cells may contribute to the patient variation in myelosuppression that occurs after treatment with microtubule-binding drugs.

MAX mutations cause hereditary and sporadic pheochromocytoma and paraganglioma

PURPOSE

Pheochromocytomas (PCC) and paragangliomas (PGL) are genetically heterogeneous neural crest-derived neoplasms. Recently we identified germline mutations in a new tumor suppressor susceptibility gene, MAX (MYC-associated factor X), which predisposes carriers to PCC. How MAX mutations contribute to PCC/PGL and associated phenotypes remain unclear. This study aimed to examine the prevalence and associated phenotypic features of germline and somatic MAX mutations in PCC/PGL.

DESIGN

We sequenced MAX in 1,694 patients with PCC or PGL (without mutations in other major susceptibility genes) from 17 independent referral centers. We screened for large deletions/duplications in 1,535 patients using a multiplex PCR-based method. Somatic mutations were searched for in tumors from an additional 245 patients. The frequency and type of MAX mutation was assessed overall and by clinical characteristics.

RESULTS

Sixteen MAX pathogenic mutations were identified in 23 index patients. All had adrenal tumors, including 13 bilateral or multiple PCCs within the same gland (P < 0.001), 15.8% developed additional tumors at thoracoabdominal sites, and 37% had familial antecedents. Age at diagnosis was lower (P = 0.001) in MAX mutation carriers compared with nonmutated cases. Two patients (10.5%) developed metastatic disease. A mutation affecting MAX was found in five tumors, four of them confirmed as somatic (1.65%). MAX tumors were characterized by substantial increases in normetanephrine, associated with normal or minor increases in metanephrine.

CONCLUSIONS

Germline mutations in MAX are responsible for 1.12% of PCC/PGL in patients without evidence of other known mutations and should be considered in the genetic work-up of these patients.

The miR-200 family controls beta-tubulin III expression and is associated with paclitaxel-based treatment response and progression-free survival in ovarian cancer patients

Ovarian cancer remains one of the leading causes of cancer deaths. Thus, new biomarkers predictive of response to the standard paclitaxel-carboplatin treatment are needed to improve chemotherapy strategies. MicroRNAs have the potential to modify drug outcomes. Based on this, we have demonstrated in this study that patients with a high expression of the miR-200 family show low levels of β-tubulin class III in ovarian carcinoma. In addition, we have established the clinical relevance of these microRNAs for ovarian cancer patients' treatment response and survival. In a well-characterized series of 72 ovarian carcinomas, the expressions of miR-141, miR-200a, miR-200b, miR-200c, and miR-429 were quantified by quantitative reverse transcription-PCR, and the protein content of β-tubulin isotypes I, II, and III was determined by immunohistochemistry. The relationship between these microRNAs, β-tubulin expression, response to paclitaxel-based treatment, progression-free survival (PFS) and overall survival was determined. While isotype I had constant high levels, protein expression of β-tubulins II and III was mutually exclusive. Low tumoral miR-200 expression was significantly associated with high β-tubulin III protein content (P values range, 0.047-<0.0001), and patients without complete response (CR) had lower miR-200c levels than patients with CR (hazard ratio (HR)=1.43, 95% confidence interval (CI)=1.02-1.99, P=0.037, multivariate analysis). Additionally, low miR-200 family expression had a trend toward poor PFS (HR>2.0, P values 0.051, 0.054, and 0.079 for miR-200c, miR-141, and miR-429 respectively, multivariate analysis). In conclusion, miR-200 family members affect the final β-tubulin III protein content of ovarian carcinomas. Furthermore, these microRNAs might constitute the biomarkers of response to paclitaxel-based treatments and relapse/progression of advanced stage ovarian carcinoma patients.

Research resource: Transcriptional profiling reveals different pseudohypoxic signatures in SDHB and VHL-related pheochromocytomas

The six major genes involved in hereditary susceptibility for pheochromocytoma (PCC)/paraganglioma (PGL) (RET, VHL, NF1, SDHB, SDHC, and SDHD) have been recently integrated into the same neuronal apoptotic pathway where mutations in any of these genes lead to cell death. In this model, prolyl hydroxylase 3 (EglN3) abrogation plays a pivotal role, but the molecular mechanisms underlying its inactivation are currently unknown. The aim of the study was to decipher specific alterations associated with the different genetic classes of PCCs/PGLs. With this purpose, 84 genetically characterized tumors were analyzed by means of transcriptional profiling. The analysis revealed a hypoxia-inducible factor (HIF)-related signature common to succinate dehydrogenase (SDH) and von Hippel-Lindau (VHL) tumors, that differentiated them from RET and neurofibromatosis type 1 cases. Both canonical HIF-1α and HIF-2α target genes were overexpressed in the SDH/VHL cluster, suggesting that a global HIF deregulation accounts for this common profile. Nevertheless, when we compared VHL tumors with SDHB cases, which often exhibit a malignant behavior, we found that HIF-1α target genes showed a predominant activation in the VHL PCCs. Expression data from 67 HIF target genes was sufficient to cluster SDHB and VHL tumors into two different groups, demonstrating different pseudo-hypoxic signatures. In addition, VHL-mutated tumors showed an unexpected overexpression of EglN3 mRNA that did not lead to significantly different EglN3 protein levels. These findings pave the way for more specific therapeutic approaches for malignant PCCs/PGLs management based on the patient's genetic alteration.

Allelic variant at -79 (C>T) in CDKN1B (p27Kip1) confers an increased risk of thyroid cancer and alters mRNA levels

The aim of this study is to assess if common genetic variants located in the CDKN1B locus, coding for the cell cycle inhibitor p27(Kip1), are involved in thyroid cancer susceptibility. Based on the literature and functional predictions, we selected three polymorphisms within the CDKN1B gene (rs2066827 (T326G, V109G), rs34330 (-79C>T) and rs36228499 (-838C>A)) to perform the first case-control study in thyroid cancer involving this locus. We had 649 Spanish patients with sporadic thyroid cancer and 385 healthy representative controls available. Luciferase reporter gene assays, real-time quantitative reverse transcription-PCR and immunoblot experiments were carried out to demonstrate the putative effect of the associated variant. The polymorphism rs34330 (-79C>T) was identified as a risk factor for developing the follicular variant of papillary thyroid carcinoma (FVPTC), fitting a recessive model (odds ratio=2.12; 95% confidence interval=1.09-4.15; P value=0.023). The risk allele (T) of this single nucleotide polymorphism led to a lower transcription rate in cells transfected with a luciferase reporter driven by the polymorphic p27(Kip1) promoter (P value <0.001). This effect was observed in -79TT genotype control carriers, who showed a tendency towards lower CDKN1B mRNA levels in lymphocytes, as well as at the protein level. This is the first study that identifies CDKN1B as a low-penetrance gene in thyroid cancer, and specifically in FVPTC subtype. We propose a reduced CDKN1B gene transcription depending on the genotype of the -79C>T (rs34330) variant as a novel mechanism underlying p27(Kip1) downregulation.

The variant rs1867277 in FOXE1 gene confers thyroid cancer susceptibility through the recruitment of USF1/USF2 transcription factors

In order to identify genetic factors related to thyroid cancer susceptibility, we adopted a candidate gene approach. We studied tag- and putative functional SNPs in genes involved in thyroid cell differentiation and proliferation, and in genes found to be differentially expressed in thyroid carcinoma. A total of 768 SNPs in 97 genes were genotyped in a Spanish series of 615 cases and 525 controls, the former comprising the largest collection of patients with this pathology from a single population studied to date. SNPs in an LD block spanning the entire FOXE1 gene showed the strongest evidence of association with papillary thyroid carcinoma susceptibility. This association was validated in a second stage of the study that included an independent Italian series of 482 patients and 532 controls. The strongest association results were observed for rs1867277 (OR[per-allele] = 1.49; 95%CI = 1.30–1.70; P = 5.9×10⁻⁹). Functional assays of rs1867277 (NM_004473.3:c.−283G>A) within the FOXE1 5′ UTR suggested that this variant affects FOXE1 transcription. DNA-binding assays demonstrated that, exclusively, the sequence containing the A allele recruited the USF1/USF2 transcription factors, while both alleles formed a complex in which DREAM/CREB/αCREM participated. Transfection studies showed an allele-dependent transcriptional regulation of FOXE1. We propose a FOXE1 regulation model dependent on the rs1867277 genotype, indicating that this SNP is a causal variant in thyroid cancer susceptibility. Our results constitute the first functional explanation for an association identified by a GWAS and thereby elucidate a mechanism of thyroid cancer susceptibility. They also attest to the efficacy of candidate gene approaches in the GWAS era.

Although follicular cell-derived thyroid cancer has an important genetic component, efforts in identifying major susceptibility genes have not been successful. Probably this is due to the complex nature of this disease that involves both genetic and environmental factors, as well as the interaction between them, which could be ultimately modulating the individual susceptibility. In this study, focused on genes carefully selected by their biological relation with the disease, and using more than 1,000 cases and 1,000 representative controls from two independent Caucasian populations, we demonstrate that FOXE1 is associated with Papillary Thyroid Cancer susceptibility. Functional assays prove that rs1867277 behaves as a genetic causal variant that regulates FOXE1 expression through a complex transcription factor network. This approach constitutes a successful approximation to define thyroid cancer risk genes related to individual susceptibility, and identifies FOXE1 as a key factor for its development.