key: cord-0307982-psf09obg authors: King, Kelly; Larsen, Brendan B.; Gryseels, Sophie; Richet, Cécile; Kraberger, Simona; Jackson, Robert; Worobey, Michael; Harrison, Joseph S.; Varsani, Arvind; Doorslaer, Koenraad Van title: Co-evolutionary analysis suggests a role for TLR9 in papillomavirus restriction date: 2021-04-18 journal: bioRxiv DOI: 10.1101/2021.04.17.440006 sha: 7775d73a7ee82d085602beb94a7fda03a08fc9b4 doc_id: 307982 cord_uid: psf09obg Upon infection, DNA viruses can be sensed by pattern recognition receptors (PRRs) leading to the activation of type I and III interferons, aimed at blocking infection. Therefore, viruses must inhibit these signaling pathways, avoid being detected, or both. Papillomavirus virions are trafficked from early endosomes to the Golgi apparatus and wait for the onset of mitosis to complete nuclear entry. This unique subcellular trafficking strategy avoids detection by cytoplasmic PRRs, a property that may contribute to establishment of infection. However, as the capsid uncoats within acidic endosomal compartments, the viral DNA may be exposed to detection by toll-like receptor (TLR) 9. In this study we characterize two new papillomaviruses from bats and use molecular archeology to demonstrate that their genomes altered their nucleotide composition to avoid detection by TLR9, providing evidence that TLR9 acts as a PRR during papillomavirus infection. Furthermore, we demonstrate that TLR9, like other components of the innate immune system, is under evolutionary selection in bats, providing the first direct evidence for co-evolution between papillomaviruses and their hosts. reduced CpG content (Warren, Van Doorslaer, et al. 2015 ; Upadhyay and Vivekanandan 2015) . This overall dinucleotide depletion confounds the ability to demonstrate the cause of this 105 depletion. Bats serve as reservoirs for many viruses and and have served as the source of well-documented 107 cross-species transmission events of pathogens responsible for a myriad of epidemics, the most 108 notable being severe acute respiratory syndrome coronavirus 1, middle east respiratory syndrome To address this question, we determined the genomes of two novel bat papillomaviruses from 124 Tadarida brasiliensis (TbraPV2, TbraPV3). Taxonomically, bats are classified into two suborders; 125 Yinpterochiroptera and Yangochiroptera. By comparing the genomes of papillomaviruses 126 associated with bats from either suborder, we demonstrate that TLR9 target motifs are 127 significantly depleted and impact papillomavirus evolution. Furthermore, we extend existing data showing that Yangochiroptera TLR9 is evolving under diversifying selection. This argues that papillomavirus genomes are evolving in response to a host change. This is the first direct evidence 130 of PVs evolving in response to host evolutionary changes, thus providing direct evidence for co-131 evolution. Also, these data argue that TLR9 is a restriction factor for papillomavirus infection. In support of this notion, these novel bat papillomaviruses cluster with other previously described Chiropteran papillomaviruses ( Figure 1A) Figure 1A ). In addition 171 to 6 species of Chiroptera, the subtree contains 16 host species classified in 5 mammalian orders. We wanted to compare the evolutionary history of these diverse viruses to their hosts. Due to 173 intra-host divergence and niche adaptation, papillomaviruses infecting the same host can be 174 found in multiple phylogenetic tree clades. To ensure that viruses with a similar tissue tropism 175 and evolutionary history are compared, we extracted a subtree from the maximum likelihood tree 176 (Figure 2 ) (Smeele et al. 2018 ). This clade contains the newly identified TbraPV2 and TbraPV3 177 embedded within the largest monophyletic Chiroptera papillomavirus clade (red arrow in Figure 178 1A). We used a tanglegram to address our hypothesis of virus-host co-evolution (Figure 2A) To formally quantify the degree of co-speciation, we focused on two datasets. First, we used the 186 phylogenetic tree shown in Figure 2 . Because it was previously shown that members of the Lambdapapillomavirus genus co-evolve with their hosts (Rector et al. 2007 ), we also tested a 188 smaller subtree to avoid skewing the results (indicated with the red arrow in Figure 2 shown in Figure 2 . A ratio close to 1 indicates that a dinucleotide is seen in the sequence as often 211 as expected based on each sequence's nucleotide composition. Values lower than 1 suggest that 212 a dinucleotide is depleted. While the ApC, ApT, GpT, TpA, and TpC ratios are significantly lower 213 than 1 (one-sample t-test p-value < 0.001), we observed the most significant decrease in the CpG 214 dinucleotide ratio (Figure 3) . The median CpG content for these evolutionarily related viruses is 215 0.46. However, we noticed that the distribution has a long tail towards even lower values, suggesting that some viruses have a further reduced CpG content (Figure 3 ). When we plotted the CpG ratio for each virus on a phylogenetic tree ( Figure This reduction in the CpG ratio could be due to an overall lower GC content. We compared the 228 CpG ratio to total genomic GC content ( Figure 4C ). This analysis demonstrates that there is no 229 correlation between the decreased CpG ratio and the total GC content (linear regression: R 2 = 230 0.008, p-value = 0.27), and the reduced CpG ratio is not simply due to a lower GC content. The viruses infecting bats in the Yangochiroptera suborder have a reduced CpG content, raising 232 the possibility that the host species is influencing the CpG ratio and evolutionary trajectory of 233 these viruses. C5. CpG depletion alters codon usage without changing amino acid composition. The CpG dinucleotide is present in 8 codons coding for five different amino acids. Therefore, reducing CpG dinucleotides is expected to lead to a bias in codon usage or amino acid We determined codon usage tables for the non-overlapping coding sequences for each of the 247 papillomaviruses in the subtree described in Figure 2 . These codon usage tables were compared 248 using the Emboss 'codcmp' tool to calculate codon usage differences. The more diverse the 249 codon usage, the larger the differences between both tables. This analysis shows that compared 250 to each other, the Yangochiroptera papillomaviruses codon usage is more similar than when other 251 viruses are compared ( Figure 5A ). This suggests that the reduction in CpG leads to a more for episodic diversifying selection on the branches leading to the Yangochiroptera ( Figure 6B ). Finally, FEL identified 7 sites under diversifying selection within the Yangochiroptera. We mapped a subset of these residues as well as residues previously identified to be under The Arg at position 76 is much larger than the canonical His at this position, which will presumably Figure 7A ). This initial analysis indicates that ACGT, GCGT, TCGA, and TCGT are diminished in Yangochiroptera. We calculated the average tetramer ratio for each 307 group of viruses to normalize for the differences in overall CpG content between Yangochiroptera 308 and other viruses, as in Figure 4 . We compared the proportion of Yang to Yin, Yang to other, and Yin to other ( Figure 7B) However, the reason for this depletion is unclear. Of note, in mammalian genomes, CpGs are 331 rare outside of so-called CpG islands (Illingworth and Bird 2009 ). This is believed to be mainly 332 due to the observation that methylated CpGs are prone to deamination, resulting in C à T 333 mutations, leading to a depletion of CpG sites in the mammalian genomes over evolutionary time. Our data in Figure 3 show an increase in TpG and CpA. However, this increase does not appear 335 to be of the same magnitude as the dramatic reduction in CpG seen in the same dataset. The zinc-finger antiviral protein (ZAP) acts as a broad-spectrum antiviral restriction protein that Therefore, it is likely that the reduction in DNMT1 levels will lead to a loss of methylation on the 376 viral genomes destined for packaging and infection of the new tissue. Co-evolution alongside their hosts has been suggested to be an essential factor in the evolution The co-evolution theory would predict that the virus would need to adapt when the host evolves 434 a new skill. Therefore, as TLR9 evolves new functionalities, the virus would need to respond to 435 preserve the balance between virus and host. Indeed, our data suggest that as Yangochiroptera TLR9 is undergoing diversifying selection, papillomavirus genomes infecting these bats further depleted their CpG content, specifically in the context of a known TLR9 PAMP. This is the first 438 direct evidence of co-evolution between this family of viruses and their hosts. In the phylogenetic analysis, the viruses that infect Yinpterochiroptera and Yangochiroptera, 440 respectively, are not monophyletic but rather are present in three mixed clades (Figure 2 ; EsPV1, 441 EsPV3, and RfPV1; MscP2 and EhPV1; TbraPV1-3). This suggests that these three main clades 442 diverged before the ancestor of Yinpterochiroptera and Yangochiroptera split over 65 million years 443 ago. As these ancestral viruses co-evolved with the Yangochiroptera hosts, they selected for loss 444 of CpG. This occurred at least three separate times in the evolution of Yangochiroptera viruses. This strongly argues against a founder effect but in favor of recurring co-evolutionary interactions. Bats were captured in mist nets set over water sources, extracted from the nets and put in brown 485 paper bags. Bats were held in bags for 20 minutes, removed, and then measured and weighed. Individuals were identified to species in the field using metrics such as forearm length and weight. Feces and urine were collected from the bag or swabbed directly off the bat using a PurFlock 488 0.14" Ultrafine swab (Puritan, Guilford, Maine). All feces and urine samples were put into tubes containing 0.5 mL buffer consisting of 1x PBS and 50% Glycerol. These samples were held on 490 ice until returning to the lab where they were stored in a -80°C freezer. All applicable international, The tetramer content for each genome was calculated as described for the dinucleotides. To To test whether any of the tested tetramers depletions are statistically significant, we randomly 520 shuffled each viral genome. To ensure that each randomly shuffled sequence would maintain the 521 same dinucleotide ratio as the original sequence, we used the Altschul An additional subtree was extracted to minimize the impact of the genus Lamdapapillomavirus. The viral types included in this smaller dataset are underlined in Figure 2 . The amino acid composition for each sequence was calculated as described (Carugo 2008 and-whisker plots with the outliers (colored circles) identified using Tukey's method. RSCU values for each amino acid were compared using a two-way ANOVA with Tukey's posthoc (7) Omikronpapillomaviridae (5) Dyopipapillomaviridae (2) Omegapapillomaviridae (5) Deltapapillomaviridae (17) Epsilonpapillomaviridae (4) Pipapillomaviridae (8) Taupapillomaviridae (13) Xipapillomaviridae (15) Chipapillomaviridae (13) Dyozetapapillomaviridae (2) Lambdapapillomaviridae (12) Etapapillomaviridae (2) Treisepsilonpapillomaviridae (2) Zetapapillomaviridae (2) Dyorhopapillomaviridae (3) Dyoiotapapillomaviridae (3) Iotapapillomaviridae (4) Dyonupapillomaviridae (3) Dyokappapapillomaviridae (5) Rhopapillomaviridae (4) Dyoomikronpapillomaviridae (4) RtiPV2 EaPV1 Pathogen Recognition and Innate 618 Immunity HPV Entry into Cells Recognition of Double-Stranded RNA and Activation of NF-ΚB by Toll-like Receptor 624 3 Basic Local 626 Alignment Search Tool The Evolution of Flight in Bats: A Novel 629 Hypothesis A Central Region in the Minor Capsid Protein Papillomaviruses Facilitates Viral Genome Tethering and Membrane Penetration for 634 Mitotic Nuclear Entry Large Scale RNAi Reveals the 638 Requirement of Nuclear Envelope Breakdown for Nuclear Import of Human 639 Papillomaviruses PACo: A Novel 642 Procrustes Application to Cophylogenetic Analysis Novel Insights Into Immune Systems of Bats SPAdes: A New Genome Assembly 649 Algorithm and Its Applications to Single-Cell Sequencing ATM Deficiency and Oxidative Stress: A New Dimension of Defective 652 Response to DNA Damage The DsRNA Binding 655 Site of Human Toll-like Receptor 3 Classification of Papillomaviruses (PVs) Based on 659 189 PV Types and Proposal of Taxonomic Amendments Hemagglutinin Protein of Wild-Type Measles Virus 664 Activates Toll-Like Receptor 2 Signaling New Insights for Development of a Safe and Protective RSV Vaccine Envelope Glycoproteins B and H Are Necessary for TLR2 Activation in Permissive 671 Cells Trimmomatic: A Flexible Trimmer 674 for Illumina Sequence Data Viral Evasion and Subversion of Pattern-680 Recognition Receptor Signalling Bats as 'Special' Reservoirs for Emerging 683 Zoonotic Pathogens The Papillomavirus Major 686 Capsid Protein L1 The Ancient Evolutionary History of 690 Polyomaviruses Epigenetic Regulation of Human 693 Papillomavirus Transcription in the Productive Virus Life Cycle BBMap: A Fast, Accurate, Splice-Aware Aligner L2/VDNA Complex across the Limiting Membrane Requires the Onset of Mitosis Subcellular Trafficking of the Papillomavirus Genome during Initial 703 Infection: The Remarkable Abilities of Minor Capsid Protein L2 Recurrent Loss-of-Function Mutations Reveal Costs 707 to OAS1 Antiviral Activity in Primates Amino Acid Composition and Protein Dimension Kallikrein-8 Proteolytically Processes Human Papillomaviruses in the Extracellular 714 Space To Facilitate Entry into Host Cells CpG Usage in RNA Viruses: Data and Hypotheses Crystal Structure of Human Toll-Like Receptor 3 (TLR3) Ectodomain Structural RNA Has Lower Folding Energy than Random RNA of the Same 722 Establishment of 724 Expression Concepts of Papillomavirus Entry into Host 728 Cells Identification of a Role for the Trans-Golgi Network in Human 732 Papillomavirus 16 Pseudovirus Infection Innate Antiviral Responses by Means of TLR7-Mediated Recognition of 735 Single-Stranded RNA Human Papillomavirus Major Capsid Protein L1 Remains 739 Associated with the Incoming Viral Genome throughout the Entry Process Cruising the Cellular Highways: How Human Papillomavirus Travels from the Surface 743 to the Nucleus The 747 Evolution of Bat Nucleic Acid-Sensing Toll-like Receptors EBV 751 Latent Membrane Protein 1 Is a Negative Regulator of TLR9 High and Low Levels of Cottontail Rabbit 755 Papillomavirus E2 Protein Generate Opposite Effects on Gene Expression The DNA Methyltransferases Associate with HP1 and the SUV39H1 Histone 758 Methyltransferase Inhibition of Retroviral RNA Production by ZAP, a CCCH-Type Zinc Finger 761 Protein The World Goes Bats: 763 Living Longer and Tolerating Viruses Quantifying the Phylodynamic Forces Driving 767 Papillomavirus Evolution Multiple Evolutionary Mechanisms Drive Papillomavirus 771 Diversification TLR3 Deletion Limits 775 Mortality and Disease Severity Due to Phlebovirus Infection Patterns of 778 Evolution and Host Gene Mimicry in Influenza and Other RNA Viruses Cross Talk between Leishmania Donovani CpG DNA and Toll-like 783 Receptor 9: An Immunoinformatics Approach Mechanism and Function of a Newly Identified 786 CpG DNA Motif in Human Primary B Cells TLR9 Expression and Function Is 790 Abolished by the Cervical Cancer-Associated Human Papillomavirus Type 16 The Human Papillomavirus Type 16 794 E7 Oncoprotein Induces a Transcriptional Repressor Complex on the Toll-like Receptor A Metaanalysis of Bat Phylogenetics and Positive Selection 799 Based on Genomes and Transcriptomes from 18 Species TimeTree: A Public Knowledge-Base 802 of Divergence Times among Organisms Tree of Life Reveals Clock-like Speciation and Diversification Species-Specific Recognition of Single-Stranded RNA via Toll-like Receptor 7 808 and 8 Small Anti-Viral Compounds Activate Immune Cells via the TLR7 MyD88-Dependent 812 Signaling Pathway Enzymatic and Structural Insights for Substrate Specificity of a Family of Jumonji 815 Paco: Implementing Procrustean Approach to Cophylogeny in 819 R CpG Islands -'A Rough Guide Selective Evolution of Toll-like Receptors 3, 7, 8, and 9 in Bats Epigenetics of Human Papillomaviruses Subversion of the Innate 830 Immune System by a Retrovirus Human TLR7 or TLR8 Independently Confer 834 Responsiveness to the Antiviral Compound R-848 MAFFT: A Novel Method for Rapid Multiple Sequence Alignment Based on 837 Fast Fourier Transform Alignment Innate Immune Recognition of Viral Infection Methylation Patterns of Papillomavirus DNA, Its Influence on E2 Function, and 845 Implications in Viral Infection Building (Viral) Phylogenetic Trees Using 848 a Maximum Likelihood Approach Not So Different After All: A 851 Comparison of Methods for Detecting Amino Acid Sites Under Selection CpG Motifs in Bacterial 855 DNA Trigger Direct B-Cell Activation TimeTree: A 858 Resource for Timelines, Timetrees, and Divergence Times PartitionFinder 2: New Methods for Selecting Partitioned Models of Evolution for 862 Molecular and Morphological Phylogenetic Analyses Suv39h-Mediated Histone H3 Lysine 9 Methylation Directs DNA 867 Methylation to Major Satellite Repeats at Pericentric Heterochromatin Phylogenomic Analyses of Bat Subordinal Relationships 870 Based on Transcriptome Data Characterizing and Comparing Phylogenies from 873 Their Laplacian Spectrum Human Cytomegalovirus Evades ZAP Detection by Suppressing CpG Dinucleotides in 878 the Major Immediate Early 1 Gene Genome-Wide SiRNA Screen Identifies the 882 Retromer as a Cellular Entry Factor for Human Papillomavirus Efficient Comparative Phylogenetics on Large 886 Trees Corepressor-Dependent Silencing of Chromosomal Regions Encoding 889 Neuronal Genes Molecular Mechanism of RNA Recognition by Zinc-Finger Antiviral Protein The Papillomavirus E2 Proteins Recognition of Microorganisms and Activation of the Immune 897 Response Creating the CIPRES Science 899 Gateway for Inference of Large Phylogenetic Trees Evolution of Overlapping Genes Syncytial Virus Activates Innate Immunity through Toll-Like Receptor 2 Transcriptional Repression by the Methyl-CpG-Binding Protein MeCP2 Involves 910 a Histone Deacetylase Complex Structural Basis of CpG and 914 Inhibitory DNA Recognition by Toll-like Receptor 9 Bat Flight and Zoonotic Viruses The TLR3/TICAM-1 Pathway Is Mandatory for 922 Innate Immune Responses to Poliovirus Infection Structural Basis for Recognition of H3K4 Methylation 926 Status by the DNA Methyltransferase 3A ATRX-DNMT3-DNMT3L Domain UV Radiation Activates Toll-Like Receptor 9 930 Expression in Primary Human Keratinocytes, an Event Inhibited by Human 931 Papillomavirus 38 E6 and E7 Oncoproteins Downregulation of Toll-Like 935 Receptor 9 Expression by Beta Human Papillomavirus 38 and Implications for Cell Cycle 936 Ape 5.0: An Environment for Modern 939 Phylogenetics and Evolutionary Analyses in R Novel Circoviruses Detected in Feces of Sonoran Felids Zinc-Finger Antiviral Protein (ZAP) Is a Restriction Factor for 947 Replication of Modified Vaccinia Virus Ankara (MVA) in Human Cells Species-Specific Minimal Sequence Motif for 952 Oligodeoxyribonucleotides Activating Mouse TLR9 Direct 956 Binding of Retromer to Human Papillomavirus Type 16 Minor Capsid Protein L2 Mediates Endosome Exit during Viral Infection Histone 961 Modifications in Papillomavirus Virion Minichromosomes Maximum-Likelihood Trees for Large Alignments Murine Retroviruses Activate B Cells via Interaction with Toll-like Receptor 4 Ancient Papillomavirus-Host Co-Speciation in 973 Phytools: An R Package for Phylogenetic Comparative Biology (and 975 Other Things): Phytools: R Package Cleavage 978 of the Papillomavirus Minor Capsid Protein, L2, at a Furin Consensus Site Is Necessary 979 for Infection Understanding the Relationship between DNA 982 Methylation and Histone Lysine Methylation Ancient Adaptive Evolution of 986 the Primate Antiviral DNA-Editing Enzyme APOBEC3G Entry of Human Papillomavirus Type 16 990 by Actin-Dependent, Clathrin-and Lipid Raft-Independent Endocytosis 993 "DNMT1 Maintains Progenitor Function in Self-Renewing Somatic Tissue The Critical DNA Flanking Sequences of a CpG 997 but Not the 6 Base CpG Motif, Can Be Replaced with RNA 998 without Quantitative or Qualitative Changes in Toll-like Receptor 9-Mediated Activity Hyun Jin Kwun, Iris Cornet, Djamel Saidj The T Antigen Locus of Merkel Cell Polyomavirus 1002 Downregulates Human Toll-like Receptor 9 Expression Codon Usage in Yeast: 1005 Cluster Analysis Clearly Differentiates Highly and Lowly Expressed Genes Adaptive 1008 Evolution of Energy Metabolism Genes and the Origin of Flight in Bats Diverse Papillomaviruses Less Is More: An Adaptive Branch-Site Random Effects 1017 Model for Efficient Detection of Episodic Diversifying Selection Signatures of Host-Pathogen Evolutionary Conflict Reveal MISTR-A 1021 Conserved MItochondrial STress Response Network Integrated Fossil and Molecular Data Reconstruct Bat Echolocation RAxML Version 8: A Tool for Phylogenetic Analysis and Post-1028 Analysis of Large Phylogenies Dose-Dependent Regulation of the Early Promoter of Human 1031 Papillomavirus Type 18 by the Viral E2 Protein Significance of Nucleotide Sequence 1034 Alignments: A Method for Random Sequence Permutation That Preserves Dinucleotide 1035 and Codon Usage Sp100 Colocalizes with HPV Replication Foci and Restricts the Productive Stage 1039 of the Infectious Cycle Immune System Modulation and Viral 1042 Persistence in Bats: Understanding Viral Spillover CG Dinucleotide Suppression Enables 1046 Antiviral Defence Targeting Non-Self RNA Novel Bat Adenoviruses with Low G+C Content Shed 1050 New Light on the Evolution of Adenoviruses Influenza Virus Evolution, Host Adaptation, 1053 and Pandemic Formation Microbat Paraphyly and the Convergent Evolution of a Key Innovation 1057 in Old World Rhinolophoid Microbats DNA Binding and 1060 Bending by the Human Papillomavirus Type 16 E2 Protein. Recognition of an Extended 1061 Binding Site The BPV1-E2 Trans-Acting Protein Can Be Either an 1064 Activator or a Repressor of the HPV18 Regulatory Region Toll-like Receptors Regulation of Viral Infection and 1067 Disease☆ Pattern Recognition Receptors and the Innate Immune Response to Viral 1071 Infection Phylogenomic Analyses Elucidate the Evolutionary Relationships of Bats Attenuation of CGAS/STING Activity during Mitosis Vesicular Trafficking Permits Evasion of CGAS/STING 1080 Surveillance during Initial Human Papillomavirus Infection Depletion of CpG Dinucleotides in 1084 Papillomaviruses and Polyomaviruses: A Role for Divergent Evolutionary Pressures Evolution of the Papillomaviridae ICTV Virus Taxonomy Profile: 1091 Papillomaviridae The Launch of an International Animal 1094 A Major Update to the Papillomavirus Sequence Database Molecular Archeological Evidence in 1100 Support of the Repeated Loss of a Papillomavirus Gene The 1104 Papillomavirus Episteme: A Central Resource for Papillomavirus Sequence Data and 1105 Analysis Classification of Papillomaviruses Hepatitis B Virus Impairs TLR9 Expression and 1112 Function in Plasmacytoid Dendritic Cells Differential Methylation of 1115 the HPV 16 Upstream Regulatory Region during Epithelial Differentiation and 1116 Edited by Torbjorn Ramqvist Evidence for SARS-CoV-2 1120 Related Coronaviruses Circulating in Bats and Pangolins in Southeast Asia L2, the Minor Capsid Protein of 1123 Treeio: An R Package for Phylogenetic Tree Input and Output 1126 with Richly Annotated and Associated Data Role of the Host Restriction Factor APOBEC3 on Papillomavirus 1130 Evolution APOBEC3A Functions as a 1133 Restriction Factor of Human Papillomavirus Double-Stranded RNA Is Produced by Positive-Strand RNA Viruses and 1137 DNA Viruses but Not in Detectable Amounts by Negative-Strand RNA Viruses Distribution, Silencing Potential and Evolutionary 1141 Impact of Promoter DNA Methylation in the Human Genome RELAX: Detecting Relaxed Selection in a Phylogenetic Framework Origin and Evolution of Papillomavirus 1148 (Onco)Genes and Genomes Codon Usage Bias and the Evolution of Influenza A Viruses. Codon Usage Biases of 1153 Influenza Virus Dampened STING-Dependent Interferon Activation in Bats Oligodeoxyribonucleotides Rescue Mature Spleen B Cells from Spontaneous Apoptosis 1160 and Promote Cell Cycle Entry Two Methods 1163 for Mapping and Visualizing Associated Data on Phylogeny Using Ggtree Comparative Analysis of Bat Genomes 1167 Provides Insight into the Evolution of Flight and Immunity Cutting Edge: Overlapping Functions 1171 of TLR7 and TLR9 for Innate Defense against a Herpesvirus Infection test. Significance is indicated as shown in the legend.