Search

The study of insular populations was key in the development of evolutionary theory. The successful colonisation of an island depends on the geographic context, and specific characteristics of the organism and the island, but also on stochastic processes. As a result, apparently identical islands may harbour populations with contrasting histories. Here, we use whole genome sequences of 65 barn owls to investigate the patterns of inbreeding and genetic diversity of insular populations in the eastern Mediterranean Sea. We focus on Crete and Cyprus, islands with similar size, climate and distance to mainland, that provide natural replicates for a comparative analysis of the impacts of microevolutionary processes on isolated populations. We show that barn owl populations from each island have a separate origin, Crete being genetically more similar to other Greek islands and mainland Greece, and Cyprus more similar to the Levant. Further, our data show that their respective demographic histories following colonisation were also distinct. On the one hand, Crete harbours a small population and maintains very low levels of gene flow with neighbouring populations. This has resulted in low genetic diversity, strong genetic drift, increased relatedness in the population and remote inbreeding. Cyprus, on the other hand, appears to maintain enough gene flow with the mainland to avoid such an outcome. Our study provides a comparative population genomic analysis of the effects of neutral processes on a classical island-mainland model system. It provides empirical evidence for the role of stochastic processes in determining the fate of diverging isolated populations.

The climate fluctuations of the Quaternary shaped the movement of species in and out of glacial refugia. In Europe, the majority of species followed one of the described traditional postglacial recolonization routes from the southern peninsulas towards the north. Like most organisms, barn owls are assumed to have colonized the British Isles by crossing over Doggerland, a land bridge that connected Britain to northern Europe. However, while they are dark rufous in northern Europe, barn owls in the British Isles are conspicuously white, a contrast that could suggest selective forces are at play on the islands. However, analysis of known candidate genes involved in colouration found no signature of selection. Instead, using whole genome sequences and species distribution modelling, we found that owls colonised the British Isles soon after the last glaciation, directly from a white coloured refugium in the Iberian Peninsula, before colonising northern Europe. They would have followed a yet unknown post-glacial colonization route to the Isles over a westwards path of suitable habitat in now submerged land in the Bay of Biscay, thus not crossing Doggerland. As such, they inherited the white colour of their Iberian founders and maintained it through low gene flow with the mainland that prevents the import of rufous alleles. Thus, we contend that neutral processes likely explain this contrasting white colour compared to continental owls. With the barn owl being a top predator, we expect future research will show this unanticipated route was used by other species from its paleo community.

Examination of genetic polymorphisms in outbred wild-living species provides insights into the evolution of complex systems. In higher vertebrates, the proopiomelanocortin (POMC) precursor gives rise to α-, β-, and γ-melanocyte-stimulating hormones (MSH), which are involved in numerous physiological aspects. Genetic defects in POMC are linked to metabolic disorders in humans and animals. In the present study, we undertook an evolutionary genetic approach complemented with biochemistry to investigate the functional consequences of genetic polymorphisms in the POMC system of free-living outbred barn owl species (family Tytonidae) at the molecular level. Our phylogenetic studies revealed a striking correlation between a loss-of-function H9P mutation in the β-MSH receptor-binding motif and an extension of a poly-serine stretch in γ3-MSH to ≥7 residues that arose in the barn owl group 6-8 MYA ago. We found that extension of the poly-serine stretches in the γ-MSH locus affects POMC precursor processing, increasing γ3-MSH production at the expense of γ2-MSH and resulting in an overall reduction of γ-MSH signaling, which may be part of a negative feedback mechanism. Extension of the γ3-MSH poly-serine stretches ≥7 further markedly increases peptide hormone stability in plasma, which is conserved in humans, and is likely relevant to its endocrine function. In sum, our phylogenetic analysis of POMC in wild living owls uncovered a H9P β-MSH mutation subsequent to serine extension in γ3-MSH to 7 residues, which was then followed by further serine extension. The linked MSH mutations highlight the genetic plasticity enabled by the modular design of the POMC gene.

New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single‐Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations.

Knowledge of how and why secondary sexual characters are associated with sex hormones is important to understand their signalling function. Such a link can occur if i) testosterone participates in the elaboration of sex-traits, ii) the display of an ornament triggers behavioural response in conspecifics that induce a rise in testosterone, or iii) genes implicated in the elaboration of a sex-trait pleiotropically regulate testosterone physiology. To evaluate the origin of the co-variation between melanism and testosterone, we measured this hormone and the expression of enzymes involved in its metabolism in feathers of barn owl (Tyto alba) nestlings at the time of melanogenesis and in adults outside the period of melanogenesis. Male nestlings displaying smaller black feather spots had higher levels of circulating testosterone, potentially suggesting that testosterone could block the production of eumelanin pigments, or that genes involved in the production of small spots pleiotropically regulate testosterone production. In contrast, the enzyme 5α-reductase, that metabolizes testosterone to DHT, was more expressed in feathers of reddish-brown than light-reddish nestlings. This is consistent with the hypothesis that testosterone might be involved in the expression of reddish-brown pheomelanic pigments. In breeding adults, male barn owls displaying smaller black spots had higher levels of circulating testosterone, whereas in females the opposite result was detected during the rearing period, but not during incubation. The observed sex- and age-specific co-variations between black spottiness and testosterone in nestling and adult barn owls may not result from testosterone-dependent melanogenesis, but from melanogenic genes pleiotropically regulating testosterone, or from colour-specific life history strategies that influence testosterone levels.

