High-throughput RNA isoform sequencing using programmed cDNA concatenation

Full-length RNA-sequencing methods using long-read technologies can capture complete transcript isoforms, but their throughput is limited. We introduce multiplexed arrays isoform sequencing (MAS-ISO-seq), a technique for programmably concatenating complementary DNAs (cDNAs) into molecules optimal for long-read sequencing, increasing the throughput >15-fold to nearly 40 million cDNA reads per run on the Sequel IIe sequencer. When applied to single-cell RNA sequencing of tumor-infiltrating T cells, MAS-ISO-seq demonstrated a 12- to 32-fold increase in the discovery of differentially spliced genes.

CXCL9:SPP1 macrophage polarity identifies a network of cellular programs that control human cancers

Tumor microenvironments (TMEs) influence cancer progression but are complex and often differ between patients. Considering that microenvironment variations may reveal rules governing intratumoral cellular programs and disease outcome, we focused on tumor-to-tumor variation to examine 52 head and neck squamous cell carcinomas. We found that macrophage polarity-defined by CXCL9 and SPP1 (CS) expression but not by conventional M1 and M2 markers-had a noticeably strong prognostic association. CS macrophage polarity also identified a highly coordinated network of either pro- or antitumor variables, which involved each tumor-associated cell type and was spatially organized. We extended these findings to other cancer indications. Overall, these results suggest that, despite their complexity, TMEs coordinate coherent responses that control human cancers and for which CS macrophage polarity is a relevant yet simple variable.

Longitudinal characterization of circulating neutrophils uncovers phenotypes associated with severity in hospitalized COVID-19 patients

Mechanisms of neutrophil involvement in severe coronavirus disease 2019 (COVID-19) remain incompletely understood. Here, we collect longitudinal blood samples from 306 hospitalized COVID-19+ patients and 86 controls and perform bulk RNA sequencing of enriched neutrophils, plasma proteomics, and high-throughput antibody profiling to investigate relationships between neutrophil states and disease severity. We identify dynamic switches between six distinct neutrophil subtypes. At days 3 and 7 post-hospitalization, patients with severe disease display a granulocytic myeloid-derived suppressor cell-like gene expression signature, while patients with resolving disease show a neutrophil progenitor-like signature. Humoral responses are identified as potential drivers of neutrophil effector functions, with elevated severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-specific immunoglobulin G1 (IgG1)-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirm that while patient-derived IgG antibodies induce phagocytosis in healthy donor neutrophils, IgA antibodies predominantly induce neutrophil cell death. Overall, our study demonstrates a dysregulated myelopoietic response in severe COVID-19 and a potential role for IgA-dominant responses contributing to mortality.

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.

Longitudinal characterization of circulating neutrophils uncovers distinct phenotypes associated with disease severity in hospitalized COVID-19 patients

Multiple studies have identified an association between neutrophils and COVID-19 disease severity; however, the mechanistic basis of this association remains incompletely understood. Here we collected 781 longitudinal blood samples from 306 hospitalized COVID-19 ⁺ patients, 78 COVID-19 ^âˆ' acute respiratory distress syndrome patients, and 8 healthy controls, and performed bulk RNA-sequencing of enriched neutrophils, plasma proteomics, cfDNA measurements and high throughput antibody profiling assays to investigate the relationship between neutrophil states and disease severity or death. We identified dynamic switches between six distinct neutrophil subtypes using non-negative matrix factorization (NMF) clustering. At days 3 and 7 post-hospitalization, patients with severe disease had an enrichment of a granulocytic myeloid derived suppressor cell-like state gene expression signature, while non-severe patients with resolved disease were enriched for a progenitor-like immature neutrophil state signature. Severe disease was associated with gene sets related to neutrophil degranulation, neutrophil extracellular trap (NET) signatures, distinct metabolic signatures, and enhanced neutrophil activation and generation of reactive oxygen species (ROS). We found that the majority of patients had a transient interferon-stimulated gene signature upon presentation to the emergency department (ED) defined here as Day 0, regardless of disease severity, which persisted only in patients who subsequently died. Humoral responses were identified as potential drivers of neutrophil effector functions, as enhanced antibody-dependent neutrophil phagocytosis and reduced NETosis was associated with elevated SARS-CoV-2-specific IgG1-to-IgA1 ratios in plasma of severe patients who survived. In vitro experiments confirmed that while patient-derived IgG antibodies mostly drove neutrophil phagocytosis and ROS production in healthy donor neutrophils, patient-derived IgA antibodies induced a predominant NETosis response. Overall, our study demonstrates neutrophil dysregulation in severe COVID-19 and a potential role for IgA-dominant responses in driving neutrophil effector functions in severe disease and mortality.