ENTRAIN: integrating trajectory inference and gene regulatory networks with spatial data to co-localize the receptor-ligand interactions that specify cell fate.

MOTIVATION: Cell fate is commonly studied by profiling the gene expression of single cells to infer developmental trajectories based on expression similarity, RNA velocity, or statistical mechanical properties. However, current approaches do not recover microenvironmental signals from the cellular niche that drive a differentiation trajectory. RESULTS: We resolve this with environment-aware trajectory inference (ENTRAIN), a computational method that integrates trajectory inference methods with ligand-receptor pair gene regulatory networks to identify extracellular signals and evaluate their relative contribution towards a differentiation trajectory. The output from ENTRAIN can be superimposed on spatial data to co-localize cells and molecules in space and time to map cell fate potentials to cell-cell interactions. We validate and benchmark our approach on single-cell bone marrow and spatially resolved embryonic neurogenesis datasets to identify known and novel environmental drivers of cellular differentiation. AVAILABILITY AND IMPLEMENTATION: ENTRAIN is available as a public package at https://github.com/theimagelab/entrain and can be used on both single-cell and spatially resolved datasets.

Circulating Tumour DNA Biomarkers Associated with Outcomes in Metastatic Prostate Cancer Treated with Lutetium-177-PSMA-617.

Background: Lutetium-177-prostate-specific membrane antigen- 617 (Lu-PSMA) is an effective therapy for metastatic castration-resistant prostate cancer (mCRPC). However, treatment responses are heterogeneous despite stringent positron emission tomography (PET)-based imaging selection criteria. Molecularly based biomarkers have potential to refine patient selection and optimise outcomes. Objective: To identify circulating tumour DNA (ctDNA) features associated with treatment outcomes for men treated with Lu-PSMA. Design setting and participants: ctDNA from men treated with Lu-PSMA in combination with idronoxil for progressive mCRPC were analysed using an 85-gene customised sequencing assay. ctDNA fractions, molecular profiles, and the presence of alterations in aggressive-variant prostate cancer (AVPC) genes were analysed at baseline, cycle 3 and at disease progression. Intervention: Men received Lu-PSMA with idronoxil every 6 wk for up to six cycles. Outcome measurements and statistical analysis: Baseline and exit PSMA and fluorodeoxyglucose PET/computed tomography (CT) imaging was conducted at baseline and study exit. Single-photon emission CT (SPECT) scans were performed 24 h after Lu-PSMA. Blood samples were collected at baseline,cycle 3 and at disease progression. Cox proportional-hazards models were used to assess associations and derive hazard ratios (HRs) and confidence intervals (CIs) for associations between molecular factors, imaging features, and clinical outcomes. Results and limitations: Sixty samples from 32 men were sequenced (32 at baseline, 24 at cycle 3, four from patients with disease progression); two samples (baseline, on-treatment) from one individual were excluded from analysis owing to poor quality of the baseline sequencing data. Alterations in AVPC genes were associated with shorter prostate-specific antigen (PSA) progression-free survival (PFS) and overall survival (OS) in univariate (HR 3.4, 95% CI 1.5-7.7; p = 0.0036; and HR 3.3, 95% CI 1.4-7.7; p = 0.0063, respectively) and multivariate analyses (HR 4.8, 95% CI 1.8-13; p = 0.0014; and HR 4.1, 95% CI 1.6-11; p = 0.004). Conclusions: ctDNA alterations in AVPC genes were associated with shorter PSA PFS and OS among men treated with Lu-PSMA and intermittent idronoxil. These candidate molecular biomarkers warrant further study to determine whether they have predictive value and potential to guide synergistic combination strategies to enhance outcomes for men treated with Lu-PSMA for mCRPC. Patient summary: Certain DNA/gene changes detected in the blood of men with advanced prostate cancer were associated with shorter benefit from lutetium PSMA, a targeted radioactive therapy. This information may be useful in determining which men may benefit most from this treatment, but additional research is needed.

Neutrophil Conversion to a Tumor-Killing Phenotype Underpins Effective Microbial Therapy.

