Search

Inhibition of B-cell receptor (BCR) signaling through the BTK inhibitor, ibrutinib, has generated a remarkable response in mantle cell lymphoma (MCL). However, approximately one third of patients do not respond well to the drug, and disease relapse on ibrutinib is nearly universal. Alternative therapeutic strategies aimed to prevent and overcome ibrutinib resistance are needed. We compared and contrasted the effects of selinexor, a selective inhibitor of nuclear export, with ibrutinib in six MCL cell lines that display differential intrinsic sensitivity to ibrutinib. We found that selinexor had a broader antitumor activity in MCL than ibrutinib. MCL cell lines resistant to ibrutinib remained sensitive to selinexor. We showed that selinexor induced apoptosis/cell-cycle arrest and XPO-1 knockdown also retarded cell growth. Furthermore, downregulation of the NFκB gene signature, as opposed to BCR signature, was a common feature that underlies the response of MCL to both selinexor and ibrutinib. Meanwhile, unaltered NFκB was associated with ibrutinib resistance. Mechnistically, selinexor induced nuclear retention of IκB that was accompanied by the reduction of DNA-binding activity of NFκB, suggesting that NFκB is trapped in an inhibitory complex. Coimmunoprecipitation confirmed that p65 of NFκB and IκB were physically associated. In primary MCL tumors, we further demonstrated that the number of cells with IκB nuclear retention was linearly correlated with the degree of apoptosis. Our data highlight the role of NFκB pathway in drug response to ibrutinib and selinexor and show the potential of using selinexor to prevent and overcome intrinsic ibrutinib resistance through NFκB inhibition.

Suppl. Fig.1. Selinexor changes cell cycle progression in both ibrutinib-sensitive and -resistant MCL cells. Suppl. Fig.2. Selinexor retains IκBα and NF-κB subunits P65/P50 in the nuclei of Granta-519 cells. Suppl. Fig. 3. Selinexor retains IκBα and NF-κB subunits P65 in the nuclei.

