Longitudinal Evaluation of Circulating Tumor DNA Using Sensitive Amplicon-Based Next-Generation Sequencing to Identify Resistance Mechanisms to Immune Checkpoint Inhibitors for Advanced Urothelial Carcinoma.

Serial evaluation of circulating tumor DNA may allow noninvasive assessment of drivers of resistance to immune checkpoint inhibitors (ICIs) in advanced urothelial cancer (aUC). We used a novel, amplicon-based next-generation sequencing assay to identify genomic alterations (GAs) pre- and post-therapy in 39 patients with aUC receiving ICI and 6 receiving platinum-based chemotherapy (PBC). One or more GA was seen in 95% and 100% of pre- and post-ICI samples, respectively, commonly in TP53 (54% and 54%), TERT (49% and 59%), and BRCA1/BRCA2 (33% and 33%). Clearance of ≥1 GA was seen in 7 of 9 patients responding to ICI, commonly in TP53 (n = 4), PIK3CA (n = 2), and BRCA1/BRCA2 (n = 2). A new GA was seen in 17 of 20 patients progressing on ICI, frequently in BRCA1/BRCA2 (n = 6), PIK3CA (n = 3), and TP53 (n = 3), which seldom emerged in patients receiving PBC. These findings highlight the potential for longitudinal circulating tumor DNA evaluation in tracking response and resistance to therapy.

RNA-SeQC 2: efficient RNA-seq quality control and quantification for large cohorts.

SUMMARY: Post-sequencing quality control is a crucial component of RNA sequencing (RNA-seq) data generation and analysis, as sample quality can be affected by sample storage, extraction and sequencing protocols. RNA-seq is increasingly applied to cohorts ranging from hundreds to tens of thousands of samples in size, but existing tools do not readily scale to these sizes, and were not designed for a wide range of sample types and qualities. Here, we describe RNA-SeQC 2, an efficient reimplementation of RNA-SeQC (DeLuca et al., 2012) that adds multiple metrics designed to characterize sample quality across a wide range of RNA-seq protocols. AVAILABILITY AND IMPLEMENTATION: The command-line tool, documentation and C++ source code are available at the GitHub repository https://github.com/getzlab/rnaseqc. Code and data for reproducing the figures in this paper are available at https://github.com/getzlab/rnaseqc2-paper. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Let-7 represses Nr6a1 and a mid-gestation developmental program in adult fibroblasts.

MicroRNAs (miRNAs) are critical to proliferation, differentiation, and development. Here, we characterize gene expression in murine Dicer-null adult mesenchymal stem cell lines, a fibroblast cell type. Loss of Dicer leads to derepression of let-7 targets at levels that exceed 10-fold to 100-fold with increases in transcription. Direct and indirect targets of this miRNA belong to a mid-gestation embryonic program that encompasses known oncofetal genes as well as oncogenes not previously associated with an embryonic state. Surprisingly, this mid-gestation program represents a distinct period that occurs between the pluripotent state of the inner cell mass at embryonic day 3.5 (E3.5) and the induction of let-7 upon differentiation at E10.5. Within this mid-gestation program, we characterize the let-7 target Nr6a1, an embryonic transcriptional repressor that regulates gene expression in adult fibroblasts following miRNA loss. In total, let-7 is required for the continual suppression of embryonic gene expression in adult cells, a mechanism that may underlie its tumor-suppressive function.

Genome-wide impact of a recently expanded microRNA cluster in mouse.

Variations in microRNA (miRNA) gene and/or target repertoire are likely to be key drivers of phenotypic differences between species. To better understand these changes, we developed a computational method that identifies signatures of species-specific target site gain and loss associated with miRNA acquisition. Interestingly, several of the miRNAs implicated in mouse 3' UTR evolution derive from a single rapidly expanded rodent-specific miRNA cluster. Located in the intron of Sfmbt2, a maternally imprinted polycomb gene, these miRNAs (referred to as the Sfmbt2 cluster) are expressed in both embryonic stem cells and the placenta. One abundant miRNA from the cluster, miR-467a, functionally overlaps with the mir-290-295 cluster in promoting growth and survival of mouse embryonic stem cells. Predicted novel targets of the remaining cluster members are enriched in pathways regulating cell survival. Two relevant species-specific target candidates, Lats2 and Dedd2, were validated in cultured cells. We suggest that the rapid evolution of the Sfmbt2 cluster may be a result of intersex conflict for growth regulation in early mammalian development and could provide a general model for the genomic response to acquisition of miRNAs and similar regulatory factors.

