Comparison of antiviral responses in two bat species reveals conserved and divergent innate immune pathways.

Bats host a range of disease-causing viruses without displaying clinical symptoms. The mechanisms behind this are a continuous source of interest. Here, we studied the antiviral response in the Egyptian fruit bat and Kuhl's pipistrelle, representing two subordinal clades. We profiled the antiviral response in fibroblasts using RNA sequencing and compared bat with primate and rodent responses. Both bats upregulate similar genes; however, a subset of these genes is transcriptionally divergent between them. These divergent genes also evolve rapidly in sequence, have specific promoter architectures, and are associated with programs underlying tolerance and resistance. Finally, we characterized antiviral genes that expanded in bats, with duplicates diverging in sequence and expression. Our study reveals a largely conserved antiviral program across bats and points to a set of genes that rapidly evolve through multiple mechanisms. These can contribute to bat adaptation to viral infection and provide directions to understanding the mechanisms behind it.

Mapping interindividual dynamics of innate immune response at single-cell resolution.

Common genetic variants across individuals modulate the cellular response to pathogens and are implicated in diverse immune pathologies, yet how they dynamically alter the response upon infection is not well understood. Here, we triggered antiviral responses in human fibroblasts from 68 healthy donors, and profiled tens of thousands of cells using single-cell RNA-sequencing. We developed GASPACHO (GAuSsian Processes for Association mapping leveraging Cell HeterOgeneity), a statistical approach designed to identify nonlinear dynamic genetic effects across transcriptional trajectories of cells. This approach identified 1,275 expression quantitative trait loci (local false discovery rate 10%) that manifested during the responses, many of which were colocalized with susceptibility loci identified by genome-wide association studies of infectious and autoimmune diseases, including the OAS1 splicing quantitative trait locus in a COVID-19 susceptibility locus. In summary, our analytical approach provides a unique framework for delineation of the genetic variants that shape a wide spectrum of transcriptional responses at single-cell resolution.

Tissue memory relies on stem cell priming in distal undamaged areas.

Epithelial cells that participated in wound repair elicit a more efficient response to future injuries, which is believed to be locally restricted. Here we show that cell adaptation resulting from a localized tissue damage has a wide spatial impact at a scale not previously appreciated. We demonstrate that a specific stem cell population, distant from the original injury, originates long-lasting wound memory progenitors residing in their own niche. Notably, these distal memory cells have not taken part in the first healing but become intrinsically pre-activated through priming. This cell state, maintained at the chromatin and transcriptional level, leads to an enhanced wound repair that is partially recapitulated through epigenetic perturbation. Importantly wound memory has long-term harmful consequences, exacerbating tumourigenesis. Overall, we show that sub-organ-scale adaptation to injury relies on spatially organized memory-dedicated progenitors, characterized by an actionable cell state that establishes an epigenetic field cancerization and predisposes to tumour onset.

Promoter evolution of mammalian gene duplicates.

BACKGROUND: Gene duplication is thought to be a central process in evolution to gain new functions. The factors that dictate gene retention following duplication as well paralog gene divergence in sequence, expression and function have been extensively studied. However, relatively little is known about the evolution of promoter regions of gene duplicates and how they influence gene duplicate divergence. Here, we focus on promoters of paralog genes, comparing their similarity in sequence, in the sets of transcription factors (TFs) that bind them, and in their overall promoter architecture. RESULTS: We observe that promoters of recent duplications display higher sequence similarity between them and that sequence similarity rapidly declines between promoters of more ancient paralogs. In contrast, similarity in cis-regulation, as measured by the set of TFs that bind promoters of both paralogs, does not simply decrease with time from duplication and is instead related to promoter architecture-paralogs with CpG Islands (CGIs) in their promoters share a greater fraction of TFs, while CGI-less paralogs are more divergent in their TF binding set. Focusing on recent duplication events and partitioning them by their duplication mechanism enables us to uncover promoter properties associated with gene retention, as well as to characterize the evolution of promoters of newly born genes: In recent retrotransposition-mediated duplications, we observe asymmetry in cis-regulation of paralog pairs: Retrocopy genes are lowly expressed and their promoters are bound by fewer TFs and are depleted of CGIs, in comparison with the original gene copy. Furthermore, looking at recent segmental duplication regions in primates enable us to compare successful retentions versus loss of duplicates, showing that duplicate retention is associated with fewer TFs and with CGI-less promoter architecture. CONCLUSIONS: In this work, we profiled promoters of gene duplicates and their inter-paralog divergence. We also studied how their characteristics are associated with duplication time and duplication mechanism, as well as with the fate of these duplicates. These results underline the importance of cis-regulatory mechanisms in shaping the evolution of new genes and their fate following duplication.

