ArCH: improving the performance of clonal hematopoiesis variant calling and interpretation.

MOTIVATION: The acquisition of somatic mutations in hematopoietic stem and progenitor stem cells with resultant clonal expansion, termed clonal hematopoiesis (CH), is associated with increased risk of hematologic malignancies and other adverse outcomes. CH is generally present at low allelic fractions, but clonal expansion and acquisition of additional mutations leads to hematologic cancers in a small proportion of individuals. With high depth and high sensitivity sequencing, CH can be detected in most adults and its clonal trajectory mapped over time. However, accurate CH variant calling is challenging due to the difficulty in distinguishing low frequency CH mutations from sequencing artifacts. The lack of well-validated bioinformatic pipelines for CH calling may contribute to lack of reproducibility in studies of CH. RESULTS: Here, we developed ArCH, an Artifact filtering Clonal Hematopoiesis variant calling pipeline for detecting single nucleotide variants and short insertions/deletions by combining the output of four variant calling tools and filtering based on variant characteristics and sequencing error rate estimation. ArCH is an end-to-end cloud-based pipeline optimized to accept a variety of inputs with customizable parameters adaptable to multiple sequencing technologies, research questions, and datasets. Using deep targeted sequencing data generated from six acute myeloid leukemia patient tumor: normal dilutions, 31 blood samples with orthogonal validation, and 26 blood samples with technical replicates, we show that ArCH improves the sensitivity and positive predictive value of CH variant detection at low allele frequencies compared to standard application of commonly used variant calling approaches. AVAILABILITY AND IMPLEMENTATION: The code for this workflow is available at: https://github.com/kbolton-lab/ArCH.

Identification of canine circulating miRNAs as tumor biospecific markers using Next-Generation Sequencing and Q-RT-PCR.

Delay in cancer diagnosis often results in metastasis and an inability to successfully treat the tumor. The use of broadly cancer-specific biomarkers at an early stage may improve cancer treatment and staging. This study has explored circulatory exosomal miRNAs as potential diagnostic biomarkers to identify cancer patients. Secretory exosomal miRNAs were isolated from 13 canine cancer cell lines (lymphoma, mast cell tumor, histiocytic cell line, osteosarcoma, melanoma, and breast tumor) and were sequenced by Next-Generation sequencing (NGS). We have identified 6 miRNAs (cfa-miR-9, -1841, -1306, -345, -132, and -26b) by NGS that were elevated in all cancer cell types. The miRNAs identified by NGS were then examined by Q-RT-PCR. The PCR data demonstrated similar expression patterns to those seen with NGS but provided fold differences that were much lower than those seen for NGS. Cfa-miR-9 was found to be the most consistently elevated miRNA in NGS and PCR, making it the most likely miRNA to prove diagnostic. In this study, we have demonstrated that it is possible to identify exosomal miRNAs with elevated secretion across multiple tumor types that could be used as circulatory diagnostic biomarkers for liquid biopsy in the future.

Differential Expression of miRNAs in Hypoxia ("HypoxamiRs") in Three Canine High-Grade Glioma Cell Lines.

Dogs with spontaneous high-grade gliomas increasingly are being proposed as useful large animal pre-clinical models for the human disease. Hypoxia is a critical microenvironmental condition that is common in both canine and human high-grade gliomas and drives increased angiogenesis, chemo- and radioresistance, and acquisition of a stem-like phenotype. Some of this effect is mediated by the hypoxia-induced expression of microRNAs, small (~22 nucleotides long), non-coding RNAs that can modulate gene expression through interference with mRNA translation. Using an in vitro model with three canine high-grade glioma cell lines (J3T, SDT3G, and G06A) exposed to 72 h of 1.5% oxygen vs. standard 20% oxygen, we examined the global "hypoxamiR" profile using small RNA-Seq and performed pathway analysis for targeted genes using both Panther and NetworkAnalyst. Important pathways include many that are well-established as being important in glioma biology, general cancer biology, hypoxia, angiogenesis, immunology, and stem-ness, among others. This work provides the first examination of the effect of hypoxia on miRNA expression in the context of canine glioma, and highlights important similarities with the human disease.

Estimated number of adult survivors of childhood cancer in United States with cancer-predisposing germline variants.

