A specific and portable gene expression program underlies antigen archiving by lymphatic endothelial cells.

Antigens from protein subunit vaccination traffic from the tissue to the draining lymph node, either passively via the lymph or carried by dendritic cells at the local injection site. Lymph node (LN) lymphatic endothelial cells (LEC) actively acquire and archive foreign antigens, and archived antigen can be released during subsequent inflammatory stimulus to improve immune responses. Here, we answer questions about how LECs achieve durable antigen archiving and whether there are transcriptional signatures associated with LECs containing high levels of antigen. We used single cell sequencing in dissociated LN tissue to quantify antigen levels in LEC and dendritic cell populations at multiple timepoints after immunization, and used machine learning to define a unique transcriptional program within archiving LECs that can predict LEC archiving capacity in independent data sets. Finally, we validated this modeling, showing we could predict antigen archiving from a transcriptional dataset of CHIKV infected mice and demonstrated in vivo the accuracy of our prediction. Collectively, our findings establish a unique transcriptional program in LECs that promotes antigen archiving that can be translated to other systems.

Immunization-induced antigen archiving enhances local memory CD8+ T cell responses following an unrelated viral infection.

Antigens from viruses or immunizations can persist or are archived in lymph node stromal cells such as lymphatic endothelial cells (LEC) and fibroblastic reticular cells (FRC). Here, we find that, during the time frame of antigen archiving, LEC apoptosis caused by a second, but unrelated, innate immune stimulus such as vaccina viral infection or CpG DNA administration resulted in cross-presentation of archived antigens and boosted memory CD8 + T cells specific to the archived antigen. In contrast to "bystander" activation associated with unrelated infections, the memory CD8 + T cells specific to the archived antigen from the immunization were significantly higher than memory CD8 + T cells of a different antigen specificity. Finally, the boosted memory CD8 + T cells resulted in increased protection against Listeria monocytogenes expressing the antigen from the immunization, but only for the duration that the antigen was archived. These findings outline an important mechanism by which lymph node stromal cell archived antigens, in addition to bystander activation, can augment memory CD8 + T cell responses during repeated inflammatory insults.

Chikungunya virus infection disrupts lymph node lymphatic endothelial cell composition and function via MARCO.

Infection with chikungunya virus (CHIKV) causes disruption of draining lymph node (dLN) organization, including paracortical relocalization of B cells, loss of the B cell-T cell border, and lymphocyte depletion that is associated with infiltration of the LN with inflammatory myeloid cells. Here, we found that, during the first 24 hours of infection, CHIKV RNA accumulated in MARCO-expressing lymphatic endothelial cells (LECs) in both the floor and medullary LN sinuses. The accumulation of viral RNA in the LN was associated with a switch to an antiviral and inflammatory gene expression program across LN stromal cells, and this inflammatory response - including recruitment of myeloid cells to the LN - was accelerated by CHIKV-MARCO interactions. As CHIKV infection progressed, both floor and medullary LECs diminished in number, suggesting further functional impairment of the LN by infection. Consistent with this idea, antigen acquisition by LECs, a key function of LN LECs during infection and immunization, was reduced during pathogenic CHIKV infection.

Systematic analysis of nonsense variants uncovers peptide release rate as a novel modifier of nonsense-mediated mRNA decay efficiency.

Nonsense variants underlie many genetic diseases. The phenotypic impact of nonsense variants is determined by Nonsense-mediated mRNA decay (NMD), which degrades transcripts with premature termination codons (PTCs). NMD activity varies across transcripts and cellular contexts via poorly understood mechanisms. Here, by leveraging human genetic datasets, we uncover that the amino acid preceding the PTC dramatically affects NMD activity in human cells. We find that glycine codons in particular support high levels of NMD and are enriched before PTCs but depleted before normal termination codons (NTCs). Gly-PTC enrichment is most pronounced in human genes that tolerate loss-of-function variants. This suggests a strong biological impact for Gly-PTC in ensuring robust elimination of potentially toxic truncated proteins from non-essential genes. Biochemical assays revealed that the peptide release rate during translation termination is highly dependent on the identity of the amino acid preceding the stop codon. This release rate is the most critical feature determining NMD activity across our massively parallel reporter assays. Together, we conclude that NMD activity is significantly modulated by the "window of opportunity" offered by translation termination kinetics. Integrating the window of opportunity model with the existing framework of NMD would enable more accurate nonsense variant interpretation in the clinic.

