Ubiquitylation of lipopolysaccharide by RNF213 during bacterial infection.

Ubiquitylation is a widespread post-translational protein modification in eukaryotes and marks bacteria that invade the cytosol as cargo for antibacterial autophagy1-3. The identity of the ubiquitylated substrate on bacteria is unknown. Here we show that the ubiquitin coat on Salmonella that invade the cytosol is formed through the ubiquitylation of a non-proteinaceous substrate, the lipid A moiety of bacterial lipopolysaccharide (LPS), by the E3 ubiquitin ligase ring finger protein 213 (RNF213). RNF213 is a risk factor for moyamoya disease4,5, which is a progressive stenosis of the supraclinoid internal carotid artery that causes stroke (especially in children)6,7. RNF213 restricts the proliferation of cytosolic Salmonella and is essential for the generation of the bacterial ubiquitin coat, both directly (through the ubiquitylation of LPS) and indirectly (through the recruitment of LUBAC, which is a downstream E3 ligase that adds M1-linked ubiquitin chains onto pre-existing ubiquitin coats8). In cells that lack RNF213, bacteria do not attract ubiquitin-dependent autophagy receptors or induce antibacterial autophagy. The ubiquitylation of LPS on Salmonella that invade the cytosol requires the dynein-like core of RNF213, but not its RING domain. Instead, ubiquitylation of LPS relies on an RZ finger in the E3 shell. We conclude that ubiquitylation extends beyond protein substrates and that ubiquitylation of LPS triggers cell-autonomous immunity, and we postulate that non-proteinaceous substances other than LPS may also become ubiquitylated.

Cryo-EM reveals the complex architecture of dynactin's shoulder region and pointed end.

Dynactin is a 1.1 MDa complex that activates the molecular motor dynein for ultra-processive transport along microtubules. In order to do this, it forms a tripartite complex with dynein and a coiled-coil adaptor. Dynactin consists of an actin-related filament whose length is defined by its flexible shoulder domain. Despite previous cryo-EM structures, the molecular architecture of the shoulder and pointed end of the filament is still poorly understood due to the lack of high-resolution information in these regions. Here we combine multiple cryo-EM datasets and define precise masking strategies for particle signal subtraction and 3D classification. This overcomes domain flexibility and results in high-resolution maps into which we can build the shoulder and pointed end. The unique architecture of the shoulder securely houses the p150 subunit and positions the four identical p50 subunits in different conformations to bind dynactin's filament. The pointed end map allows us to build the first structure of p62 and reveals the molecular basis for cargo adaptor binding to different sites at the pointed end.

The Molecular Basis for Ubiquitin and Ubiquitin-like Specificities in Bacterial Effector Proteases.

Pathogenic bacteria rely on secreted effector proteins to manipulate host signaling pathways, often in creative ways. CE clan proteases, specific hydrolases for ubiquitin-like modifications (SUMO and NEDD8) in eukaryotes, reportedly serve as bacterial effector proteins with deSUMOylase, deubiquitinase, or, even, acetyltransferase activities. Here, we characterize bacterial CE protease activities, revealing K63-linkage-specific deubiquitinases in human pathogens, such as Salmonella, Escherichia, and Shigella, as well as ubiquitin/ubiquitin-like cross-reactive enzymes in Chlamydia, Rickettsia, and Xanthomonas. Five crystal structures, including ubiquitin/ubiquitin-like complexes, explain substrate specificities and redefine relationships across the CE clan. Importantly, this work identifies novel family members and provides key discoveries among previously reported effectors, such as the unexpected deubiquitinase activity in Xanthomonas XopD, contributed by an unstructured ubiquitin binding region. Furthermore, accessory domains regulate properties such as subcellular localization, as exemplified by a ubiquitin-binding domain in Salmonella Typhimurium SseL. Our work both highlights and explains the functional adaptations observed among diverse CE clan proteins.

Structure and assembly pathway of the ribosome quality control complex.

