First evidence of independent pseudogenization of toll-like receptor 5 in passerine birds

Recurrent Independent Pseudogenization Events of the Sperm Fertilization Gene ZP3r in Apes and Monkeys

Many reproductive proteins show signatures of rapid evolution through sequence divergence and duplication. These features of reproductive genes may complicate the detection of orthologs across taxa, making it difficult to connect studies in model systems to human biology. In mice, ZP3r/sp56 is a binding partner to the egg coat protein ZP3 and may mediate induction of the acrosome reaction, a crucial step in fertilization. In rodents, ZP3r, as a member of the Regulators of Complement Activation cluster, is surrounded by paralogs, some of which have been shown to be evolving under positive selection. Although primate egg coats also contain ZP3, sequence divergence paired with paralogous relationships with neighboring genes has complicated the accurate identification of the human ZP3r ortholog. Here, we phylogenetically and syntenically resolve that the human ortholog of ZP3r is the pseudogene C4BPAP1. We investigate the evolution of this gene within primates. We observe independent pseudogenization events of ZP3r in all Apes with the exception of Orangutans, and independent pseudogenization events in many monkey species. ZP3r in both primates that retain ZP3r and in rodents contains positively selected sites. We hypothesize that redundant mechanisms mediate ZP3 recognition in mammals and ZP3r's relative importance to ZP recognition varies across species.

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.

Divergent evolution drives high diversity of toll-like receptors (TLRs) in passerine birds: Buntings and finches

Toll-like receptors (TLRs) form a key component of animal innate immunity, being responsible for recognition of conserved microbial structures. As such, TLRs may be subject to diversifying and balancing selection, which maintains allelic variation both within and between populations. However, most research on TLRs in non-model avian species is focused on bottlenecked populations with depleted genetic variation. Here, we assessed variation at the extracellular domains of three TLR genes (TLR1LA, TLR3, TLR4) across eleven species from two passerine families of buntings (Emberizidae) and finches (Fringillidae), all having large breeding population sizes (millions of individuals). We found extraordinary TLR polymorphism in our study taxa, with >100 alleles detected at TLR1LA and TLR4 across species and high haplotype diversity (>0.75) in several species. Despite recent species divergence, no nucleotide allelic variants were shared between species, suggesting rapid TLR evolution. Higher variation at TLR1LA and TLR4 than TLR3 was associated with a stronger signal of diversifying selection, as measured with nucleotide substitutions rates and the number of positively selected sites (PSS). Structural protein modelling of TLRs showed that some PSS detected within TLR1LA and TLR4 were previously recognized as functionally important sites or were located in their proximity, possibly affecting ligand recognition. Furthermore, we identified PSS responsible for major surface electrostatic charge clustering, which may indicate their adaptive importance. Our study provides compelling evidence for the divergent evolution of TLR genes in buntings and finches and indicates that high TLR variation may be adaptively maintained via diversifying selection acting on functional ligand binding sites.

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.

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.

Cannabinoid receptor 2 evolutionary gene loss makes parrots more susceptible to neuroinflammation

In vertebrates, cannabinoids modulate neuroimmune interactions through two cannabinoid receptors (CNRs) conservatively expressed in the brain (CNR1, syn. CB1) and in the periphery (CNR2, syn. CB2). Our comparative genomic analysis indicates several evolutionary losses in the CNR2 gene that is involved in immune regulation. Notably, we show that the CNR2 gene pseudogenized in all parrots (Psittaciformes). This CNR2 gene loss occurred because of chromosomal rearrangements. Our positive selection analysis suggests the absence of any specific molecular adaptations in parrot CNR1 that would compensate for the CNR2 loss in the modulation of the neuroimmune interactions. Using transcriptomic data from the brains of birds with experimentally induced sterile inflammation we highlight possible functional effects of such a CNR2 gene loss. We compare the expression patterns of CNR and neuroinflammatory markers in CNR2-deficient parrots (represented by the budgerigar, Melopsittacus undulatus and five other parrot species) with CNR2-intact passerines (represented by the zebra finch, Taeniopygia guttata). Unlike in passerines, stimulation with lipopolysaccharide resulted in neuroinflammation in the parrots linked with a significant upregulation of expression in proinflammatory cytokines (including interleukin 1 beta (IL1B) and 6 (IL6)) in the brain. Our results indicate the functional importance of the CNR2 gene loss for increased sensitivity to brain inflammation.

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.

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.

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

