Identification of liver cancer progenitors whose malignant progression depends on autocrine IL-6 signaling.

Hepatocellular carcinoma (HCC) is a slowly developing malignancy postulated to evolve from premalignant lesions in chronically damaged livers. However, it was never established that premalignant lesions actually contain tumor progenitors that give rise to cancer. Here, we describe isolation and characterization of HCC progenitor cells (HcPCs) from different mouse HCC models. Unlike fully malignant HCC, HcPCs give rise to cancer only when introduced into a liver undergoing chronic damage and compensatory proliferation. Although HcPCs exhibit a similar transcriptomic profile to bipotential hepatobiliary progenitors, the latter do not give rise to tumors. Cells resembling HcPCs reside within dysplastic lesions that appear several months before HCC nodules. Unlike early hepatocarcinogenesis, which depends on paracrine IL-6 production by inflammatory cells, due to upregulation of LIN28 expression, HcPCs had acquired autocrine IL-6 signaling that stimulates their in vivo growth and malignant progression. This may be a general mechanism that drives other IL-6-producing malignancies.

Systematic analysis of binding of transcription factors to noncoding variants.

Many sequence variants have been linked to complex human traits and diseases1, but deciphering their biological functions remains challenging, as most of them reside in noncoding DNA. Here we have systematically assessed the binding of 270 human transcription factors to 95,886 noncoding variants in the human genome using an ultra-high-throughput multiplex protein-DNA binding assay, termed single-nucleotide polymorphism evaluation by systematic evolution of ligands by exponential enrichment (SNP-SELEX). The resulting 828 million measurements of transcription factor-DNA interactions enable estimation of the relative affinity of these transcription factors to each variant in vitro and evaluation of the current methods to predict the effects of noncoding variants on transcription factor binding. We show that the position weight matrices of most transcription factors lack sufficient predictive power, whereas the support vector machine combined with the gapped k-mer representation show much improved performance, when assessed on results from independent SNP-SELEX experiments involving a new set of 61,020 sequence variants. We report highly predictive models for 94 human transcription factors and demonstrate their utility in genome-wide association studies and understanding of the molecular pathways involved in diverse human traits and diseases.

Enhancer release and retargeting activates disease-susceptibility genes.

The functional engagement between an enhancer and its target promoter ensures precise gene transcription1. Understanding the basis of promoter choice by enhancers has important implications for health and disease. Here we report that functional loss of a preferred promoter can release its partner enhancer to loop to and activate an alternative promoter (or alternative promoters) in the neighbourhood. We refer to this target-switching process as 'enhancer release and retargeting'. Genetic deletion, motif perturbation or mutation, and dCas9-mediated CTCF tethering reveal that promoter choice by an enhancer can be determined by the binding of CTCF at promoters, in a cohesin-dependent manner-consistent with a model of 'enhancer scanning' inside the contact domain. Promoter-associated CTCF shows a lower affinity than that at chromatin domain boundaries and often lacks a preferred motif orientation or a partnering CTCF at the cognate enhancer, suggesting properties distinct from boundary CTCF. Analyses of cancer mutations, data from the GTEx project and risk loci from genome-wide association studies, together with a focused CRISPR interference screen, reveal that enhancer release and retargeting represents an overlooked mechanism that underlies the activation of disease-susceptibility genes, as exemplified by a risk locus for Parkinson's disease (NUCKS1-RAB7L1) and three loci associated with cancer (CLPTM1L-TERT, ZCCHC7-PAX5 and PVT1-MYC).

Glioblastoma cellular cross-talk converges on NF-κB to attenuate EGFR inhibitor sensitivity.

