Yann X. C. Bourgeois

Extremely reduced dispersal and gene flow in an island bird

The Réunion grey white-eye, Zosterops borbonicus, a passerine bird endemic to Réunion Island in the Mascarene archipelago, represents an extreme case of microgeographical plumage colour variation in birds, with four distinct colour forms occupying different parts of this small island (2512u2009km2). To understand whether such population differentiation may reflect low levels of dispersal and gene flow at a very small spatial scale, we examined population structure and gene flow by analysing variation at 11 microsatellite loci among four geographically close localities (<26u2009km apart) sampled within the distribution range of one of the colour forms, the brown-headed brown form. Our results revealed levels of genetic differentiation that are exceptionally high for birds at such a small spatial scale. This strong population structure appears to reflect low levels of historical and contemporary gene flow among populations, unless very close geographically (<10u2009km). Thus, we suggest that the Réunion grey white-eye shows an extremely reduced propensity to disperse, which is likely to be related to behavioural processes.

HYBRIDIZATION AND BARRIERS TO GENE FLOW IN AN ISLAND BIRD RADIATION

While reinforcement may play a role in all major modes of speciation, relatively little is known about the timescale over which species hybridize without evolving complete reproductive isolation. Birds have high potential for hybridization, and islands provide simple settings for uncovering speciation and hybridization patterns. Here we develop a phylogenetic hypothesis for a phenotypically diverse radiation of finch‐like weaver‐birds (Foudia) endemic to the western Indian Ocean islands. We find that unlike Darwins finches, each island‐endemic Foudia population is a monophyletic entity for which speciation can be considered complete. In explaining the only exceptions—mismatches between taxonomy, mitochondrial, and nuclear data—phylogenetic and coalescent methods support introgressive hybridization rather than incomplete lineage sorting. Human introductions of known timing of one island‐endemic species, to all surrounding archipelagos provide two fortuitous experiments; (1) population sampling at known times in recent evolutionary history, (2) bringing allopatric lineages of an island radiation into secondary contact. Our results put a minimum time bound on introgression (235 years), and support hybridization between species in natural close contact (parapatry), but not between those in natural allopatry brought into contact by human introduction. Time in allopatry, rather than in sympatry, appears key in the reproductive isolation of Foudia species.

Mass production of SNP markers in a nonmodel passerine bird through RAD sequencing and contig mapping to the zebra finch genome

Here, we present an adaptation of restriction‐site‐associated DNA sequencing (RAD‐seq) to the Illumina HiSeq2000 technology that we used to produce SNP markers in very large quantities at low cost per unit in the Réunion grey white‐eye (Zosterops borbonicus), a nonmodel passerine bird species with no reference genome. We sequenced a set of six pools of 18–25 individuals using a single sequencing lane. This allowed us to build around 600 000 contigs, among which at least 386 000 could be mapped to the zebra finch (Taeniopygia guttata) genome. This yielded more than 80 000 SNPs that could be mapped unambiguously and are evenly distributed across the genome. Thus, our approach provides a good illustration of the high potential of paired‐end RAD sequencing of pooled DNA samples combined with comparative assembly to the zebra finch genome to build large contigs and characterize vast numbers of informative SNPs in nonmodel passerine bird species in a very efficient and cost‐effective way.

The role of selection and historical factors in driving population differentiation along an elevational gradient in an island bird

Adaptation to local environmental conditions and the range dynamics of populations can influence evolutionary divergence along environmental gradients. Thus, it is important to investigate patterns of both phenotypic and genetic variations among populations to reveal the respective roles of these two types of factors in driving population differentiation. Here, we test for evidence of phenotypic and genetic structure across populations of a passerine bird (Zosterops borbonicus) distributed along a steep elevational gradient on the island of Réunion. Using 11 microsatellite loci screened in 401 individuals from 18 localities distributed along the gradient, we found that genetic differentiation occurred at two spatial levels: (i) between two main population groups corresponding to highland and lowland areas, respectively, and (ii) within each of these two groups. In contrast, several morphological traits varied gradually along the gradient. Comparison of neutral genetic differentiation (FST) and phenotypic differentiation (PST) showed that PST largely exceeds FST at several morphological traits, which is consistent with a role for local adaptation in driving morphological divergence along the gradient. Overall, our results revealed an area of secondary contact midway up the gradient between two major, cryptic, population groups likely diverged in allopatry. Remarkably, local adaptation has shaped phenotypic differentiation irrespective of population history, resulting in different patterns of variation along the elevational gradient. Our findings underscore the importance of understanding both historical and selective factors when trying to explain variation along environmental gradients.