New genomic tools open doors to study ecology, evolution, and population genomics of wild animals. For the Barn owl species complex, a cosmopolitan nocturnal raptor, a very fragmented draft genome was assembled for the American species (Tyto furcata pratincola) (Jarvis et al. 2014). To improve the genome, we assembled de novo Illumina and Pacific Biosciences (PacBio) long reads sequences of its European counterpart (Tyto alba alba). This genome assembly of 1.219 Gbp comprises 21,509 scaffolds and results in a N50 of 4,615,526 bp. BUSCO (Universal Single‐Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). By mapping the reads of the female American barn owl to the male European barn owl reads, we detected several structural variants and identified 70 Mbp of the Z chromosome. The barn owl scaffolds were further mapped to the chromosomes of the zebra finch. In addition, the completeness of the European barn owl genome is demonstrated with 94 of 128 proteins missing in the chicken genome retrieved in the European barn owl transcripts. This improved genome will help future barn owl population genomic investigations., We assemble the genome of the European barn owl (Tyto alba alba) of 1.219 Gbp by mixing Illumina (117×) and Pacific Biosciences (PacBio) long reads (12×) technologies. BUSCO (Universal Single‐Copy Orthologs) analysis revealed an assembly completeness of 94.8% with only 1.8% of the genes missing out of 4,915 avian orthologs searched, a proportion similar to that found in the genomes of the zebra finch (Taeniopygia guttata) or the collared flycatcher (Ficedula albicollis). In addition, we identify the Z chromosome' scaffolds, 73% of proteins thought to be missing in the chicken.

The melanocortin-1 receptor (MC1R) gene influences coloration by altering the expression of genes acting downstream in the melanin synthesis. MC1R belongs to the melanocortin system, a genetic network coding for the ligands that regulate MC1R and other melanocortin receptors controlling different physiological and behavioural traits. The impact of MC1R variants on these regulatory melanocortin genes was never considered, even though MC1R mutations could alter the influence of these genes on coloration (e.g. by decreasing MC1R response to melanocortin ligands). Using barn owl growing feathers, we investigated the differences between MC1R genotypes in the (co)expression of six melanocortin and nine melanogenic-related genes and in the association between melanocortin gene expression and phenotype (feather pheomelanin content). Compared to the MC1R rufous allele, responsible for reddish coloration, the white allele was not only associated with an expected lower expression of melanogenic-related genes (TYR, TYRP1, OCA2, SLC45A2, KIT, DCT) but also with a lower MC1R expression and a higher expression of ASIP, the MC1R antagonist. More importantly, the expression of PCSK2, responsible for the maturation of the MC1R agonist, α-melanocyte-stimulating hormone, was positively related to pheomelanin content in MC1R white homozygotes but not in individuals carrying the MC1R rufous allele. These findings indicate that MC1R mutations not only alter the expression of melanogenic-related genes but also the association between coloration and the expression of melanocortin genes upstream of MC1R. This suggests that MC1R mutations can modulate the regulation of coloration by the pleiotropic melanocortin genes, potentially decoupling the often-observed associations between coloration and other phenotypes.

Glucocorticoid hormones are important intermediates between an organism and its environment. They enable an organism to adjust its behavioural and physiological processes in response to environmental changes by binding to mineralocorticoid receptors (MR) and glucocorticoid receptors (GR) expressed in many tissues, including the integument. The regulation of glucocorticoids co-varies with melanin-based colouration in numerous species, an association that might result from pleiotropic effects of genes in the melanocortin system and evolve within a signalling context. Most studies have focused on the circulating levels of glucocorticoids disregarding the receptors that mediate their action, and that might partly account for the covariation between the regulation of stress and melanin-based colouration. We investigated the association of the expression levels of GR and MR genes with melanin-based colouration in the growing feathers of nestling barn owls (Tyto alba). We also explored the association between GR and MR expression levels and the expression of genes related to the melanocortin system and melanogenesis to better understand the origin of the link between the expression of receptors to which corticosterone binds and melanin-based colouration. Nestling barn owls displaying larger eumelanic black feather spots expressed GR and MR at lower levels than smaller-spotted individuals. However, we found that the expression of the GR and MR genes was positively rather than negatively correlated with the expression of genes involved in the deposition of melanin pigments at the time we sampled the nestlings. This provides mixed evidence of the association between melanin-based traits and MR and GR gene expression. The finding that the expression of GR and MR was positively associated with the expression of the PCSK2 gene (encoding one of the protein convertase responsible for the production of hormone peptide ACTH and α-MSH) suggests that the melanocortin system may be implicated in the establishment of the covariation between melanin-based colour and the expression of receptors to which glucocorticoids bind. However, further studies investigating the expression of melanin-based traits with stress-related endpoints at different time points of feather development will be necessary to understand better the proximate mechanism linking melanin-based traits with stress.