PURPOSE: To estimate the absolute number of adult survivors of childhood cancer in the U.S. population who carry a pathogenic or likely pathogenic variant in a cancer predisposition gene. METHODS: Using the Surveillance, Epidemiology, and End Results (SEER) Program, we estimated the number of childhood cancer survivors on December 31, 2016 for each childhood cancer diagnosis, multiplied this by the proportion of carriers of pathogenic/likely pathogenic variants in the St. Jude Lifetime Cohort (SJLIFE) study, and projected the resulting number onto the U.S. RESULTS: Based on genome sequence data, 11.8% of 2450 SJLIFE participants carry a pathogenic/likely pathogenic variant in one of 156 cancer predisposition genes. Given this information, we estimate that 21 800 adult survivors of childhood cancer in the United States carry a pathogenic/likely pathogenic variant in one of these genes. The highest estimated absolute number of variant carriers are among survivors of central nervous system tumors (n = 4300), particularly astrocytoma (n = 1800) and other gliomas (n = 1700), acute lymphoblastic leukemia (n = 4300), and retinoblastoma (n = 3500). The most frequently mutated genes are RB1 (n = 3000), NF1 (n = 2300), and BRCA2 (n = 800). CONCLUSION: Given the increasing number of childhood cancer survivors in the United States, clinicians should counsel survivors regarding their potential genetic risk, consider referral for genetic counseling and testing, and, as appropriate, implement syndrome-specific cancer surveillance or risk-reducing measures.

The whole-genome landscape of Burkitt lymphoma subtypes.

Burkitt lymphoma (BL) is an aggressive, MYC-driven lymphoma comprising 3 distinct clinical subtypes: sporadic BLs that occur worldwide, endemic BLs that occur predominantly in sub-Saharan Africa, and immunodeficiency-associated BLs that occur primarily in the setting of HIV. In this study, we comprehensively delineated the genomic basis of BL through whole-genome sequencing (WGS) of 101 tumors representing all 3 subtypes of BL to identify 72 driver genes. These data were additionally informed by CRISPR screens in BL cell lines to functionally annotate the role of oncogenic drivers. Nearly every driver gene was found to have both coding and non-coding mutations, highlighting the importance of WGS for identifying driver events. Our data implicate coding and non-coding mutations in IGLL5, BACH2, SIN3A, and DNMT1. Epstein-Barr virus (EBV) infection was associated with higher mutation load, with type 1 EBV showing a higher mutational burden than type 2 EBV. Although sporadic and immunodeficiency-associated BLs had similar genetic profiles, endemic BLs manifested more frequent mutations in BCL7A and BCL6 and fewer genetic alterations in DNMT1, SNTB2, and CTCF. Silencing mutations in ID3 were a common feature of all 3 subtypes of BL. In vitro, mass spectrometry-based proteomics demonstrated that the ID3 protein binds primarily to TCF3 and TCF4. In vivo knockout of ID3 potentiated the effects of MYC, leading to rapid tumorigenesis and tumor phenotypes consistent with those observed in the human disease.

Analysis of error profiles in deep next-generation sequencing data.

BACKGROUND: Sequencing errors are key confounding factors for detecting low-frequency genetic variants that are important for cancer molecular diagnosis, treatment, and surveillance using deep next-generation sequencing (NGS). However, there is a lack of comprehensive understanding of errors introduced at various steps of a conventional NGS workflow, such as sample handling, library preparation, PCR enrichment, and sequencing. In this study, we use current NGS technology to systematically investigate these questions. RESULTS: By evaluating read-specific error distributions, we discover that the substitution error rate can be computationally suppressed to 10-5 to 10-4, which is 10- to 100-fold lower than generally considered achievable (10-3) in the current literature. We then quantify substitution errors attributable to sample handling, library preparation, enrichment PCR, and sequencing by using multiple deep sequencing datasets. We find that error rates differ by nucleotide substitution types, ranging from 10-5 for A>C/T>G, C>A/G>T, and C>G/G>C changes to 10-4 for A>G/T>C changes. Furthermore, C>T/G>A errors exhibit strong sequence context dependency, sample-specific effects dominate elevated C>A/G>T errors, and target-enrichment PCR led to ~ 6-fold increase of overall error rate. We also find that more than 70% of hotspot variants can be detected at 0.1 ~ 0.01% frequency with the current NGS technology by applying in silico error suppression. CONCLUSIONS: We present the first comprehensive analysis of sequencing error sources in conventional NGS workflows. The error profiles revealed by our study highlight new directions for further improving NGS analysis accuracy both experimentally and computationally, ultimately enhancing the precision of deep sequencing.

