Evidence of positive selection at codon sites localized in extracellular domains of mammalian CC motif chemokine receptor proteins

Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events

Chemokine receptors are an important determinant for the infectiousness of different pathogens, which are able to target the host cells by binding to the extracellular domains of these proteins. This is the mechanism of infection of HIV-1, among other concerning human diseases. Over the past years, it has been shown that two chemokine receptors, CCR2 and CCR5, have been shaped by events of gene conversion in different mammalian lineages, which has been linked to a possible selective advantage against pathogens. Here, by taking advantage of available bat genomes, we present the first insight of CCR2 and CCR5 evolution within the Chiroptera order. In total, four independent events of recombination between CCR2 and CCR5 were detected: two in a single species, Miniopterus natalensis; one in two species from the Rhinolophoidea superfamily; and one in four species from the Pteropodidae family. The regions affected by the gene conversions were generally extensive and always encompassed extracellular domains. Overall, we demonstrate that CCR2 and CCR5 have been subject to extensive gene conversion in multiple species of bats. Considering that bats are known to be large reservoirs of virus in nature, these results might indicate that chimeric CCR2-CCR5 genes might grant some bat species a selective advantage against viruses that rely in the extracellular portions of either CCR2 or CCR5 as gateways into the cell.

Genome-wide identification and analysis of chemokine receptor superfamily in miiuy croaker, Miichthys miiuy

The chemokine receptor (ChemR) superfamily, which is divided into 4 subfamilies (CXCR, CCR, XCR, and CX3CR), is the main receptors of chemokines in innate immune responses. In the current study, we have identified 27 ChemRs in miiuy croaker: 13 CCR genes, 11 CXCR genes, and 3 XCR genes. Multiple characteristics of these genes, including phylogeny, gene structures, conserved motifs, chromosome locations, evolutionary mechanism, and expression levels upon the bacterial challenge were analyzed. Gene structure and location analysis showed that all ChemR genes contain fewer introns (≤4) and they are unevenly distributed on the 12 chromosomes. And the XCR subfamily of miiuy croaker don't have the DRY motif of ChemR. Phylogenetic and synteny analysis showed that these genes experienced tandem and segmental duplication event in several species, and tandem duplication might be the main expansion way in miiuy croaker. The major ChemRs of each orthologous group in vertebrates were selected for molecular evolution analysis, the results of which indicated that compared with vertebrates, ChemRs of teleost fishes may have a relatively high evolutionary dynamic. In addition, a total of 21 positively selected codons were detected in vertebrate ChemRs under Model 8. RNA-Seq analysis and qRT-PCR verification demonstrated that CXCR3.2, CXCR5, and XCR1 genes were up-regulated significantly upon the Vibrio harveyi infection. These results provide valuable information for investigating the evolutionary relationships of chemokine receptor superfamily in miiuy croaker and laid the basis for further functional analysis.

Unraveling the evolutionary origin of ELR motif using fish CXC chemokine CXCL8

Chemokines are chemotactic proteins involved in host defense through the migration of immune-regulatory cells to the site of infection. Interleukin-8 (CXCL8/IL8) is the most studied "ELR-CXC chemokine/neutrophil activating chemokine (NAC) that regulate neutrophil trafficking during infections and inflammation by binding to its cognate G-protein coupled receptors CXCR1/CXCR2. The "ELR" motif of NAC chemokines is essential for the CXCR1/CXCR2 receptor activation. In order to understand the evolutionary origin of "ELR" motif in the CXC chemokines, a thorough evolutionary study of CXCL8 gene from various fishes and primates was performed. Phylogenetic analysis revealed that the CXCL8 gene can be classified into four distinct lineages (CXCL8-L1a, CXCL8-L1b, CXCL8-L2, and CXCL8-L3), where CXCL8-L1a is the fastest evolving lineage and CXCL8-L3 is the slowest. Selection analysis suggested that The "ELR/DLR" motif containing branches (gadoid and coelacanth) are positively selected. The probable evolutionary trend of "ELR" motif suggested that this motif in ancestor CXCL8 is evolved from the GGR of Lamprey (Agnatha), followed by duplication giving rise to two main motifs in CXCL8 "NXH" in L3 lineage and "ELR/DLR" in L1a/L1b lineages. Although, structural analysis suggested that the overall topology of the CXCL8 proteins is similar, differences do exist at the individual structural elements among the members of different lineages. Functional distance analysis suggested that the CXCL8-L3 lineage is more distant compared to the CXCL8-L1a and L1b lineages from the inferred ancestor. Functional divergence analysis between different lineages suggested that most of the selected residues are important for receptor or glycosaminoglycan binding. Such a functional diversification can be attributed to the novel set of functions adopted by CXCL8 in various species.

