Genome sequencing reveals loci under artificial selection that underlie disease phenotypes in the laboratory rat.

Large numbers of inbred laboratory rat strains have been developed for a range of complex disease phenotypes. To gain insights into the evolutionary pressures underlying selection for these phenotypes, we sequenced the genomes of 27 rat strains, including 11 models of hypertension, diabetes, and insulin resistance, along with their respective control strains. Altogether, we identified more than 13 million single-nucleotide variants, indels, and structural variants across these rat strains. Analysis of strain-specific selective sweeps and gene clusters implicated genes and pathways involved in cation transport, angiotensin production, and regulators of oxidative stress in the development of cardiovascular disease phenotypes in rats. Many of the rat loci that we identified overlap with previously mapped loci for related traits in humans, indicating the presence of shared pathways underlying these phenotypes in rats and humans. These data represent a step change in resources available for evolutionary analysis of complex traits in disease models.

Detecting Förster resonance energy transfer in living cells by conventional and spectral flow cytometry.

Assays based on Förster resonance energy transfer (FRET) can be used to study many processes in cell biology. Although this is most often done with microscopy for fluorescence detection, we report two ways to measure FRET in living cells by flow cytometry. Using a conventional flow cytometer and the "3-cube method" for intensity-based calculation of FRET efficiency, we measured the enzymatic activity of specific kinases in cells expressing a genetically-encoded reporter. For both AKT and protein kinase A, the method measured kinase activity in time-course, dose-response, and kinetic assays. Using the Cytek Aurora spectral flow cytometer, which applies linear unmixing to emission measured in multiple wavelength ranges, FRET from the same reporters was measured with greater single-cell precision, in real time and in the presence of other fluorophores. Results from gene-knockout studies suggested that spectral flow cytometry might enable the sorting of cells on the basis of FRET. The methods we present provide convenient and flexible options for using FRET with flow cytometry in studies of cell biology.

Subtype-specific and co-occurring genetic alterations in B-cell non-Hodgkin lymphoma.

B-cell non-Hodgkin lymphoma (B-NHL) encompasses multiple clinically and phenotypically distinct subtypes of malignancy with unique molecular etiologies. Common subtypes of B-NHL, such as diffuse large B-cell lymphoma, have been comprehensively interrogated at the genomic level, but rarer subtypes, such as mantle cell lymphoma, remain less extensively characterized. Furthermore, multiple B-NHL subtypes have thus far not been comprehensively compared using the same methodology to identify conserved or subtype-specific patterns of genomic alterations. Here, we employed a large targeted hybrid-capture sequencing approach encompassing 380 genes to interrogate the genomic landscapes of 685 B-NHL tumors at high depth, including diffuse large B-cell lymphoma, mantle cell lymphoma, follicular lymphoma, and Burkitt lymphoma. We identified conserved hallmarks of B-NHL that were deregulated in the majority of tumors from each subtype, including frequent genetic deregulation of the ubiquitin proteasome system. In addition, we identified subtype-specific patterns of genetic alterations, including clusters of co-occurring mutations and DNA copy number alterations. The cumulative burden of mutations within a single cluster were more discriminatory of B-NHL subtypes than individual mutations, implicating likely patterns of genetic cooperation that contribute to disease etiology. We therefore provide the first cross-sectional analysis of mutations and DNA copy number alterations across major B-NHL subtypes and a framework of co-occurring genetic alterations that deregulate genetic hallmarks and likely cooperate in lymphomagenesis.

Genome variation and conserved regulation identify genomic regions responsible for strain specific phenotypes in rat.