During ribosome-associated quality control, stalled ribosomes are split into subunits and the 60S-housed nascent polypeptides are poly-ubiquitinated by Listerin. How this low-abundance ubiquitin ligase targets rare stall-generated 60S among numerous empty 60S is unknown. Here, we show that Listerin specificity for nascent chain-60S complexes depends on nuclear export mediator factor (NEMF). The 3.6 Å cryo-EM structure of a nascent chain-containing 60S-Listerin-NEMF complex revealed that NEMF makes multiple simultaneous contacts with 60S and peptidyl-tRNA to sense nascent chain occupancy. Structural and mutational analyses showed that ribosome-bound NEMF recruits and stabilizes Listerin's N-terminal domain, while Listerin's C-terminal RWD domain directly contacts the ribosome to position the adjacent ligase domain near the nascent polypeptide exit tunnel. Thus, highly specific nascent chain targeting by Listerin is imparted by the avidity gained from a multivalent network of context-specific individually weak interactions, highlighting a new principle of client recognition during protein quality control.

Live imaging RNAi screen reveals genes essential for meiosis in mammalian oocytes.

During fertilization, an egg and a sperm fuse to form a new embryo. Eggs develop from oocytes in a process called meiosis. Meiosis in human oocytes is highly error-prone, and defective eggs are the leading cause of pregnancy loss and several genetic disorders such as Down's syndrome. Which genes safeguard accurate progression through meiosis is largely unclear. Here we develop high-content phenotypic screening methods for the systematic identification of mammalian meiotic genes. We targeted 774 genes by RNA interference within follicle-enclosed mouse oocytes to block protein expression from an early stage of oocyte development onwards. We then analysed the function of several genes simultaneously by high-resolution imaging of chromosomes and microtubules in live oocytes and scored each oocyte quantitatively for 50 phenotypes, generating a comprehensive resource of meiotic gene function. The screen generated an unprecedented annotated data set of meiotic progression in 2,241 mammalian oocytes, which allowed us to analyse systematically which defects are linked to abnormal chromosome segregation during meiosis, identifying progression into anaphase with misaligned chromosomes as well as defects in spindle organization as risk factors. This study demonstrates how high-content screens can be performed in oocytes, and allows systematic studies of meiosis in mammals.

Gene discovery by chemical mutagenesis and whole-genome sequencing in Dictyostelium.

Whole-genome sequencing is a useful approach for identification of chemical-induced lesions, but previous applications involved tedious genetic mapping to pinpoint the causative mutations. We propose that saturation mutagenesis under low mutagenic loads, followed by whole-genome sequencing, should allow direct implication of genes by identifying multiple independent alleles of each relevant gene. We tested the hypothesis by performing three genetic screens with chemical mutagenesis in the social soil amoeba Dictyostelium discoideum Through genome sequencing, we successfully identified mutant genes with multiple alleles in near-saturation screens, including resistance to intense illumination and strong suppressors of defects in an allorecognition pathway. We tested the causality of the mutations by comparison to published data and by direct complementation tests, finding both dominant and recessive causative mutations. Therefore, our strategy provides a cost- and time-efficient approach to gene discovery by integrating chemical mutagenesis and whole-genome sequencing. The method should be applicable to many microbial systems, and it is expected to revolutionize the field of functional genomics in Dictyostelium by greatly expanding the mutation spectrum relative to other common mutagenesis methods.

Terpene synthase genes in eukaryotes beyond plants and fungi: Occurrence in social amoebae.

Terpenes are structurally diverse natural products involved in many ecological interactions. The pivotal enzymes for terpene biosynthesis, terpene synthases (TPSs), had been described only in plants and fungi in the eukaryotic domain. In this report, we systematically analyzed the genome sequences of a broad range of nonplant/nonfungus eukaryotes and identified putative TPS genes in six species of amoebae, five of which are multicellular social amoebae from the order of Dictyosteliida. A phylogenetic analysis revealed that amoebal TPSs are evolutionarily more closely related to fungal TPSs than to bacterial TPSs. The social amoeba Dictyostelium discoideum was selected for functional study of the identified TPSs. D. discoideum grows as a unicellular organism when food is abundant and switches from vegetative growth to multicellular development upon starvation. We found that expression of most D. discoideum TPS genes was induced during development. Upon heterologous expression, all nine TPSs from D. discoideum showed sesquiterpene synthase activities. Some also exhibited monoterpene and/or diterpene synthase activities. Direct measurement of volatile terpenes in cultures of D. discoideum revealed essentially no emission at an early stage of development. In contrast, a bouquet of terpenes, dominated by sesquiterpenes including ß-barbatene and (E,E)-α-farnesene, was detected at the middle and late stages of development, suggesting a development-specific function of volatile terpenes in D. discoideum. The patchy distribution of TPS genes in the eukaryotic domain and the evidence for TPS function in D. discoideum indicate that the TPS genes mediate lineage-specific adaptations.

