A comprehensive SARS-CoV-2-human protein-protein interactome reveals COVID-19 pathobiology and potential host therapeutic targets.

Studying viral-host protein-protein interactions can facilitate the discovery of therapies for viral infection. We use high-throughput yeast two-hybrid experiments and mass spectrometry to generate a comprehensive SARS-CoV-2-human protein-protein interactome network consisting of 739 high-confidence binary and co-complex interactions, validating 218 known SARS-CoV-2 host factors and revealing 361 novel ones. Our results show the highest overlap of interaction partners between published datasets and of genes differentially expressed in samples from COVID-19 patients. We identify an interaction between the viral protein ORF3a and the human transcription factor ZNF579, illustrating a direct viral impact on host transcription. We perform network-based screens of >2,900 FDA-approved or investigational drugs and identify 23 with significant network proximity to SARS-CoV-2 host factors. One of these drugs, carvedilol, shows clinical benefits for COVID-19 patients in an electronic health records analysis and antiviral properties in a human lung cell line infected with SARS-CoV-2. Our study demonstrates the value of network systems biology to understand human-virus interactions and provides hits for further research on COVID-19 therapeutics.

Confined migration induces heterochromatin formation and alters chromatin accessibility.

During migration, cells often squeeze through small constrictions, requiring extensive deformation. We hypothesized that nuclear deformation associated with such confined migration could alter chromatin organization and function. By studying cells migrating through microfluidic devices that mimic interstitial spaces in vivo, we found that confined migration results in increased H3K9me3 and H3K27me3 heterochromatin marks that persist for days. This "confined migration-induced heterochromatin" (CMiH) was distinct from heterochromatin formation during migration initiation. Confined migration decreased chromatin accessibility at intergenic regions near centromeres and telomeres, suggesting heterochromatin spreading from existing sites. Consistent with the overall decrease in accessibility, global transcription was decreased during confined migration. Intriguingly, we also identified increased accessibility at promoter regions of genes linked to chromatin silencing, tumor invasion, and DNA damage response. Inhibiting CMiH reduced migration speed, suggesting that CMiH promotes confined migration. Together, our findings indicate that confined migration induces chromatin changes that regulate cell migration and other functions.

The Drosophila BEAF insulator protein interacts with the polybromo subunit of the PBAP chromatin remodeling complex.

The Drosophila Boundary Element-Associated Factor of 32 kDa (BEAF) binds in promoter regions of a few thousand mostly housekeeping genes. BEAF is implicated in both chromatin domain boundary activity and promoter function, although molecular mechanisms remain elusive. Here, we show that BEAF physically interacts with the polybromo subunit (Pbro) of PBAP, a SWI/SNF-class chromatin remodeling complex. BEAF also shows genetic interactions with Pbro and other PBAP subunits. We examine the effect of this interaction on gene expression and chromatin structure using precision run-on sequencing and micrococcal nuclease sequencing after RNAi-mediated knockdown in cultured S2 cells. Our results are consistent with the interaction playing a subtle role in gene activation. Fewer than 5% of BEAF-associated genes were significantly affected after BEAF knockdown. Most were downregulated, accompanied by fill-in of the promoter nucleosome-depleted region and a slight upstream shift of the +1 nucleosome. Pbro knockdown caused downregulation of several hundred genes and showed a correlation with BEAF knockdown but a better correlation with promoter-proximal GAGA factor binding. Micrococcal nuclease sequencing supports that BEAF binds near housekeeping gene promoters while Pbro is more important at regulated genes. Yet there is a similar general but slight reduction of promoter-proximal pausing by RNA polymerase II and increase in nucleosome-depleted region nucleosome occupancy after knockdown of either protein. We discuss the possibility of redundant factors keeping BEAF-associated promoters active and masking the role of interactions between BEAF and the Pbro subunit of PBAP in S2 cells. We identify Facilitates Chromatin Transcription (FACT) and Nucleosome Remodeling Factor (NURF) as candidate redundant factors.

A comprehensive SARS-CoV-2-human protein-protein interactome network identifies pathobiology and host-targeting therapies for COVID-19.