Large Scale Analyses and Visualization of Adaptive Amino Acid Changes Projects

When changes at few amino acid sites are the target of selection, adaptive amino acid changes in protein sequences can be identified using maximum-likelihood methods based on models of codon substitution (such as codeml). Although such methods have been employed numerous times using a variety of different organisms, the time needed to collect the data and prepare the input files means that tens or hundreds of coding regions are usually analyzed. Nevertheless, the recent availability of flexible and easy to use computer applications that collect relevant data (such as BDBM) and infer positively selected amino acid sites (such as ADOPS), means that the entire process is easier and quicker than before. However, the lack of a batch option in ADOPS, here reported, still precludes the analysis of hundreds or thousands of sequence files. Given the interest and possibility of running such large-scale projects, we have also developed a database where ADOPS projects can be stored. Therefore, this study also presents the B+ database, which is both a data repository and a convenient interface that looks at the information contained in ADOPS projects without the need to download and unzip the corresponding ADOPS project file. The ADOPS projects available at B+ can also be downloaded, unzipped, and opened using the ADOPS graphical interface. The availability of such a database ensures results repeatability, promotes data reuse with significant savings on the time needed for preparing datasets, and effortlessly allows further exploration of the data contained in ADOPS projects.

Evolution of toll-like receptors in the context of terrestrial ungulates and cetaceans diversification

Background

Toll-like receptors (TLRs) are the frontline actors in the innate immune response to various pathogens and are expected to be targets of natural selection in species adapted to habitats with contrasting pathogen burdens. The recent publication of genome sequences of giraffe and okapi together afforded the opportunity to examine the evolution of selected TLRs in broad range of terrestrial ungulates and cetaceans during their complex habitat diversification. Through direct sequence comparisons and standard evolutionary approaches, the extent of nucleotide and protein sequence diversity in seven Toll-like receptors (TLR2, TLR3, TLR4, TLR5, TLR7, TLR9 and TLR10) between giraffe and closely related species was determined. In addition, comparison of the patterning of key TLR motifs and domains between giraffe and related species was performed. The quantification of selection pressure and divergence on TLRs among terrestrial ungulates and cetaceans was also performed.

Results

Sequence analysis shows that giraffe has 94–99% nucleotide identity with okapi and cattle for all TLRs analyzed. Variations in the number of Leucine-rich repeats were observed in some of TLRs between giraffe, okapi and cattle. Patterning of key TLR domains did not reveal any significant differences in the domain architecture among giraffe, okapi and cattle. Molecular evolutionary analysis for selection pressure identifies positive selection on key sites for all TLRs examined suggesting that pervasive evolutionary pressure has taken place during the evolution of terrestrial ungulates and cetaceans. Analysis of positively selected sites showed some site to be part of Leucine-rich motifs suggesting functional relevance in species-specific recognition of pathogen associated molecular patterns. Notably, clade analysis reveals significant selection divergence between terrestrial ungulates and cetaceans in viral sensing TLR3. Mapping of giraffe TLR3 key substitutions to the structure of the receptor indicates that at least one of giraffe altered sites coincides with TLR3 residue known to play a critical role in receptor signaling activity.

Conclusion

There is overall structural conservation in TLRs among giraffe, okapi and cattle indicating that the mechanism for innate immune response utilizing TLR pathways may not have changed very much during the evolution of these species. However, a broader phylogenetic analysis revealed signatures of adaptive evolution among terrestrial ungulates and cetaceans, including the observed selection divergence in TLR3. This suggests that long term ecological dynamics has led to species-specific innovation and functional variation in the mechanisms mediating innate immunity in terrestrial ungulates and cetaceans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0901-7) contains supplementary material, which is available to authorized users.