Ectonucleoside Triphosphate Diphosphohydrolase-3 Antibody Targets Adult Human Pancreatic ß Cells for In Vitro and In Vivo Analysis.

Identification of cell-surface markers specific to human pancreatic ß cells would allow in vivo analysis and imaging. Here we introduce a biomarker, ectonucleoside triphosphate diphosphohydrolase-3 (NTPDase3), that is expressed on the cell surface of essentially all adult human ß cells, including those from individuals with type 1 or type 2 diabetes. NTPDase3 is expressed dynamically during postnatal human pancreas development, appearing first in acinar cells at birth, but several months later its expression declines in acinar cells while concurrently emerging in islet ß cells. Given its specificity and membrane localization, we utilized an NTPDase3 antibody for purification of live human ß cells as confirmed by transcriptional profiling, and, in addition, for in vivo imaging of transplanted human ß cells. Thus, NTPDase3 is a cell-surface biomarker of adult human ß cells, and the antibody directed to this protein should be a useful new reagent for ß cell sorting, in vivo imaging, and targeting.

Diverse Long RNAs Are Differentially Sorted into Extracellular Vesicles Secreted by Colorectal Cancer Cells.

The regulation and functional roles of secreted coding and long noncoding RNAs (lncRNAs; >200 nt) are largely unknown. We previously showed that mutant KRAS colorectal cancer (CRC) cells release extracellular vesicles (EVs) containing distinct proteomes, microRNAs (miRNAs), and circular RNAs. Here, we comprehensively identify diverse classes of CRC extracellular long RNAs secreted in EVs and demonstrate differential export of specific RNAs. Distinct noncoding RNAs, including antisense transcripts and transcripts derived from pseudogenes, are enriched in EVs compared to cellular profiles. We detected strong enrichment of Rab13 in mutant KRAS EVs and demonstrate functional delivery of Rab13 mRNA to recipient cells. To assay functional transfer of lncRNAs, we implemented a CRISPR/Cas9-based RNA-tracking system to monitor delivery to recipient cells. We show that gRNAs containing export signals from secreted RNAs can be transferred from donor to recipient cells. Our data support the existence of cellular mechanisms to selectively export diverse classes of RNA.

Genetic Risk for Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer.

Purpose Childhood cancer survivors are at increased risk of subsequent neoplasms (SNs), but the germline genetic contribution is largely unknown. We assessed the contribution of pathogenic/likely pathogenic (P/LP) mutations in cancer predisposition genes to their SN risk. Patients and Methods Whole-genome sequencing (30-fold) was performed on samples from childhood cancer survivors who were ≥ 5 years since initial cancer diagnosis and participants in the St Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Germline mutations in 60 genes known to be associated with autosomal dominant cancer predisposition syndromes with moderate to high penetrance were classified by their pathogenicity according to the American College of Medical Genetics and Genomics guidelines. Relative rates (RRs) and 95% CIs of SN occurrence by mutation status were estimated using multivariable piecewise exponential regression stratified by radiation exposure. Results Participants were 3,006 survivors (53% male; median age, 35.8 years [range, 7.1 to 69.8 years]; 56% received radiotherapy), 1,120 SNs were diagnosed among 439 survivors (14.6%), and 175 P/LP mutations were identified in 5.8% (95% CI, 5.0% to 6.7%) of survivors. Mutations were associated with significantly increased rates of breast cancer (RR, 13.9; 95% CI, 6.0 to 32.2) and sarcoma (RR, 10.6; 95% CI, 4.3 to 26.3) among irradiated survivors and with increased rates of developing any SN (RR, 4.7; 95% CI, 2.4 to 9.3), breast cancer (RR, 7.7; 95% CI, 2.4 to 24.4), nonmelanoma skin cancer (RR, 11.0; 95% CI, 2.9 to 41.4), and two or more histologically distinct SNs (RR, 18.6; 95% CI, 3.5 to 99.3) among nonirradiated survivors. Conclusion The findings support referral of all survivors for genetic counseling for potential clinical genetic testing, which should be prioritized for nonirradiated survivors with any SN and for those with breast cancer or sarcoma in the field of prior irradiation.