In glioblastoma (GBM), heterogeneous expression of amplified and mutated epidermal growth factor receptor (EGFR) presents a substantial challenge for the effective use of EGFR-directed therapeutics. Here we demonstrate that heterogeneous expression of the wild-type receptor and its constitutively active mutant form, EGFRvIII, limits sensitivity to these therapies through an interclonal communication mechanism mediated by interleukin-6 (IL-6) cytokine secreted from EGFRvIII-positive tumor cells. IL-6 activates a NF-κB signaling axis in a paracrine and autocrine manner, leading to bromodomain protein 4 (BRD4)-dependent expression of the prosurvival protein survivin (BIRC5) and attenuation of sensitivity to EGFR tyrosine kinase inhibitors (TKIs). NF-κB and survivin are coordinately up-regulated in GBM patient tumors, and functional inhibition of either protein or BRD4 in in vitro and in vivo models restores sensitivity to EGFR TKIs. These results provide a rationale for improving anti-EGFR therapeutic efficacy through pharmacological uncoupling of a convergence point of NF-κB-mediated survival that is leveraged by an interclonal circuitry mechanism established by intratumoral mutational heterogeneity.

Deciphering the genetic architecture and ethnographic distribution of IRD in three ethnic populations by whole genome sequence analysis.

Patients with inherited retinal dystrophies (IRDs) were recruited from two understudied populations: Mexico and Pakistan as well as a third well-studied population of European Americans to define the genetic architecture of IRD by performing whole-genome sequencing (WGS). Whole-genome analysis was performed on 409 individuals from 108 unrelated pedigrees with IRDs. All patients underwent an ophthalmic evaluation to establish the retinal phenotype. Although the 108 pedigrees in this study had previously been examined for mutations in known IRD genes using a wide range of methodologies including targeted gene(s) or mutation(s) screening, linkage analysis and exome sequencing, the gene mutations responsible for IRD in these 108 pedigrees were not determined. WGS was performed on these pedigrees using Illumina X10 at a minimum of 30X depth. The sequence reads were mapped against hg19 followed by variant calling using GATK. The genome variants were annotated using SnpEff, PolyPhen2, and CADD score; the structural variants (SVs) were called using GenomeSTRiP and LUMPY. We identified potential causative sequence alterations in 61 pedigrees (57%), including 39 novel and 54 reported variants in IRD genes. For 57 of these pedigrees the observed genotype was consistent with the initial clinical diagnosis, the remaining 4 had the clinical diagnosis reclassified based on our findings. In seven pedigrees (12%) we observed atypical causal variants, i.e. unexpected genotype(s), including 4 pedigrees with causal variants in more than one IRD gene within all affected family members, one pedigree with intrafamilial genetic heterogeneity (different affected family members carrying causal variants in different IRD genes), one pedigree carrying a dominant causative variant present in pseudo-recessive form due to consanguinity and one pedigree with a de-novo variant in the affected family member. Combined atypical and large structural variants contributed to about 20% of cases. Among the novel mutations, 75% were detected in Mexican and 50% found in European American pedigrees and have not been reported in any other population while only 20% were detected in Pakistani pedigrees and were not previously reported. The remaining novel IRD causative variants were listed in gnomAD but were found to be very rare and population specific. Mutations in known IRD associated genes contributed to pathology in 63% Mexican, 60% Pakistani and 45% European American pedigrees analyzed. Overall, contribution of known IRD gene variants to disease pathology in these three populations was similar to that observed in other populations worldwide. This study revealed a spectrum of mutations contributing to IRD in three populations, identified a large proportion of novel potentially causative variants that are specific to the corresponding population or not reported in gnomAD and shed light on the genetic architecture of IRD in these diverse global populations.

In heart failure reactivation of RNA-binding proteins is associated with the expression of 1,523 fetal-specific isoforms.