Toll-like receptors (TLRs) play an important role for the innate immune system by detecting pathogen-associated molecular patterns. TLR5 encodes the major extracellular receptor for bacterial flagellin and frequently evolves under positive selection, consistent with coevolutionary arms races between the host and pathogens. Furthermore, TLR5 is inactivated in several vertebrates and a TLR5 stop codon polymorphism is widespread in human populations. Here, we analyzed the genomes of 120 mammals and discovered that TLR5 is convergently lost in four independent lineages, comprising guinea pigs, Yangtze river dolphin, pinnipeds, and pangolins. Validated inactivating mutations, absence of protein-coding transcript expression, and relaxed selection on the TLR5 remnants confirm these losses. PCR analysis further confirmed the loss of TLR5 in the pinniped stem lineage. Finally, we show that TLR11, encoding a second extracellular flagellin receptor, is also absent in these four lineages. Independent losses of TLR5 and TLR11 suggest that a major pathway for detecting flagellated bacteria is not essential for different mammals and predicts an impaired capacity to sense extracellular flagellin.

Evolutionary History of the Toll-Like Receptor Gene Family across Vertebrates

Adaptation to a wide range of pathogenic environments is a major aspect of the ecological adaptations of vertebrates during evolution. Toll-like receptors (TLRs) are ancient membrane-bound sensors in animals and are best known for their roles in detecting and defense against invading pathogenic microorganisms. To understand the evolutionary history of the vertebrate TLR gene family, we first traced the origin of single-cysteine cluster TLRs that share the same protein architecture with vertebrate TLRs in early-branching animals and then analyzed all members of the TLR family in over 200 species covering all major vertebrate clades. Our results indicate that although the emergence of single-cysteine cluster TLRs predates the separation of bilaterians and cnidarians, most vertebrate TLR members originated shortly after vertebrate emergence. Phylogenetic analyses divided 1,726 vertebrate TLRs into 8 subfamilies, and TLR3 may represent the most ancient subfamily that emerged before the branching of deuterostomes. Our analysis reveals that purifying selection predominated in the evolution of all vertebrate TLRs, with mean dN/dS (ω) values ranging from 0.082 for TLR21 in birds to 0.434 for TLR11 in mammals. However, we did observe patterns of positive selection acting on specific codons (527 of 60,294 codons across all vertebrate TLRs, 8.7‰), which are significantly concentrated in ligand-binding extracellular domains and suggest host–pathogen coevolutionary interactions. Additionally, we found stronger positive selection acting on nonviral compared with viral TLRs, indicating the more essential nonredundant function of viral TLRs in host immunity. Taken together, our findings provide comprehensive insight into the complex evolutionary processes of the vertebrate TLR gene family, involving gene duplication, pseudogenization, purification, and positive selection.

Avian Toll-like receptor allelic diversity far exceeds human polymorphism: an insight from domestic chicken breeds

Immune genes show remarkable levels of adaptive variation shaped by pathogen-mediated selection. Compared to humans, however, population polymorphism in animals has been understudied. To provide an insight into immunogenetic diversity in birds, we sequenced complete protein-coding regions of all Toll-like receptor (TLR) genes with direct orthology between mammals and birds (TLR3, TLR4, TLR5 and TLR7) in 110 domestic chickens from 25 breeds and compared their variability with a corresponding human dataset. Chicken TLRs (chTLRs) exhibit on average nine-times higher nucleotide diversity than human TLRs (hTLRs). Increased potentially functional non-synonymous variability is found in chTLR ligand-binding ectodomains. While we identified seven sites in chTLRs under positive selection and found evidence for convergence between alleles, no selection or convergence was detected in hTLRs. Up to six-times more alleles were identified in fowl (70 chTLR4 alleles vs. 11 hTLR4 alleles). In chTLRs, high numbers of alleles are shared between the breeds and the allelic frequencies are more equal than in hTLRs. These differences may have an important impact on infectious disease resistance and host-parasite co-evolution. Though adaptation through high genetic variation is typical for acquired immunity (e.g. MHC), our results show striking levels of intraspecific polymorphism also in poultry innate immune receptors.

Toll-Like Receptor Evolution in Birds: Gene Duplication, Pseudogenization, and Diversifying Selection

