Signatures of diversifying selection and convergence acting on passerine Toll‐like receptor 4 in an evolutionary context

Positive selection acting on Toll‐like receptors (TLRs) has been recently investigated to reveal evolutionary mechanisms of host–pathogen molecular co‐adaptation. Much of this research, however, has focused mainly on the identification of sites predicted to be under positive selection, bringing little insight into the functional differences and similarities among species and a limited understanding of convergent evolution in the innate immune molecules. In this study, we provide evidence of phenotypic variability in the avian TLR4 ligand‐binding region (LBR), the direct interface between host and pathogen molecular structures. We show that 55 passerine species vary substantially in the distribution of electrostatic potential on the surface of the receptor, and based on these distinct patterns, we identified four species clusters. Seven of the 34 evolutionarily nonconservative and positively selected residues correspond topologically to sites previously identified as being important for lipopolysaccharide, lipid IVa or MD‐2 binding. Five of these positions codetermine the identity of the charge clusters. Groups of species that host‐related communities of pathogens were predicted to cluster based on their TLR4 LBR charge. Despite some evidence for convergence among taxa, there were no clear associations between the TLR4 LBR charge distribution and any of the general ecological characteristics compared (migration, latitudinal distribution and diet). Closely related species, however, mostly belonged to the same surface charge cluster indicating that phylogenetic constraints are key determinants shaping TLR4 adaptive evolution. Our results suggest that host innate immune evolution is consistent with Fahrenholz's rule on the cospeciation of hosts and their parasites.     

The evolution of bat nucleic acid‐sensing Toll‐like receptors

We characterized the nucleic acid‐sensing Toll‐like receptors (TLR) of a New World bat species, the common vampire bat (Desmodus rotundus), and through a comparative molecular evolutionary approach searched for general adaptation patterns among the nucleic acid‐sensing TLRs of eight different bats species belonging to three families (Pteropodidae, Vespertilionidae and Phyllostomidae). We found that the bat TLRs are evolving slowly and mostly under purifying selection and that the divergence pattern of such receptors is overall congruent with the species tree, consistent with the evolution of many other mammalian nuclear genes. However, the chiropteran TLRs exhibited unique mutations fixed in ligand‐binding sites, some of which involved nonconservative amino acid changes and/or targets of positive selection. Such changes could potentially modify protein function and ligand‐binding properties, as some changes were predicted to alter nucleic acid binding motifs in TLR 9. Moreover, evidence for episodic diversifying selection acting specifically upon the bat lineage and sublineages was detected. Thus, the long‐term adaptation of chiropterans to a wide variety of environments and ecological niches with different pathogen profiles is likely to have shaped the evolution of the bat TLRs in an order‐specific manner. The observed evolutionary patterns provide evidence for potential functional differences between bat and other mammalian TLRs in terms of resistance to specific pathogens or recognition of nucleic acids in general.     

Positive selection drives adaptive diversification of the 4‐coumarate: CoA ligase (4CL) gene in angiosperms

Lignin and flavonoids play a vital role in the adaption of plants to a terrestrial environment. 4‐Coumarate: coenzyme A ligase (4CL) is a key enzyme of general phenylpropanoid metabolism which provides the precursors for both lignin and flavonoids biosynthesis. However, very little is known about how such essential enzymatic functions evolve and diversify. Here, we analyze 4CL sequence variation patterns in a phylogenetic framework to further identify the evolutionary forces that lead to functional divergence. The results reveal that lignin‐biosynthetic 4CLs are under positive selection. The majority of the positively selected sites are located in the substrate‐binding pocket and the catalytic center, indicating that nonsynonymous substitutions might contribute to the functional evolution of 4CLs for lignin biosynthesis. The evolution of 4CLs involved in flavonoid biosynthesis is constrained by purifying selection and maintains the ancestral role of the protein in response to biotic and abiotic factors. Overall, our results demonstrate that protein sequence evolution via positive selection is an important evolutionary force driving adaptive diversification in 4CL proteins in angiosperms. This diversification is associated with adaption to a terrestrial environment.     

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0072-6) contains supplementary material, which is available to authorized users.     

How phenotypic convergence arises in experimental evolution

