Hybrid Gene Origination Creates Human-Virus Chimeric Proteins during Infection.

RNA viruses are a major human health threat. The life cycles of many highly pathogenic RNA viruses like influenza A virus (IAV) and Lassa virus depends on host mRNA, because viral polymerases cleave 5'-m7G-capped host transcripts to prime viral mRNA synthesis ("cap-snatching"). We hypothesized that start codons within cap-snatched host transcripts could generate chimeric human-viral mRNAs with coding potential. We report the existence of this mechanism of gene origination, which we named "start-snatching." Depending on the reading frame, start-snatching allows the translation of host and viral "untranslated regions" (UTRs) to create N-terminally extended viral proteins or entirely novel polypeptides by genetic overprinting. We show that both types of chimeric proteins are made in IAV-infected cells, generate T cell responses, and contribute to virulence. Our results indicate that during infection with IAV, and likely a multitude of other human, animal and plant viruses, a host-dependent mechanism allows the genesis of hybrid genes.

Genomic Decoding of Neuronal Depolarization by Stimulus-Specific NPAS4 Heterodimers.

Cells regulate gene expression in response to salient external stimuli. In neurons, depolarization leads to the expression of inducible transcription factors (ITFs) that direct subsequent gene regulation. Depolarization encodes both a neuron's action potential (AP) output and synaptic inputs, via excitatory postsynaptic potentials (EPSPs). However, it is unclear if distinct types of electrical activity can be transformed by an ITF into distinct modes of genomic regulation. Here, we show that APs and EPSPs in mouse hippocampal neurons trigger two spatially segregated and molecularly distinct induction mechanisms that lead to the expression of the ITF NPAS4. These two pathways culminate in the formation of stimulus-specific NPAS4 heterodimers that exhibit distinct DNA binding patterns. Thus, NPAS4 differentially communicates increases in a neuron's spiking output and synaptic inputs to the nucleus, enabling gene regulation to be tailored to the type of depolarizing activity along the somato-dendritic axis of a neuron.

Transcription Elongation Can Affect Genome 3D Structure.

How transcription affects genome 3D organization is not well understood. We found that during influenza A (IAV) infection, rampant transcription rapidly reorganizes host cell chromatin interactions. These changes occur at the ends of highly transcribed genes, where global inhibition of transcription termination by IAV NS1 protein causes readthrough transcription for hundreds of kilobases. In these readthrough regions, elongating RNA polymerase II disrupts chromatin interactions by inducing cohesin displacement from CTCF sites, leading to locus decompaction. Readthrough transcription into heterochromatin regions switches them from the inert (B) to the permissive (A) chromatin compartment and enables transcription factor binding. Data from non-viral transcription stimuli show that transcription similarly affects cohesin-mediated chromatin contacts within gene bodies. Conversely, inhibition of transcription elongation allows cohesin to accumulate at previously transcribed intragenic CTCF sites and to mediate chromatin looping and compaction. Our data indicate that transcription elongation by RNA polymerase II remodels genome 3D architecture.

Primate-specific ORF0 contributes to retrotransposon-mediated diversity.

LINE-1 retrotransposons are fast-evolving mobile genetic entities that play roles in gene regulation, pathological conditions, and evolution. Here, we show that the primate LINE-1 5'UTR contains a primate-specific open reading frame (ORF) in the antisense orientation that we named ORF0. The gene product of this ORF localizes to promyelocytic leukemia-adjacent nuclear bodies. ORF0 is present in more than 3,000 loci across human and chimpanzee genomes and has a promoter and a conserved strong Kozak sequence that supports translation. By virtue of containing two splice donor sites, ORF0 can also form fusion proteins with proximal exons. ORF0 transcripts are readily detected in induced pluripotent stem (iPS) cells from both primate species. Capped and polyadenylated ORF0 mRNAs are present in the cytoplasm, and endogenous ORF0 peptides are identified upon proteomic analysis. Finally, ORF0 enhances LINE-1 mobility. Taken together, these results suggest a role for ORF0 in retrotransposon-mediated diversity.

Position-dependent function of human sequence-specific transcription factors.

