Napabucasin overcomes cisplatin resistance in ovarian germ cell tumor-derived cell line by inhibiting cancer stemness.

BACKGROUND: Cisplatin resistance of ovarian yolk sac tumors (oYST) is a clinical challenge due to dismal patient prognosis, even though the disease is extremely rare. We investigated potential association between cisplatin resistance and cancer stem cell (CSC) markers in chemoresistant oYST cells and targeting strategies to overcome resistance in oYST. METHODS: Chemoresistant cells were derived from chemosensitive human oYST cells by cultivation in cisplatin in vitro. Derivative cells were characterized by chemoresistance, functional assays, flow cytometry, gene expression and protein arrays focused on CSC markers. RNAseq, methylation and microRNA profiling were performed. Quail chorioallantoic membranes (CAM) with implanted oYST cells were used to analyze the micro-tumor extent and interconnection with the CAM. Tumorigenicity in vivo was determined on immunodeficient mouse model. Chemoresistant cells were treated by inhibitors intefering with the CSC properties to examine the chemosensitization to cisplatin. RESULTS: Long-term cisplatin exposure resulted in seven-fold higher IC50 value in resistant cells, cross-resistance to oxaliplatin and carboplatin, and increased migratory capacity, invasiveness and tumorigenicity, associated with hypomethylation of differentially methylated genes/promotors. Resistant cells exhibited increased expression of prominin-1 (CD133), ATP binding cassette subfamily G member 2 (ABCG2), aldehyde dehydrogenase 3 isoform A1 (ALDH3A1), correlating with reduced gene and promoter methylation, as well as increased expression of ALDH1A3 and higher overall ALDH enzymatic activity, rendering them cross-resistant to DEAB, disulfiram and napabucasin. Salinomycin and tunicamycin were significantly more toxic to resistant cells. Pretreatment with napabucasin resensitized the cells to cisplatin and reduced their tumorigenicity in vivo. CONCLUSIONS: The novel chemoresistant cells represent unique model of refractory oYST. CSC markers are associated with cisplatin resistance being possible targets in chemorefractory oYST.

Molecular heterogeneity and early metastatic clone selection in testicular germ cell cancer development.

BACKGROUND: Testicular germ cell cancer (TGCC), being the most frequent malignancy in young Caucasian males, is initiated from an embryonic germ cell. This study determines intratumour heterogeneity to unravel tumour progression from initiation until metastasis. METHODS: In total, 42 purified samples of four treatment-resistant nonseminomatous (NS) TGCC were investigated, including the precursor germ cell neoplasia in situ (GCNIS) and metastatic specimens, using whole-genome and targeted sequencing. Their evolution was reconstructed. RESULTS: Intratumour molecular heterogeneity did not correspond to the supposed primary tumour histological evolution. Metastases after systemic treatment could be derived from cancer stem cells not identified in the primary cancer. GCNIS mostly lacked the molecular marks of the primary NS and comprised dominant clones that failed to progress. A BRCA-like mutational signature was observed without evidence for direct involvement of BRCA1 and BRCA2 genes. CONCLUSIONS: Our data strongly support the hypothesis that NS is initiated by whole-genome duplication, followed by chromosome copy number alterations in the cancer stem cell population, and accumulation of low numbers of somatic mutations, even in therapy-resistant cases. These observations of heterogeneity at all stages of tumourigenesis should be considered when treating patients with GCNIS-only disease, or with clinically overt NS.

Imprints and DPPA3 are bypassed during pluripotency- and differentiation-coupled methylation reprogramming in testicular germ cell tumors.

