Convergent Losses of TLR5 Suggest Altered Extracellular Flagellin Detection in Four Mammalian Lineages

Natural selection directing molecular evolution in vertebrate viral sensors

Diseases caused by pathogens contribute to molecular adaptations in host immunity. Variety of viral pathogens challenging animal immunity can drive positive selection diversifying receptors recognising the infections. However, whether distinct virus sensing systems differ across animals in their evolutionary modes remains unclear. Our review provides a comparative overview of natural selection shaping molecular evolution in vertebrate viral-binding pattern recognition receptors (PRRs). Despite prevailing negative selection arising from the functional constraints, multiple lines of evidence now suggest diversifying selection in the Toll-like receptors (TLRs), NOD-like receptors (NLRs), RIG-I-like receptors (RLRs) and oligoadenylate synthetases (OASs). In several cases, location of the positively selected sites in the ligand-binding regions suggests effects on viral detection although experimental support is lacking. Unfortunately, in most other PRR families including the AIM2-like receptor family, C-type lectin receptors (CLRs), and cyclic GMP-AMP synthetase studies characterising their molecular evolution are rare, preventing comparative insight. We indicate shared characteristics of the viral sensor evolution and highlight priorities for future research.

An RNA-Seq analysis of coronavirus in the skin of the Pangolin

Protection of the Critically Endangered East Asian Pangolin species is hampered by the vulnerability of captive individuals to infection. Studies have previously shown the pangolin to have a unique pseudogenisation of many immunity genes (including IFNE, IFIH1, cGAS, STING, TLR5, and TLR11), and we suspected that these losses could account for this vulnerability. Here we used RNA-Seq data to show the effect of these gene losses on the transcriptional response to a viral skin infection in a deceased pangolin. This virus is very closely related to the one causing the current COVID-19 pandemic in the human population (SARS-CoV2), and we found the most upregulated pathway was the same one previously identified in the lungs of SARS-CoV2-infected humans. As predicted, we found that the pathways downstream of the lost genes were not upregulated. For example, the pseudogenised interferon epsilon (IFNE) is known to be particularly important in epithelial immunity, and we show that interferon-related responses were not upregulated in the infected pangolin skin. We suggest that the pangolin's innate gene pseudogenisation is indeed likely to be responsible for the animal's vulnerability to infection.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Dynamic Evolution of Avian RNA Virus Sensors: Repeated Loss of RIG-I and RIPLET

Retinoic acid-inducible gene I (RIG-I) and melanoma differentiation-associated protein 5 (MDA5) are key RNA virus sensors belonging to the RIG-I-like receptor (RLR) family. The activation of the RLR inflammasome leads to the establishment of antiviral state, mainly through interferon-mediated signaling. The evolutionary dynamics of RLRs has been studied mainly in mammals, where rare cases of RLR gene losses were described. By in silico screening of avian genomes, we previously described two independent disruptions of MDA5 in two bird orders. Here, we extend this analysis to approximately 150 avian genomes and report 16 independent evolutionary events of RIG-I inactivation. Interestingly, in almost all cases, these inactivations are coupled with genetic disruptions of RIPLET/RNF135, an ubiquitin ligase RIG-I regulator. Complete absence of any detectable RIG-I sequences is unique to several galliform species, including the domestic chicken (Gallus gallus). We further aimed to determine compensatory evolution of MDA5 in RIG-I-deficient species. While we were unable to show any specific global pattern of adaptive evolution in RIG-I-deficient species, in galliforms, the analyses of positive selection and surface charge distribution support the hypothesis of some compensatory evolution in MDA5 after RIG-I loss. This work highlights the dynamic nature of evolution in bird RNA virus sensors.

Host control and the evolution of cooperation in host microbiomes

Humans, and many other species, are host to diverse symbionts. It is often suggested that the mutual benefits of host-microbe relationships can alone explain cooperative evolution. Here, we evaluate this hypothesis with evolutionary modelling. Our model predicts that mutual benefits are insufficient to drive cooperation in systems like the human microbiome, because of competition between symbionts. However, cooperation can emerge if hosts can exert control over symbionts, so long as there are constraints that limit symbiont counter evolution. We test our model with genomic data of two bacterial traits monitored by animal immune systems. In both cases, bacteria have evolved as predicted under host control, tending to lose flagella and maintain butyrate production when host-associated. Moreover, an analysis of bacteria that retain flagella supports the evolution of host control, via toll-like receptor 5, which limits symbiont counter evolution. Our work puts host control mechanisms, including the immune system, at the centre of microbiome evolution.