Patterns of transcriptional activity are encoded in our genome through regulatory elements such as promoters or enhancers that, paradoxically, contain similar assortments of sequence-specific transcription factor (TF) binding sites1-3. Knowledge of how these sequence motifs encode multiple, often overlapping, gene expression programs is central to understanding gene regulation and how mutations in non-coding DNA manifest in disease4,5. Here, by studying gene regulation from the perspective of individual transcription start sites (TSSs), using natural genetic variation, perturbation of endogenous TF protein levels and massively parallel analysis of natural and synthetic regulatory elements, we show that the effect of TF binding on transcription initiation is position dependent. Analysing TF-binding-site occurrences relative to the TSS, we identified several motifs with highly preferential positioning. We show that these patterns are a combination of a TF's distinct functional profiles-many TFs, including canonical activators such as NRF1, NFY and Sp1, activate or repress transcription initiation depending on their precise position relative to the TSS. As such, TFs and their spacing collectively guide the site and frequency of transcription initiation. More broadly, these findings reveal how similar assortments of TF binding sites can generate distinct gene regulatory outcomes depending on their spatial configuration and how DNA sequence polymorphisms may contribute to transcription variation and disease and underscore a critical role for TSS data in decoding the regulatory information of our genome.

B cell super-enhancers and regulatory clusters recruit AID tumorigenic activity.

The antibody gene mutator activation-induced cytidine deaminase (AID) promiscuously damages oncogenes, leading to chromosomal translocations and tumorigenesis. Why nonimmunoglobulin loci are susceptible to AID activity is unknown. Here, we study AID-mediated lesions in the context of nuclear architecture and the B cell regulome. We show that AID targets are not randomly distributed across the genome but are predominantly grouped within super-enhancers and regulatory clusters. Unexpectedly, in these domains, AID deaminates active promoters and eRNA(+) enhancers interconnected in some instances over megabases of linear chromatin. Using genome editing, we demonstrate that 3D-linked targets cooperate to recruit AID-mediated breaks. Furthermore, a comparison of hypermutation in mouse B cells, AID-induced kataegis in human lymphomas, and translocations in MEFs reveals that AID damages different genes in different cell types. Yet, in all cases, the targets are predominantly associated with topological complex, highly transcribed super-enhancers, demonstrating that these compartments are key mediators of AID recruitment.

Functional landscape of SARS-CoV-2 cellular restriction.

A deficient interferon (IFN) response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection has been implicated as a determinant of severe coronavirus disease 2019 (COVID-19). To identify the molecular effectors that govern IFN control of SARS-CoV-2 infection, we conducted a large-scale gain-of-function analysis that evaluated the impact of human IFN-stimulated genes (ISGs) on viral replication. A limited subset of ISGs were found to control viral infection, including endosomal factors inhibiting viral entry, RNA binding proteins suppressing viral RNA synthesis, and a highly enriched cluster of endoplasmic reticulum (ER)/Golgi-resident ISGs inhibiting viral assembly/egress. These included broad-acting antiviral ISGs and eight ISGs that specifically inhibited SARS-CoV-2 and SARS-CoV-1 replication. Among the broad-acting ISGs was BST2/tetherin, which impeded viral release and is antagonized by SARS-CoV-2 Orf7a protein. Overall, these data illuminate a set of ISGs that underlie innate immune control of SARS-CoV-2/SARS-CoV-1 infection, which will facilitate the understanding of host determinants that impact disease severity and offer potential therapeutic strategies for COVID-19.

The selection and function of cell type-specific enhancers.

The human body contains several hundred cell types, all of which share the same genome. In metazoans, much of the regulatory code that drives cell type-specific gene expression is located in distal elements called enhancers. Although mammalian genomes contain millions of potential enhancers, only a small subset of them is active in a given cell type. Cell type-specific enhancer selection involves the binding of lineage-determining transcription factors that prime enhancers. Signal-dependent transcription factors bind to primed enhancers, which enables these broadly expressed factors to regulate gene expression in a cell type-specific manner. The expression of genes that specify cell type identity and function is associated with densely spaced clusters of active enhancers known as super-enhancers. The functions of enhancers and super-enhancers are influenced by, and affect, higher-order genomic organization.