Distinct gene programs underpinning disease tolerance and resistance in influenza virus infection.

When challenged with an invading pathogen, the host-defense response is engaged to eliminate the pathogen (resistance) and to maintain health in the presence of the pathogen (disease tolerance). However, the identification of distinct molecular programs underpinning disease tolerance and resistance remained obscure. We exploited transcriptional and physiological monitoring across 33 mouse strains, during in vivo influenza virus infection, to identify two host-defense gene programs-one is associated with hallmarks of disease tolerance and the other with hallmarks of resistance. Both programs constitute generic responses in multiple mouse and human cell types. Our study describes the organizational principles of these programs and validates Arhgdia as a regulator of disease-tolerance states in epithelial cells. We further reveal that the baseline disease-tolerance state in peritoneal macrophages is associated with the pathophysiological response to injury and infection. Our framework provides a paradigm for the understanding of disease tolerance and resistance at the molecular level.

Rapidly evolving viral motifs mostly target biophysically constrained binding pockets of host proteins.

Evolutionary changes in host-virus interactions can alter the course of infection, but the biophysical and regulatory constraints that shape interface evolution remain largely unexplored. Here, we focus on viral mimicry of host-like motifs that allow binding to host domains and modulation of cellular pathways. We observe that motifs from unrelated viruses preferentially target conserved, widely expressed, and highly connected host proteins, enriched with regulatory and essential functions. The interface residues within these host domains are more conserved and bind a larger number of cellular proteins than similar motif-binding domains that are not known to interact with viruses. In contrast, rapidly evolving viral-binding human proteins form few interactions with other cellular proteins and display high tissue specificity, and their interfaces have few inter-residue contacts. Our results distinguish between conserved and rapidly evolving host-virus interfaces and show how various factors limit host capacity to evolve, allowing for efficient viral subversion of host machineries.

Principles of dengue virus evolvability derived from genotype-fitness maps in human and mosquito cells.

Dengue virus (DENV) cycles between mosquito and mammalian hosts. To examine how DENV populations adapt to these different host environments, we used serial passage in human and mosquito cell lines and estimated fitness effects for all single-nucleotide variants in these populations using ultra-deep sequencing. This allowed us to determine the contributions of beneficial and deleterious mutations to the collective fitness of the population. Our analysis revealed that the continuous influx of a large burden of deleterious mutations counterbalances the effect of rare, host-specific beneficial mutations to shape the path of adaptation. Beneficial mutations preferentially map to intrinsically disordered domains in the viral proteome and cluster to defined regions in the genome. These phenotypically redundant adaptive alleles may facilitate host-specific DENV adaptation. Importantly, the evolutionary constraints described in our simple system mirror trends observed across DENV and Zika strains, indicating it recapitulates key biophysical and biological constraints shaping long-term viral evolution.

Viruses are constantly evolving as a result of mutations in their genetic material and environmental pressures. Viruses switching between insects and mammals face unique evolutionary pressures because they must retain their ability to infect both types of organisms. Yet, the mutations in a virus that may be beneficial in an insect may be different from the ones that may be beneficial in a mammal. Mutations in one host may be even harmful in the other. To learn more about how such viruses thrive as they switch between hosts, Dolan, Taguwa et al. studied the dengue virus, which causes over 390 million infections and over 10,000 deaths each year around the globe. They compared the mutations that occurred as the virus multiplied in human and mosquito cells grown in a laboratory. In the experiments, they used a method called ultra-deep RNA sequencing to identify every change that occurred in the genetic material of the virus each time it multiplied. They determined whether the mutations were beneficial or harmful based on whether they became more common ­ suggesting they helped the virus survive ­ or whether they did not persist because they were likely harmful or even fatal to the virus. The experiments showed that many harmful mutations constantly occur in the virus, in both human and mosquito cells. Beneficial changes happen rarely, and those that do are usually only helpful in one type of cell. Fatal mutations tended to occur in parts of the genetic material that encodes regions in the viral proteins that must remain the same. These structural elements appear to be essential to the virus's survival and unable to undergo change, which makes them good targets for antiviral drugs or vaccines. The techniques used in the study may be useful for investigating other viruses and for understanding the evolutionary constraints on viruses more generally. This may help scientists develop antiviral drugs or vaccines that will remain effective even as viruses continue to evolve and mutate.