Physical interactions between viral and host proteins are responsible for almost all aspects of the viral life cycle and the host's immune response. Studying viral-host protein-protein interactions is thus crucial for identifying strategies for treatment and prevention of viral infection. Here, we use high-throughput yeast two-hybrid and affinity purification followed by mass spectrometry to generate a comprehensive SARS-CoV-2-human protein-protein interactome network consisting of both binary and co-complex interactions. We report a total of 739 high-confidence interactions, showing the highest overlap of interaction partners among published datasets as well as the highest overlap with genes differentially expressed in samples (such as upper airway and bronchial epithelial cells) from patients with SARS-CoV-2 infection. Showcasing the utility of our network, we describe a novel interaction between the viral accessory protein ORF3a and the host zinc finger transcription factor ZNF579 to illustrate a SARS-CoV-2 factor mediating a direct impact on host transcription. Leveraging our interactome, we performed network-based drug screens for over 2,900 FDA-approved/investigational drugs and obtained a curated list of 23 drugs that had significant network proximities to SARS-CoV-2 host factors, one of which, carvedilol, showed promising antiviral properties. We performed electronic health record-based validation using two independent large-scale, longitudinal COVID-19 patient databases and found that carvedilol usage was associated with a significantly lowered probability (17%-20%, P < 0.001) of obtaining a SARS-CoV-2 positive test after adjusting various confounding factors. Carvedilol additionally showed anti-viral activity against SARS-CoV-2 in a human lung epithelial cell line [half maximal effective concentration (EC 50 ) value of 4.1 µM], suggesting a mechanism for its beneficial effect in COVID-19. Our study demonstrates the value of large-scale network systems biology approaches for extracting biological insight from complex biological processes.

Roles of transposable elements in the regulation of mammalian transcription.

Transposable elements (TEs) comprise about half of the mammalian genome. TEs often contain sequences capable of recruiting the host transcription machinery, which they use to express their own products and promote transposition. However, the regulatory sequences carried by TEs may affect host transcription long after the TEs have lost the ability to transpose. Recent advances in genome analysis and engineering have facilitated systematic interrogation of the regulatory activities of TEs. In this Review, we discuss diverse mechanisms by which TEs contribute to transcription regulation. Notably, TEs can donate enhancer and promoter sequences that influence the expression of host genes, modify 3D chromatin architecture and give rise to novel regulatory genes, including non-coding RNAs and transcription factors. We discuss how TEs spur regulatory evolution and facilitate the emergence of genetic novelties in mammalian physiology and development. By virtue of their repetitive and interspersed nature, TEs offer unique opportunities to dissect the effects of mutation and genomic context on the function and evolution of cis-regulatory elements. We argue that TE-centric studies hold the key to unlocking general principles of transcription regulation and evolution.

Spatiotemporal localization of proteins in mycobacteria.

Although prokaryotic organisms lack traditional organelles, they must still organize cellular structures in space and time, challenges that different species solve differently. To systematically define the subcellular architecture of mycobacteria, we perform high-throughput imaging of a library of fluorescently tagged proteins expressed in Mycobacterium smegmatis and develop a customized computational pipeline, MOMIA and GEMATRIA, to analyze these data. Our results establish a spatial organization network of over 700 conserved mycobacterial proteins and reveal a coherent localization pattern for many proteins of known function, including those in translation, energy metabolism, cell growth and division, as well as proteins of unknown function. Furthermore, our pipeline exploits morphologic proxies to enable a pseudo-temporal approximation of protein localization and identifies previously uncharacterized cell-cycle-dependent dynamics of essential mycobacterial proteins. Collectively, these data provide a systems perspective on the subcellular organization of mycobacteria and provide tools for the analysis of bacteria with non-standard growth characteristics.

Fluorescence Imaging-Based Discovery of Membrane Domain-Associated Proteins in Mycobacterium smegmatis.

Mycobacteria spatially organize their plasma membrane, and many enzymes involved in envelope biosynthesis associate with a membrane compartment termed the intracellular membrane domain (IMD). The IMD is concentrated in the polar regions of growing cells and becomes less polarized under nongrowing conditions. Because mycobacteria elongate from the poles, the observed polar localization of the IMD during growth likely supports the localized biosynthesis of envelope components. While we have identified more than 300 IMD-associated proteins by proteomic analyses, only a few of these have been verified by independent experimental methods. Furthermore, some IMD-associated proteins may have escaped proteomic identification and remain to be identified. Here, we visually screened an arrayed library of 523 Mycobacterium smegmatis strains, each producing a Dendra2-FLAG-tagged recombinant protein. We identified 29 fusion proteins that showed polar fluorescence patterns characteristic of IMD proteins. Twenty of these had previously been suggested to localize to the IMD based on proteomic data. Of the nine remaining IMD candidate proteins, three were confirmed by biochemical methods to be associated with the IMD. Taken together, this new colocalization strategy is effective in verifying the IMD association of proteins found by proteomic analyses while facilitating the discovery of additional IMD-associated proteins. IMPORTANCE The intracellular membrane domain (IMD) is a membrane subcompartment found in Mycobacterium smegmatis cells. Proteomic analysis of purified IMD identified more than 300 proteins, including enzymes involved in cell envelope biosynthesis. However, proteomics on its own is unlikely to detect every IMD-associated protein because of technical and biological limitations. Here, we describe fluorescent protein colocalization as an alternative, independent approach. Using a combination of fluorescence microscopy, proteomics, and subcellular fractionation, we identified three new proteins associated with the IMD. Such a robust method to rigorously define IMD proteins will benefit future investigations to decipher the synthesis, maintenance, and functions of this membrane domain and help delineate a more general mechanism of subcellular protein localization in mycobacteria.