BACKGROUND: The genomes of laboratory rat strains are characterised by a mosaic haplotype structure caused by their unique breeding history. These mosaic haplotypes have been recently mapped by extensive sequencing of key strains. Comparison of genomic variation between two closely related rat strains with different phenotypes has been proposed as an effective strategy for the discovery of candidate strain-specific regions involved in phenotypic differences. We developed a method to prioritise strain-specific haplotypes by integrating genomic variation and genomic regulatory data predicted to be involved in specific phenotypes. Specifically, we aimed to identify genomic regions associated with Metabolic Syndrome (MetS), a disorder of energy utilization and storage affecting several organ systems. RESULTS: We compared two Lyon rat strains, Lyon Hypertensive (LH) which is susceptible to MetS, and Lyon Low pressure (LL), which is susceptible to obesity as an intermediate MetS phenotype, with a third strain (Lyon Normotensive, LN) that is resistant to both MetS and obesity. Applying a novel metric, we ranked the identified strain-specific haplotypes using evolutionary conservation of the occupancy three liver-specific transcription factors (HNF4A, CEBPA, and FOXA1) in five rodents including rat. Consideration of regulatory information effectively identified regions with liver-associated genes and rat orthologues of human GWAS variants related to obesity and metabolic traits. We attempted to find possible causative variants and compared them with the candidate genes proposed by previous studies. In strain-specific regions with conserved regulation, we found a significant enrichment for published evidence to obesity-one of the metabolic symptoms shown by the Lyon strains-amongst the genes assigned to promoters with strain-specific variation. CONCLUSIONS: Our results show that the use of functional regulatory conservation is a potentially effective approach to select strain-specific genomic regions associated with phenotypic differences among Lyon rats and could be extended to other systems.

Genomic structure of nucleotide diversity among Lyon rat models of metabolic syndrome.

BACKGROUND: The metabolic syndrome (MetS), a complex disorder involving hypertension, obesity, dyslipidemia and insulin resistance, is a major risk factor for heart disease, stroke, and diabetes. The Lyon Hypertensive (LH), Lyon Normotensive (LN) and Lyon Low-pressure (LL) rats are inbred strains simultaneously derived from a common outbred Sprague Dawley colony by selection for high, normal, and low blood pressure, respectively. Further studies found that LH is a MetS susceptible strain, while LN is resistant and LL has an intermediate phenotype. Whole genome sequencing determined that, while the strains are phenotypically divergent, they are nearly 98% similar at the nucleotide level. Using the sequence of the three strains, we applied an approach that harnesses the distribution of Observed Strain Differences (OSD), or nucleotide diversity, to distinguish genomic regions of identity-by-descent (IBD) from those with divergent ancestry between the three strains. This information was then used to fine-map QTL identified in a cross between LH and LN rats in order to identify candidate genes causing the phenotypes. RESULTS: We identified haplotypes that, in total, contain at least 95% of the identifiable polymorphisms between the Lyon strains that are likely of differing ancestral origin. By intersecting the identified haplotype blocks with Quantitative Trait Loci (QTL) previously identified in a cross between LH and LN strains, the candidate QTL regions have been narrowed by 78%. Because the genome sequence has been determined, we were further able to identify putative functional variants in genes that are candidates for causing the QTL. CONCLUSIONS: Whole genome sequence analysis between the LH, LN, and LL strains identified the haplotype structure of these three strains and identified candidate genes with sequence variants predicted to affect gene function. This approach, merged with additional integrative genetics approaches, will likely lead to novel mechanisms underlying complex disease and provide new drug targets and therapies.

Smart Start: Rituximab, Lenalidomide, and Ibrutinib in Patients With Newly Diagnosed Large B-Cell Lymphoma.