Reactivation of fetal-specific genes and isoforms occurs during heart failure. However, the underlying molecular mechanisms and the extent to which the fetal program switch occurs remains unclear. Limitations hindering transcriptome-wide analyses of alternative splicing differences (i.e. isoform switching) in cardiovascular system (CVS) tissues between fetal, healthy adult and heart failure have included both cellular heterogeneity across bulk RNA-seq samples and limited availability of fetal tissue for research. To overcome these limitations, we have deconvoluted the cellular compositions of 996 RNA-seq samples representing heart failure, healthy adult (heart and arteria), and fetal-like (iPSC-derived cardiovascular progenitor cells) CVS tissues. Comparison of the expression profiles revealed that reactivation of fetal-specific RNA-binding proteins (RBPs), and the accompanied re-expression of 1,523 fetal-specific isoforms, contribute to the transcriptome differences between heart failure and healthy adult heart. Of note, isoforms for 20 different RBPs were among those that reverted in heart failure to the fetal-like expression pattern. We determined that, compared with adult-specific isoforms, fetal-specific isoforms encode proteins that tend to have more functions, are more likely to harbor RBP binding sites, have canonical sequences at their splice sites, and contain typical upstream polypyrimidine tracts. Our study suggests that compared with healthy adult, fetal cardiac tissue requires stricter transcriptional regulation, and that during heart failure reversion to this stricter transcriptional regulation occurs. Furthermore, we provide a resource of cardiac developmental stage-specific and heart failure-associated genes and isoforms, which are largely unexplored and can be exploited to investigate novel therapeutics for heart failure.

Ultra-Sharp Nanowire Arrays Natively Permeate, Record, and Stimulate Intracellular Activity in Neuronal and Cardiac Networks.

We report innovative scalable, vertical, ultra-sharp nanowire arrays that are individually addressable to enable long-term, native recordings of intracellular potentials. Stable amplitudes of intracellular potentials from 3D tissue-like networks of neurons and cardiomyocytes are obtained. Individual electrical addressability is necessary for high-fidelity intracellular electrophysiological recordings. This study paves the way toward predictive, high-throughput, and low-cost electrophysiological drug screening platforms.

Systems genomics in age-related macular degeneration.

Genomic studies in age-related macular degeneration (AMD) have identified genetic variants that account for the majority of AMD risk. An important next step is to understand the functional consequences and downstream effects of the identified AMD-associated genetic variants. Instrumental for this next step are 'omics' technologies, which enable high-throughput characterization and quantification of biological molecules, and subsequent integration of genomics with these omics datasets, a field referred to as systems genomics. Single cell sequencing studies of the retina and choroid demonstrated that the majority of candidate AMD genes identified through genomic studies are expressed in non-neuronal cells, such as the retinal pigment epithelium (RPE), glia, myeloid and choroidal cells, highlighting that many different retinal and choroidal cell types contribute to the pathogenesis of AMD. Expression quantitative trait locus (eQTL) studies in retinal tissue have identified putative causal genes by demonstrating a genetic overlap between gene regulation and AMD risk. Linking genetic data to complement measurements in the systemic circulation has aided in understanding the effect of AMD-associated genetic variants in the complement system, and supports that protein QTL (pQTL) studies in plasma or serum samples may aid in understanding the effect of genetic variants and pinpointing causal genes in AMD. A recent epigenomic study fine-mapped AMD causal variants by determing regulatory regions in RPE cells differentiated from induced pluripotent stem cells (iPSC-RPE). Another approach that is being employed to pinpoint causal AMD genes is to produce synthetic DNA assemblons representing risk and protective haplotypes, which are then delivered to cellular or animal model systems. Pinpointing causal genes and understanding disease mechanisms is crucial for the next step towards clinical translation. Clinical trials targeting proteins encoded by the AMD-associated genomic loci C3, CFB, CFI, CFH, and ARMS2/HTRA1 are currently ongoing, and a phase III clinical trial for C3 inhibition recently showed a modest reduction of lesion growth in geographic atrophy. The EYERISK consortium recently developed a genetic test for AMD that allows genotyping of common and rare variants in AMD-associated genes. Polygenic risk scores (PRS) were applied to quantify AMD genetic risk, and may aid in predicting AMD progression. In conclusion, genomic studies represent a turning point in our exploration of AMD. The results of those studies now serve as a driving force for several clinical trials. Expanding to omics and systems genomics will further decipher function and causality from the associations that have been reported, and will enable the development of therapies that will lessen the burden of AMD.

PI3Kγ is a molecular switch that controls immune suppression.