Testicular germ cell tumors (TGCTs) share germline ancestry but diverge phenotypically and clinically as seminoma (SE) and nonseminoma (NSE), the latter including the pluripotent embryonal carcinoma (EC) and its differentiated derivatives, teratoma (TE), yolk sac tumor (YST), and choriocarcinoma. Epigenomes from TGCTs may illuminate reprogramming in both normal development and testicular tumorigenesis. Herein we investigate pure-histological forms of 130 TGCTs for conserved and subtype-specific DNA methylation, including analysis of relatedness to pluripotent stem cell (ESC, iPSC), primordial germ cell (PGC), and differentiated somatic references. Most generally, TGCTs conserve PGC-lineage erasure of maternal and paternal genomic imprints and DPPA3 (also known as STELLA); however, like ESCs, TGCTs show focal recurrent imprinted domain hypermethylation. In this setting of shared physiologic erasure, NSEs harbor a malignancy-associated hypermethylation core, akin to that of a diverse cancer compendium. Beyond these concordances, we found subtype epigenetic homology with pluripotent versus differentiated states. ECs demonstrate a striking convergence of both CpG and CpH (non-CpG) methylation with pluripotent states; the pluripotential methyl-CpH signature crosses species boundaries and is distinct from neuronal methyl-CpH. EC differentiation to TE and YST entails reprogramming toward the somatic state, with loss of methyl-CpH but de novo methylation of pluripotency loci such as NANOG Extreme methyl-depletion among SE reflects the PGC methylation nadir. Adjacent to TGCTs, benign testis methylation profiles are determined by spermatogenetic proficiency measured by Johnsen score. In sum, TGCTs share collective entrapment in a PGC-like state of genomic-imprint and DPPA3 erasure, recurrent hypermethylation of cancer-associated targets, and subtype-dependent pluripotent, germline, or somatic methylation.

DMRforPairs: identifying differentially methylated regions between unique samples using array based methylation profiles.

BACKGROUND: Array based methylation profiling is a cost-effective solution to study the association between genome methylation and human disease & development. Available tools to analyze the Illumina Infinium HumanMethylation450 BeadChip focus on comparing methylation levels per locus. Other tools combine multiple probes into a range, identifying differential methylated regions (DMRs). These tools all require groups of samples to compare. However, comparison of unique, individual samples is essential in situations where larger sample sizes are not possible. RESULTS: DMRforPairs was designed to compare regional methylation status between unique samples. It identifies probe dense genomic regions and quantifies/tests their (difference in) methylation level between the samples. As a proof of concept, DMRforPairs is applied to public data from four human cell lines: two lymphoblastoid cell lines from healthy individuals and the cancer cell lines A431 and MCF7 (including 2 technical replicates each). DMRforPairs identified an increasing number of DMRs related to the sample phenotype when biological similarity of the samples decreased. DMRs identified by DMRforPairs were related to the biological origin of the cell lines. CONCLUSION: To our knowledge, DMRforPairs is the first tool to identify and visualize relevant and significant differentially methylated regions between unique samples.

miMsg: a target enrichment algorithm for predicted miR-mRNA interactions based on relative ranking of matched expression data.

MOTIVATION: Algorithms predicting microRNA (miR)-mRNA interactions generate high numbers of possible interactions, many of which might be non-existent or irrelevant in a certain biological context. It is desirable to develop a transparent, user-friendly, unbiased tool to enrich miR-mRNA predictions. RESULTS: The miMsg algorithm uses matched miR/mRNA expression data to enrich miR-mRNA predictions. It grades interactions by the number, magnitude and significance of misplacements in the combined ranking profiles of miR/mRNA expression assessed over multiple biological samples. miMsg requires minimal user input and makes no statistical assumptions. It identified 921 out of 56 262 interactions as top scoring and significant in an actual germ cell cancer dataset. Twenty-eight miR-mRNA pairs were deemed of highest interest based on ranking by miMsg and supported by current knowledge about validated interactions and biological function. To conclude, miMsg is an effective algorithm to reduce a high number of predicted interactions to a small set of high confidence interactions for further study. AVAILABILITY AND IMPLEMENTATION: Matlab source code and datasets available at www.martinrijlaarsdam.nl/mimsg . SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Protein pathway activation mapping reveals molecular networks associated with antiestrogen resistance in breast cancer cell lines.