A latent pro-survival function for the mir-290-295 cluster in mouse embryonic stem cells.

MicroRNAs (miRNAs) post-transcriptionally regulate the expression of thousands of distinct mRNAs. While some regulatory interactions help to maintain basal cellular functions, others are likely relevant in more specific settings, such as response to stress. Here we describe such a role for the mir-290-295 cluster, the dominant miRNA cluster in mouse embryonic stem cells (mESCs). Examination of a target list generated from bioinformatic prediction, as well as expression data following miRNA loss, revealed strong enrichment for apoptotic regulators, two of which we validated directly: Caspase 2, the most highly conserved mammalian caspase, and Ei24, a p53 transcriptional target. Consistent with these predictions, mESCs lacking miRNAs were more likely to initiate apoptosis following genotoxic exposure to gamma irradiation or doxorubicin. Knockdown of either candidate partially rescued this pro-apoptotic phenotype, as did transfection of members of the mir-290-295 cluster. These findings were recapitulated in a specific mir-290-295 deletion line, confirming that they reflect miRNA functions at physiological levels. In contrast to the basal regulatory roles previously identified, the pro-survival phenotype shown here may be most relevant to stressful gestations, where pro-oxidant metabolic states induce DNA damage. Similarly, this cluster may mediate chemotherapeutic resistance in a neoplastic context, making it a useful clinical target.

Comprehensive multiplexed autoantibody profiling of patients with advanced urothelial cancer.

BACKGROUND: Comprehensive profiling of autoantibodies (AAbs) in metastatic urothelial cancer (mUC) has not been performed to date. This may aid in diagnosis of UC, uncover novel therapeutic targets in this disease as well as identify associations between AAbs and response and toxicity to systemic therapies. METHODS: We used serum from patients with mUC collected prior to and after systemic therapy (immune checkpoint inhibitor (ICI) or platinum-based chemotherapy (PBC)) at Dana-Farber Cancer Institute. 38 age-matched and sex-matched healthy controls (HCs) from healthy blood donors were also evaluated. The SeroTag immuno-oncology discovery array (Oncimmune) was used, with quantification of the AAb reactivity toward 1132 antigens. Bound AAbs were detected using an anti-immunoglobulin G-specific detection antibody conjugated to the fluorescent reporter dye phycoerythrin. The AAb reactivity was reported as the median fluorescence intensity for each color and sample using a Luminex FlexMAP3D analyzer. Clinical outcomes of interest included radiographic response and development of immune-related adverse events (irAEs). Significance analysis of microarray was used to compare mUC versus HC and radiographic response. Associations with irAE were evaluated using a logistic regression model. P<0.05 was considered statistically significant. RESULTS: 66 patients were included with a median age of 68 years; 54 patients (82%) received ICI and 12 patients (18%) received PBC. Compared with HCs, AAbs against the cancer/testis antigens (CTAG1B, CTAG2, MAGEB18), HSPA1A, TP53, KRAS, and FGFR3 were significantly elevated in patients with mUC. AAbs against BRCA2, TP53, and CTNBB1 were associated with response, and those against BICD2 and UACA were associated with resistance to ICI therapy. AAbs against MITF, CDH3, and KDM4A were associated with development of irAEs in patient who received an ICI. A higher variance in pre-to-post treatment fold change in AAb levels was seen in patients treated with ICI versus PBC and was associated with response to ICI. CONCLUSIONS: This is the first report of comprehensive AAb profiling of patients with mUC and identified key AAbs that were elevated in patients with mUC versus HCs as well as AAbs associated with therapeutic response to ICI. These findings are hypothesis generating and further mechanistic studies evaluating humoral immunity in UC are required.

TGF-ß in the microenvironment induces a physiologically occurring immune-suppressive senescent state.

TGF-ß induces senescence in embryonic tissues. Whether TGF-ß in the hypoxic tumor microenvironment (TME) induces senescence in cancer and how the ensuing senescence-associated secretory phenotype (SASP) remodels the cellular TME to influence immune checkpoint inhibitor (ICI) responses are unknown. We show that TGF-ß induces a deeper senescent state under hypoxia than under normoxia; deep senescence correlates with the degree of E2F suppression and is marked by multinucleation, reduced reentry into proliferation, and a distinct 14-gene SASP. Suppressing TGF-ß signaling in tumors in an immunocompetent mouse lung cancer model abrogates endogenous senescent cells and suppresses the 14-gene SASP and immune infiltration. Untreated human lung cancers with a high 14-gene SASP display immunosuppressive immune infiltration. In a lung cancer clinical trial of ICIs, elevated 14-gene SASP is associated with increased senescence, TGF-ß and hypoxia signaling, and poor progression-free survival. Thus, TME-induced senescence may represent a naturally occurring state in cancer, contributing to an immune-suppressive phenotype associated with immune therapy resistance.