A signalling cascade of IL-33 to IL-13 regulates metaplasia in the mouse stomach.

OBJECTIVE: Alternatively activated macrophages (M2) are associated with the progression of spasmolytic polypeptide-expressing metaplasia (SPEM) in the stomach. However, the precise mechanism(s) and critical mediators that induce SPEM are unknown. DESIGN: To determine candidate genes important in these processes, macrophages from the stomach corpus of mice with SPEM (DMP-777-treated) or advanced SPEM (L635-treated) were isolated and RNA sequenced. Effects on metaplasia development after acute parietal cell loss induced by L635 were evaluated in interleukin (IL)-33, IL-33 receptor (ST2) and IL-13 knockout (KO) mice. RESULTS: Profiling of metaplasia-associated macrophages in the stomach identified an M2a-polarised macrophage population. Expression of IL-33 was significantly upregulated in macrophages associated with advanced SPEM. L635 induced metaplasia in the stomachs of wild-type mice, but not in the stomachs of IL-33 and ST2 KO mice. While IL-5 and IL-9 were not required for metaplasia induction, IL-13 KO mice did not develop metaplasia in response to L635. Administration of IL-13 to ST2 KO mice re-established the induction of metaplasia following acute parietal cell loss. CONCLUSIONS: Metaplasia induction and macrophage polarisation after parietal cell loss is coordinated through a cytokine signalling network of IL-33 and IL-13, linking a combined response to injury by both intrinsic mucosal mechanisms and infiltrating M2 macrophages.

Age-dependent human ß cell proliferation induced by glucagon-like peptide 1 and calcineurin signaling.

Inadequate pancreatic ß cell function underlies type 1 and type 2 diabetes mellitus. Strategies to expand functional cells have focused on discovering and controlling mechanisms that limit the proliferation of human ß cells. Here, we developed an engraftment strategy to examine age-associated human islet cell replication competence and reveal mechanisms underlying age-dependent decline of ß cell proliferation in human islets. We found that exendin-4 (Ex-4), an agonist of the glucagon-like peptide 1 receptor (GLP-1R), stimulates human ß cell proliferation in juvenile but not adult islets. This age-dependent responsiveness does not reflect loss of GLP-1R signaling in adult islets, since Ex-4 treatment stimulated insulin secretion by both juvenile and adult human ß cells. We show that the mitogenic effect of Ex-4 requires calcineurin/nuclear factor of activated T cells (NFAT) signaling. In juvenile islets, Ex-4 induced expression of calcineurin/NFAT signaling components as well as target genes for proliferation-promoting factors, including NFATC1, FOXM1, and CCNA1. By contrast, expression of these factors in adult islet ß cells was not affected by Ex-4 exposure. These studies reveal age-dependent signaling mechanisms regulating human ß cell proliferation, and identify elements that could be adapted for therapeutic expansion of human ß cells.

Comparative analysis of chimeric ZFP-, TALE- and Cas9-piggyBac transposases for integration into a single locus in human cells.

Integrating DNA delivery systems hold promise for many applications including treatment of diseases; however, targeted integration is needed for improved safety. The piggyBac (PB) transposon system is a highly active non-viral gene delivery system capable of integrating defined DNA segments into host chromosomes without requiring homologous recombination. We systematically compared four different engineered zinc finger proteins (ZFP), four transcription activator-like effector proteins (TALE), CRISPR associated protein 9 (SpCas9) and the catalytically inactive dSpCas9 protein fused to the amino-terminus of the transposase enzyme designed to target the hypoxanthine phosphoribosyltransferase (HPRT) gene located on human chromosome X. Chimeric transposases were evaluated for expression, transposition activity, chromatin immunoprecipitation at the target loci, and targeted knockout of the HPRT gene in human cells. One ZFP-PB and one TALE-PB chimera demonstrated notable HPRT gene targeting. In contrast, Cas9/dCas9-PB chimeras did not result in gene targeting. Instead, the HPRT locus appeared to be protected from transposon integration. Supplied separately, PB permitted highly efficient isolation of Cas9-mediated knockout of HPRT, with zero transposon integrations in HPRT by deep sequencing. In summary, these tools may allow isolation of 'targeted-only' cells, be utilized to protect a genomic locus from transposon integration, and enrich for Cas9-mutated cells.