PURPOSE: Chemoimmunotherapy for patients with newly diagnosed diffuse large B-cell lymphoma (DLBCL) is largely unchanged for decades. Both preclinical models and clinical data suggest the combination of lenalidomide and ibrutinib may have synergy in DLBCL, particularly in the non-germinal center B-cell-like subset. METHODS: We enrolled 60 patients with newly diagnosed non-germinal center B-cell-like DLBCL in this investigator-initiated, single-arm phase II trial of rituximab, lenalidomide, and ibrutinib (RLI) with the sequential addition of chemotherapy (ClinicalTrials.gov identifier: NCT02636322). Patients were treated with rituximab 375 mg/m2 intravenous once on day 1, lenalidomide 25 mg once per day on days 1-10, and ibrutinib 560 mg once daily continuously of each 21-day cycle (RLI). After two cycles, standard chemotherapy was added to RLI for six additional cycles. The primary end points were overall response rate (ORR) after two cycles of RLI alone and complete response rate after completion of RLI with chemotherapy. In evaluable samples, circulating tumor DNA and DLBCL90 assays were performed. RESULTS: The median age was 63.5 years (range, 29-83 years) with 28% age 70 years or older. The revised international prognostic index identified 42% as high risk, and 62% were double expressor of MYC and BCL2 protein. The ORR after two cycles of RLI was 86.2%, and the complete response rate at the end of RLI-chemotherapy was 94.5%. With a median follow-up of 31 months, the progression-free survival and overall survival were at 91.3% and 96.6% at 2 years, respectively. CONCLUSION: Smart Start is the first study, to our knowledge, to treat newly diagnosed DLBCL with a targeted therapy combination before chemotherapy. RLI produced a high ORR, and RLI with chemotherapy resulted in durable responses. This establishes the potential for developing biologically driven and noncytotoxic first-line therapies for DLBCL.

Safety and activity of pembrolizumab in combination with rituximab in relapsed or refractory follicular lymphoma.

PD-1 blockade enhances the function of antitumor T cells and antibody-dependent, cell-mediated cytotoxicity (ADCC) of NK cells. In a single-center, open-label, phase 2 trial, we tested the combination of pembrolizumab, an anti-PD-1 monoclonal antibody, and rituximab, an anti-CD20 monoclonal antibody that induces ADCC, in 30 patients with follicular lymphoma (FL) with rituximab-sensitive disease who had relapsed after ≥1 prior therapy. Pembrolizumab was administered at 200 mg IV every 3 weeks for up to 16 cycles, and rituximab was given at 375 mg/m2 IV weekly for 4 weeks in cycle 1 only. The most common grade 3/4 adverse events (AEs) were liver enzyme abnormalities (3%), diarrhea (3%), nausea (3%), aseptic meningitis (3%), and pancreatitis (3%). Low-grade immune-related AEs were reported in 80% of patients, including diarrhea (43%), liver enzyme abnormalities (33%), thyroid dysfunction (27%), and rash (23%). Grade 3 or 4 immune-related AEs occurred in 13% of the patients. Treatment-related AEs led to discontinuation in 6 (20%) patients. The overall response rate (primary end point) was 67%, and the complete response (CR) rate was 50%. Median progression-free survival (PFS) was 12.6 months (95% confidence interval, 8.2-27.6), the 3-year overall survival rate was 97%, and 23% of patients were in remission at a median follow-up of 35 months. The presence of a high CD8+ T-effector score at baseline in the tumor was associated with induction of a CR and improved PFS. In this single-arm, phase 2 study, the combination of pembrolizumab and rituximab demonstrates favorable efficacy and safety profile in relapsed FL. This trial is registered at www.clinicaltrials.gov as #NCT02446457.

Follicular Lymphoma Microenvironment Characteristics Associated with Tumor Cell Mutations and MHC Class II Expression.

Follicular lymphoma (FL) is a B-cell malignancy with a complex tumor microenvironment that is rich in nonmalignant immune cells. We applied single-cell RNA sequencing to characterize the diverse tumor and immune cell populations of FL and identified major phenotypic subsets of FL T cells, including a cytotoxic CD4 T-cell population. We characterized four major FL subtypes with differential representation or relative depletion of distinct T-cell subsets. By integrating exome sequencing, we observed that somatic mutations are associated with, but not definitive for, reduced MHC expression on FL cells. In turn, expression of MHCII genes by FL cells was associated with significant differences in the proportions and targetable immunophenotypic characteristics of T cells. This provides a classification framework of the FL microenvironment in association with FL genotypes and MHC expression, and informs different potential immunotherapeutic strategies based upon tumor cell MHCII expression. SIGNIFICANCE: We have characterized the FL-infiltrating T cells, identified cytotoxic CD4 T cells as an important component that is associated with tumor cell-intrinsic characteristics, and identified sets of targetable immune checkpoints on T cells that differed from FLs with normal versus low MHC expression. See related commentary by Melnick, p. 374. This article is highlighted in the In This Issue feature, p. 369.