The GATA transcription factor gene gtaG is required for terminal differentiation in Dictyostelium.

The GATA transcription factor GtaG is conserved in Dictyostelids and essential for terminal differentiation in Dictyostelium discoideum, but its function is not well understood. Here we show that gtaG is expressed in prestalk cells at the anterior region of fingers and in the extending stalk during culmination. The gtaG- phenotype is cell-autonomous in prestalk cells and non-cell-autonomous in prespore cells. Transcriptome analyses reveal that GtaG regulates prestalk gene expression during cell differentiation before culmination and is required for progression into culmination. GtaG-dependent genes include genetic suppressors of the Dd-STATa-defective phenotype as well as Dd-STATa target-genes, including extra cellular matrix genes. We show that GtaG may be involved in the production of two culmination-signaling molecules, cyclic di-GMP and the spore differentiation factor SDF-1 and that addition of c-di-GMP rescues the gtaG- culmination and spore formation deficiencies. We propose that GtaG is a regulator of terminal differentiation that functions in concert with Dd-STATa and controls culmination through regulating c-di-GMP and SDF-1 production in prestalk cells.

Allorecognition, via TgrB1 and TgrC1, mediates the transition from unicellularity to multicellularity in the social amoeba Dictyostelium discoideum.

The social amoeba Dictyostelium discoideum integrates into a multicellular organism when individual starving cells aggregate and form a mound. The cells then integrate into defined tissues and develop into a fruiting body that consists of a stalk and spores. Aggregation is initially orchestrated by waves of extracellular cyclic adenosine monophosphate (cAMP), and previous theory suggested that cAMP and other field-wide diffusible signals mediate tissue integration and terminal differentiation as well. Cooperation between cells depends on an allorecognition system comprising the polymorphic adhesion proteins TgrB1 and TgrC1. Binding between compatible TgrB1 and TgrC1 variants ensures that non-matching cells segregate into distinct aggregates prior to terminal development. Here, we have embedded a small number of cells with incompatible allotypes within fields of developing cells with compatible allotypes. We found that compatibility of the allotype encoded by the tgrB1 and tgrC1 genes is required for tissue integration, as manifested in cell polarization, coordinated movement and differentiation into prestalk and prespore cells. Our results show that the molecules that mediate allorecognition in D. discoideum also control the integration of individual cells into a unified developing organism, and this acts as a gating step for multicellularity.

Proto-genes and de novo gene birth.

Novel protein-coding genes can arise either through re-organization of pre-existing genes or de novo. Processes involving re-organization of pre-existing genes, notably after gene duplication, have been extensively described. In contrast, de novo gene birth remains poorly understood, mainly because translation of sequences devoid of genes, or 'non-genic' sequences, is expected to produce insignificant polypeptides rather than proteins with specific biological functions. Here we formalize an evolutionary model according to which functional genes evolve de novo through transitory proto-genes generated by widespread translational activity in non-genic sequences. Testing this model at the genome scale in Saccharomyces cerevisiae, we detect translation of hundreds of short species-specific open reading frames (ORFs) located in non-genic sequences. These translation events seem to provide adaptive potential, as suggested by their differential regulation upon stress and by signatures of retention by natural selection. In line with our model, we establish that S. cerevisiae ORFs can be placed within an evolutionary continuum ranging from non-genic sequences to genes. We identify ~1,900 candidate proto-genes among S. cerevisiae ORFs and find that de novo gene birth from such a reservoir may be more prevalent than sporadic gene duplication. Our work illustrates that evolution exploits seemingly dispensable sequences to generate adaptive functional innovation.

