Identification of Circular RNAs using RNA Sequencing.

Circular RNAs (circRNAs) are a class of non-coding RNAs involved in functions including micro-RNA (miRNA) regulation, mediation of protein-protein interactions, and regulation of parental gene transcription. In classical next generation RNA sequencing (RNA-seq), circRNAs are typically overlooked as a result of poly-A selection during construction of mRNA libraries, or are found at very low abundance, and are therefore difficult to isolate and detect. Here, a circRNA library construction protocol was optimized by comparing library preparation kits, pre-treatment options and various total RNA input amounts. Two commercially available whole transcriptome library preparation kits, with and without RNase R pre-treatment, and using variable amounts of total RNA input (1 to 4 µg), were tested. Lastly, multiple tissue types; including liver, lung, lymph node, and pancreas; as well as multiple brain regions; including the cerebellum, inferior parietal lobe, middle temporal gyrus, occipital cortex, and superior frontal gyrus; were compared to evaluate circRNA abundance across tissue types. Analysis of the generated RNA-seq data using six different circRNA detection tools (find_circ, CIRI, Mapsplice, KNIFE, DCC, and CIRCexplorer) revealed that a stranded total RNA library preparation kit with RNase R pre-treatment and 4 µg RNA input is the optimal method for identifying the highest relative number of circRNAs. Consistent with previous findings, the highest enrichment of circRNAs was observed in brain tissues compared to other tissue types.

Multiscale, multimodal analysis of tumor heterogeneity in IDH1 mutant vs wild-type diffuse gliomas.

