Growth inhibition of pancreatic cancer by targeted delivery of gemcitabine via fucoidan-coated pH-sensitive liposomes.

Fucoidan-coated pH sensitive liposomes were designed for targeted delivery of gemcitabine (FU-GEM PSL) to treat pancreatic cancer (PC). FU-GEM PSL had a particle size of 175.3 ± 4.9 nm, zeta potential of -19.0 ± 3.7 mV, encapsulation efficiency (EE) of 74.05 ± 0.17 %, and drug loading (DL) of 21.27 ± 0.05 %. Cell experiments in vitro showed that FU-GEM PSL could increase the release of GEM and drug concentration, and could inhibit tumor cell proliferation by affecting the cell cycle. FU-GEM PSL entered cells through macropinocytosis and caveolin-mediated endocytosis to exert effects. Meanwhile, the expression of P-selectin was detected in human tissues, demonstrating the feasibility of targeting FU. Moreover, combined with animal experiments in vivo, FU-GEM PSL could inhibit the development of PC. Furthermore, anti-tumor experiments in vivo carried on BALB/c mice indicated that FU-GEM PSL had tumor suppression abilities and safety. Therefore, FU-GEM PSL is a promising formulation for PC therapy.

Triple-negative Breast Carcinoma With Apocrine and Histiocytoid Features: A Clinicopathologic and Molecular Study of 18 Cases.

Triple-negative breast cancer (TNBC) is a heterogenous group of tumors. Most TNBCs are high-grade aggressive tumors, but a minority of TNBCs are not high grade, with relatively indolent behavior and specific morphologic and molecular features. We performed a clinicopathologic and molecular assessment of 18 non-high-grade TNBCs with apocrine and/or histiocytoid features. All were grade I or II with low Ki-67 (≤20%). Thirteen (72%) showed apocrine features, and 5 (28%) showed histiocytoid and lobular features. In all, 17/18 expressed the androgen receptor, and 13/13 expressed gross cystic disease fluid protein 15. Four (22.2%) patients were treated with neoadjuvant chemotherapy, but none achieved a pathologic complete response. In all, 2/18 patients (11%) had lymph node metastasis at the time of surgery. None of the cases had a recurrence or disease-specific death, with an average follow-up time of 38 months. Thirteen cases were profiled by targeted capture-based next-generation DNA sequencing. Genomic alterations (GAs) were most significant for PI3K-PKB/Akt pathway (69%) genes, including PIK3R1 (23%), PIK3CA (38%), and PTEN (23%), and RTK-RAS pathway (62%) including FGFR4 (46%) and ERBB2 (15%). TP53 GA was seen in only 31% of patients. Our findings support those on high-grade TNBCs with apocrine and/or histiocytoid features as a clinicopathologic and genetically distinct subgroup of TNBC. They can be defined by features including tubule formation, rare mitosis, low Ki-67 (≤20%), triple-negative status, expression of androgen receptor and/or gross cystic disease fluid protein 15, and GA in the PI3K-PKB/Akt and/or RTK-RAS pathway. These tumors are not sensitive to chemotherapy but have favorable clinical behavior. Tumor subtype definitions are the first step to implementing future trial designs to select these patients.

Genomic Complexity Predicts Resistance to Endocrine Therapy and CDK4/6 Inhibition in Hormone Receptor-Positive (HR+)/HER2-Negative Metastatic Breast Cancer.

PURPOSE: Clinical biomarkers to identify patients unlikely to benefit from CDK4/6 inhibition (CDK4/6i) in combination with endocrine therapy (ET) are lacking. We implemented a comprehensive circulating tumor DNA (ctDNA) analysis to identify genomic features for predicting and monitoring treatment resistance. EXPERIMENTAL DESIGN: ctDNA was isolated from 216 plasma samples collected from 51 patients with hormone receptor-positive (HR+)/HER2-negative (HER2-) metastatic breast cancer (MBC) on a phase II trial of palbociclib combined with letrozole or fulvestrant (NCT03007979). Boosted whole-exome sequencing (WES) was performed at baseline and clinical progression to evaluate genomic alterations, mutational signatures, and blood tumor mutational burden (bTMB). Low-pass whole-genome sequencing was performed at baseline and serial timepoints to assess blood copy-number burden (bCNB). RESULTS: High bTMB and bCNB were associated with lack of clinical benefit and significantly shorter progression-free survival (PFS) compared with patients with low bTMB or low bCNB (all P < 0.05). Dominant APOBEC signatures were detected at baseline exclusively in cases with high bTMB (5/13, 38.5%) versus low bTMB (0/37, 0%; P = 0.0006). Alterations in ESR1 were enriched in samples with high bTMB (P = 0.0005). There was a high correlation between bTMB determined by WES and bTMB determined using a 600-gene panel (R = 0.98). During serial monitoring, an increase in bCNB score preceded radiographic progression in 12 of 18 (66.7%) patients. CONCLUSIONS: Genomic complexity detected by noninvasive profiling of bTMB and bCNB predicted poor outcomes in patients treated with ET and CDK4/6i and identified early disease progression before imaging. Novel treatment strategies including immunotherapy-based combinations should be investigated in this population.