An inter-species protein-protein interaction network across vast evolutionary distance.

In cellular systems, biophysical interactions between macromolecules underlie a complex web of functional interactions. How biophysical and functional networks are coordinated, whether all biophysical interactions correspond to functional interactions, and how such biophysical-versus-functional network coordination is shaped by evolutionary forces are all largely unanswered questions. Here, we investigate these questions using an "inter-interactome" approach. We systematically probed the yeast and human proteomes for interactions between proteins from these two species and functionally characterized the resulting inter-interactome network. After a billion years of evolutionary divergence, the yeast and human proteomes are still capable of forming a biophysical network with properties that resemble those of intra-species networks. Although substantially reduced relative to intra-species networks, the levels of functional overlap in the yeast-human inter-interactome network uncover significant remnants of co-functionality widely preserved in the two proteomes beyond human-yeast homologs. Our data support evolutionary selection against biophysical interactions between proteins with little or no co-functionality. Such non-functional interactions, however, represent a reservoir from which nascent functional interactions may arise.

Leaps and lulls in the developmental transcriptome of Dictyostelium discoideum.

BACKGROUND: Development of the soil amoeba Dictyostelium discoideum is triggered by starvation. When placed on a solid substrate, the starving solitary amoebae cease growth, communicate via extracellular cAMP, aggregate by tens of thousands and develop into multicellular organisms. Early phases of the developmental program are often studied in cells starved in suspension while cAMP is provided exogenously. Previous studies revealed massive shifts in the transcriptome under both developmental conditions and a close relationship between gene expression and morphogenesis, but were limited by the sampling frequency and the resolution of the methods. RESULTS: Here, we combine the superior depth and specificity of RNA-seq-based analysis of mRNA abundance with high frequency sampling during filter development and cAMP pulsing in suspension. We found that the developmental transcriptome exhibits mostly gradual changes interspersed by a few instances of large shifts. For each time point we treated the entire transcriptome as single phenotype, and were able to characterize development as groups of similar time points separated by gaps. The grouped time points represented gradual changes in mRNA abundance, or molecular phenotype, and the gaps represented times during which many genes are differentially expressed rapidly, and thus the phenotype changes dramatically. Comparing developmental experiments revealed that gene expression in filter developed cells lagged behind those treated with exogenous cAMP in suspension. The high sampling frequency revealed many genes whose regulation is reproducibly more complex than indicated by previous studies. Gene Ontology enrichment analysis suggested that the transition to multicellularity coincided with rapid accumulation of transcripts associated with DNA processes and mitosis. Later development included the up-regulation of organic signaling molecules and co-factor biosynthesis. Our analysis also demonstrated a high level of synchrony among the developing structures throughout development. CONCLUSIONS: Our data describe D. discoideum development as a series of coordinated cellular and multicellular activities. Coordination occurred within fields of aggregating cells and among multicellular bodies, such as mounds or migratory slugs that experience both cell-cell contact and various soluble signaling regimes. These time courses, sampled at the highest temporal resolution to date in this system, provide a comprehensive resource for studies of developmental gene expression.

Exploring GPCR signaling pathway networks as cancer therapeutic targets.

GPCR signaling can contribute to establishing the tumor microenvironment and influence the progression and metabolism of tumors. Arora et al.1 describe a systems-level approach to investigate the patterns of co-expression of GPCR signaling pathway networks across diverse tumors and identify network components that correlate with patient-survival data across different cancer types.

Dynamic conformational changes of a tardigrade group-3 late embryogenesis abundant protein modulate membrane biophysical properties.