The mean phenotypic effects of a discovered variant help to predict major aspects of the evolution and inheritance of a phenotype. However, differences in the phenotypic variance associated to distinct genotypes are often overlooked despite being suggestive of processes that largely influence phenotypic evolution, such as interactions between the genotypes with the environment or the genetic background. We present empirical evidence for a mutation at the melanocortin-1-receptor gene, a major vertebrate coloration gene, affecting phenotypic variance in the barn owl, Tyto alba. The white MC1R allele, which associates with whiter plumage coloration, also associates with a pronounced phenotypic and additive genetic variance for distinct color traits. Contrarily, the rufous allele, associated with a rufous coloration, relates to a lower phenotypic and additive genetic variance, suggesting that this allele may be epistatic over other color loci. Variance differences between genotypes entailed differences in the strength of phenotypic and genetic associations between color traits, suggesting that differences in variance also alter the level of integration between traits. This study highlights that addressing variance differences of genotypes in wild populations provides interesting new insights into the evolutionary mechanisms and the genetic architecture underlying the phenotype.

The melanocortin system is implicated in the expression of many phenotypic traits. Activation of the melanocortin MC(1) receptor by melanocortin hormones induces the production of brown/black eumelanic pigments, while activation of the four other melanocortin receptors affects other physiological and behavioural functions including stress response, energy homeostasis, anti-inflammatory and sexual activity, aggressiveness and resistance to oxidative stress. We recently proposed the hypothesis that some melanocortin-physiological and -behavioural traits are correlated within individuals. This hypothesis predicts that the degree of eumelanin production may, in some cases, be associated with the regulation of glucocorticoids, immunity, resistance to oxidative stress, energy homeostasis, sexual activity, and aggressiveness. A review of the zoological literature and detailed experimental studies in a free-living population of barn owls (Tyto alba) showed that indeed melanic coloration is often correlated with the predicted physiological and behavioural traits. Support for predictions of the hypothesis that covariations between coloration and other phenotypic traits stem from pleiotropic effects of the melanocortin system raises a number of theoretical and empirical issues from evolutionary and pharmacological point of views.

Sexual selection theory predicts that males advertise quality by displaying extravagant ornaments. By contrast, whether phenotypic variation in females has a signalling function remains an open question. Here, to our knowledge, we provide the first evidence that a female plumage trait can signal fluctuating asymmetry in the offspring. We experimentally demonstrate in wild barn owls (Tyto alba) that the extent to which females display black spots on their plumage does not only signal offspring parasite resistance as shown in a previous study but also developmental homeostasis in the offspring. A greater number of spotted females produced offspring that had more symmetrical feathers during the period of growth. Males, that pair non-randomly with respect to female plumage spottiness therefore appear to gain substantial benefits by mating with heavily spotted females. Genetic variation in plumage spottiness is nevertheless maintained as the covariation between offspring body mass and mother plumage spottiness varies annually depending on environmental conditions.

Variants of the melanocortin-1 receptor (MC1R) gene result in abrupt, naturally selected colour morphs. These genetic variants may differentially affect sexual dimorphism if one morph is naturally selected in the two sexes but another morph is naturally or sexually selected only in one of the two sexes (e.g. to confer camouflage in reproductive females or confer mating advantage in males). Therefore, the balance between natural and sexual selections can differ between MC1R variants, as suggest studies showing interspecific correlations between sexual dimorphism and the rate of nonsynonymous vs. synonymous amino acid substitutions at the MC1R. Surprisingly, how MC1R is related to within-species sexual dimorphism, and thereby to sex-specific selection, has not yet been investigated. We tackled this issue in the barn owl (Tyto alba), a species showing pronounced variation in the degree of reddish pheomelanin-based coloration and in the number and size of black feather spots. We found that a valine (V)-to-isoleucine (I) substitution at position 126 explains up to 30% of the variation in the three melanin-based colour traits and in feather melanin content. Interestingly, MC1R genotypes also differed in the degree of sexual colour dimorphism, with individuals homozygous for the II MC1R variant being 2 times redder and 2.5 times less sexually dimorphic than homozygous individuals for the VV MC1R variant. These findings support that MC1R interacts with the expression of sexual dimorphism and suggest that a gene with major phenotypic effects and weakly influenced by variation in body condition can participate in sex-specific selection processes.