Macrophages play critical, but opposite, roles in acute and chronic inflammation and cancer. In response to pathogens or injury, inflammatory macrophages express cytokines that stimulate cytotoxic T cells, whereas macrophages in neoplastic and parasitic diseases express anti-inflammatory cytokines that induce immune suppression and may promote resistance to T cell checkpoint inhibitors. Here we show that macrophage PI 3-kinase γ controls a critical switch between immune stimulation and suppression during inflammation and cancer. PI3Kγ signalling through Akt and mTor inhibits NFκB activation while stimulating C/EBPß activation, thereby inducing a transcriptional program that promotes immune suppression during inflammation and tumour growth. By contrast, selective inactivation of macrophage PI3Kγ stimulates and prolongs NFκB activation and inhibits C/EBPß activation, thus promoting an immunostimulatory transcriptional program that restores CD8+ T cell activation and cytotoxicity. PI3Kγ synergizes with checkpoint inhibitor therapy to promote tumour regression and increased survival in mouse models of cancer. In addition, PI3Kγ-directed, anti-inflammatory gene expression can predict survival probability in cancer patients. Our work thus demonstrates that therapeutic targeting of intracellular signalling pathways that regulate the switch between macrophage polarization states can control immune suppression in cancer and other disorders.

Genomics Links Inflammation With Neurocognitive Impairment in Children Living With Human Immunodeficiency Virus Type-1.

BACKGROUND: We identified host single-nucleotide variants (SNVs) associated with neurocognitive impairment (NCI) in perinatally HIV-infected (PHIV) children. METHODS: Whole-exome sequencing (WES) was performed on 217 PHIV with cognitive score for age (CSA) < 70 and 247 CSA ≥ 70 (discovery cohort [DC]). SNVs identified in DC were evaluated in 2 validation cohorts (VC). Logistic regression was used to estimate adjusted odds ratios (ORs) for NCI. A human microglia NLRP3 inflammasome assay characterized the role of identified genes. RESULTS: Twenty-nine SNVs in 24 genes reaching P ≤ .002 and OR ≥ 1.5 comparing CSA < 70 to CSA ≥ 70 were identified in the DC, of which 3 SNVs were identified in VCs for further study. Combining the 3 cohorts, SNV in CCRL2 (rs3204849) was associated with decreased odds of NCI (P < .0001); RETREG1/FAM134B (rs61733811) and YWHAH (rs73884247) were associated with increased risk of NCI (P < .0001 and P < .001, respectively). Knockdown of CCRL2 led to decreased microglial release of IL-1ß following exposure to ssRNA40 while knockdown of RETREG1 and YWHAH resulted in increased IL-1ß release. CONCLUSIONS: Using WES and 2 VCs, and gene silencing of microglia we identified 3 genetic variants associated with NCI and inflammation in HIV-infected children.

Detection and validation of novel mutations in MERTK in a simplex case of retinal degeneration using WGS and hiPSC-RPEs model.

Inherited retinal degenerations (IRDs) are a group of genetically heterogeneous conditions with a broad phenotypic heterogeneity. Here, we report detection and validation of the underlying cause of progressive retinal degeneration in a nuclear family of European descent with a single affected individual. Whole genome sequencing of the proband and her unaffected sibling identified a novel intron 8 donor splice site variant (c.1296 + 1G>A) and a novel 731 base pair deletion encompassing exon 9 (Chr2:g.112751488_112752218 del) resulting in c.1297_1451del; p.K433_G484fsTer3 in the Mer tyrosine kinase protooncogene (MERTK), which is highly expressed in the retinal pigment epithelium (RPE). The proband carried both variants in the heterozygous state, which segregated with disease in the pedigree. These MERTK variants are predicted to result in the defective splicing of exon 8 and loss of exon 9 respectively. To evaluate the impact of these novel variants, peripheral blood mononuclear cells of the proband and her parents were reprogrammed to humaninduced pluripotent stem cell (hiPSC) lines, which were subsequently differentiated to hiPSC-RPE. Analysis of the proband's hiPSC-RPE revealed the absence of both MERTK transcript and its respective protein as well as abnormal phagocytosis when compared with the parental hiPSC-RPE. In summary, whole genome sequencing identified novel compound heterozygous variants in MERTK as the underlying cause of progressive retinal degeneration in a simplex case. Further, analysis using an hiPSC-RPE model established the functional impact of novel MERTK mutations and revealed the potential mechanism underlying pathology in the proband.

