Sites of transcription initiation drive mRNA isoform selection.

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

The key features of SARS-CoV-2 leader and NSP1 required for viral escape of NSP1-mediated repression.

SARS-CoV-2, responsible for the ongoing global pandemic, must overcome a conundrum faced by all viruses. To achieve its own replication and spread, it simultaneously depends on and subverts cellular mechanisms. At the early stage of infection, SARS-CoV-2 expresses the viral nonstructural protein 1 (NSP1), which inhibits host translation by blocking the mRNA entry tunnel on the ribosome; this interferes with the binding of cellular mRNAs to the ribosome. Viral mRNAs, on the other hand, overcome this blockade. We show that NSP1 enhances expression of mRNAs containing the SARS-CoV-2 leader. The first stem-loop (SL1) in the viral leader is both necessary and sufficient for this enhancement mechanism. Our analysis pinpoints specific residues within SL1 (three cytosine residues at the positions 15, 19, and 20) and another within NSP1 (R124), which are required for viral evasion, and thus might present promising drug targets. We target SL1 with the antisense oligo (ASO) to efficiently and specifically down-regulate SARS-CoV-2 mRNA. Additionally, we carried out analysis of a functional interactome of NSP1 using BioID and identified components of antiviral defense pathways. Our analysis therefore suggests a mechanism by which NSP1 inhibits the expression of host genes while enhancing that of viral RNA. This analysis helps reconcile conflicting reports in the literature regarding the mechanisms by which the virus avoids NSP1 silencing.

Transcriptomic profiling of SARS-CoV-2 infected human cell lines identifies HSP90 as target for COVID-19 therapy.

Detailed knowledge of the molecular biology of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is crucial for understanding of viral replication, host responses, and disease progression. Here, we report gene expression profiles of three SARS-CoV- and SARS-CoV-2-infected human cell lines. SARS-CoV-2 elicited an approximately two-fold higher stimulation of the innate immune response compared to SARS-CoV in the human epithelial cell line Calu-3, including induction of miRNA-155. Single-cell RNA sequencing of infected cells showed that genes induced by virus infections were broadly upregulated, whereas interferon beta/lambda genes, a pro-inflammatory cytokines such as IL-6, were expressed only in small subsets of infected cells. Temporal analysis suggested that transcriptional activities of interferon regulatory factors precede those of nuclear factor κB. Lastly, we identified heat shock protein 90 (HSP90) as a protein relevant for the infection. Inhibition of the HSP90 activity resulted in a reduction of viral replication and pro-inflammatory cytokine expression in primary human airway epithelial cells.

Lymphoma Driver Mutations in the Pathogenic Evolution of an Iconic Human Autoantibody.

Pathogenic autoantibodies arise in many autoimmune diseases, but it is not understood how the cells making them evade immune checkpoints. Here, single-cell multi-omics analysis demonstrates a shared mechanism with lymphoid malignancy in the formation of public rheumatoid factor autoantibodies responsible for mixed cryoglobulinemic vasculitis. By combining single-cell DNA and RNA sequencing with serum antibody peptide sequencing and antibody synthesis, rare circulating B lymphocytes making pathogenic autoantibodies were found to comprise clonal trees accumulating mutations. Lymphoma driver mutations in genes regulating B cell proliferation and V(D)J mutation (CARD11, TNFAIP3, CCND3, ID3, BTG2, and KLHL6) were present in rogue B cells producing the pathogenic autoantibody. Antibody V(D)J mutations conferred pathogenicity by causing the antigen-bound autoantibodies to undergo phase transition to insoluble aggregates at lower temperatures. These results reveal a pre-neoplastic stage in human lymphomagenesis and a cascade of somatic mutations leading to an iconic pathogenic autoantibody.

B-cell receptor reconstruction from single-cell RNA-seq with VDJPuzzle.