Inhibition of translation initiation factor eIF4a inactivates heat shock factor 1 (HSF1) and exerts anti-leukemia activity in AML.

Eukaryotic initiation factor 4A (eIF4A), the enzymatic core of the eIF4F complex essential for translation initiation, plays a key role in the oncogenic reprogramming of protein synthesis, and thus is a putative therapeutic target in cancer. As important component of its anticancer activity, inhibition of translation initiation can alleviate oncogenic activation of HSF1, a stress-inducible transcription factor that enables cancer cell growth and survival. Here, we show that primary acute myeloid leukemia (AML) cells exhibit the highest transcript levels of eIF4A1 compared to other cancer types. eIF4A inhibition by the potent and specific compound rohinitib (RHT) inactivated HSF1 in these cells, and exerted pronounced in vitro and in vivo anti-leukemia effects against progenitor and leukemia-initiating cells, especially those with FLT3-internal tandem duplication (ITD). In addition to its own anti-leukemic activity, genetic knockdown of HSF1 also sensitized FLT3-mutant AML cells to clinical FLT3 inhibitors, and this synergy was conserved in FLT3 double-mutant cells carrying both ITD and tyrosine kinase domain mutations. Consistently, the combination of RHT and FLT3 inhibitors was highly synergistic in primary FLT3-mutated AML cells. Our results provide a novel therapeutic rationale for co-targeting eIF4A and FLT3 to address the clinical challenge of treating FLT3-mutant AML.

Characteristics of anti-CD19 CAR T cell infusion products associated with efficacy and toxicity in patients with large B cell lymphomas.

Autologous chimeric antigen receptor (CAR) T cell therapies targeting CD19 have high efficacy in large B cell lymphomas (LBCLs), but long-term remissions are observed in less than half of patients, and treatment-associated adverse events, such as immune effector cell-associated neurotoxicity syndrome (ICANS), are a clinical challenge. We performed single-cell RNA sequencing with capture-based cell identification on autologous axicabtagene ciloleucel (axi-cel) anti-CD19 CAR T cell infusion products to identify transcriptomic features associated with efficacy and toxicity in 24 patients with LBCL. Patients who achieved a complete response by positron emission tomography/computed tomography at their 3-month follow-up had three-fold higher frequencies of CD8 T cells expressing memory signatures than patients with partial response or progressive disease. Molecular response measured by cell-free DNA sequencing at day 7 after infusion was significantly associated with clinical response (P = 0.008), and a signature of CD8 T cell exhaustion was associated (q = 2.8 × 10-149) with a poor molecular response. Furthermore, a rare cell population with monocyte-like transcriptional features was associated (P = 0.0002) with high-grade ICANS. Our results suggest that heterogeneity in the cellular and molecular features of CAR T cell infusion products contributes to variation in efficacy and toxicity after axi-cel therapy in LBCL, and that day 7 molecular response might serve as an early predictor of CAR T cell efficacy.

Selective Inhibition of HDAC3 Targets Synthetic Vulnerabilities and Activates Immune Surveillance in Lymphoma.

CREBBP mutations are highly recurrent in B-cell lymphomas and either inactivate its histone acetyltransferase (HAT) domain or truncate the protein. Herein, we show that these two classes of mutations yield different degrees of disruption of the epigenome, with HAT mutations being more severe and associated with inferior clinical outcome. Genes perturbed by CREBBP mutation are direct targets of the BCL6-HDAC3 onco-repressor complex. Accordingly, we show that HDAC3-selective inhibitors reverse CREBBP-mutant aberrant epigenetic programming, resulting in: (i) growth inhibition of lymphoma cells through induction of BCL6 target genes such as CDKN1A and (ii) restoration of immune surveillance due to induction of BCL6-repressed IFN pathway and antigen-presenting genes. By reactivating these genes, exposure to HDAC3 inhibitors restored the ability of tumor-infiltrating lymphocytes to kill DLBCL cells in an MHC class I and II-dependent manner, and synergized with PD-L1 blockade in a syngeneic model in vivo. Hence, HDAC3 inhibition represents a novel mechanism-based immune epigenetic therapy for CREBBP-mutant lymphomas. SIGNIFICANCE: We have leveraged the molecular characterization of different types of CREBBP mutations to define a rational approach for targeting these mutations through selective inhibition of HDAC3. This represents an attractive therapeutic avenue for targeting synthetic vulnerabilities in CREBBP-mutant cells in tandem with promoting antitumor immunity.This article is highlighted in the In This Issue feature, p. 327.