Previously, we have identified a panel of breast cancer antiestrogen resistance (BCAR) genes. Several of these genes have clinical relevance because mRNA or protein levels associate with tamoxifen resistance or tumor aggressiveness. We postulated that changes in activation status of protein signaling networks induced by BCAR genes may provide better insight into the mechanisms underlying antiestrogen resistance. Key signal transduction pathways were analyzed for changes in activation or expression using reverse-phase protein microarrays probed with 78 antibodies against signaling proteins with known roles in tumorigenesis. We used ZR-75-1-derived cell lines transduced with AKT1, AKT2, BCAR1, BCAR3, BCAR4, EGFR, GRB7, HRAS, HRAS(v12) or HEF1 and MCF7-derived cell lines transduced with BCAR3, BCAR4 or EGFR. In the antiestrogen-resistant cell lines, we observed increased phosphorylation of several pathways involved in cell proliferation and survival. All tamoxifen-resistant cell lines contained high levels of phosphorylated AKT and its biochemically linked substrates Forkhead box O1/3. The activation of ERBB2, ERBB3 and the downstream modulators focal adhesion kinase and SHC were activated in cells with overexpression of BCAR4. Remarkable differences were observed for the levels of activated AMPK alpha1, cyclins, STAT5, STAT6, ERK1/2 and BCL2. The comparison of the cell signaling networks in estrogen-dependent and -independent cell lines revealed biochemically linked kinase-substrate markers that comprised systemically activated signaling pathways involved in tamoxifen resistance. Our results show that this model provides insights into the molecular and cellular mechanisms of breast cancer progression and antiestrogen resistance. This knowledge may help the development of novel targeted treatments.

Chromosome 3p25.3 Gain Is Associated With Cisplatin Resistance and Is an Independent Predictor of Poor Outcome in Male Malignant Germ Cell Tumors.

PURPOSE: Cisplatin is the main systemic treatment modality for male type II germ cell tumors (GCTs). Although generally very effective, 5%-10% of patients suffer from cisplatin-resistant disease. Identification of the driving mechanisms of resistance will enable improved risk stratification and development of alternative treatments. METHODS: We developed and characterized cisplatin-resistant GCT cell line models and compared their molecular characteristics with patient samples with cisplatin resistance and/or a poor clinical outcome. Subsequently, the association between the overlapping genetic features and clinical data was assessed. Finally, we used Cox regression to determine the prognostic relevance of these features within the currently used risk classification. RESULTS: Gain of chromosome 3p25.3 was detected in all cisplatin-resistant cell lines, and copy number of this region correlated with the level of resistance (R = 0.96, P = 1.5e-04). Gain of this region was detected at low frequencies in primary tumors and at higher frequencies in relapsed and/or cisplatin-resistant tumors. Chromosome 3p25.3 gain was associated with shorter progression-free survival and overall survival, with the strongest association observed in nonseminomas excluding pure teratomas. 3p25.3 gain was more frequently observed in tumors with yolk sac tumor histology and predicted adverse outcome independent of the International Germ Cell Cancer Collaborative Group risk classification and the presence of TP53/MDM2 alterations. CONCLUSION: On the basis of both in vitro analyses and clinical data, we found 3p25.3 to be strongly associated with cisplatin resistance and poor clinical outcome in male type II GCTs. Using genomic profiling, 3p25.3 status could help to improve risk stratification in male patients with type II GCT. Further characterization of this locus and underlying mechanisms of resistance is warranted to guide development of novel treatment approaches for cisplatin-resistant disease.

Characterization of BCAR4, a novel oncogene causing endocrine resistance in human breast cancer cells.

Resistance to the antiestrogen tamoxifen remains a major problem in the management of estrogen receptor-positive breast cancer. Knowledge on the resistance mechanisms is needed to develop more effective therapies. Breast cancer antiestrogen resistance 4 (BCAR4) was identified in a functional screen for genes involved in tamoxifen resistance. BCAR4 is expressed in 27% of primary breast tumors. In patients treated with tamoxifen for metastized disease high BCAR4 mRNA levels are associated with reduced clinical benefit and progression-free survival. Regarding tumor aggressiveness high BCAR4 mRNA levels are associated with a shorter metastasis free survival and overall survival. In the present study, we investigated the role of BCAR4 in endocrine resistance. Forced expression of BCAR4 in human ZR-75-1 and MCF7 breast cancer cells resulted in cell proliferation in the absence of estrogen and in the presence of various antiestrogens. Inhibition of estrogen receptor 1 (ESR1) expression with small interfering RNA (siRNA), implied that the BCAR4-induced mechanism of resistance is independent of ESR1. Highly conserved BCAR4 homologues of rhesus monkey, green monkey, and the less conserved common marmoset gene induced tamoxifen-resistant cell proliferation, in contrast to the distant BCAR4 homologues of bovine and rabbit. Injection of BCAR4-expressing ZR-75-1 cells into nude mice resulted in rapidly growing tumors. In silico analysis showed that BCAR4 mRNA is highly expressed in human placenta and oocyte, and absent in other normal tissues. In conclusion, BCAR4 is a strong transforming gene causing estrogen-independent growth and antiestrogen resistance, and induces tumor formation in vivo. Due to its restricted expression, BCAR4 may be a good target for treating antiestrogen-resistant breast cancer.