Motivation: The B-cell receptor (BCR) performs essential functions for the adaptive immune system including recognition of pathogen-derived antigens. The vast repertoire and adaptive variation of BCR sequences due to V(D)J recombination and somatic hypermutation necessitates single-cell characterization of BCR sequences. Single-cell RNA sequencing presents the opportunity for simultaneous capture of paired BCR heavy and light chains and the transcriptomic signature. Results: We developed VDJPuzzle, a novel bioinformatic tool that reconstructs productive, full-length B-cell receptor sequences of both heavy and light chains and extract somatic mutations on the VDJ region. VDJPuzzle successfully reconstructed BCRs from 100% (n=117) human and 96.5% (n=200) murine B cells. The reconstructed BCRs were successfully validated with single-cell Sanger sequencing. Availability and implementation: VDJPuzzle is available at https://bitbucket.org/kirbyvisp/vdjpuzzle2. Supplementary information: Supplementary data are available at Bioinformatics online.

Principles of long noncoding RNA evolution derived from direct comparison of transcriptomes in 17 species.

The inability to predict long noncoding RNAs from genomic sequence has impeded the use of comparative genomics for studying their biology. Here, we develop methods that use RNA sequencing (RNA-seq) data to annotate the transcriptomes of 16 vertebrates and the echinoid sea urchin, uncovering thousands of previously unannotated genes, most of which produce long intervening noncoding RNAs (lincRNAs). Although in each species, >70% of lincRNAs cannot be traced to homologs in species that diverged >50 million years ago, thousands of human lincRNAs have homologs with similar expression patterns in other species. These homologs share short, 5'-biased patches of sequence conservation nested in exonic architectures that have been extensively rewired, in part by transposable element exonization. Thus, over a thousand human lincRNAs are likely to have conserved functions in mammals, and hundreds beyond mammals, but those functions require only short patches of specific sequences and can tolerate major changes in gene architecture.

Sequencing the cap-snatching repertoire of H1N1 influenza provides insight into the mechanism of viral transcription initiation.

The influenza polymerase cleaves host RNAs ∼10-13 nucleotides downstream of their 5' ends and uses this capped fragment to prime viral mRNA synthesis. To better understand this process of cap snatching, we used high-throughput sequencing to determine the 5' ends of A/WSN/33 (H1N1) influenza mRNAs. The sequences provided clear evidence for nascent-chain realignment during transcription initiation and revealed a strong influence of the viral template on the frequency of realignment. After accounting for the extra nucleotides inserted through realignment, analysis of the capped fragments indicated that the different viral mRNAs were each prepended with a common set of sequences and that the polymerase often cleaved host RNAs after a purine and often primed transcription on a single base pair to either the terminal or penultimate residue of the viral template. We also developed a bioinformatic approach to identify the targeted host transcripts despite limited information content within snatched fragments and found that small nuclear RNAs and small nucleolar RNAs contributed the most abundant capped leaders. These results provide insight into the mechanism of viral transcription initiation and reveal the diversity of the cap-snatched repertoire, showing that noncoding transcripts as well as mRNAs are used to make influenza mRNAs.

Global analyses of the effect of different cellular contexts on microRNA targeting.

MicroRNA (miRNA) regulation clearly impacts animal development, but the extent to which development-with its resulting diversity of cellular contexts-impacts miRNA regulation is unclear. Here, we compared cohorts of genes repressed by the same miRNAs in different cell lines and tissues and found that target repertoires were largely unaffected, with secondary effects explaining most of the differential responses detected. Outliers resulting from differential direct targeting were often attributable to alternative 3' UTR isoform usage that modulated the presence of miRNA sites. More inclusive examination of alternative 3' UTR isoforms revealed that they influence ∼10% of predicted targets when comparing any two cell types. Indeed, considering alternative 3' UTR isoform usage improved prediction of targeting efficacy significantly beyond the improvements observed when considering constitutive isoform usage. Thus, although miRNA targeting is remarkably consistent in different cell types, considering the 3' UTR landscape helps predict targeting efficacy and explain differential regulation that is observed.

Extensive alternative polyadenylation during zebrafish development.

