SEX CHROMOSOME LINKAGE OF MATE PREFERENCE AND COLOR SIGNAL MAINTAINS ASSORTATIVE MATING BETWEEN INTERBREEDING FINCH MORPHS

Assortative mating is a key aspect in the speciation process because it is important for both initial divergence and maintenance of distinct species. However, it remains a challenge to explain how assortative mating evolves when diverging populations are undergoing gene flow (e.g., during hybridization). Here I experimentally test how assortative mating is maintained with frequent gene flow between diverged head‐color morphs of the Gouldian finch (Erythrura gouldiae). Contrary to the predominant view on the development of sexual preferences in birds, cross‐fostered offspring did not imprint on the phenotype of their conspecific (red or black morphs) or heterospecific (Bengalese finch) foster parents. Instead, the mating preferences of F₁ and F₂ intermorph‐hybrids are consistent with inheritance on the Z chromosomes, which are also the location for genes controlling color expression and the genes causing low fitness of intermorph‐hybrids. Genetic associations between color signal and preference loci on the sex chromosomes may prevent recombination from breaking down these associations when the morphs interbreed, helping to maintain assortative mating in the face of gene flow. Although sex linkage of reproductively isolating traits is theoretically expected to promote speciation, social and ecological constraints may enforce frequent interbreeding between the morphs, thus preventing complete reproductive isolation.     

Sex chromosome inversions enforce reproductive isolation across an avian hybrid zone

Across hybrid zones, the sex chromosomes are often more strongly differentiated than the autosomes. This is regularly attributed to the greater frequency of reproductive incompatibilities accumulating on sex chromosomes and their exposure in the heterogametic sex. Working within an avian hybrid zone, we explore the possibility that chromosome inversions differentially accumulate on the Z chromosome compared to the autosomes and thereby contribute to Z chromosome differentiation. We analyse the northern Australian hybrid zone between two subspecies of the long‐tailed finch (Poephila acuticauda), first described based on differences in bill colour, using reduced‐representation genomic sequencing for 293 individuals over a 1,530‐km transect. Autosomal differentiation between subspecies is minimal. In contrast, 75% of the Z chromosome is highly differentiated and shows a steep genomic cline, which is displaced 350 km to the west of the cline in bill colour. Differentiation is associated with two or more putative chromosomal inversions, each predominating in one subspecies. If inversions reduce recombination between hybrid incompatibilities, they are selectively favoured and should therefore accumulate in hybrid zones. We argue that this predisposes inversions to differentially accumulate on the Z chromosome. One genomic region affecting bill colour is on the Z, but the main candidates are on chromosome 8. This and the displacement of the bill colour and Z chromosome cline centres suggest that bill colour has not strongly contributed to inversion accumulation. Based on cline width, however, the Z chromosome and bill colour both contribute to reproductive isolation established between this pair of subspecies.     

The colour of paternity: extra‐pair paternity in the wild Gouldian finch does not appear to be driven by genetic incompatibility between morphs

In socially monogamous species, individuals can use extra‐pair paternity and offspring sex allocation as adaptive strategies to ameliorate costs of genetic incompatibility with their partner. Previous studies on domesticated Gouldian finches (Erythrura gouldiae) demonstrated a genetic incompatibility between head colour morphs, the effects of which are more severe in female offspring. Domesticated females use differential sex allocation, and extra‐pair paternity with males of compatible head colour, to reduce fitness costs associated with incompatibility in mixed‐morph pairings. However, laboratory studies are an oversimplification of the complex ecological factors experienced in the wild and may only reflect the biology of a domesticated species. This study aimed to examine the patterns of parentage and sex ratio bias with respect to colour pairing combinations in a wild population of the Gouldian finch. We utilized a novel PCR assay that allowed us to genotype the morph of offspring before the morph phenotype develops and to explore bias in morph paternity and selection at the nest. Contrary to previous findings in the laboratory, we found no effect of pairing combinations on patterns of extra‐pair paternity, offspring sex ratio or selection on morphs in nestlings. In the wild, the effect of morph incompatibility is likely much smaller, or absent, than was observed in the domesticated birds. Furthermore, the previously studied domesticated population is genetically differentiated from the wild population, consistent with the effects of domestication. It is possible that the domestication process fostered the emergence (or enhancement) of incompatibility between colour morphs previously demonstrated in the laboratory.     