Prevalence of the variant allele rs61764370 T>G in the 3'UTR of KRAS among Dutch BRCA1, BRCA2 and non-BRCA1/BRCA2 breast cancer families.

Recently, a variant allele in the 3'UTR of the KRAS gene (rs61764370 T>G) was shown to be associated with an increased risk for developing non-small cell lung cancer, as well as ovarian cancer, and was most enriched in ovarian cancer patients from hereditary breast and ovarian cancer families. This functional variant has been shown to disrupt a let-7 miRNA binding site leading to increased expression of KRAS in vitro. In the current study, we have genotyped this KRAS-variant in breast cancer index cases from 268 BRCA1 families, 89 BRCA2 families, 685 non-BRCA1/BRCA2 families, and 797 geographically matched controls. The allele frequency of the KRAS-variant was found to be increased among patients with breast cancer from BRCA1, but not BRCA2 or non-BRCA1/BRCA2 families as compared to controls. As BRCA1 carriers mostly develop ER-negative breast cancers, we also examined the variant allele frequency among indexes from non-BRCA1/BRCA2 families with ER-negative breast cancer. The prevalence of the KRAS-variant was, however, not significantly increased as compared to controls, suggesting that the variant allele not just simply associates with ER-negative breast cancer. Subsequent expansion of the number of BRCA1 carriers with breast cancer by including other family members in addition to the index cases resulted in loss of significance for the association between the variant allele and mutant BRCA1 breast cancer. In this same cohort, the KRAS-variant did not appear to modify breast cancer risk for BRCA1 carriers. Importantly, results from the current study suggest that KRAS-variant frequencies might be increased among BRCA1 carriers, but solid proof requires confirmation in a larger cohort of BRCA1 carriers.

The substrate domain of BCAR1 is essential for anti-estrogen-resistant proliferation of human breast cancer cells.

To unravel the mechanisms underlying failure of endocrine therapy of breast cancer, we have previously executed a functional genetic screen and identified the adaptor protein BCAR1 to be causative for tamoxifen resistance. As a consequence of the manifold of interactions with other proteins, we characterized the contribution of individual protein domains of BCAR1 to anti-estrogen-resistant proliferation of human breast cancer cells. We took advantage of the observation that the closely related family member HEF1 was unable to support long-term anti-estrogen-resistant cell proliferation. Chimerical proteins containing defined domains of BCAR1 and HEF1 were evaluated for anti-estrogen-resistant growth. Exchange of the SH3 and C-terminal domains did not modify the capacity to support cell proliferation. Full support of anti-estrogen resistant proliferation was observed for chimerical molecules containing the central part of BCAR1. The bi-partite SRC-binding site or the Serine-rich domain did not explain the differential capacity of BCAR1. These findings indicate that the differences between BCAR1 and HEF1 with respect to support of anti-estrogen resistance reside in the substrate domain which contains multiple sites for tyrosine phosphorylation. The crucial interactions required for anti-estrogen resistance occur within the substrate domain of BCAR1. Further deciphering of these interactions may resolve the growth regulatory mechanism and provide an explanation for the observation that primary tumors with high levels of BCAR1 are likely to fail on tamoxifen therapy. This information may also help to devise alternative personalized treatment strategies with improved outcome for breast cancer patients.

Comparing genome-scale DNA methylation and CNV marks between adult human cultured ITGA6+ testicular cells and seminomas to assess in vitro genomic stability.