The post-transcriptional fate of messenger RNAs (mRNAs) is largely dictated by their 3' untranslated regions (3' UTRs), which are defined by cleavage and polyadenylation (CPA) of pre-mRNAs. We used poly(A)-position profiling by sequencing (3P-seq) to map poly(A) sites at eight developmental stages and tissues in the zebrafish. Analysis of over 60 million 3P-seq reads substantially increased and improved existing 3' UTR annotations, resulting in confidently identified 3' UTRs for >79% of the annotated protein-coding genes in zebrafish. mRNAs from most zebrafish genes undergo alternative CPA, with those from more than a thousand genes using different dominant 3' UTRs at different stages. These included one of the poly(A) polymerase genes, for which alternative CPA reinforces its repression in the ovary. 3' UTRs tend to be shortest in the ovaries and longest in the brain. Isoforms with some of the shortest 3' UTRs are highly expressed in the ovary, yet absent in the maternally contributed RNAs of the embryo, perhaps because their 3' UTRs are too short to accommodate a uridine-rich motif required for stability of the maternal mRNA. At 2 h post-fertilization, thousands of unique poly(A) sites appear at locations lacking a typical polyadenylation signal, which suggests a wave of widespread cytoplasmic polyadenylation of mRNA degradation intermediates. Our insights into the identities, formation, and evolution of zebrafish 3' UTRs provide a resource for studying gene regulation during vertebrate development.

Biodegradation of polyester polyurethane by endophytic fungi.

Bioremediation is an important approach to waste reduction that relies on biological processes to break down a variety of pollutants. This is made possible by the vast metabolic diversity of the microbial world. To explore this diversity for the breakdown of plastic, we screened several dozen endophytic fungi for their ability to degrade the synthetic polymer polyester polyurethane (PUR). Several organisms demonstrated the ability to efficiently degrade PUR in both solid and liquid suspensions. Particularly robust activity was observed among several isolates in the genus Pestalotiopsis, although it was not a universal feature of this genus. Two Pestalotiopsis microspora isolates were uniquely able to grow on PUR as the sole carbon source under both aerobic and anaerobic conditions. Molecular characterization of this activity suggests that a serine hydrolase is responsible for degradation of PUR. The broad distribution of activity observed and the unprecedented case of anaerobic growth using PUR as the sole carbon source suggest that endophytes are a promising source of biodiversity from which to screen for metabolic properties useful for bioremediation.

Comparison and calibration of transcriptome data from RNA-Seq and tiling arrays.

BACKGROUND: Tiling arrays have been the tool of choice for probing an organism's transcriptome without prior assumptions about the transcribed regions, but RNA-Seq is becoming a viable alternative as the costs of sequencing continue to decrease. Understanding the relative merits of these technologies will help researchers select the appropriate technology for their needs. RESULTS: Here, we compare these two platforms using a matched sample of poly(A)-enriched RNA isolated from the second larval stage of C. elegans. We find that the raw signals from these two technologies are reasonably well correlated but that RNA-Seq outperforms tiling arrays in several respects, notably in exon boundary detection and dynamic range of expression. By exploring the accuracy of sequencing as a function of depth of coverage, we found that about 4 million reads are required to match the sensitivity of two tiling array replicates. The effects of cross-hybridization were analyzed using a "nearest neighbor" classifier applied to array probes; we describe a method for determining potential "black list" regions whose signals are unreliable. Finally, we propose a strategy for using RNA-Seq data as a gold standard set to calibrate tiling array data. All tiling array and RNA-Seq data sets have been submitted to the modENCODE Data Coordinating Center. CONCLUSIONS: Tiling arrays effectively detect transcript expression levels at a low cost for many species while RNA-Seq provides greater accuracy in several regards. Researchers will need to carefully select the technology appropriate to the biological investigations they are undertaking. It will also be important to reconsider a comparison such as ours as sequencing technologies continue to evolve.

The RNA polymerase "switch region" is a target for inhibitors.

The alpha-pyrone antibiotic myxopyronin (Myx) inhibits bacterial RNA polymerase (RNAP). Here, through a combination of genetic, biochemical, and structural approaches, we show that Myx interacts with the RNAP "switch region"--the hinge that mediates opening and closing of the RNAP active center cleft--to prevent interaction of RNAP with promoter DNA. We define the contacts between Myx and RNAP and the effects of Myx on RNAP conformation and propose that Myx functions by interfering with opening of the RNAP active-center cleft during transcription initiation. We further show that the structurally related alpha-pyrone antibiotic corallopyronin (Cor) and the structurally unrelated macrocyclic-lactone antibiotic ripostatin (Rip) function analogously to Myx. The RNAP switch region is distant from targets of previously characterized RNAP inhibitors, and, correspondingly, Myx, Cor, and Rip do not exhibit crossresistance with previously characterized RNAP inhibitors. The RNAP switch region is an attractive target for identification of new broad-spectrum antibacterial therapeutic agents.