Candidate Gene Analysis Suggests Untapped Genetic Complexity in Melanin-Based Pigmentation in Birds

Studies on melanin-based color variation in a context of natural selection have provided a wealth of information on the link between phenotypic and genetic variation. Here, we evaluated associations between melanic plumage patterns and genetic polymorphism in the Réunion grey white-eye (Zosterops borbonicus), a species in which mutations on MC1R do not seem to play any role in explaining melanic variation. This species exhibits 5 plumage color variants that can be grouped into 3 color forms which occupy discrete geographic regions in the lowlands of Réunion, and a fourth high-elevation form which comprises 2 color morphs (grey and brown) and represents a true color polymorphism. We conducted a comprehensive survey of sequence variation in 96 individuals at a series of 7 candidate genes other than MC1R that have been previously shown to influence melanin-based color patterns in vertebrates, including genes that have rarely been studied in a wild bird species before: POMC, Agouti, TYR, TYRP1, DCT, Corin, and SLC24A5 Of these 7 genes, 2 (Corin and TYRP1) displayed an interesting shift in allele frequencies between lowland and highland forms and a departure from mutation-drift equilibrium consistent with balancing selection in the polymorphic highland form only. Sequence variation at Agouti, a gene frequently involved in melanin-based pigmentation patterning, was not associated with color forms or morphs. Thus, we suggest that functionally important changes in loci other than those classically studied are involved in the color polymorphism exhibited by the Réunion grey white-eye and possibly many other nonmodel species.

The genetic basis of resistance and matching-allele interactions of a host-parasite system

Negative frequency-dependent selection (NFDS) is an evolutionary mechanism suggested to govern host-parasite coevolution and the maintenance of genetic diversity at host resistance loci, such as the vertebrate MHC and R-genes in plants. Matching-allele interactions of hosts and parasites that prevent the emergence of host and parasite genotypes that are universally resistant and infective are a genetic mechanism predicted to underpin NFDS. The underlying genetics of matching-allele interactions are unknown even in host-parasite systems with empirical support for coevolution by NFDS, as is the case for the planktonic crustacean Daphnia magna and the bacterial pathogen Pasteuria ramosa. We fine-map one locus associated with D. magna resistance to P. ramosa and genetically characterize two haplotypes of the Pasteuria resistance (PR-) locus using de novo genome and transcriptome sequencing. Sequence comparison of PR-locus haplotypes finds dramatic structural polymorphisms between PR-locus haplotypes including a large portion of each haplotype being composed of non-homologous sequences resulting in haplotypes differing in size by 66 kb. The high divergence of PR-locus haplotypes suggest a history of multiple, diverse and repeated instances of structural mutation events and restricted recombination. Annotation of the haplotypes reveals striking differences in gene content. In particular, a group of glycosyltransferase genes that is present in the susceptible but absent in the resistant haplotype. Moreover, in natural populations, we find that the PR-locus polymorphism is associated with variation in resistance to different P. ramosa genotypes, pointing to the PR-locus polymorphism as being responsible for the matching-allele interactions that have been previously described for this system. Our results conclusively identify a genetic basis for the matching-allele interaction observed in a coevolving host-parasite system and provide a first insight into its molecular basis.

A novel locus on chromosome 1 underlies the evolution of a melanic plumage polymorphism in a wild songbird