Autologous transplantation of spermatogonial stem cells is a promising new avenue to restore fertility in infertile recipients. Expansion of the initial spermatogonial stem cell pool through cell culturing is a necessary step to obtain enough cells for effective repopulation of the testis after transplantation. Since in vitro propagation can lead to (epi-)genetic mutations and possibly malignant transformation of the starting cell population, we set out to investigate genome-wide DNA methylation status in uncultured and cultured primary testicular ITGA6+ sorted cells and compare them with germ cell tumor samples of the seminoma subtype. Seminomas displayed a severely global hypomethylated profile, including loss of genomic imprinting, which we did not detect in cultured primary testicular ITGA6+ cells. Differential methylation analysis revealed altered regulation of gamete formation and meiotic processes in cultured primary testicular ITGA6+ cells but not in seminomas. The pivotal POU5F1 marker was hypomethylated in seminomas but not in uncultured or cultured primary testicular ITGA6+ cells, which is reflected in the POU5F1 mRNA expression levels. Lastly, seminomas displayed a number of characteristic copy number variations that were not detectable in primary testicular ITGA6+ cells, either before or after culture. Together, the data show a distinct DNA methylation patterns in cultured primary testicular ITGA6+ cells that does not resemble the pattern found in seminomas, but also highlight the need for more sensitive methods to fully exclude the presence of malignant cells after culture and to further study the epigenetic events that take place during in vitro culture.

Functional identification of genes causing estrogen independence of human breast cancer cells.

Endocrine treatment of breast cancer is widely applied and effective. However, in advanced disease cases, the tumors will eventually progress into an estrogen-independent and therapy-resistant phenotype. To elucidate the molecular mechanisms underlying this endocrine therapy failure, we applied retroviral insertion mutagenesis to identify the main genes conferring estrogen independence to human breast cancer cells. Estrogen-dependent ZR-75-1 cells were infected with replication-defective retroviruses followed by selection with the anti-estrogen 4-hydroxy-tamoxifen. In the resulting panel of 79 tamoxifen-resistant cell lines, the viral integrations were mapped within the human genome. Genes located in the immediate proximity of the retroviral integration sites were characterized for altered expression and their capacity to confer anti-estrogen resistance when transfected into breast cancer cells. Out of 15 candidate BCAR (breast cancer anti-estrogen resistance) genes, seven (AKT1, AKT2, BCAR1, BCAR3, EGFR, GRB7, and TRERF1/BCAR2) were shown to directly underlie estrogen independence. Our results show that insertion mutagenesis is a powerful tool to identify BCAR loci, which may provide insights into the molecular and cellular mechanisms of breast tumor progression and therapy resistance thereby offering novel targets for the development of tailor-made therapeutical and prevention strategies.

TSC22D1 and PSAP predict clinical outcome of tamoxifen treatment in patients with recurrent breast cancer.

Purpose Two genes, TSC22 domain family, member 1 (TSC22D1) and prosaposin (PSAP) were identified in an in vitro functional screen for genes having a causative role in tamoxifen resistance. These genes were also present in our previously established 81-gene signature for resistance to first-line tamoxifen therapy. The aim of this study was to investigate the predictive value of these genes for tamoxifen therapy failure in patients with recurrent breast cancer. Experimental Design The mRNA levels of TSC22D1 and PSAP were analyzed by quantitative real-time polymerase chain reaction (qRT-PCR) in 223 estrogen receptor-positive primary breast tumors of patients with recurrent disease treated with first-line tamoxifen therapy. The main objective of this study was the length of progression-free survival (PFS). Results High mRNA levels of TSC22D1 and PSAP were significantly associated with shorter PFS and both were independent of the traditional predictive factors (HR = 1.30, 95% CI = 1.04-1.64 P = 0.023; and HR = 1.40, 95% CI = 1.03-1.88, P = 0.029, respectively). In multivariate analysis, patients with high mRNA levels of both genes associated significantly with no clinical benefit (OR = 0.19, 95% CI = 0.06-0.62, P = 0.006) and had the shortest PFS (HR = 2.05, 95% CI = 1.29-3.25, P = 0.002). Conclusion These results confirm our previous in vitro and tumor-related findings and are indicative for the failure of tamoxifen treatment in breast-cancer patients. Both TSC22D1 and PSAP are associated with clinical outcome and may have a functional role in therapy resistance.

Downregulation of SIAH2, an ubiquitin E3 ligase, is associated with resistance to endocrine therapy in breast cancer.

