Genome-wide evolution analysis reveals low CpG contents of fast-evolving genes and identifies antiviral microRNAs

Agarwal, 2015, Predicting effective microRNA target sites in mammalian mRNAs, eLife, 4, 10.7554/eLife.05005 Alinejad-Rokny, 2016, Source of CpG depletion in the HIV-1 genome, Mol. Biol. Evol., 33, 3205, 10.1093/molbev/msw205 Amaral, 2017, miRNA profiling of human naive CD4 T cells links miR-34c-5p to cell activation and HIV replication, EMBO J., 36, 346, 10.15252/embj.201694335 Atkinson, 2014, The influence of CpG and UpA dinucleotide frequencies on RNA virus replication and characterization of the innate cellular pathways underlying virus attenuation and enhanced replication, Nucleic Acids Res., 42, 4527, 10.1093/nar/gku075 Barreiro, 2009, From evolutionary genetics to human immunology: how selection shapes host defence genes, Nat. Rev. Genet., 11, 17, 10.1038/nrg2698 Bergelson, 1994, Decay-accelerating factor (CD55), a glycosylphosphatidylinositol-anchored complement regulatory protein, is a receptor for several echoviruses, Proc. Natl. Acad. Sci. U. S. A., 91, 6245, 10.1073/pnas.91.13.6245 Blanchette, 2004, Aligning multiple genomic sequences with the threaded blockset aligner, Genome Res., 14, 708, 10.1101/gr.1933104 Blanco-Melo, 2012, Intrinsic cellular defenses against human immunodeficiency viruses, Immunity, 37, 399, 10.1016/j.immuni.2012.08.013 Boliar, 2019, Inhibition of the lncRNA SAF drives activation of apoptotic effector caspases in HIV-1-infected human macrophages, Proc. Natl. Acad. Sci. U. S. A., 116, 7431, 10.1073/pnas.1818662116 Bruscella, 2017, Viruses and miRNAs: more friends than foes, Front. Microbiol., 8, 824, 10.3389/fmicb.2017.00824 Burns, 2009, Genetic inactivation of poliovirus infectivity by increasing the frequencies of CpG and UpA dinucleotides within and across synonymous capsid region codons, J. Virol., 83, 9957, 10.1128/JVI.00508-09 Casanova, 2013, Immunology taught by human genetics, Cold Spring Harbor Symp. Quant. Biol., 78, 157, 10.1101/sqb.2013.78.019968 Chamary, 2006, Hearing silence: non-neutral evolution at synonymous sites in mammals, Nat. Rev. Genet., 7, 98, 10.1038/nrg1770 Cheng, 2013, CpG usage in RNA viruses: data and hypotheses, PLoS One, 8, e74109, 10.1371/journal.pone.0074109 Compton, 2012, The host restriction factor APOBEC3G and retroviral Vif protein coevolve due to ongoing genetic conflict, Cell Host Microbe, 11, 91, 10.1016/j.chom.2011.11.010 Compton, 2013, Host gene evolution traces the evolutionary history of ancient primate lentiviruses, Philos. Trans. R. Soc. Lond. B Biol. Sci., 368, 20120496, 10.1098/rstb.2012.0496 Daub, 2017, Detection of pathways affected by positive selection in primate lineages ancestral to humans, Mol. Biol. Evol., 34, 1391, 10.1093/molbev/msx083 Daugherty, 2012, Rules of engagement: molecular insights from host-virus arms races, Annu. Rev. Genet., 46, 677, 10.1146/annurev-genet-110711-155522 Deaton, 2011, CpG islands and the regulation of transcription, Genes Dev., 25, 1010, 10.1101/gad.2037511 Dörig, 1993, The human CD46 molecule is a receptor for measles virus (Edmonston strain), Cell, 75, 295, 10.1016/0092-8674(93)80071-L Duggal, 2012, Evolutionary conflicts between viruses and restriction factors shape immunity, Nat. Rev. Immunol., 12, 687, 10.1038/nri3295 Dutheil, 2014, MafFilter: a highly flexible and extensible multiple genome alignment files processor, BMC Genomics, 15, 53, 10.1186/1471-2164-15-53 Fumagalli, 2011, Signatures of environmental genetic adaptation pinpoint pathogens as the main selective oressure through human evolution, PLoS Genet., 7, e1002355, 10.1371/journal.pgen.1002355 Gaggar, 2003, CD46 is a cellular receptor for group B adenoviruses, Nat. Med., 9, 1408, 10.1038/nm952 Gaunt, 2016, Elevation of CpG frequencies in influenza A genome attenuates pathogenicity but enhances host response to infection, eLife, 5, e12735, 10.7554/eLife.12735 Girardi, 2018, On the importance of host microRNAs during viral infection, Front. Genet., 9, 439, 10.3389/fgene.2018.00439 Goldman, 1994, A codon-based model of nucleotide substitution for protein-coding DNA sequences, Mol. Biol. Evol., 11, 725 Goymer, 2007, Synonymous mutations break their silence, Nat. Rev. Genet., 8, 92 Greenbaum, 2014, Quantitative theory of entropic forces acting on constrained nucleotide sequences applied to viruses, Proc. Natl. Acad. Sci. U. S. A., 111, 5054, 10.1073/pnas.1402285111 Greve, 1989, The major human rhinovirus receptor is ICAM-1, Cell, 56, 839, 10.1016/0092-8674(89)90688-0 Hamilton, 1990, Sexual reproduction as an adaptation to resist parasites (a review), Proc. Natl. Acad. Sci. U. S. A., 87, 3566, 10.1073/pnas.87.9.3566 Hodgkinson, 2011, Variation in the mutation rate across mammalian genomes, Nat. Rev. Genet., 12, 756, 10.1038/nrg3098 Hoover, 2015, Global survey of variation in a human olfactory receptor gene reveals signatures of non-neutral evolution, Chem. Senses, 40, 481, 10.1093/chemse/bjv030 Jiang, 2007, Features and trend of loss of promoter-associated CpG islands in the human and mouse genomes, Mol. Biol. Evol., 24, 1991, 10.1093/molbev/msm128 Jiang, 2018, Self-recognition of an inducible host lncRNA by RIG-I feedback restricts innate immune response, Cell, 173, 906, 10.1016/j.cell.2018.03.064 Karlin, 1994, Why is CpG suppressed in the genomes of virtually all small eukaryotic viruses but not in those of large eukaryotic viruses?, J. Virol., 68, 2889, 10.1128/JVI.68.5.2889-2897.1994 Karlin, 1997, Compositional differences within and between eukaryotic genomes, Proc. Natl. Acad. Sci. U. S. A., 94, 10227, 10.1073/pnas.94.19.10227 Karlsson, 2014, Natural selection and infectious disease in human populations, Nat. Rev. Genet., 15, 379, 10.1038/nrg3734 Karnauchow, 1996, The HeLa cell receptor for enterovirus 70 is decay-accelerating factor (CD55), J. Virol., 70, 5143, 10.1128/JVI.70.8.5143-5152.1996 Kimchi-Sarfaty, 2007, A “silent” polymorphism in the MDR1 gene changes substrate specificity, Science, 315, 525, 10.1126/science.1135308 Kimura, 2013, Stabilization of human interferon-α1 mRNA by its antisense RNA, Cell. Mol. Life Sci., 70, 1451, 10.1007/s00018-012-1216-x Kimura, 2015, Interferon-alpha competing endogenous RNA network antagonizes microRNA-1270, Cell. Mol. Life Sci., 72, 2749, 10.1007/s00018-015-1875-5 Kirchhoff, 2010, Immune evasion and counteraction of restriction factors by HIV-1 and other primate lentiviruses, Cell Host Microbe, 8, 55, 10.1016/j.chom.2010.06.004 Knox, 2018, A functional SNP in the 3′-UTR of TAP2 gene interacts with microRNA hsa-miR-1270 to suppress the gene expression, Environ. Mol. Mutagen., 59, 134, 10.1002/em.22159 Kondratowicz, 2011, T-cell immunoglobulin and mucin domain 1 (TIM-1) is a receptor for Zaire Ebolavirus and Lake Victoria Marburgvirus, Proc. Natl. Acad. Sci. U. S. A., 108, 8426, 10.1073/pnas.1019030108 Kosakovsky Pond, 2005, Not so different after all: a comparison of methods for detecting amino acid sites under selection, Mol. Biol. Evol., 22, 1208, 10.1093/molbev/msi105 Kosiol, 2008, Patterns of positive selection in six mammalian genomes, PLoS Genet., 4, e1000144, 10.1371/journal.pgen.1000144 Kotzin, 2019, The long noncoding RNA Morrbid regulates CD8 T cells in response to viral infection, Proc. Natl. Acad. Sci. U. S. A., 116, 11916, 10.1073/pnas.1819457116 Kozomara, 2018, miRBase: from microRNA sequences to function, Nucleic Acids Res., 47, D155, 10.1093/nar/gky1141 Laguette, 2012, Evolutionary and functional analyses of the interaction between the myeloid restriction factor SAMHD1 and the lentiviral Vpx protein, Cell Host Microbe, 11, 205, 10.1016/j.chom.2012.01.007 Larsen, 2004, Inhibition of infectious pancreatic necrosis virus replication by atlantic salmon Mx1 protein, J. Virol., 78, 7938, 10.1128/JVI.78.15.7938-7944.2004 Lasky, 1987, Delineation of a region of the human immunodeficiency virus type 1 gp120 glycoprotein critical for interaction with the CD4 receptor, Cell, 50, 975, 10.1016/0092-8674(87)90524-1 Li, 2014, The long noncoding RNA THRIL regulates TNFα expression through its interaction with hnRNPL, Proc. Natl. Acad. Sci. U. S. A., 111, 1002, 10.1073/pnas.1313768111 Li, 2016, Long noncoding RNA NRON contributes to HIV-1 latency by specifically inducing tat protein degradation, Nat. Commun., 7, 11730, 10.1038/ncomms11730 Li, 1993, Unbiased estimation of the rates of synonymous and nonsynonymous substitution, J. Mol. Evol., 36, 96, 10.1007/BF02407308 Li, 2017, TRIM25 enhances the antiviral action of zinc-finger antiviral protein (ZAP), PLoS Pathog., 13, e1006145, 10.1371/journal.ppat.1006145 Li, 1985, A new method for estimating synonymous and nonsynonymous rates of nucleotide substitution considering the relative likelihood of nucleotide and codon changes, Mol. Biol. Evol., 2, 150 Lim, 2012, The ability of primate lentiviruses to degrade the monocyte restriction factor SAMHD1 preceded the birth of the viral accessory protein Vpx, Cell Host Microbe, 11, 194, 10.1016/j.chom.2012.01.004 Lim, 2010, Ancient adaptive evolution of tetherin shaped the functions of Vpu and Nef in human immunodeficiency virus and primate lentiviruses, J. Virol., 84, 7124, 10.1128/JVI.00468-10 Lin, 2019, The long noncoding RNA Lnczc3h7a promotes a TRIM25-mediated RIG-I antiviral innate immune response, Nat. Immunol., 20, 812, 10.1038/s41590-019-0379-0 Lindblad-Toh, 2011, A high-resolution map of human evolutionary constraint using 29 mammals, Nature, 478, 476, 10.1038/nature10530 Liu, 2019, Proteomic profiling of HIV-1 infection of human CD4(+) T cells identifies PSGL-1 as an HIV restriction factor, Nat. Microbiol., 4, 813, 10.1038/s41564-019-0372-2 Lodge, 2017, Host microRNAs-221 and -222 inhibit HIV-1 entry in macrophages by targeting the CD4 viral receptor, Cell Rep., 21, 141, 10.1016/j.celrep.2017.09.030 Londrigan, 2011, N-linked glycosylation facilitates sialic acid-independent attachment and entry of influenza A viruses into cells expressing DC-SIGN or L-SIGN, J. Virol., 85, 2990, 10.1128/JVI.01705-10 Love, 2014, Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2, Genome Biol., 15, 550, 10.1186/s13059-014-0550-8 Ma, 2016, miRNA-mRNA interaction network in non-small cell lung cancer, Interdiscip. Sci. Comput. Life Sci., 8, 209, 10.1007/s12539-015-0117-8 Maenner, 2013, ATP-dependent roX RNA remodeling by the helicase maleless enables specific association of MSL proteins, Mol. Cell, 51, 174, 10.1016/j.molcel.2013.06.011 Makova, 2015, The effects of chromatin organization on variation in mutation rates in the genome, Nat. Rev. Genet., 16, 213, 10.1038/nrg3890 Malik, 2006, Genetic conflicts during meiosis and the evolutionary origins of centromere complexity, Behav. Ecol. Sociobiol., 34, 569 Malim, 2008, HIV-1 accessory proteins—ensuring viral survival in a hostile environment, Cell Host Microbe, 3, 388, 10.1016/j.chom.2008.04.008 Mangeat, 2003, Broad antiretroviral defence by human APOBEC3G through lethal editing of nascent reverse transcripts, Nature, 424, 99, 10.1038/nature01709 Martin, 2011, Cutadapt removes adapter sequences from high-throughput sequencing reads, EMBnet. J., 17, 10, 10.14806/ej.17.1.200 McNatt, 2009, Species-specific activity of HIV-1 Vpu and positive selection of tetherin transmembrane domain variants, PLoS Pathog., 5, e1000300, 10.1371/journal.ppat.1000300 Mitchell, 2012, Evolution-guided identification of antiviral specificity determinants in the broadly acting interferon-induced innate immunity factor MxA, Cell Host Microbe, 12, 598, 10.1016/j.chom.2012.09.005 Miyata, 1980, Molecular evolution of mRNA: a method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application, J. Mol. Evol., 16, 23, 10.1007/BF01732067 Mourtada-Maarabouni, 2008, Growth arrest in human T-cells is controlled by the non-coding RNA growth-arrest-specific transcript 5 (GAS5), J. Cell Sci., 121, 939, 10.1242/jcs.024646 Muller, 1994, Functional role of type I and type II interferons in antiviral defense, Science, 264, 1918, 10.1126/science.8009221 Muse, 1994, A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome, Mol. Biol. Evol., 11, 715 Nei, 1986, Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions, Mol. Biol. Evol., 3, 418 Neil, 2008, Tetherin inhibits retrovirus release and is antagonized by HIV-1 Vpu, Nature, 451, 425, 10.1038/nature06553 Nguyen, 2018, Potential role of microRNAs in the regulation of antiviral responses to influenza infection, Front. Immunol., 9, 1541, 10.3389/fimmu.2018.01541 Nickel, 2007, Human PAML browser: a database of positive selection on human genes using phylogenetic methods, Nucleic Acids Res., 36, D800, 10.1093/nar/gkm764 Nielsen, 2005, A scan for positively selected genes in the genomes of humans and chimpanzees, PLoS Biol., 3, e170, 10.1371/journal.pbio.0030170 Nishimura, 2009, Human P-selectin glycoprotein ligand-1 is a functional receptor for enterovirus 71, Nat. Med., 15, 794, 10.1038/nm.1961 Ogembo, 2013, Human complement receptor type 1/CD35 is an Epstein-Barr virus receptor, Cell Rep., 3, 371, 10.1016/j.celrep.2013.01.023 Pal, 2007, Coevolution with viruses drives the evolution of bacterial mutation rates, Nature, 450, 1079, 10.1038/nature06350 Paradis, 2004, APE: analyses of phylogenetics and evolution in R language, Bioinformatics, 20, 289, 10.1093/bioinformatics/btg412 Perez, 2013, Insertion of a GFP reporter gene in influenza virus, Curr. Protoc. Microbiol., 15G, 1 Pheasant, 2007, Raising the estimate of functional human sequences, Genome Res., 17, 1245, 10.1101/gr.6406307 Pohl, 2014, bwtool: a tool for bigWig files, Bioinformatics, 30, 1618, 10.1093/bioinformatics/btu056 Pollard, 2010, Detection of nonneutral substitution rates on mammalian phylogenies, Genome Res., 20, 110, 10.1101/gr.097857.109 Pollard, 2006, An RNA gene expressed during cortical development evolved rapidly in humans, Nature, 443, 167, 10.1038/nature05113 Quintana-Murci, 2019, Human immunology through the lens of evolutionary genetics, Cell, 177, 184, 10.1016/j.cell.2019.02.033 Ramani, 2018, PhastWeb: a web interface for evolutionary conservation scoring of multiple sequence alignments using phastCons and phyloP, Bioinformatics, 35, 2320, 10.1093/bioinformatics/bty966 Rothenburg, 2008, Rapid evolution of protein kinase PKR alters sensitivity to viral inhibitors, Nat. Struct. Mol. Biol., 16, 63, 10.1038/nsmb.1529 Santoro, 1999, CD46 is a cellular receptor for human herpesvirus 6, Cell, 99, 817, 10.1016/S0092-8674(00)81678-5 Sawyer, 2005, Positive selection of primate TRIM5α identifies a critical species-specific retroviral restriction domain, Proc. Natl. Acad. Sci. U. S. A., 102, 2832, 10.1073/pnas.0409853102 Saxonov, 2006, A genome-wide analysis of CpG dinucleotides in the human genome distinguishes two distinct classes of promoters, Proc. Natl. Acad. Sci. U. S. A., 103, 1412, 10.1073/pnas.0510310103 Schwefel, 2015, Molecular determinants for recognition of divergent SAMHD1 proteins by the lentiviral accessory protein Vpx, Cell Host Microbe, 17, 489, 10.1016/j.chom.2015.03.004 Short, 2006, Members of adenovirus species B utilize CD80 and CD86 as cellular attachment receptors, Virus Res., 122, 144, 10.1016/j.virusres.2006.07.009 Shultz, 2019, Immune genes are hotspots of shared positive selection across birds and mammals, eLife, 8, e41815, 10.7554/eLife.41815 Siepel, 2005, Evolutionarily conserved elements in vertebrate, insect, worm, and yeast genomes, Genome Res., 15, 1034, 10.1101/gr.3715005 Sironi, 2015, Evolutionary insights into host-pathogen interactions from mammalian sequence data, Nat. Rev. Genet., 16, 224, 10.1038/nrg3905 Staunton, 1989, A cell adhesion molecule, ICAM-1, is the major surface receptor for rhinoviruses, Cell, 56, 849, 10.1016/0092-8674(89)90689-2 Stremlau, 2004, The cytoplasmic body component TRIM5α restricts HIV-1 infection in Old World monkeys, Nature, 427, 848, 10.1038/nature02343 Takata, 2017, CG dinucleotide suppression enables antiviral defence targeting non-self RNA, Nature, 550, 124, 10.1038/nature24039 Tambyah, 2016, MicroRNA expression in blood of dengue patients, Ann. Clin. Biochem., 53, 466, 10.1177/0004563215604001 Trobaugh, 2017, MicroRNA regulation of RNA virus replication and pathogenesis, Trends Mol. Med., 23, 80, 10.1016/j.molmed.2016.11.003 Tulloch, 2014, RNA virus attenuation by codon pair deoptimisation is an artefact of increases in CpG/UpA dinucleotide frequencies, eLife, 3, e04531, 10.7554/eLife.04531 van der Lee, 2017, Genome-scale detection of positive selection in nine primates predicts human-virus evolutionary conflicts, Nucleic Acids Res., 45, 10634, 10.1093/nar/gkx704 Vitti, 2013, Detecting natural selection in genomic data, Annu. Rev. Genet., 47, 97, 10.1146/annurev-genet-111212-133526 Wang, 2010, Viral induction of the zinc finger antiviral protein is IRF3-dependent but NF-κB-independent, J. Biol. Chem., 285, 6080, 10.1074/jbc.M109.054486 Wang, 2017, An interferon-independent lncRNA promotes viral replication by modulating cellular metabolism, Science, 358, 1051, 10.1126/science.aao0409 Wang, 2014, The STAT3-binding long noncoding RNA lnc-DC controls human dendritic cell differentiation, Science, 344, 310, 10.1126/science.1251456 Wang, 2012, Changes of microRNA expression profiles in Vero cells induced by HSV-2 LAT overexpression, J. South. Med. Univ., 32, 1440 Wong, 2004, Detecting selection in noncoding regions of nucleotide sequences, Genetics, 167, 949, 10.1534/genetics.102.010959 Yan, 2012, Intrinsic antiviral immunity, Nat. Immunol., 13, 214, 10.1038/ni.2229 Yang, 1998, Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution, Mol. Biol. Evol., 15, 568, 10.1093/oxfordjournals.molbev.a025957 Yang, 2007, PAML 4: phylogenetic analysis by maximum likelihood, Mol. Biol. Evol., 24, 1586, 10.1093/molbev/msm088 Yi, 2019, MicroRNA-1270 modulates papillary thyroid cancer cell development by regulating SCAI, Biomed. Pharmacother., 109, 2357, 10.1016/j.biopha.2018.08.150 Zerbino, 2015, The ensembl regulatory build, Genome Biol., 16, 56, 10.1186/s13059-015-0621-5 Zhao, 2014, GPC5, a tumor suppressor, is regulated by miR-620 in lung adenocarcinoma, Mol Med Rep, 9, 2540, 10.3892/mmr.2014.2092 Zhao, 2007, Methylation-dependent transition rates are dependent on local sequence lengths and genomic regions, Mol. Biol. Evol., 24, 23, 10.1093/molbev/msl156 Zheng, 2017, TRIM25 is required for the antiviral activity of zinc finger antiviral protein, J. Virol., 91, 10.1128/JVI.00088-17