Improving nanopore read accuracy with the R2C2 method enables the sequencing of highly multiplexed full-length single-cell cDNA.

High-throughput short-read sequencing has revolutionized how transcriptomes are quantified and annotated. However, while Illumina short-read sequencers can be used to analyze entire transcriptomes down to the level of individual splicing events with great accuracy, they fall short of analyzing how these individual events are combined into complete RNA transcript isoforms. Because of this shortfall, long-distance information is required to complement short-read sequencing to analyze transcriptomes on the level of full-length RNA transcript isoforms. While long-read sequencing technology can provide this long-distance information, there are issues with both Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT) long-read sequencing technologies that prevent their widespread adoption. Briefly, PacBio sequencers produce low numbers of reads with high accuracy, while ONT sequencers produce higher numbers of reads with lower accuracy. Here, we introduce and validate a long-read ONT-based sequencing method. At the same cost, our Rolling Circle Amplification to Concatemeric Consensus (R2C2) method generates more accurate reads of full-length RNA transcript isoforms than any other available long-read sequencing method. These reads can then be used to generate isoform-level transcriptomes for both genome annotation and differential expression analysis in bulk or single-cell samples

Depletion of Hemoglobin Transcripts and Long-Read Sequencing Improves the Transcriptome Annotation of the Polar Bear (Ursus maritimus)

Transcriptome studies evaluating whole blood and tissues are often confounded by overrepresentation of highly abundant transcripts. These abundant transcripts are problematic, as they compete with and prevent the detection of rare RNA transcripts, obscuring their biological importance. This issue is more pronounced when using long-read sequencing technologies for isoform-level transcriptome analysis, as they have relatively lower throughput compared to short-read sequencers. As a result, long-read based transcriptome analysis is prohibitively expensive for non-model organisms. While there are off-the-shelf kits available for select model organisms capable of depleting highly abundant transcripts for alpha (HBA) and beta (HBB) hemoglobin, they are unsuitable for non-model organisms. To address this, we have adapted the recent CRISPR/Cas9-based depletion method (depletion of abundant sequences by hybridization) for long-read full-length cDNA sequencing approaches that we call Long-DASH. Using a recombinant Cas9 protein with appropriate guide RNAs, full-length hemoglobin transcripts can be depleted in vitro prior to performing any short- and long-read sequencing library preparations. Using this method, we sequenced depleted full-length cDNA in parallel using both our Oxford Nanopore Technology (ONT) based R2C2 long-read approach, as well as the Illumina short-read based Smart-seq2 approach. To showcase this, we have applied our methods to create an isoform-level transcriptome from whole blood samples derived from three polar bears (Ursus maritimus). Using Long-DASH, we succeeded in depleting hemoglobin transcripts and generated deep Smart-seq2 Illumina datasets and 3.8 million R2C2 full-length cDNA consensus reads. Applying Long-DASH with our isoform identification pipeline, Mandalorion, we discovered ∼6,000 high-confidence isoforms and a number of novel genes. This indicates that there is a high diversity of gene isoforms within U. maritimus not yet reported. This reproducible and straightforward approach has not only improved the polar bear transcriptome annotations but will serve as the foundation for future efforts to investigate transcriptional dynamics within the 19 polar bear subpopulations around the Arctic

Circadian Oscillations of Protein-Coding and Regulatory RNAs in a Highly Dynamic Mammalian Liver Epigenome

SummaryIn the mouse liver, circadian transcriptional rhythms are necessary for metabolic homeostasis. Whether dynamic epigenomic modifications are associated with transcript oscillations has not been systematically investigated. We found that several antisense RNA, lincRNA, and microRNA transcripts also showed circadian oscillations in adult mouse livers. Robust transcript oscillations often correlated with rhythmic histone modifications in promoters, gene bodies, or enhancers, although promoter DNA methylation levels were relatively stable. Such integrative analyses identified oscillating expression of an antisense transcript (asPer2) to the gene encoding the circadian oscillator component Per2. Robust transcript oscillations often accompanied rhythms in multiple histone modifications and recruitment of multiple chromatin-associated clock components. Coupling of cycling histone modifications with nearby oscillating transcripts thus established a temporal relationship between enhancers, genes, and transcripts on a genome-wide scale in a mammalian liver. The results offer a framework for understanding the dynamics of metabolism, circadian clock, and chromatin modifications involved in metabolic homeostasis

Horizontal Transmission and Recombination Maintain forever Young Bacterial Symbiont Genomes

Bacterial symbionts bring a wealth of functions to the associations they participate in, but by doing so, they endanger the genes and genomes underlying these abilities. When bacterial symbionts become obligately associated with their hosts, their genomes are thought to decay towards an organelle-like fate due to decreased homologous recombination and inef- ficient selection. However, numerous associations exist that counter these expectations, especially in marine environments, possibly due to ongoing horizontal gene flow. Despite extensive theoretical treatment, no empirical study thus far has connected these underlying population genetic processes with long-term evolutionary outcomes. By sampling marine chemosynthetic bacterial-bivalve endosymbioses that range from primarily vertical to strictly horizontal transmission, we tested this canonical theory. We found that transmission mode strongly predicts homologous recombination rates, and that exceedingly low recombination rates are associated with moderate genome degradation in the marine symbionts with nearly strict vertical transmission. Nonetheless, even the most degraded marine endosym- biont genomes are occasionally horizontally transmitted and are much larger than their ter- restrial insect symbiont counterparts. Therefore, horizontal transmission and recombination enable efficient natural selection to maintain intermediate symbiont genome sizes and sub- stantial functional genetic variation

Harmonics of Circadian Gene Transcription in Mammals