PURPOSE: In our microarray analysis we observed that Seven-in-Absentia Homolog 2 (SIAH2) levels were low in estrogen receptor (ER) positive breast tumors of patients resistant to first-line tamoxifen therapy. The aim of this study was to evaluate SIAH2 for its (a) predictive/prognostic value, and (b) functional role in endocrine therapy resistance. PATIENTS AND METHODS: SIAH2 expression was measured with quantitative Real-Time-PCR (qRT-PCR) in 1205 primary breast tumor specimens and related to disease outcome. The functional role of SIAH2 was determined in human breast cancer cell lines ZR-75-1, ZR/HERc, and MCF7. Cell lines were treated with estrogen (E2), anti-estrogen ICI164.384 or epidermal growth factor (EGF). Moreover, MCF7 was treated with ICI164.384 after silencing SIAH2 expression. RESULTS: SIAH2 was not prognostic in 603 lymph node negative patients who had not received adjuvant systemic therapy. In multivariate analysis of ER-positive tumors of 235 patients with recurrent disease, SIAH2 as continuous variable, significantly predicted first-line tamoxifen treatment failure (OR = 1.48; P = 0.05) and progression-free survival (PFS) (HR = 0.79; P = 0.007). Furthermore, in primary breast cancer patients treated with adjuvant tamoxifen, SIAH2 predicted metastasis-free survival (MFS) (HR = 0.73; P = 0.005). In vitro experiments showed that SIAH2 silencing in MCF7 cells resulted in resistance to ICI164.384-treatment when compared with mock silenced cells (P = 0.008). Interestingly, in ZR cells transfected with EGFR (ZR/HERc), SIAH2 expression was induced by E2 but downregulated by EGF. CONCLUSION: In primary breast tumor specimens as well as in vitro low SIAH2 levels associated with resistance to endocrine therapy. Moreover, SIAH2 expression showed an opposite regulation by E2 and EGF.

Cell-free MicroRNA miR-505-3p in Graft Preservation Fluid Is an Independent Predictor of Delayed Graft Function After Kidney Transplantation.

BACKGROUND: Delayed graft function (DGF), a common complication after transplantation of deceased donor kidneys, affects both short- and long-term outcomes. Currently available biomarkers during graft preservation lack sensitivity in predicting risk for DGF. The aim of this study is to identify cell-free micro ribonucleic acid (miRNA) biomarkers in graft preservation fluid predictive of DGF after kidney transplantation. METHODS: Vascular bed preservation fluid was collected from 48 kidney grafts from donation after circulatory death (DCD) or donation after brain death (DBD) donors. miRNA profiles were determined by polymerase chain reaction (PCR) array (n = 8) and validated by reverse transcription and quantitative PCR (n = 40). Graft function posttransplantation was defined as immediate good function (IF) or DGF. RESULTS: A total of 223 miRNAs fulfilled the preset parameters (Ct < 40 in 3 or more samples) and were included in the analysis. Thirty-two miRNAs were significantly different between DGF and IF kidney grafts (P < 0.05) but, after correction for multiple testing, only miR-505-3p remained significant. The significant association of high miR-505-3p levels with DGF was confirmed in an independent validation cohort using conventional reverse transcription and quantitative PCR detection. Multivariate analyses showed miR-505-3p as an independent predictor for DGF (odds ratio, 1.12; P = 0.028). If stratified for donor type, miR-505-3p levels remained significantly different between IF and DGF in DCD grafts (P < 0.01), but not in DBD grafts. Receiver operating characteristic curve analysis showed a high sensitivity and specificity (area under the curve, 0.833). CONCLUSIONS: In DCD grafts, high levels of miR-505-3p in preservation fluid are associated with increased risk of DGF after kidney transplantation. Further study is required to confirm the utility of cell-free miR-505-3p as prognostic biomarker for DGF.

Identification and Validation Model for Informative Liquid Biopsy-Based microRNA Biomarkers: Insights from Germ Cell Tumor In Vitro, In Vivo and Patient-Derived Data.

Liquid biopsy-based biomarkers, such as microRNAs, represent valuable tools for patient management, but often do not make it to integration in the clinic. We aim to explore issues impeding this transition, in the setting of germ cell tumors, for which novel biomarkers are needed. We describe a model for identifying and validating clinically relevant microRNAs for germ cell tumor patients, using both in vitro, in vivo (mouse model) and patient-derived data. Initial wide screening of candidate microRNAs is performed, followed by targeted profiling of potentially relevant biomarkers. We demonstrate the relevance of appropriate (negative) controls, experimental conditions (proliferation), and issues related to sample origin (serum, plasma, cerebral spinal fluid) and pre-analytical variables (hemolysis, contaminants, temperature), all of which could interfere with liquid biopsy-based studies and their conclusions. Finally, we show the value of our identification model in a specific scenario, contradicting the presumed role of miR-375 as marker of teratoma histology in liquid biopsy setting. Our findings indicate other putative microRNAs (miR-885-5p, miR-448 and miR-197-3p) fulfilling this clinical need. The identification model is informative to identify the best candidate microRNAs to pursue in a clinical setting.