A high-throughput molecular data resource for cutaneous neurofibromas.

Neurofibromatosis type 1 (NF1) is a genetic disorder with a range of clinical manifestations such as widespread growth of benign tumours called neurofibromas, pain, learning disorders, bone deformities, vascular abnormalities and even malignant tumours. With the establishment of the Children's Tumour Foundation biobank, neurofibroma samples can now be collected directly from patients to be analysed by the larger scientific community. This work describes a pilot study to characterize one class of neurofibroma, cutaneous neurofibromas, by molecularly profiling of ~40 cutaneous neurofibromas collected from 11 individual patients. Data collected from each tumour includes (1) SNP Arrays, (2) Whole genome sequencing (WGS) and (3) RNA-Sequencing. These data are now freely available for further analysis at http://www.synapse.org/cutaneousNF.

Enteropathy-associated T cell lymphoma subtypes are characterized by loss of function of SETD2.

Enteropathy-associated T cell lymphoma (EATL) is a lethal, and the most common, neoplastic complication of celiac disease. Here, we defined the genetic landscape of EATL through whole-exome sequencing of 69 EATL tumors. SETD2 was the most frequently silenced gene in EATL (32% of cases). The JAK-STAT pathway was the most frequently mutated pathway, with frequent mutations in STAT5B as well as JAK1, JAK3, STAT3, and SOCS1 We also identified mutations in KRAS, TP53, and TERT Type I EATL and type II EATL (monomorphic epitheliotropic intestinal T cell lymphoma) had highly overlapping genetic alterations indicating shared mechanisms underlying their pathogenesis. We modeled the effects of SETD2 loss in vivo by developing a T cell-specific knockout mouse. These mice manifested an expansion of γδ T cells, indicating novel roles for SETD2 in T cell development and lymphomagenesis. Our data render the most comprehensive genetic portrait yet of this uncommon but lethal disease and may inform future classification schemes.

Cell-Based Systems Biology Analysis of Human AS03-Adjuvanted H5N1 Avian Influenza Vaccine Responses: A Phase I Randomized Controlled Trial.

BACKGROUND: Vaccine development for influenza A/H5N1 is an important public health priority, but H5N1 vaccines are less immunogenic than seasonal influenza vaccines. Adjuvant System 03 (AS03) markedly enhances immune responses to H5N1 vaccine antigens, but the underlying molecular mechanisms are incompletely understood. OBJECTIVE AND METHODS: We compared the safety (primary endpoint), immunogenicity (secondary), gene expression (tertiary) and cytokine responses (exploratory) between AS03-adjuvanted and unadjuvanted inactivated split-virus H5N1 influenza vaccines. In a double-blinded clinical trial, we randomized twenty adults aged 18-49 to receive two doses of either AS03-adjuvanted (n = 10) or unadjuvanted (n = 10) H5N1 vaccine 28 days apart. We used a systems biology approach to characterize and correlate changes in serum cytokines, antibody titers, and gene expression levels in six immune cell types at 1, 3, 7, and 28 days after the first vaccination. RESULTS: Both vaccines were well-tolerated. Nine of 10 subjects in the adjuvanted group and 0/10 in the unadjuvanted group exhibited seroprotection (hemagglutination inhibition antibody titer > 1:40) at day 56. Within 24 hours of AS03-adjuvanted vaccination, increased serum levels of IL-6 and IP-10 were noted. Interferon signaling and antigen processing and presentation-related gene responses were induced in dendritic cells, monocytes, and neutrophils. Upregulation of MHC class II antigen presentation-related genes was seen in neutrophils. Three days after AS03-adjuvanted vaccine, upregulation of genes involved in cell cycle and division was detected in NK cells and correlated with serum levels of IP-10. Early upregulation of interferon signaling-related genes was also found to predict seroprotection 56 days after first vaccination. CONCLUSIONS: Using this cell-based systems approach, novel mechanisms of action for AS03-adjuvanted pandemic influenza vaccination were observed. TRIAL REGISTRATION: ClinicalTrials.gov NCT01573312.

Genomic regulation of invasion by STAT3 in triple negative breast cancer.

