All coronaviruses known to have recently emerged as human pathogens probably originated in bats
. Here we use a single experimental platform based on immunodeficient mice implanted with human lung ...tissue (hereafter, human lung-only mice (LoM)) to demonstrate the efficient in vivo replication of severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East respiratory syndrome coronavirus (MERS-CoV) and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), as well as two endogenous SARS-like bat coronaviruses that show potential for emergence as human pathogens. Virus replication in this model occurs in bona fide human lung tissue and does not require any type of adaptation of the virus or the host. Our results indicate that bats contain endogenous coronaviruses that are capable of direct transmission to humans. Our detailed analysis of in vivo infection with SARS-CoV-2 in human lung tissue from LoM showed a predominant infection of human lung epithelial cells, including type-2 pneumocytes that are present in alveoli and ciliated airway cells. Acute infection with SARS-CoV-2 was highly cytopathic and induced a robust and sustained type-I interferon and inflammatory cytokine and chemokine response. Finally, we evaluated a therapeutic and pre-exposure prophylaxis strategy for SARS-CoV-2 infection. Our results show that therapeutic and prophylactic administration of EIDD-2801-an oral broad-spectrum antiviral agent that is currently in phase II/III clinical trials-markedly inhibited SARS-CoV-2 replication in vivo, and thus has considerable potential for the prevention and treatment of COVID-19.
Alignment is the first step in most RNA-seq analysis pipelines, and the accuracy of downstream analyses depends heavily on it. Unlike most steps in the pipeline, alignment is particularly amenable to ...benchmarking with simulated data. We performed a comprehensive benchmarking of 14 common splice-aware aligners for base, read, and exon junction-level accuracy and compared default with optimized parameters. We found that performance varied by genome complexity, and accuracy and popularity were poorly correlated. The most widely cited tool underperforms for most metrics, particularly when using default settings.
Although the genome is generally thought to be transcriptionally silent during mitosis, technical limitations have prevented sensitive mapping of transcription during mitosis and mitotic exit. Thus, ...the means by which the interphase expression pattern is transduced to daughter cells have been unclear. We used 5-ethynyluridine to pulse-label transcripts during mitosis and mitotic exit and found that many genes exhibit transcription during mitosis, as confirmed with fluorescein isothiocyanate–uridine 5′-triphosphate labeling, RNA fluorescence in situ hybridization, and quantitative reverse transcription polymerase chain reaction. The first round of transcription immediately after mitosis primarily activates genes involved in the growth and rebuilding of daughter cells, rather than cell type–specific functions. We propose that the cell’s transcription pattern is largely retained at a low level through mitosis, whereas the amplitude of transcription observed in interphase is reestablished during mitotic exit.
The mechanisms underlying the formation of acyl protein modifications remain poorly understood. By investigating the reactivity of endogenous acyl-CoA metabolites, we found a class of acyl-CoAs that ...undergo intramolecular catalysis to form reactive intermediates that non-enzymatically modify proteins. Based on this mechanism, we predicted, validated, and characterized a protein modification: 3-hydroxy-3-methylglutaryl(HMG)-lysine. In a model of altered HMG-CoA metabolism, we found evidence of two additional protein modifications: 3-methylglutaconyl(MGc)-lysine and 3-methylglutaryl(MG)-lysine. Using quantitative proteomics, we compared the "acylomes" of two reactive acyl-CoA species, namely HMG-CoA and glutaryl-CoA, which are generated in different pathways. We found proteins that are uniquely modified by each reactive metabolite, as well as common proteins and pathways. We identified the tricarboxylic acid cycle as a pathway commonly regulated by acylation and validated malate dehydrogenase as a key target. These data uncover a fundamental relationship between reactive acyl-CoA species and proteins and define a new regulatory paradigm in metabolism.
High-throughput RNA sequencing is an increasingly accessible method for studying gene structure and activity on a genome-wide scale. A critical step in RNA-seq data analysis is the alignment of ...partial transcript reads to a reference genome sequence. To assess the performance of current mapping software, we invited developers of RNA-seq aligners to process four large human and mouse RNA-seq data sets. In total, we compared 26 mapping protocols based on 11 programs and pipelines and found major performance differences between methods on numerous benchmarks, including alignment yield, basewise accuracy, mismatch and gap placement, exon junction discovery and suitability of alignments for transcript reconstruction. We observed concordant results on real and simulated RNA-seq data, confirming the relevance of the metrics employed. Future developments in RNA-seq alignment methods would benefit from improved placement of multimapped reads, balanced utilization of existing gene annotation and a reduced false discovery rate for splice junctions.
Influenza is a leading cause of respiratory mortality and morbidity. While inflammation is essential for fighting infection, a balance of anti-viral defense and host tolerance is necessary for ...recovery. Circadian rhythms have been shown to modulate inflammation. However, the importance of diurnal variability in the timing of influenza infection is not well understood. Here we demonstrate that endogenous rhythms affect survival in influenza infection. Circadian control of influenza infection is mediated by enhanced inflammation as proven by increased cellularity in bronchoalveolar lavage (BAL), pulmonary transcriptomic profile and histology and is not attributable to viral burden. Better survival is associated with a time dependent preponderance of NK and NKT cells and lower proportion of inflammatory monocytes in the lung. Further, using a series of genetic mouse mutants, we elucidate cellular mechanisms underlying circadian gating of influenza infection.
Full-length isoform quantification from RNA-Seq is a key goal in transcriptomics analyses and has been an area of active development since the beginning. The fundamental difficulty stems from the ...fact that RNA transcripts are long, while RNA-Seq reads are short. Here we use simulated benchmarking data that reflects many properties of real data, including polymorphisms, intron signal and non-uniform coverage, allowing for systematic comparative analyses of isoform quantification accuracy and its impact on differential expression analysis. Genome, transcriptome and pseudo alignment-based methods are included; and a simple approach is included as a baseline control. Salmon, kallisto, RSEM, and Cufflinks exhibit the highest accuracy on idealized data, while on more realistic data they do not perform dramatically better than the simple approach. We determine the structural parameters with the greatest impact on quantification accuracy to be length and sequence compression complexity and not so much the number of isoforms. The effect of incomplete annotation on performance is also investigated. Overall, the tested methods show sufficient divergence from the truth to suggest that full-length isoform quantification and isoform level DE should still be employed selectively.
Defective viral genomes of the copy-back type (cbDVGs) are the primary initiators of the antiviral immune response during infection with respiratory syncytial virus (RSV) both in vitro and in vivo. ...However, the mechanism governing cbDVG generation remains unknown, thereby limiting our ability to manipulate cbDVG content in order to modulate the host response to infection. Here we report a specific genomic signal that mediates the generation of a subset of RSV cbDVG species. Using a customized bioinformatics tool, we identified regions in the RSV genome frequently used to generate cbDVGs during infection. We then created a minigenome system to validate the function of one of these sequences and to determine if specific nucleotides were essential for cbDVG generation at that position. Further, we created a recombinant virus unable to produce a subset of cbDVGs due to mutations introduced in this sequence. The identified sequence was also found as a site for cbDVG generation during natural RSV infections, and common cbDVGs originated at this sequence were found among samples from various infected patients. These data demonstrate that sequences encoded in the viral genome determine the location of cbDVG formation and, therefore, the generation of cbDVGs is not a stochastic process. These findings open the possibility of genetically manipulating cbDVG formation to modulate infection outcome.
Eukaryotic circadian clocks include transcriptional/translational feedback loops that drive 24-h rhythms of transcription. These transcriptional rhythms underlie oscillations of protein abundance, ...thereby mediating circadian rhythms of behavior, physiology, and metabolism. Numerous studies over the last decade have used microarrays to profile circadian transcriptional rhythms in various organisms and tissues. Here we use RNA sequencing (RNA-seq) to profile the circadian transcriptome of Drosophila melanogaster brain from wild-type and period-null clock-defective animals. We identify several hundred transcripts whose abundance oscillates with 24-h periods in either constant darkness or 12 h light/dark diurnal cycles, including several noncoding RNAs (ncRNAs) that were not identified in previous microarray studies. Of particular interest are U snoRNA host genes (Uhgs), a family of diurnal cycling noncoding RNAs that encode the precursors of more than 50 box-C/D small nucleolar RNAs, key regulators of ribosomal biogenesis. Transcriptional profiling at the level of individual exons reveals alternative splice isoforms for many genes whose relative abundances are regulated by either period or circadian time, although the effect of circadian time is muted in comparison to that of period. Interestingly, period loss of function significantly alters the frequency of RNA editing at several editing sites, suggesting an unexpected link between a key circadian gene and RNA editing. We also identify tens of thousands of novel splicing events beyond those previously annotated by the modENCODE Consortium, including several that affect key circadian genes. These studies demonstrate extensive circadian control of ncRNA expression, reveal the extent of clock control of alternative splicing and RNA editing, and provide a novel, genome-wide map of splicing in Drosophila brain.
A critical task in high-throughput sequencing is aligning millions of short reads to a reference genome. Alignment is especially complicated for RNA sequencing (RNA-Seq) because of RNA splicing. A ...number of RNA-Seq algorithms are available, and claim to align reads with high accuracy and efficiency while detecting splice junctions. RNA-Seq data are discrete in nature; therefore, with reasonable gene models and comparative metrics RNA-Seq data can be simulated to sufficient accuracy to enable meaningful benchmarking of alignment algorithms. The exercise to rigorously compare all viable published RNA-Seq algorithms has not been performed previously.
We developed an RNA-Seq simulator that models the main impediments to RNA alignment, including alternative splicing, insertions, deletions, substitutions, sequencing errors and intron signal. We used this simulator to measure the accuracy and robustness of available algorithms at the base and junction levels. Additionally, we used reverse transcription-polymerase chain reaction (RT-PCR) and Sanger sequencing to validate the ability of the algorithms to detect novel transcript features such as novel exons and alternative splicing in RNA-Seq data from mouse retina. A pipeline based on BLAT was developed to explore the performance of established tools for this problem, and to compare it to the recently developed methods. This pipeline, the RNA-Seq Unified Mapper (RUM), performs comparably to the best current aligners and provides an advantageous combination of accuracy, speed and usability.
The RUM pipeline is distributed via the Amazon Cloud and for computing clusters using the Sun Grid Engine (http://cbil.upenn.edu/RUM).
ggrant@pcbi.upenn.edu; epierce@mail.med.upenn.edu
The RNA-Seq sequence reads described in the article are deposited at GEO, accession GSE26248.