Developmental cell programs are co-opted in inflammatory skin disease.

The skin confers biophysical and immunological protection through a complex cellular network established early in embryonic development. We profiled the transcriptomes of more than 500,000 single cells from developing human fetal skin, healthy adult skin, and adult skin with atopic dermatitis and psoriasis. We leveraged these datasets to compare cell states across development, homeostasis, and disease. Our analysis revealed an enrichment of innate immune cells in skin during the first trimester and clonal expansion of disease-associated lymphocytes in atopic dermatitis and psoriasis. We uncovered and validated in situ a reemergence of prenatal vascular endothelial cell and macrophage cellular programs in atopic dermatitis and psoriasis lesional skin. These data illustrate the dynamism of cutaneous immunity and provide opportunities for targeting pathological developmental programs in inflammatory skin diseases.

Single-Cell RNA Sequencing Reveals a Dynamic Stromal Niche That Supports Tumor Growth.

Here, using single-cell RNA sequencing, we examine the stromal compartment in murine melanoma and draining lymph nodes (LNs) at points across tumor development, providing data at http://www.teichlab.org/data/. Naive lymphocytes from LNs undergo activation and clonal expansion within the tumor, before PD1 and Lag3 expression, while tumor-associated myeloid cells promote the formation of a suppressive niche. We identify three temporally distinct stromal populations displaying unique functional signatures, conserved across mouse and human tumors. Whereas "immune" stromal cells are observed in early tumors, "contractile" cells become more prevalent at later time points. Complement component C3 is specifically expressed in the immune population. Its cleavage product C3a supports the recruitment of C3aR+ macrophages, and perturbation of C3a and C3aR disrupts immune infiltration, slowing tumor growth. Our results highlight the power of scRNA-seq to identify complex interplays and increase stromal diversity as a tumor develops, revealing that stromal cells acquire the capacity to modulate immune landscapes from early disease.

Cardelino: computational integration of somatic clonal substructure and single-cell transcriptomes.