Breast cancer is a heterogeneous disease comprised of four molecular subtypes defined by whether the tumor-originating cells are luminal or basal epithelial cells. Breast cancers arising from the luminal mammary duct often express estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth receptor 2 (HER2). Tumors expressing ER and/or PR are treated with anti-hormonal therapies, while tumors overexpressing HER2 are targeted with monoclonal antibodies. Immunohistochemical detection of ER, PR, and HER2 receptors/proteins is a critical step in breast cancer diagnosis and guided treatment. Breast tumors that do not express these proteins are known as "triple negative breast cancer" (TNBC) and are typically basal-like. TNBCs are the most aggressive subtype, with the highest mortality rates and no targeted therapy, so there is a pressing need to identify important TNBC tumor regulators. The signal transducer and activator of transcription 3 (STAT3) transcription factor has been previously implicated as a constitutively active oncogene in TNBC. However, its direct regulatory gene targets and tumorigenic properties have not been well characterized. By integrating RNA-seq and ChIP-seq data from 2 TNBC tumors and 5 cell lines, we discovered novel gene signatures directly regulated by STAT3 that were enriched for processes involving inflammation, immunity, and invasion in TNBC. Functional analysis revealed that STAT3 has a key role regulating invasion and metastasis, a characteristic often associated with TNBC. Our findings suggest therapies targeting STAT3 may be important for preventing TNBC metastasis.

A comparison of microRNA expression profiles from splenic hemangiosarcoma, splenic nodular hyperplasia, and normal spleens of dogs.

BACKGROUND: Splenic masses are common in older dogs; yet diagnosis preceding splenectomy and histopathology remains elusive. MicroRNAs (miRNAs) are short, non-coding RNAs that play a role in post-transcriptional regulation, and differential expression of miRNAs between normal and tumor tissue has been used to diagnose neoplastic diseases. The objective of this study was to determine differential expression of miRNAs by use of RNA-sequencing in canine spleens that were histologically confirmed as hemangiosarcoma, nodular hyperplasia, or normal. RESULTS: Twenty-two miRNAs were found to be differentially expressed in hemangiosarcoma samples (4 between hemangiosarcoma and both nodular hyperplasia and normal spleen and 18 between hemangiosarcoma and normal spleen only). In particular, mir-26a, mir-126, mir-139, mir-140, mir-150, mir-203, mir-424, mir-503, mir-505, mir-542, mir-30e, mir-33b, mir-365, mir-758, mir-22, and mir-452 are of interest in the pathogenesis of hemangiosarcoma. CONCLUSIONS: Findings of this study confirm the hypothesis that miRNA expression profiles are different between canine splenic hemangiosarcoma, nodular hyperplasia, and normal spleens. A large portion of the differentially expressed miRNAs have roles in angiogenesis, with an additional group of miRNAs being dysregulated in vascular disease processes. Two other miRNAs have been implicated in cancer pathways such as PTEN and cell cycle checkpoints. The finding of multiple miRNAs with roles in angiogenesis and vascular disease is important, as hemangiosarcoma is a tumor of endothelial cells, which are driven by angiogenic stimuli. This study shows that miRNA dysregulation is a potential player in the pathogenesis of canine splenic hemangiosarcoma.

Circular RNAs are down-regulated in KRAS mutant colon cancer cells and can be transferred to exosomes.

Recent studies have shown that circular RNAs (circRNAs) are abundant, widely expressed in mammals, and can display cell-type specific expression. However, how production of circRNAs is regulated and their precise biological function remains largely unknown. To study how circRNAs might be regulated during colorectal cancer progression, we used three isogenic colon cancer cell lines that differ only in KRAS mutation status. Cellular RNAs from the parental DLD-1 cells that contain both wild-type and G13D mutant KRAS alleles and isogenically-matched derivative cell lines, DKO-1 (mutant KRAS allele only) and DKs-8 (wild-type KRAS allele only) were analyzed using RNA-Seq. We developed a bioinformatics pipeline to identify and evaluate circRNA candidates from RNA-Seq data. Hundreds of high-quality circRNA candidates were identified in each cell line. Remarkably, circRNAs were significantly down-regulated at a global level in DLD-1 and DKO-1 cells compared to DKs-8 cells, indicating a widespread effect of mutant KRAS on circRNA abundance. This finding was confirmed in two independent colon cancer cell lines HCT116 (KRAS mutant) and HKe3 (KRAS WT). In all three cell lines, circRNAs were also found in secreted extracellular-vesicles, and circRNAs were more abundant in exosomes than cells. Our results suggest that circRNAs may serve as promising cancer biomarkers.