Vision-related convergent gene losses reveal SERPINE3's unknown role in the eye

Despite decades of research, knowledge about the genes that are important for development and function of the mammalian eye and are involved in human eye disorders remains incomplete. During mammalian evolution, mammals that naturally exhibit poor vision or regressive eye phenotypes have independently lost many eye-related genes. This provides an opportunity to predict novel eye-related genes based on specific evolutionary gene loss signatures. Building on these observations, we performed a genome-wide screen across 49 mammals for functionally uncharacterized genes that are preferentially lost in species exhibiting lower visual acuity values. The screen uncovered several genes, including SERPINE3, a putative serine proteinase inhibitor. A detailed investigation of 381 additional mammals revealed that SERPINE3 is independently lost in 18 lineages that typically do not primarily rely on vision, predicting a vision-related function for this gene. To test this, we show that SERPINE3 has the highest expression in eyes of zebrafish and mouse. In the zebrafish retina, serpine3 is expressed in Müller glia cells, a cell type essential for survival and maintenance of the retina. A CRISPR-mediated knockout of serpine3 in zebrafish resulted in alterations in eye shape and defects in retinal layering. Furthermore, two human polymorphisms that are in linkage with SERPINE3 are associated with eye-related traits. Together, these results suggest that SERPINE3 has a role in vertebrate eyes. More generally, by integrating comparative genomics with experiments in model organisms, we show that screens for specific phenotype-associated gene signatures can predict functions of uncharacterized genes.

Parallel decay of vision genes in subterranean water beetles

In the framework of neutral theory of molecular evolution, genes specific to the development and function of eyes in subterranean animals living in permanent darkness are expected to evolve by relaxed selection, ultimately becoming pseudogenes. However, definitive empirical evidence for the role of neutral processes in the loss of vision over evolutionary time remains controversial. In previous studies, we characterized an assemblage of independently-evolved water beetle (Dytiscidae) species from a subterranean archipelago in Western Australia, where parallel vision and eye loss have occurred. Using a combination of transcriptomics and exon capture, we present evidence of parallel coding sequence decay, resulting from the accumulation of frameshift mutations and premature stop codons, in eight phototransduction genes (arrestins, opsins, ninaC and transient receptor potential channel genes) in 32 subterranean species in contrast to surface species, where these genes have open reading frames. Our results provide strong evidence to support neutral evolutionary processes as a major contributing factor to the loss of phototransduction genes in subterranean animals, with the ultimate fate being the irreversible loss of a light detection system.

Decay of TRPV3 as the genomic trace of epidermal structure changes in the land-to-sea transition of mammals

The epidermis plays an indispensable barrier function in animals. Some species have evolved unique epidermal structures to adapt to different environments. Aquatic and semi-aquatic mammals (cetaceans, manatees, and hippopotamus) are good models to study the evolution of epidermal structures because of their exceptionally thickened stratum spinosum, the lack of stratum granulosum, and the parakeratotic stratum corneum. This study aimed to analyze an upstream regulatory gene transient receptor potential cation channel, subfamily V, member 3 (TRPV3) of epidermal differentiation so as to explore the association between TRPV3 evolution and epidermal changes in mammals. Inactivating mutations were detected in almost all the aquatic cetaceans and several terrestrial mammals. Relaxed selective pressure was examined in the cetacean lineages with inactivated TRPV3, which might contribute to its exceptionally thickened stratum spinosum as the significant thickening of stratum spinosum in TRPV3 knock-out mouse. However, functional TRPV3 may exist in several terrestrial mammals due to their strong purifying selection, although they have "inactivating mutations." Further, for intact sequences, relaxed selective constraints on the TRPV3 gene were also detected in aquatic cetaceans, manatees, and semi-aquatic hippopotamus. However, they had intact TRPV3, suggesting that the accumulation of inactivating mutations might have lagged behind the relaxed selective pressure. The results of this study revealed the decay of TRPV3 being the genomic trace of epidermal development in aquatic and semi-aquatic mammals. They provided insights into convergently evolutionary changes of epidermal structures during the transition from the terrestrial to the aquatic environment.

