# Maximum-likelihood approaches reveal signatures of positive selection in IL genes in mammals

Fabiana Neves, Joana Abrantes, John W Steinke, Pedro J Esteves

<figure><img src="image_1.png" /></figure>

## To cite this version:

Fabiana Neves, Joana Abrantes, John W Steinke, Pedro J Esteves. Maximum-likelihood approaches reveal signatures of positive selection in IL genes in mammals. Innate Immunity, 2013, 20 (2), pp.184 - 191. 10.1177/1753425913486687. hal-04969724

HAL Id: hal-04969724

https://hal.science/hal-04969724v1

Submitted on 27 Feb 2025

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

SAGE

# Maximum-likelihood approaches reveal signatures of positive selection in IL genes in mammals

Fabiana Neves<fnref n="1" />,<fnref n="2" />, Joana Abrantes<fnref n="1" />,<fnref n="3" />, John W Steinke<fnref n="4" /> and Pedro J Esteves<fnref n="1" />,<fnref n="5" />

## Abstract

ILs are part of the immune system and are involved in multiple biological activities. ILs have been shown to evolve under positive selection; however, little information exists regarding which codons are specifically selected. By using different codon-based maximum-likelihood (ML) approaches, signatures of positive selection in mammalian ILs were searched for. Sequences of 46 ILs were retrieved from publicly available databases of mammalian genomes to detect signatures of positive selection in individual codons. Evolutionary analyses were conducted under two ML frameworks, the HyPhy package implemented in the Data Monkey Web Server and CODEML implemented in PAML. Signatures of positive selection were found in 28 ILs: IL-1A and B; IL-2, IL-4 to IL-10, IL-12A and B; IL-14 to IL-17A and C; IL-18, IL-20 to IL-22, IL-25, IL-26, IL-27B, IL-31, IL-34, IL-36A; and G. Codons under positive selection varied between 1 and 15. No evidence of positive selection was detected in IL-13; IL-17B and F; IL-19, IL-23, IL-24, IL-27A; or IL-29. Most mammalian ILs have sites evolving under positive selection, which may be explained by the multitude of biological processes in which ILs are enrolled. The results obtained raise hypotheses concerning the ILs functions, which should be pursued by using mutagenesis and crystallographic approaches.

## Keywords

ILs, IL receptors, maximum likelihood, positive selection

Date received: 15 October 2012; revised: 8 January 2013; accepted: 27 January 2013

## Introduction

Different cells of the immune system are able to secrete regulatory proteins, namely cytokines (cks) in response to a variety of stimuli. In mammals, cks produced as a part of the innate immune response have the ability to influence the extent and nature of the adaptive immunity response, and are thus crucial for many aspects of the immune response.<fnref n="1" />⁻⁶ When cks are secreted by leukocytes and act on other cells they are called IL.<fnref n="3" />,⁷⁻¹⁰ ILs are polypeptides of low molecular mass involved in several biological activities, including immunity, inflammation, inflammatory diseases, hematopoiesis, oncogenesis, neurogenesis and fertility, among many others. Each IL is normally involved in a multiple biological processes.⁹,¹¹⁻¹³ Currently, 37 ILs have been identified in mammals, some with a variable number of variants. Classification of ILs has been based on sequence homology, receptor chain similarities, activity, structural or functional features, leading

to complex and intricate classifications.⁸,⁹,¹⁴⁻²⁰ The activities of ILs are possible through the binding of these polypeptides to specific cell-surface receptors, components of which can be shared by several ILs, which are able to transmit intracellular signals through different signaling pathways.<fnref n="1" />,¹¹,¹³,²¹ Sometimes, for optimal function, interaction between complementary

## Corresponding author

Pedro J. Esteves, CIBIO-UP, Centro de Investigação em Biodiversidade e Recursos Genéticos, Universidade do Porto, Campus Agrário de Vairão, Rua Padre Armando Quintas, no. 7, 4485-661 Vairão, Portugal.
Email: pjesteves@cibio.up.pt

<footnote n="1">¹CIBIO/UP – Centro de Investigação em Biodiversidade e Recursos Genéticos, Universidade do Porto, Porto, Portugal</footnote>

<footnote n="2">²UMIB/UP – Unidade Multidisciplinar de Investigação Biomédica, Universidade do Porto, Porto, Portugal</footnote>

<footnote n="3">³INSERM, U892, Université de Nantes, Nantes, France</footnote>

<footnote n="4">⁴University of Virginia, Charlottesville, VA, USA</footnote>

<footnote n="5">⁵CITS – Centro de Investigação em Tecnologias de Saúde, CESPU, Gandra, Portugal</footnote>
ILs is necessary.²² These interactions can be synergistic, additive or antagonistic, and exhibit both negative and positive regulatory effects.⁹,¹¹,²²⁻²⁴