Connecting the Dots: Therapy-Induced Senescence and a Tumor-Suppressive Immune Microenvironment.

BACKGROUND: Tumor cell senescence is a common outcome of anticancer therapy. Here we investigated how therapy-induced senescence (TIS) affects tumor-infiltrating leukocytes (TILs) and the efficacy of immunotherapy in melanoma. METHODS: Tumor senescence was induced by AURKA or CDK4/6 inhibitors (AURKAi, CDK4/6i). Transcriptomes of six mouse tumors with differential response to AURKAi were analyzed by RNA sequencing, and TILs were characterized by flow cytometry. Chemokine RNA and protein expression were determined by quantitative real-time polymerase chain reaction and enzyme-linked immunosorbent assay. Therapeutic response was queried in immunodeficient mice, in mice with CCL5-deficient tumors, and in mice cotreated with CD137 agonist to activate TILs. CCL5 expression in reference to TIS and markers of TILs was studied in human melanoma tumors using patient-derived xenografts (n = 3 patients, n = 3 mice each), in AURKAi clinical trial samples (n = 3 patients, before/after therapy), and in The Cancer Genome Atlas (n = 278). All statistical tests were two-sided. RESULTS: AURKAi response was associated with induction of the immune transcriptome (P = 3.5 x 10-29) while resistance inversely correlated with TIL numbers (Spearman r = -0.87, P < .001). AURKAi and CDK4/6i promoted the recruitment of TILs by inducing CCL5 secretion in melanoma cells (P ≤ .005) in an NF-κB-dependent manner. Therapeutic response to AURKAi was impaired in immunodeficient compared with immunocompetent mice (0% vs 67% tumors regressed, P = .01) and in mice bearing CCL5-deficient vs control tumors (P = .61 vs P = .02); however, AURKAi response was greatly enhanced in mice also receiving T-cell-activating immunotherapy (P < .001). In human tumors, CCL5 expression was also induced by AURKAi (P ≤ .02) and CDK4/6i (P = .01) and was associated with increased immune marker expression (P = 1.40 x 10-93). CONCLUSIONS: Senescent melanoma cells secret CCL5, which promotes recruitment of TILs. Combining TIS with immunotherapy that enhances tumor cell killing by TILs is a promising novel approach to improve melanoma outcomes.

Age and Nursing Affect the Neonatal Porcine Uterine Transcriptome.

The lactocrine hypothesis for maternal programming of neonatal development was proposed to describe a mechanism through which milk-borne bioactive factors, delivered from mother to nursing offspring, could affect development of tissues, including the uterus. Porcine uterine development, initiated before birth, is completed postnatally. However, age- and lactocrine-sensitive elements of the neonatal porcine uterine developmental program are undefined. Here, effects of age and nursing on the uterine transcriptome for 48 h from birth (Postnatal Day [PND] = 0) were identified using RNA sequencing (RNAseq). Uterine tissues were obtained from neonatal gilts (n = 4 per group) within 1 h of birth and before feeding (PND 0), or 48 h after nursing ad libitum (PND 2N) or feeding a commercial milk replacer (PND 2R). RNAseq analysis revealed differentially expressed genes (DEGs) associated with both age (PND 2N vs. PND 0; 3283 DEGs) and nursing on PND 2 (PND 2N vs PND 2R; 896 DEGs). Expression of selected uterine genes was validated using quantitative real-time PCR. Bioinformatic analyses revealed multiple biological processes enriched in response to both age and nursing, including cell adhesion, morphogenesis, and cell-cell signaling. Age-sensitive pathways also included estrogen receptor-alpha and hedgehog signaling cascades. Lactocrine-sensitive processes in nursed gilts included those involved in response to wounding, the plasminogen activator network and coagulation. Overall, RNAseq analysis revealed comprehensive age- and nursing-related transcriptomic differences in the neonatal porcine uterus and identified novel pathways and biological processes regulating uterine development.