Innate Immune Landscape in Early Lung Adenocarcinoma by Paired Single-Cell Analyses.

To guide the design of immunotherapy strategies for patients with early stage lung tumors, we developed a multiscale immune profiling strategy to map the immune landscape of early lung adenocarcinoma lesions to search for tumor-driven immune changes. Utilizing a barcoding method that allows a simultaneous single-cell analysis of the tumor, non-involved lung, and blood cells, we provide a detailed immune cell atlas of early lung tumors. We show that stage I lung adenocarcinoma lesions already harbor significantly altered T cell and NK cell compartments. Moreover, we identified changes in tumor-infiltrating myeloid cell (TIM) subsets that likely compromise anti-tumor T cell immunity. Paired single-cell analyses thus offer valuable knowledge of tumor-driven immune changes, providing a powerful tool for the rational design of immune therapies. VIDEO ABSTRACT.

Host-Protozoan Interactions Protect from Mucosal Infections through Activation of the Inflammasome.

While conventional pathogenic protists have been extensively studied, there is an underappreciated constitutive protist microbiota that is an integral part of the vertebrate microbiome. The impact of these species on the host and their potential contributions to mucosal immune homeostasis remain poorly studied. Here, we show that the protozoan Tritrichomonas musculis activates the host epithelial inflammasome to induce IL-18 release. Epithelial-derived IL-18 promotes dendritic cell-driven Th1 and Th17 immunity and confers dramatic protection from mucosal bacterial infections. Along with its role as a "protistic" antibiotic, colonization with T. musculis exacerbates the development of T-cell-driven colitis and sporadic colorectal tumors. Our findings demonstrate a novel mutualistic host-protozoan interaction that increases mucosal host defenses at the cost of an increased risk of inflammatory disease.

A neoantigen fitness model predicts tumour response to checkpoint blockade immunotherapy.

Checkpoint blockade immunotherapies enable the host immune system to recognize and destroy tumour cells. Their clinical activity has been correlated with activated T-cell recognition of neoantigens, which are tumour-specific, mutated peptides presented on the surface of cancer cells. Here we present a fitness model for tumours based on immune interactions of neoantigens that predicts response to immunotherapy. Two main factors determine neoantigen fitness: the likelihood of neoantigen presentation by the major histocompatibility complex (MHC) and subsequent recognition by T cells. We estimate these components using the relative MHC binding affinity of each neoantigen to its wild type and a nonlinear dependence on sequence similarity of neoantigens to known antigens. To describe the evolution of a heterogeneous tumour, we evaluate its fitness as a weighted effect of dominant neoantigens in the subclones of the tumour. Our model predicts survival in anti-CTLA-4-treated patients with melanoma and anti-PD-1-treated patients with lung cancer. Importantly, low-fitness neoantigens identified by our method may be leveraged for developing novel immunotherapies. By using an immune fitness model to study immunotherapy, we reveal broad similarities between the evolution of tumours and rapidly evolving pathogens.

Clinical and Genomic Characterization of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS CoV-2) Infections in mRNA Vaccinated Health Care Personnel in New York City.

BACKGROUND: Vaccine-induced clinical protection against severe acute respiratory syndrome coronavirus 2 (SARS CoV-2) variants is an evolving target. There are limited genomic level data on SARS CoV-2 breakthrough infections and vaccine effectiveness (VE) since the global spread of the B.1.617.2 (Delta) variant. METHODS: In a retrospective study from 1 November 2020 to 31 August 2021, divided as pre-Delta and Delta-dominant periods, laboratory-confirmed SARS CoV-2 infections among healthcare personnel (HCP) at a large tertiary cancer center in New York City were examined to compare the weekly infection rate-ratio in vaccinated, partially vaccinated, and unvaccinated HCP. We describe the clinical and genomic epidemiologic features of post-vaccine infections to assess for selection of variants of concern (VOC)/variants of interest (VOI) in the early post-vaccine period and impact of B.1.617.2 (Delta) variant domination on VE. RESULTS: Among 13658 HCP in our cohort, 12379 received at least 1 dose of a messenger RNA (mRNA) vaccine. In the pre-Delta period overall VE was 94.5%. Whole genome sequencing (WGS) of 369 isolates in the pre-Delta period did not reveal a clade bias for VOC/VOI specific to post-vaccine infections. VE in the Delta dominant phase was 75.6%. No hospitalizations occurred among vaccinated HCP in the entire study period, compared to 17 hospitalizations and 1 death among unvaccinated HCP. CONCLUSIONS: Findings show high VE among HCP in New York City in the pre-Delta phase, with moderate decline in VE post-Delta emergence. SARS CoV-2 clades were similarly distributed among vaccinated and unvaccinated infected HCP without apparent clustering during the pre-Delta period of diverse clade circulation. Strong vaccine protection against hospitalization was maintained through the entire study period.

