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.

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.

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.

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.

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.

Recovery and analysis of transcriptome subsets from pooled single-cell RNA-seq libraries.

Single-cell RNA sequencing (scRNA-seq) methods generate sparse gene expression profiles for thousands of single cells in a single experiment. The information in these profiles is sufficient to classify cell types by distinct expression patterns but the high complexity of scRNA-seq libraries often prevents full characterization of transcriptomes from individual cells. To extract more focused gene expression information from scRNA-seq libraries, we developed a strategy to physically recover the DNA molecules comprising transcriptome subsets, enabling deeper interrogation of the isolated molecules by another round of DNA sequencing. We applied the method in cell-centric and gene-centric modes to isolate cDNA fragments from scRNA-seq libraries. First, we resampled the transcriptomes of rare, single megakaryocytes from a complex mixture of lymphocytes and analyzed them in a second round of DNA sequencing, yielding up to 20-fold greater sequencing depth per cell and increasing the number of genes detected per cell from a median of 1313 to 2002. We similarly isolated mRNAs from targeted T cells to improve the reconstruction of their VDJ-rearranged immune receptor mRNAs. Second, we isolated CD3D mRNA fragments expressed across cells in a scRNA-seq library prepared from a clonal T cell line, increasing the number of cells with detected CD3D expression from 59.7% to 100%. Transcriptome resampling is a general approach to recover targeted gene expression information from single-cell RNA sequencing libraries that enhances the utility of these costly experiments, and may be applicable to the targeted recovery of molecules from other single-cell assays.

Genetic bypass of essential RNA repair enzymes in budding yeast.

RNA repair enzymes catalyze rejoining of an RNA molecule after cleavage of phosphodiester linkages. RNA repair in budding yeast is catalyzed by two separate enzymes that process tRNA exons during their splicing and HAC1 mRNA exons during activation of the unfolded protein response (UPR). The RNA ligase Trl1 joins 2',3'-cyclic phosphate and 5'-hydroxyl RNA fragments, creating a phosphodiester linkage with a 2'-phosphate at the junction. The 2'-phosphate is removed by the 2'-phosphotransferase Tpt1. We bypassed the essential functions of TRL1 and TPT1 in budding yeast by expressing "prespliced," intronless versions of the 10 normally intron-containing tRNAs, indicating this repair pathway does not have additional essential functions. Consistent with previous studies, expression of intronless tRNAs failed to rescue the growth of cells with deletions in components of the SEN complex, implying an additional essential role for the splicing endonuclease. The trl1Δ and tpt1Δ mutants accumulate tRNA and HAC1 splicing intermediates indicative of RNA repair defects and are hypersensitive to drugs that inhibit translation. Failure to induce the unfolded protein response in trl1Δ cells grown with tunicamycin is lethal owing to their inability to ligate HAC1 after its cleavage by Ire1. In contrast, tpt1Δ mutants grow in the presence of tunicamycin despite reduced accumulation of spliced HAC1 mRNA. We optimized a PCR-based method to detect RNA 2'-phosphate modifications and show they are present on ligated HAC1 mRNA. These RNA repair mutants enable new studies of the role of RNA repair in cellular physiology.

Dynamic temperature-sensitive A-to-I RNA editing in the brain of a heterothermic mammal during hibernation.

RNA editing diversifies genomically encoded information to expand the complexity of the transcriptome. In ectothermic organisms, including Drosophila and Cephalopoda, where body temperature mirrors ambient temperature, decreases in environmental temperature lead to increases in A-to-I RNA editing and cause amino acid recoding events that are thought to be adaptive responses to temperature fluctuations. In contrast, endothermic mammals, including humans and mice, typically maintain a constant body temperature despite environmental changes. Here, A-to-I editing primarily targets repeat elements, rarely results in the recoding of amino acids, and plays a critical role in innate immune tolerance. Hibernating ground squirrels provide a unique opportunity to examine RNA editing in a heterothermic mammal whose body temperature varies over 30°C and can be maintained at 5°C for many days during torpor. We profiled the transcriptome in three brain regions at six physiological states to quantify RNA editing and determine whether cold-induced RNA editing modifies the transcriptome as a potential mechanism for neuroprotection at low temperature during hibernation. We identified 5165 A-to-I editing sites in 1205 genes with dynamically increased editing after prolonged cold exposure. The majority (99.6%) of the cold-increased editing sites are outside of previously annotated coding regions, 82.7% lie in SINE-derived repeats, and 12 sites are predicted to recode amino acids. Additionally, A-to-I editing frequencies increase with increasing cold-exposure, demonstrating that ADAR remains active during torpor. Our findings suggest that dynamic A-to-I editing at low body temperature may provide a neuroprotective mechanism to limit aberrant dsRNA accumulation during torpor in the mammalian hibernator.