Genomic and transcriptomic analysis of checkpoint blockade response in advanced non-small cell lung cancer.

Anti-PD-1/PD-L1 agents have transformed the treatment landscape of advanced non-small cell lung cancer (NSCLC). To expand our understanding of the molecular features underlying response to checkpoint inhibitors in NSCLC, we describe here the first joint analysis of the Stand Up To Cancer-Mark Foundation cohort, a resource of whole exome and/or RNA sequencing from 393 patients with NSCLC treated with anti-PD-(L)1 therapy, along with matched clinical response annotation. We identify a number of associations between molecular features and outcome, including (1) favorable (for example, ATM altered) and unfavorable (for example, TERT amplified) genomic subgroups, (2) a prominent association between expression of inducible components of the immunoproteasome and response and (3) a dedifferentiated tumor-intrinsic subtype with enhanced response to checkpoint blockade. Taken together, results from this cohort demonstrate the complexity of biological determinants underlying immunotherapy outcomes and reinforce the discovery potential of integrative analysis within large, well-curated, cancer-specific cohorts.

Role of tumor infiltrating lymphocytes and spatial immune heterogeneity in sensitivity to PD-1 axis blockers in non-small cell lung cancer.

BACKGROUND: Tumor infiltrating lymphocytes (TILs) reflect adaptive antitumor immune responses in cancer and are generally associated with favorable prognosis. However, the relationships between TILs subsets and their spatial arrangement with clinical benefit from immune checkpoint inhibitors (ICI) in non-small cell lung cancer (NSCLC) remains less explored. METHODS: We used multiplexed quantitative immunofluorescence panels to determine the association of major TILs subpopulations, CD8+ cytotoxic T cells, CD4+ helper T cells and CD20+ B cells, and T cell exhaustion markers, programmed cell death protein-1 (PD-1),lymphocyte-activation gene 3 (LAG-3) and T cell immunoglobulin mucin-3 (TIM-3) with outcomes in a multi-institutional cohort of baseline tumor samples from 179 patients with NSCLC treated with ICI. The analysis of full-face tumor biopsies including numerous fields of view allowed a detailed spatial analysis and assessment of tumor immune heterogeneity using a multiparametric quadratic entropy metric (Rao's Q Index (RQI)). RESULTS: TILs were preferentially located in the stromal tissue areas surrounding tumor-cell nests and CD8+ T cells were the most abundant subset. Higher density of stromal CD8+ cytotoxic T cells was significantly associated with longer survival, and this effect was more prominent in programmed death ligand-1 (PD-L1) positive cases. The role of baseline T cell infiltration to stratify PD-L1 expressing cases was confirmed measuring the T cell receptor-burden in an independent NSCLC cohort studied with whole-exome DNA sequencing. High levels of LAG-3 on T cells or elevated RQI heterogeneity index were associated with worse survival in the cohort. CONCLUSION: Baseline T cell density and T cell exhaustion marker expression can stratify outcomes in PD-L1 positive patients with NSCLC treated with ICI. Spatial immune heterogeneity can be measured using the RQI and is associated with survival in NSCLC.

Association of High Tumor Mutation Burden in Non-Small Cell Lung Cancers With Increased Immune Infiltration and Improved Clinical Outcomes of PD-L1 Blockade Across PD-L1 Expression Levels.

Importance: Although tumor mutation burden (TMB) has been explored as a potential biomarker of immunotherapy efficacy in solid tumors, there still is a lack of consensus about the optimal TMB threshold that best discriminates improved outcomes of immune checkpoint inhibitor therapy among patients with non-small cell lung cancer (NSCLC). Objectives: To determine the association between increasing TMB levels and immunotherapy efficacy across clinically relevant programmed death ligand-1 (PD-L1) levels in patients with NSCLC. Design, Setting, and Participants: This multicenter cohort study included patients with advanced NSCLC treated with immunotherapy who received programmed cell death-1 (PD-1) or PD-L1 inhibition in the Dana-Farber Cancer Institute (DFCI), Memorial Sloan Kettering Cancer Center (MSKCC), and in the Stand Up To Cancer (SU2C)/Mark Foundation data sets. Clinicopathological and genomic data were collected from patients between September 2013 and September 2020. Data analysis was performed from November 2021 to February 2022. Exposures: Treatment with PD-1/PD-L1 inhibition without chemotherapy. Main Outcomes and Measures: Association of TMB levels with objective response rate (ORR), progression-free survival (PFS), and overall survival (OS). Results: In the entire cohort of 1552 patients with advanced NSCLC who received PD-1/PD-L1 blockade, the median (range) age was 66 (22-92) years, 830 (53.5%) were women, and 1347 (86.8%) had cancer with nonsquamous histologic profile. A regression tree modeling ORR as a function of TMB identified 2 TMB groupings in the discovery cohort (MSKCC), defined as low TMB (≤19.0 mutations per megabase) and high TMB (>19.0 mutations per megabase), which were associated with increasing improvements in ORR, PFS, and OS in the discovery cohort and in 2 independent cohorts (DFCI and SU2C/Mark Foundation). These TMB levels also were associated with significant improvements in outcomes of immunotherapy in each PD-L1 tumor proportion score subgroup of less than 1%, 1% to 49%, and 50% or higher. The ORR to PD-1/PD-L1 inhibition was as high as 57% in patients with high TMB and PD-L1 expression 50% or higher and as low as 8.7% in patients with low TMB and PD-L1 expression less than 1%. Multiplexed immunofluorescence and transcriptomic profiling revealed that high TMB levels were associated with increased CD8-positive, PD-L1-positive T-cell infiltration, increased PD-L1 expression on tumor and immune cells, and upregulation of innate and adaptive immune response signatures. Conclusions and Relevance: These findings suggest that increasing TMB levels are associated with immune cell infiltration and an inflammatory T-cell-mediated response, resulting in increased sensitivity to PD-1/PD-L1 blockade in NSCLC across PD-L1 expression subgroups.

Combined tumor and immune signals from genomes or transcriptomes predict outcomes of checkpoint inhibition in melanoma.

Immune checkpoint blockade (CPB) improves melanoma outcomes, but many patients still do not respond. Tumor mutational burden (TMB) and tumor-infiltrating T cells are associated with response, and integrative models improve survival prediction. However, integrating immune/tumor-intrinsic features using data from a single assay (DNA/RNA) remains underexplored. Here, we analyze whole-exome and bulk RNA sequencing of tumors from new and published cohorts of 189 and 178 patients with melanoma receiving CPB, respectively. Using DNA, we calculate T cell and B cell burdens (TCB/BCB) from rearranged TCR/Ig sequences and find that patients with TMBhigh and TCBhigh or BCBhigh have improved outcomes compared to other patients. By combining pairs of immune- and tumor-expressed genes, we identify three gene pairs associated with response and survival, which validate in independent cohorts. The top model includes lymphocyte-expressed MAP4K1 and tumor-expressed TBX3. Overall, RNA or DNA-based models combining immune and tumor measures improve predictions of melanoma CPB outcomes.

TRPV4 channels mediate cyclic strain-induced endothelial cell reorientation through integrin-to-integrin signaling.

Cyclic mechanical strain produced by pulsatile blood flow regulates the orientation of endothelial cells lining blood vessels and influences critical processes such as angiogenesis. Mechanical stimulation of stretch-activated calcium channels is known to mediate this reorientation response; however, the molecular basis remains unknown. Here, we show that cyclically stretching capillary endothelial cells adherent to flexible extracellular matrix substrates activates mechanosensitive TRPV4 (transient receptor potential vanilloid 4) ion channels that, in turn, stimulate phosphatidylinositol 3-kinase-dependent activation and binding of additional beta1 integrin receptors, which promotes cytoskeletal remodeling and cell reorientation. Inhibition of integrin activation using blocking antibodies and knock down of TRPV4 channels using specific small interfering RNA suppress strain-induced capillary cell reorientation. Thus, mechanical forces that physically deform extracellular matrix may guide capillary cell reorientation through a strain-dependent "integrin-to-integrin" signaling mechanism mediated by force-induced activation of mechanically gated TRPV4 ion channels on the cell surface.

Predictive Role of Computed Tomography Texture Analysis in Patients with Metastatic Urothelial Cancer Treated with Programmed Death-1 and Programmed Death-ligand 1 Inhibitors.

BACKGROUND: Reliable biomarkers to predict the response of metastatic urothelial cancer (mUC) to programmed death-1 and programmed death-ligand 1 (PD-1/PD-L1) inhibitors are being investigated. Texture analysis represents tumor heterogeneity and may serve as a predictor of response in mUC. OBJECTIVE: To assess the predictive ability of computed tomography (CT) texture analysis for progression-free survival (PFS) in patients with mUC treated with PD-1/PD-L1 inhibitors. DESIGN, SETTING, AND PARTICIPANTS: Forty-two postplatinum patients with mUC treated with PD-1/PD-L1 inhibitors from 2013 to 2018, including those with measurable disease per RECIST 1.1 who had contrast-enhanced baseline or first follow-up CT within 3mo after starting treatment, were included. PFS was calculated based on serial follow-up CT scans. Eleven patients with follow-up of <12mo without progression were excluded. Texture features of measurable lesions on baseline and first follow-up CT were extracted using commercially available software (TexRAD; Feedback Plc, Cambridge, UK) using different spatial scaling factors (0, 2-6). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Stepwise logistic regression analysis was conducted to identify patients with PFS <12mo, and performance was assessed using receiver operator characteristic curves. RESULTS AND LIMITATIONS: Of 31 included patients, 18 had PFS <12mo. Twenty-five baseline CT and 29 first follow-up CT scans met the inclusion criteria. In patients with PFS <12mo, entropy and mean were higher on first follow-up CT (p=0.02 and p=0.005, respectively). A predictive model including mean and entropy on first follow-up CT yielded 95% sensitivity, 80% specificity, 90% positive predictive value, 89% negative predictive value, and 90% accuracy (area under the curve=0.963) to identify patients with PFS <12mo. Limitations include retrospective nature and small sample size. CONCLUSIONS: CT texture analysis can help predict early progression with high accuracy soon after starting PD-1/PD-L1 inhibitors. Studies investigating the correlation of texture analysis with survival endpoints may help validate texture analysis as a biomarker of PD-1/PD-L1 inhibitors' treatment response. PATIENT SUMMARY: Computed tomography texture analysis can help predict durability of response in patients with metastatic urothelial cancer early during treatment with programmed death-1 and programmed death-ligand 1 (PD-1/PD-L1) inhibitors.

Polycythemia and Thrombocytosis.

Myeloproliferative neoplasms (MPNs) are diseases of excess cell proliferation from bone marrow precursors. Two classic MPNs, polycythemia vera (PV) and essential thrombocytosis (ET), are conditions of excess proliferation of red blood cells and platelets, respectively. Although PV and ET involve different cells in the myeloid lineage, their clinical presentations have shared features, consistent with overlapping mutations in growth factor signaling. The management of both diseases involves minimizing the risk of thrombotic and hemorrhagic complications. Both PV and ET can progress to myelofibrosis or acute myeloid leukemia, portending a poor prognosis. MPNs can also present as primary myelofibrosis.

Proliferation and tumorigenesis of a murine sarcoma cell line in the absence of DICER1.

MicroRNAs are a class of short ~22 nucleotide RNAs predicted to regulate nearly half of all protein coding genes, including many involved in basal cellular processes and organismal development. Although a global reduction in miRNAs is commonly observed in various human tumors, complete loss has not been documented, suggesting an essential function for miRNAs in tumorigenesis. Here we present the finding that transformed or immortalized Dicer1 null somatic cells can be isolated readily in vitro, maintain the characteristics of DICER1-expressing controls and remain stably proliferative. Furthermore, Dicer1 null cells from a sarcoma cell line, though depleted of miRNAs, are competent for tumor formation. Hence, miRNA levels in cancer may be maintained in vivo by a complex stabilizing selection in the intratumoral environment.

Computational model for early cardiac looping.

Looping is a vital event during early cardiac morphogenesis, as the initially straight heart tube bends and twists into a curved tube, laying out the basic pattern of the future four-chambered heart. Despite intensive study for almost a century, the biophysical mechanisms that drive this process are not well understood. To explore a recently proposed hypothesis for looping, we constructed a finite element model for the embryonic chick heart during the first phase of looping, called c-looping. The model includes the main structures of the early heart (heart tube, omphalomesenteric veins, and dorsal mesocardium), and the analysis features realistic three-dimensional geometry, nonlinear passive and active material properties, and anisotropic growth. As per our earlier hypothesis for c-looping, actin-based morpho-genetic processes (active cell shape change, cytoskeletal contraction, and cell migration) are simulated in specific regions of the model. The model correctly predicts the initial gross morphological shape changes of the heart, as well as distributions of morphogenetic stresses and strains measured in embryonic chick hearts. The model was tested further in studies that perturbed normal cardiac morphogenesis. The model, taken together with the new experimental data, supports our hypothesis for the mechanisms that drive early looping.