Supplemental Figures 1 and 2 (S1. Measurement of autophagic vesicles in PBMC and S2 Massively parallel sequencing and outcome.

A, Schematic describing the primary B-cell lymphoma PDX SCID/NSG-hu mouse model. Six-to 8-week-old male NSG or CB-17 SCID mice were implanted with human fetal bone. Approximately 4 to 6 weeks following bone implantation, one injection of 5 Ã- 106 patient tumor cells were administered directly into the human fetal bone chip as the first generation (G1). The levels of circulating human β2M in mouse serum were used to monitor tumor engraftment and growth in G1 mice. Once tumor burden was detected, the pieces of tumor with bone chip were subcutaneously transplanted into NSG mice as the second generation (G2). After G2, the consistent growths of tumor in NSG mice were acquired. Single-agent targeted therapy, drug-resistant treatment, or combination therapy, were administered in the indicated generations. B, Human β2M reflected the tumor burdens of 16 B-cell lymphomas in PDX-G1 mice. After 3-4 weeks of tumor inoculation, mouse sera were collected from the peripheral blood of tail vein, and the level of human β2M was determined by ELISA.

Purpose:Preclinical studies suggest SYK and JAK contribute to tumor-intrinsic and microenvironment-derived survival signals. The pharmacodynamics of cerdulatinib, a dual SYK/JAK inhibitor, and associations with tumor response were investigated.Patients and Methods:In a phase I dose-escalation study in adults with relapsed/refractory B-cell malignancies, cerdulatinib was administered orally to sequential dose-escalation cohorts using once-daily or twice-daily schedules. The study enrolled 8 patients with chronic lymphocytic leukemia (CLL)/small lymphocytic lymphoma (SLL), 13 with follicular lymphoma, 16 with diffuse large B-cell lymphoma (DLBCL), and 6 with mantle cell lymphoma. Correlation of tumor response with pharmacodynamic markers was determined in patients with meaningful clinical responses.Results:Following cerdulatinib administration, complete SYK and JAK pathway inhibition was achieved in whole blood of patients at tolerated exposures. Target inhibition correlated with serum cerdulatinib concentration, and IC50 values against B-cell antigen receptor (BCR), IL2, IL4, and IL6 signaling pathways were 0.27 to 1.11 μmol/L, depending on the phosphorylation event. Significant correlations were observed between SYK and JAK pathway inhibition and tumor response. Serum inflammation markers were reduced by cerdulatinib, and several significantly correlated with tumor response. Diminished expression of CD69 and CD86 (B-cell activation markers), CD5 (negative regulator of BCR signaling), and enhanced expression of CXCR4 were observed in 2 patients with CLL, consistent with BCR and IL4 suppression and loss of proliferative capacity.Conclusions:Cerdulatinib potently and selectively inhibited SYK/JAK signaling at tolerated exposures in patients with relapsed/refractory B-cell malignancies. The extent of target inhibition in whole-blood assays and suppression of inflammation correlated with tumor response. (ClinicalTrials.gov ID:NCT01994382).

Supplementary Table S3. OncoPlus Gene List and Identified Variants

PT8 was clinically primary ibrutinib resistant. PT8-MCL PDX was generated and passaged. G3 mice were treated with vehicle control, IBN, and CFZ (P < 0.01, CFZ versus vehicle control or IBN, n = 5).

Supplementary Table S2. The RPPA antibody list

Representative tumor mass and histology of tumor series passages of PT1 MCL and PT2 Burkitt's lymphoma (BL). B-C, PDX tumors invaded the mouse spleen and formed splenomegaly. D, Patient abdomen PET/CT images of PT1 and PT5 showed splenomegaly (left panel). Autopsy of PT1 G3 NSG mouse and PT5 G1 SCID-hu mouse also showed splenomegaly (middle panel), and H&E staining of mouse spleens in PDXs indicated extensive lymphoma involvement (right panel).

Purpose: Patients with B-cell lymphomas often relapse after frontline therapy, and novel therapies are urgently needed to provide long-term remission. We established B-cell lymphoma patient-derived xenograft (PDX) models to assess their ability to mimic tumor biology and to identify B-cell lymphoma patient treatment options.Experimental Design: We established the PDX models from 16 patients with diffuse large B-cell lymphoma, mantle cell lymphoma, follicular lymphoma, marginal zone lymphoma, or Burkitt lymphoma by inoculating the patient tumor cells into a human bone chip implanted into mice. We subjected the PDX models to histopathologic and phenotypical examination, sequencing, and drug efficacy analysis. Primary and acquired resistance to ibrutinib, an oral covalent inhibitor of Bruton tyrosine kinase, were investigated to elucidate the mechanisms underlying ibrutinib resistance and to identify drug treatments to overcome resistance.Results: The PDXs maintained the same biological, histopathologic, and immunophenotypical features, retained similar genetic mutations, and produced comparable drug responses with the original patient tumors. In the acquired ibrutinib-resistant PDXs, PLC-γ2, p65, and Src were downregulated; however, a PI3K signaling pathway member was upregulated. Inactivation of the PI3K pathway with the inhibitor idelalisib in combination with ibrutinib significantly inhibited the growth of the ibrutinib-resistant tumors. Furthermore, we used a PDX model derived from a clinically ibrutinib-relapsed patient to evaluate various therapeutic choices, ultimately eliminating the tumor cells in the patient's peripheral blood.Conclusions: Our results demonstrate that the B-cell lymphoma PDX model is an effective system to predict and personalize therapies and address therapeutic resistance in B-cell lymphoma patients. Clin Cancer Res; 23(15); 4212–23. ©2017 AACR.

Supplemental Appendix (Table S1, Figure S1, Figure S2, Figure S3)

Purpose:Everolimus inhibits the mTOR, activating cytoprotective autophagy. Hydroxychloroquine inhibits autophagy. On the basis of preclinical data demonstrating synergistic cytotoxicity when mTOR inhibitors are combined with an autophagy inhibitor, we launched a clinical trial of combined everolimus and hydroxychloroquine, to determine its safety and activity in patients with clear-cell renal cell carcinoma (ccRCC).Patients and Methods:Three centers conducted a phase I/II trial of everolimus 10 mg daily and hydroxychloroquine in patients with advanced ccRCC. The objectives were to determine the MTD of hydroxychloroquine with daily everolimus, and to estimate the rate of 6-month progression-free survival (PFS) in patients with ccRCC receiving everolimus/hydroxychloroquine after 1–3 prior treatment regimens. Correlative studies to identify patient subpopulations that achieved the most benefit included population pharmacokinetics, measurement of autophagosomes by electron microscopy, and next-generation tumor sequencing.Results:No dose-limiting toxicity was observed in the phase I trial. The recommended phase II dose of hydroxychloroquine 600 mg twice daily with everolimus was identified. Disease control [stable disease + partial response (PR)] occurred in 22 of 33 (67%) evaluable patients. PR was observed in 2 of 33 patients (6%). PFS ≥ 6 months was achieved in 15 of 33 (45%) of patients who achieved disease control.Conclusions:Combined hydroxychloroquine 600 mg twice daily with 10 mg daily everolimus was tolerable. The primary endpoint of >40% 6-month PFS rate was met. Hydroxychloroquine is a tolerable autophagy inhibitor in future RCC or other trials.

Freshly isolated lymphoma cells from PT2-BL-G6 and G7 PDXs. Cell viability was tested by CellTiter-Glo luminescent cell viability assay after 48-hour incubation with the indicated drug treatment. The in vitro drug and dose information for IBN and ABT-199 are listed in Supplementary Table S1. IBN, ibrutinib; PT, patient; BL, Burkitt's lymphoma.

Include supplementary Figure S1-S4 legends; Supplementary Table S1 and S4; and list Supplementary Table S2 and S3 titles. Supplementary Table S1. The panel of drugs used in vitro and in vivo with their doses and duration. Supplementary Table S4. The mean passage time per generation for 5 generations

Genomewide loss-of-function (LOF) screening using CRISPR (clustered regularly interspaced short palindromic repeats) has facilitated the discovery of novel gene functions across diverse physiological and pathophysiological systems. A challenge with conventional genomewide CRISPR/Cas9 libraries is the unwieldy size (60,000 to 120,000 constructs), which is resource intensive and prohibitive in some experimental contexts. One solution to streamlining CRISPR screening is by multiplexing two or more guides per gene on a single construct, which enables functional redundancy to compensate for suboptimal gene knockout by individual guides. In this regard, AsCas12a (Cpf1) and its derivatives, e.g., enhanced AsCas12a (enAsCas12a), have enabled multiplexed guide arrays to be specifically and efficiently processed for genome editing. Prior studies have established that multiplexed CRISPR/Cas12a libraries perform comparably to the larger equivalent CRISPR/Cas9 libraries, yet the most efficient CRISPR/Cas12a library design remains unresolved. Here, we demonstrate that CRISPR/Cas12a genomewide LOF screening performed optimally with three guides arrayed per gene construct and could be adapted to robotic cell culture without noticeable differences in screen performance. Thus, the conclusions from this study provide novel insight to streamlining genomewide LOF screening using CRISPR/Cas12a and robotic cell culture.

RNA as an effector of biological functions often adopts secondary and tertiary structural folds. Plasmodium falciparum is a deadly human pathogen responsible for the devastating disease called malaria. In this study, we measured the differential accumulation of RNA secondary structures in coding and noncoding transcripts from the asexual developmental cycle in P. falciparum in human red blood cells. Our comprehensive analysis, combining high-throughput nuclease mapping of RNA structures by duplex RNA-seq, immunoaffinity purification and RNA analysis, collectively measured differentially base-paired RNA regions during the parasite development. Our mapping data not only aligned to a diverse pool of RNAs with known structures but also enabled us to identify new structural RNA regions in the malaria genome. On average, ~71% of the genes with secondary structures are found to be protein coding mRNAs. Mapping pattern of these base-paired RNAs corresponded to all parts of protein-coding mRNAs, including 5’ UTR, CDS and 3’ UTR. In addition to histone family genes which are known to form secondary structures in their mRNAs, transcripts from genes which are important for transcriptional and post-transcriptional control, such as unique plant-like transcription factor family, ApiAP2, DNA/RNA binding protein family, Alba, ribosomal proteins and eukaryotic initiation factors involved in translational control and the ones important for RBC invasion and cytoadherence also show strong accumulation of duplex RNA reads in various asexual stages. Intriguingly, our study determined a positive relationship between mRNA structural contents and translation efficiency in P. falciparum asexual blood stages, suggesting an essential role of RNA structural changes in malaria gene expression programs.

SLAM-associated protein (SAP) is an adaptor molecule that facilitates critical effector functions in immune cells, and its deficiency causes X-linked lymphoproliferative disease type 1 in which effector responses directed against EBV are severely compromised. The primary objective of this study was to phenotypically and functionally characterize a rare, CD8 T cell–restricted bimodal SAP expression pattern observed in healthy, human donors with the widely used 1C9-SAP mAb clone. We initially observed this pattern during the clinical validation of our flow cytometry–based assay to diagnose X-linked lymphoproliferative disease type 1 in our laboratory. For this validation study, we used multiparameter flow cytometry to identify cytosolic SAP expression in lymphocyte subsets, and CD8 T cells from the donors displaying the rare SAP expression pattern mentioned above were separately further evaluated by intracellular cytokine and CD107a staining to examine polyfunctionality following PMA/ionomycin and HLA class I allele–restricted EBV peptide epitope–induced T cell activation. Our data revealed that SAP 1C9-hi CD8 T cells clearly displayed higher polyfunctional responses versus SAP 1C9-lo CD8 T cells following PMA/ionomycin stimulation. Furthermore, polyfunctional EBV-specific CD8 T cell responses segregated with the SAP 1C9-hi CD8 T cells and not the SAP 1C9-lo CD8 T cells. Additionally, and rather intriguingly, short- and long-term T cell stimulation selectively diminished the signal for the 1C9-hi subset. Overall, our data suggest that although rare, this unique SAP expression pattern merits further evaluation as it has the potential to provide some insight into fundamental processes as they might relate to host–pathogen dynamics.

In silico approaches for next-generation sequencing (NGS) data modeling have utility in the clinical laboratory as a tool for clinical assay validation. In silico NGS data can take a variety of forms, including pure simulated data or manipulated data files in which variants are inserted into existing data files. In silico data enable simulation of a range of variants that may be difficult to obtain from a single physical sample. Such data allow laboratories to more accurately test the performance of clinical bioinformatics pipelines without sequencing additional cases. For example, clinical laboratories may use in silico data to simulate low variant allele fraction variants to test the analytical sensitivity of variant calling software or simulate a range of insertion/deletion sizes to determine the performance of insertion/deletion calling software. In this article, the Working Group reviews the different types of in silico data with their strengths and limitations, methods to generate in silico data, and how data can be used in the clinical molecular diagnostic laboratory. Survey data indicate how in silico NGS data are currently being used. Finally, potential applications for which in silico data may become useful in the future are presented.

Importance Postmortem genetic testing of young individuals with sudden death has previously identified pathogenic gene variants. However, prior studies primarily considered highly penetrant monogenic variants, often without detailed decedent and family clinical information. Objective To assess genotype and phenotype risk in a diverse cohort of young decedents with sudden death and their families. Design, Setting, and Participants Pathological and whole-genome sequence analysis was conducted in a cohort referred from a national network of medical examiners. Cases were accrued prospectively from May 2015 to March 2019 across 24 US states. Analysis began September 2016 and ended November 2020. Exposures Evaluation of autopsy and clinical data integrated with whole-genome sequence data and family member evaluation. Results A total of 103 decedents (mean [SD] age at death, 23.7 [11.9] years; age range, 1-44 years), their surviving family members, and 140 sex- and genetic ancestry–matched controls were analyzed. Among 103 decedents, autopsy and clinical data review categorized 36 decedents with postmortem diagnoses, 23 decedents with findings of uncertain significance, and 44 with sudden unexplained death. Pathogenic/likely pathogenic (P/LP) genetic variants in arrhythmia or cardiomyopathy genes were identified in 13 decedents (12.6%). A multivariable analysis including decedent phenotype, ancestry, and sex demonstrated that younger decedents had a higher burden of P/LP variants and select variants of uncertain significance (effect size, −1.64;P = .001). These select, curated variants of uncertain significance in cardiac genes were more common in decedents than controls (83 of 103 decedents [86%] vs 100 of 140 controls [71%];P = .005), and decedents harbored more rare cardiac variants than controls (2.3 variants per individual vs 1.8 in controls;P = .006). Genetic testing of 31 parent-decedent trios and 14 parent-decedent dyads revealed 8 transmitted P/LP variants and 1 de novo P/LP variant. Incomplete penetrance was present in 6 of 8 parents who transmitted a P/LP variant. Conclusions and Relevance Whole-genome sequencing effectively identified P/LP variants in cases of sudden death in young individuals, implicating both arrhythmia and cardiomyopathy genes. Genomic analyses and familial phenotype association suggest potentially additive, oligogenic risk mechanisms for sudden death in this cohort.

Background: Black patients have been historically underrepresented in studies investigating molecular patterns in non-small cell lung cancer (NSCLC). We aimed to investigate differences in actionable mutations among patients at our urban, diverse medical center. Results: 146 patients were included (59 black, 76 white, 7 Asian, 3 Hispanic, 1 mixed). 35 patients had a targetable mutation. Seven black patients (11.8%) had a targetable mutation compared to 28 non-black patients (32.2%, p = 0.005). 15 black patients had PD-L1 expression ≥50% compared to 19 non-black (25.4% vs 21.8%, p = 0.69). Black patients had a higher TMB compared to non-black (15.3 mutations/Mb compared to 11.5 mutations/Mb, p = 0.001). In a multivariate analysis, TMB was driven by smoking (p < 0.01), without any additive interaction in black patients who smoke (p = 0.8). Conclusion: NSCLC tumors from black patients had a higher TMB and were less likely to carry a targetable mutation. The higher TMB seen was driven by a higher prevalence of smoking among black patients in our study, which may not reflect nationwide trends. Our results serve as a proof of concept that differences in molecular markers exist between black and non-black patients, and that these differences may impact the treatment options available to black patients. Methods: Retrospective chart review of patients with a diagnosis of NSCLC who underwent both PD-L1 testing and massively parallel sequencing (UCM-OncoPlus) was conducted. We examined whether high PD-L1 expression, tumor mutational burden (TMB), and presence of targetable mutations (EGFR, BRAF, ERBB2, RET or ALK translocations, ROS1 rearrangements) occur at different frequencies in tumors from black patients compared to non-black patients.

Renal cell carcinoma (RCC) with leiomyomatous stroma is a provisional category of RCC in the 2016 World Health Organization Classification of Tumors of the Urinary System. Microscopic examination of hematoxylin and eosin-stained sections reveals this entity to be well-circumscribed with tubulopapillary growth of cells with clear cytoplasm in a background of leiomyomatous stroma. Herein we describe the genetic features of 15 University of Chicago Medical Center archived cases with hematoxylin and eosin histology matching the provisional diagnosis. Immunohistochemical (IHC) stains revealed 1/15 of these tumors to be clear cell renal cell carcinoma (ccRCC) and 6/15 to be clear cell papillary renal cell carcinoma (ccpRCC), demonstrating the morphologic overlap with these discrete known entities. Interestingly 3/6 of the ccpRCCs had chromosome 18 gain suggesting there may be novel specific genetic changes in ccpRCC with leiomyomatous stroma. Of the remaining 8 tumors with IHC staining patterns that do not fit either ccRCC or ccpRCC only 3 of these had mutations in the recently described TCEB1 gene with concurrent monosomy of chromosome 8. These 3 cases had a somewhat unique IHC pattern that possibly could separate them from the 5 other non-ccRCC/non-ccpRCC cases. This descriptive study, although small, demonstrates the difficulty in applying the current World Health Organization provisional criteria at a single institution with suggestion of an immunohistochemcial panel that may assist in the diagnosis of TCEB1-mutated RCC with leiomyomatous stroma.

Purpose: Everolimus inhibits the mTOR, activating cytoprotective autophagy. Hydroxychloroquine inhibits autophagy. On the basis of preclinical data demonstrating synergistic cytotoxicity when mTOR inhibitors are combined with an autophagy inhibitor, we launched a clinical trial of combined everolimus and hydroxychloroquine, to determine its safety and activity in patients with clear-cell renal cell carcinoma (ccRCC). Patients and Methods: Three centers conducted a phase I/II trial of everolimus 10 mg daily and hydroxychloroquine in patients with advanced ccRCC. The objectives were to determine the MTD of hydroxychloroquine with daily everolimus, and to estimate the rate of 6-month progression-free survival (PFS) in patients with ccRCC receiving everolimus/hydroxychloroquine after 1–3 prior treatment regimens. Correlative studies to identify patient subpopulations that achieved the most benefit included population pharmacokinetics, measurement of autophagosomes by electron microscopy, and next-generation tumor sequencing. Results: No dose-limiting toxicity was observed in the phase I trial. The recommended phase II dose of hydroxychloroquine 600 mg twice daily with everolimus was identified. Disease control [stable disease + partial response (PR)] occurred in 22 of 33 (67%) evaluable patients. PR was observed in 2 of 33 patients (6%). PFS ≥ 6 months was achieved in 15 of 33 (45%) of patients who achieved disease control. Conclusions: Combined hydroxychloroquine 600 mg twice daily with 10 mg daily everolimus was tolerable. The primary endpoint of >40% 6-month PFS rate was met. Hydroxychloroquine is a tolerable autophagy inhibitor in future RCC or other trials.

The 2016 WHO classification of brain tumors represents a major step towards the integration of molecular data into pathologic diagnoses. Our institution has included massively parallel sequencing technology in the diagnostic work-up of all gliomas since January 2016. The utilized platform successfully identifies copy number variations, individual gene mutations, small insertions and deletions, and selected gene fusions. Herein, we retrospectively review the first 51 glial tumor samples run for clinical purposes using the UCM-OncoPlus platform, a 1213 gene targeted hybrid-capture next generation sequencing (NGS) panel. NGS paired with histomorphology and clinical data allowed for reliable, comprehensive, and cost-effective classification of all the analyzed gliomas (51/51) with minimal tissue required and without the need for additional testing. In addition to detecting all diagnostically relevant mutations according to the 2016 WHO system, our data suggest a large NGS-based platform may improve the accuracy of classifying gliomas beyond the 2016 WHO system, to provide truly personalized diagnostics. Furthermore, this methodology assists in classifying histologically challenging or clinically unusual cases. And, finally, the versatile nature of this testing methodology allows for near effortless expansion as new therapeutic targets and prognostic markers are discovered.

Lack of reliable reference samples containing different mutations of interest across large sets of disease-relevant loci limits the extensive validation clinical next-generation sequencing (NGS) assays and their associated bioinformatics pipelines. Herein, we have generated a publicly available, highly flexible tool, in silico Mutator (insiM), to introduce point mutations, insertions, deletions, and duplications of any size into real data sets of amplicon-based or hybrid-capture NGS assays. insiM accepts an alignment file along with target territory and produces paired-end FASTQ files containing specified mutations via modification of original sequencing reads. Mutant signal is, thus, generated within the context of existing real-world data to most closely mimic assay performance. Resulting files may then be passed through the assay's bioinformatics pipeline to assist with assay/bioinformatics validation and to identify performance gaps in detection. To establish the basic functionality of the software, a series of simulation experiments with varying mutation types, sizes, and allele frequencies were performed across the entire clinical territory of hybrid-capture and amplicon-based clinical assays developed at The University of Chicago. This work demonstrates the utility of insiM as a supplementary tool during the validation of an NGS assay's bioinformatics pipeline.

Classification and characterization of neuronal types are critical for understanding their function and dysfunction. Neuronal classification schemes typically rely on measurements of electrophysiological, morphological, and molecular features, but aligning such datasets has been challenging. Here, we present a unified classification of mouse retinal ganglion cells (RGCs), the sole retinal output neurons. We used visually-evoked responses to classify 1859 mouse RGCs into 42 types. We also obtained morphological or transcriptomic data from subsets and used these measurements to align the functional classification to publicly available morphological and transcriptomic data sets. We created an online database that allows users to browse or download the data and to classify RGCs from their light responses using a machine learning algorithm. This work provides a resource for studies of RGCs, their upstream circuits in the retina, and their projections in the brain, and establishes a framework for future efforts in neuronal classification and open data distribution.

Plasmodium falciparum is a deadly human pathogen responsible for the devastating disease called malaria. In this study, we measured the differential accumulation of RNA secondary structures in coding and non-coding transcripts from the asexual developmental cycle in P. falciparum in human red blood cells. Our comprehensive analysis that combined high-throughput nuclease mapping of RNA structures by duplex RNA-seq, SHAPE-directed RNA structure validation, immunoaffinity purification and characterization of antisense RNAs collectively measured differentially base-paired RNA regions throughout the parasite’s asexual RBC cycle. Our mapping data not only aligned to a diverse pool of RNAs with known structures but also enabled us to identify new structural RNA regions in the malaria genome. On average, approximately 71% of the genes with secondary structures are found to be protein coding mRNAs. The mapping pattern of these base-paired RNAs corresponded to all regions of mRNAs, including the 5ʹ UTR, CDS and 3ʹ UTR as well as the start and stop codons. Histone family genes which are known to form secondary structures in their mRNAs and transcripts from genes which are important for transcriptional and post-transcriptional control, such as the unique plant-like transcription factor family, ApiAP2, DNA-/RNA-binding protein, Alba3 and proteins important for RBC invasion and malaria cytoadherence also showed strong accumulation of duplex RNA reads in various asexual stages in P. falciparum. Intriguingly, our study determined stage-specific, dynamic relationships between mRNA structural contents and translation efficiency in P. falciparum asexual blood stages, suggesting an essential role of RNA structural changes in malaria gene expression programs. Abbreviations: CDS: Coding Sequence; DNA: Deoxyribonucleic Acid; dsRNA: double-stranded RNA; IDC: Intra-erythrocytic Developmental Cycle (IDC); m6A: N6-methyladenosine; mRNA: Messenger RNA; ncRNA: Non-coding RNA; RBC: Red Blood cells; RBP: RNA-Binding Protein; REC: Relative Expression Counts; RNA-seq: RNA-sequencing; RNA: Ribonucleic Acid; RNP: Ribonucleoprotein; RPKM: Reads Per Kilobase of transcript Per Million; rRNA: Ribosomal RNA 16. RUFs: RNAs of Unknown Function; SHAPE: Selective 2’-hydroxyl acylation analysed by primer extension; snoRNA: Small Nucleolar RNA; snRNA: Small Nuclear RNA; SRP-RNA: Signal Recognition Particle RNA; ssRNA: (Single-stranded RNA); TE: Translation Efficiency; tRNA: transfer RNA; UTR: Untranslated Region

Systematic implementation of bioinformatics resources for next generation sequencing (NGS)-based clinical testing is an arduous undertaking. One of the key challenges involves developing an ecosystem of information technology infrastructure for enabling scalable and reproducible bioinformatics services that is resilient and secure for handling genetic and protected health information, often embedded in an existing non-bioinformatics-oriented infrastructure. Container technology provides an ideal and infrastructure-agnostic solution for molecular laboratories developing and using bioinformatics pipelines, whether on-premise or using the cloud. A container is a technology that provides a consistent computational environment and enables reproducibility, scalability, and security when developing NGS bioinformatics analysis pipelines. Containers can increase the bioinformatics team's productivity by automating and simplifying the maintenance of complex bioinformatics resources, as well as facilitate validation, version control, and documentation necessary for clinical laboratory regulatory compliance. Although there is increasing popularity in adopting containers for developing NGS bioinformatics pipelines, there is wide variability and inconsistency in the usage of containers that may result in suboptimal performance and potentially compromise the security and privacy of protected health information. In this article, the authors highlight the current state and provide best or recommended practices for building, using containers in NGS bioinformatics solutions in a clinical setting with focus on scalability, optimization, maintainability, and data security.

Transcriptional analysis can be utilized to reconcile variants of uncertain significance, particularly those predicted to impact splicing. Laboratory analysis of the predicted mRNA transcript may allow inference of the in vivo impact of the variant and aid prediction of its clinical significance. We present a patient with classical features of primary ciliary dyskinesia (PCD) who was identified to have compound heterozygous variants in the DNAH11 gene (c.10691 + 2T > C, c.13523_13543dup21) via trio whole-exome sequencing in 2013. These variants were originally classified as Mutation and Likely Mutation. However, these variants were downgraded to variants of uncertain significance (VUSs) during reanalysis in 2016 because of uncertainty that they caused a loss of function of the gene. c.10691 + 2T > C is predicted to abrogate the canonical splice site and lead to the skipping of exon 65, but the adjoining of exon 64 and exon 66 in the DNAH11 transcript preserves the reading frame of the resultant protein. c.13523_13543dup21 is located in the last exon of the DNAH11 coding sequence, upstream of the canonical stop codon, which suggests a reduced likelihood to trigger nonsense-mediated decay (NMD). Transcriptional analysis was performed to characterize the impact of the variants, resulting in reclassification of c.10691 + 2T > C to Likely Pathogenic by providing evidence that it results in a deleterious effect and subsequent downstream reclassification of c.13523_13543dup21 to Likely Pathogenic as well. Our case illustrates the potential impact of transcriptional analysis on variant resolution, supporting its usage on variants that exert an unpredictable effect on splicing.

Introduction: Genetic testing after sudden death in the young can identify pathogenic cardiac gene variants. Hypothesis: Genomic methods, coupled with phenotype evaluation, reveal non-Mendalian risks. Methods: We conducted clinical analysis and whole genome sequencing on 103 decedents aged 1-40 (mean age at death 23.7 years) accrued prospectively from 2015 to 2019 across 22 states. Postmortem pathological findings were classified as: known cardiac disorders, findings of uncertain significance (FUS), or sudden unexplained death (SUD, indicating no postmortem pathological diagnosis). Parental DNA and clinical data were obtained where possible. Variants were classified by an independent clinical genetic laboratory. Results: Among the 103 decedents, 34 had a postmortem clinical diagnosis, 23 had FUS, and 46 were classified as SUD. Pathogenic/likely pathogenic (P/LP) variants in arrhythmia or cardiomyopathy genes were identified in 17 (16.5%) decedents. The distribution of P/LP variants was not associated with age at death (OR 1.01 [0.97, 1.05], p=0.54); however, a multivariable analysis including decedent phenotype, ancestry and sex demonstrated that younger decedents had a higher burden of curated P/LP/VUS variants (effect size -1.5, p=0.0019). DNA from 31 parent-decedent trios and 14 parent-decedent dyads revealed 9 transmitted P/LP variants and 1 de novo P/LP variant. More than half of parents transmitting a P/LP variant (5/9) did not have clinical findings associated with the genotype. Conclusions: Whole genome sequencing effectively revealed P/LP variants in cases of sudden death in the young, implicating both arrhythmia and cardiomyopathy genes. In addition, both genotype and phenotype analyses suggest additional non-Mendelian risk mechanisms.

In follicular thyroid neoplasms without invasion, a diagnosis of atypical adenoma (AA) (follicular tumor of uncertain malignant potential) may be rendered if atypical features (indefinite capsular/vascular invasion, necrosis, solid growth, increased mitoses) are present. This study compares clinical, histologic, and molecular features of patients with AAs (n=31), nonmetastatic follicular thyroid carcinoma (nmFTC) (n=18), and metastatic follicular thyroid carcinoma (mFTC) (n=38). Patients with mFTC were older. Mitotic activity in areas of solid growth was greatest in mFTC (P=0.05). Oncocytic tumors tended to show solid growth (P=0.04). The presence or frequency of capsular and/or vascular invasion was not different between nmFTC and mFTC. TERT promoter mutations were higher in patients with mFTC (50%) than nmFTC (25%) and AA (10%) (P=0.02). TERT promoter mutation was associated with necrosis (P=0.01) and solid growth plus increased mitoses (P=0.03). Necrosis and TERT promoter mutations were identified in all groups, most frequently in mFTC. The combination of solid growth with increased mitoses, necrosis, and TERT promoter mutation was only seen in follicular carcinomas. Poorly differentiated features, vascular invasion, and TERT promoter mutation correlated with metastasis in FTC. Given the low frequency of necrosis and TERT promoter mutation in AAs, close clinical follow-up is recommended in patients with these findings, especially if additional atypical features (such as solid growth plus mitoses) are present.

Inhibition of B-cell receptor (BCR) signaling through the BTK inhibitor, ibrutinib, has generated a remarkable response in mantle cell lymphoma (MCL). However, approximately one third of patients do not respond well to the drug, and disease relapse on ibrutinib is nearly universal. Alternative therapeutic strategies aimed to prevent and overcome ibrutinib resistance are needed. We compared and contrasted the effects of selinexor, a selective inhibitor of nuclear export, with ibrutinib in six MCL cell lines that display differential intrinsic sensitivity to ibrutinib. We found that selinexor had a broader antitumor activity in MCL than ibrutinib. MCL cell lines resistant to ibrutinib remained sensitive to selinexor. We showed that selinexor induced apoptosis/cell-cycle arrest and XPO-1 knockdown also retarded cell growth. Furthermore, downregulation of the NFκB gene signature, as opposed to BCR signature, was a common feature that underlies the response of MCL to both selinexor and ibrutinib. Meanwhile, unaltered NFκB was associated with ibrutinib resistance. Mechnistically, selinexor induced nuclear retention of IκB that was accompanied by the reduction of DNA-binding activity of NFκB, suggesting that NFκB is trapped in an inhibitory complex. Coimmunoprecipitation confirmed that p65 of NFκB and IκB were physically associated. In primary MCL tumors, we further demonstrated that the number of cells with IκB nuclear retention was linearly correlated with the degree of apoptosis. Our data highlight the role of NFκB pathway in drug response to ibrutinib and selinexor and show the potential of using selinexor to prevent and overcome intrinsic ibrutinib resistance through NFκB inhibition.

Cribriform-morular variant of papillary thyroid carcinoma (CMVPTC) is usually an inherited malignancy and may be a presenting indicator of familial adenomatous polyposis syndrome although it may occasionally be sporadic. Known CMVPTC mutations include adenomatous polyposis coli ( APC) and β-catenin ( CTNNB1) genes. Despite its malignant classification, CMVPTC is considered to be a well-differentiated thyroid tumor with a generally good behavior. In contrast, poorly differentiated thyroid carcinoma is an aggressive tumor. We report a case of CMVPTC with poorly differentiated features in a young female without phenotypic features of familial adenomatous polyposis but with known germline alterations of the APC gene. High throughput sequencing showed germline chromosome 5q deletion encompassing the APC gene in all components with additional unique genetic alterations in the somatic components. A single nucleotide substitution (c.1548+1G>A, NM_000038.5) located one base pair downstream of exon 12 of the APC gene was identified in the CMVPTC component, and a pathogenic frameshift deletion in exon 14 of APC (c.3642del, p.Ser1214Argfs*51, NM_000038.5) was identified in the poorly differentiated thyroid carcinoma component. No other cancer-associated genes were identified by our techniques. Our case represents a rare phenomenon of poorly differentiated features in association with CMVPTC. To our knowledge, ours is the only such report of poorly differentiated features arising in association with an inherited CMVPTC.

The oligometastasis hypothesis suggests a spectrum of metastatic virulence where some metastases are limited in extent and curable with focal therapies. A subset of patients with metastatic colorectal cancer achieves prolonged survival after resection of liver metastases consistent with oligometastasis. Here we define three robust subtypes of de novo colorectal liver metastasis through integrative molecular analysis. Patients with metastases exhibiting MSI-independent immune activation experience the most favorable survival. Subtypes with adverse outcomes demonstrate VEGFA amplification in concert with (i) stromal, mesenchymal, and angiogenic signatures, or (ii) exclusive NOTCH1 and PIK3C2B mutations with E2F/MYC activation. Molecular subtypes complement clinical risk stratification to distinguish low-risk, intermediate-risk, and high-risk patients with 10-year overall survivals of 94%, 45%, and 19%, respectively. Our findings provide a framework for integrated classification and treatment of metastasis and support the biological basis of curable oligometastatic colorectal cancer. These concepts may be applicable to many patients with metastatic cancer., The oligometastasis hypothesis suggests certain metastases are limited in extent and curable with focal therapies. Here they identify three integrated molecular subtypes of colorectal cancer liver metastasis, which complement clinical risk stratification to distinguish the subset of oligometastatic patients.

Context.—Proposed noninvasive follicular thyroid neoplasm with papillary-like nuclear features (NIFTPs), formerly noninvasive encapsulated papillary carcinoma, follicular variant (PTC-FV), is an indolent tumor with follicular growth and frequent RAS mutations.Objective.—To detect histologic and molecular differences separating NIFTP from follicular adenomas (FAs) and invasive carcinomas, particularly papillary carcinomas with extensive follicular growth (PTC-EFGs) and invasive encapsulated PTC-FV (IE-PTC-FV).Design.—Sixty-one tumors were reviewed histologically and reclassified into 32 NIFTPs (52%), 4 IE-PTC-FVs (7%), 14 PTC-EFGs (23%), and 11 FAs (18%). Next-generation sequencing for mutations in 50 genes was performed. Clinical outcomes were recorded.Results.—The NIFTPs and FAs were well circumscribed and unencapsulated. The FAs had bland nuclei, whereas the NIFTPs showed at least 2 of 3 (67%; sufficient) nuclear features (enlargement, irregular contours, chromatin clearing). The IE-PTC-FVs had follicular growth, sufficient nuclear features, and extensive capsular invasion. The PTC-EFGs had a median of 5% papillae with intrathyroidal invasion (broad-based, sclerotic, or small follicle growth patterns); intranuclear pseudoinclusions were present only in PTC-EFGs (9 of 14; 64%). Mutations included RAS in 20 of the 32 NIFTPs (62%), 4 of the 11 FAs (36%), and 3 of the 4 IE-PTC-FVs (75%); BRAF K601E in 1 NIFTP (3%); BRAF V600E in 5 PTC-EFGs (36%). No NIFTPs or FAs recurred or metastasized. All 4 IE-PTC-FVs (100%) had hematogenous metastasis. Two PTC-EFGs (14%) had lymphatic metastasis.Conclusions.—The morphologic similarity and RAS mutations in FAs, NIFTPs, and IE-PTC-FVs supports the genetic similarity of those follicular neoplasms in contrast to the unique presence of BRAF V600E mutations in PTC-EFGs. Using strict diagnostic criteria supported by molecular testing, tumors with extensive follicular growth can be classified into follicular type or RAS-like (FA, NIFTP, IE-PTC-FV) versus papillary type or BRAF V600E–like (PTC-EFG).

Next-generation sequencing (NGS)–based targeted gene capture panels are used to profile hematopoietic malignancies to guide prognostication and treatment decisions. Because these panels include genes associated with hereditary hematopoietic malignancies (HHMs), we hypothesized that these panels could identify pathogenic germ line variants in malignant cells, thereby identifying patients at risk for HHMs. In total, pathogenic or likely pathogenic variants in ANKRD26 , CEBPA , DDX41 , ETV6 , GATA2 , RUNX1 , or TP53 were identified in 74 (21%) of 360 patients. Germ line tissue was available for 24 patients with 25 pathogenic or likely pathogenic variants with variant allele frequencies >0.4. Six (24%) of these 25 variants were of germ line origin. Three DDX41 variants, 2 GATA2 variants, and a TP53 variant previously implicated in Li-Fraumeni syndrome were of germ line origin. No likely pathogenic/pathogenic germ line variants possessed variant allele frequencies DDX41 and GATA2 are especially likely to be of germ line origin. Thus, tumor-based panels may augment, but should not replace, comprehensive germ line–based testing and counseling.

Ibrutinib has generated remarkable responses in patients with chronic lymphocytic leukemia (CLL), including those with an unfavorable cytogenetic profile. However, patients develop resistance, with poor outcomes and no established treatment options. Mutations in BTK and PLCG2 have emerged as main mechanisms of drug resistance, but not all patients carry these mutations. Further understanding of mechanisms of resistance is urgently needed and will support rational development of new therapeutic strategies. To that end, we characterized the genomic profiles of serial samples from 9 patients with ibrutinib-relapsed disease, including 6 who had Richter transformation. Mutations, indels, copy-number aberrations, and loss of heterozygosity were assessed using next-generation sequencing and single-nucleotide polymorphism array. We found that 18p deletion (del(18p)), together with del(17p)/TP53 mutations, was present in 5 of 9 patients before ibrutinib therapy. In addition to BTKC481 , we identified BTKT316A , a structurally novel mutation located in the SH2 domain of BTK. Minor BTK clones with low allele frequencies were captured in addition to major BTK clones. Although TP53 loss predisposes patients for relapse, clone size of TP53 loss may diminish during disease progression while mutant BTK clone expands. In patients who had Richter transformation, we found that the transformed cells were clonal descendants of circulating leukemia cells but continued to undergo evolution and drifts. Surprisingly, transformed lymphoma cells in tissue may acquire a different BTK mutation from that in the CLL leukemia cells. Collectively, these results provide insights into clonal evolution underlying ibrutinib relapse and prompt further investigation on genomic abnormalities that have clinical application potential.