O Cdc7 kinase where art thou?

Although Cdc7 protein kinase is important for regulating DNA replication in all eukaryotes and is a target for cancer therapy, it has never been localized in cells. Recently, a novel molecular genomic method used by our laboratory to localize Cdc7 to regions of chromosomes. Originally, mutations in the CDC7 gene were found in the classic cdc mutant collection of Hartwell et al. (Genetics 74:267-286, 1973). The CDC7 gene was found to encode a protein kinase called DDK that has been studied for many years, establishing its precise role in the initiation of DNA replication at origins. Recently, clinical studies are underway with DDK inhibitors against DDK in cancer patients. However, the conundrum is that Cdc7 has never been detected at origins of replication even though many studies have suggested it should be there. We used "Calling Card" system in which DNA binding proteins are localized to the genome via retrotransposon insertion and deep-sequencing methods. We have shown that Cdc7 localizes at many regions of the genome and was enriched at functional origins of replication. These results are consistent with DDK's role in many additional genomic processes including mutagenesis, chromatid cohesion, and meiotic recombination. Thus, the main conclusion from our studies is that Cdc7 kinase is found at many locations in the genome, but is enriched at functional origins of replication. Furthermore, we propose that application of the Calling Card system to other eukaryotes should be useful in identification of functional origins in other eukaryotic cells.

Ribose-seq: global mapping of ribonucleotides embedded in genomic DNA.

Abundant ribonucleotide incorporation in DNA during replication and repair has profound consequences for genome stability, but the global distribution of ribonucleotide incorporation is unknown. We developed ribose-seq, a method for capturing unique products generated by alkaline cleavage of DNA at embedded ribonucleotides. High-throughput sequencing of these fragments in DNA from the yeast Saccharomyces cerevisiae revealed widespread ribonucleotide distribution, with a strong preference for cytidine and guanosine, and identified hotspots of ribonucleotide incorporation in nuclear and mitochondrial DNA. Ribonucleotides were primarily incorporated on the newly synthesized leading strand of nuclear DNA and were present upstream of (G+C)-rich tracts in the mitochondrial genome. Ribose-seq is a powerful tool for the systematic profiling of ribonucleotide incorporation in genomic DNA.

High resolution mapping of modified DNA nucleobases using excision repair enzymes.

The incorporation and creation of modified nucleobases in DNA have profound effects on genome function. We describe methods for mapping positions and local content of modified DNA nucleobases in genomic DNA. We combined in vitro nucleobase excision with massively parallel DNA sequencing (Excision-seq) to determine the locations of modified nucleobases in genomic DNA. We applied the Excision-seq method to map uracil in E. coli and budding yeast and discovered significant variation in uracil content, wherein uracil is excluded from the earliest and latest replicating regions of the genome, possibly driven by changes in nucleotide pool composition. We also used Excision-seq to identify sites of pyrimidine dimer formation induced by UV light exposure, where the method could distinguish between sites of cyclobutane and 6-4 photoproduct formation. These UV mapping data enabled analysis of local sequence bias around pyrimidine dimers and suggested a preference for an adenosine downstream from 6-4 photoproducts. The Excision-seq method is broadly applicable for high precision, genome-wide mapping of modified nucleobases with cognate repair enzymes.

Global analysis of RNA cleavage by 5'-hydroxyl RNA sequencing.

RNA cleavage by some endoribonucleases and self-cleaving ribozymes produces RNA fragments with 5'-hydroxyl (5'-OH) and 2',3'-cyclic phosphate termini. To identify 5'-OH RNA fragments produced by these cleavage events, we exploited the unique ligation mechanism of Escherichia coli RtcB RNA ligase to attach an oligonucleotide linker to RNAs with 5'-OH termini, followed by steps for library construction and analysis by massively parallel DNA sequencing. We applied the method to RNA from budding yeast and captured known 5'-OH fragments produced by tRNA Splicing Endonuclease (SEN) during processing of intron-containing pre-tRNAs and by Ire1 cleavage of HAC1 mRNA following induction of the unfolded protein response (UPR). We identified numerous novel 5'-OH fragments derived from mRNAs: some 5'-OH mRNA fragments were derived from single, localized cleavages, while others were likely produced by multiple, distributed cleavages. Many 5'-OH fragments derived from mRNAs were produced upstream of codons for highly electrostatic peptides, suggesting that the fragments may be generated by co-translational mRNA decay. Several 5'-OH RNA fragments accumulated during the induction of the UPR, some of which share a common sequence motif that may direct cleavage of these mRNAs. This method enables specific capture of 5'-OH termini and complements existing methods for identifying RNAs with 2',3'-cyclic phosphate termini.