Glioma is recognized to be a highly heterogeneous CNS malignancy, whose diverse cellular composition and cellular interactions have not been well characterized. To gain new clinical- and biological-insights into the genetically-bifurcated IDH1 mutant (mt) vs wildtype (wt) forms of glioma, we integrated data from protein, genomic and MR imaging from 20 treatment-naïve glioma cases and 16 recurrent GBM cases. Multiplexed immunofluorescence (MxIF) was used to generate single cell data for 43 protein markers representing all cancer hallmarks, Genomic sequencing (exome and RNA (normal and tumor) and magnetic resonance imaging (MRI) quantitative features (protocols were T1-post, FLAIR and ADC) from whole tumor, peritumoral edema and enhancing core vs equivalent normal region were also collected from patients. Based on MxIF analysis, 85,767 cells (glioma cases) and 56,304 cells (GBM cases) were used to generate cell-level data for 24 biomarkers. K-means clustering was used to generate 7 distinct groups of cells with divergent biomarker profiles and deconvolution was used to assign RNA data into three classes. Spatial and molecular heterogeneity metrics were generated for the cell data. All features were compared between IDH mt and IDHwt patients and were finally combined to provide a holistic/integrated comparison. Protein expression by hallmark was generally lower in the IDHmt vs wt patients. Molecular and spatial heterogeneity scores for angiogenesis and cell invasion also differed between IDHmt and wt gliomas irrespective of prior treatment and tumor grade; these differences also persisted in the MR imaging features of peritumoral edema and contrast enhancement volumes. A coherent picture of enhanced angiogenesis in IDHwt tumors was derived from multiple platforms (genomic, proteomic and imaging) and scales from individual proteins to cell clusters and heterogeneity, as well as bulk tumor RNA and imaging features. Longer overall survival for IDH1mt glioma patients may reflect mutation-driven alterations in cellular, molecular, and spatial heterogeneity which manifest in discernable radiological manifestations.

E6201, an intravenous MEK1 inhibitor, achieves an exceptional response in BRAF V600E-mutated metastatic malignant melanoma with brain metastases.

Malignant melanoma (MM) exhibits a high propensity for central nervous system dissemination with ~50% of metastatic MM patients developing brain metastases (BM). Targeted therapies and immune checkpoint inhibitors have improved overall survival for MM patients with BM. However, responses are usually of short duration and new agents that effectively penetrate the blood brain barrier (BBB) are needed. Here, we report a MM patient with BM who experienced an exceptional response to E6201, an ATP-competitive MEK1 inhibitor, on a Phase 1 study, with ongoing near-complete response and overall survival extending beyond 8 years. Whole exome and transcriptome sequencing revealed a high mutational burden tumor (22 mutations/Megabase) with homozygous BRAF V600E mutation. Correlative preclinical studies demonstrated broad activity for E6201 across BRAF V600E mutant melanoma cell lines and effective BBB penetration in vivo. Together, these results suggest that E6201 may represent a potential new treatment option for BRAF-mutant MM patients with BM.

Identification of therapeutic targets in chordoma through comprehensive genomic and transcriptomic analyses.

Chordoma is a rare, orphan cancer arising from embryonal precursors of bone. Surgery and radiotherapy (RT) provide excellent local control, often at the price of significant morbidity because of the structures involved and the need for relatively high doses of RT; however, recurrence remains high. Although our understanding of the genetic changes that occur in chordoma is evolving rapidly, this knowledge has yet to translate into treatments. We performed comprehensive DNA (paired tumor/normal whole-exome and shallow whole-genome) and RNA (tumor whole-transcriptome) next-generation sequencing analyses of archival sacral and clivus chordoma specimens. Incorporation of transcriptomic data enabled the identification of gene overexpression and expressed DNA alterations, thus providing additional support for potential therapeutic targets. In three patients, we identified alterations that may be amenable to off-label FDA-approved treatments for other tumor types. These alterations include FGFR1 overexpression (ponatinib, pazopanib) and copy-number duplication of CDK4 (palbociclib) and ERBB3 (gefitinib). In a third patient, germline DNA demonstrated predicted pathogenic changes in CHEK2 and ATM, which may have predisposed the patient to developing chordoma at a young age and may also be associated with potential sensitivity to PARP inhibitors because of homologous recombination repair deficiency. Last, in the fourth patient, a missense mutation in IGF1R was identified, suggesting potential activity for investigational anti-IGF1R strategies. Our findings demonstrate that chordoma patients present with aberrations in overlapping pathways. These results provide support for targeting the IGF1R/FGFR/EGFR and CDK4/6 pathways as treatment strategies for chordoma patients. This study underscores the value of comprehensive genomic and transcriptomic analysis in the development of rational, individualized treatment plans for chordoma.

Circular RNA expression and regulatory network prediction in posterior cingulate astrocytes in elderly subjects.

BACKGROUND: Circular RNAs (circRNAs) are a novel class of endogenous, non-coding RNAs that form covalently closed continuous loops and that are both highly conserved and abundant in the mammalian brain. A role for circRNAs in sponging microRNAs (miRNAs) has been proposed, but the circRNA-miRNA-mRNA interaction networks in human brain cells have not been defined. Therefore, we identified circRNAs in RNA sequencing data previously generated from astrocytes microdissected from the posterior cingulate (PC) of Alzheimer's disease (AD) patients (N = 10) and healthy elderly controls (N = 10) using four circRNA prediction algorithms - CIRI, CIRCexplorer, find_circ and KNIFE. RESULTS: Overall, utilizing these four tools, we identified a union of 4438 unique circRNAs across all samples, of which 70.3% were derived from exonic regions. Notably, the widely reported CDR1as circRNA was detected in all samples across both groups by find_circ. Given the putative miRNA regulatory function of circRNAs, we identified potential miRNA targets of circRNAs, and further, delineated circRNA-miRNA-mRNA networks using in silico methods. Pathway analysis of the genes regulated by these miRNAs identified significantly enriched immune response pathways, which is consistent with the known function of astrocytes as immune sensors in the brain. CONCLUSIONS: In this study, we performed circRNA detection on cell-specific transcriptomic data and identified potential circRNA-miRNA-mRNA regulatory networks in PC astrocytes. Given the known function of astrocytes in cerebral innate immunity and our identification of significantly enriched immune response pathways, the circRNAs we identified may be associated with such key functions. While we did not detect recurrent differentially expressed circRNAs in the context of healthy controls or AD, we report for the first time circRNAs and their potential regulatory impact in a cell-specific and region-specific manner in aged subjects. These predicted regulatory network and pathway analyses may help provide new insights into transcriptional regulation in the brain.

Oligomeric amyloid ß preferentially targets neuronal and not glial mitochondrial-encoded mRNAs.

INTRODUCTION: Our laboratories have demonstrated that accumulation of oligomeric amyloid ß (OAß) in neurons is an essential step leading to OAß-mediated mitochondrial dysfunction. METHODS: Alzheimer's disease (AD) and matching control hippocampal neurons, astrocytes, and microglia were isolated by laser-captured microdissection from the same subjects, followed by whole-transcriptome sequencing. Complementary in vitro work was performed in OAß-treated differentiated SH-SY5Y, followed by the use of a novel CoQ10 analogue for protection. This compound is believed to be effective both in suppressing reactive oxygen species and also functioning in mitochondrial electron transport. RESULTS: We report decreases in the same mitochondrial-encoded mRNAs in Alzheimer's disease laser-captured CA1 neurons and in OAß-treated SH-SY5Y cells, but not in laser-captured microglia and astrocytes. Pretreatment with a novel CoQ10 analogue, protects neuronal mitochondria from OAß-induced mitochondrial changes. DISCUSSION: Similarity of expression changes in neurons from Alzheimer's disease brain and neuronal cells treated with OAß, and the effect of a CoQ10 analogue on the latter, suggests a pretreatment option to prevent OAß toxicity, long before the damage is apparent.

Whole Genome Library Construction for Next Generation Sequencing.

With the rapid evolution of genomics technologies over the past decade, whole genome sequencing (WGS) has become an increasingly accessible tool in biomedical research. WGS applications include analysis of genomic DNA from single individuals, multiple related family members, and tumor/normal samples from the same patient in the context of oncology. A number of different modalities are available for performing WGS; this chapter focuses on wet lab library construction procedures for complex short insert WGS libraries using the KAPA Hyper Prep Kit (Kapa Biosystems), and includes a discussion of appropriate quality control measures for sequencing on the Illumina HiSeq2000 platform. Additional modifications to the protocol for long insert WGS library construction, to assess structural alterations and copy number changes, are also described.

Whole Exome Library Construction for Next Generation Sequencing.

Whole exome sequencing (WES) is a DNA sequencing strategy that provides a survey of base substitutions across coding genomic locations and other regions of interest. As the coding portion of the genome encompasses only 1-2% of the entire genome, this approach represents a more cost-effective strategy to detect DNA alterations that may alter protein function, compared to whole genome sequencing. Although the research community has and is currently delineating the functional implications of sequence changes in noncoding regions of the genome, WES is a currently available assay that provides valuable information for both discovery research and precision medicine applications. In this chapter, we present a WES library preparation protocol using the KAPA Hyper Prep Kit with Agilent SureSelect Human All Exon V5+UTR probes that demonstrates high DNA-to-library conversion efficiency for sequencing on the Illumina HiSeq platform.

Laser-captured microglia in the Alzheimer's and Parkinson's brain reveal unique regional expression profiles and suggest a potential role for hepatitis B in the Alzheimer's brain.

Expression array data from dozens of laboratories, including our own, show significant changes in expression of many genes in Alzheimer's disease (AD) patients compared with normal controls. These data typically rely on brain homogenates, and information about transcripts specific to microglia and other central nervous system (CNS) cell types, which far outnumber microglia-specific transcripts, is lost. We therefore used single-cell laser capture methods to assess the full range of microglia-specific expression changes that occur in different brain regions (substantia nigra and hippocampus CA1) and disease states (AD, Parkinson's disease, and normal controls). Two novel pathways, neuronal repair and viral processing were identified. Based on KEGG analysis (Kyoto Encyclopedia of Genes and Genomes, a collection of biological pathways), one of the most significant viruses was hepatitis B virus (HBV) (false discovery rate < 0.00000001). Immunohistochemical analysis using HBV-core antibody in HBV-positive control, amnestic mild cognitive impairment, and HBV-positive AD cases show increased HBV immunoreactivity as disease pathology increases. These results are the first, to our knowledge, to show regional differences in human microglia. In addition, these data reveal new functions for microglia and suggest a novel risk factor for AD.