Targetable genetic alterations of TCF4 (E2-2) drive immunoglobulin expression in diffuse large B cell lymphoma.

The activated B cell (ABC-like) subtype of diffuse large B cell lymphoma (DLBCL) is characterized by chronic activation of signaling initiated by immunoglobulin µ (IgM). By analyzing the DNA copy number profiles of 1000 DLBCL tumors, we identified gains of 18q21.2 as the most frequent genetic alteration in ABC-like DLBCL. Using integrative analysis of matched gene expression profiling data, we found that the TCF4 (E2-2) transcription factor gene was the target of these alterations. Overexpression of TCF4 in ABC-like DLBCL cell lines led to its occupancy on immunoglobulin (IGHM) and MYC gene enhancers and increased expression of these genes at the transcript and protein levels. Inhibition of TCF4 activity with dominant-negative constructs was synthetically lethal to ABC-like DLBCL cell lines harboring TCF4 DNA copy gains, highlighting these gains as an attractive potential therapeutic target. Furthermore, the TCF4 gene was one of the top BRD4-regulated genes in DLBCL cell lines. BET proteolysis-targeting chimera (PROTAC) ARV771 extinguished TCF4, MYC, and IgM expression and killed ABC-like DLBCL cells in vitro. In DLBCL xenograft models, ARV771 treatment reduced tumor growth and prolonged survival. This work highlights a genetic mechanism for promoting immunoglobulin signaling in ABC-like DLBCL and provides a functional rationale for the use of BET inhibitors in this disease.

Mitochondrial ClpP-Mediated Proteolysis Induces Selective Cancer Cell Lethality.

The mitochondrial caseinolytic protease P (ClpP) plays a central role in mitochondrial protein quality control by degrading misfolded proteins. Using genetic and chemical approaches, we showed that hyperactivation of the protease selectively kills cancer cells, independently of p53 status, by selective degradation of its respiratory chain protein substrates and disrupts mitochondrial structure and function, while it does not affect non-malignant cells. We identified imipridones as potent activators of ClpP. Through biochemical studies and crystallography, we show that imipridones bind ClpP non-covalently and induce proteolysis by diverse structural changes. Imipridones are presently in clinical trials. Our findings suggest a general concept of inducing cancer cell lethality through activation of mitochondrial proteolysis.

Predictive Gene Signatures Determine Tumor Sensitivity to MDM2 Inhibition.

Early clinical trials using murine double minute 2 (MDM2) inhibitors demonstrated proof-of-concept of p53-induced apoptosis by MDM2 inhibition in cancer cells; however, not all wild-type TP53 tumors are sensitive to MDM2 inhibition. Therefore, more potent inhibitors and biomarkers predictive of tumor sensitivity are needed. The novel MDM2 inhibitor DS-3032b is 10-fold more potent than the first-generation inhibitor nutlin-3a. TP53 mutations were predictive of resistance to DS-3032b, and allele frequencies of TP53 mutations were negatively correlated with sensitivity to DS-3032b. However, sensitivity to DS-3032b of TP53 wild-type tumors varied greatly. We thus used two methods to create predictive gene signatures. First, by comparing sensitivity to MDM2 inhibition with basal mRNA expression profiles in 240 cancer cell lines, a 175-gene signature was defined and validated in patient-derived tumor xenograft models and ex vivo human acute myeloid leukemia (AML) cells. Second, an AML-specific 1,532-gene signature was defined by performing random forest analysis with cross-validation using gene expression profiles of 41 primary AML samples. The combination of TP53 mutation status with the two gene signatures provided the best positive predictive values (81% and 82%, compared with 62% for TP53 mutation status alone). In addition, the top-ranked 50 genes selected from the AML-specific 1,532-gene signature conserved high predictive performance, suggesting that a more feasible size of gene signature can be generated through this method for clinical implementation. Our model is being tested in ongoing clinical trials of MDM2 inhibitors.Significance: This study demonstrates that gene expression profiling combined with TP53 mutational status predicts antitumor effects of MDM2 inhibitors in vitro and in vivoCancer Res; 78(10); 2721-31. ©2018 AACR.