POSTZYGOTIC GENETIC INCOMPATIBILITY BETWEEN SYMPATRIC COLOR MORPHS

Alternative genetically determined color morphs within a population or species are believed to successfully interbreed within a population. However, the occurrence of prezygotic or ecological selection in a number of polymorphic systems may lead to nonrandom mating and prevent genetic morphs from fully interbreeding. Here we show that postzygotic incompatibility significantly limits gene flow between the sympatric red and black color morphs of the Gouldian finch ( Erythrura gouldiae ). Using a balanced within-female experimental design, in which individuals were forced to breed in pure and mixed morph crosses, we found large inviability effects (>30%) in offspring resulting from genetically mixed genotypes. The consistent mortality effects across different stages of development (e.g., prehatching, juvenile, adulthood), unconfounded by environmentally derived parental effects or social environments, reveal an underlying genetic incompatibility between different genotypes. Furthermore, mortality in mixed morph genotypes was particularly severe (43.6%) for the heterogametic sex (daughters), which is consistent with Haldane's rule predicted for postzygotic incompatibilities between hybridizing species. This significant, but incomplete, postzygotic isolation suggests that the sympatric morphs may represent transient stages in the speciation–hybridization process.     

The mutation for medullary thyroid carcinoma with parathyroid tumors (MTC with PTs) is closely linked to the centromeric region of chromosome 10.

Two new morphs (F and G) detected by the centromeric alpha satellite probe p alpha 10RP8 and D10Z1 in HinfI digests are linked to the PstI polymorphisms of D10Z1, confirming their chromosome 10 location. The F and G morphs were in strong linkage disequilibrium with each other but were in weak linkage disequilibrium with the A and B morphs defined in PstI digests. Data for haplotypes formed by using the A and F morphs improved the lod score for linkage between the disease locus for multiple endocrine neoplasia type 2A (MEN2A) and D10Z1 (Z = 14.06 at theta = 0) in the six large families studied by Wu et al. Furthermore, the locus that codes for a distinct phenotype, medullary thyroid carcinoma (MTC) with parathyroid tumors (PTs) and no pheochromocytomas (PHEOs) (referred to as MTC with PTs), in one of the families was closely linked to two markers, D10Z1 and RBP3, with lodscores of 2.86 and 3.54, respectively, at theta = 0. A possible allelic association was noted between disease phenotypes and centromeric haplotypes. The phenotype MTC and PHEOs with and without PTs was associated with the same relatively common centromeric haplotype (A + B-F-G-) in the four families in which all four morphs could be determined, while the phenotype MTC with PTs was associated with the rare centromeric haplotype (A-B-F-G+) in one family.     

Morph-dependent resource acquisition and fitness in a polymorphic bird

Understanding genetic colour polymorphism has proved a major challenge, both in terms of the underlying genetic mechanisms and the evolutionarily forces maintaining such genetic variation. In this context, genetic differences in aggression or competitive-related traits may covary with the expression of alternative phenotypes, and affect the evolutionary stability and maintenance of colour polymorphisms. Genetic red and black head-colour morphs of the Gouldian finch (Erythrura gouldiae) co-occur in temporally and geographically stable frequencies in sympatric populations. Gouldian finches are obligate cavity-nesters with highly specific preferences for nest-site morphometry that directly affect reproductive success. Because intra- and interspecific competition for high quality nest-sites is prevalent, and fitness is directly related to nest-site quality, we investigated the relative access (and consequences for reproductive success) of alternative morphs to this critical limiting resource in the wild. Red males defended higher quality nest-sites, and overcame greater levels of nest-site competition against conspecifics and superior heterospecific competitors than black males. Red-headed males also produced more fledglings (especially with red-headed females) than black-headed males, independent of nest-site quality. Finally, the independent (positive) effect of nest-site quality on reproductive success was confirmed. Such competitive asymmetries are important to relative selection among coexisting morphs, and are likely to contribute to the maintenance of alternative sympatric colour-morphs in wild populations.     