Chikungunya virus infection disrupts lymph node lymphatic endothelial cell composition and function via MARCO.

Infection with chikungunya virus (CHIKV) causes disruption of draining lymph node (dLN) organization, including paracortical relocalization of B cells, loss of the B cell-T cell border, and lymphocyte depletion that is associated with infiltration of the LN with inflammatory myeloid cells. Here, we find that during the first 24 h of infection, CHIKV RNA accumulates in MARCO-expressing lymphatic endothelial cells (LECs) in both the floor and medullary LN sinuses. The accumulation of viral RNA in the LN was associated with a switch to an antiviral and inflammatory gene expression program across LN stromal cells, and this inflammatory response, including recruitment of myeloid cells to the LN, was accelerated by CHIKV-MARCO interactions. As CHIKV infection progressed, both floor and medullary LECs diminished in number, suggesting further functional impairment of the LN by infection. Consistent with this idea, we find that antigen acquisition by LECs, a key function of LN LECs during infection and immunization, was reduced during pathogenic CHIKV infection.

Vaccine adjuvant-elicited CD8+ T cell immunity is co-dependent on T-bet and FOXO1.

T-bet and FOXO1 are transcription factors canonically associated with effector and memory T cell fates, respectively. During an infectious response, these factors direct the development of CD8+ T cell fates, where T-bet deficiency leads to ablation of only short-lived effector cells, while FOXO1 deficiency results in selective loss of memory. In contrast, following adjuvanted subunit vaccination in mice, both effector- and memory-fated T cells are compromised in the absence of either T-bet or FOXO1. Thus, unlike responses to challenge with Listeria monocytogenes, productive CD8+ T cell responses to adjuvanted vaccination require coordinated regulation of FOXO1 and T-bet transcriptional programs. Single-cell RNA sequencing analysis confirms simultaneous T-bet, FOXO1, and TCF1 transcriptional activity in vaccine-elicited, but not infection-elicited, T cells undergoing clonal expansion. Collectively, our data show that subunit vaccine adjuvants elicit T cell responses dependent on transcription factors associated with effector and memory cell fates.

Nanopore sequencing of internal 2'-PO4 modifications installed by RNA repair.

Ligation by plant and fungal RNA ligases yields an internal 2'-phosphate group on each RNA ligation product. In budding yeast, this covalent mark occurs at the splice junction of two targets of ligation: intron-containing tRNAs and the messenger RNA HAC1 The repertoire of RNA molecules repaired by RNA ligation has not been explored due to a lack of unbiased approaches for identifying RNA ligation products. Here, we define several unique signals produced by 2'-phosphorylated RNAs during nanopore sequencing. A 2'-phosphate at the splice junction of HAC1 mRNA inhibits 5' → 3' degradation, enabling detection of decay intermediates in yeast RNA repair mutants by nanopore sequencing. During direct RNA sequencing, intact 2'-phosphorylated RNAs on HAC1 and tRNAs produce diagnostic changes in nanopore current properties and base calling features, including stalls produced as the modified RNA translocates through the nanopore motor protein. These approaches enable directed studies to identify novel RNA repair events in other contexts.

Modification mapping by nanopore sequencing.