Selective Landscapes in newt Immune Genes Inferred from Patterns of Nucleotide Variation

Host–pathogen interactions may result in either directional selection or in pressure for the maintenance of polymorphism at the molecular level. Hence signatures of both positive and balancing selection are expected in immune genes. Because both overall selective pressure and specific targets may differ between species, large-scale population genomic studies are useful in detecting functionally important immune genes and comparing selective landscapes between taxa. Such studies are of particular interest in amphibians, a group threatened worldwide by emerging infectious diseases. Here, we present an analysis of polymorphism and divergence of 634 immune genes in two lineages of Lissotriton newts: L. montandoni and L. vulgaris graecus. Variation in newt immune genes has been shaped predominantly by widespread purifying selection and strong evolutionary constraint, implying long-term importance of these genes for functioning of the immune system. The two evolutionary lineages differ in the overall strength of purifying selection which can partially be explained by demographic history but may also signal differences in long-term pathogen pressure. The prevalent constraint notwithstanding, 23 putative targets of positive selection and 11 putative targets of balancing selection were identified. The latter were detected by composite tests involving the demographic model and further validated in independent population samples. Putative targets of balancing selection encode proteins which may interact closely with pathogens but include also regulators of immune response. The identified candidates will be useful for testing whether genes affected by balancing selection are more prone to interspecific introgression than other genes in the genome.

Evidence of positive selection associated with placental loss in tiger sharks

BACKGROUND

All vertebrates initially feed their offspring using yolk reserves. In some live-bearing species these yolk reserves may be supplemented with extra nutrition via a placenta. Sharks belonging to the Carcharhinidae family are all live-bearing, and with the exception of the tiger shark (Galeocerdo cuvier), develop placental connections after exhausting yolk reserves. Phylogenetic relationships suggest the lack of placenta in tiger sharks is due to secondary loss. This represents a dramatic shift in reproductive strategy, and is likely to have left a molecular footprint of positive selection within the genome.

RESULTS

We sequenced the transcriptome of the tiger shark and eight other live-bearing shark species. From this data we constructed a time-calibrated phylogenetic tree estimating the tiger shark lineage diverged from the placental carcharhinids approximately 94 million years ago. Along the tiger shark lineage, we identified five genes exhibiting a signature of positive selection. Four of these genes have functions likely associated with brain development (YWHAE and ARL6IP5) and sexual reproduction (VAMP4 and TCTEX1D2).

CONCLUSIONS

Our results indicate the loss of placenta in tiger sharks may be associated with subsequent adaptive changes in brain development and sperm production.

Evolution of CCL11: genetic characterization in lagomorphs and evidence of positive and purifying selection in mammals

The interactions between chemokines and their receptors are crucial for differentiation and activation of inflammatory cells. CC chemokine ligand 11 (CCL11) binds to CCR3 and to CCR5 that in leporids underwent gene conversion with CCR2. Here, we genetically characterized CCL11 in lagomorphs (leporids and pikas). All lagomorphs have a potentially functional CCL11, and the Pygmy rabbit has a mutation in the stop codon that leads to a longer protein. Other mammals also have mutations at the stop codon that result in proteins with different lengths. By employing maximum likelihood methods, we observed that, in mammals, CCL11 exhibits both signatures of purifying and positive selection. Signatures of purifying selection were detected in sites important for receptor binding and activation. Of the three sites detected as under positive selection, two were located close to the stop codon. Our results suggest that CCL11 is functional in all lagomorphs, and that the signatures of purifying and positive selection in mammalian CCL11 probably reflect the protein's biological roles.

Signatures of natural selection at the FTO (fat mass and obesity associated) locus in human populations

Background and aims

Polymorphisms in the first intron of FTO have been robustly replicated for associations with obesity. In the Sorbs, a Slavic population resident in Germany, the strongest effect on body mass index (BMI) was found for a variant in the third intron of FTO (rs17818902). Since this may indicate population specific effects of FTO variants, we initiated studies testing FTO for signatures of selection in vertebrate species and human populations.

Methods