Evolutionary convergence is a core issue in the study of adaptive evolution, as well as a highly debated topic at present. Few studies have analyzed this issue using a “real‐time” or evolutionary trajectory approach. Do populations that are initially differentiated converge to a similar adaptive state when experiencing a common novel environment? Drosophila subobscura populations founded from different locations and years showed initial differences and variation in evolutionary rates in several traits during short‐term (～20 generations) laboratory adaptation. Here, we extend that analysis to 40 more generations to analyze (1) how differences in evolutionary dynamics among populations change between shorter and longer time spans, and (2) whether evolutionary convergence occurs after 60 generations of evolution in a common environment. We found substantial variation in longer term evolutionary trajectories and differences between short‐ and longer term evolutionary dynamics. Although we observed pervasive patterns of convergence toward the character values of long‐established populations, populations still remain differentiated for several traits at the final generations analyzed. This pattern might involve transient divergence, as we report in some cases, indicating that more generations should lead to final convergence. These findings highlight the importance of longer term studies for understanding convergent evolution.     

Homology modeling of human Toll‐like receptors TLR7, 8, and 9 ligand‐binding domains

Toll‐like receptors (TLRs) play a key role in the innate immune system. The TLR7, 8, and 9 compose a family of intracellularly localized TLRs that signal in response to pathogen‐derived nucleic acids. So far, there are no crystallographic structures for TLR7, 8, and 9. For this reason, their ligand‐binding mechanisms are poorly understood. To enable first predictions of the receptor–ligand interaction sites, we developed three‐dimensional structures for the leucine‐rich repeat ectodomains of human TLR7, 8, and 9 based on homology modeling. To achieve a high sequence similarity between targets and templates, structural segments from all known TLR ectodomain structures (human TLR1/2/3/4 and mouse TLR3/4) were used as candidate templates for the modeling. The resulting models support previously reported essential ligand‐binding residues. They also provide a basis to identify three potential receptor dimerization mechanisms. Additionally, potential ligand‐binding residues are identified using combined procedures. We suggest further investigations of these residues through mutation experiments. Our modeling approach can be extended to other members of the TLR family or other repetitive proteins.     

Tracing the evolution of novel features of human Toll‐like receptor 4

Toll‐like receptor 4 (TLR4) is a critical innate immune protein that activates inflammation in response to extracellular cues. Much of the work to understand how the protein works in humans has been done using mouse models. Although human and mouse TLR4 have many shared features, they have also diverged significantly since their last common ancestor, acquiring 277 sequence differences. Functional differences include the extent of ligand‐independent activation, whether lipid IVa acts as an antagonist or agonist, and the relative species cross‐compatibility of their MD‐2 cofactor. We set out to understand the evolutionary origins for these functional differences between human and mouse TLR4. Using a combination of phylogenetics, ancestral sequence reconstruction, and functional characterization, we found that evolutionary changes to the human TLR4, rather than changes to the mouse TLR4, were largely responsible for these functional changes. Human TLR4 repressed ancestral ligand‐independent activity and gained antagonism to lipid IVa. Additionally, mutations to the human TLR4 cofactor MD‐2 led to lineage‐specific incompatibility between human and opossum TLR4 complex members. These results were surprising, as mouse TLR4 has acquired many more mutations than human TLR4 since their last common ancestor. Our work has polarized this set of transitions and sets up work to study the mechanistic underpinnings for the evolution of new functions in TLR4.     

Convergent selection pressures drive the evolution of rhodopsin kinetics at high altitudes via nonparallel mechanisms

Convergent evolution in response to similar selective pressures is a well‐known phenomenon in evolutionary biology. Less well understood is how selection drives convergence in protein function, and the underlying mechanisms by which this can be achieved. Here, we investigate functional convergence in the visual system of two distantly related lineages of high‐altitude adapted Andean and Himalayan catfishes. Statistical analyses revealed in the two high‐altitude lineages, a parallel acceleration of evolutionary rates in rhodopsin, the dim‐light visual pigment. However, the elevated rates were found to be accompanied by substitutions at different sites in the protein. Experiments substituting Andean‐ or Himalayan‐specific residues significantly accelerated the kinetic rates of rhodopsin, destabilizing the ligand‐bound forms. As found in cold‐adapted enzymes, this phenotype likely compensates for a cold‐induced decrease in kinetic rates, properties of rhodopsin mediating rod sensitivity and visual performance. Our study suggests that molecular convergence in protein function can be driven by parallel shifts in evolutionary rates but via nonparallel molecular mechanisms. Signatures of natural selection may therefore be a powerful guide for identifying complex instances of functional convergence across a wider range of protein systems.     