Next generation sequencing (NGS) has provided biologists with an unprecedented view into biological processes and their regulation over the past 2 decades, fueling a wave of development of high throughput methods based on short read DNA and RNA sequencing. For nucleic acid modifications, NGS has been coupled with immunoprecipitation, chemical treatment, enzymatic treatment, and/or the use of reverse transcriptase enzymes with fortuitous activities to enrich for and to identify covalent modifications of RNA and DNA. However, the majority of nucleic acid modifications lack commercial monoclonal antibodies, and mapping techniques that rely on chemical or enzymatic treatments to manipulate modification signatures add additional technical complexities to library preparation. Moreover, such approaches tend to be specific to a single class of RNA or DNA modification, and generate only indirect readouts of modification status. Third generation sequencing technologies such as the commercially available "long read" platforms from Pacific Biosciences and Oxford Nanopore Technologies are an attractive alternative for high throughput detection of nucleic acid modifications. While the former can indirectly sense modified nucleotides through changes in the kinetics of reverse transcription reactions, nanopore sequencing can in principle directly detect any nucleic acid modification that produces a signal distortion as the nucleic acid passes through a nanopore sensor embedded within a charged membrane. To date, more than a dozen endogenous DNA and RNA modifications have been interrogated by nanopore sequencing, as well as a number of synthetic nucleic acid modifications used in metabolic labeling, structure probing, and other emerging applications. This review is intended to introduce the reader to nanopore sequencing and key principles underlying its use in direct detection of nucleic acid modifications in unamplified DNA or RNA samples, and outline current approaches for detecting and quantifying nucleic acid modifications by nanopore sequencing. As this technology matures, we anticipate advances in both sequencing chemistry and analysis methods will lead to rapid improvements in the identification and quantification of these epigenetic marks.

Neoplastic and immune single-cell transcriptomics define subgroup-specific intra-tumoral heterogeneity of childhood medulloblastoma.

BACKGROUND: Medulloblastoma (MB) is a heterogeneous disease in which neoplastic cells and associated immune cells contribute to disease progression. We aimed to determine the influence of neoplastic and immune cell diversity on MB biology in patient samples and animal models. METHODS: To better characterize cellular heterogeneity in MB we used single-cell RNA sequencing, immunohistochemistry, and deconvolution of transcriptomic data to profile neoplastic and immune populations in patient samples and animal models across childhood MB subgroups. RESULTS: Neoplastic cells cluster primarily according to individual sample of origin which is influenced by chromosomal copy number variance. Harmony alignment reveals novel MB subgroup/subtype-associated subpopulations that recapitulate neurodevelopmental processes, including photoreceptor and glutamatergic neuron-like cells in molecular subgroups GP3 and GP4, and a specific nodule-associated neuronally differentiated subpopulation in the sonic hedgehog subgroup. We definitively chart the spectrum of MB immune cell infiltrates, which include subpopulations that recapitulate developmentally related neuron-pruning and antigen-presenting myeloid cells. MB cellular diversity matching human samples is mirrored in subgroup-specific mouse models of MB. CONCLUSIONS: These findings provide a clearer understanding of the diverse neoplastic and immune cell subpopulations that constitute the MB microenvironment.

MARCO+ lymphatic endothelial cells sequester arthritogenic alphaviruses to limit viremia and viral dissemination.

Viremia in the vertebrate host is a major determinant of arboviral reservoir competency, transmission efficiency, and disease severity. However, immune mechanisms that control arboviral viremia are poorly defined. Here, we identify critical roles for the scavenger receptor MARCO in controlling viremia during arthritogenic alphavirus infections in mice. Following subcutaneous inoculation, arthritogenic alphavirus particles drain via the lymph and are rapidly captured by MARCO+ lymphatic endothelial cells (LECs) in the draining lymph node (dLN), limiting viral spread to the bloodstream. Upon reaching the bloodstream, alphavirus particles are cleared from the circulation by MARCO-expressing Kupffer cells in the liver, limiting viremia and further viral dissemination. MARCO-mediated accumulation of alphavirus particles in the draining lymph node and liver is an important host defense mechanism as viremia and viral tissue burdens are elevated in MARCO-/- mice and disease is more severe. In contrast to prior studies implicating a key role for lymph node macrophages in limiting viral dissemination, these findings exemplify a previously unrecognized arbovirus-scavenging role for lymphatic endothelial cells and improve our mechanistic understanding of viremia control during arthritogenic alphavirus infection.