Centroid based clustering of high throughput sequencing reads based on n-mer counts.

BACKGROUND: Many problems in computational biology require alignment-free sequence comparisons. One of the common tasks involving sequence comparison is sequence clustering. Here we apply methods of alignment-free comparison (in particular, comparison using sequence composition) to the challenge of sequence clustering. RESULTS: We study several centroid based algorithms for clustering sequences based on word counts. Study of their performance shows that using k-means algorithm with or without the data whitening is efficient from the computational point of view. A higher clustering accuracy can be achieved using the soft expectation maximization method, whereby each sequence is attributed to each cluster with a specific probability. We implement an open source tool for alignment-free clustering. It is publicly available from github: https://github.com/luscinius/afcluster. CONCLUSIONS: We show the utility of alignment-free sequence clustering for high throughput sequencing analysis despite its limitations. In particular, it allows one to perform assembly with reduced resources and a minimal loss of quality. The major factor affecting performance of alignment-free read clustering is the length of the read.

Cancer cells co-evolve with retrotransposons to mitigate viral mimicry.

Overexpression of repetitive elements is an emerging hallmark of human cancers 1 . Diverse repeats can mimic viruses by replicating within the cancer genome through retrotransposition, or presenting pathogen-associated molecular patterns (PAMPs) to the pattern recognition receptors (PRRs) of the innate immune system 2-5 . Yet, how specific repeats affect tumor evolution and shape the tumor immune microenvironment (TME) in a pro- or anti-tumorigenic manner remains poorly defined. Here, we integrate whole genome and total transcriptome data from a unique autopsy cohort of multiregional samples collected in pancreatic ductal adenocarcinoma (PDAC) patients, into a comprehensive evolutionary analysis. We find that more recently evolved S hort I nterspersed N uclear E lements (SINE), a family of retrotransposable repeats, are more likely to form immunostimulatory double-strand RNAs (dsRNAs). Consequently, younger SINEs are strongly co-regulated with RIG-I like receptor associated type-I interferon genes but anti-correlated with pro-tumorigenic macrophage infiltration. We discover that immunostimulatory SINE expression in tumors is regulated by either L ong I nterspersed N uclear E lements 1 (LINE1/L1) mobility or ADAR1 activity in a TP53 mutation dependent manner. Moreover, L1 retrotransposition activity tracks with tumor evolution and is associated with TP53 mutation status. Altogether, our results suggest pancreatic tumors actively evolve to modulate immunogenic SINE stress and induce pro-tumorigenic inflammation. Our integrative, evolutionary analysis therefore illustrates, for the first time, how dark matter genomic repeats enable tumors to co-evolve with the TME by actively regulating viral mimicry to their selective advantage.

Prolonged SARS-CoV-2 Infection in Patients with Lymphoid Malignancies.

Coronavirus disease 2019 (COVID-19) infection results in both acute mortality and persistent and/or recurrent disease in patients with hematologic malignancies, but the drivers of persistent infection in this population are unknown. We found that B-cell lymphomas were at particularly high risk for persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) positivity. Further analysis of these patients identified discrete risk factors for initial disease severity compared with disease chronicity. Active therapy and diminished T-cell counts were drivers of acute mortality in COVID-19-infected patients with lymphoma. Conversely, B cell-depleting therapy was the primary driver of rehospitalization for COVID-19. In patients with persistent SARS-CoV-2 positivity, we observed high levels of viral entropy consistent with intrahost viral evolution, particularly in patients with impaired CD8+ T-cell immunity. These results suggest that persistent COVID-19 infection is likely to remain a risk in patients with impaired adaptive immunity and that additional therapeutic strategies are needed to enable viral clearance in this high-risk population. SIGNIFICANCE: We describe the largest cohort of persistent symptomatic COVID-19 infection in patients with lymphoid malignancies and identify B-cell depletion as the key immunologic driver of persistent infection. Furthermore, we demonstrate ongoing intrahost viral evolution in patients with persistent COVID-19 infection, particularly in patients with impaired CD8+ T-cell immunity.This article is highlighted in the In This Issue feature, p. 1.

Impaired humoral immunity is associated with prolonged COVID-19 despite robust CD8 T cell responses.