Integrated genomic analyses reveal frequent TERT aberrations in acral melanoma.

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver BRAF/NRAS/NF1 mutations. In contrast with CM, we observed PAK1 copy gains in 15% of patients, and somatic TERT translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of TERT in ALM tumorigenesis and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.

Clinical Implementation of Integrated Genomic Profiling in Patients with Advanced Cancers.

DNA focused panel sequencing has been rapidly adopted to assess therapeutic targets in advanced/refractory cancer. Integrated Genomic Profiling (IGP) utilising DNA/RNA with tumour/normal comparisons in a Clinical Laboratory Improvement Amendments (CLIA) compliant setting enables a single assay to provide: therapeutic target prioritisation, novel target discovery/application and comprehensive germline assessment. A prospective study in 35 advanced/refractory cancer patients was conducted using CLIA-compliant IGP. Feasibility was assessed by estimating time to results (TTR), prioritising/assigning putative therapeutic targets, assessing drug access, ascertaining germline alterations, and assessing patient preferences/perspectives on data use/reporting. Therapeutic targets were identified using biointelligence/pathway analyses and interpreted by a Genomic Tumour Board. Seventy-five percent of cases harboured 1-3 therapeutically targetable mutations/case (median 79 mutations of potential functional significance/case). Median time to CLIA-validated results was 116 days with CLIA-validation of targets achieved in 21/22 patients. IGP directed treatment was instituted in 13 patients utilising on/off label FDA approved drugs (n = 9), clinical trials (n = 3) and single patient IND (n = 1). Preliminary clinical efficacy was noted in five patients (two partial response, three stable disease). Although barriers to broader application exist, including the need for wider availability of therapies, IGP in a CLIA-framework is feasible and valuable in selection/prioritisation of anti-cancer therapeutic targets.

Next-generation profiling to identify the molecular etiology of Parkinson dementia.