Cell-level metadata are indispensable for documenting single-cell sequencing datasets.

Single-cell RNA sequencing (scRNA-seq) provides an unprecedented view of cellular diversity of biological systems. However, across the thousands of publications and datasets generated using this technology, we estimate that only a minority (<25%) of studies provide cell-level metadata information containing identified cell types and related findings of the published dataset. Metadata omission hinders reproduction, exploration, validation, and knowledge transfer and is a common problem across journals, data repositories, and publication dates. We encourage investigators, reviewers, journals, and data repositories to improve their standards and ensure proper documentation of these valuable datasets.

Molecular tracking devices quantify antigen distribution and archiving in the murine lymph node.

The detection of foreign antigens in vivo has relied on fluorescent conjugation or indirect read-outs such as antigen presentation. In our studies, we found that these widely used techniques had several technical limitations that have precluded a complete picture of antigen trafficking or retention across lymph node cell types. To address these limitations, we developed a 'molecular tracking device' to follow the distribution, acquisition, and retention of antigen in the lymph node. Utilizing an antigen conjugated to a nuclease-resistant DNA tag, acting as a combined antigen-adjuvant conjugate, and single-cell mRNA sequencing, we quantified antigen abundance in the lymph node. Variable antigen levels enabled the identification of caveolar endocytosis as a mechanism of antigen acquisition or retention in lymphatic endothelial cells. Thus, these molecular tracking devices enable new approaches to study dynamic tissue dissemination of antigen-adjuvant conjugates and identify new mechanisms of antigen acquisition and retention at cellular resolution in vivo.

The lymphatic system is a network of ducts that transports fluid, proteins, and immune cells from different organs around the body. Lymph nodes provide pit stops at hundreds of points along this network where immune cells reside, and lymph fluid can be filtered and cleaned. When pathogens, such as viruses or bacteria, enter the body during an infection, fragments of their proteins can get swept into the lymph nodes. These pathogenic proteins or protein fragments activate resident immune cells and kickstart the immune response. Vaccines are designed to mimic this process by introducing isolated pathogenic proteins in a controlled way to stimulate similar immune reactions in lymph nodes. Once an infection has been cleared by the immune system, or a vaccination has triggered the immune system, most pathogenic proteins get cleared away. However, a small number of pathogenic proteins remain in the lymph nodes to enable immune cells to respond more strongly and quickly the next time they see the same pathogen. Yet it is largely unclear how much protein remains for training and how or where it is all stored. Current techniques are not sensitive or long-lived enough to accurately detect and track these small protein deposits over time. Walsh, Sheridan, Lucas, et al. have addressed this problem by developing biological tags that can be attached to the pathogenic proteins so they can be traced. These tags were designed so the body cannot easily break them down, helping them last as long as the proteins they are attached to. Walsh, Sheridan, Lucas et al. tested whether vaccinating mice with the tagged proteins allowed the proteins to be tracked. The method they used was designed to identify individual cell types based on their genetic information along with the tag. This allowed them to accurately map the complex network of cells involved in storing and retrieving archived protein fragments, as well as those involved in training new immune cells to recognize them. These results provide important insights into the protein archiving system that is involved in enhancing immune memory. This may help guide the development of new vaccination strategies that can manipulate how proteins are archived to establish more durable immune protection. The biological tags developed could also be used to track therapeutic proteins, allowing scientists to determine how long cancer drugs, antibody therapies or COVID19 anti-viral agents remain in the body. This information could then be used by doctors to plan specific and personalized treatment timetables for patients.