First, we analyzed the coding region of 35 vertebrate FTO orthologs with Phylogenetic Analysis by Maximum Likelihood (PAML, ω = dN/dS) to screen for signatures of selection among species. Second, we investigated human population (Europeans/CEU, Yoruba/YRI, Chinese/CHB, Japanese/JPT, Sorbs) SNP data for footprints of selection using DnaSP version 4.5 and the Haplotter/PhaseII. Finally, using ConSite we compared transcription factor (TF) binding sites at sequences harbouring FTO SNPs in intron three.

Results

PAML analyses revealed strong conservation in coding region of FTO (ω<1). Sliding-window results from population genetic analyses provided highly significant (p<0.001) signatures for balancing selection specifically in the third intron (e.g. Tajima’s D in Sorbs = 2.77). We observed several alterations in TF binding sites, e.g. TCF3 binding site introduced by the rs17818902 minor allele.

Conclusion

Population genetic analysis revealed signatures of balancing selection at the FTO locus with a prominent signal in intron three, a genomic region with strong association with BMI in the Sorbs. Our data support the hypothesis that genes associated with obesity may have been under evolutionary selective pressure.

The contribution of natural selection to present-day susceptibility to chronic inflammatory and autoimmune disease

Chronic inflammatory and autoimmune diseases have been the focus of many genome-wide association studies (GWAS) because they represent a significant cause of illness and morbidity, and many are heritable. Almost a decade of GWAS studies suggests that the pathological inflammation associated with these diseases is controlled by a limited number of networked immune system genes. Chronic inflammatory and autoimmune diseases are enigmatic from an evolutionary perspective because they exert a negative affect on reproductive fitness. The persistence of these conditions may be partially explained by the important roles the implicated immune genes play in pathogen defense and other functions thought to be under strong natural selection in humans. The evolutionary reasons for chronic inflammatory and autoimmune disease persistence and uneven distribution across populations are the focus of this review.

Tracing evolutionary relicts of positive selection on eight malaria-related immune genes in mammals

Plasmodium-induced malaria widely infects primates and other mammals. Multiple past studies have revealed that positive selection could be the main evolutionary force triggering the genetic diversity of anti-malaria resistance-associated genes in human or primates. However, researchers focused most of their attention on the infra-generic and intra-specific genome evolution rather than analyzing the complete evolutionary history of mammals. Here we extend previous research by testing the evolutionary link of natural selection on eight candidate genes associated with malaria resistance in mammals. Three of the eight genes were detected to be affected by recombination, including TNF-α, iNOS and DARC. Positive selection was detected in the rest five immunogenes multiple times in different ancestral lineages of extant species throughout the mammalian evolution. Signals of positive selection were exposed in four malaria-related immunogenes in primates: CCL2, IL-10, HO1 and CD36. However, selection signals of G6PD have only been detected in non-primate eutherians. Significantly higher evolutionary rates and more radical amino acid replacement were also detected in primate CD36, suggesting its functional divergence from other eutherians. Prevalent positive selection throughout the evolutionary trajectory of mammalian malaria-related genes supports the arms race evolutionary hypothesis of host genetic response of mammalian immunogenes to infectious pathogens.

Simulation of genome-wide evolution under heterogeneous substitution models and complex multispecies coalescent histories

Genomic evolution can be highly heterogeneous. Here, we introduce a new framework to simulate genome-wide sequence evolution under a variety of substitution models that may change along the genome and the phylogeny, following complex multispecies coalescent histories that can include recombination, demographics, longitudinal sampling, population subdivision/species history, and migration. A key aspect of our simulation strategy is that the heterogeneity of the whole evolutionary process can be parameterized according to statistical prior distributions specified by the user. We used this framework to carry out a study of the impact of variable codon frequencies across genomic regions on the estimation of the genome-wide nonsynonymous/synonymous ratio. We found that both variable codon frequencies across genes and rate variation among sites and regions can lead to severe underestimation of the global dN/dS values. The program SGWE—Simulation of Genome-Wide Evolution—is freely available from http://code.google.com/p/sgwe-project/, including extensive documentation and detailed examples.

Characterization of the CCR3 and CCR9 genes in miiuy croaker and different selection pressures imposed on different domains between mammals and teleosts