Author Correction: Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Ultrasonic disruption extraction, characterization and bioactivities of polysaccharides from wild Armillaria mellea.

An efficient ultrasonic disruption extraction (UDE) of polysaccharides from Armillariamellea (AM) were optimized by response surface methodology (RSM). Under optimum conditions: ultrasonic power 915.00 W, temperature 69.27 °C and time 39.13 min, the crude polysaccharides (AMPs) yield was 19.5%. Two purified fractions AMPs-1-1 and AMPs-2-1 were obtained through anion-exchange and gel chromatography. AMPs-1-1 was a heteropolysaccharide with average molecular weights of 1.23 × 105 Da, and composed of Glc, Gal and GlcA with mole percentages of 89.06%, 9.59% and 1.34%, respectively, owning a backbone structure of (1→)-ß-d-Glcp, (1 â†’ 3,6)-α-d-Glcp and (1 â†’ 3)-ß-d-Glcp residues. AMPs-2-1 was a heteropolysaccharide with average molecular weights of 6.76 × 104 Da, and composed of Glc, Gal, GlcA and Man with mole percentages of 65.28%, 22.87%, 2.87% and 8.98%, containing a main backbone chain of (1 â†’ 3,6)-α-d-Glcp and (1 â†’ 6)-ß-d-Glcp residues. AMPs-2-1 possessed obviously antioxidant activities in terms of stronger scavenging activity against DPPH· and ABTS+⋅, higher FRAP and ORAC value than AMPs-1-1. AMPs-2-1 could promote splenocyte lymphocytes and RAW264.7 macrophages proliferation and enhanced the phagocytosis of macrophages, exhibited significant immunomodulatory activities. These results suggested that UDE is an effective extract technology, and AMPs-2-1 could be explored as potential natural antioxidants and immunomodulatory agents.

High prevalence of focal and multi-focal somatic genetic variants in the human brain.

Somatic mutations during stem cell division are responsible for several cancers. In principle, a similar process could occur during the intense cell proliferation accompanying human brain development, leading to the accumulation of regionally distributed foci of mutations. Using dual platform >5000-fold depth sequencing of 102 genes in 173 adult human brain samples, we detect and validate somatic mutations in 27 of 54 brains. Using a mathematical model of neurodevelopment and approximate Bayesian inference, we predict that macroscopic islands of pathologically mutated neurons are likely to be common in the general population. The detected mutation spectrum also includes DNMT3A and TET2 which are likely to have originated from blood cell lineages. Together, these findings establish developmental mutagenesis as a potential mechanism for neurodegenerative disorders, and provide a novel mechanism for the regional onset and focal pathology in sporadic cases.

Cloche is a bHLH-PAS transcription factor that drives haemato-vascular specification.

Vascular and haematopoietic cells organize into specialized tissues during early embryogenesis to supply essential nutrients to all organs and thus play critical roles in development and disease. At the top of the haemato-vascular specification cascade lies cloche, a gene that when mutated in zebrafish leads to the striking phenotype of loss of most endothelial and haematopoietic cells and a significant increase in cardiomyocyte numbers. Although this mutant has been analysed extensively to investigate mesoderm diversification and differentiation and continues to be broadly used as a unique avascular model, the isolation of the cloche gene has been challenging due to its telomeric location. Here we used a deletion allele of cloche to identify several new cloche candidate genes within this genomic region, and systematically genome-edited each candidate. Through this comprehensive interrogation, we succeeded in isolating the cloche gene and discovered that it encodes a PAS-domain-containing bHLH transcription factor, and that it is expressed in a highly specific spatiotemporal pattern starting during late gastrulation. Gain-of-function experiments show that it can potently induce endothelial gene expression. Epistasis experiments reveal that it functions upstream of etv2 and tal1, the earliest expressed endothelial and haematopoietic transcription factor genes identified to date. A mammalian cloche orthologue can also rescue blood vessel formation in zebrafish cloche mutants, indicating a highly conserved role in vertebrate vasculogenesis and haematopoiesis. The identification of this master regulator of endothelial and haematopoietic fate enhances our understanding of early mesoderm diversification and may lead to improved protocols for the generation of endothelial and haematopoietic cells in vivo and in vitro.