The adaptive function of melanin-based coloration is a long-standing debate. A recent genetic model suggested that pleiotropy could account for covariations between pigmentation, behaviour, morphology, physiology and life history traits. We explored whether the expression levels of genes belonging to the melanocortin system (MC1R, POMC, PC1/3, PC2 and the antagonist ASIP), which have many pleiotropic effects, are associated with melanogenesis (through variation in the expression of the genes MITF, SLC7A11, TYR, TYRP1) and in turn melanin-based coloration. We considered the tawny owl (Strix aluco) because individuals vary continuously from light to dark reddish, and thus, colour variation is likely to stem from differences in the levels of gene expression. We measured gene expression in feather bases collected in nestlings at the time of melanin production. As expected, the melanocortin system was associated with the expression of melanogenic genes and pigmentation. Offspring of darker reddish fathers expressed PC1/3 to lower levels but tended to express PC2 to higher levels. The convertase enzyme PC1/3 cleaves the POMC prohormone to obtain ACTH, while the convertase enzyme PC2 cleaves ACTH to produce α-melanin-stimulating hormone (α-MSH). ACTH regulates glucocorticoids, hormones that modulate stress responses, while α-MSH induces eumelanogenesis. We therefore conclude that the melanocortin system, through the convertase enzymes PC1/3 and PC2, may account for part of the interindividual variation in melanin-based coloration in nestling tawny owls. Pleiotropy may thus account for the covariation between phenotypic traits involved in social interactions (here pigmentation) and life history, morphology, behaviour and physiology.

Knowledge of the hormonal pathway controlling genotype-specific norms of reaction would shed light on the ecological factors to which each genotype is adapted. Environmentally mediated changes in the sign and magnitude of covariations between heritable melanin-based colouration and fitness components are frequent, revealing that extreme melanin-based phenotypes can display different physiological states depending on the environment. Yet, the hormonal mechanism underlying this phenomenon is poorly understood. One novel hypothesis proposes that these covariations stem from pleiotropic effects of the melanocortin system. Melanocortins are post-translationally modified bioactive peptides derived from the POMC prohormone that are involved in melanogenesis, anti-inflammation, energy homeostasis and stress responses. Thus, differential regulation of fitness components in relation to environmental factors by pale and dark melanic individuals may be due to colour-specific regulation of the POMC prohormone. Accordingly, we found that the degree of reddish melanic colouration was negatively correlated with blood circulating levels of the POMC prohormone in female tawny owls (Strix aluco) rearing a brood for which the size was experimentally reduced, but not when enlarged, and in females located in rich but not in poor territories. Our findings support the hypothesis that the widespread links between melanin-based colouration and fitness components may be mediated, at least in part, by the melanocortin system.

Summary 1 Melanin pigments provide the most widespread source of coloration in vertebrates, but the adaptive function of such traits remains poorly known. 2 In a wild population of tawny owls (Strix aluco), we investigated the relationships between plumage coloration, which varies continuously from dark to pale reddish, and the strength and cost of an induced immune response. 3 The degree of reddishness in tawny owl feather colour was positively correlated with the concentration of phaeomelanin and eumelanin pigments, and plumage coloration was highly heritable (h2 = 0·93). No carotenoids were detected in the feathers. 4 In mothers, the degree of melanin-based coloration was associated with antibody production against a vaccine, with dark reddish females maintaining a stronger level of antibody for a longer period of time compared to pale reddish females, but at a cost in terms of greater loss of body mass. 5 A cross-fostering experiment showed that, independent of maternal coloration, foster chicks reared by vaccinated mothers were lighter than those reared by nonvaccinated mothers. Hence, even though dark reddish mothers suffered a stronger immune cost than pale reddish mothers, this asymmetric cost was not translated to offspring growth. 6 Our study suggests that different heritable melanin-based colorations are associated with alternative strategies to resist parasite attacks, with dark reddish individuals investing more resources towards the humoral immune response than lightly reddish conspecifics.

In vertebrates, melanin-based coloration is often associated with variation in physiological and behavioural traits. We propose that this association stems from pleiotropic effects of the genes regulating the synthesis of brown to black eumelanin. The most important regulators are the melanocortin 1 receptor and its ligands, the melanocortin agonists and the agouti-signalling protein antagonist. On the basis of the physiological and behavioural functions of the melanocortins, we predict five categories of traits correlated with melanin-based coloration. A review of the literature indeed reveals that, as predicted, darker wild vertebrates are more aggressive, sexually active and resistant to stress than lighter individuals. Pleiotropic effects of the melanocortins might thus account for the widespread covariance between melanin-based coloration and other phenotypic traits in vertebrates.