Adaptation and cryptic pseudogenization in penguin Toll-like Receptors

Penguins (Sphenisciformes) are an iconic order of flightless, diving seabirds distributed across a large latitudinal range in the Southern Hemisphere. The extensive area over which penguins are endemic is likely to have fostered variation in pathogen pressure, which in turn will have imposed differential selective pressures on the penguin immune system. At the front line of pathogen detection and response, the Toll-like receptors (TLRs) provide insight into host evolution in the face of microbial challenge. TLRs respond to conserved pathogen-associated molecular patterns and are frequently found to be under positive selection, despite retaining specificity for defined agonist classes. We undertook a comparative immunogenetics analysis of TLRs for all penguin species and found evidence of adaptive evolution that was largely restricted to the cell surface-expressed TLRs, with evidence of positive selection at, or near, key agonist-binding sites in TLR1B, TLR4, and TLR5. Intriguingly, TLR15, which is activated by fungal products, appeared to have been pseudogenized multiple times in the Eudyptes spp., but a full-length form was present as a rare haplotype at the population level. However, in vitro analysis revealed that even the full-length form of Eudyptes TLR15 was nonfunctional, indicating an ancestral cryptic pseudogenization prior to its eventual disruption multiple times in the Eudyptes lineage. This unusual pseudogenization event could provide an insight into immune adaptation to fungal pathogens such as Aspergillus, which is responsible for significant mortality in wild and captive bird populations.

Reconstruction of evolutionary changes in fat and toxin consumption reveals associations with gene losses in mammals: a case study for the lipase inhibitor PNLIPRP1 and the xenobiotic receptor NR1I3

The inactivation of ancestral protein-coding genes (gene loss) can be associated with phenotypic modifications. Within placental mammals, repeated losses of PNLIPRP1 (gene inhibiting fat digestion) occurred preferentially in strictly herbivorous species, while repeated NR1I3 losses (gene involved in detoxification) occurred preferentially in strictly carnivorous species. It was hypothesized that lower fat contents of herbivorous diets and lower toxin contents of carnivorous diets cause relaxed selection pressure on these genes resulting in the accumulation of mutations and ultimately to convergent gene losses. However, since herbivorous and carnivorous diets differ vastly in their composition, a fine-grained analysis is required for hypothesis testing. We generated a trait matrix recording diet and semi-quantitative estimates of fat and toxin consumption for 52 placental species. By including data from 31 fossil taxa, we reconstructed the ancestral diets in major lineages (grundplan reconstruction). We found support that PNLIPRP1 loss is primarily associated with low levels of fat intake and not simply with herbivory/carnivory. In particular, PNLIPRP1 loss also occurred in carnivorous lineages feeding on a fat-poor diet, suggesting that the loss of this gene may be beneficial for occupying ecological niches characterized by fat-poor food resources. Similarly, we demonstrated that carnivorous species are indeed less exposed to diet-related toxins suggesting that the loss of NR1I3 and related genes (NR1I2, UGT1A6) resulted from relaxed selection pressure. This study illustrates the need of detailed phenotype studies to obtain a deeper understanding of factors underlying gene losses and to progress in understanding genomic causes of phenotypic variation in mammals.

NOD2 and reproduction-associated NOD-like receptors have been lost during the evolution of pangolins

NOD-like receptors (NLRs) are sensors of pathogen-associated molecular patterns with critical roles in the control of immune responses and programmed cell death. Recent studies have revealed inter-species differences in mammalian innate immune genes and a particular degeneration of nucleic acid sensing pathways in pangolins, which are currently investigated as potential hosts for zoonotic pathogens. Here, we used comparative genomics to determine which NLR genes are conserved or lost in pangolins and related mammals. We show that NOD2, which is implicated in sensing bacterial muramyl dipeptide and viral RNA, is a pseudogene in pangolins, but not in any other mammalian species investigated. NLRC4 and NAIP are absent in pangolins and canine carnivorans, suggesting convergent loss of cytoplasmic sensing of bacterial flagellin in these taxa. Among NLR family pyrin domain containing proteins (NLRPs), skin barrier-related NLRP10 has been lost in pangolins after the evolutionary divergence from Carnivora. Strikingly, pangolins lack all NLRPs associated with reproduction (germ cells and embryonic development) in other mammals, i.e., NLRP2, 4, 5, 7, 8, 9, 11, 13, and 14. Taken together, our study shows a massive degeneration of NLR genes in pangolins and suggests that these endangered mammals may have unique adaptations of innate immunity and reproductive cell biology.