Contrasting patterns of introgression at X-linked loci across the hybrid zone between subspecies of the European rabbit (Oryctolagus cuniculus)

Hybrid zones provide an excellent opportunity for studying the consequences of genetic changes between closely related taxa. Here we investigate patterns of genetic variability and gene flow at four X-linked loci within and between the two subspecies of European rabbit (Oryctolagus cuniculus cuniculus and O. c. algirus). Two of these genes are located near the centromere and two are located near the telomeres. We observed a deep split in the genealogy of each gene with the root located along the deepest branch in each case, consistent with the evolution of these subspecies in allopatry. The two centromeric loci showed low levels of variability, high levels of linkage disequilibrium, and little introgression between subspecies. In contrast, the two telomeric loci showed high levels of variability, low levels of linkage disequilibrium, and considerable introgression between subspecies. These data are consistent with suppression of recombination near the centromere of the rabbit X chromosome. These observations support a view of speciation where genomic incompatibilities at different loci in the genome create localized differences in levels of gene flow between nascent species.     

A high‐density linkage map enables a second‐generation collared flycatcher genome assembly and reveals the patterns of avian recombination rate variation and chromosomal evolution

Detailed linkage and recombination rate maps are necessary to use the full potential of genome sequencing and population genomic analyses. We used a custom collared flycatcher 50 K SNP array to develop a high‐density linkage map with 37 262 markers assigned to 34 linkage groups in 33 autosomes and the Z chromosome. The best‐order map contained 4215 markers, with a total distance of 3132 cM and a mean genetic distance between markers of 0.12 cM . Facilitated by the array being designed to include markers from most scaffolds, we obtained a second‐generation assembly of the flycatcher genome that approaches full chromosome sequences (N50 super‐scaffold size 20.2 Mb and with 1.042 Gb (of 1.116 Gb) anchored to and mostly ordered and oriented along chromosomes). We found that flycatcher and zebra finch chromosomes are entirely syntenic but that inversions at mean rates of 1.5–2.0 event (6.6–7.5 Mb) per My have changed the organization within chromosomes, rates high enough for inversions to potentially have been involved with many speciation events during avian evolution. The mean recombination rate was 3.1 cM /Mb and correlated closely with chromosome size, from 2 cM /Mb for chromosomes >100 Mb to >10 cM /Mb for chromosomes <10 Mb. This size dependence seemed entirely due to an obligate recombination event per chromosome; if 50 cM was subtracted from the genetic lengths of chromosomes, the rate per physical unit DNA was constant across chromosomes. Flycatcher recombination rate showed similar variation along chromosomes as chicken but lacked the large interior recombination deserts characteristic of zebra finch chromosomes.     

Nucleotide Variation, Linkage Disequilibrium and Founder-Facilitated Speciation in Wild Populations of the Zebra Finch (Taeniopygia guttata)

The zebra finch has long been an important model system for the study of vocal learning, vocal production, and behavior. With the imminent sequencing of its genome, the zebra finch is now poised to become a model system for population genetics. Using a panel of 30 noncoding loci, we characterized patterns of polymorphism and divergence among wild zebra finch populations. Continental Australian populations displayed little population structure, exceptionally high levels of nucleotide diversity (π = 0.010), a rapid decay of linkage disequilibrium (LD), and a high population recombination rate (ρ ≈ 0.05), all of which suggest an open and fluid genomic background that could facilitate adaptive variation. By contrast, substantial divergence between the Australian and Lesser Sunda Island populations (K_ST = 0.193), reduced genetic diversity (π = 0.002), and higher levels of LD in the island population suggest a strong but relatively recent founder event, which may have contributed to speciation between these populations as envisioned under founder-effect speciation models. Consistent with this hypothesis, we find that under a simple quantitative genetic model both drift and selection could have contributed to the observed divergence in six quantitative traits. In both Australian and Lesser Sundas populations, diversity in Z-linked loci was significantly lower than in autosomal loci. Our analysis provides a quantitative framework for studying the role of selection and drift in shaping patterns of molecular evolution in the zebra finch genome.     