How immune dysregulation affects recovery from COVID-19 infection in patients with cancer remains unclear. We analyzed cellular and humoral immune responses in 103 patients with prior COVID-19 infection, more than 20% of whom had delayed viral clearance. Delayed clearance was associated with loss of antibodies to nucleocapsid and spike proteins with a compensatory increase in functional T cell responses. High-dimensional analysis of peripheral blood samples demonstrated increased CD8+ effector T cell differentiation and a broad but poorly converged COVID-specific T cell receptor (TCR) repertoire in patients with prolonged disease. Conversely, patients with a CD4+ dominant immunophenotype had a lower incidence of prolonged disease and exhibited a deep and highly select COVID-associated TCR repertoire, consistent with effective viral clearance and development of T cell memory. These results highlight the importance of B cells and CD4+ T cells in promoting durable SARS-CoV-2 clearance and the significance of coordinated cellular and humoral immunity for long-term disease control.

Y RNAs are conserved endogenous RIG-I ligands across RNA virus infection and are targeted by HIV-1.

Pattern recognition receptors (PRRs) protect against microbial invasion by detecting specific molecular patterns found in pathogens and initiating an immune response. Although microbial-derived PRR ligands have been extensively characterized, the contribution and relevance of endogenous ligands to PRR activation remains overlooked. Here, we characterize the landscape of endogenous ligands that engage RIG-I-like receptors (RLRs) upon infection by different RNA viruses. In each infection, several RNAs transcribed by RNA polymerase III (Pol3) specifically engaged RLRs, particularly the family of Y RNAs. Sensing of Y RNAs was dependent on their mimicking of viral secondary structure and their 5'-triphosphate extremity. Further, we found that HIV-1 triggered a VPR-dependent downregulation of RNA triphosphatase DUSP11 in vitro and in vivo, inducing a transcriptome-wide change of cellular RNA 5'-triphosphorylation that licenses Y RNA immunogenicity. Overall, our work uncovers the contribution of endogenous RNAs to antiviral immunity and demonstrates the importance of this pathway in HIV-1 infection.

Satellite repeat RNA expression in epithelial ovarian cancer associates with a tumor-immunosuppressive phenotype.

Aberrant expression of viral-like repeat elements is a common feature of epithelial cancers, and the substantial diversity of repeat species provides a distinct view of the cancer transcriptome. Repeatome profiling across ovarian, pancreatic, and colorectal cell lines identifies distinct clustering independent of tissue origin that is seen with coding gene analysis. Deeper analysis of ovarian cancer cell lines demonstrated that human satellite II (HSATII) satellite repeat expression was highly associated with epithelial-mesenchymal transition (EMT) and anticorrelated with IFN-response genes indicative of a more aggressive phenotype. SATII expression - and its correlation with EMT and anticorrelation with IFN-response genes - was also found in ovarian cancer RNA-Seq data and was associated with significantly shorter survival in a second independent cohort of patients with ovarian cancer. Repeat RNAs were enriched in tumor-derived extracellular vesicles capable of stimulating monocyte-derived macrophages, demonstrating a mechanism that alters the tumor microenvironment with these viral-like sequences. Targeting of HSATII with antisense locked nucleic acids stimulated IFN response and induced MHC I expression in ovarian cancer cell lines, highlighting a potential strategy of modulating the repeatome to reestablish antitumor cell immune surveillance.

Reverse Transcriptase Inhibition Disrupts Repeat Element Life Cycle in Colorectal Cancer.

Altered RNA expression of repetitive sequences and retrotransposition are frequently seen in colorectal cancer, implicating a functional importance of repeat activity in cancer progression. We show the nucleoside reverse transcriptase inhibitor 3TC targets activities of these repeat elements in colorectal cancer preclinical models with a preferential effect in p53-mutant cell lines linked with direct binding of p53 to repeat elements. We translate these findings to a human phase II trial of single-agent 3TC treatment in metastatic colorectal cancer with demonstration of clinical benefit in 9 of 32 patients. Analysis of 3TC effects on colorectal cancer tumorspheres demonstrates accumulation of immunogenic RNA:DNA hybrids linked with induction of interferon response genes and DNA damage response. Epigenetic and DNA-damaging agents induce repeat RNAs and have enhanced cytotoxicity with 3TC. These findings identify a vulnerability in colorectal cancer by targeting the viral mimicry of repeat elements. SIGNIFICANCE: Colorectal cancers express abundant repeat elements that have a viral-like life cycle that can be therapeutically targeted with nucleoside reverse transcriptase inhibitors (NRTI) commonly used for viral diseases. NRTIs induce DNA damage and interferon response that provide a new anticancer therapeutic strategy. This article is highlighted in the In This Issue feature, p. 1397.