Physiologic RNA targets and refined sequence specificity of coronavirus EndoU.

Coronavirus EndoU inhibits dsRNA-activated antiviral responses; however, the physiologic RNA substrates of EndoU are unknown. In this study, we used mouse hepatitis virus (MHV)-infected bone marrow-derived macrophage (BMM) and cyclic phosphate cDNA sequencing to identify the RNA targets of EndoU. EndoU targeted viral RNA, cleaving the 3' side of pyrimidines with a strong preference for U ↓ A and C ↓ A sequences (endoY ↓ A). EndoU-dependent cleavage was detected in every region of MHV RNA, from the 5' NTR to the 3' NTR, including transcriptional regulatory sequences (TRS). Cleavage at two CA dinucleotides immediately adjacent to the MHV poly(A) tail suggests a mechanism to suppress negative-strand RNA synthesis and the accumulation of viral dsRNA. MHV with EndoU (EndoUmut) or 2'-5' phosphodiesterase (PDEmut) mutations provoked the activation of RNase L in BMM, with corresponding cleavage of RNAs by RNase L. The physiologic targets of EndoU are viral RNA templates required for negative-strand RNA synthesis and dsRNA accumulation. Coronavirus EndoU cleaves U ↓ A and C ↓ A sequences (endoY ↓ A) within viral (+) strand RNA to evade dsRNA-activated host responses.

Single-Cell RNA Sequencing of Childhood Ependymoma Reveals Neoplastic Cell Subpopulations That Impact Molecular Classification and Etiology.

Ependymoma (EPN) is a brain tumor commonly presenting in childhood that remains fatal in most children. Intra-tumoral cellular heterogeneity in bulk-tumor samples significantly confounds our understanding of EPN biology, impeding development of effective therapy. We, therefore, use single-cell RNA sequencing, histology, and deconvolution to catalog cellular heterogeneity of the major childhood EPN subgroups. Analysis of PFA subgroup EPN reveals evidence of an undifferentiated progenitor subpopulation that either differentiates into subpopulations with ependymal cell characteristics or transitions into a mesenchymal subpopulation. Histological analysis reveals that progenitor and mesenchymal subpopulations co-localize in peri-necrotic zones. In conflict with current classification paradigms, relative PFA subpopulation proportions are shown to determine bulk-tumor-assigned subgroups. We provide an interactive online resource that facilitates exploration of the EPN single-cell dataset. This atlas of EPN cellular heterogeneity increases understanding of EPN biology.

clustifyr: an R package for automated single-cell RNA sequencing cluster classification.

Assignment of cell types from single-cell RNA sequencing (scRNA-seq) data remains a time-consuming and error-prone process. Current packages for identity assignment use limited types of reference data and often have rigid data structure requirements. We developed the clustifyr R package to leverage several external data types, including gene expression profiles to assign likely cell types using data from scRNA-seq, bulk RNA-seq, microarray expression data, or signature gene lists. We benchmark various parameters of a correlation-based approach and implement gene list enrichment methods. clustifyr is a lightweight and effective cell-type assignment tool developed for compatibility with various scRNA-seq analysis workflows. clustifyr is publicly available at https://github.com/rnabioco/clustifyr.

Simultaneous measurement of biochemical phenotypes and gene expression in single cells.

Methods to measure heterogeneity among cells are rapidly transforming our understanding of biology but are currently limited to molecular abundance measurements. We developed an approach to simultaneously measure biochemical activities and mRNA abundance in single cells to understand the heterogeneity of DNA repair activities across thousands of human lymphocytes, identifying known and novel cell-type-specific DNA repair phenotypes.

Dynamic RNA Regulation in the Brain Underlies Physiological Plasticity in a Hibernating Mammal.