No evidence for Z‐chromosome rearrangements between the pied flycatcher and the collared flycatcher as judged by gene‐based comparative genetic maps

Revealing the genetic basis of reproductive isolation is fundamental for understanding the speciation process. Chromosome speciation models propose a role for chromosomal rearrangements in promoting the build up of reproductive isolation between diverging populations and empirical data from several animal and plant taxa support these models. The pied flycatcher and the collared flycatcher are two closely related species that probably evolved reproductive isolation during geographical separation in Pleistocene glaciation refugia. Despite the short divergence time and current hybridization, these two species demonstrate a high degree of intrinsic post‐zygotic isolation and previous studies have shown that traits involved in mate choice and hybrid viability map to the Z‐chromosome. Could rearrangements of the Z‐chromosome between the species explain their reproductive isolation? We developed high coverage Z‐chromosome linkage maps for both species, using gene‐based markers and large‐scale SNP genotyping. Best order maps contained 57–62 gene markers with an estimated average density of one every 1–1.5 Mb. We estimated the recombination rates in flycatcher Z‐chromosomes to 1.1–1.3 cM/Mb. A comparison of the maps of the two species revealed extensive co‐linearity with no strong evidence for chromosomal rearrangements. This study does therefore not provide support the idea that sex chromosome rearrangements have caused the relatively strong post‐zygotic reproductive isolation between these two Ficedula species.     

Interspecific crossing and genetic mapping reveal intrinsic genomic incompatibility between two Senecio species that form a hybrid zone on Mount Etna,Sicily

Studies of hybridizing species can reveal much about the genetic basis and maintenance of species divergence in the face of gene flow. Here we report a genetic segregation and linkage analysis conducted on F₂ progeny of a reciprocal cross between Senecio aethnensis and S. chrysanthemifolius that form a hybrid zone on Mount Etna, Sicily, aimed at determining the genetic basis of intrinsic hybrid barriers between them. Significant transmission ratio distortion (TRD) was detected at 34 (∼27%) of 127 marker loci located in nine distinct clusters across seven of the ten linkage groups detected, indicating genomic incompatibility between the species. TRD at these loci could not be attributed entirely to post-zygotic selective loss of F₂ individuals that failed to germinate or flower (16.7%). At four loci tests indicated that pre-zygotic events, such as meiotic drive in F₁ parents or gametophytic selection, contributed to TRD. Additional tests revealed that cytonuclear incompatibility contributed to TRD at five loci, Bateson–Dobzhansky–Muller (BDM) incompatibilities involving epistatic interactions between loci contributed to TRD at four loci, and underdominance (heterozygote disadvantage) was a possible cause of TRD at one locus. Major chromosomal rearrangements were probably not a cause of interspecific incompatibility at the scale that could be examined with current map marker density. Intrinsic genomic incompatibility between S. aethnensis and S. chrysanthemifolius revealed by TRD across multiple genomic regions in early-generation hybrids is likely to impact the genetic structure of the natural hybrid zone on Mount Etna by limiting introgression and promoting divergence across the genome.     

Red dominates black: agonistic signalling among head morphs in the colour polymorphic Gouldian finch

Recent sexual selection studies on the evolution of bird colouration have mainly focused on signals with a high level of condition-dependent variation, with much less attention given to colour traits whose expression is genetically controlled. Here, we experimentally tested the relative importance of a genetic colour polymorphism in determining male dominance in the Gouldian finch (Erythrura gouldiae), a species displaying three completely discrete but naturally co-occurring genetically inherited phenotypes; yellow-, red- (carotenoid) and black-headed (melanin) morphs. First, in staged dominance contests between unfamiliar birds of different head morphs, red-headed males dominated black-headed males, both of which dominated the yellow-headed birds. Second, within morphs, the intensity and size of the strongly ultraviolet-blue collar determined the outcome of these contests, and among the red-headed males, redder males dominated less chromatic birds. Lastly, when the dominance signal of red-headed birds was experimentally destabilized (i.e. blackened or reddened), naturally red-headed morphs continued to dominate both the black-and yellow-headed morphs. Together, these results suggest that intrinsic dominance-related behavioural differences between the three colour morphs, which are likely to influence the relative fitness of each morph, contribute to the complex selective patterns maintaining these three discrete phenotypes in relatively stable frequencies in wild populations.     