Molecular evolution of the avian growth hormone gene and comparison with its mammalian counterpart

The molecular evolution of all available avian growth hormone (GH) gene sequences was investigated using both maximum-likelihood and parsimony methods, and the patterns compared to those found in mammals. In contrast to the rapid bursts of evolution observed for mammalian GH, the evolutionary rate of the avian GH mature peptide appears to have been more constant. However several positively selected sites were identified at functionally important positions in the avian signal peptide by the site-specific likelihood method. This implies that sequence variation in the avian GH signal peptide may be adaptive, although more conservative parsimony methods failed to confirm this. Nevertheless, the differing patterns of avian and mammalian GH signal peptide molecular evolution are consistent with the apparently differing roles of GH in controlling growth in these taxonomic groups and support the hypothesis that signal peptide sequence variation may in fact be the basis for increased functional complexity.     

Up‐regulated TLR4 in cardiomyocytes exacerbates heart failure after long‐term myocardial infarction

It remains unclear whether and how cardiomyocytes contribute to the inflammation in chronic heart failure (CHF). We recently reviewed the capacity of cardiomyocytes to initiate inflammation, by means of expressing certain immune receptors such as toll‐like receptors (TLRs) that respond to pathogen‐ and damage‐associated molecular patterns (PAMP and DAMP). Previous studies observed TLR4‐mediated inflammation within days of myocardial infarction (MI). This study examined TLR4 expression and function in cardiomyocytes of failing hearts after 4 weeks of MI in rats. The increases of TLR4 mRNA and proteins, as well as inflammatory cytokine production, were observed in both the infarct and remote myocardium. Enhanced immunostaining for TLR4 was observed in cardiomyocytes but not infiltrating leucocytes. The injection of lentivirus shRNA against TLR4 into the infarcted heart decreased inflammatory cytokine production and improved heart function in vivo. Accordingly, in cardiomyocytes isolated from CHF hearts, increases of TLR4 mRNA and proteins were detected. More robust binding of TLR4 with lipopolysaccharide (LPS), a PAMP ligand for TLR4, and heat shock protein 60 (HSP60), a DAMP ligand for TLR4, was observed in CHF cardiomyocytes under a confocal microscope. The maximum binding capacity (B_max) of TLR4 was increased for LPS and HSP60, whereas the binding affinity (Kd) was not significantly changed. Furthermore, both LPS and HSP60 induced more robust production of inflammatory cytokines in CHF cardiomyocytes, which was reduced by TLR4‐blocking antibodies. We conclude that the expression, ligand‐binding capacity and pro‐inflammatory function of cardiomyocyte TLR4 are up‐regulated after long‐term MI, which promote inflammation and exacerbate heart failure.     

Impact of temperature shifts on the joint evolution of seed dormancy and size

Seed dormancy and size are two important life‐history traits that interplay as adaptation to varying environmental settings. As evolution of both traits involves correlated selective pressures, it is of interest to comparatively investigate the evolution of the two traits jointly as well as independently. We explore evolutionary trajectories of seed dormancy and size using adaptive dynamics in scenarios of deterministic or stochastic temperature variations. Ecological dynamics usually result in unbalanced population structures, and temperature shifts or fluctuations of high magnitude give rise to more balanced ecological structures. When only seed dormancy evolves, it is counter‐selected and temperature shifts hasten this evolution. Evolution of seed size results in the fixation of a given strategy and evolved seed size decreases when seed dormancy is lowered. When coevolution is allowed, evolutionary variations are reduced while the speed of evolution becomes faster given temperature shifts. Such coevolution scenarios systematically result in reduced seed dormancy and size and similar unbalanced population structures. We discuss how this may be linked to the system stability. Dormancy is counter‐selected because population dynamics lead to stable equilibrium, while small seeds are selected as the outcome of size‐number trade‐offs. Our results suggest that unlike random temperature variation between generations, temperature shifts with high magnitude can considerably alter population structures and accelerate life‐history evolution. This study increases our understanding of plant evolution and persistence in the context of climate changes.     

Pathogen richness and abundance predict patterns of adaptive major histocompatibility complex variation in insular amphibians