Breakthrough infections by SARS-CoV-2 variants boost cross-reactive hybrid immune responses in mRNA-vaccinated Golden Syrian hamsters.

Hybrid immunity (vaccination + natural infection) to SARS-CoV-2 provides superior protection to re-infection. We performed immune profiling studies during breakthrough infections in mRNA-vaccinated hamsters to evaluate hybrid immunity induction. The mRNA vaccine, BNT162b2, was dosed to induce binding antibody titers against ancestral spike, but inefficient serum virus neutralization of ancestral SARS-CoV-2 or variants of concern (VoCs). Vaccination reduced morbidity and controlled lung virus titers for ancestral virus and Alpha but allowed breakthrough infections in Beta, Delta and Mu-challenged hamsters. Vaccination primed for T cell responses that were boosted by infection. Infection back-boosted neutralizing antibody responses against ancestral virus and VoCs. Hybrid immunity resulted in more cross-reactive sera, reflected by smaller antigenic cartography distances. Transcriptomics post-infection reflects both vaccination status and disease course and suggests a role for interstitial macrophages in vaccine-mediated protection. Therefore, protection by vaccination, even in the absence of high titers of neutralizing antibodies in the serum, correlates with recall of broadly reactive B- and T-cell responses.

Yeast derlin Dfm1 employs a chaperone-like function to resolve misfolded membrane protein stress.

Protein aggregates are a common feature of diseased and aged cells. Membrane proteins comprise a quarter of the proteome, and yet, it is not well understood how aggregation of membrane proteins is regulated and what effects these aggregates can have on cellular health. We have determined in yeast that the derlin Dfm1 has a chaperone-like activity that influences misfolded membrane protein aggregation. We establish that this function of Dfm1 does not require recruitment of the ATPase Cdc48 and it is distinct from Dfm1's previously identified function in dislocating misfolded membrane proteins from the endoplasmic reticulum (ER) to the cytosol for degradation. Additionally, we assess the cellular impacts of misfolded membrane proteins in the absence of Dfm1 and determine that misfolded membrane proteins are toxic to cells in the absence of Dfm1 and cause disruptions to proteasomal and ubiquitin homeostasis.

Discovery of SARS-CoV-2 antiviral drugs through large-scale compound repurposing.

The emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in 2019 has triggered an ongoing global pandemic of the severe pneumonia-like disease coronavirus disease 2019 (COVID-19)1. The development of a vaccine is likely to take at least 12-18 months, and the typical timeline for approval of a new antiviral therapeutic agent can exceed 10 years. Thus, repurposing of known drugs could substantially accelerate the deployment of new therapies for COVID-19. Here we profiled a library of drugs encompassing approximately 12,000 clinical-stage or Food and Drug Administration (FDA)-approved small molecules to identify candidate therapeutic drugs for COVID-19. We report the identification of 100 molecules that inhibit viral replication of SARS-CoV-2, including 21 drugs that exhibit dose-response relationships. Of these, thirteen were found to harbour effective concentrations commensurate with probable achievable therapeutic doses in patients, including the PIKfyve kinase inhibitor apilimod2-4 and the cysteine protease inhibitors MDL-28170, Z LVG CHN2, VBY-825 and ONO 5334. Notably, MDL-28170, ONO 5334 and apilimod were found to antagonize viral replication in human pneumocyte-like cells derived from induced pluripotent stem cells, and apilimod also demonstrated antiviral efficacy in a primary human lung explant model. Since most of the molecules identified in this study have already advanced into the clinic, their known pharmacological and human safety profiles will enable accelerated preclinical and clinical evaluation of these drugs for the treatment of COVID-19.

The Dfm1 Derlin Is Required for ERAD Retrotranslocation of Integral Membrane Proteins.