Evolution of mouse circadian enhancers from transposable elements.

BACKGROUND: Transposable elements are increasingly recognized as a source of cis-regulatory variation. Previous studies have revealed that transposons are often bound by transcription factors and some have been co-opted into functional enhancers regulating host gene expression. However, the process by which transposons mature into complex regulatory elements, like enhancers, remains poorly understood. To investigate this process, we examined the contribution of transposons to the cis-regulatory network controlling circadian gene expression in the mouse liver, a well-characterized network serving an important physiological function. RESULTS: ChIP-seq analyses reveal that transposons and other repeats contribute ~ 14% of the binding sites for core circadian regulators (CRs) including BMAL1, CLOCK, PER1/2, and CRY1/2, in the mouse liver. RSINE1, an abundant murine-specific SINE, is the only transposon family enriched for CR binding sites across all datasets. Sequence analyses and reporter assays reveal that the circadian regulatory activity of RSINE1 stems from the presence of imperfect CR binding motifs in the ancestral RSINE1 sequence. These motifs matured into canonical motifs through point mutations after transposition. Furthermore, maturation occurred preferentially within elements inserted in the proximity of ancestral CR binding sites. RSINE1 also acquired motifs that recruit nuclear receptors known to cooperate with CRs to regulate circadian gene expression specifically in the liver. CONCLUSIONS: Our results suggest that the birth of enhancers from transposons is predicated both by the sequence of the transposon and by the cis-regulatory landscape surrounding their genomic integration site.

Membrane-partitioned cell wall synthesis in mycobacteria.

Many antibiotics target the assembly of cell wall peptidoglycan, an essential, heteropolymeric mesh that encases most bacteria. In rod-shaped bacteria, cell wall elongation is spatially precise yet relies on limited pools of lipid-linked precursors that generate and are attracted to membrane disorder. By tracking enzymes, substrates, and products of peptidoglycan biosynthesis in Mycobacterium smegmatis, we show that precursors are made in plasma membrane domains that are laterally and biochemically distinct from sites of cell wall assembly. Membrane partitioning likely contributes to robust, orderly peptidoglycan synthesis, suggesting that these domains help template peptidoglycan synthesis. The cell wall-organizing protein DivIVA and the cell wall itself promote domain homeostasis. These data support a model in which the peptidoglycan polymer feeds back on its membrane template to maintain an environment conducive to directional synthesis. Our findings are applicable to rod-shaped bacteria that are phylogenetically distant from M. smegmatis, indicating that horizontal compartmentalization of precursors may be a general feature of bacillary cell wall biogenesis.

Recurrent evolution of vertebrate transcription factors by transposase capture.

Genes with novel cellular functions may evolve through exon shuffling, which can assemble novel protein architectures. Here, we show that DNA transposons provide a recurrent supply of materials to assemble protein-coding genes through exon shuffling. We find that transposase domains have been captured-primarily via alternative splicing-to form fusion proteins at least 94 times independently over the course of ~350 million years of tetrapod evolution. We find an excess of transposase DNA binding domains fused to host regulatory domains, especially the Krüppel-associated box (KRAB) domain, and identify four independently evolved KRAB-transposase fusion proteins repressing gene expression in a sequence-specific fashion. The bat-specific KRABINER fusion protein binds its cognate transposons genome-wide and controls a network of genes and cis-regulatory elements. These results illustrate how a transcription factor and its binding sites can emerge.

Pioneer-like factor GAF cooperates with PBAP (SWI/SNF) and NURF (ISWI) to regulate transcription.

Transcriptionally silent genes must be activated throughout development. This requires nucleosomes be removed from promoters and enhancers to allow transcription factor (TF) binding and recruitment of coactivators and RNA polymerase II (Pol II). Specialized pioneer TFs bind nucleosome-wrapped DNA to perform this chromatin opening by mechanisms that remain incompletely understood. Here, we show that GAGA factor (GAF), a Drosophila pioneer-like factor, functions with both SWI/SNF and ISWI family chromatin remodelers to allow recruitment of Pol II and entry to a promoter-proximal paused state, and also to promote Pol II's transition to productive elongation. We found that GAF interacts with PBAP (SWI/SNF) to open chromatin and allow Pol II to be recruited. Importantly, this activity is not dependent on NURF as previously proposed; however, GAF also synergizes with NURF downstream from this process to ensure efficient Pol II pause release and transition to productive elongation, apparently through its role in precisely positioning the +1 nucleosome. These results demonstrate how a single sequence-specific pioneer TF can synergize with remodelers to activate sets of genes. Furthermore, this behavior of remodelers is consistent with findings in yeast and mice, and likely represents general, conserved mechanisms found throughout eukarya.