Activation of hedgehog signaling associates with early disease progression in chronic lymphocytic leukemia.

Targeted sequencing of 103 leukemia-associated genes in leukemia cells from 841 treatment-naive patients with chronic lymphocytic leukemia (CLL) identified 89 (11%) patients as having CLL cells with mutations in genes encoding proteins that putatively are involved in hedgehog (Hh) signaling. Consistent with this finding, there was a significant association between the presence of these mutations and the expression of GLI1 (χ2 test, P < .0001), reflecting activation of the Hh pathway. However, we discovered that 38% of cases without identified mutations also were GLI1+ Patients with GLI1+ CLL cells had a shorter median treatment-free survival than patients with CLL cells lacking expression of GLI1 independent of IGHV mutation status. We found that GANT61, a small molecule that can inhibit GLI1, was highly cytotoxic for GLI1+ CLL cells relative to that of CLL cells without GLI1. Collectively, this study shows that a large proportion of patients have CLL cells with activated Hh signaling, which is associated with early disease progression and enhanced sensitivity to inhibition of GLI1.

Genomic and transcriptomic association studies identify 16 novel susceptibility loci for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality. To advance our understanding of the biology contributing to VTE, we conducted a genome-wide association study (GWAS) of VTE and a transcriptome-wide association study (TWAS) based on imputed gene expression from whole blood and liver. We meta-analyzed GWAS data from 18 studies for 30 234 VTE cases and 172 122 controls and assessed the association between 12 923 718 genetic variants and VTE. We generated variant prediction scores of gene expression from whole blood and liver tissue and assessed them for association with VTE. Mendelian randomization analyses were conducted for traits genetically associated with novel VTE loci. We identified 34 independent genetic signals for VTE risk from GWAS meta-analysis, of which 14 are newly reported associations. This included 11 newly associated genetic loci (C1orf198, PLEK, OSMR-AS1, NUGGC/SCARA5, GRK5, MPHOSPH9, ARID4A, PLCG2, SMG6, EIF5A, and STX10) of which 6 replicated, and 3 new independent signals in 3 known genes. Further, TWAS identified 5 additional genetic loci with imputed gene expression levels differing between cases and controls in whole blood (SH2B3, SPSB1, RP11-747H7.3, RP4-737E23.2) and in liver (ERAP1). At some GWAS loci, we found suggestive evidence that the VTE association signal for novel and previously known regions colocalized with expression quantitative trait locus signals. Mendelian randomization analyses suggested that blood traits may contribute to the underlying risk of VTE. To conclude, we identified 16 novel susceptibility loci for VTE; for some loci, the association signals are likely mediated through gene expression of nearby genes.

A large-scale exome array analysis of venous thromboembolism.

Although recent Genome-Wide Association Studies have identified novel associations for common variants, there has been no comprehensive exome-wide search for low-frequency variants that affect the risk of venous thromboembolism (VTE). We conducted a meta-analysis of 11 studies comprising 8,332 cases and 16,087 controls of European ancestry and 382 cases and 1,476 controls of African American ancestry genotyped with the Illumina HumanExome BeadChip. We used the seqMeta package in R to conduct single variant and gene-based rare variant tests. In the single variant analysis, we limited our analysis to the 64,794 variants with at least 40 minor alleles across studies (minor allele frequency [MAF] ~0.08%). We confirmed associations with previously identified VTE loci, including ABO, F5, F11, and FGA. After adjusting for multiple testing, we observed no novel significant findings in single variant or gene-based analysis. Given our sample size, we had greater than 80% power to detect minimum odds ratios greater than 1.5 and 1.8 for a single variant with MAF of 0.01 and 0.005, respectively. Larger studies and sequence data may be needed to identify novel low-frequency and rare variants associated with VTE risk.