LIN28A Modulates Splicing and Gene Expression Programs in Breast Cancer Cells.

LIN28 is an evolutionarily conserved RNA-binding protein with critical functions in developmental timing and cancer. However, the molecular mechanisms underlying LIN28's oncogenic properties are yet to be described. RNA-protein immunoprecipitation coupled with genome-wide sequencing (RIP-Seq) analysis revealed significant LIN28 binding within 843 mRNAs in breast cancer cells. Many of the LIN28-bound mRNAs are implicated in the regulation of RNA and cell metabolism. We identify heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1), a protein with multiple roles in mRNA metabolism, as a LIN28-interacting partner. Subsequently, we used a custom computational method to identify differentially spliced gene isoforms in LIN28 and hnRNP A1 small interfering RNA (siRNA)-treated cells. The results reveal that these proteins regulate alternative splicing and steady-state mRNA expression of genes implicated in aspects of breast cancer biology. Notably, cells lacking LIN28 undergo significant isoform switching of the ENAH gene, resulting in a decrease in the expression of the ENAH exon 11a isoform. The expression of ENAH isoform 11a has been shown to be elevated in breast cancers that express HER2. Intriguingly, analysis of publicly available array data from the Cancer Genome Atlas (TCGA) reveals that LIN28 expression in the HER2 subtype is significantly different from that in other breast cancer subtypes. Collectively, our data suggest that LIN28 may regulate splicing and gene expression programs that drive breast cancer subtype phenotypes.

mRNA-Seq of single prostate cancer circulating tumor cells reveals recapitulation of gene expression and pathways found in prostate cancer.

Circulating tumor cells (CTC) mediate metastatic spread of many solid tumors and enumeration of CTCs is currently used as a prognostic indicator of survival in metastatic prostate cancer patients. Some evidence suggests that it is possible to derive additional information about tumors from expression analysis of CTCs, but the technical difficulty of isolating and analyzing individual CTCs has limited progress in this area. To assess the ability of a new generation of MagSweeper to isolate intact CTCs for downstream analysis, we performed mRNA-Seq on single CTCs isolated from the blood of patients with metastatic prostate cancer and on single prostate cancer cell line LNCaP cells spiked into the blood of healthy donors. We found that the MagSweeper effectively isolated CTCs with a capture efficiency that matched the CellSearch platform. However, unlike CellSearch, the MagSweeper facilitates isolation of individual live CTCs without contaminating leukocytes. Importantly, mRNA-Seq analysis showed that the MagSweeper isolation process did not have a discernible impact on the transcriptional profile of single LNCaPs isolated from spiked human blood, suggesting that any perturbations caused by the MagSweeper process on the transcriptional signature of isolated cells are modest. Although the RNA from patient CTCs showed signs of significant degradation, consistent with reports of short half-lives and apoptosis amongst CTCs, transcriptional signatures of prostate tissue and of cancer were readily detectable with single CTC mRNA-Seq. These results demonstrate that the MagSweeper provides access to intact CTCs and that these CTCs can potentially supply clinically relevant information.

Monitoring chronic lymphocytic leukemia progression by whole genome sequencing reveals heterogeneous clonal evolution patterns.

Chronic lymphocytic leukemia is characterized by relapse after treatment and chemotherapy resistance. Similarly, in other malignancies leukemia cells accumulate mutations during growth, forming heterogeneous cell populations that are subject to Darwinian selection and may respond differentially to treatment. There is therefore a clinical need to monitor changes in the subclonal composition of cancers during disease progression. Here, we use whole-genome sequencing to track subclonal heterogeneity in 3 chronic lymphocytic leukemia patients subjected to repeated cycles of therapy. We reveal different somatic mutation profiles in each patient and use these to establish probable hierarchical patterns of subclonal evolution, to identify subclones that decline or expand over time, and to detect founder mutations. We show that clonal evolution patterns are heterogeneous in individual patients. We conclude that genome sequencing is a powerful and sensitive approach to monitor disease progression repeatedly at the molecular level. If applied to future clinical trials, this approach might eventually influence treatment strategies as a tool to individualize and direct cancer treatment.

Highly sensitive profiling of CD44+/CD24- breast cancer stem cells by combining global mRNA amplification and next generation sequencing: evidence for a hyperactive PI3K pathway.