Coordination between Drosophila Arc1 and a specific population of brain neurons regulates organismal fat.

The brain plays a critical yet incompletely understood role in regulating organismal fat. We performed a neuronal silencing screen in Drosophila larvae to identify brain regions required to maintain proper levels of organismal fat. When used to modulate synaptic activity in specific brain regions, the enhancer-trap driver line E347 elevated fat upon neuronal silencing, and decreased fat upon neuronal activation. Unbiased sequencing revealed that Arc1 mRNA levels increase upon E347 activation. We had previously identified Arc1 mutations in a high-fat screen. Here we reveal metabolic changes in Arc1 mutants consistent with a high-fat phenotype and an overall shift toward energy storage. We find that Arc1-expressing cells neighbor E347 neurons, and manipulating E347 synaptic activity alters Arc1 expression patterns. Elevating Arc1 expression in these cells decreased fat, a phenocopy of E347 activation. Finally, loss of Arc1 prevented the lean phenotype caused by E347 activation, suggesting that Arc1 activity is required for E347 control of body fat. Importantly, neither E347 nor Arc1 manipulation altered energy-related behaviors. Our results support a model wherein E347 neurons induce Arc1 in specific neighboring cells to prevent excess fat accumulation.

Improvements to the HITS-CLIP protocol eliminate widespread mispriming artifacts.

BACKGROUND: High-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP) allows for high resolution, genome-wide mapping of RNA-binding proteins. This methodology is frequently used to validate predicted targets of microRNA binding, as well as direct targets of other RNA-binding proteins. Hence, the accuracy and sensitivity of binding site identification is critical. RESULTS: We found that substantial mispriming during reverse transcription results in the overrepresentation of sequences complementary to the primer used for reverse transcription. Up to 45 % of peaks in publicly available HITS-CLIP libraries are attributable to this mispriming artifact, and the majority of libraries have detectable levels of mispriming. We also found that standard techniques for validating microRNA-target interactions fail to differentiate between artifactual peaks and physiologically relevant peaks. CONCLUSIONS: Here, we present a modification to the HITS-CLIP protocol that effectively eliminates this artifact and improves the sensitivity and complexity of resulting libraries.

A homolog of lariat-debranching enzyme modulates turnover of branched RNA.

Turnover of the branched RNA intermediates and products of pre-mRNA splicing is mediated by the lariat-debranching enzyme Dbr1. We characterized a homolog of Dbr1 from Saccharomyces cerevisiae, Drn1/Ygr093w, that has a pseudo-metallophosphodiesterase domain with primary sequence homology to Dbr1 but lacks essential active site residues found in Dbr1. Whereas loss of Dbr1 results in lariat-introns failing broadly to turnover, loss of Drn1 causes low levels of lariat-intron accumulation. Conserved residues in the Drn1 C-terminal CwfJ domains, which are not present in Dbr1, are required for efficient intron turnover. Drn1 interacts with Dbr1, components of the Nineteen Complex, U2 snRNA, branched intermediates, and products of splicing. Drn1 enhances debranching catalyzed by Dbr1 in vitro, but does so without significantly improving the affinity of Dbr1 for branched RNA. Splicing carried out in in vitro extracts in the absence of Drn1 results in an accumulation of branched splicing intermediates and products released from the spliceosome, likely due to less active debranching, as well as the promiscuous release of cleaved 5'-exon. Drn1 enhances Dbr1-mediated turnover of lariat-intermediates and lariat-intron products, indicating that branched RNA turnover is regulated at multiple steps during splicing.

RNase L targets distinct sites in influenza A virus RNAs.

UNLABELLED: Influenza A virus (IAV) infections are influenced by type 1 interferon-mediated antiviral defenses and by viral countermeasures to these defenses. When IAV NS1 protein is disabled, RNase L restricts virus replication; however, the RNAs targeted for cleavage by RNase L under these conditions have not been defined. In this study, we used deep-sequencing methods to identify RNase L cleavage sites within host and viral RNAs from IAV PR8ΔNS1-infected A549 cells. Short hairpin RNA knockdown of RNase L allowed us to distinguish between RNase L-dependent and RNase L-independent cleavage sites. RNase L-dependent cleavage sites were evident at discrete locations in IAV RNA segments (both positive and negative strands). Cleavage in PB2, PB1, and PA genomic RNAs suggests that viral RNPs are susceptible to cleavage by RNase L. Prominent amounts of cleavage mapped to specific regions within IAV RNAs, including some areas of increased synonymous-site conservation. Among cellular RNAs, RNase L-dependent cleavage was most frequent at precise locations in rRNAs. Our data show that RNase L targets specific sites in both host and viral RNAs to restrict influenza virus replication when NS1 protein is disabled. IMPORTANCE: RNase L is a critical component of interferon-regulated and double-stranded-RNA-activated antiviral host responses. We sought to determine how RNase L exerts its antiviral activity during influenza virus infection. We enhanced the antiviral activity of RNase L by disabling a viral protein, NS1, that inhibits the activation of RNase L. Then, using deep-sequencing methods, we identified the host and viral RNAs targeted by RNase L. We found that RNase L cleaved viral RNAs and rRNAs at very precise locations. The direct cleavage of IAV RNAs by RNase L highlights an intimate battle between viral RNAs and an antiviral endonuclease.