A number of intrinsically disordered proteins (IDPs) encoded in stress-tolerant organisms, such as tardigrade, can confer fitness advantage and abiotic stress tolerance when heterologously expressed. Tardigrade-specific disordered proteins including the cytosolic-abundant heat-soluble proteins are proposed to confer stress tolerance through vitrification or gelation, whereas evolutionarily conserved IDPs in tardigrades may contribute to stress tolerance through other biophysical mechanisms. In this study, we characterized the mechanism of action of an evolutionarily conserved, tardigrade IDP, HeLEA1, which belongs to the group-3 late embryogenesis abundant (LEA) protein family. HeLEA1 homologs are found across different kingdoms of life. HeLEA1 is intrinsically disordered in solution but shows a propensity for helical structure across its entire sequence. HeLEA1 interacts with negatively charged membranes via dynamic disorder-to-helical transition, mainly driven by electrostatic interactions. Membrane interaction of HeLEA1 is shown to ameliorate excess surface tension and lipid packing defects. HeLEA1 localizes to the mitochondrial matrix when expressed in yeast and interacts with model membranes mimicking inner mitochondrial membrane. Yeast expressing HeLEA1 shows enhanced tolerance to hyperosmotic stress under nonfermentative growth and increased mitochondrial membrane potential. Evolutionary analysis suggests that although HeLEA1 homologs have diverged their sequences to localize to different subcellular organelles, all homologs maintain a weak hydrophobic moment that is characteristic of weak and reversible membrane interaction. We suggest that such dynamic and weak protein-membrane interaction buffering alterations in lipid packing could be a conserved strategy for regulating membrane properties and represent a general biophysical solution for stress tolerance across the domains of life.

Viral perturbations of host networks reflect disease etiology.

Many human diseases, arising from mutations of disease susceptibility genes (genetic diseases), are also associated with viral infections (virally implicated diseases), either in a directly causal manner or by indirect associations. Here we examine whether viral perturbations of host interactome may underlie such virally implicated disease relationships. Using as models two different human viruses, Epstein-Barr virus (EBV) and human papillomavirus (HPV), we find that host targets of viral proteins reside in network proximity to products of disease susceptibility genes. Expression changes in virally implicated disease tissues and comorbidity patterns cluster significantly in the network vicinity of viral targets. The topological proximity found between cellular targets of viral proteins and disease genes was exploited to uncover a novel pathway linking HPV to Fanconi anemia.

Systematic assessment of prognostic molecular features across cancers.

Precision oncology promises accurate prediction of disease trajectories by utilizing molecular features of tumors. We present a systematic analysis of the prognostic potential of diverse molecular features across large cancer cohorts. We find that the mRNA expression of biologically coherent sets of genes (modules) is substantially more predictive of patient survival than single-locus genomic and transcriptomic aberrations. Extending our analysis beyond existing curated gene modules, we find a large novel class of highly prognostic DNA/RNA cis-regulatory modules associated with dynamic gene expression within cancers. Remarkably, in more than 82% of cancers, modules substantially improve survival stratification compared with conventional clinical factors and prominent genomic aberrations. The prognostic potential of cancer modules generalizes to external cohorts better than conventionally used single-gene features. Finally, a machine-learning framework demonstrates the combined predictive power of multiple modules, yielding prognostic models that perform substantially better than existing histopathological and clinical factors in common use.

Cellular allostatic load is linked to increased energy expenditure and accelerated biological aging.

Stress triggers anticipatory physiological responses that promote survival, a phenomenon termed allostasis. However, the chronic activation of energy-dependent allostatic responses results in allostatic load, a dysregulated state that predicts functional decline, accelerates aging, and increases mortality in humans. The energetic cost and cellular basis for the damaging effects of allostatic load have not been defined. Here, by longitudinally profiling three unrelated primary human fibroblast lines across their lifespan, we find that chronic glucocorticoid exposure increases cellular energy expenditure by ∼60%, along with a metabolic shift from glycolysis to mitochondrial oxidative phosphorylation (OxPhos). This state of stress-induced hypermetabolism is linked to mtDNA instability, non-linearly affects age-related cytokines secretion, and accelerates cellular aging based on DNA methylation clocks, telomere shortening rate, and reduced lifespan. Pharmacologically normalizing OxPhos activity while further increasing energy expenditure exacerbates the accelerated aging phenotype, pointing to total energy expenditure as a potential driver of aging dynamics. Together, our findings define bioenergetic and multi-omic recalibrations of stress adaptation, underscoring increased energy expenditure and accelerated cellular aging as interrelated features of cellular allostatic load.