OBJECTIVE: We sought to determine the underlying cortical gene expression changes associated with Parkinson dementia using a next-generation RNA sequencing approach. METHODS: In this study, we used RNA sequencing to evaluate differential gene expression and alternative splicing in the posterior cingulate cortex from neurologically normal control patients, patients with Parkinson disease, and patients with Parkinson disease with dementia. RESULTS: Genes overexpressed in both disease states were involved with an immune response, whereas shared underexpressed genes functioned in signal transduction or as components of the cytoskeleton. Alternative splicing analysis produced a pattern of immune and RNA-processing disturbances. CONCLUSIONS: Genes with the greatest degree of differential expression did not overlap with genes exhibiting significant alternative splicing activity. Such variation indicates the importance of broadening expression studies to include exon-level changes because there can be significant differential splicing activity with potential structural consequences, a subtlety that is not detected when examining differential gene expression alone, or is underrepresented with probe-limited array technology.

Evidence Suggesting That Discontinuous Dosing of ALK Kinase Inhibitors May Prolong Control of ALK+ Tumors.

The anaplastic lymphoma kinase (ALK) is chromosomally rearranged in a subset of certain cancers, including 2% to 7% of non-small cell lung cancers (NSCLC) and ∼70% of anaplastic large cell lymphomas (ALCL). The ALK kinase inhibitors crizotinib and ceritinib are approved for relapsed ALK(+) NSCLC, but acquired resistance to these drugs limits median progression-free survival on average to ∼10 months. Kinase domain mutations are detectable in 25% to 37% of resistant NSCLC samples, with activation of bypass signaling pathways detected frequently with or without concurrent ALK mutations. Here we report that, in contrast to NSCLC cells, drug-resistant ALCL cells show no evidence of bypassing ALK by activating alternate signaling pathways. Instead, drug resistance selected in this setting reflects upregulation of ALK itself. Notably, in the absence of crizotinib or ceritinib, we found that increased ALK signaling rapidly arrested or killed cells, allowing a prolonged control of drug-resistant tumors in vivo with the administration of discontinuous rather than continuous regimens of drug dosing. Furthermore, even when drug resistance mutations were detected in the kinase domain, overexpression of the mutant ALK was toxic to tumor cells. We confirmed these findings derived from human ALCL cells in murine pro-B cells that were transformed to cytokine independence by ectopic expression of an activated NPM-ALK fusion oncoprotein. In summary, our results show how ALK activation functions as a double-edged sword for tumor cell viability, with potential therapeutic implications.

Personalized treatment of Sézary syndrome by targeting a novel CTLA4:CD28 fusion.

Matching molecularly targeted therapies with cancer subtype-specific gene mutations is revolutionizing oncology care. However, for rare cancers this approach is problematic due to the often poor understanding of the disease's natural history and phenotypic heterogeneity, making treatment of these cancers a particularly unmet medical need in clinical oncology. Advanced Sézary syndrome (SS), an aggressive, exceedingly rare variant of cutaneous T-cell lymphoma (CTCL) is a prototypical example of a rare cancer. Through whole genome and RNA sequencing (RNA-seq) of a SS patient's tumor we discovered a highly expressed gene fusion between CTLA4 (cytotoxic T lymphocyte antigen 4) and CD28 (cluster of differentiation 28), predicting a novel stimulatory molecule on the surface of tumor T cells. Treatment with the CTLA4 inhibitor ipilimumab resulted in a rapid clinical response. Our findings suggest a novel driver mechanism for SS, and cancer in general, and exemplify an emerging model of cancer treatment using exploratory genomic analysis to identify a personally targeted treatment option when conventional therapies are exhausted.

Alzheimer's disease is associated with altered expression of genes involved in immune response and mitochondrial processes in astrocytes.

Alzheimer's disease (AD) is characterized by deficits in cerebral metabolic rates of glucose in the posterior cingulate (PC) and precuneus in AD subjects, and in APOEε4 carriers, decades before the onset of measureable cognitive deficits. However, the cellular and molecular basis of this phenotype remains to be clarified. Given the roles of astrocytes in energy storage and brain immunity, we sought to characterize the transcriptome of AD PC astrocytes. Cells were laser capture microdissected from AD (n = 10) and healthy elderly control (n = 10) subjects for RNA sequencing. We generated >5.22 billion reads and compared sequencing data between controls and AD patients. We identified differentially expressed mitochondria-related genes including TRMT61B, FASTKD2, and NDUFA4L2, and using pathway and weighted gene coexpression analyses, we identified differentially expressed immune response genes. A number of these genes, including CLU, C3, and CD74, have been implicated in beta amyloid generation or clearance. These data provide key insights into astrocyte-specific contributions to AD, and we present this data set as a publicly available resource.