Contribution of independent and pleiotropic genetic effects in the metabolic syndrome in a hypertensive rat.

Hypertension is a major risk factor for cardiovascular disease, Type 2 diabetes, and end organ failure, and is often found concomitant with disorders characteristic of the Metabolic Syndrome (MetS), including obesity, dyslipidemia, and insulin resistance. While the associated features often occur together, the pathway(s) or mechanism(s) linking hypertension in MetS are not well understood. Previous work determined that genetic variation on rat chromosome 17 (RNO17) contributes to several MetS-defining traits (including hypertension, obesity, and dyslipidemia) in the Lyon Hypertensive (LH) rat, a genetically determined MetS model. We hypothesized that at least some of the traits on RNO17 are controlled by a single gene with pleiotropic effects. To address this hypothesis, consomic and congenic strains were developed, whereby a defined fragment of RNO17 from the LH rat was substituted with the control Lyon Normotensive (LN) rat, and MetS phenotypes were measured in the resultant progeny. Compared to LH rats, LH-17LN consomic rats have significantly reduced body weight, blood pressure, and lipid profiles. A congenic strain (LH-17LNc), with a substituted fragment at the distal end of RNO17 (17q12.3; 74-97 Mb; rn4 assembly), showed differences from the LH rat in blood pressure and serum total cholesterol and triglycerides. Interestingly, there was no difference in body weight between the LH-17LNc and the parental LH rat. These data indicate that blood pressure and serum lipids are regulated by a gene(s) in the distal congenic interval, and could be due to pleiotropy. The data also indicate that body weight is not determined by the same gene(s) at this locus. Interestingly, only two small haplotypes spanning a total of approximately 0.5 Mb differ between the LH and LN genomes in the congenic interval. Genes in these haplotypes are strong candidate genes for causing dyslipidemia in the LH rat. Overall, MetS, even in a simplified genetic model such as the LH-17LN rat, is likely due to both independent and pleiotropic gene effects.

Systems biology with high-throughput sequencing reveals genetic mechanisms underlying the metabolic syndrome in the Lyon hypertensive rat.

BACKGROUND: The metabolic syndrome (MetS) is a collection of co-occurring complex disorders including obesity, hypertension, dyslipidemia, and insulin resistance. The Lyon hypertensive and Lyon normotensive rats are models of MetS sensitivity and resistance, respectively. To identify genetic determinants and mechanisms underlying MetS, an F2 intercross between Lyon hypertensive and Lyon normotensive was comprehensively studied. METHODS AND RESULTS: Multidimensional data were obtained including genotypes of 1536 single-nucleotide polymorphisms, 23 physiological traits, and >150 billion nucleotides of RNA-seq reads from the livers of F2 intercross offspring and parental rats. Phenotypic and expression quantitative trait loci (eQTL) were mapped. Application of systems biology methods identified 17 candidate MetS genes. Several putative causal cis-eQTL were identified corresponding with phenotypic QTL loci. We found an eQTL hotspot on rat chromosome 17 that is causally associated with multiple MetS-related traits and found RGD1562963, a gene regulated in cis by this eQTL hotspot, as the most likely eQTL driver gene directly affected by genetic variation between Lyon hypertensive and Lyon normotensive rats. CONCLUSIONS: Our study sheds light on the intricate pathogenesis of MetS and demonstrates that systems biology with high-throughput sequencing is a powerful method to study the pathogenesis of complex genetic diseases.