Fibroblast growth factor receptor 4 predicts failure on tamoxifen therapy in patients with recurrent breast cancer.

Tamoxifen treatment of estrogen-dependent breast cancer ultimately loses its effectiveness due to the development of resistance. From a functional screen for identifying genes responsible for tamoxifen resistance in human ZR-75-1 breast cancer cells, fibroblast growth factor (FGF) 17 was recovered. The aim of this exploratory study was to assess the predictive value of FGF17 and the receptors FGFR1-4 for the type of response to tamoxifen treatment (clinical benefit) and the duration of progression-free survival (PFS) in patients with recurrent breast cancer. mRNA levels of FGF17 and FGFR1-4 were quantified by real-time reverse transcriptase PCR in 285 estrogen receptor-positive breast carcinomas with clinical follow-up. All patients had recurrent disease and were treated with tamoxifen as first-line systemic therapy for local or distant relapse. FGF17 and FGFR1-3 mRNA levels had no significant predictive value for this group of patients. However, high FGFR4 mRNA levels analyzed as a continuous log-transformed variable predicted poor clinical benefit (odds ratio=1.22; P=0.009) and shorter PFS (hazard ratio=1.18; P<0.001). In addition, in multivariable analysis, the predictive value of FGFR4 was independent from the traditional predictive factors. Our analyses show that FGFR4 may play a role in the biological response of the tumor to tamoxifen treatment. In addition, as altered expression of FGF17 causes tamoxifen resistance in vitro, the FGF signaling pathway could be a valuable target in the treatment of breast cancer patients resistant to endocrine treatment.

TargetClone: A multi-sample approach for reconstructing subclonal evolution of tumors.

Most tumors are composed of a heterogeneous population of subclones. A more detailed insight into the subclonal evolution of these tumors can be helpful to study progression and treatment response. Problematically, tumor samples are typically very heterogeneous, making deconvolving individual tumor subclones a major challenge. To overcome this limitation, reducing heterogeneity, such as by means of microdissections, coupled with targeted sequencing, is a viable approach. However, computational methods that enable reconstruction of the evolutionary relationships require unbiased read depth measurements, which are commonly challenging to obtain in this setting. We introduce TargetClone, a novel method to reconstruct the subclonal evolution tree of tumors from single-nucleotide polymorphism allele frequency and somatic single-nucleotide variant measurements. Furthermore, our method infers copy numbers, alleles and the fraction of the tumor component in each sample. TargetClone was specifically designed for targeted sequencing data obtained from microdissected samples. We demonstrate that our method obtains low error rates on simulated data. Additionally, we show that our method is able to reconstruct expected trees in a testicular germ cell cancer and ovarian cancer dataset. The TargetClone package including tree visualization is written in Python and is publicly available at https://github.com/UMCUGenetics/targetclone.

The MicroRNA-371 Family as Plasma Biomarkers for Monitoring Undifferentiated and Potentially Malignant Human Pluripotent Stem Cells in Teratoma Assays.

Predicting developmental potency and risk of posttransplantation tumor formation by human pluripotent stem cells (hPSCs) and their derivatives largely rely on classical histological analysis of teratomas. Here, we investigated whether an assay based on microRNAs (miRNA) in blood plasma is able to detect potentially malignant elements. Several hPSCs and human malignant germ cell tumor (hGCT) lines were investigated in vitro and in vivo after mouse xenografting. The multiple conventional hPSC lines generated mature teratomas, while xenografts from induced hPSCs (hiPSCs) with reactivated reprogramming transgenes and hGCT lines contained undifferentiated and potentially malignant components. The presence of these elements was reflected in the mRNA and miRNA profiles of the xenografts with OCT3/4 mRNA and the miR-371 and miR-302 families readily detectable. miR-371 family members were also identified in mouse plasma faithfully reporting undifferentiated elements in the xenografts. This study demonstrated that undifferentiated and potentially malignant cells could be detected in vivo.