Hibernation is a physiological and behavioral phenotype that minimizes energy expenditure. Hibernators cycle between profound depression and rapid hyperactivation of multiple physiological processes, challenging our concept of mammalian homeostasis. How the hibernator orchestrates and survives these extremes while maintaining cell to organismal viability is unknown. Here, we enhance the genome integrity and annotation of a model hibernator, the 13-lined ground squirrel. Our new assembly brings this genome to near chromosome-level contiguity and adds thousands of previously unannotated genes. These new genomic resources were used to identify 6,505 hibernation-related, differentially-expressed and processed transcripts using RNA-seq data from three brain regions in animals whose physiological status was precisely defined using body temperature telemetry. A software tool, squirrelBox, was developed to foster further data analyses and visualization. SquirrelBox includes a comprehensive toolset for rapid visualization of gene level and cluster group dynamics, sequence scanning of k-mer and domains, and interactive exploration of gene lists. Using these new tools and data, we deconvolute seasonal from temperature-dependent effects on the brain transcriptome during hibernation for the first time, highlighting the importance of carefully timed samples for studies of differential gene expression in hibernation. The identified genes include a regulatory network of RNA binding proteins that are dynamic in hibernation along with the composition of the RNA pool. In addition to passive effects of temperature, we provide evidence for regulated transcription and RNA turnover during hibernation. Significant alternative splicing, largely temperature dependent, also occurs during hibernation. These findings form a crucial first step and provide a roadmap for future work toward defining novel mechanisms of tissue protection and metabolic depression that may 1 day be applied toward improving human health.

Chronic Liver Disease in Humans Causes Expansion and Differentiation of Liver Lymphatic Endothelial Cells.

Liver lymphatic vessels support liver function by draining interstitial fluid, cholesterol, fat, and immune cells for surveillance in the liver draining lymph node. Chronic liver disease is associated with increased inflammation and immune cell infiltrate. However, it is currently unknown if or how lymphatic vessels respond to increased inflammation and immune cell infiltrate in the liver during chronic disease. Here we demonstrate that lymphatic vessel abundance increases in patients with chronic liver disease and is associated with areas of fibrosis and immune cell infiltration. Using single-cell mRNA sequencing and multi-spectral immunofluorescence analysis we identified liver lymphatic endothelial cells and found that chronic liver disease results in lymphatic endothelial cells (LECs) that are in active cell cycle with increased expression of CCL21. Additionally, we found that LECs from patients with NASH adopt a transcriptional program associated with increased IL13 signaling. Moreover, we found that oxidized low density lipoprotein, associated with NASH pathogenesis, induced the transcription and protein production of IL13 in LECs both in vitro and in a mouse model. Finally, we show that oxidized low density lipoprotein reduced the transcription of PROX1 and decreased lymphatic stability. Together these data indicate that LECs are active participants in the liver, expanding in an attempt to maintain tissue homeostasis. However, when inflammatory signals, such as oxidized low density lipoprotein are increased, as in NASH, lymphatic function declines and liver homeostasis is impeded.

Multiple decay events target HAC1 mRNA during splicing to regulate the unfolded protein response.

In the unfolded protein response (UPR), stress in the endoplasmic reticulum (ER) activates a large transcriptional program to increase ER folding capacity. During the budding yeast UPR, Ire1 excises an intron from the HAC1 mRNA and the exon products of cleavage are ligated, and the translated protein induces hundreds of stress-response genes. Using cells with mutations in RNA repair and decay enzymes, we show that phosphorylation of two different HAC1 splicing intermediates is required for their degradation by the 5'→3' exonuclease Xrn1 to enact opposing effects on the UPR. We also found that ligated but 2'-phosphorylated HAC1 mRNA is cleaved, yielding a decay intermediate with both 5'- and 2'-phosphates at its 5'-end that inhibit 5'→3' decay and suggesting that Ire1 degrades incompletely processed HAC1. These decay events expand the scope of RNA-based regulation in the budding yeast UPR and have implications for the control of the metazoan UPR.