Owing to the host-pathogen co-evolution, the immune system and its genes are constantly evolving, being under constant pressure and selection for adaptation, and where advantageous mutations are highly favored and deleterious mutations are quickly eliminated.⁵,²⁵ ILs such as IL-2, IL-3, IL-4, IL-5, IL-13, IL-23A, IL-28A, IL-28B and IL-29¹,²,⁴,²⁶⁻³¹ have been identified as some of the immune system genes under positive selection in different mammals. Despite all of these studies, most do not specify the codons positively selected,²,⁴,²⁷,²⁹⁻³¹ with the exception of IL-4²⁶ in which 15 residues were detected as being under positive selection and located on sites responsible for binding to receptors. With these results in mind, we have extended the search for signatures of positive selection to include mammalian ILs 1-37 by using different maximum-likelihood approaches (ML).

# Materials and methods

## Sequences

The sequences of the mammalian ILs used in the analyses were retrieved from GenBank (http://www.ncbi.nlm.nih.gov/), Ensembl (http://useast.ensembl.org/index.html) and UniProt (http://www.uniprot.org/). For each IL, amino acid residues were numbered from the first human methionine residue, with signal peptides and propeptide amino acids included in the numbering. The number of sequences retrieved for each IL ranged from 14 to 47 (some of the IL genes were not found in all mammalian genomes) and included representatives of most mammalian groups (e.g. Artiodactyla, Carnivores, Lagomorphs, Primates, Rodents, etc). The identification of the species used for each IL and the accession numbers are listed in Appendix 1 (Supplementary material). Each of the 37 ILs were aligned using ClustalW, as implemented in the software program BioEdit version 7.1.3,³² and adjusted manually.

## Codon-based analyses of positive selection

According to Poon and collaborators,⁶ and in order to reliably identify codons under positive selection, only ILs that were represented in at least 10 species were analyzed. Hence, IL-17D, IL-28A and B, IL-32, IL-36B and IL-37 were not considered. IL-3 and IL-33 were also excluded as the alignments produced were not reliable and prone to affect the detection of positive selection by leading to false predictions.³³ Signatures of positive selection are inferred if the ratio of non-synonymous substitutions per non-synonymous substitutions site ($d_N$) over synonymous substitutions per

synonymous sites ($d_S$) is statistically higher than the value observed under neutrality, 1. Here, to detect such signatures in individual codons of mammalian IL sequences, the different $d_N/d_S$ ratios ($\omega$) were compared using two ML frameworks, the HyPhy package implemented in the Data Monkey Web Server (http://www.datamonkey.org/)⁶ and CODEML implemented in PAML version 4,³⁴,³⁵ being considered in the analysis the results where $d_N/d_S$ ratios were significantly higher than 1.

In the Data Monkey Web Server, the best-fitting nucleotide substitution model was first determined using the automated tool available. Sequences of the IL genes were analyzed under three available models, single likelihood ancestor counting (SLAC), fixed-effect likelihood (FEL) and random effect likelihood (REL). The SLAC model is based on the reconstruction of ancestral sequences and counts the number of $d_S$ and $d_N$ changes at each codon position of the phylogeny. FEL estimates ratios of $d_N$ to $d_S$ changes for each site in an alignment. REL uses a flexible distribution and allows $d_S$ and $d_N$ to vary across sites independently.³⁶⁻³⁸

In CODEML two opposing models, M7 and M8, were compared using likelihood ratio tests. While M7 assumes that $\omega$ ratios are distributed among sites according to a beta distribution allowing codons to evolve neutrally or under purifying selection, M8 is an extension of the M7 model with an extra class of sites with an independent $\omega$ ratio freely estimated from the data allowing positive selection. M7 versus M8 were compared by computing twice the difference in the natural logs of the likelihoods ($2\Delta\ln L$). The value obtained was used in a likelihood ratio test along with the degrees of freedom (2) and allowed the rejection of the null model ($P < 0.05$). Amino acids detected as under positive selection in M8 were identified using the Bayes Empirical Bayes approach (BEB), with posterior probability > 95%. BEB is the preferred approach because it accounts sampling errors in the ML.³⁴,³⁵,³⁹⁻⁴¹ For each gene, a neighbor-joining tree was constructed in MEGA⁵⁴² as working topology with selected options p-distance as the substitution model and complete deletion to gaps/missing data treatment.

For a more conservative approach, and as used previously,⁴³,⁴⁴ only sites detected to be under positive selection in more than one ML method were considered.

# Results

Positive selection was found in 28 of the 46 ILs: IL-1 (A and B), IL-2, IL-4, IL-5, IL-6, IL-7, IL-8, IL-9, IL-10, IL-12 (A and B), IL-14, IL-15, IL-16, IL-17 (A and C), IL-18, IL-20, IL-21, IL-22, IL-25, IL-26, IL-27B, IL-31, IL-34 and IL-36 (A and G), as shown in Table 1. Codons under positive selection detected by
<table><caption>**Table 1.** Phylogenetic tests of positive selection.ª</caption><thead><tr><th>Gene</th><th>No. of aa residues</th><th>No. of species</th><th>In L M7</th><th>In L M8</th><th>2ΔlnL</th><th>PAML</th><th>SLAC<sup>b</sup></th><th>FEL<sup>c</sup></th><th>REL<sup>d</sup></th><th>Total no. of sites</th><th>% of sites</th></tr></thead><tbody><tr><td>IL-1 A</td><td>825</td><td>33</td><td>-10,987, 711</td><td>-10,939, 544</td><td>96,335*</td><td><u>104, 122, 132, 147, 148, 173, 174, 175, 178, 180, 190, 227, 230, 248, 258</u></td><td><u>173, 174</u></td><td>66, <u>173, 174, 180, 188, 215, 227, 230, 265</u></td><td>54, 99, <u>104, 116, 122, 125, 127, 129, 132, 139, 158, 162, 173, 174, 175, 179, 180, 215, 227, 230, 248, 250</u></td><td>11</td><td>1.333</td></tr><tr><td>IL-1 B</td><td>801</td><td>37</td><td>-9578, 964</td><td>-9565, 064</td><td>27,800*</td><td>79, 150, 205</td><td>8, <u>84</u></td><td>8, 36, 71, 84, 93, 96, 104, 129, 153</td><td>84</td><td>2</td><td>0.25</td></tr><tr><td>IL-2</td><td>513</td><td>47</td><td>-5429, 515</td><td>-5417, 492</td><td>24,045*</td><td><u>100, 101, 124</u></td><td>-</td><td>28, 39, 42, 58, 69</td><td>43, 46, 66, 68, 91, 98, <u>100, 101, 102, 115, 122, 124, 152</u></td><td>3</td><td>0.585</td></tr><tr><td>IL-4</td><td>465</td><td>44</td><td>-6436, 977</td><td>-6424, 860</td><td>24,235*</td><td>84, <u>105, 106, 120</u></td><td>84</td><td>5, 59, <u>105, 120, 128, 147</u></td><td><u>84, 105, 106, 120</u></td><td>4</td><td>0.86</td></tr><tr><td>IL-5</td><td>408</td><td>36</td><td>-4996, 583</td><td>-4994, 044</td><td>5078</td><td>-</td><td>-</td><td>4, 93, 96, <u>100</u></td><td><u>100</u></td><td>1</td><td>0.245</td></tr><tr><td>IL-6</td><td>657</td><td>42</td><td>-9824, 485</td><td>-9805, 607</td><td>37,755*</td><td><u>46, 55, 62, 65, 75, 120, 141, 145, 149, 158, 161, 162, 177</u></td><td>8, <u>145</u></td><td>8, <u>162, 168, 174, 193</u></td><td>30, 38, <u>46, 55, 56, 62, 75, 76, 96, 113, 120, 141, 145, 149, 153, 158, 161, 162, 168, 171, 173, 174, 177</u></td><td>15</td><td>2.283</td></tr><tr><td>IL-7</td><td>540</td><td>21</td><td>-3887, 223</td><td>-3861, 408</td><td>51,630*</td><td><u>47, 51, 54, 60, 92, 102, 113, 126, 153, 172, 174</u></td><td>92, <u>144</u></td><td>38, 43, <u>51, 53, 55, 92, 105, 106, 144</u></td><td>30, 38, 46, 47, <u>51, 53, 54, 55, 92, 101, 105, 106, 113, 143, 144</u></td><td>11</td><td>2.037</td></tr><tr><td>IL-8</td><td>297</td><td>33</td><td>-3256, 753</td><td>-3252, 791</td><td>7924*</td><td>71</td><td>-</td><td>22, <u>30</u></td><td>22, <u>30, 64, 71, 74</u></td><td>3</td><td>1010</td></tr><tr><td>IL-9</td><td>423</td><td>20</td><td>-4237, 433</td><td>-4217, 874</td><td>39,118*</td><td><u>5, 24, 25, 84, 86</u></td><td>84</td><td><u>79, 84, 106</u></td><td>5, 20, 22, 24, 25, 28, 31, 35, 49, <u>79, 82, 84, 86, 90, 98, 106, 131, 136, 140</u></td><td>7</td><td>1.655</td></tr><tr><td>IL-10</td><td>534</td><td>44</td><td>-6188, 435</td><td>-6184, 527</td><td>7815*</td><td>26</td><td><u>18, 63, 71, 121</u></td><td><u>18, 63, 71, 121</u></td><td><u>18, 20, 26, 35, 63, 68, 71, 121</u></td><td>5</td><td>0.936</td></tr><tr><td>IL-12A</td><td>669</td><td>27</td><td>-5853, 066</td><td>-5847, 159</td><td>11,813*</td><td><u>161</u></td><td>-</td><td>24, 35, 39, 91, 154</td><td>39, <u>161, 172</u></td><td>2</td><td>0.299</td></tr><tr><td>IL-12B</td><td>1032</td><td>33</td><td>-9173, 438</td><td>-9170, 760</td><td>5357</td><td>-</td><td><u>83, 147</u></td><td><u>83, 147</u></td><td>21, <u>83, 182, 185, 188, 216, 289, 302, 306</u></td><td>1</td><td>0.097</td></tr><tr><td>IL-14</td><td>1698</td><td>18</td><td>-8069, 857</td><td>-8069, 766</td><td>0.183</td><td>-</td><td>-</td><td><u>501</u></td><td><u>501</u></td><td>1</td><td>0.059</td></tr></tbody></table>

(continued)
<table><caption>Table I. Continued.</caption><thead><tr><th>Gene</th><th>No. of aa residues</th><th>No. of species</th><th>In L M7</th><th>In L M8</th><th>2Δln L</th><th>PAML</th><th>SLAC<sup>b</sup></th><th>FEL<sup>c</sup></th><th>REL<sup>d</sup></th><th>Total no. of sites</th><th>% of sites</th></tr></thead><tbody><tr><td>IL-15</td><td>483</td><td>29</td><td>-4003, 633</td><td>-3987, 300</td><td>32,666*</td><td>51, 105, 106, 108</td><td>-</td><td>4, 59, 116</td><td>51, 52, 55, 59, 62, 96, 105, 106, 108, 116, 152</td><td>6</td><td>1.242</td></tr><tr><td>IL-16</td><td>4095</td><td>17</td><td>-27,560, 331</td><td>-27,552, 931</td><td>14,800*</td><td>622</td><td>335</td><td>118, 188, 335, 633, 708, 734, 745, 801, 930, 1075</td><td>118, 188, 335, 705, 745, 801, 930, 1075</td><td>8</td><td>0.195</td></tr><tr><td>IL-17A</td><td>483</td><td>27</td><td>-5623, 567</td><td>-5618, 970</td><td>9193*</td><td>-</td><td>18, 98</td><td>18, 98</td><td>-</td><td>2</td><td>0.414</td></tr><tr><td>IL-17C</td><td>600</td><td>18</td><td>-4833, 872</td><td>-4829, 855</td><td>8034*</td><td>1</td><td>11, 17</td><td>11, 17</td><td>11, 17, 55</td><td>2</td><td>0.333</td></tr><tr><td>IL-18</td><td>588</td><td>24</td><td>-4474, 258</td><td>-4458, 874</td><td>30,768*</td><td>4, 5, 29, 45, 75, 105, 176, 192</td><td>180</td><td>6, 32, 37, 45, 159, 180</td><td>4, 6, 45, 75, 82, 105, 180, 192</td><td>7</td><td>1.19</td></tr><tr><td>IL-20</td><td>525</td><td>25</td><td>-4976, 365</td><td>-4974, 835</td><td>3060</td><td>-</td><td>-</td><td>56, 62</td><td>56, 62</td><td>2</td><td>0.381</td></tr><tr><td>IL-21</td><td>456</td><td>19</td><td>-3201, 764</td><td>-3191, 105</td><td>21,319*</td><td>98, 119</td><td>86, 112</td><td>86, 112</td><td>86, 98, 112, 119</td><td>4</td><td>0.877</td></tr><tr><td>IL-22</td><td>543</td><td>25</td><td>-5167, 215</td><td>-5159, 805</td><td>14,821*</td><td>86, 173</td><td>53</td><td>8, 53, 84</td><td>3, 8, 53, 65, 84</td><td>3</td><td>0.552</td></tr><tr><td>IL-25</td><td>510</td><td>20</td><td>-4452, 843</td><td>-4441, 182</td><td>23,323*</td><td>1, 2, 137</td><td>-</td><td>53, 137</td><td>36, 40, 55</td><td>1</td><td>0.196</td></tr><tr><td>IL-26</td><td>504</td><td>20</td><td>-3693, 866</td><td>-3687, 927</td><td>11,877*</td><td>2</td><td>2</td><td>2</td><td>-</td><td>1</td><td>0.198</td></tr><tr><td>IL-27B</td><td>708</td><td>16</td><td>-5613, 682</td><td>-5608, 707</td><td>9949*</td><td>-</td><td>114</td><td>114, 139</td><td>-</td><td>1</td><td>0.141</td></tr><tr><td>IL-31</td><td>489</td><td>14</td><td>-4093, 760</td><td>-4063, 715</td><td>60,091*</td><td>31, 33, 38, 67, 103, 133</td><td>-</td><td>30, 33, 38, 43, 56, 67, 141</td><td>31, 33, 38, 56, 66, 98, 133</td><td>6</td><td>1.227</td></tr><tr><td>IL-34</td><td>717</td><td>17</td><td>-4918, 901</td><td>-4917, 438</td><td>2926</td><td>-</td><td>134</td><td>134</td><td>134</td><td>1</td><td>0.139</td></tr><tr><td>IL-36A</td><td>477</td><td>18</td><td>-4805, 698</td><td>-4766, 449</td><td>78,496*</td><td>2, 4, 5, 8, 9, 12, 38, 133</td><td>-</td><td>8, 133</td><td>1, 2, 3, 5, 7, 8, 9, 12, 71, 133, 156</td><td>6</td><td>1.258</td></tr><tr><td>IL-36G</td><td>471</td><td>21</td><td>-4966, 855</td><td>-4916, 882</td><td>99,946*</td><td>2, 3, 5, 6, 104, 140</td><td>2, 6, 104</td><td>2, 6, 104, 141</td><td>2, 5, 6, 18, 26, 41, 52, 55, 69, 73, 91, 104, 128, 144, 148</td><td>4</td><td>0.849</td></tr></tbody></table>

<footnote><sup>a</sup>Codons identified by more than one ML method are underlined. Codons located in close vicinity of N-glycosylation sites are shaded and those codons related with disulfide bonds are boxed.</footnote>

<footnote><sup>b,c</sup>Codons with P-values < 0.05.</footnote>

<footnote><sup>d</sup>Codons with Bayes factor > 50.</footnote>

<footnote>*P < 0.05.</footnote>
at least two ML approaches ranged between 1 and 15, with IL-6 being the IL with the highest percentage of positively selected codons (2.28% with 15 codons under positive selection) followed by IL-7 (2.04% and 11 codons under positive selection). The smallest percentages of positively selected codons were 0.06% in IL-14 and 0.10% in IL-12B, each with 1 codon. No positively selected codons were detected in IL-13, IL-17B and F, IL-19, IL-23, IL-24, IL-27A and IL-29.

Alterations in coding regions of genes that lead to amino acids substitutions can induce changes in protein conformation. These changes may be conservative or radical, and may alter physiochemical properties of the proteins, such as charge and polarity.³³ From our physiochemical study of the amino acids under positive selection (see Appendix 2 in Supplementary material) we observed that for the majority of codons under positive selection changes in charge and polarity occur. For example, in IL-36G, amino acid positions 2 and 5 correspond to 13 and 12 amino acid possibilities, respectively, i.e. the different amino acids that can be found in all the species analyzed.

## Discussion

Mammalian genes involved in the immune response are among the most rapidly evolving genes²⁵,⁴⁵,⁴⁶ as they include proteins with biological activities designed to protect the host (antibacterial, antiviral, antifungal or antiparasitic). As positive selection is likely correlated with sites of important activity, an effort was made to try and verify if the positions under positive selection have any particular function that led to their evolution. For most ILs with sites under positive selection, a correlation between those sites and interaction with other molecules, such as receptors, pathogens and binding proteins (BP), was found. It is known that sites where interaction of the ILs with their receptors occurs are essential for their function.²⁶,⁴⁷ To verify a cause-effect relationship between ILs and their receptors, we extended our positive selection analyses to IL receptors (ILR), in particular focusing on sites that interact with ILs at sites detected as being under positive selection (Table 2). From all ILs where positive selection was detected, only in IL-4 and IL-18 was that correlation found. For IL-4, and by focusing on the residues under positive selection known to interact with IL-4 receptors, mainly IL-4Rα and IL-13Rα (see Figure 1), our results are in agreement with previous studies.⁴⁸⁻⁵¹ Indeed, of the IL-4 residues detected under positive selection, Glu84 is located between sites of interaction with IL-4Rα and IL-13Rα, while Arg105 and Phe106 are located in interaction sites of IL-4 with IL-4Rα. No codons under positive selection were detected at the sites of IL-4:γc interaction, as suggested by others.¹,²⁶ Of the two residues of IL-4Rα under positive selection, only Ser95 has a direct interaction with IL-4 (Asn113).

IL-4 positively selected sites Arg105 and Phe106 are located near sites where IL-4 interacts with IL-4Rα.²⁶,⁴⁹,⁵² For IL-13Rα, none of the seven residues detected as under positive selection interacted with Glu84 of IL-4, but the positively selected ILe254 from IL-13Rα interacted with IL-4. The sites detected are in accordance with those described by Koyanagi et al.,²⁶ but not with those detected by O'Connel et al.¹ One explanation for this discordance is the number of species used, which differs between studies. For IL-18, of the seven residues detected as under positive selection, only Leu45 is located near the site of interaction with IL-18Rα (Lys44) and IL-18BP (Glu42). The signatures of positive selection in these sites might also be involved in modulating signaling intensity.²⁶

In their evolution, pathogens have evolved an arsenal of immune-evading strategies, which include antagonists for the host immune-related proteins²⁸ that target critical sites for these proteins' functions. Thus, to escape from antagonists, those targeting sites need to continuously evolve so that the IL-ILR interaction can still be functional. Supporting this concept, sites under positive selection for IL-4 have been suggested to be associated with escaping from pathogen-encoded antagonistic proteins.²⁶ BP are polypeptides that have the ability to interact with corresponding proteins and neutralize them. In the ILs studied, there are many examples of such proteins. Indeed, IL-18BP has the ability to bind to IL-18 with high affinity preventing IL-18-IL-18R interaction and, therefore, neutralizing their biological activities. Some authors have demonstrated that IL-18BP prevents LPS-induced IFN-γ production.⁵³⁻⁵⁵ In some models, it has been described that blocking IL-18 results in a reduction in disease severity in injuries related to IL-18 increased production, namely injuries of heart, kidneys, liver, arthritis, etc.⁵⁶ By studying IL-18, Leu45 which was detected as under positive selection is located close to the site of interaction of IL-18 (Glu42) with IL-18BP. In addition, when we studied IL-18BP, we detected that Cys131 is under positive selection next to the site of interaction of IL-18BP (Lys130) with IL-18.

Nevertheless, for some ILs, the evidence of positive selection is not so intelligible. When pathogens or stress are encountered, the initial response is via the innate immune system with activation of the inflammasome. Interestingly, when examining Table 1 two of the genes with the largest number of positively selected residues are IL-1A and IL-6. IL-1A is the initial cytokine produced upon inflammasome activation, which directly leads to increased expression of IL-6.⁹ Given their critical importance in the early activation of the immune response, it is tempting to speculate that there has been selective pressure to maintain (rather than to change) structure and function of these molecules from their early precursors. However, knowledge on the functional role of IL-1A and IL-6 derives mostly from
<table><caption>**Table 2.** Phylogenetic tests of positive selection for receptors and binding proteins.ª</caption><thead><tr><th>Gene</th><th>No. of aa residues</th><th>No. of Species</th><th>In L M7</th><th>In L M8</th><th>2Δln L</th><th>PAML</th><th>SLAC<sup>b</sup></th><th>FEL<sup>c</sup></th><th>REL<sup>d</sup></th><th>Total no. of sites</th><th>% of sites</th></tr></thead><tbody><tr><td>IL-4Rα</td><td>2529</td><td>10</td><td>-13,357, 029</td><td>-13346,832</td><td>20,393*</td><td><u>95</u></td><td>-</td><td>471, 479, 733, 764</td><td>67, 68, 93, 95, 529, 590, 677, 689, 764</td><td>2</td><td>0.079</td></tr><tr><td>IL-13Rα</td><td>1296</td><td>19</td><td>-9056, 699</td><td>-9038, 203</td><td>36,992*</td><td><u>254, 276</u></td><td><u>254</u></td><td><u>83, 196, 225, 237, 254, 259</u></td><td>16, 19, 25, 31, 83, 171, 196, 225, 235, 237, 239, 242, 254, 259, 272, 274, 276, 289, 334, 335</td><td>7</td><td>0.540</td></tr><tr><td>IL-18BP</td><td>573</td><td>21</td><td>-5034, 333</td><td>-5016, 838</td><td>34,991*</td><td><u>54, 60, 126, 131, 156, 176</u></td><td><u>11</u></td><td><u>9, 11, 24, 60, 61</u></td><td><u>11, 60, 63, 126, 131, 153, 156</u></td><td>5</td><td>0.873</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Positively selected sites in the three-dimensional structures of IL-4-IL-4Rα-IL-13Rα (marked in yellow). IL-4 appears in pink, IL-4Rα in blue and IL-13Rα in brown.</figcaption></figure>

studies in humans with almost nothing known for other mammals. It is possible that, in other species, these ILs might have evolved for other tasks, leading to the observed divergence in sequence across species.

Codons under positive selection were also detected within sites of N-glycosylation or where disulfide bonds exist. Glycosylation is considered to be important for protein folding, oligomerization, intrinsic stability, solubility, capacity to diffuse throughout the organism, interaction with cell surface receptors and subsequent biological activity.⁵⁷⁻⁵⁹ This modification is an effective way to generate diversity and modulate protein properties due to inherent structural variations of the glycans.⁵⁷,⁵⁹,⁶⁰ For IL-1A, IL-9, IL-22, IL-25 and IL-31, residues under positive selection were found in close proximity to sites where N-glycosylation occurs (Table 1). Disulfide bonds play an important role in folding, stability and function of the protein, being, according to some authors,⁶¹,⁶² associated with functional differentiation of the proteins. These bonds are thought to be well conserved in proteins.⁶¹ Positive selection was detected at sites where disulfide bonds exist for IL-2 and IL-17A (Table 1) and may compromise their structure, stability and function. The changes in charge and polarity observed at the positively selected codons are likely to impose changes in the protein conformation with consequences at the protein function.

## Conclusions

Positive selection may function to maintain a host response to pathogens, diseases and environmental

<footnote>ªCodons identified by more than one ML method are underlined.</footnote>

<footnote>ᵇCodons with P-values < 0.05.</footnote>

<footnote>ᶜCodons with Bayes factor > 50.</footnote>

<footnote>*P < 0.05.</footnote>
conditions. Here, we have detected positive selection in 28 ILs with some of the identified codons associated with critical functions of these proteins. Several reasons might underlie our results, including the specific biological functions in which IL-s are involved and the maintenance of critical structural features of the ILs. However, the limited knowledge of the role of specific codons in IL functions for most mammalian species hampers the complete understanding of our observations. Further functional and structural studies by using mutagenesis and crystallographic approaches should be performed for full comprehension of the role of the observed variation in mammalian ILs.

## Funding

The Portuguese Foundation for Science and Technology (FCT-Portugal) supported the doctoral fellowship of Fabiana Neves (SFRH/BD/81916/2011), the post-doctoral fellowships of Joana Abrantes (SFRH/BPD/73512/2010) and Pedro J. Esteves (SPRH/BPD/27021/2006). This study was performed under the project conceded by FCT (PTDC/BIA-BEC/103158/2008).

## References

1. O'Connell MJ and McInerney JO. Gamma chain receptor interleukins: evidence for positive selection driving the evolution of cell-to-cell communicators in the mammalian immune system. *J Mol Evol* 2005; 61: 608-619.

2. Zelus D, Robinson-Rechavi M, Delacre M, et al. Fast evolution of interleukin-2 in mammals and positive selection in ruminants. *J Mol Evol* 2000; 51: 234-244.

3. Kaiser P, Rothwell L, Avery S and Balu S. Evolution of the interleukins. *Develop Comp Immunol* 2004; 28: 375-394.

4. Zhang J and Nei M. Positive selection in the evolution of mammalian interleukin-2 genes. *Mol Biol Evol* 2000; 17: 1413-1416.

5. Brocker C, Thompson D, Matsumoto A, et al. Evolutionary divergence and functions of the human interleukin (IL) gene family. *Hum Genomics* 2010; 5: 30-55.

6. Poon AF, Frost SD and Pond SL. Detecting signatures of selection from DNA sequences using Datamonkey. *Methods Mol Biol* 2009; 537: 163-183.

7. Schrader JW. Interleukin is as interleukin does. *J Immunol Methods* 2003; 276: 1-3.

8. Beadling C and Slifka MK. Regulation of innate and adaptive immune responses by the related cytokines IL-12, IL-23, and IL-27. *Arch Immunol Ther Exp (Warsz)* 2006; 54: 15-24.

9. Akdis M, Burgler S, Crameri R, et al. Interleukins, from 1 to 37, and interferon-gamma: receptors, functions, and roles in diseases. *J Allergy Clin Immunol* 2011; 127: 701-721 e1-e70.

10. Ozaki K and Leonard WJ. Cytokine and cytokine receptor pleiotropy and redundancy. *J Biol Chem* 2002; 277: 29355-29358.

11. Afzal N, Tahir R and Jahan S. Cytokines: an ever expanding area. *Biol Biomed Rep* 2012; 2: 37-43.

12. Heinrich PC, Behrmann I, Haan S, et al. Principles of interleukin (IL)-6-type cytokine signalling and its regulation. *Biochem J* 2003; 374: 1-20.

13. Ishihara K and Hirano T. Molecular basis of the cell specificity of cytokine action. *Biochim Biophys Acta* 2002; 1592: 281-296.

14. Borish LC and Steinke JW. 2. Cytokines and chemokines. *J Allergy Clin Immunol* 2003; 111: S460-S475.

15. Commins SP, Borish L and Steinke JW. Immunologic messenger molecules: cytokines, interferons, and chemokines. *J Allergy Clin Immunol* 2010; 125: S53-S72.

16. Steinke JW and Borish L. 3. Cytokines and chemokines. *J Allergy Clin Immunol* 2006; 117: S441-S445.

17. Sabat R. IL-10 family of cytokines. *Cytokine Growth Factor Rev* 2010; 21: 315-324.

18. Krause CD and Pestka S. Evolution of the Class 2 cytokines and receptors, and discovery of new friends and relatives. *Pharmacol Ther* 2005; 106: 299-346.

19. Leonard WJ and Lin JX. Cytokine receptor signaling pathways. *J Clin Immunol* 2000; 105: 877-888.

20. Secombes CJ, Wang T and Bird S. The interleukins of fish. *Develop Immunol* 2011; 35: 1336-1345.

21. Fumagalli M, Pozzoli U, Cagliani R, et al. Parasites represent a major selective force for interleukins genes and shape the genetic predisposition to autoimmune conditions. *J Exp Med* 2009; 206: 1395-1408.

22. Commins S, Steinke JW and Borish L. The extended IL-10 superfamily: IL-10, IL-19, IL-20, IL-22, IL-24, IL-26, IL-28, and IL-29. *J Allergy Clin Immunol* 2008; 121: 1108-1111.

23. Feghali CA and Wright TM. Cytokines in acute and chronic inflammation. *Front Biosci* 1997; 2: d12-d26.

24. Haddad JJ. Cytokines and related receptor-mediated signaling pathways. *Biochem Biophys Res Commun* 2002; 297: 700-713.

25. Barreiro LB and Quintana-Murci L. From evolutionary genetics to human immunology: how selection shapes host defence genes. *Nat Rev Genet* 2010; 11: 17-30.

26. Koyanagi M, Kerns JA, Chung L, et al. Diversifying selection and functional analysis of interleukin-4 suggests antagonism-driven evolution at receptor-binding interfaces. *BMC Evol Biol* 2010; 10: 223.

27. Makalowski W and Boguski MS. Synonymous and nonsynonymous substitution distances are correlated in mouse and rat genes. *J Mol Evol* 1998; 47: 119-121.

28. PIL-lai MR and Bix M. Evolution of IL4 and pathogen antagonism. *Growth Factors* 2011; 29: 153-160.

29. Manry J, Laval G, Patin E, et al. Evolutionary genetic dissection of human interferons. *J Exp Med* 2011; 208: 2747-2759.

30. Tindall EA and Hayes VM. Comprehensive sequence analysis of the human IL23A gene defines new variation content and high rate of evolutionary conservation. *DNA Res* 2010; 17: 117-122.

31. Zhou G, Zhai Y, Dong X, et al. Haplotype structure and evidence for positive selection at the human IL13 locus. *Mol Biol Evol* 2004; 21: 29-35.

32. Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. *Nucl Acids Symp Ser* 1999; 41: 95-98.

33. Majewski J and Ott J. Amino acid substitutions in the human genome: evolutionary implications of single nucleotide polymorphisms. *Gene* 2003; 305: 167-173.

34. Yang Z. PAML: a program package for phylogenetic analysis by maximum likelihood. *Comput Appl Biosci* 1997; 13: 555-556.

35. Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. *Mol Biol Evol* 2007; 24: 1586-1591.

36. Kosakovsky Pond SL and Frost SD. Not so different after all: a comparison of methods for detecting amino acid sites under selection. *Mol Biol Evol* 2005; 22: 1208-1222.

37. Pond SL and Frost SD. Datamonkey: rapid detection of selective pressure on individual sites of codon alignments. *Bioinformatics* 2005; 21: 2531-2533.

38. Poon AF, Frost SD and Pond SL. Detecting signatures of selection from DNA sequences using Datamonkey. *Methods Mol Biol* 2009; 537: 163-183.

39. Bielawski JP and Yang Z. Maximum likelihood methods for detecting adaptive evolution after gene duplication. *J Struct Funct Genomics* 2003; 3: 201-212.

40. Yang Z. Inference of selection from multiple species alignments. *Curr Opin Genet Dev* 2002; 12: 688-694.
41. Yang Z, Nielsen R, Goldman N and Pedersen AM. Codon-substitution models for heterogeneous selection pressure at amino acid sites. *Genetics* 2000; 155: 431–449.

42. Tamura K, Peterson D, Peterson N, et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. *Mol Biol Evol* 2011; 28: 2731–2739.

43. Areal H, Abrantes J and Esteves PJ. Signatures of positive selection in Toll-like receptor (TLR) genes in mammals. *BMC Evol Biol* 2011; 11: 368.

44. Wlasiuk G and Nachman MW. Adaptation and constraint at Toll-like receptors in primates. *Mol Biol Evol* 2010; 27: 2172–2186.

45. Kosiol C, Vinar T, da Fonseca RR, et al. Patterns of positive selection in six Mammalian genomes. *PLoS Genet* 2008; 4: e1000144.

46. Ferrer-Admetlla A, Bosch E, Sikora M, et al. Balancing selection is the main force shaping the evolution of innate immunity genes. *J Immunol* 2008; 181: 1315–1322.

47. Smith KA. The structure of IL2 bound to the three chains of the IL2 receptor and how signaling occurs. *Med Immunol* 2006; 5: 3.

48. Hage T, Sebald W and Reinemer P. Crystal structure of the interleukin-4/receptor alpha chain complex reveals a mosaic binding interface. *Cell* 1999; 97: 271–281.

49. LaPorte SL, Juo ZS, Vaclavikova J, et al. Molecular and structural basis of cytokine receptor pleiotropy in the interleukin-4/13 system. *Cell* 2008; 132: 259–272.

50. Zarlenga DS, Dawson H, Kringel H, et al. Molecular cloning of the swine IL-4 receptor alpha and IL-13 receptor 1-chains: effects of experimental *Toxoplasma gondii*, *Ascaris suum* and *Trichuris suis* infections on tissue mRNA levels. *Vet Immunol Immunopathol* 2004; 101: 223–234.

51. Zhang JL, Simeonowa I, Wang Y and Sebald W. The high-affinity interaction of human IL-4 and the receptor alpha chain is

constituted by two independent binding clusters. *J Mol Biol* 2002; 315: 399–407.

52. Mueller TD, Zhang JL, Sebald W and Duschl A. Structure, binding, and antagonists in the IL-4/IL-13 receptor system. *Biochim Biophys Acta* 2002; 1592: 237–250.

53. Dinarello CA. Targeting interleukin 18 with interleukin 18 binding protein. *Ann Rheum Dis* 2000; 59(Suppl. 1): i17–i20.

54. He Z, Lu L, Altmann C, et al. Interleukin-18 binding protein transgenic mice are protected against ischemic acute kidney injury. *Am J Physiol Renal Physiol* 2008; 295: F1414–F1421.

55. Mallat Z, Corbaz A, Scoazec A, et al. Interleukin-18/interleukin-18 binding protein signaling modulates atherosclerotic lesion development and stability. *Circ Res* 2001; 89: E41–E45.

56. Dinarello CA. Interleukin 1 and interleukin 18 as mediators of inflammation and the aging process. *Am J Clin Nutr* 2006; 83: 447S–455S.

57. Shental-Bechor D and Levy Y. Effect of glycosylation on protein folding: a close look at thermodynamic stabilization. *Proc Natl Acad Sci USA* 2008; 105: 8256–8261.

58. Szabo TG, Palotai R, Antal P, et al. Critical role of glycosylation in determining the length and structure of T cell epitopes. *Immunome Res* 2009; 5: 4.

59. Waetzig GH, Chalaris A, Rosenstiel P, et al. N-linked glycosylation is essential for the stability but not the signaling function of the interleukin-6 signal transducer glycoprotein 130. *J Biol Chem* 2010; 285: 1781–1789.

60. Chamorey AL, Magne N, Pivot X and MIL-ano G. Impact of glycosylation on the effect of cytokines. A special focus on oncology. *Eur Cytokine Netw* 2002; 13: 154–160.

61. Li XQ, Zhang T and Donnelly D. Selective loss of cysteine residues and disulphide bonds in a potato proteinase inhibitor II family. *PLoS One* 2011; 6: e18615.

62. Fass D. Disulfide bonding in protein biophysics. *Annu Rev Biophys* 2012; 41: 63–79.