Fibrinogen gamma gene rs2066865 and risk of cancer-related venous thromboembolism.

Venous thromboembolism (VTE) is a frequent complication in patients with cancer. Homozygous carriers of the fibrinogen gamma gene (FGG) rs2066865 have a moderately increased risk of VTE, but the effect of the FGG variant in cancer is unknown. We aimed to investigate the effect of the FGG variant and active cancer on the risk of VTE. Cases with incident VTE (n=640) and a randomly selected age-weighted sub-cohort (n=3,734) were derived from a population-based cohort (the Tromsø study). Cox-regression was used to estimate hazard ratios (HR) with 95% confidence intervals (CI) for VTE according to categories of cancer and FGG In those without cancer, homozygosity at the FGG variant was associated with a 70% (HR 1.7, 95% CI: 1.2-2.3) increased risk of VTE compared to non-carriers. Cancer patients homozygous for the FGG variant had a two-fold (HR 2.0, 95% CI: 1.1-3.6) higher risk of VTE than cancer patients without the variant. Moreover, the six-months cumulative incidence of VTE among cancer patients was 6.4% (95% CI: 3.5-11.6) in homozygous carriers of FGG and 3.1% (95% CI: 2.3-4.7) in those without risk alleles. A synergistic effect was observed between rs2066865 and active cancer on the risk of VTE (synergy index: 1.81, 95% CI: 1.02-3.21, attributable proportion: 0.43, 95% CI: 0.11-0.74). In conclusion, homozygosity at the FGG variant and active cancer yielded a synergistic effect on the risk of VTE.

A mutation in IFT43 causes non-syndromic recessive retinal degeneration.

The aim of this work is to identify the molecular cause of autosomal recessive early onset retinal degeneration in a consanguineous pedigree. Seventeen members of a four-generation Pakistani family were recruited and underwent a detailed ophthalmic examination. Exomes of four affected and two unaffected individuals were sequenced. Variants were filtered using exomeSuite to identify rare potentially pathogenic variants in genes expressed in the retina and/or brain and consistent with the pattern of inheritance. Effect of the variant observed in the gene Intraflagellar Transport Protein 43 (IFT43) was studied by heterologous expression in mIMCD3 and MDCK cells. Expression and sub-cellular localization of IFT43 in the retina and transiently transfected cells was examined by RT-PCR, western blot analysis, and immunohistochemistry. Affected members were diagnosed with early onset non-syndromic progressive retinal degeneration and the presence of bone spicules distributed throughout the retina at younger ages while the older affected members showed severe central choroidal atrophy. Whole-exome sequencing analysis identified a novel homozygous c.100 G > A change in IFT43 segregating with retinal degeneration and not present in ethnicity-matched controls. Immunostaining showed IFT43 localized in the photoreceptors, and to the tip of the cilia in transfected mIMCD3 and MDCK cells. The cilia in mIMCD3 and MDCK cells expressing mutant IFT43 were found to be significantly shorter (P < 0.001) than cells expressing wild-type IFT43. Our studies identified a novel homozygous mutation in the ciliary protein IFT43 as the underlying cause of recessive inherited retinal degeneration. This is the first report demonstrating the involvement of IFT43 in retinal degeneration.

IFT88 mutations identified in individuals with non-syndromic recessive retinal degeneration result in abnormal ciliogenesis.

