Marie-Ka Tilak, Constance Colombier, Juan Carlos Illera, Martim Melo, Claire Doutrelant, Claire Loiseau, Mathilde Barthe, Thibault Leroy, Benoit Nabholz, Rita Covas, Institut des Sciences de l'Evolution de Montpellier (UMR ISEM), Centre de Coopération Internationale en Recherche Agronomique pour le Développement (Cirad)-École Pratique des Hautes Études (EPHE), Université Paris sciences et lettres (PSL)-Université Paris sciences et lettres (PSL)-Institut de recherche pour le développement [IRD] : UR226-Centre National de la Recherche Scientifique (CNRS)-Université de Montpellier (UM), Centre d’Ecologie Fonctionnelle et Evolutive (CEFE), Université Paul-Valéry - Montpellier 3 (UPVM)-École Pratique des Hautes Études (EPHE), Université Paris sciences et lettres (PSL)-Université Paris sciences et lettres (PSL)-Centre National de la Recherche Scientifique (CNRS)-Institut de Recherche pour le Développement (IRD [France-Sud])-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE)-Institut Agro Montpellier, Institut national d'enseignement supérieur pour l'agriculture, l'alimentation et l'environnement (Institut Agro)-Institut national d'enseignement supérieur pour l'agriculture, l'alimentation et l'environnement (Institut Agro)-Université de Montpellier (UM), Research Centre in Biodiversity and Genetic Resources (CIBIO-InBio), Program in Genomics, Biodiversity and Land Planning (BIOPOLIS), FitzPatrick Institute of African Ornithology, University of Cape Town-DST-NRF Centre of Excellence, Centro de Investigacão em Biodiversidade e Recursos Genéticos (CIBIO), Universidade do Porto = University of Porto, Consejo Superior de Investigaciones Científicas [Madrid] (CSIC), Institut de Recherche en Horticulture et Semences (IRHS), Université d'Angers (UA)-Institut National de Recherche pour l’Agriculture, l’Alimentation et l’Environnement (INRAE)-Institut Agro Rennes Angers, Institut national d'enseignement supérieur pour l'agriculture, l'alimentation et l'environnement (Institut Agro)-Institut national d'enseignement supérieur pour l'agriculture, l'alimentation et l'environnement (Institut Agro), CIBIO-InBIO - Universidade do Porto, Institut Universitaire de France (IUF), Ministère de l'Education nationale, de l’Enseignement supérieur et de la Recherche (M.E.N.E.S.R.), Labex CeMEB (project ISLAND IMMUNITY), the National Geographic Society (Grant/Award Number:W251-12), the British Ecological Society (Grant/Award Number: 369/4558), the Portuguese Foundation for Science and Technology under the PTDC/BIA-EVL/29390/2017 'DEEP' Research Project, CIBIO (UIDB/50027/2021), the European Regional Development Fund (Ref.: PGC2018-097575-B-I00, European Union’s Horizon 2020 research and innovation programme under grant agreement 854248., and ANR-14-CE02-0002,BirdIslandGenomic,Influence de la taille des populations sur l'évolution du génome: le cas des Oiseaux endémiques de l'île de La Réunion(2014)
Shared ecological conditions encountered by species that colonize islands often lead to the evolution of convergent phenotypes, commonly referred to as “island syndrome”. Reduced immune functions have been previously proposed to be part of this syndrome, as a consequence of the reduced diversity of pathogens on island ecosystems. According to this hypothesis, immune genes are expected to exhibit genomic signatures of relaxed selection pressure in island species. In this study, we used comparative genomic methods to study immune genes in island species (N = 20) and their mainland relatives (N = 14). We gathered public data as well as generated new data on innate (TLR: Toll-Like Receptors, BD: Beta Defensins) and acquired immune genes (MHC: Major Histocompatibility Complex classes I and II), but also on hundreds of genes with various immune functions. As a control, we used a set of 97 genes, not known to be involved in immune functions based on the literature, to account for the increased drift effects of the lower effective population sizes in island species. We used synonymous and non-synonymous variants to estimate the selection pressure acting on immune genes. We found that BDs and TLRs have higher ratios of non-synonymous over synonymous polymorphisms (Pn/Ps) than randomly selected control genes, suggesting that they evolve under a different selection regime. However, simulations show that this is unlikely to be explained by ongoing positive selection or balancing selection. For the MHC genes, which evolve under balancing selection, we used simulations to estimate the impact of population size variation. We found a significant effect of drift on immune genes of island species leading to a reduction in genetic diversity and efficacy of selection. However, the intensity of relaxed selection was not significantly different from control genes, except for MHC class II genes. These genes exhibit a significantly higher level of non-synonymous loss of polymorphism than expected assuming only drift and evolution under frequency dependent selection, possibly due to a reduction of extracellular parasite communities on islands. Overall, our results showed that demographic effects lead to a decrease in the immune functions of island species, but the relaxed selection that is expected to be caused by a reduced parasite pressure may only occur in some categories of immune genes.