The identification of the factors responsible for genetic variation and differentiation at adaptive loci can provide important insights into the evolutionary process and is crucial for the effective management of threatened species. We studied the impact of environmental viral richness and abundance on functional diversity and differentiation of the MHC class Ia locus in populations of the black‐spotted pond frog (Pelophylax nigromaculatus), an IUCN‐listed species, on 24 land‐bridge islands of the Zhoushan Archipelago and three nearby mainland sites. We found a high proportion of private MHC alleles in mainland and insular populations, corresponding to 32 distinct functional supertypes, and strong positive selection on MHC antigen‐binding sites in all populations. Viral pathogen diversity and abundance were reduced at island sites relative to the mainland, and islands housed distinctive viral communities. Standardized MHC diversity at island sites exceeded that found at neutral microsatellites, and the representation of key functional supertypes was positively correlated with the abundance of specific viruses in the environment (Frog virus 3 and Ambystoma tigrinum virus). These results indicate that pathogen‐driven diversifying selection can play an important role in maintaining functionally important MHC variation following island isolation, highlighting the importance of considering functionally important genetic variation and host–pathogen associations in conservation planning and management.     

Positive selection underlies the species‐specific binding of Plasmodium falciparum RH5 to human basigin

Plasmodium falciparum, the causative agent of the deadliest form of malaria, is a member of the Laverania subgenus, which includes ape‐infecting parasites. P. falciparum is thought to have originated in gorillas, although infection is now restricted to humans. Laverania parasites display remarkable host‐specificity, which is partially mediated by the interaction between parasite ligands and host receptors. We analyse the evolution of BSG (basigin) and GYPA (glycophorin A) in primates/hominins, as well as of their Plasmodium‐encoded ligands, PfRH5 and PfEBA175. We show that, in primates, positive selection targeted two sites in BSG (F27 and H102), both involved in PfRH5 binding. A population genetics–phylogenetics approach detected the strongest selection for the gorilla lineage: one of the positively selected sites (K191) is a major determinant of PfRH5 binding affinity. Analysis of RH5 genes indicated episodic selection on the P. falciparum branch; the positively selected W447 site is known to stabilize the interaction with human basigin. Conversely, we detect no selection in the receptor‐binding region of EBA175 in the P. falciparum lineage. Its host receptor, GYPA, shows evidence of positive selection in all hominid lineages; selected codons include glycosylation sites that modulate PfEBA175 binding affinity. Data herein provide an evolutionary explanation for species‐specific binding of the PfRH5‐BSG ligand–receptor pair and support the hypothesis that positive selection at these genes drove the host shift leading to the emergence of P. falciparum as a human pathogen.     

Malagasy cichlids differentially limit impacts of body shape evolution on oral jaw functional morphology

Patterns of trait covariation, such as integration and modularity, are vital factors that influence the evolution of vertebrate body plans. In functional systems, decoupling of morphological modules buffers functional change in one trait by reducing correlated variation with another. However, for complex morphologies with many‐to‐one mapping of form to function (MTOM), resistance to functional change may also be achieved by constraining morphological variation within a functionally stable region of morphospace. For this research, we used geometric morphometrics to evaluate the evolution of body shape and its relationship with jaw functional morphology in two independent radiations of endemic Malagasy cichlid (Teleostei: Cichlidae). Our results suggested that the two subfamilies used different strategies to mitigate impacts of body shape variation on a metric of jaw function, maxillary kinematic transmission (MKT): (1) modularity between cranial and postcranial morphologies, and (2) integration of body and jaw evolution, with jaw morphologies varying in a manner that limits change in MKT. This research shows that, unlike modularity, MTOM allows traits to retain strong evolutionary covariation while still reducing impacts on functionality. These results suggest that MTOM, and its influence on the evolution of correlated traits, is likely much more widespread than is currently understood.     

Imprint of evolutionary conservation and protein structure variation on the binding function of protein tyrosine kinases