The inflammatory microenvironment of solid tumors creates a protumorigenic milieu that resembles chronic inflammation akin to a subverted wound healing response. Here, we investigated the effect of converting the tumor microenvironment from a chronically inflamed state to one of acute microbial inflammation by injecting microbial bioparticles directly into tumors. Intratumoral microbial bioparticle injection led to rapid and dramatic changes in the tumor immune composition, the most striking of which was a substantial increase in the presence of activated neutrophils. In situ photoconversion and intravital microscopy indicated that tumor neutrophils transiently switched from sessile producers of VEGF to highly motile neutrophils that clustered to make neutrophil-rich domains in the tumor. The neutrophil clusters remodeled tumor tissue and repressed tumor growth. Single-cell transcriptional analysis of microbe-stimulated neutrophils showed a profound shift in gene expression towards heightened activation and antimicrobial effector function. Microbe-activated neutrophils also upregulated chemokines known to regulate neutrophil and CD8+ T-cell recruitment. Microbial therapy also boosted CD8+ T-cell function and enhanced the therapeutic benefit of checkpoint inhibitor therapy in tumor-bearing mice and provided protection in a model of tumor recurrence. These data indicate that one of the major effector mechanisms of microbial therapy is the conversion of tumor neutrophils from a wound healing to an acutely activated cytotoxic phenotype, highlighting a rationale for broader deployment of microbial therapy in the treatment of solid cancers. SIGNIFICANCE: Intratumoral injection of microbial bioparticles stimulates neutrophil antitumor functions, suggesting pathways for optimizing efficacy of microbial therapies and paving the way for their broader utilization in the clinic.

Quantitative ctDNA Detection in Hepatoblastoma: Implications for Precision Medicine.

Hepatoblastoma is characterized by driver mutations in CTNNB1, making it an attractive biomarker for a liquid biopsy approach utilizing circulating tumor DNA (ctDNA). This prospective observational study sought to ascertain the feasibility of ctDNA detection in patients with hepatoblastoma and explore its associations with established clinical indicators and biomarkers, including serum Alpha-fetoprotein (AFP). We obtained 38 plasma samples and 17 tumor samples from 20 patients with hepatoblastoma. These samples were collected at various stages: 10 at initial diagnosis, 17 during neoadjuvant chemotherapy, 6 post-operatively, and 5 at disease recurrence. Utilizing a bespoke sequencing assay we developed called QUENCH, we identified single nucleotide variants and deletions in CTNNB1 ctDNA. Our study demonstrated the capability to quantitate ctDNA down to a variant allele frequency of 0.3%, achieving a sensitivity of 90% for patients at initial diagnosis, and a specificity of 100% at the patient level. Notably, ctDNA positivity correlated with tumor burden, and ctDNA levels exhibited associations with macroscopic residual disease and treatment response. Our findings provide evidence for the utility of quantitative ctDNA detection in hepatoblastoma management. Given the distinct detection targets, ctDNA and AFP-based stratification and monitoring approaches could synergize to enhance clinical decision-making. Further research is needed to elucidate the interplay between ctDNA and AFP and determine the optimal clinical applications for both methods in risk stratification and residual disease detection.

Embryonic lethality and defective mammary gland development of activator-function impaired conditional knock-in Erbb3 V943R mice.

ERBB3 is a pseudokinase domain-containing member of the ERBB family of receptor tyrosine kinases (RTKs). Following ligand binding, ERBB receptors homo- or hetero-dimerize, leading to a head-to-tail arrangement of the intracellular kinase domains, where the "receiver" kinase domain of one ERBB is activated by the "activator" domain of the other ERBB in the dimer. In ERBB3, a conserved valine at codon 943 (V943) in the kinase C-terminal domain has been shown to be important for its function as an "activator" kinase in vitro. Here we report a knock-in mouse model where we have modified the endogenous Erbb3 allele to allow for tissue-specific conditional expression of Erbb3 V943R (Erbb3 CKI-V943R ). Additionally, we generated an Erbb3 D850N (Erbb3 CKI-D850N ) conditional knock-in mouse model where the conserved aspartate in the DFG motif of the pseudokinase domain was mutated to abolish any potential residual kinase activity. While Erbb3 D850N/D850N animals developed normally, homozygous Erbb3 V943R/V943R expression during development resulted in embryonic lethality. Further, tissue specific expression of Erbb3 V943R/V943R in the mammary gland epithelium following its activation using MMTV-Cre resulted in delayed elongation of the ductal network during puberty. Single-cell RNA-seq analysis of Erbb3 V943R/V943R mammary glands showed a reduction in a specific subset of fibrinogen-producing luminal epithelial cells.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