We performed next generation sequencing- and microarray-based gene expression profiling of CD44(+)/CD24(-)/CD45(-) breast CSCs (cancer stem cells) isolated from primary ERα-positive breast cancer. By combining semi-automated dissociation of human tumor tissue, magnetic cell sorting and cDNA amplification less than 500 CSCs were required for transcriptome analyses. Besides overexpressing genes involved in maintenance of stemness, the CSCs showed higher levels of genes that drive the PI3K pathway, including EGFR, HB-EGF, PDGFRA/B, PDGF, MET, PIK3CA, PIK3R1 and PIK3R2. This suggests that, in CSCs of ERα-positive breast cancer, the PI3K pathway which is involved in endocrine resistance is hyperactive.

Full-length mRNA-Seq from single-cell levels of RNA and individual circulating tumor cells.

Genome-wide transcriptome analyses are routinely used to monitor tissue-, disease- and cell type­specific gene expression, but it has been technically challenging to generate expression profiles from single cells. Here we describe a robust mRNA-Seq protocol (Smart-Seq) that is applicable down to single cell levels. Compared with existing methods, Smart-Seq has improved read coverage across transcripts, which enhances detailed analyses of alternative transcript isoforms and identification of single-nucleotide polymorphisms. We determined the sensitivity and quantitative accuracy of Smart-Seq for single-cell transcriptomics by evaluating it on total RNA dilution series. We found that although gene expression estimates from single cells have increased noise, hundreds of differentially expressed genes could be identified using few cells per cell type. Applying Smart-Seq to circulating tumor cells from melanomas, we identified distinct gene expression patterns, including candidate biomarkers for melanoma circulating tumor cells. Our protocol will be useful for addressing fundamental biological problems requiring genome-wide transcriptome profiling in rare cells.

Genome sequencing and analysis of the Tasmanian devil and its transmissible cancer.

The Tasmanian devil (Sarcophilus harrisii), the largest marsupial carnivore, is endangered due to a transmissible facial cancer spread by direct transfer of living cancer cells through biting. Here we describe the sequencing, assembly, and annotation of the Tasmanian devil genome and whole-genome sequences for two geographically distant subclones of the cancer. Genomic analysis suggests that the cancer first arose from a female Tasmanian devil and that the clone has subsequently genetically diverged during its spread across Tasmania. The devil cancer genome contains more than 17,000 somatic base substitution mutations and bears the imprint of a distinct mutational process. Genotyping of somatic mutations in 104 geographically and temporally distributed Tasmanian devil tumors reveals the pattern of evolution and spread of this parasitic clonal lineage, with evidence of a selective sweep in one geographical area and persistence of parallel lineages in other populations.

Integrated analysis of gene expression, CpG island methylation, and gene copy number in breast cancer cells by deep sequencing.

We used deep sequencing technology to profile the transcriptome, gene copy number, and CpG island methylation status simultaneously in eight commonly used breast cell lines to develop a model for how these genomic features are integrated in estrogen receptor positive (ER+) and negative breast cancer. Total mRNA sequence, gene copy number, and genomic CpG island methylation were carried out using the Illumina Genome Analyzer. Sequences were mapped to the human genome to obtain digitized gene expression data, DNA copy number in reference to the non-tumor cell line (MCF10A), and methylation status of 21,570 CpG islands to identify differentially expressed genes that were correlated with methylation or copy number changes. These were evaluated in a dataset from 129 primary breast tumors. Gene expression in cell lines was dominated by ER-associated genes. ER+ and ER- cell lines formed two distinct, stable clusters, and 1,873 genes were differentially expressed in the two groups. Part of chromosome 8 was deleted in all ER- cells and part of chromosome 17 amplified in all ER+ cells. These loci encoded 30 genes that were overexpressed in ER+ cells; 9 of these genes were overexpressed in ER+ tumors. We identified 149 differentially expressed genes that exhibited differential methylation of one or more CpG islands within 5 kb of the 5' end of the gene and for which mRNA abundance was inversely correlated with CpG island methylation status. In primary tumors we identified 84 genes that appear to be robust components of the methylation signature that we identified in ER+ cell lines. Our analyses reveal a global pattern of differential CpG island methylation that contributes to the transcriptome landscape of ER+ and ER- breast cancer cells and tumors. The role of gene amplification/deletion appears to more modest, although several potentially significant genes appear to be regulated by copy number aberrations.

Global analysis of estrogen receptor beta binding to breast cancer cell genome reveals an extensive interplay with estrogen receptor alpha for target gene regulation.