The innate immune system can recognize non-self through pattern recognition receptors and provides a first line of antimicrobial host defense. Thus innate immunity plays a very important role in resistance against major bacterial disease in vertebrates. In the innate immune responses, the chemokine receptors act as the main receptors of chemokines which are released at the sites of infection, inflammation and injury. In this study, the Miichthys miiuy CCR3 and CCR9 genes were cloned and characterized, showing that MIMI-CCR3 possesses a highly conserved DRYLA motif similar to that of other fishes. MIMI-CCR3 and CCR9 were ubiquitously expressed in all tested tissues and the expressions were significantly up-regulated after infection with Vibrio anguillarum except that of CCR9 in spleen. Evolutionary analysis showed that all the ancestral lineages of CCR3 and CCR9 in mammals and teleosts underwent positive selection, indicating that the ancestor of terrestrial animals further evolved to adapt to terrestrial environments and the continuous intrusion of microbes stimulated the evolution of CCR genes in the ancestor of teleost. Multiple ML methods were used to detect the robust candidates for sites under positive selection. In total, 11 and 8 positively selected sites were found in the subsets of current mammal and teleost CCR3 genes, and 3 and 2 sites were detected in the subsets of current mammals and teleosts in CCR9. Interestingly, for mammal CCR3 and CCR9 genes, the robust candidates of positively selected sites were mainly located in the extracellular domains which were the ligand binding and pathogen interaction regions. For teleost CCR3 and CCR9 genes, the positively selected sites were not only located in the extracellular domains but also in the C-terminal and intracellular domains, indicating mammals and teleosts experienced different selection pressures upon their N-terminus, C-terminus and intracellular loops of CCRs.

Directional Darwinian Selection in proteins

Background

Molecular evolution is a very active field of research, with several complementary approaches, including dN/dS, HON90, MM01, and others. Each has documented strengths and weaknesses, and no one approach provides a clear picture of how natural selection works at the molecular level. The purpose of this work is to present a simple new method that uses quantitative amino acid properties to identify and characterize directional selection in proteins.

Methods

Inferred amino acid replacements are viewed through the prism of a single physicochemical property to determine the amount and direction of change caused by each replacement. This allows the calculation of the probability that the mean change in the single property associated with the amino acid replacements is equal to zero (H₀: μ = 0; i.e., no net change) using a simple two-tailed t-test.

Results

Example data from calanoid and cyclopoid copepod cytochrome oxidase subunit I sequence pairs are presented to demonstrate how directional selection may be linked to major shifts in adaptive zones, and that convergent evolution at the whole organism level may be the result of convergent protein adaptations.

Conclusions

Rather than replace previous methods, this new method further complements existing methods to provide a holistic glimpse of how natural selection shapes protein structure and function over evolutionary time.

Evolution of vertebrate immunity

All multicellular organisms protect themselves against pathogens using sophisticated immune defenses. Functionally interconnected humoral and cellular facilities maintain immune homeostasis in the absence of overt infection and regulate the initiation and termination of immune responses directed against pathogens. Immune responses of invertebrates, such as flies, are innate and usually stereotyped; those of vertebrates, encompassing species as diverse as jawless fish and humans, are additionally adaptive, enabling more rapid and efficient immune reactivity upon repeated encounters with a pathogen. Many of the attributes historically defining innate and adaptive immunity are in fact common to both, blurring their functional distinction and emphasizing shared ancestry and co-evolution. These findings provide indications of the evolutionary forces underlying the origin of somatic diversification of antigen receptors and contribute to our understanding of the complex phenotypes of human immune disorders. Moreover, informed by phylogenetic considerations and inspired by improved knowledge of functional networks, new avenues emerge for innovative therapeutic strategies.

Species-specific chitin-binding module 18 expansion in the amphibian pathogen Batrachochytrium dendrobatidis

Batrachochytrium dendrobatidis is the causative agent of chytridiomycosis, which is considered one of the driving forces behind the worldwide decline in populations of amphibians. As a member of the phylum Chytridiomycota, B. dendrobatidis has diverged significantly to emerge as the only pathogen of adult vertebrates. Such shifts in lifestyle are generally accompanied by various degrees of genomic modifications, yet neither its mode of pathogenicity nor any factors associated with it have ever been identified. Presented here is the identification and characterization of a unique expansion of the carbohydrate-binding module family 18 (CBM18), specific to B. dendrobatidis. CBM (chitin-binding module) expansions have been likened to the evolution of pathogenicity in a variety of fungus species, making this expanded group a prime candidate for the identification of potential pathogenicity factors. Furthermore, the CBM18 expansions are confined to three categories of genes, each having been previously implicated in host-pathogen interactions. These correlations highlight this specific domain expansion as a potential key player in the mode of pathogenicity in this unique fungus. The expansion of CBM18 in B. dendrobatidis is exceptional in its size and diversity compared to other pathogenic species of fungi, making this genomic feature unique in an evolutionary context as well as in pathogenicity.