Inducers of the endothelial cell barrier identified through chemogenomic screening in genome-edited hPSC-endothelial cells.

The blood-retina barrier and blood-brain barrier (BRB/BBB) are selective and semipermeable and are critical for supporting and protecting central nervous system (CNS)-resident cells. Endothelial cells (ECs) within the BRB/BBB are tightly coupled, express high levels of Claudin-5 (CLDN5), a junctional protein that stabilizes ECs, and are important for proper neuronal function. To identify novel CLDN5 regulators (and ultimately EC stabilizers), we generated a CLDN5-P2A-GFP stable cell line from human pluripotent stem cells (hPSCs), directed their differentiation to ECs (CLDN5-GFP hPSC-ECs), and performed flow cytometry-based chemogenomic library screening to measure GFP expression as a surrogate reporter of barrier integrity. Using this approach, we identified 62 unique compounds that activated CLDN5-GFP. Among them were TGF-ß pathway inhibitors, including RepSox. When applied to hPSC-ECs, primary brain ECs, and retinal ECs, RepSox strongly elevated barrier resistance (transendothelial electrical resistance), reduced paracellular permeability (fluorescein isothiocyanate-dextran), and prevented vascular endothelial growth factor A (VEGFA)-induced barrier breakdown in vitro. RepSox also altered vascular patterning in the mouse retina during development when delivered exogenously. To determine the mechanism of action of RepSox, we performed kinome-, transcriptome-, and proteome-profiling and discovered that RepSox inhibited TGF-ß, VEGFA, and inflammatory gene networks. In addition, RepSox not only activated vascular-stabilizing and barrier-establishing Notch and Wnt pathways, but also induced expression of important tight junctions and transporters. Taken together, our data suggest that inhibiting multiple pathways by selected individual small molecules, such as RepSox, may be an effective strategy for the development of better BRB/BBB models and novel EC barrier-inducing therapeutics.

Heterogeneity of Satellite Cells Implicates DELTA1/NOTCH2 Signaling in Self-Renewal.

How satellite cells and their progenitors balance differentiation and self-renewal to achieve sustainable tissue regeneration is not well understood. A major roadblock to understanding satellite cell fate decisions has been the difficulty of studying this process in vivo. By visualizing expression dynamics of myogenic transcription factors during early regeneration in vivo, we identify the time point at which cells undergo decisions to differentiate or self-renew. Single-cell RNA sequencing reveals heterogeneity of satellite cells, including a subpopulation enriched in Notch2 receptor expression, during both muscle homeostasis and regeneration. Furthermore, we reveal that differentiating cells express the Dll1 ligand. Using antagonistic antibodies, we demonstrate that the DLL1 and NOTCH2 signaling pair is required for satellite cell self-renewal. Thus, differentiating cells provide the self-renewing signal during regeneration, enabling proportional regeneration in response to injury while maintaining the satellite cell pool. These findings have implications for therapeutic control of muscle regeneration.

Restricted epitope specificity determined by variable region germline segment pairing in rodent antibody repertoires.

Antibodies from B-cell clonal lineages share sequence and structural properties as well as epitope specificity. Clonally unrelated antibodies can similarly share sequence and specificity properties and are said to be convergent. Convergent antibody responses against several antigens have been described in humans and mice and include different classes of shared sequence features. In particular, some antigens and epitopes can induce convergent responses of clonally unrelated antibodies with restricted heavy (VH) and light (VL) chain variable region germline segment usage without similarity in the heavy chain third complementarity-determining region (CDR H3), a critical specificity determinant. Whether these V germline segment-restricted responses reflect a general epitope specificity restriction of antibodies with shared VH/VL pairing is not known. Here, we investigated this question by determining patterns of antigen binding competition between clonally unrelated antigen-specific rat antibodies from paired-chain deep sequencing datasets selected based solely on VH/VL pairing. We found that antibodies with shared VH/VL germline segment pairings but divergent CDR H3 sequences almost invariably have restricted epitope specificity indicated by shared binding competition patterns. This epitope restriction included 82 of 85 clonally unrelated antibodies with 13 different VH/VL pairings binding in 8 epitope groups in 2 antigens. The corollary that antibodies with shared VH/VL pairing and epitope-restricted binding can accommodate widely divergent CDR H3 sequences was confirmed by in vitro selection of variants of anti-human epidermal growth factor receptor 2 antibodies known to mediate critical antigen interactions through CDR H3. Our results show that restricted epitope specificity determined by VH/VL germline segment pairing is a general property of rodent antigen-specific antibodies.