Endoplasmic reticulum (ER)-associated degradation (ERAD) removes misfolded proteins from the ER membrane and lumen by the ubiquitin-proteasome pathway. Retrotranslocation of ubiquitinated substrates to the cytosol is a universal feature of ERAD that requires the Cdc48 AAA-ATPase. Despite intense efforts, the mechanism of ER exit, particularly for integral membrane (ERAD-M) substrates, has remained unclear. Using a self-ubiquitinating substrate (SUS), which undergoes normal retrotranslocation independently of known ERAD factors, and the new SPOCK (single plate orf compendium kit) micro-library to query all yeast genes, we found the rhomboid derlin Dfm1 was required for retrotranslocation of both HRD and DOA ERAD pathway integral membrane substrates. Dfm1 recruited Cdc48 to the ER membrane with its unique SHP motifs, and it catalyzed substrate extraction through its conserved rhomboid motifs. Surprisingly, dfm1Δ can undergo rapid suppression, restoring wild-type ERAD-M. This unexpected suppression explained earlier studies ruling out Dfm1, and it revealed an ancillary ERAD-M retrotranslocation pathway requiring Hrd1.

Heterochromatin-Encoded Satellite RNAs Induce Breast Cancer.

Heterochromatic repetitive satellite RNAs are extensively transcribed in a variety of human cancers, including BRCA1 mutant breast cancer. Aberrant expression of satellite RNAs in cultured cells induces the DNA damage response, activates cell cycle checkpoints, and causes defects in chromosome segregation. However, the mechanism by which satellite RNA expression leads to genomic instability is not well understood. Here we provide evidence that increased levels of satellite RNAs in mammary glands induce tumor formation in mice. Using mass spectrometry, we further show that genomic instability induced by satellite RNAs occurs through interactions with BRCA1-associated protein networks required for the stabilization of DNA replication forks. Additionally, de-stabilized replication forks likely promote the formation of RNA-DNA hybrids in cells expressing satellite RNAs. These studies lay the foundation for developing novel therapeutic strategies that block the effects of non-coding satellite RNAs in cancer cells.

Triclosan administration to humanized UDP-glucuronosyltransferase 1 neonatal mice induces UGT1A1 through a dependence on PPARα and ATF4.

Triclosan (TCS) is an antimicrobial toxicant found in a myriad of consumer products and has been detected in human tissues, including breastmilk. We have evaluated the impact of lactational TCS on UDP-glucuronosyltransferase 1A1 (UGT1A1) expression and bilirubin metabolism in humanized UGT1 (hUGT1) neonatal mice. In hUGT1 mice, expression of the hepatic UGT1A1 gene is developmentally delayed resulting in elevated total serum bilirubin (TSB) levels. We found that newborn hUGT1 mice breastfed or orally treated with TCS presented lower TSB levels along with induction of hepatic UGT1A1. Lactational and oral treatment by gavage with TCS leads to the activation of hepatic nuclear receptors constitutive androstane receptor (CAR), peroxisome proliferator-activated receptor alpha (PPARα), and stress sensor, activating transcription factor 4 (ATF4). When CAR-deficient hUGT1 mice (hUGT1/Car-/-) were treated with TCS, TSB levels were reduced with a robust induction of hepatic UGT1A1, leaving us to conclude that CAR is not tied to UGT1A1 induction. Alternatively, when PPARα-deficient hUGT1 mice (hUGT1/Pparα-/-) were treated with TCS, hepatic UGT1A1 was not induced. Additionally, we had previously demonstrated that TCS is a potent inducer of ATF4, a transcriptional factor linked to the integrated stress response. When ATF4 was deleted in liver of hUGT1 mice (hUGT1/Atf4ΔHep) and these mice treated with TCS, we observed superinduction of hepatic UGT1A1. Oxidative stress genes in livers of hUGT1/Atf4ΔHep treated with TCS were increased, suggesting that ATF4 protects liver from excessive oxidative stress. The increase oxidative stress may be associated with superinduction of UGT1A1. The expression of ATF4 in neonatal hUGT1 hepatic tissue may play a role in the developmental repression of UGT1A1.

MHC class II proteins mediate cross-species entry of bat influenza viruses.