Amphibian populations are declining worldwide at an unprecedented rate. Although various factors are thought to contribute to this phenomenon, chytridiomycosis has been identified as one of the leading causes. This deadly fungal disease is cause by Batrachochytrium dendrobatidis, a chytrid fungus species unique in its pathogenicity and, furthermore, its specificity to amphibians. Despite more than two decades of research, the biology of this fungus species and its deadly interaction with amphibians had been notoriously difficult to unravel. Due to the alarming rate of worldwide spread and associated decline in amphibian populations, it is imperative to incorporate novel genomic and genetic techniques into the study of this species. In this study, we present the first reported potential pathogenicity factors in B. dendrobatidis. In silico studies such as this allow us to identify putative targets for more specific molecular analyses, furthering our hope for the control of this pathogen.

Insight into pattern of codon biasness and nucleotide base usage in serotonin receptor gene family from different mammalian species

5-HT (5-Hydroxy-tryptamine) or serotonin receptors are found both in central and peripheral nervous system as well as in non-neuronal tissues. In the animal and human nervous system, serotonin produces various functional effects through a variety of membrane bound receptors. In this study, we focus on 5-HT receptor family from different mammals and examined the factors that account for codon and nucleotide usage variation. A total of 110 homologous coding sequences from 11 different mammalian species were analyzed using relative synonymous codon usage (RSCU), correspondence analysis (COA) and hierarchical cluster analysis together with nucleotide base usage frequency of chemically similar amino acid codons. The mean effective number of codon (ENc) value of 37.06 for 5-HT(6) shows very high codon bias within the family and may be due to high selective translational efficiency. The COA and Spearman's rank correlation reveals that the nucleotide compositional mutation bias as the major factors influencing the codon usage in serotonin receptor genes. The hierarchical cluster analysis suggests that gene function is another dominant factor that affects the codon usage bias, while species is a minor factor. Nucleotide base usage was reported using Goldman, Engelman, Stietz (GES) scale reveals the presence of high uracil (>45%) content at functionally important hydrophobic regions. Our in silico approach will certainly help for further investigations on critical inference on evolution, structure, function and gene expression aspects of 5-HT receptors family which are potential antipsychotic drug targets.

Evidence for positive selection in the extracellular domain of human cytomegalovirus encoded G protein-coupled receptor US28

The human cytomegalovirus (HCMV)-encoded chemokine receptor US28 is also a seven-transmembrane G-protein coupled receptor, whose signaling pathway is known for its involvement in host immune system evasion. HCMV infection can result in serious disease in immunocompromised individuals and is also linked to atherosclerosis and cardiovascular disease. Identifying amino acid residues that play a crucial role in successful viral adaptation in response to the host's immune defense is critical for effective drug design. In this study maximum likelihood-based codon substitution analyses were carried out to determine whether any codon of US28 has evolved adaptively. If the rate of nonsynonymous (dn) to the rate of synonymous (ds) nucleotide substitutions (ω = dn/ds) is greater than one, the codon is said to be under positive selection, indicating adaptive evolution. Although the overall ω for US28 gene was 0.154, indicating that most codon sites were subject to strong purifying selection, five codon sites are under strong positive selection. Three (E18D/L, D19A/E/G, and R267K/Q) of these positively selected sites are located in extracellular domains, the domains that play a crucial role for successful viral adaptation in response to the host's immune defense. The C-terminal (R329Q/W) and the fifth transmembrane domain (V190I), each have one positively selected site. These results suggest that relative to the extracellular domains, amino acid residues present in intracellular domains are more selectively constrained. A few amino acid residues in extracellular domains of US28 evolved more rapidly, presumably due to positive selection pressure resulting from ligand-binding and pathogen interactions of extracellular domains.