Peripheral T cell expansion predicts tumour infiltration and clinical response.

Despite the resounding clinical success in cancer treatment of antibodies that block the interaction of PD1 with its ligand PDL11, the mechanisms involved remain unknown. A major limitation to understanding the origin and fate of T cells in tumour immunity is the lack of quantitative information on the distribution of individual clonotypes of T cells in patients with cancer. Here, by performing deep single-cell sequencing of RNA and T cell receptors in patients with different types of cancer, we survey the profiles of various populations of T cells and T cell receptors in tumours, normal adjacent tissue, and peripheral blood. We find clear evidence of clonotypic expansion of effector-like T cells not only within the tumour but also in normal adjacent tissue. Patients with gene signatures of such clonotypic expansion respond best to anti-PDL1 therapy. Notably, expanded clonotypes found in the tumour and normal adjacent tissue can also typically be detected in peripheral blood, which suggests a convenient approach to patient identification. Analyses of our data together with several external datasets suggest that intratumoural T cells, especially in responsive patients, are replenished with fresh, non-exhausted replacement cells from sites outside the tumour, suggesting continued activity of the cancer immunity cycle in these patients, the acceleration of which may be associated with clinical response.

Integration of eQTL and a Single-Cell Atlas in the Human Eye Identifies Causal Genes for Age-Related Macular Degeneration.

Age-related macular degeneration (AMD) is a leading cause of vision loss. To better understand disease pathogenesis and identify causal genes in GWAS loci for AMD risk, we present a comprehensive database of human retina and retinal pigment epithelium (RPE). Our database comprises macular and non-macular RNA sequencing (RNA-seq) profiles from 129 donors, a genome-wide expression quantitative trait loci (eQTL) dataset that includes macula-specific retina and RPE/choroid, and single-nucleus RNA-seq (NucSeq) from human retina and RPE with subtype resolution from more than 100,000 cells. Using NucSeq, we find enriched expression of AMD candidate genes in RPE cells. We identify 15 putative causal genes for AMD on the basis of co-localization of genetic association signals for AMD risk and eye eQTL, including the genes TSPAN10 and TRPM1. These results demonstrate the value of our human eye database for elucidating genetic pathways and potential therapeutic targets for ocular diseases.

The Indian cobra reference genome and transcriptome enables comprehensive identification of venom toxins.

Snakebite envenoming is a serious and neglected tropical disease that kills ~100,000 people annually. High-quality, genome-enabled comprehensive characterization of toxin genes will facilitate development of effective humanized recombinant antivenom. We report a de novo near-chromosomal genome assembly of Naja naja, the Indian cobra, a highly venomous, medically important snake. Our assembly has a scaffold N50 of 223.35 Mb, with 19 scaffolds containing 95% of the genome. Of the 23,248 predicted protein-coding genes, 12,346 venom-gland-expressed genes constitute the 'venom-ome' and this included 139 genes from 33 toxin families. Among the 139 toxin genes were 19 'venom-ome-specific toxins' (VSTs) that showed venom-gland-specific expression, and these probably encode the minimal core venom effector proteins. Synthetic venom reconstituted through recombinant VST expression will aid in the rapid development of safe and effective synthetic antivenom. Additionally, our genome could serve as a reference for snake genomes, support evolutionary studies and enable venom-driven drug discovery.

S100a4 upregulation in Pik3caH1047R;Trp53R270H;MMTV-Cre-driven mammary tumors promotes metastasis.