Whole genome sequencing (WGS) was performed to identify the variants responsible for inherited retinal degeneration (IRD) in a Caucasian family. Segregation analysis of selected rare variants with pathogenic potential identified a set of compound heterozygous changes p.Arg266*:c.796C>T and p.Ala568Thr:c.1702G>A in the intraflagellar transport protein-88 (IFT88) gene segregating with IRD. Expression of IFT88 with the p.Arg266* and p.Ala568Thr mutations in mIMDC3 cells by transient transfection and in HeLa cells by introducing the mutations using CRISPR-cas9 system suggested that both mutations result in the formation of abnormal ciliary structures. The introduction of the IFT88 p.Arg266* variant in the homozygous state in HeLa cells by CRISPR-Cas9 genome-editing revealed that the mutant transcript undergoes nonsense-mediated decay leading to a significant depletion of IFT88 transcript. Additionally, abnormal ciliogenesis was observed in these cells. These observations suggest that the rare and unique combination of IFT88 alleles observed in this study provide insight into the physiological role of IFT88 in humans and the likely mechanism underlying retinal pathology in the pedigree with IRD.

High-level ROR1 associates with accelerated disease progression in chronic lymphocytic leukemia.

ROR1 is an oncoembryonic orphan receptor found on chronic lymphocytic leukemia (CLL) B cells, but not on normal postpartum tissues. ROR1 is a receptor for Wnt5a that may complex with TCL1, a coactivator of AKT that is able to promote development of CLL. We found the CLL cells of a few patients expressed negligible ROR1 (ROR1Neg), but expressed TCL1A at levels comparable to those of samples that expressed ROR1 (ROR1Pos). Transcriptome analyses revealed that ROR1Neg cases generally could be distinguished from those that were ROR1Pos in unsupervised gene-expression clustering analysis. Gene-set enrichment analyses demonstrated that ROR1Neg CLL had lower expression and activation of AKT signaling pathways relative to ROR1Pos CLL, similar to what was noted for leukemia that respectively developed in TCL1 vs ROR1xTCL1 transgenic mice. In contrast to its effect on ROR1Pos CLL, Wnt5a did not enhance the proliferation, chemotaxis, or survival of ROR1Neg CLL. We examined the CLL cells from 1568 patients, which we randomly assigned to a training or validation set of 797 or 771 cases, respectively. Using recursive partitioning, we defined a threshold for ROR1 surface expression that could segregate samples of the training set into ROR1-Hi vs ROR1-Lo subgroups that differed significantly in their median treatment-free survival (TFS). Using this threshold, we found that ROR1-Hi cases had a significantly shorter median TFS and overall survival than ROR1-Lo cases in the validation set. These data demonstrate that expression of ROR1 may promote leukemia-cell activation and survival and enhance disease progression in patients with CLL.

Systematic transcriptome analysis reveals tumor-specific isoforms for ovarian cancer diagnosis and therapy.

Tumor-specific molecules are needed across diverse areas of oncology for use in early detection, diagnosis, prognosis and therapy. Large and growing public databases of transcriptome sequencing data (RNA-seq) derived from tumors and normal tissues hold the potential of yielding tumor-specific molecules, but because the data are new they have not been fully explored for this purpose. We have developed custom bioinformatic algorithms and used them with 296 high-grade serous ovarian (HGS-OvCa) tumor and 1,839 normal RNA-seq datasets to identify mRNA isoforms with tumor-specific expression. We rank prioritized isoforms by likelihood of being expressed in HGS-OvCa tumors and not in normal tissues and analyzed 671 top-ranked isoforms by high-throughput RT-qPCR. Six of these isoforms were expressed in a majority of the 12 tumors examined but not in 18 normal tissues. An additional 11 were expressed in most tumors and only one normal tissue, which in most cases was fallopian or colon. Of the 671 isoforms, the topmost 5% (n = 33) ranked based on having tumor-specific or highly restricted normal tissue expression by RT-qPCR analysis are enriched for oncogenic, stem cell/cancer stem cell, and early development loci--including ETV4, FOXM1, LSR, CD9, RAB11FIP4, and FGFRL1. Many of the 33 isoforms are predicted to encode proteins with unique amino acid sequences, which would allow them to be specifically targeted for one or more therapeutic strategies--including monoclonal antibodies and T-cell-based vaccines. The systematic process described herein is readily and rapidly applicable to the more than 30 additional tumor types for which sufficient amounts of RNA-seq already exist.