The circadian clock is a molecular and cellular oscillator found in most mammalian tissues that regulates rhythmic physiology and behavior. Numerous investigations have addressed the contribution of circadian rhythmicity to cellular, organ, and organismal physiology. We recently developed a method to look at transcriptional oscillations with unprecedented precision and accuracy using high-density time sampling. Here, we report a comparison of oscillating transcription from mouse liver, NIH3T3, and U2OS cells. Several surprising observations resulted from this study, including a 100-fold difference in the number of cycling transcripts in autonomous cellular models of the oscillator versus tissues harvested from intact mice. Strikingly, we found two clusters of genes that cycle at the second and third harmonic of circadian rhythmicity in liver, but not cultured cells. Validation experiments show that 12-hour oscillatory transcripts occur in several other peripheral tissues as well including heart, kidney, and lungs. These harmonics are lost ex vivo, as well as under restricted feeding conditions. Taken in sum, these studies illustrate the importance of time sampling with respect to multiple testing, suggest caution in use of autonomous cellular models to study clock output, and demonstrate the existence of harmonics of circadian gene expression in the mouse

Inducible Ablation of Melanopsin-Expressing Retinal Ganglion Cells Reveals Their Central Role in Non-Image Forming Visual Responses

Rod/cone photoreceptors of the outer retina and the melanopsin-expressing retinal ganglion cells (mRGCs) of the inner retina mediate non-image forming visual responses including entrainment of the circadian clock to the ambient light, the pupillary light reflex (PLR), and light modulation of activity. Targeted deletion of the melanopsin gene attenuates these adaptive responses with no apparent change in the development and morphology of the mRGCs. Comprehensive identification of mRGCs and knowledge of their specific roles in image-forming and non-image forming photoresponses are currently lacking. We used a Cre-dependent GFP expression strategy in mice to genetically label the mRGCs. This revealed that only a subset of mRGCs express enough immunocytochemically detectable levels of melanopsin. We also used a Cre-inducible diphtheria toxin receptor (iDTR) expression approach to express the DTR in mRGCs. mRGCs develop normally, but can be acutely ablated upon diphtheria toxin administration. The mRGC-ablated mice exhibited normal outer retinal function. However, they completely lacked non-image forming visual responses such as circadian photoentrainment, light modulation of activity, and PLR. These results point to the mRGCs as the site of functional integration of the rod/cone and melanopsin phototransduction pathways and as the primary anatomical site for the divergence of image-forming and non-image forming photoresponses in mammals

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification

The Long-read RNA-Seq Genome Annotation Assessment Project (LRGASP) Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. The consortium generated over 427 million long-read sequences from cDNA and direct RNA datasets, encompassing human, mouse, and manatee species, using different protocols and sequencing platforms. These data were utilized by developers to address challenges in transcript isoform detection and quantification, as well as de novo transcript isoform identification. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. When aiming to detect rare and novel transcripts or when using reference-free approaches, incorporating additional orthogonal data and replicate samples are advised. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis

Acute Liver Injury Is Independent of B Cells or Immunoglobulin M

Acute liver injury is a clinically important pathology and results in the release of Danger Associated Molecular Patterns, which initiate an immune response. Withdrawal of the injurious agent and curtailing any pathogenic secondary immune response may allow spontaneous resolution of injury. The role B cells and Immunoglobulin M (IgM) play in acute liver injury is largely unknown and it was proposed that B cells and/or IgM would play a significant role in its pathogenesis.Tissue from 3 models of experimental liver injury (ischemia-reperfusion injury, concanavalin A hepatitis and paracetamol-induced liver injury) and patients transplanted following paracetamol overdose were stained for evidence of IgM deposition. Mice deficient in B cells (and IgM) were used to dissect out the role B cells and/or IgM played in the development or resolution of injury. Serum transfer into mice lacking IgM was used to establish the role IgM plays in injury.Significant deposition of IgM was seen in the explanted livers of patients transplanted following paracetamol overdose as well as in 3 experimental models of acute liver injury (ischemia-reperfusion injury, concanavalin A hepatitis and paracetamol-induced liver injury). Serum transfer into IgM-deficient mice failed to reconstitute injury (p = 0.66), despite successful engraftment of IgM. Mice deficient in both T and B cells (RAG1-/-) mice (p<0.001), but not B cell deficient (μMT) mice (p = 0.93), were significantly protected from injury. Further interrogation with T cell deficient (CD3εKO) mice confirmed that the T cell component is a key mediator of sterile liver injury. Mice deficient in B cells and IgM mice did not have a significant delay in resolution following acute liver injury.IgM deposition appears to be common feature of both human and murine sterile liver injury. However, neither IgM nor B cells, play a significant role in the development of or resolution from acute liver injury. T cells appear to be key mediators of injury. In conclusion, the therapeutic targeting of IgM or B cells (e.g. with Rituximab) would have limited benefit in protecting patients from acute liver injury

Systematic assessment of long-read RNA-seq methods for transcript identification and quantification

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis

Single-cell isoform analysis in human immune cells

High-throughput single-cell analysis today is facilitated by protocols like the 10X Genomics platform or Drop-Seq which generate cDNA pools in which the origin of a transcript is encoded at its 5' or 3' end. Here, we used R2C2 to sequence and demultiplex 12 million full-length cDNA molecules generated by the 10X Genomics platform from ~3000 peripheral blood mononuclear cells. We use these reads, independent from Illumina data, to identify B cell, T cell, and monocyte clusters and generate isoform-level transcriptomes for cells and cell types. Finally, we extract paired adaptive immune receptor sequences unique to each T and B cell