MOTIVATION: According to the models of divergent molecular evolution, the evolvability of new protein function may depend on the induction of new phenotypic traits by a small number of mutations of the binding site residues. Evolutionary relationships between protein kinases are often employed to infer inhibitor binding profiles from sequence analysis. However, protein kinases binding profiles may display inhibitor selectivity within a given kinase subfamily, while exhibiting cross-activity between kinases that are phylogenetically remote from the prime target. The emerging insights into kinase function and evolution combined with a rapidly growing number of publically available crystal structures of protein kinases complexes have motivated structural bioinformatics analysis of sequence-structure relationships in determining the binding function of protein tyrosine kinases. RESULTS: In silico profiling of Imatinib mesylate and PD-173955 kinase inhibitors with protein tyrosine kinases is conducted on kinome scale by using evolutionary analysis and fingerprinting inhibitor-protein interactions with the panel of all publically available protein tyrosine kinases crystal structures. We have found that sequence plasticity of the binding site residues alone may not be sufficient to enable protein tyrosine kinases to readily evolve novel binding activities with inhibitors. While evolutionary signal derived solely from the tyrosine kinase sequence conservation can not be readily translated into the ligand binding phenotype, the proposed structural bioinformatics analysis can discriminate a functionally relevant kinase binding signal from a simple phylogenetic relationship. The results of this work reveal that protein conformational diversity is intimately linked with sequence plasticity of the binding site residues in achieving functional adaptability of protein kinases towards specific drug binding. This study offers a plausible molecular rationale to the experimental binding profiles of the studied kinase inhibitors and provides a theoretical basis for constructing functionally relevant kinase binding trees.     

Small‐scale spatial variation in infection risk shapes the evolution of a Borrelia resistance gene in wild rodents

Spatial variation in pathogen‐mediated selection is predicted to influence the evolutionary trajectory of host populations and lead to spatial variation in their immunogenetic composition. However, to date few studies have been able to directly link small‐scale spatial variation in infection risk to host immune gene evolution in natural, nonhuman populations. Here, we use a natural rodent–Borrelia system to test for associations between landscape‐level spatial variation in Borrelia infection risk along replicated elevational gradients in the Swiss Alps and Toll‐like receptor 2 (TLR2) evolution, a candidate gene for Borrelia resistance, across bank vole (Myodes glareolus) populations. We found that Borrelia infection risk (i.e., the product of Borrelia prevalence in questing ticks and the average tick load of voles at a sampling site) was spatially variable and significantly negatively associated with elevation. Across sampling sites, Borrelia prevalence in bank voles was significantly positively associated with Borrelia infection risk along the elevational clines. We observed a significant association between naturally occurring TLR2 polymorphisms in hosts and their Borrelia infection status. The TLR2 variant associated with a reduced likelihood of Borrelia infection was most common in rodent populations at lower elevations that face a high Borrelia infection risk, and its frequency changed in accordance with the change in Borrelia infection risk along the elevational clines. These results suggest that small‐scale spatial variation in parasite‐mediated selection affects the immunogenetic composition of natural host populations, providing a striking example that the microbial environment shapes the evolution of the host's immune system in the wild.     

Phenotypic integration: studying the ecology and evolution of complex phenotypes

Phenotypic integration refers to the study of complex patterns of covariation among functionally related traits in a given organism. It has been investigated throughout the 20th century, but has only recently risen to the forefront of evolutionary ecological research. In this essay, I identify the reasons for this late flourishing of studies on integration, and discuss some of the major areas of current endeavour: the interplay of adaptation and constraints, the genetic and molecular bases of integration, the role of phenotypic plasticity, macroevolutionary studies of integration, and statistical and conceptual issues in the study of the evolution of complex phenotypes. I then conclude with a brief discussion of what I see as the major future directions of research on phenotypic integration and how they relate to our more general quest for the understanding of phenotypic evolution within the neo‐Darwinian framework. I suggest that studying integration provides a particularly stimulating and truly interdisciplinary convergence of researchers from fields as disparate as molecular genetics, developmental biology, evolutionary ecology, palaeontology and even philosophy of science.     

Toll-like receptor 5 forms asymmetric dimers in the absence of flagellin

The structure of full-length human TLR5 determined by electron microscopy single-particle image reconstruction at 26 Å resolution shows that TLR5 forms an asymmetric homodimer via ectodomain interactions. The structure shows that like TLR9, TLR5 dimerizes in the absence of ligand. The asymmetry of the dimer suggests that TLR5 may recognize two flagellin molecules cooperatively to establish an optimal flagellin response threshold. A TLR5 homology model was generated and fitted into the electron microscopy structure. All seven predicted N-linked glycosylation sites are exposed on the molecular surface, away from the dimer interface. Glycosylation at the first five sites was confirmed by tandem mass spectrometry. Two aspartate residues proposed to interact with flagellin (Asp294 and Asp366) are sterically occluded by a glycan at position 342. In contrast, the central region of the ectodomains near the dimer interface is unobstructed by glycans. Ligand binding in this region would be consistent with the ligand binding sites of other TLRs.