The RtcB RNA ligase is an essential component of the metazoan unfolded protein response.

RNA ligation can regulate RNA function by altering RNA sequence, structure and coding potential. For example, the function of XBP1 in mediating the unfolded protein response requires RNA ligation, as does the maturation of some tRNAs. Here, we describe a novel in vivo model in Caenorhabditis elegans for the conserved RNA ligase RtcB and show that RtcB ligates the xbp-1 mRNA during the IRE-1 branch of the unfolded protein response. Without RtcB, protein stress results in the accumulation of unligated xbp-1 mRNA fragments, defects in the unfolded protein response, and decreased lifespan. RtcB also ligates endogenous pre-tRNA halves, and RtcB mutants have defects in growth and lifespan that can be bypassed by expression of pre-spliced tRNAs. In addition, animals that lack RtcB have defects that are independent of tRNA maturation and the unfolded protein response. Thus, RNA ligation by RtcB is required for the function of multiple endogenous target RNAs including both xbp-1 and tRNAs. RtcB is uniquely capable of performing these ligation functions, and RNA ligation by RtcB mediates multiple essential processes in vivo.

Ribonuclease L and metal-ion-independent endoribonuclease cleavage sites in host and viral RNAs.

Ribonuclease L (RNase L) is a metal-ion-independent endoribonuclease associated with antiviral and antibacterial defense, cancer and lifespan. Despite the biological significance of RNase L, the RNAs cleaved by this enzyme are poorly defined. In this study, we used deep sequencing methods to reveal the frequency and location of RNase L cleavage sites within host and viral RNAs. To make cDNA libraries, we exploited the 2', 3'-cyclic phosphate at the end of RNA fragments produced by RNase L and other metal-ion-independent endoribonucleases. We optimized and validated 2', 3'-cyclic phosphate cDNA synthesis and Illumina sequencing methods using viral RNAs cleaved with purified RNase L, viral RNAs cleaved with purified RNase A and RNA from uninfected and poliovirus-infected HeLa cells. Using these methods, we identified (i) discrete regions of hepatitis C virus and poliovirus RNA genomes that were profoundly susceptible to RNase L and other single-strand specific endoribonucleases, (ii) RNase L-dependent and RNase L-independent cleavage sites within ribosomal RNAs (rRNAs) and (iii) 2', 3'-cyclic phosphates at the ends of 5S rRNA and U6 snRNA. Monitoring the frequency and location of metal-ion-independent endoribonuclease cleavage sites within host and viral RNAs reveals, in part, how these enzymes contribute to health and disease.

Genome-wide analysis of miRNA-mRNA interactions in marrow stromal cells.

Regulation of hematopoietic stem cell proliferation, lineage commitment, and differentiation in adult vertebrates requires extrinsic signals provided by cells in the marrow microenvironment (ME) located within the bone marrow. Both secreted and cell-surface bound factors critical to this regulation have been identified, yet control of their expression by cells within the ME has not been addressed. Herein we hypothesize that microRNAs (miRNAs) contribute to their controlled expression. MiRNAs are small noncoding RNAs that bind to target mRNAs and downregulate gene expression by either initiating mRNA degradation or preventing peptide translation. Testing the role of miRNAs in downregulating gene expression has been difficult since conventional techniques used to define miRNA-mRNA interactions are indirect and have high false-positive and negative rates. In this report, a genome-wide biochemical technique (high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation or HITS-CLIP) was used to generate unbiased genome-wide maps of miRNA-mRNA interactions in two critical cellular components of the marrow ME: marrow stromal cells and bone marrow endothelial cells. Analysis of these datasets identified miRNAs as direct regulators of JAG1, WNT5A, MMP2, and VEGFA; four factors that are important to ME function. Our results show the feasibility and utility of unbiased genome-wide biochemical techniques in dissecting the role of miRNAs in regulation of complex tissues such as the marrow ME.