Bulk and single-cell DNA sequencing has enabled reconstructing clonal substructures of somatic tissues from frequency and cooccurrence patterns of somatic variants. However, approaches to characterize phenotypic variations between clones are not established. Here we present cardelino (https://github.com/single-cell-genetics/cardelino), a computational method for inferring the clonal tree configuration and the clone of origin of individual cells assayed using single-cell RNA-seq (scRNA-seq). Cardelino flexibly integrates information from imperfect clonal trees inferred based on bulk exome-seq data, and sparse variant alleles expressed in scRNA-seq data. We apply cardelino to a published cancer dataset and to newly generated matched scRNA-seq and exome-seq data from 32 human dermal fibroblast lines, identifying hundreds of differentially expressed genes between cells from different somatic clones. These genes are frequently enriched for cell cycle and proliferation pathways, indicating a role for cell division genes in somatic evolution in healthy skin.

Gene expression variability across cells and species shapes innate immunity.

As the first line of defence against pathogens, cells mount an innate immune response, which varies widely from cell to cell. The response must be potent but carefully controlled to avoid self-damage. How these constraints have shaped the evolution of innate immunity remains poorly understood. Here we characterize the innate immune response's transcriptional divergence between species and variability in expression among cells. Using bulk and single-cell transcriptomics in fibroblasts and mononuclear phagocytes from different species, challenged with immune stimuli, we map the architecture of the innate immune response. Transcriptionally diverging genes, including those that encode cytokines and chemokines, vary across cells and have distinct promoter structures. Conversely, genes that are involved in the regulation of this response, such as those that encode transcription factors and kinases, are conserved between species and display low cell-to-cell variability in expression. We suggest that this expression pattern, which is observed across species and conditions, has evolved as a mechanism for fine-tuned regulation to achieve an effective but balanced response.

Neutralizing Gatad2a-Chd4-Mbd3/NuRD Complex Facilitates Deterministic Induction of Naive Pluripotency.

Mbd3, a member of nucleosome remodeling and deacetylase (NuRD) co-repressor complex, was previously identified as an inhibitor for deterministic induced pluripotent stem cell (iPSC) reprogramming, where up to 100% of donor cells successfully complete the process. NuRD can assume multiple mutually exclusive conformations, and it remains unclear whether this deterministic phenotype can be attributed to a specific Mbd3/NuRD subcomplex. Moreover, since complete ablation of Mbd3 blocks somatic cell proliferation, we aimed to explore functionally relevant alternative ways to neutralize Mbd3-dependent NuRD activity. We identify Gatad2a, a NuRD-specific subunit, whose complete deletion specifically disrupts Mbd3/NuRD repressive activity on the pluripotency circuitry during iPSC differentiation and reprogramming without ablating somatic cell proliferation. Inhibition of Gatad2a facilitates deterministic murine iPSC reprogramming within 8 days. We validate a distinct molecular axis, Gatad2a-Chd4-Mbd3, within Mbd3/NuRD as being critical for blocking reestablishment of naive pluripotency and further highlight signaling-dependent and post-translational modifications of Mbd3/NuRD that influence its interactions and assembly.

Nonspecific yet decisive: Ubiquitination can affect the native-state dynamics of the modified protein.

Ubiquitination is one of the most common post-translational modifications of proteins, and mediates regulated protein degradation among other cellular processes. A fundamental question regarding the mechanism of protein ubiquitination is whether and how ubiquitin affects the biophysical nature of the modified protein. For some systems, it was shown that the position of ubiquitin within the attachment site is quite flexible and ubiquitin does not specifically interact with its substrate. Nevertheless, it was revealed that polyubiquitination can decrease the thermal stability of the modified protein in a site-specific manner because of alterations of the thermodynamic properties of the folded and unfolded states. In this study, we used detailed atomistic simulations to focus on the molecular effects of ubiquitination on the native structure of the modified protein. As a model, we used Ubc7, which is an E2 enzyme whose in vivo ubiquitination process is well characterized and known to lead to degradation. We found that, despite the lack of specific direct interactions between the ubiquitin moiety and Ubc7, ubiquitination decreases the conformational flexibility of certain regions of the substrate Ubc7 protein, which reduces its entropy and thus destabilizes it. The strongest destabilizing effect was observed for systems in which Lys48-linked tetra-ubiquitin was attached to sites used for in vivo degradation. These results reveal how changes in the configurational entropy of the folded state may modulate the stability of the protein's native state. Overall, our results imply that ubiquitination can modify the biophysical properties of the attached protein in the folded state and that, in some proteins, different ubiquitination sites will lead to different biophysical outcomes. We propose that this destabilizing effect of polyubiquitin on the substrate is linked to the functions carried out by the modification, and in particular, regulatory control of protein half-life through proteasomal degradation.

Rapid evolution of virus sequences in intrinsically disordered protein regions.

Nodamura Virus (NoV) is a nodavirus originally isolated from insects that can replicate in a wide variety of hosts, including mammals. Because of their simplicity and ability to replicate in many diverse hosts, NoV, and the Nodaviridae in general, provide a unique window into the evolution of viruses and host-virus interactions. Here we show that the C-terminus of the viral polymerase exhibits extreme structural and evolutionary flexibility. Indeed, fewer than 10 positively charged residues from the 110 amino acid-long C-terminal region of protein A are required to support RNA1 replication. Strikingly, this region can be replaced by completely unrelated protein sequences, yet still produce a functional replicase. Structure predictions, as well as evolutionary and mutational analyses, indicate that the C-terminal region is structurally disordered and evolves faster than the rest of the viral proteome. Thus, the function of an intrinsically unstructured protein region can be independent of most of its primary sequence, conferring both functional robustness and sequence plasticity on the protein. Our results provide an experimental explanation for rapid evolution of unstructured regions, which enables an effective exploration of the sequence space, and likely function space, available to the virus.

Use of host-like peptide motifs in viral proteins is a prevalent strategy in host-virus interactions.

Viruses interact extensively with host proteins, but the mechanisms controlling these interactions are not well understood. We present a comprehensive analysis of eukaryotic linear motifs (ELMs) in 2,208 viral genomes and reveal that viruses exploit molecular mimicry of host-like ELMs to possibly assist in host-virus interactions. Using a statistical genomics approach, we identify a large number of potentially functional ELMs and observe that the occurrence of ELMs is often evolutionarily conserved but not uniform across virus families. Some viral proteins contain multiple types of ELMs, in striking similarity to complex regulatory modules in host proteins, suggesting that ELMs may act combinatorially to assist viral replication. Furthermore, a simple evolutionary model suggests that the inherent structural simplicity of ELMs often enables them to tolerate mutations and evolve quickly. Our findings suggest that ELMs may allow fast rewiring of host-virus interactions, which likely assists rapid viral evolution and adaptation to diverse environments.

NADH binds and stabilizes the 26S proteasomes independent of ATP.

The 26S proteasome is the end point of the ubiquitin- and ATP-dependent degradation pathway. The 26S proteasome complex (26S PC) integrity and function has been shown to be highly dependent on ATP and its homolog nucleotides. We report here that the redox molecule NADH binds the 26S PC and is sufficient in maintaining 26S PC integrity even in the absence of ATP. Five of the 19S proteasome complex subunits contain a putative NADH binding motif (GxGxxG) including the AAA-ATPase subunit, Psmc1 (Rpt2). We demonstrate that recombinant Psmc1 binds NADH via the GxGxxG motif. Introducing the ΔGxGxxG Psmc1 mutant into cells results in reduced NADH-stabilized 26S proteasomes and decreased viability following redox stress induced by the mitochondrial inhibitor rotenone. The newly identified NADH binding of 26S proteasomes advances our understanding of the molecular mechanisms of protein degradation and highlights a new link between protein homeostasis and the cellular metabolic/redox state.

Deterministic direct reprogramming of somatic cells to pluripotency.

Somatic cells can be inefficiently and stochastically reprogrammed into induced pluripotent stem (iPS) cells by exogenous expression of Oct4 (also called Pou5f1), Sox2, Klf4 and Myc (hereafter referred to as OSKM). The nature of the predominant rate-limiting barrier(s) preventing the majority of cells to successfully and synchronously reprogram remains to be defined. Here we show that depleting Mbd3, a core member of the Mbd3/NuRD (nucleosome remodelling and deacetylation) repressor complex, together with OSKM transduction and reprogramming in naive pluripotency promoting conditions, result in deterministic and synchronized iPS cell reprogramming (near 100% efficiency within seven days from mouse and human cells). Our findings uncover a dichotomous molecular function for the reprogramming factors, serving to reactivate endogenous pluripotency networks while simultaneously directly recruiting the Mbd3/NuRD repressor complex that potently restrains the reactivation of OSKM downstream target genes. Subsequently, the latter interactions, which are largely depleted during early pre-implantation development in vivo, lead to a stochastic and protracted reprogramming trajectory towards pluripotency in vitro. The deterministic reprogramming approach devised here offers a novel platform for the dissection of molecular dynamics leading to establishing pluripotency at unprecedented flexibility and resolution.

Stabilization of a protein conferred by an increase in folded state entropy.

Entropic stabilization of native protein structures typically relies on strategies that serve to decrease the entropy of the unfolded state. Here we report, using a combination of experimental and computational approaches, on enhanced thermodynamic stability conferred by an increase in the configurational entropy of the folded state. The enhanced stability is observed upon modifications of a loop region in the enzyme acylphosphatase and is achieved despite significant enthalpy losses. The modifications that lead to increased stability, as well as those that result in destabilization, however, strongly compromise enzymatic activity, rationalizing the preservation of the native loop structure even though it does not provide the protein with maximal stability or kinetic foldability.

Modulation of folding kinetics of repeat proteins: interplay between intra- and interdomain interactions.

Repeat proteins have unique elongated structures that, unlike globular proteins, are quite modular. Despite their simple one-dimensional structure, repeat proteins exhibit intricate folding behavior with a complexity similar to that of globular proteins. Therefore, repeat proteins allow one to quantify fundamental aspects of the biophysics of protein folding. One important feature of repeat proteins is the interfaces between the repeating units. In particular, the distribution of stabilities within and between the repeats was previously suggested to affect their folding characteristics. In this study, we explore how the interface affects folding kinetics and cooperativity by investigating two families of repeat proteins, namely, the Ankyrin and tetratricopeptide repeat proteins, which differ in the number of interfacial contacts that are formed between their units as well as in their folding behavior. By using simple topology-based models, we show that modulating the energetic strength of the interface relative to that of the repeat itself can drastically change the protein stability, folding rate, and cooperativity. By further dissecting the interfacial contacts into several subsets, we isolated the effects of each of these groups on folding kinetics. Our study highlights the importance of interface connectivity in determining the folding behavior.