OxPhos defects cause hypermetabolism and reduce lifespan in cells and in patients with mitochondrial diseases.

Patients with primary mitochondrial oxidative phosphorylation (OxPhos) defects present with fatigue and multi-system disorders, are often lean, and die prematurely, but the mechanistic basis for this clinical picture remains unclear. By integrating data from 17 cohorts of patients with mitochondrial diseases (n = 690) we find evidence that these disorders increase resting energy expenditure, a state termed hypermetabolism. We examine this phenomenon longitudinally in patient-derived fibroblasts from multiple donors. Genetically or pharmacologically disrupting OxPhos approximately doubles cellular energy expenditure. This cell-autonomous state of hypermetabolism occurs despite near-normal OxPhos coupling efficiency, excluding uncoupling as a general mechanism. Instead, hypermetabolism is associated with mitochondrial DNA instability, activation of the integrated stress response (ISR), and increased extracellular secretion of age-related cytokines and metabokines including GDF15. In parallel, OxPhos defects accelerate telomere erosion and epigenetic aging per cell division, consistent with evidence that excess energy expenditure accelerates biological aging. To explore potential mechanisms for these effects, we generate a longitudinal RNASeq and DNA methylation resource dataset, which reveals conserved, energetically demanding, genome-wide recalibrations. Taken together, these findings highlight the need to understand how OxPhos defects influence the energetic cost of living, and the link between hypermetabolism and aging in cells and patients with mitochondrial diseases.

Functional genetic screen identifies ITPR3/calcium/RELB axis as a driver of colorectal cancer metastatic liver colonization.

Metastatic colonization is the primary cause of death from colorectal cancer (CRC). We employed genome-scale in vivo short hairpin RNA (shRNA) screening and validation to identify 26 promoters of CRC liver colonization. Among these genes, we identified a cluster that contains multiple targetable genes, including ITPR3, which promoted liver-metastatic colonization and elicited similar downstream gene expression programs. ITPR3 is a caffeine-sensitive inositol 1,4,5-triphosphate (IP3) receptor that releases calcium from the endoplasmic reticulum and enhanced metastatic colonization by inducing expression of RELB, a transcription factor that is associated with non-canonical NF-κB signaling. Genetic, cell biological, pharmacologic, and clinical association studies revealed that ITPR3 and RELB drive CRC colony formation by promoting cell survival upon substratum detachment or hypoxic exposure. RELB was sufficient to drive colonization downstream of ITPR3. Our findings implicate the ITPR3/calcium/RELB axis in CRC metastatic colony formation and uncover multiple clinico-pathologically associated targetable proteins as drivers of CRC metastatic colonization.

A multi-omics longitudinal aging dataset in primary human fibroblasts with mitochondrial perturbations.

Aging is a process of progressive change. To develop biological models of aging, longitudinal datasets with high temporal resolution are needed. Here we report a multi-omics longitudinal dataset for cultured primary human fibroblasts measured across their replicative lifespans. Fibroblasts were sourced from both healthy donors (n = 6) and individuals with lifespan-shortening mitochondrial disease (n = 3). The dataset includes cytological, bioenergetic, DNA methylation, gene expression, secreted proteins, mitochondrial DNA copy number and mutations, cell-free DNA, telomere length, and whole-genome sequencing data. This dataset enables the bridging of mechanistic processes of aging as outlined by the "hallmarks of aging", with the descriptive characterization of aging such as epigenetic age clocks. Here we focus on bridging the gap for the hallmark mitochondrial metabolism. Our dataset includes measurement of healthy cells, and cells subjected to over a dozen experimental manipulations targeting oxidative phosphorylation (OxPhos), glycolysis, and glucocorticoid signaling, among others. These experiments provide opportunities to test how cellular energetics affect the biology of cellular aging. All data are publicly available at our webtool: https://columbia-picard.shinyapps.io/shinyapp-Lifespan_Study/.