Toll-like receptors (TLRs) are key sensor molecules in vertebrates triggering initial phases of immune responses to pathogens. The avian TLR family typically consists of ten receptors, each adapted to distinct ligands. To understand the complex evolutionary history of each avian TLR, we analyzed all members of the TLR family in the whole genome assemblies and target sequence data of 63 bird species covering all major avian clades. Our results indicate that gene duplication events most probably occurred in TLR1 before synapsids diversified from sauropsids. Unlike mammals, ssRNA-recognizing TLR7 has duplicated independently in several avian taxa, while flagellin-sensing TLR5 has pseudogenized multiple times in bird phylogeny. Our analysis revealed stronger positive, diversifying selection acting in TLR5 and the three-domain TLRs (TLR10 [TLR1A], TLR1 [TLR1B], TLR2A, TLR2B, TLR4) that face the extracellular space and bind complex ligands than in single-domain TLR15 and endosomal TLRs (TLR3, TLR7, TLR21). In total, 84 out of 306 positively selected sites were predicted to harbor substitutions dramatically changing the amino acid physicochemical properties. Furthermore, 105 positively selected sites were located in the known functionally relevant TLR regions. We found evidence for convergent evolution acting between birds and mammals at 54 of these sites. Our comparative study provides a comprehensive insight into the evolution of avian TLR genetic variability. Besides describing the history of avian TLR gene gain and gene loss, we also identified candidate positions in the receptors that have been likely shaped by direct molecular host-pathogen coevolutionary interactions and most probably play key functional roles in birds.

A High-Quality, Long-Read De Novo Genome Assembly to Aid Conservation of Hawaii's Last Remaining Crow Species

Genome-level data can provide researchers with unprecedented precision to examine the causes and genetic consequences of population declines, which can inform conservation management. Here, we present a high-quality, long-read, de novo genome assembly for one of the world’s most endangered bird species, the ʻAlalā (Corvus hawaiiensis; Hawaiian crow). As the only remaining native crow species in Hawaiʻi, the ʻAlalā survived solely in a captive-breeding program from 2002 until 2016, at which point a long-term reintroduction program was initiated. The high-quality genome assembly was generated to lay the foundation for both comparative genomics studies and the development of population-level genomic tools that will aid conservation and recovery efforts. We illustrate how the quality of this assembly places it amongst the very best avian genomes assembled to date, comparable to intensively studied model systems. We describe the genome architecture in terms of repetitive elements and runs of homozygosity, and we show that compared with more outbred species, the ʻAlalā genome is substantially more homozygous. We also provide annotations for a subset of immunity genes that are likely to be important in conservation management, and we discuss how this genome is currently being used as a roadmap for downstream conservation applications.

Functional expression of TLR5 of different vertebrate species and diversification in intestinal pathogen recognition

Toll-like receptor 5 (TLR5) is activated by bacterial flagellins and plays a crucial role in the first-line defence against pathogenic bacteria and in immune homeostasis, and is highly conserved in vertebrate species. However, little comparative information is available on TLR5 functionality. In this study, we compared TLR5 activation using full-length and chimeric TLR5 of various vertebrate species (human, chicken, mouse, pig, cattle). Chimeric TLR5 receptors, consisting of human transmembrane and intracellular domains, linked to extracellular domains of animal origin, were generated and expressed. The comparison of chimeric TLR5s and their full-length counterparts revealed significant functional disparities. While porcine and chicken full-length TLR5s showed a strongly reduced functionality in human cells, all chimeric receptors were functional when challenged with TLR5 ligand Salmonella FliC. Using chimeric receptors as a tool allowed for the identification of ectodomain-dependent activation potential and partially host species-specific differences in response to various enteric bacterial strains and their purified flagellins. We conclude that both the extra- and intracellular determinants of TLR5 receptors are crucial for compatibility with the species expression background and hence for proper receptor functionality. TLR5 receptors with a common intracellular domain provide a useful system to investigate bacteria- and host-specific differences in receptor activation.

Purifying selection and concerted evolution of RNA-sensing toll-like receptors in migratory waders

