key: cord-0298416-zgcs1ac9 authors: Sorouri, Mahsa; Chang, Tyron; Jesudhasan, Palmy; Pinkham, Chelsea; Elde, Nels C.; Hancks, Dustin C. title: MISTR: A conserved MItochondrial STress Response network revealed by signatures of evolutionary conflict date: 2020-01-26 journal: bioRxiv DOI: 10.1101/2020.01.25.919811 sha: 9f35bd4aa65968346f18e04ab17db4921606e490 doc_id: 298416 cord_uid: zgcs1ac9 Host-pathogen conflicts leave genetic signatures of variation in homologous host genes. Using these “molecular scars” as a guide, we discovered a vertebrate-specific MItochondrial STress Response circuit (MISTR). MISTR proteins are associated with electron transport chain factors and activated by stress signals such as interferon-gamma and hypoxia. Upon stress, ultraconserved miRNAs downregulate MISTR1 followed by replacement with paralogs MISTR AntiViral (MISTRAV) or MISTR Hypoxia (MISTRH), depending on the insult. While cells lacking MISTR1 are more sensitive to apoptotic triggers, cells lacking MISTRAV or expressing the poxvirus-encoded vMISTRAV exhibit resistance to the same insults. Rapid evolution signatures across primate genomes for MISTR1 and MISTRAV indicate ancient and ongoing conflicts with pathogens. MISTR proteins are also found in plants, yeasts, and an algal virus indicating ancient origins and suggesting diverse means of altering mitochondrial function under stress. The discovery of MISTR circuitry highlights the use of evolution-guided studies to reveal fundamental biological processes. Innate immunity is a critical frontline host defense mechanism in response to pathogen infection. At the onset of infections in vertebrates, a set of more than 400 genes are transcriptionally upregulated 51 by interferon and thus termed interferon-stimulated genes (ISGs). ISGs display diverse functions such 52 as activation of cell death programs and recruitment of immune cells (e.g. dendritic cells (Schneider et To investigate MISTRAV biology, we generated three A549 clonal cell lines -C15∆1, C15∆2, C15∆3 -182 with distinct indels that disrupted the MISTRAV ORF using CRISPR/CAS ( Figure 3A ). A549 cells were 183 selected because: 1) MISTRAV is interferon-inducible in these cells ( Figure 1D , Figure RT-PCR indicated that C15∆1 and C15∆2 cells lack full-length (FL) MISTRAV RNA expression in 192 IFN-g treated cells at steady-state while C15∆3 cells display a fortuitous and drastic increase of the same 193 transcript ( Figure 3C ). miRNA qPCR demonstrated that C15∆1 and C15∆2, maintain miR-147b at levels 194 comparable to wild-type with expression of miR-147b in C15∆3 ~30 fold greater than WT ( Figure 3D ). Thus, C15∆1 and C15∆2 lack MISTRAV protein but maintain the miRNA, while C15∆3 lacks MISTRAV 196 with a gain-of-function of miR-147b expression. Based on MISTRAV mitochondrial localization and numerous documented connections between 198 immune responses involving cell death mediated through mitochondria, we reasoned that MISTR might 199 mediate apoptotic responses. We primed WT and KO cells with IFN-g then added the commonly used 200 activator of apoptosis, staurosporine (STS), for 16 hours followed by functional analysis ( Figure 3E ). Assays were normalized to either untreated controls or to the number of cells being tested to account for 202 differences in proliferation rates (Figure 3-figure supplement 1) . Interestingly, we observed that C15∆1 203 and C15∆2 displayed reduced sensitivity to STS, while C15∆3 showed increased sensitivity to STS 204 compared to WT cells ( Figure 3F ). Consistent with these results, we detected robust caspase-3/7 205 cleavage activity for C15∆3 compared to WT ( Figure 3G ; p<0.0001). In addition, detectable decreases in 206 caspase-3/7 activity were observed for C15∆1 and C15∆2 relative to WT treated cells ( Figure 3G ; 207 p<0.0001). The differential sensitivities of the clones to STS relative to WT were consistent with levels of 9 PARP cleavage across the clones in response to STS ( Figure 3H ). These data suggest a role for 209 MISTRAV and miR-147b in apoptosis. We define MISTRAV as an IFN-g-inducible gene ( Figure 1D ) and protein ( Figure proposed circuit had not been tested ( Figure 5H ). Our model predicts that MISTR1 is a ubiquitously expressed sensor of stress. Specific stress 369 signals induce miRNA expression leading to the downregulation of MISTR1 and its replacement by The embedded nature of miR-147b implies a step-wise molecular progression of this response. The next day, cells were transfected with 50 ng/well of the psiCHECK-2 (Promega) construct using the Confocal Images. One day following transfection of either 1 µg hMISTRAV-GFP or 1 µg vMISTRAV-GFP, HeLa cells were fixed and stained with mitoTracker Red (Thermo). Images were taken in the confocal 533 microscopy core at the University of Utah. The location of pre-mir-147b is marked by the blue box below predicted protein-coding mRNAs. Predicting effective microRNA target sites 571 in mammalian mRNAs OASes and STING: Adaptive Evolution in Concert Detecting 688 Individual Sites Subject to Episodic Diversifying Selection Unconstrained Bayesian AppRoximation for Inferring Selection IL-10 encoded by viruses: a remarkable example of independent 696 acquisition of a cellular gene by viruses and its subsequent evolution in the viral genome A Mitochondrial Protein Compendium Elucidates Complex I 701 Disease Biology A Molecular Phylogeny of Living Primates UCSF Chimera-A visualization system for exploratory research and analysis NDUFA4 (Renamed COXFA4) Is a Cytochrome-712 c Oxidase Subunit NDUFA4 Mutations Underlie Dysfunction of a Cytochrome c 716 Oxidase Subunit Linked to Human Neurological Disease Ancient Adaptive Evolution of the Primate Antiviral 719 DNA-Editing Enzyme APOBEC3G Interferon-Stimulated Genes: A Complex Web 722 of Host Defenses Interferon-stimulated genes: roles in viral pathogenesis A diverse range of gene products are effectors of the type I interferon antiviral response Pan-viral specificity of IFN-induced genes 733 reveals new roles for cGAS in innate immunity Evolutionarily conserved elements in vertebrate, insect, 737 worm, and yeast genomes Nucleotide sequences related to the transforming gene 740 of avian sarcoma virus are present in DNA of uninfected vertebrates Induction of the Mitochondrial 745 NDUFA4L2 Protein by HIF-1α Decreases Oxygen Consumption by Inhibiting Complex I Activity SWISS-MODEL: homology modelling of protein structures 750 and complexes Modern Web Application for Characterizing Selective and Other Evolutionary Processes miRDB: an online resource for microRNA target prediction and 757 functional annotations PAML 4: Phylogenetic Analysis by Maximum Likelihood