Transposons take remote control.

A family of retroviral-like elements in the human genome has a pervasive influence on gene expression.

Direct cell-cell contact activates SigM to express the ESX-4 secretion system in Mycobacterium smegmatis.

Conjugal cell-cell contact between strains of Mycobacterium smegmatis induces the esxUT transcript, which encodes the putative primary substrates of the ESAT-6 secretion system 4 (ESX-4) secretion system. This recipient response was required for conjugal transfer of chromosomal DNA from the donor strain. Here we show that the extracytoplasmic σ factor, SigM, is a cell contact-dependent activator of ESX-4 expression and is required for conjugal transfer of DNA in the recipient strain. The SigM regulon includes genes outside the seven-gene core esx4 locus that we show are also required for conjugation, and we show that some of these SigM-induced proteins likely function through ESX-4. A fluorescent reporter revealed that SigM is specifically activated in recipient cells in direct contact with donor cells. Coculture RNA-seq experiments indicated that SigM regulon induction occurred early and before transconjugants are detected. This work supports a model wherein donor contact with the recipient cell surface inactivates the transmembrane anti-SigM, thereby releasing SigM. Free SigM induces an extended ESX-4 secretion system, resulting in changes that facilitate chromosomal transfer. The contact-dependent inactivation of an extracytoplasmic σ-factor that tightly controls ESX-4 activity suggests a mechanism dedicated to detect, and appropriately respond to, external stimuli from mycobacteria.

Application of Distributive Conjugal DNA Transfer in Mycobacterium smegmatis To Establish a Genome-Wide Synthetic Genetic Array.

Genetic redundancy can obscure phenotypic effects of single-gene mutations. Two individual mutations may be viable separately but are lethal when combined, thus synthetically linking the two gene products in an essential process. Synthetic genetic arrays (SGAs), in which defined mutations are combined, provide a powerful approach to identify novel genetic interactions and redundant pathways. A genome-scale SGA can offer an initial assignment of function to hypothetical genes by uncovering interactions with known genes or pathways. Here, we take advantage of the chromosomal conjugation system of Mycobacterium smegmatis to combine individual donor and recipient mutations on a genome-wide scale. We demonstrated the feasibility of a high-throughput mycobacterial SGA (mSGA) screen by using mutants of esx3, fxbA, and recA as query genes, which were combined with an arrayed library of transposon mutants by conjugation. The mSGA identified interacting genes that we had predicted and, most importantly, identified novel interacting genes-encoding both proteins and a noncoding RNA (ncRNA). In combination with other molecular genetic approaches, the mSGA has great potential to both reduce the high number of conserved hypothetical protein annotations in mycobacterial genomes and further define mycobacterial pathways and gene interactions.IMPORTANCE Mycobacterium smegmatis is the model organism of choice for the study of mycobacterial pathogens, because it is a fast-growing nonpathogenic species harboring many genes that are conserved throughout mycobacteria. In this work, we describe a synthetic genetic array (mSGA) approach for M. smegmatis, which combines mutations on a genome-wide scale with high efficiency. Analysis of the double mutant strains enables the identification of interacting genes and pathways that are normally hidden by redundant biological pathways. The mSGA is a powerful genetic tool that enables functions to be assigned to the many conserved hypothetical genes found in all mycobacterial species.

Pomegranate Intake Protects Against Genomic Instability Induced by Medical X-rays In Vivo in Mice.

Ionizing radiation (IR) is a well-documented human carcinogen. The increased use of IR in medical procedures has doubled the annual radiation dose and may increase cancer risk. Genomic instability is an intermediate lesion in IR-induced cancer. We examined whether pomegranate extract (PE) suppresses genomic instability induced by x-rays. Mice were treated orally with PE and exposed to an x-ray dose of 2 Gy. PE intake suppressed x-ray-induced DNA double-strand breaks (DSBs) in peripheral blood and chromosomal damage in bone marrow. We hypothesized that PE-mediated protection against x-ray-induced damage may be due to the upregulation of DSB repair and antioxidant enzymes and/or increase in glutathione (GSH) levels. We found that expression of DSB repair genes was not altered (Nbs1 and Rad50) or was reduced (Mre11, DNA-PKcs, Ku80, Rad51, Rad52 and Brca2) in the liver of PE-treated mice. Likewise, mRNA levels of antioxidant enzymes were reduced (Gpx1, Cat, and Sod2) or were not altered (HO-1 and Sod1) as a function of PE treatment. In contrast, PE-treated mice with and without IR exposure displayed higher hepatic GSH concentrations than controls. Thus, ingestion of pomegranate polyphenols is associated with inhibition of x-ray-induced genomic instability and elevated GSH, which may reduce cancer risk.