Understanding the mechanisms responsible for phenotypic diversification within and among species ultimately rests with linking naturally occurring mutations to functionally and ecologically significant traits. Colour polymorphisms are of great interest in this context because discrete colour patterns within a population are often controlled by just a few genes in a common environment. We investigated how and why phenotypic diversity arose and persists in the Zosterops borbonicus white-eye of Reunion (Mascarene archipelago), a colour polymorphic songbird in which all highland populations contain individuals belonging to either a brown or a grey plumage morph. Using extensive phenotypic and genomic data, we demonstrate that this melanin-based colour polymorphism is controlled by a single locus on chromosome 1 with two large-effect alleles, which was not previously described as affecting hair or feather colour. Differences between colour morphs appear to rely upon complex cis-regulatory variation that either prevents the synthesis of pheomelanin in grey feathers, or increases its production in brown ones. We used coalescent analyses to show that, from a ‘brown’ ancestral population, the dominant ‘grey’ allele spread quickly once it arose from a new mutation. Since colour morphs are always found in mixture, this implies that the selected allele does not go to fixation, but instead reaches an intermediate frequency, as would be expected under balancing selection.

Investigating the Role of the Melanocortin-1 Receptor Gene in an Extreme Case of Microgeographical Variation in the Pattern of Melanin-Based Plumage Pigmentation

The Réunion grey white-eye (Zosterops borbonicus) is a single-island endemic passerine bird that exhibits striking geographically structured melanic polymorphism at a very small spatial scale. We investigated the genetic basis of this color polymorphism by testing whether the melanocortin-1 receptor (MC1R), a gene often involved in natural melanic polymorphism in birds, was associated with the observed plumage variation. Although we found three non-synonymous mutations, we detected no association between MC1R variants and color morphs, and the main amino-acid variant found in the Réunion grey white-eye was also present at high frequency in the Mauritius grey white-eye (Zosterops mauritianus), its sister species which shows no melanic polymorphism. In addition, neutrality tests and analysis of population structure did not reveal any obvious pattern of positive or balancing selection acting on MC1R. Altogether these results indicate that MC1R does not play a role in explaining the melanic variation observed in the Réunion grey white-eye. We propose that other genes such as POMC, Agouti or any other genes involved in pigment synthesis will need to be investigated in future studies if we are to understand how selection shapes complex patterns of melanin-based plumage pigmentation. Trial Registration All sequences submitted to Genbank. Accession number: JX914505 to JX914564.

Isolation and characterization of twelve polymorphic microsatellite loci for investigating an extreme case of microgeographical variation in an island bird (Zosterops borbonicus)

Twelve polymorphic microsatellite loci were characterized in order to investigate an extreme case of microgeographical variation in an island-dwelling passerine (Zosterops borbonicus). Variation was examined in two distinct natural populations of 27 individuals each. The loci displayed 5–20 alleles, with observed heterozygosities ranging between 0.33 and 0.93. All but one loci (Z16) conformed to Hardy–Weinberg equilibrium. None of the pairwise comparisons among loci showed significant linkage disequilibrium after sequential Bonferonni correction. Successful cross-species amplification in additional Zosterops species and other passerines suggests their broad applicability.

Parasitism drives host genome evolution: Insights from the Pasteuria ramosa–Daphnia magna system

Because parasitism is thought to play a major role in shaping host genomes, it has been predicted that genomic regions associated with resistance to parasites should stand out in genome scans, revealing signals of selection above the genomic background. To test whether parasitism is indeed such a major factor in host evolution and to better understand host–parasite interaction at the molecular level, we studied genome‐wide polymorphisms in 97 genotypes of the planktonic crustacean Daphnia magna originating from three localities across Europe. Daphnia magna is known to coevolve with the bacterial pathogen Pasteuria ramosa for which host genotypes (clonal lines) are either resistant or susceptible. Using association mapping, we identified two genomic regions involved in resistance to P. ramosa, one of which was already known from a previous QTL analysis. We then performed a naïve genome scan to test for signatures of positive selection and found that the two regions identified with the association mapping further stood out as outliers. Several other regions with evidence for selection were also found, but no link between these regions and phenotypic variation could be established. Our results are consistent with the hypothesis that parasitism is driving host genome evolution.