BACKGROUND: Estrogen receptors alpha (ERα) and beta (ERß) are transcription factors (TFs) that mediate estrogen signaling and define the hormone-responsive phenotype of breast cancer (BC). The two receptors can be found co-expressed and play specific, often opposite, roles, with ERß being able to modulate the effects of ERα on gene transcription and cell proliferation. ERß is frequently lost in BC, where its presence generally correlates with a better prognosis of the disease. The identification of the genomic targets of ERß in hormone-responsive BC cells is thus a critical step to elucidate the roles of this receptor in estrogen signaling and tumor cell biology. RESULTS: Expression of full-length ERß in hormone-responsive, ERα-positive MCF-7 cells resulted in a marked reduction in cell proliferation in response to estrogen and marked effects on the cell transcriptome. By ChIP-Seq we identified 9702 ERß and 6024 ERα binding sites in estrogen-stimulated cells, comprising sites occupied by either ERß, ERα or both ER subtypes. A search for TF binding matrices revealed that the majority of the binding sites identified comprise one or more Estrogen Response Element and the remaining show binding matrixes for other TFs known to mediate ER interaction with chromatin by tethering, including AP2, E2F and SP1. Of 921 genes differentially regulated by estrogen in ERß+ vs ERß- cells, 424 showed one or more ERß site within 10 kb. These putative primary ERß target genes control cell proliferation, death, differentiation, motility and adhesion, signal transduction and transcription, key cellular processes that might explain the biological and clinical phenotype of tumors expressing this ER subtype. ERß binding in close proximity of several miRNA genes and in the mitochondrial genome, suggests the possible involvement of this receptor in small non-coding RNA biogenesis and mitochondrial genome functions. CONCLUSIONS: Results indicate that the vast majority of the genomic targets of ERß can bind also ERα, suggesting that the overall action of ERß on the genome of hormone-responsive BC cells depends mainly on the relative concentration of both ERs in the cell.

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis.

Monozygotic or 'identical' twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4(+) lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among approximately 3.6 million single nucleotide polymorphisms (SNPs) or approximately 0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of approximately 19,000 genes in CD4(+) T cells. Only 2 to 176 differences in the methylation of approximately 2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to approximately 800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Estrogen receptor alpha controls a gene network in luminal-like breast cancer cells comprising multiple transcription factors and microRNAs.

Luminal-like breast tumor cells express estrogen receptor alpha (ERalpha), a member of the nuclear receptor family of ligand-activated transcription factors that controls their proliferation, survival, and functional status. To identify the molecular determinants of this hormone-responsive tumor phenotype, a comprehensive genome-wide analysis was performed in estrogen stimulated MCF-7 and ZR-75.1 cells by integrating time-course mRNA expression profiling with global mapping of genomic ERalpha binding sites by chromatin immunoprecipitation coupled to massively parallel sequencing, microRNA expression profiling, and in silico analysis of transcription units and receptor binding regions identified. All 1270 genes that were found to respond to 17beta-estradiol in both cell lines cluster in 33 highly concordant groups, each of which showed defined kinetics of RNA changes. This hormone-responsive gene set includes several direct targets of ERalpha and is organized in a gene regulation cascade, stemming from ligand-activated receptor and reaching a large number of downstream targets via AP-2gamma, B-cell activating transcription factor, E2F1 and 2, E74-like factor 3, GTF2IRD1, hairy and enhancer of split homologue-1, MYB, SMAD3, RARalpha, and RXRalpha transcription factors. MicroRNAs are also integral components of this gene regulation network because miR-107, miR-424, miR-570, miR-618, and miR-760 are regulated by 17beta-estradiol along with other microRNAs that can target a significant number of transcripts belonging to one or more estrogen-responsive gene clusters.

Chimeric transcript discovery by paired-end transcriptome sequencing.

Recurrent gene fusions are a prevalent class of mutations arising from the juxtaposition of 2 distinct regions, which can generate novel functional transcripts that could serve as valuable therapeutic targets in cancer. Therefore, we aim to establish a sensitive, high-throughput methodology to comprehensively catalog functional gene fusions in cancer by evaluating a paired-end transcriptome sequencing strategy. Not only did a paired-end approach provide a greater dynamic range in comparison with single read based approaches, but it clearly distinguished the high-level "driving" gene fusions, such as BCR-ABL1 and TMPRSS2-ERG, from potential lower level "passenger" gene fusions. Also, the comprehensiveness of a paired-end approach enabled the discovery of 12 previously undescribed gene fusions in 4 commonly used cell lines that eluded previous approaches. Using the paired-end transcriptome sequencing approach, we observed read-through mRNA chimeras, tissue-type restricted chimeras, converging transcripts, diverging transcripts, and overlapping mRNA transcripts. Last, we successfully used paired-end transcriptome sequencing to detect previously undescribed ETS gene fusions in prostate tumors. Together, this study establishes a highly specific and sensitive approach for accurately and comprehensively cataloguing chimeras within a sample using paired-end transcriptome sequencing.