BACKGROUND: PIK3CA mutations are frequent in human breast cancer. Pik3caH1047R mutant expression in mouse mammary gland promotes tumorigenesis. TP53 mutations co-occur with PIK3CA mutations in human breast cancers. We previously generated a conditionally activatable Pik3caH1047R;MMTV-Cre mouse model and found a few malignant sarcomatoid (spindle cell) carcinomas that had acquired spontaneous dominant-negative Trp53 mutations. METHODS: A Pik3caH1047R;Trp53R270H;MMTV-Cre double mutant mouse breast cancer model was generated. Tumors were characterized by histology, marker analysis, transcriptional profiling, single-cell RNA-seq, and bioinformatics. Cell lines were developed from mutant tumors and used to identify and confirm genes involved in metastasis. RESULTS: We found Pik3caH1047R and Trp53R270H cooperate in driving oncogenesis in mammary glands leading to a shorter latency than either alone. Double mutant mice develop multiple histologically distinct mammary tumors, including adenocarcinoma and sarcomatoid (spindle cell) carcinoma. We found some tumors to be invasive and a few metastasized to the lung and/or the lymph node. Single-cell RNA-seq analysis of the tumors identified epithelial, stromal, myeloid, and T cell groups. Expression analysis of the metastatic tumors identified S100a4 as a top candidate gene associated with metastasis. Metastatic tumors contained a much higher percentage of epithelial-mesenchymal transition (EMT)-signature positive and S100a4-expressing cells. CRISPR/CAS9-mediated knockout of S100a4 in a metastatic tumor-derived cell line disrupted its metastatic potential indicating a role for S100a4 in metastasis. CONCLUSIONS: Pik3caH1047R;Trp53R270H;MMTV-Cre mouse provides a preclinical model to mimic a subtype of human breast cancers that carry both PIK3CA and TP53 mutations. It also allows for understanding the cooperation between the two mutant genes in tumorigenesis. Our model also provides a system to study metastasis and develop therapeutic strategies for PIK3CA/TP53 double-positive cancers. S100a4 found involved in metastasis in this model can be a potential diagnostic and therapeutic target.

Massively parallel single-cell B-cell receptor sequencing enables rapid discovery of diverse antigen-reactive antibodies.

Obtaining full-length antibody heavy- and light-chain variable regions from individual B cells at scale remains a challenging problem. Here we use high-throughput single-cell B-cell receptor sequencing (scBCR-seq) to obtain accurately paired full-length variable regions in a massively parallel fashion. We sequenced more than 250,000 B cells from rat, mouse and human repertoires to characterize their lineages and expansion. In addition, we immunized rats with chicken ovalbumin and profiled antigen-reactive B cells from lymph nodes of immunized animals. The scBCR-seq data recovered 81% (n = 56/69) of B-cell lineages identified from hybridomas generated from the same set of B cells subjected to scBCR-seq. Importantly, scBCR-seq identified an additional 710 candidate lineages not recovered as hybridomas. We synthesized, expressed and tested 93 clones from the identified lineages and found that 99% (n = 92/93) of the clones were antigen-reactive. Our results establish scBCR-seq as a powerful tool for antibody discovery.

An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors.

By leveraging tumorgraft (patient-derived xenograft) RNA-sequencing data, we developed an empirical approach, DisHet, to dissect the tumor microenvironment (eTME). We found that 65% of previously defined immune signature genes are not abundantly expressed in renal cell carcinoma (RCC) and identified 610 novel immune/stromal transcripts. Using eTME, genomics, pathology, and medical record data involving >1,000 patients, we established an inflamed pan-RCC subtype (IS) enriched for regulatory T cells, natural killer cells, TH1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. IS is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. The IS subtype correlated with systemic manifestations of inflammation such as thrombocytosis and anemia, which are enigmatic predictors of poor prognosis. Furthermore, IS was a strong predictor of poor survival. Our analyses suggest that tumor cells drive the stromal immune response. These data provide a missing link between tumor cells, the TME, and systemic factors.Significance: We undertook a novel empirical approach to dissect the renal cell carcinoma TME by leveraging tumorgrafts. The dissection and downstream analyses uncovered missing links between tumor cells, the TME, systemic manifestations of inflammation, and poor prognosis. Cancer Discov; 8(9); 1142-55. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

Comprehensive genomic analysis identifies pathogenic variants in maturity-onset diabetes of the young (MODY) patients in South India.