Molecular characterization of severe and mild cases of influenza A (H1N1) 2009 strain from Argentina.

While worldwide pandemic influenza A(H1N1) pdm case fatality rate (CFR) was 0.4%, Argentina's was 4.5%. A total of 34 strains from mild and severe cases were analyzed. A full genome sequencing was carried out on 26 of these, and a partial sequencing on the remaining eight. We observed no evidence that the high CFR can be attributed to direct virus changes. No evidence of re-assortment, mutations associated with resistance to antiviral drugs, or genetic drift that might contribute to virulence was observed. Although the mutation D225G associated with severity in the latest reports from the Ukraine and Norway is not observed among the Argentine strains, an amino acid change in the area (S206T) surrounding the HA receptor binding domain was observed, the same previously established worldwide.

An integrative approach identifies dysregulated long non-coding RNAs as microRNA decoys during nevus to melanoma transformation.

Mounting evidence supports a role for dysregulated long non-coding RNAs (lncRNA) in the development of many cancers. A recently discovered function of lncRNAs is to act as microRNA (miR) decoys or competing endogenous RNAs, which sequester specific miRs and relieve negative regulation of mRNA expression by miRs. Although a large number of non-coding RNAs are thought to function as competing endogenous RNAs, miR-sequestering lncRNAs involved in nevus to melanoma transformation remain largely unknown. In this study, we applied a bioinformatics approach to a unique dataset of benign melanocytic nevi and primary melanomas of the skin in order to fill this research gap. We modified a previously published miR target prediction algorithm, RNAhybrid, and improved its search efficiency. We reported the presence of many lncRNAs and miRs deregulated when transitioning from a senescence-like state of nevi to melanoma. We provided evidence of a relatively new and understudied mechanism of gene regulation during this process and identified for the first time lncRNAs (n = 122) that may potentially function as miR decoys as well as their target miRs during nevus to melanoma transformation. The knowledge presented here can be employed for developing biomarkers for diagnostic and risk stratification purposes.

Mutation-derived Neoantigen-specific T-cell Responses in Multiple Myeloma.

PURPOSE: Somatic mutations in cancer cells can give rise to novel protein sequences that can be presented by antigen-presenting cells as neoantigens to the host immune system. Tumor neoantigens represent excellent targets for immunotherapy, due to their specific expression in cancer tissue. Despite the widespread use of immunomodulatory drugs and immunotherapies that recharge T and NK cells, there has been no direct evidence that neoantigen-specific T-cell responses are elicited in multiple myeloma. EXPERIMENTAL DESIGN: Using next-generation sequencing data we describe the landscape of neo-antigens in 184 patients with multiple myeloma and successfully validate neoantigen-specific T cells in patients with multiple myeloma and support the feasibility of neoantigen-based therapeutic vaccines for use in cancers with intermediate mutational loads such as multiple myeloma. RESULTS: In this study, we demonstrate an increase in neoantigen load in relapsed patients with multiple myeloma as compared with newly diagnosed patients with multiple myeloma. Moreover, we identify shared neoantigens across multiple patients in three multiple myeloma oncogenic driver genes (KRAS, NRAS, and IRF4). Next, we validate neoantigen T-cell response and clonal expansion in correlation with clinical response in relapsed patients with multiple myeloma. This is the first study to experimentally validate the immunogenicity of predicted neoantigens from next-generation sequencing in relapsed patients with multiple myeloma. CONCLUSIONS: Our findings demonstrate that somatic mutations in multiple myeloma can be immunogenic and induce neoantigen-specific T-cell activation that is associated with antitumor activity in vitro and clinical response in vivo. Our results provide the foundation for using neoantigen targeting strategies such as peptide vaccines in future trials for patients with multiple myeloma.

Temporal and Spatial Heterogeneity of Host Response to SARS-CoV-2 Pulmonary Infection.

The relationship of SARS-CoV-2 lung infection and severity of pulmonary disease is not fully understood. We analyzed autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter- and intra- patient heterogeneity of pulmonary virus infection. There was a spectrum of high and low virus cases that was associated with duration of disease and activation of interferon pathway genes. Using a digital spatial profiling platform, the virus corresponded to distinct spatial expression of interferon response genes and immune checkpoint genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

Temporal and spatial heterogeneity of host response to SARS-CoV-2 pulmonary infection.