Frequency-dependent physiological trade-offs between competing colour morphs

Evolutionary theory suggests that alternative colour morphs (i.e. genetically controlled phenotypes) may derive similar fitness under frequency-dependent selection. Here we experimentally demonstrate opposing effects of frequency-dependent social environments on plasma hormone levels (testosterone and corticosterone) and immune function between red- and black-headed male morphs of the Gouldian finch (Erythrura gouldiae). Red-headed males are highly sensitive to changes in the social environment, especially towards the relative density of their own aggressive morph, exhibiting high stress responses and immunosuppression in socially competitive environments. In contrast, the non-aggressive black-headed males follow a more passive strategy that appears to buffer them against social stresses. The differential effect of hormones on aggressive behaviour and immune performance reinforces the contrasting behavioural strategies employed by these colour morphs, and highlights the importance of the social environment in determining the individual basis of behavioural and physiological responses.     

Compelling evidence that a single nucleotide substitution in TYRP1 is responsible for coat-colour polymorphism in a free-living population of Soay sheep

Identifying the genes that underlie phenotypic variation in natural populations is a central objective of evolutionary genetics. Here, we report the identification of the gene and causal mutation underlying coat colour variation in a free-living population of Soay sheep (Ovis aries). We targeted tyrosinase-related protein 1 (TYRP1), a positional candidate gene based on previous work that mapped the Coat colour locus to an approximately 15cM window on chromosome 2. We identified a non-synonymous substitution in exon IV that was perfectly associated with coat colour. This polymorphism is predicted to cause the loss of a cysteine residue that is highly evolutionarily conserved and likely to be of functional significance. We eliminated the possibility that this association is due to the presence of strong linkage disequilibrium with an unknown regulatory mutation by demonstrating that there is no difference in relative TYRP1 expression between colour morphs. Analysis of this putative causal mutation in a complex pedigree of more than 500 sheep revealed almost perfect co-segregation with coat colour (chi2-test, p<0.0001, LOD=110.20), and very tight linkage between Coat colour and TYRP1 (LOD=29.50).     

A microsatellite‐based linkage map for song sparrows (Melospiza melodia)

Although linkage maps are important tools in evolutionary biology, their availability for wild populations is limited. The population of song sparrows (Melospiza melodia) on Mandarte Island, Canada, is among the more intensively studied wild animal populations. Its long‐term pedigree data, together with extensive genetic sampling, have allowed the study of a range of questions in evolutionary biology and ecology. However, the availability of genetic markers has been limited. We here describe 191 new microsatellite loci, including 160 high‐quality polymorphic autosomal, 7 Z‐linked and 1 W‐linked markers. We used these markers to construct a linkage map for song sparrows with a total sex‐averaged map length of 1731 cM and covering 35 linkage groups, and hence, these markers cover most of the 38–40 chromosomes. Female and male map lengths did not differ significantly. We then bioinformatically mapped these loci to the zebra finch (Taeniopygia guttata) genome and found that linkage groups were conserved between song sparrows and zebra finches. Compared to the zebra finch, marker order within small linkage groups was well conserved, whereas the larger linkage groups showed some intrachromosomal rearrangements. Finally, we show that as expected, recombination frequency between linked loci explained the majority of variation in gametic phase disequilibrium. Yet, there was substantial overlap in gametic phase disequilibrium between pairs of linked and unlinked loci. Given that the microsatellites described here lie on 35 of the 38–40 chromosomes, these markers will be useful for studies in this species, as well as for comparative genomics studies with other species.     