Repeated MDA5 Gene Loss in Birds: An Evolutionary Perspective

Two key cytosolic receptors belonging to the retinoic acid-inducible gene I (RIG-I)-like receptor (RLR) family sense the viral RNA-derived danger signals: RIG-I and melanoma differentiation-associated protein 5 (MDA5). Their activation establishes an antiviral state by downstream signaling that ultimately activates interferon-stimulated genes (ISGs). While in rare cases RIG-I gene loss has been detected in mammalian and avian species, most notably in the chicken, MDA5 pseudogenization has only been detected once in mammals. We have screened over a hundred publicly available avian genome sequences and describe an independent disruption of MDA5 in two unrelated avian lineages, the storks (Ciconiiformes) and the rallids (Gruiformes). The results of our RELAX analysis confirmed the absence of negative selection in the MDA5 pseudogene. In contrast to our prediction, we have shown, using multiple dN/dS-based approaches, that the MDA5 loss does not appear to have resulted in any compensatory evolution in the RIG-I gene, which may partially share its ligand-binding specificity. Together, our results indicate that the MDA5 pseudogenization may have important functional effects on immune responsiveness in these two avian clades.

Cytosolic DNA sensing through cGAS and STING is inactivated by gene mutations in pangolins

The release of DNA into the cytoplasm upon damage to the nucleus or during viral infection triggers an interferon-mediated defense response, inflammation and cell death. In human cells cytoplasmic DNA is sensed by cyclic GMP-AMP Synthase (cGAS) and Absent In Melanoma 2 (AIM2). Here, we report the identification of a “natural knockout” model of cGAS. Comparative genomics of phylogenetically diverse mammalian species showed that cGAS and its interaction partner Stimulator of Interferon Genes (STING) have been inactivated by mutations in the Malayan pangolin whereas other mammals retained intact copies of these genes. The coding sequences of CGAS and STING1 are also disrupted by premature stop codons and frame-shift mutations in Chinese and tree pangolins, suggesting that expression of these genes was lost in a common ancestor of all pangolins that lived more than 20 million years ago. AIM2 is retained in a functional form in pangolins whereas it is inactivated by mutations in carnivorans, the phylogenetic sister group of pangolins. The deficiency of cGAS and STING points to the existence of alternative mechanisms of controlling cytoplasmic DNA-associated cell damage and viral infections in pangolins.

Genomic and anatomical comparisons of skin support independent adaptation to life in water by cetaceans and hippos

The macroevolutionary transition from terra firma to obligatory inhabitance of the marine hydrosphere has occurred twice in the history of Mammalia: Cetacea and Sirenia. In the case of Cetacea (whales, dolphins, and porpoises), molecular phylogenies provide unambiguous evidence that fully aquatic cetaceans and semiaquatic hippopotamids (hippos) are each other's closest living relatives. Ancestral reconstructions suggest that some adaptations to the aquatic realm evolved in the common ancestor of Cetancodonta (Cetacea + Hippopotamidae). An alternative hypothesis is that these adaptations evolved independently in cetaceans and hippos. Here, we focus on the integumentary system and evaluate these hypotheses by integrating new histological data for cetaceans and hippos, the first genome-scale data for pygmy hippopotamus, and comprehensive genomic screens and molecular evolutionary analyses for protein-coding genes that have been inactivated in hippos and cetaceans. We identified eight skin-related genes that are inactivated in both cetaceans and hippos, including genes that are related to sebaceous glands, hair follicles, and epidermal differentiation. However, none of these genes exhibit inactivating mutations that are shared by cetaceans and hippos. Mean dates for the inactivation of skin genes in these two clades serve as proxies for phenotypic changes and suggest that hair reduction/loss, the loss of sebaceous glands, and changes to the keratinization program occurred ∼16 Ma earlier in cetaceans (∼46.5 Ma) than in hippos (∼30.5 Ma). These results, together with histological differences in the integument and prior analyses of oxygen isotopes from stem hippopotamids ("anthracotheres"), support the hypothesis that aquatic skin adaptations evolved independently in hippos and cetaceans.