The relationship of SARS-CoV-2 pulmonary infection and severity of disease is not fully understood. Here we show analysis of autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter-patient and intra-patient heterogeneity of pulmonary virus infection. There is a spectrum of high and low virus cases associated with duration of disease. High viral cases have high activation of interferon pathway genes and a predominant M1-like macrophage infiltrate. Low viral cases are more heterogeneous likely reflecting inherent patient differences in the evolution of host response, but there is consistent indication of pulmonary epithelial cell recovery based on napsin A immunohistochemistry and RNA expression of surfactant and mucin genes. Using a digital spatial profiling platform, we find the virus corresponds to distinct spatial expression of interferon response genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

A tumor-specific endogenous repetitive element is induced by herpesviruses.

Tandem satellite repeats account for 3% of the human genome. One of them, Human Satellite II (HSATII), is highly expressed in several epithelial cancers and cancer cell lines. Here we report an acute induction of HSATII RNA in human cells infected with two herpes viruses. We show that human cytomegalovirus (HCMV) IE1 and IE2 proteins cooperate to induce HSATII RNA affecting several aspects of the HCMV replication cycle, viral titers and infected-cell processes. HSATII RNA expression in tissue from two chronic HCMV colitis patients correlates with the strength of CMV antigen staining. Thus, endogenous HSATII RNA synthesis after herpesvirus infections appears to have functionally important consequences for viral replication and may provide a novel insight into viral pathogenesis. The HSATII induction seen in both infected and cancer cells suggests possible convergence upon common HSATII-based regulatory mechanisms in these seemingly disparate diseases.

Global Cancer Transcriptome Quantifies Repeat Element Polarization between Immunotherapy Responsive and T Cell Suppressive Classes.

It has been posited that anti-tumoral innate activation is driven by derepression of endogenous repeats. We compared RNA sequencing protocols to assess repeat transcriptomes in The Cancer Genome Atlas (TCGA). Although poly(A) selection efficiently detects coding genes, most non-coding genes, and limited subsets of repeats, it fails to capture overall repeat expression and co-expression. Alternatively, total RNA expression reveals distinct repeat co-expression subgroups and delivers greater dynamic changes, implying they may serve as better biomarkers of clinical outcomes. We show that endogenous retrovirus expression predicts immunotherapy response better than conventional immune signatures in one cohort yet is not predictive in another. Moreover, we find that global repeat derepression, including the HSATII satellite repeat, correlates with an immunosuppressive phenotype in colorectal and pancreatic tumors and validate in situ. In conclusion, we stress the importance of analyzing the full spectrum of repeat transcription to decode their role in tumor immunity.

FBXW7 regulates a mitochondrial transcription program by modulating MITF.

FBXW7 is well characterized as a tumor suppressor in many human cancers including melanoma; however, the mechanisms of tumor-suppressive function have not been fully elucidated. We leveraged two distinct RNA sequencing datasets: human melanoma cell lines (n = 10) with control versus silenced FBXW7 and a cohort of human melanoma tumor samples (n = 51) to define the transcriptomic fingerprint regulated by FBXW7. Here, we report that loss of FBXW7 enhances a mitochondrial gene transcriptional program that is dependent on MITF in human melanoma and confers poor patient outcomes. MITF is a lineage-specific master regulator of melanocytes and together with PGC-1alpha is a marker for melanoma subtypes with dependence for mitochondrial oxidative metabolism. We found that inactivation of FBXW7 elevates MITF protein levels in melanoma cells. In vitro studies examining loss of FBXW7 and MITF alone or in combination showed that FBXW7 is an upstream regulator for the MITF/PGC-1 signaling.

Diverse repetitive element RNA expression defines epigenetic and immunologic features of colon cancer.

There is tremendous excitement for the potential of epigenetic therapies in cancer, but the ability to predict and monitor response to these drugs remains elusive. This is in part due to the inability to differentiate the direct cytotoxic and the immunomodulatory effects of these drugs. The DNA-hypomethylating agent 5-azacitidine (AZA) has shown these distinct effects in colon cancer and appears to be linked to the derepression of repeat RNAs. LINE and HERV are two of the largest classes of repeats in the genome, and despite many commonalities, we found that there is heterogeneity in behavior among repeat subtypes. Specifically, the LINE-1 and HERV-H subtypes detected by RNA sequencing and RNA in situ hybridization in colon cancers had distinct expression patterns, which suggested that these repeats are correlated to transcriptional programs marking different biological states. We found that low LINE-1 expression correlates with global DNA hypermethylation, wild-type TP53 status, and responsiveness to AZA. HERV-H repeats were not concordant with LINE-1 expression but were found to be linked with differences in FOXP3+ Treg tumor infiltrates. Together, distinct repeat RNA expression patterns define new molecular classifications of colon cancer and provide biomarkers that better distinguish cytotoxic from immunomodulatory effects by epigenetic drugs.