Recombinational Landscape and Population Genomics of Caenorhabditis elegans

Recombination rate and linkage disequilibrium, the latter a function of population genomic processes, are the critical parameters for mapping by linkage and association, and their patterns in Caenorhabditis elegans are poorly understood. We performed high-density SNP genotyping on a large panel of recombinant inbred advanced intercross lines (RIAILs) of C. elegans to characterize the landscape of recombination and, on a panel of wild strains, to characterize population genomic patterns. We confirmed that C. elegans autosomes exhibit discrete domains of nearly constant recombination rate, and we show, for the first time, that the pattern holds for the X chromosome as well. The terminal domains of each chromosome, spanning about 7% of the genome, exhibit effectively no recombination. The RIAILs exhibit a 5.3-fold expansion of the genetic map. With median marker spacing of 61 kb, they are a powerful resource for mapping quantitative trait loci in C. elegans. Among 125 wild isolates, we identified only 41 distinct haplotypes. The patterns of genotypic similarity suggest that some presumed wild strains are laboratory contaminants. The Hawaiian strain, CB4856, exhibits genetic isolation from the remainder of the global population, whose members exhibit ample evidence of intercrossing and recombining. The population effective recombination rate, estimated from the pattern of linkage disequilibrium, is correlated with the estimated meiotic recombination rate, but its magnitude implies that the effective rate of outcrossing is extremely low, corroborating reports of selection against recombinant genotypes. Despite the low population, effective recombination rate and extensive linkage disequilibrium among chromosomes, which are techniques that account for background levels of genomic similarity, permit association mapping in wild C. elegans strains.     

Differential barrier strength and allele frequencies in hybrid zones maintained by sex-biased hybrid incompatibilities

In animals, hybrid sterility and inviability between closely related species often affect only the heterogametic sex (XY). This widespread phenomenon, known as Haldane's rule, is an early speciation event found across broad taxa, but the role of heterogametic hybrid incompatibilities, as opposed to homogametic ones, as a barrier in a speciation process remains obscure. It has been hypothesized that heterogametic incompatibility may be a more efficient mechanism in driving speciation. The population dynamics after (rather than before) the occurrence of sex-biased incompatibilities may account for Haldane's rule. In this study, a recursion model of hybrid zones was developed to investigate the differences between heterogametic and homogametic incompatibilities. The selection strengths and selection patterns of sex chromosome-linked, two-locus Bateson-Dobzhansky-Muller (BDM) genetic incompatibilities were examined. It is noted that a sex-biased hybrid incompatibility in a hybrid zone confers asymmetric and uneven impedance to gene flow. The clines of different loci in such a hybrid zone displayed diverse differentiation in their width, steepness and asymmetry. Alleles involved in the incompatibility face much stronger resistance to cross a hybrid zone. Different sex-biased BDM incompatibilities also affect the flow of neutral alleles differently. Compared to a homogametic one, heterogametic incompatibility is a weaker but more asymmetric barrier. These unique patterns of gene flow may explain uneven divergence among different genomic regions during speciation between some closely related species.     

A linkage map of the zebra finch Taeniopygia guttata provides new insights into avian genome evolution

Passeriformes are the largest order of birds and one of the most widely studied groups in evolutionary biology and ecology. Until recently genomic tools in passerines relied on chicken genomic resources. Here we report the construction and analysis of a whole-genome linkage map for the zebra finch (Taeniopygia guttata) using a 354-bird pedigree. The map contains 876 SNPs dispersed across 45 linkage groups and we found only a few instances of interchromosomal rearrangement between the zebra finch and the chicken genomes. Interestingly, there was a greater than expected degree of intrachromosomal rearrangements compared to the chicken, suggesting that gene order is not conserved within avian chromosomes. At 1068 cM the map is approximately only one quarter the length of the chicken linkage map, providing further evidence that the chicken has an unusually high recombination rate. Male and female linkage-map lengths were similar, suggesting no heterochiasmy in the zebra finch. This whole-genome map is the first for any passerine and a valuable tool for the zebra finch genome sequence project and for studies of quantitative trait loci.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.