BACKGROUND: Maturity-onset diabetes of the young (MODY) is an early-onset, autosomal dominant form of non-insulin dependent diabetes. Genetic diagnosis of MODY can transform patient management. Earlier data on the genetic predisposition to MODY have come primarily from familial studies in populations of European origin. METHODS: In this study, we carried out a comprehensive genomic analysis of 289 individuals from India that included 152 clinically diagnosed MODY cases to identify variants in known MODY genes. Further, we have analyzed exome data to identify putative MODY relevant variants in genes previously not implicated in MODY. Functional validation of MODY relevant variants was also performed. RESULTS: We found MODY 3 (HNF1A; 7.2%) to be most frequently mutated followed by MODY 12 (ABCC8; 3.3%). They together account for ~ 11% of the cases. In addition to known MODY genes, we report the identification of variants in RFX6, WFS1, AKT2, NKX6-1 that may contribute to development of MODY. Functional assessment of the NKX6-1 variants showed that they are functionally impaired. CONCLUSIONS: Our findings showed HNF1A and ABCC8 to be the most frequently mutated MODY genes in south India. Further we provide evidence for additional MODY relevant genes, such as NKX6-1, and these require further validation.

Massively parallel nanowell-based single-cell gene expression profiling.

BACKGROUND: Technological advances have enabled transcriptome characterization of cell types at the single-cell level providing new biological insights. New methods that enable simple yet high-throughput single-cell expression profiling are highly desirable. RESULTS: Here we report a novel nanowell-based single-cell RNA sequencing system, ICELL8, which enables processing of thousands of cells per sample. The system employs a 5,184-nanowell-containing microchip to capture ~1,300 single cells and process them. Each nanowell contains preprinted oligonucleotides encoding poly-d(T), a unique well barcode, and a unique molecular identifier. The ICELL8 system uses imaging software to identify nanowells containing viable single cells and only wells with single cells are processed into sequencing libraries. Here, we report the performance and utility of ICELL8 using samples of increasing complexity from cultured cells to mouse solid tissue samples. Our assessment of the system to discriminate between mixed human and mouse cells showed that ICELL8 has a low cell multiplet rate (< 3%) and low cross-cell contamination. We characterized single-cell transcriptomes of more than a thousand cultured human and mouse cells as well as 468 mouse pancreatic islets cells. We were able to identify distinct cell types in pancreatic islets, including alpha, beta, delta and gamma cells. CONCLUSIONS: Overall, ICELL8 provides efficient and cost-effective single-cell expression profiling of thousands of cells, allowing researchers to decipher single-cell transcriptomes within complex biological samples.

Recurrent Loss of NFE2L2 Exon 2 Is a Mechanism for Nrf2 Pathway Activation in Human Cancers.

The Nrf2 pathway is frequently activated in human cancers through mutations in Nrf2 or its negative regulator KEAP1. Using a cell-line-derived gene signature for Nrf2 pathway activation, we found that some tumors show high Nrf2 activity in the absence of known mutations in the pathway. An analysis of splice variants in oncogenes revealed that such tumors express abnormal transcript variants from the NFE2L2 gene (encoding Nrf2) that lack exon 2, or exons 2 and 3, and encode Nrf2 protein isoforms missing the KEAP1 interaction domain. The Nrf2 alterations result in the loss of interaction with KEAP1, Nrf2 stabilization, induction of a Nrf2 transcriptional response, and Nrf2 pathway dependence. In all analyzed cases, transcript variants were the result of heterozygous genomic microdeletions. Thus, we identify an alternative mechanism for Nrf2 pathway activation in human tumors and elucidate its functional consequences.

Prediction and Quantification of Splice Events from RNA-Seq Data.

Analysis of splice variants from short read RNA-seq data remains a challenging problem. Here we present a novel method for the genome-guided prediction and quantification of splice events from RNA-seq data, which enables the analysis of unannotated and complex splice events. Splice junctions and exons are predicted from reads mapped to a reference genome and are assembled into a genome-wide splice graph. Splice events are identified recursively from the graph and are quantified locally based on reads extending across the start or end of each splice variant. We assess prediction accuracy based on simulated and real RNA-seq data, and illustrate how different read aligners (GSNAP, HISAT2, STAR, TopHat2) affect prediction results. We validate our approach for quantification based on simulated data, and compare local estimates of relative splice variant usage with those from other methods (MISO, Cufflinks) based on simulated and real RNA-seq data. In a proof-of-concept study of splice variants in 16 normal human tissues (Illumina Body Map 2.0) we identify 249 internal exons that belong to known genes but are not related to annotated exons. Using independent RNA samples from 14 matched normal human tissues, we validate 9/9 of these exons by RT-PCR and 216/249 by paired-end RNA-seq (2 x 250 bp). These results indicate that de novo prediction of splice variants remains beneficial even in well-studied systems. An implementation of our method is freely available as an R/Bioconductor package [Formula: see text].