Zoonotic influenza A viruses of avian origin can cause severe disease in individuals, or even global pandemics, and thus pose a threat to human populations. Waterfowl and shorebirds are believed to be the reservoir for all influenza A viruses, but this has recently been challenged by the identification of novel influenza A viruses in bats1,2. The major bat influenza A virus envelope glycoprotein, haemagglutinin, does not bind the canonical influenza A virus receptor, sialic acid or any other glycan1,3,4, despite its high sequence and structural homology with conventional haemagglutinins. This functionally uncharacterized plasticity of the bat influenza A virus haemagglutinin means the tropism and zoonotic potential of these viruses has not been fully determined. Here we show, using transcriptomic profiling of susceptible versus non-susceptible cells in combination with genome-wide CRISPR-Cas9 screening, that the major histocompatibility complex class II (MHC-II) human leukocyte antigen DR isotype (HLA-DR) is an essential entry determinant for bat influenza A viruses. Genetic ablation of the HLA-DR α-chain rendered cells resistant to infection by bat influenza A virus, whereas ectopic expression of the HLA-DR complex in non-susceptible cells conferred susceptibility. Expression of MHC-II from different bat species, pigs, mice or chickens also conferred susceptibility to infection. Notably, the infection of mice with bat influenza A virus resulted in robust virus replication in the upper respiratory tract, whereas mice deficient for MHC-II were resistant. Collectively, our data identify MHC-II as a crucial entry mediator for bat influenza A viruses in multiple species, which permits a broad vertebrate tropism.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation.

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Oral arsenic administration to humanizedUDP-glucuronosyltransferase1 neonatal mice induces UGT1A1 through a dependence on Nrf2 and PXR.

Inorganic arsenic (iAs) is an environmental toxicant that can lead to severe health consequences, which can be exacerbated if exposure occurs early in development. Here, we evaluated the impact of oral iAs treatment on UDP-glucuronosyltransferase 1A1 (UGT1A1) expression and bilirubin metabolism in humanized UGT1 (hUGT1) mice. We found that oral administration of iAs to neonatal hUGT1 mice that display severe neonatal hyperbilirubinemia leads to induction of intestinal UGT1A1 and a reduction in total serum bilirubin values. Oral iAs administration accelerates neonatal intestinal maturation, an event that is directly associated with UGT1A1 induction. As a reactive oxygen species producer, oral iAs treatment activated the Keap-Nrf2 pathway in the intestinal tract and liver. When Nrf2-deficient hUGT1 mice (hUGT1/Nrf2-/-) were treated with iAs, it was shown that activated Nrf2 contributed significantly toward intestinal maturation and UGT1A1 induction. However, hepatic UGT1A1 was not induced upon iAs exposure. We previously demonstrated that the nuclear receptor PXR represses liver UGT1A1 in neonatal hUGT1 mice. When PXR was deleted in hUGT1 mice (hUGT1/Pxr-/-), derepression of UGT1A1 was evident in both liver and intestinal tissue in neonates. Furthermore, when neonatal hUGT1/Pxr-/- mice were treated with iAs, UGT1A1 was superinduced in both tissues, confirming PXR release derepressed key regulatory elements on the gene that could be activated by iAs exposure. With iAs capable of generating reactive oxygen species in both liver and intestinal tissue, we conclude that PXR deficiency in neonatal hUGT1/Pxr-/- mice allows greater access of activated transcriptional modifiers such as Nrf2 leading to superinduction of UGT1A1.

Global changes in the nuclear positioning of genes and intra- and interdomain genomic interactions that orchestrate B cell fate.

The genome is folded into domains located in compartments that are either transcriptionally inert or transcriptionally permissive. Here we used genome-wide strategies to characterize domains during B cell development. Structured interaction matrix analysis showed that occupancy by the architectural protein CTCF was associated mainly with intradomain interactions, whereas sites bound by the histone acetyltransferase p300 or the transcription factors E2A or PU.1 were associated with intra- and interdomain interactions that are developmentally regulated. We identified a spectrum of genes that switched nuclear location during early B cell development. In progenitor cells, the transcriptionally inactive locus encoding early B cell factor (Ebf1) was sequestered at the nuclear lamina, which thereby preserved their multipotency. After development into the pro-B cell stage, Ebf1 and other genes switched compartments to establish new intra- and interdomain interactions associated with a B lineage-specific transcription signature.