Migratory birds encounter a broad range of pathogens during their journeys, making them ideal models for studying immune gene evolution. Despite the potential value of these species to immunoecology and disease epidemiology, previous studies have typically focused on their adaptive immune gene repertoires. In this study, we examined the evolution of innate immune genes in three long-distance migratory waders (order Charadriiformes). We analysed two parts of the extracellular domains of two Toll-like receptors (TLR3 and TLR7) involved in virus recognition in the Sanderling (Calidris alba), Red-necked Stint (Calidris ruficollis), and Ruddy Turnstone (Arenaria interpres). Our analysis was extended to 50 avian species for which whole-genome sequences were available, including two additional waders. We found that the inferred relationships among avian TLR3 and TLR7 do not match the whole-genome phylogeny of birds. Further analyses showed that although both loci are predominantly under purifying selection, the evolution of the extracellular domain of avian TLR3 has also been driven by episodic diversifying selection. TLR7 was found to be duplicated in all five wader species and in two other orders of birds, Cuculiformes and Passeriformes. The duplication is likely to have occurred in the ancestor of each order, and the duplicated copies appear to be undergoing concerted evolution. The phylogenetic relationships of wader TLR7 matched those of the five wader species, but that of TLR3 did not. Instead, the tree inferred from TLR3 showed potential associations with the species' ecology, including migratory behaviour and exposure to pathogens. Our study demonstrates the importance of combining immunological and ecological knowledge to understand the impact of immune gene polymorphism on the evolutionary ecology of infectious diseases.

Characterisation of Toll-like receptors 4, 5 and 7 and their genetic variation in the grey partridge

Toll-like receptors (TLRs) are a cornerstone of vertebrate innate immunity. In this study, we identified orthologues of TLR4, TLR5 and TLR7 (representing both bacterial- and viral-sensing TLRs) in the grey partridge (Perdix perdix), a European Galliform game bird species. The phylogeny of all three TLR genes follows the known phylogeny of Galloanserae birds, placing grey partridge TLRs (PePeTLRs) in close proximity to their turkey and pheasant orthologues. The predicted proteins encoded by the PePeTLR genes were 843, 862-863 and 1,047 amino acids long, respectively, and clearly showed all TLR structural features. To verify functionality in these genes we mapped their tissue-expression profiles, revealing generally high PePeTLR4 and PePeTLR5 expression in the thymus and absence of PePeTLR4 and PePeTLR7 expression in the brain. Using 454 next-generation sequencing, we then assessed genetic variation within these genes for a wild grey partridge population in the Czech Republic, EU. We identified 11 nucleotide substitutions in PePeTLR4, eight in PePeTLR5 and six in PePeTLR7, resulting in four, four and three amino acid replacements, respectively. Given their locations and chemical features, most of these non-synonymous substitutions probably have a minor functional impact. As the intraspecific genetic variation of the three TLR genes was low, we assume that either negative selection or a bottleneck may have reduced TLR population variability in this species.

Protein evolution of Toll-like receptors 4, 5 and 7 within Galloanserae birds

BACKGROUND

Toll-like receptors (TLR) are essential activators of the innate part of the vertebrate immune system. In this study, we analysed the interspecific variability of three TLR (bacterial-sensing TLR4 and TLR5 and viral-sensing TLR7) within the Galloanserae bird clade, investigated their phylogeny, assessed their structural conservation and estimated site-specific selection pressures.

RESULTS

Physiochemical properties varied according to the TLR analysed, mainly with regards to the surface electrostatic potential distribution. The predicted ligand-binding features (mainly in TLR4 and TLR5) differed between the avian proteins and their fish and mammalian counterparts, but also varied within the Galloanserae birds. We identified 20 positively selected sites in the three TLR, among which several are topologically close to ligand-binding sites reported for mammalian and fish TLR. We described 26, 28 and 25 evolutionarily non-conservative sites in TLR4, TLR5 and TLR7, respectively. Thirteen of these sites in TLR4, and ten in TLR5 were located in functionally relevant regions. The variability appears to be functionally more conserved for viral-sensing TLR7 than for the bacterial-sensing TLR. Amino-acid positions 268, 270, 343, 383, 444 and 471 in TLR4 and 180, 183, 209, 216, 264, 342 and 379 in TLR5 are key candidates for further functional research.

CONCLUSIONS

Host-pathogen co-evolution has a major effect on the features of host immune receptors. Our results suggest that avian and mammalian TLR may be differentially adapted to pathogen-derived ligand recognition. We have detected signatures of positive selection even within the Galloanserae lineage. To our knowledge, this is the first study to depict evolutionary pressures on Galloanserae TLR and to estimate the validity of current knowledge on TLR function (based on mammalian and chicken models) for non-model species of this clade.