Insights into structure and dynamics of extracellular domain of Toll-like receptor 5 in Cirrhinus mrigala (mrigala): A molecular dynamics simulation approach

The toll-like receptor 5 (TLR5) is the most conserved important pattern recognition receptors (PRRs) often stimulated by bacterial flagellins and plays a major role in the first-line defense against invading pathogenic bacteria and in immune homeostasis. Experimental crystallographic studies have shown that the extracellular domain (ECD) of TLR5 recognizes flagellin of bacteria and functions as a homodimer in model organism zebrafish. However, no structural information is available on TLR5 functionality in the major carp Cirrhinus mrigala (mrigala) and its interaction with bacterial flagellins. Therefore, the present study was undertaken to unravel the structural basis of TLR5-flagellin recognition in mrigala using structural homodimeric TLR5-flagellin complex of zebrafish as reference. Integrative structural modeling and molecular dynamics simulations were employed to explore the structural and mechanistic details of TLR5 recognition. Results from structural snapshots of MD simulation revealed that TLR5 consistently formed close interactions with the three helices of the D1 domain in flagellin on its lateral side mediated by several conserved amino acids. Results from the intermolecular contact analysis perfectly substantiate with the findings of per residue-free energy decomposition analysis. The differential recognition mediated by flagellin to TLR5 in mrigala involves charged residues at the interface of binding as compared to the zebrafish complex. Overall our results shows TLR5 of mrigala involved in innate immunity specifically recognized a conserved site on flagellin which advocates the scientific community to explore host-specific differences in receptor activation.

Asymptomatic SARS-CoV-2 infection: is it all about being refractile to innate immune sensing of viral spare-parts?-Clues from exotic animal reservoirs

A vast proportion of coronavirus disease 2019 (COVID-19) individuals remain asymptomatic and can shed severe acute respiratory syndrome (SARS-CoV) type 2 virus to transmit the infection, which also explains the exponential increase in the number of COVID-19 cases globally. Furthermore, the rate of recovery from clinical COVID-19 in certain pockets of the globe is surprisingly high. Based on published reports and available literature, here, we speculated a few immunovirological mechanisms as to why a vast majority of individuals remain asymptomatic similar to exotic animal (bats and pangolins) reservoirs that remain refractile to disease development despite carrying a huge load of diverse insidious viral species, and whether such evolutionary advantage would unveil therapeutic strategies against COVID-19 infection in humans. Understanding the unique mechanisms that exotic animal species employ to achieve viral control, as well as inflammatory regulation, appears to hold key clues to the development of therapeutic versatility against COVID-19.

Pangolins Lack IFIH1/MDA5, a Cytoplasmic RNA Sensor That Initiates Innate Immune Defense Upon Coronavirus Infection

Zoonotic infections are an imminent threat to human health. Pangolins were recently identified as carriers and intermediate hosts of coronaviruses. Previous research has shown that infection with coronaviruses activates an innate immune response upon sensing of viral RNA by interferon-induced with helicase C domain 1 (IFIH1), also known as MDA5. Here, we performed a comparative genomics study of RNA sensor genes in three species of pangolins. DDX58/RIG-I, a sensor of cytoplasmic viral RNA and toll-like receptors (TLR) 3, 7, and 8, which bind RNA in endosomes, are conserved in pangolins. By contrast, IFIH1 a sensor of intracellular double-stranded RNA, has been inactivated by mutations in pangolins. Likewise, Z-DNA-binding protein (ZBP1), which senses both Z-DNA and Z-RNA, has been lost during the evolution of pangolins. These results suggest that the innate immune response to viruses differs significantly between pangolins and other mammals, including humans. We put forward the hypothesis that loss of IFIH1 and ZBP1 provided an evolutionary advantage by reducing inflammation-induced damage to host tissues and thereby contributed to a switch from resistance to tolerance of viral infections in pangolins.