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.     

Fecundity selection does not vary along a large geographical cline of trait means in a passerine bird

Local environmental and ecological conditions are commonly expected to result in local adaptation, although there are few examples of variation in phenotypic selection across continent‐wide spatial scales. We collected standardized data on selection with respect to the highly variable plumage coloration of pied flycatcher (Ficedula hypoleuca Pall.) males from 17 populations across the species' breeding range. The observed selection on multiple male coloration traits via the annual number of fledged young was generally relatively weak. The main aim of the present study, however, was to examine whether the current directional selection estimates are associated with distance to the sympatric area with the collared flycatcher (Ficedula albicollis Temminck), a sister species with which the pied flycatcher is showing character displacement. This pattern was expected because plumage traits in male pied flycatchers are changing with the distance to these areas of sympatry. However, we did not find such a pattern in current selection on coloration. There were no associations between current directional selection on ornamentation and latitude or longitude either. Interestingly, current selection on coloration traits was not associated with the observed mean plumage traits of the populations. Thus, there do not appear to be geographical gradients in current directional fecundity selection on male plumage ornamentation. The results of the present study do not support the idea that constant patterns in directional fecundity selection would play a major role in the maintenance of coloration among populations in this species. By contrast, the tendency for relatively weak mosaic‐like variation in selection among populations could reflect just a snapshot of temporally variable, potentially environment‐dependent, selection, as suggested by other studies in this system. Such fine‐grained variable selection coupled with gene flow could maintain extensive phenotypic variation across populations. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 808–827.     

Genetic mapping in a natural population of collared flycatchers (Ficedula albicollis): conserved synteny but gene order rearrangements on the avian Z chromosome

Data from completely sequenced genomes are likely to open the way for novel studies of the genetics of nonmodel organisms, in particular when it comes to the identification and analysis of genes responsible for traits that are under selection in natural populations. Here we use the draft sequence of the chicken genome as a starting point for linkage mapping in a wild bird species, the collared flycatcher - one of the most well-studied avian species in ecological and evolutionary research. A pedigree of 365 flycatchers was established and genotyped for single nucleotide polymorphisms in 23 genes selected from (and spread over most of) the chicken Z chromosome. All genes were also found to be located on the Z chromosome in the collared flycatcher, confirming conserved synteny at the level of gene content across distantly related avian lineages. This high degree of conservation mimics the situation seen for the mammalian X chromosome and may thus be a general feature in sex chromosome evolution, irrespective of whether there is male or female heterogamety. Alternatively, such unprecedented chromosomal conservation may be characteristic of most chromosomes in avian genome evolution. However, several internal rearrangements were observed, meaning that the transfer of map information from chicken to nonmodel bird species cannot always assume conserved gene orders. Interestingly, the rate of recombination on the Z chromosome of collared flycatchers was only approximately 50% that of chicken, challenging the widely held view that birds generally have high recombination rates.     

Flycatcher song in allopatry and sympatry--convergence, divergence and reinforcement

The theory of reinforcement predicts that natural selection against the production of unfit hybrids favours traits that increase assortative mating. Whether culturally inherited traits, such as bird song, can increase assortative mating by reinforcement is largely unknown. We compared songs of pied (Ficedula hypoleuca) and collared flycatchers (F. albicollis) from two hybrid zones of different ages with songs from allopatric populations. Previously, a character divergence in male plumage traits has been shown to reinforce premating isolation in sympatric flycatchers. In contrast, we find that the song of the pied flycatcher has converged towards that of the collared flycatcher (mixed singing). However, a corresponding divergence in the collared flycatcher shows that the species differences in song characters are maintained in sympatry. Genetic analyses suggest that mixed song is not caused by introgression from the collared flycatcher, but rather due to heterospecific copying. Circumstantial evidence suggests that mixed song may increase the rate of maladaptive hybridization. In the oldest hybrid zone where reinforcement on plumage traits is most pronounced, the frequency of mixed singing and hybridization is also lowest. Thus, we suggest that reinforcement has reduced the frequency of mixed singing in the pied flycatcher and caused a divergence in the song of the collared flycatcher. Whether a culturally inherited trait promotes or opposes speciation in sympatry may depend on its plasticity. The degree of plasticity may be genetically determined and accordingly under selection by reinforcement.     

Difference in plasticity of resting metabolic rate – the proximate explanation to different niche breadth in sympatric Ficedula flycatchers

Variation in relative fitness of competing recently formed species across heterogeneous environments promotes coexistence. However, the physiological traits mediating such variation in relative fitness have rarely been identified. Resting metabolic rate (RMR) is tightly associated with life history strategies, thermoregulation, diet use, and inhabited latitude and could therefore moderate differences in fitness responses to fluctuations in local environments, particularly when species have adapted to different climates in allopatry. We work in a long‐term study of collared (Ficedula albicollis) and pied flycatchers (Ficedula hypoleuca) in a recent hybrid zone located on the Swedish island of Öland in the Baltic Sea. Here, we explore whether differences in RMR match changes in relative performance of growing flycatcher nestlings across environmental conditions using an experimental approach. The fitness of pied flycatchers has previously been shown to be less sensitive to the mismatch between the peak in food abundance and nestling growth among late breeders. Here, we find that pied flycatcher nestlings have lower RMR in response to higher ambient temperatures (associated with low food availability). We also find that experimentally relaxed nestling competition is associated with an increased RMR in this species. In contrast, collared flycatcher nestlings did not vary their RMR in response to these environmental factors. Our results suggest that a more flexible nestling RMR in pied flycatchers is responsible for the better adaptation of pied flycatchers to the typical seasonal changes in food availability experienced in this hybrid zone. Generally, subtle physiological differences that have evolved when species were in allopatry may play an important role to patterns of competition, coexistence, or displacements between closely related species in secondary contact.     

Sympatric divergence and clinal variation in multiple coloration traits of Ficedula flycatchers

Geographic variation in phenotypes plays a key role in fundamental evolutionary processes such as local adaptation, population differentiation and speciation, but the selective forces behind it are rarely known. We found support for the hypothesis that geographic variation in plumage traits of the pied flycatcher Ficedula hypoleuca is explained by character displacement with the collared flycatcher Ficedula albicollis in the contact zone. The plumage traits of the pied flycatcher differed strongly from the more conspicuous collared flycatcher in a sympatric area but increased in conspicuousness with increasing distance to there. Phenotypic differentiation (P_ST) was higher than that in neutral genetic markers (F_ST), and the effect of geographic distance remained when statistically controlling for neutral genetic differentiation. This suggests that a cline created by character displacement and gene flow explains phenotypic variation across the distribution of this species. The different plumage traits of the pied flycatcher are strongly to moderately correlated, indicating that they evolve non‐independently from each other. The flycatchers provide an example of plumage patterns diverging in two species that differ in several aspects of appearance. The divergence in sympatry and convergence in allopatry in these birds provide a possibility to study the evolutionary mechanisms behind the highly divergent avian plumage patterns.     

Footprints of adaptive evolution revealed by whole Z chromosomes haplotypes in flycatchers

Detecting positive selection using genomic data is critical to understanding the role of adaptive evolution. Of particular interest in this context is sex chromosomes since they are thought to play a special role in local adaptation and speciation. We sought to circumvent the challenges associated with statistical phasing when using haplotype‐based statistics in sweep scans by benefitting from that whole chromosome haplotypes of the sex chromosomes can be obtained by resequencing of individuals of the hemizygous sex. We analyzed whole Z chromosome haplotypes from 100 females from several populations of four black and white flycatcher species (in birds, females are ZW and males ZZ). Based on integrated haplotype score (iHS) and number of segregating sites by length (nSL) statistics, we found strong and frequent haplotype structure in several regions of the Z chromosome in each species. Most of these sweep signals were population‐specific, with essentially no evidence for regions under selection shared among species. Some completed sweeps were revealed by the cross‐population extended haplotype homozygosity (XP‐EHH) statistic. Importantly, by using statistically phased Z chromosome data from resequencing of males, we failed to recover the signals of selection detected in analyses based on whole chromosome haplotypes from females; instead, what likely represent false signals of selection were frequently seen. This highlights the power issues in statistical phasing and cautions against conclusions from selection scans using such data. The detection of frequent selective sweeps on the avian Z chromosome supports a large role of sex chromosomes in adaptive evolution.     

Candidate genes for colour and vision exhibit signals of selection across the pied flycatcher (Ficedula hypoleuca) breeding range

The role of natural selection in shaping adaptive trait differentiation in natural populations has long been recognized. Determining its molecular basis, however, remains a challenge. Here, we search for signals of selection in candidate genes for colour and its perception in a passerine bird. Pied flycatcher plumage varies geographically in both its structural and pigment-based properties. Both characteristics appear to be shaped by selection. A single-locus outlier test revealed 2 of 14 loci to show significantly elevated signals of divergence. The first of these, the follistatin gene, is expressed in the developing feather bud and is found in pathways with genes that determine the structure of feathers and may thus be important in generating variation in structural colouration. The second is a gene potentially underlying the ability to detect this variation: SWS1 opsin. These two loci were most differentiated in two Spanish pied flycatcher populations, which are also among the populations that have the highest UV reflectance. The follistatin and SWS1 opsin genes thus provide strong candidates for future investigations on the molecular basis of adaptively significant traits and their co-evolution.     

Genomic analyses reveal natural selection on reproduction related genes between two closely related Populus (Salicaceae) species

Identifying the factors that cause reproductive isolation and their relative importance in species divergence is crucial to our understanding of speciation processes. In most species, natural selection is commonly considered to play a large role in driving speciation. Based on whole genome re-sequencing data from 27 Populus alba and 28 Populus adenopoda individuals, we explored the factors related to reproductive isolation of these two closely related species. The results showed that the two species diverged ~5–10 million years ago (Ma), when the Qinghai–Tibet Plateau reached a certain height and the inland climate of the Asian continent became arid. In highly differentiated genomic regions, the relative divergence (F_ST) and absolute divergence (d_xy) were significantly higher than the genomic background, θπ and shared polymorphisms decreased whereas fixed differences increased, which indicated that natural selection played a key role in the reproductive isolation of the two species. In addition, we found several genes that were related to reproduction that may be involved in explaining the reproductive isolation. Using phylogenetic trees resolved from haplotype data of Populus tomentosa and P. adenopoda, the maternal origin of P. tomentosa from P. adenopoda was likely to be located in Hubei and Chongqing Provinces.     

Tracking changes in chromosomal arrangements and their genetic content during adaptation

There is considerable evidence for an adaptive role of inversions, but how their genetic content evolves and affects the subsequent evolution of chromosomal polymorphism remains controversial. Here, we track how life‐history traits, chromosomal arrangements and 22 microsatellites, within and outside inversions, change in three replicated populations of Drosophila subobscura for 30 generations of laboratory evolution since founding from the wild. The dynamics of fitness‐related traits indicated adaptation to the new environment concomitant with directional evolution of chromosomal polymorphism. Evidence of selective changes in frequency of inversions was obtained for seven of 23 chromosomal arrangements, corroborating a role for inversions in adaptation. The evolution of linkage disequilibrium between some microsatellites and chromosomes suggested that adaptive changes in arrangements involved changes in their genetic content. Several microsatellite alleles increased in frequency more than expected by drift in targeted inversions in all replicate populations. In particular, there were signs of selection in the O₃₊₄ arrangement favouring a combination of alleles in two loci linked to the inversion and changing along with it, although the lack of linkage disequilibrium between these loci precludes epistatic selection. Seven other alleles increased in frequency within inversions more than expected by drift, but were not in linkage disequilibrium with them. Possibly these alleles were hitchhiking along with alleles under selection that were not specific to those inversions. Overall, the selection detected on the genetic content of inversions, despite limited coverage of the genome, suggests that genetic changes within inversions play an important role in adaptation.     

Parent‐of‐origin growth effects and the evolution of hybrid inviability in dwarf hamsters

Mammalian hybrids often show abnormal growth, indicating that developmental inviability may play an important role in mammalian speciation. Yet, it is unclear if this recurrent phenotype reflects a common genetic basis. Here, we describe extreme parent‐of‐origin‐dependent growth in hybrids from crosses between two species of dwarf hamsters, Phodopus campbelli and Phodopus sungorus. One cross type resulted in massive placental and embryonic overgrowth, severe developmental defects, and maternal death. Embryos from the reciprocal cross were viable and normal sized, but adult hybrid males were relatively small. These effects are strikingly similar to patterns from several other mammalian hybrids. Using comparative sequence data from dwarf hamsters and several other hybridizing mammals, we argue that extreme hybrid growth can contribute to reproductive isolation during the early stages of species divergence. Next, we tested if abnormal growth in hybrid hamsters was associated with disrupted genomic imprinting. We found no association between imprinting status at several candidate genes and hybrid growth, though two interacting genes involved in embryonic growth did show reduced expression in overgrown hybrids. Collectively, our study indicates that growth‐related hybrid inviability may play an important role in mammalian speciation but that the genetic underpinnings of these phenotypes remain unresolved.     

Contrasting patterns of polymorphism and divergence on the Z chromosome and autosomes in two Ficedula flycatcher species

In geographic areas where pied and collared flycatchers (Ficedula hypoleuca and F. albicollis) breed in sympatry, hybridization occurs, leading to gene flow (introgression) between the two recently diverged species. Notably, while such introgression is observable at autosomal loci it is apparently absent at the Z chromosome, suggesting an important role for genes on the Z chromosome in creating reproductive isolation during speciation. To further understand the role of Z-linked loci in the formation of new species, we studied genetic variation of the two species from regions where they live in allopatry. We analyzed patterns of polymorphism and divergence in introns from 9 Z-linked and 23 autosomal genes in pied and collared flycatcher males. Average variation on the Z chromosome is greatly reduced compared to neutral expectations based on autosomal diversity in both species. We also observe significant heterogeneity between patterns of polymorphism and divergence at Z-linked loci and a relative absence of polymorphisms that are shared by the two species on the Z chromosome compared to the autosomes. We suggest that these observations may indicate the action of recurrent selective sweeps on the Z chromosome during the evolution of the two species, which may be caused by sexual selection acting on Z-linked genes. Alternatively, reduced variation on the Z chromosome could result from substantially higher levels of introgression at autosomal than at Z-linked loci or from a complex demographic history, such as a population bottleneck.     

The genetic architecture of sexually selected traits in two natural populations of Drosophila montana

We investigated the genetic architecture of courtship song and cuticular hydrocarbon traits in two phygenetically distinct populations of Drosophila montana. To study natural variation in these two important traits, we analysed within-population crosses among individuals sampled from the wild. Hence, the genetic variation analysed should represent that available for natural and sexual selection to act upon. In contrast to previous between-population crosses in this species, no major quantitative trait loci (QTLs) were detected, perhaps because the between-population QTLs were due to fixed differences between the populations. Partitioning the trait variation to chromosomes suggested a broadly polygenic genetic architecture of within-population variation, although some chromosomes explained more variation in one population compared with the other. Studies of natural variation provide an important contrast to crosses between species or divergent lines, but our analysis highlights recent concerns that segregating variation within populations for important quantitative ecological traits may largely consist of small effect alleles, difficult to detect with studies of moderate power.     

Can phenotypic selection on floral traits explain the presence of enigmatic intermediate individuals in sympatric populations of Platanthera bifolia and P. chlorantha (Orchidaceae)?

Pollinators represent one of the main agents of selection on floral traits. Here, we estimated phenotypic selection on floral morphology and phenology in a sympatric population of two orchid species, Platanthera bifolia and P. chlorantha, including enigmatic individuals with intermediate column morphology (as reflected by the distance between viscidia and caudicle length, two traits involved in assortative mating and reproductive isolation among Platanthera species), but genetically indistinguishable from P. bifolia. Our aim was to clarify whether the occurrence of intermediate phenotypes could be explained by the presence of selective pressures exerted by pollinators. Simple linear and quadratic regressions together with univariate and multivariate analyses were used to evaluate the strength of directional, disruptive and stabilizing selection. We found that selection on phenotypic traits varied between groups and sex functions. Contrary to our hypothesis, selection on the viscidia distance and caudicle length appeared to be consistent in the two P. bifolia groups. Interestingly, the viscidia distance was under significant stabilizing selection through female reproductive success in intermediate individuals. Based on these results, we conclude that, despite a significant selective pressure on some phenotypic traits, the presence of individuals with intermediate phenotype is not due to selection. Stabilizing selection on distance between viscidia in intermediate individuals may suggest that assortative mating play a role in the maintenance of this phenotypic polymorphism.     

A role for migration‐linked genes and genomic islands in divergence of a songbird

Next‐generation sequencing has made it possible to begin asking questions about the process of divergence at the level of the genome. For example, recently, there has been a debate around the role of ‘genomic islands of divergence’ (i.e. blocks of outlier loci) in facilitating the process of speciation‐with‐gene‐flow. The Swainson's thrush, Catharus ustulatus, is a migratory songbird with two genetically distinct subspecies that differ in a number of traits known to be involved in reproductive isolation in birds (plumage coloration, song and migratory behaviour), despite contemporary gene flow along a secondary contact zone. Here, we use RAD‐PE sequencing to test emerging hypotheses about the process of divergence at the level of the genome and identify genes and gene regions involved in differentiation in this migratory songbird. Our analyses revealed distinct genomic islands on 15 of the 23 chromosomes and an accelerated rate of divergence on the Z chromosome, one of the avian sex chromosomes. Further, an analysis of loci linked to traits known to be involved in reproductive isolation in songbirds showed that genes linked to migration are significantly more differentiated than expected by chance, but that these genes lie primarily outside the genomic islands. Overall, our analysis supports the idea that genes linked to migration play an important role in divergence in migratory songbirds, but we find no compelling evidence that the observed genomic islands are facilitating adaptive divergence in migratory behaviour.     

Patterns of divergence across the geographic and genomic landscape of a butterfly hybrid zone associated with a climatic gradient

Hybrid zones are a valuable tool for studying the process of speciation and for identifying the genomic regions undergoing divergence and the ecological (extrinsic) and nonecological (intrinsic) factors involved. Here, we explored the genomic and geographic landscape of divergence in a hybrid zone between Papilio glaucus and Papilio canadensis. Using a genome scan of 28,417 ddRAD SNPs, we identified genomic regions under possible selection and examined their distribution in the context of previously identified candidate genes for ecological adaptations. We showed that differentiation was genomewide, including multiple candidate genes for ecological adaptations, particularly those involved in seasonal adaptation and host plant detoxification. The Z chromosome and four autosomes showed a disproportionate amount of differentiation, suggesting genes on these chromosomes play a potential role in reproductive isolation. Cline analyses of significantly differentiated genomic SNPs, and of species‐diagnostic genetic markers, showed a high degree of geographic coincidence (81%) and concordance (80%) and were associated with the geographic distribution of a climate‐mediated developmental threshold (length of the growing season). A relatively large proportion (1.3%) of the outliers for divergent selection were not associated with candidate genes for ecological adaptations and may reflect the presence of previously unrecognized intrinsic barriers between these species. These results suggest that exogenous (climate‐mediated) and endogenous (unknown) clines may have become coupled and act together to reinforce reproductive isolation. This approach of assessing divergence across both the genomic and geographic landscape can provide insight about the interplay between the genetic architecture of reproductive isolation and endogenous and exogenous selection.     

QTL analysis of white blood cell, platelet and red blood cell-related traits in an F2 intercross between Landrace and Korean native pigs

Haematological traits play important roles in disease resistance and defence functions. The objective of this study was to locate quantitative trait loci (QTL) and the associated positional candidate genes influencing haematological traits in an F₂ intercross between Landrace and Korean native pigs. Eight blood‐related traits (six erythrocyte traits, one leucocyte trait and one platelet trait) were measured in 816 F₂ progeny. All experimental animals were genotyped with 173 informative microsatellite markers located throughout the pig genome. We report that nine chromosomes harboured QTL for the baseline blood parameters: genomic regions on SSC 1, 4, 5, 6, 8, 9, 11, 13 and 17. Eight of twenty identified QTL reached genome‐wide significance. In addition, we evaluated the KIT locus, an obvious candidate gene locus affecting variation in blood‐related traits. Using dense single nucleotide polymorphism marker data on SSC 8 and the marker‐assisted association test, the strong association of the KIT locus with blood phenotypes was confirmed. In conclusion, our study identified both previously reported and novel QTL affecting baseline haematological parameters in pigs. Additionally, the positional candidate genes identified here could play an important role in elucidating the genetic architecture of haematological phenotype variation in swine and in humans.     

QTLs mapping for fruit size and shape in chromosomes 2 and 4 in pepper and a comparison of the pepper QTL map with that of tomato

Quantitative trait locus (QTL) mapping for fruit weight and shape in pepper (Capsicum spp.) was performed using C. chinense and C. frutescens introgression lines of chromosomes 2 and 4. In chromosome 2, a single major fruit-weight QTL, fw2.1, was detected in both populations that explained 62% of the trait variation. This QTL, as well as a fruit-shape QTL, fs2.1, which had a more minor effect, were localized to the tomato fruit-shape gene ovate. The cloned tomato fruit-weight QTL, fw2.2, did not play a major role in controlling fruit size variations in pepper. In chromosome 4, two fruit-weight QTLs, fw4.1 and fw4.2, were detected in the same genomic regions in both mapping populations. In addition, a single fruit-shape QTL was detected in each of the mapping populations that co-localized with one of the fruit-weight QTLs, suggesting pleiotropy or close linkage of the genes controlling size and shape. fw2.1 and fw4.2 represent major fruit-weight QTLs that are conserved in the three Capsicum species analyzed to date for fruit-size variations. Co-localization of the pepper QTLs with QTLs identified for similar traits in tomato suggests that the pepper and tomato QTLs are orthologous. Compared to fruit-shape QTLs, fruit-weight QTLs were more often conserved between pepper and tomato. This implies that different modes of selection were employed for these traits during domestication of the two Solanaceae species.S. Zygier and A. Ben Chaim contributed equally to this work.     

Ecology and genetics of speciation in Ficedula flycatchers

Birds have for long been popular study objects in speciation research. Being easy to observe in the field, they have traditionally been particularly important in studies of behavioural and ecological factors in speciation, whereas the genetic aspects of the process have been studied in other organisms, such as Drosophila. More recently, however, a stronger genetic focus has been placed on speciation research also in birds. Here, we review ecological, behavioural and genetic studies on speciation in the pied flycatcher (Ficedula hypoleuca) and the collared flycatcher (Ficedula albicollis). These well‐studied birds provide among the few proposed examples of the process of reinforcement of premating isolation, and the evidence for reinforcement is strong. They are further characterized by having strong intrinsic postzygotic barriers (female hybrid sterility), yet the two species appear to be very similar ecologically. This is in stark contrast to another well‐studied bird complex, Darwin’s finches, in which the species differ vastly in ecologically important traits but have no developmental problems arising from genetic incompatibilities, and where no evidence for reinforcement is found. In the flycatchers, sex chromosome linkage of genes affecting traits associated with both pre‐ and postzygotic barriers to gene exchange is likely to facilitate reinforcement. We discuss whether such sex‐linkage may be common in birds. The contrast between flycatchers and Darwin’s finches indicate that speciation in birds cannot always be understood mainly as a result of divergent natural selection (‘ecological speciation’), and generalizations from one system may lead us astray. We discuss to what extent insight from research on the flycatchers may point to fruitful avenues for future research on bird speciation and specifically call for a more systematic effort to simultaneously investigate ecology, behaviour and genetics of birds caught in the process of speciation.     

Use of locally weighted scatterplot smoothing (LOWESS) regression to study selection signatures in Piedmontese and Italian Brown cattle breeds

Selection is the major force affecting local levels of genetic variation in species. The availability of dense marker maps offers new opportunities for a detailed understanding of genetic diversity distribution across the animal genome. Over the last 50 years, cattle breeds have been subjected to intense artificial selection. Consequently, regions controlling traits of economic importance are expected to exhibit selection signatures. The fixation index (F_st) is an estimate of population differentiation, based on genetic polymorphism data, and it is calculated using the relationship between inbreeding and heterozygosity. In the present study, locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach were used to investigate selection signatures in two cattle breeds with different production aptitudes (dairy and beef). F_st was calculated for 42 514 SNP marker loci distributed across the genome in 749 Italian Brown and 364 Piedmontese bulls. The statistical significance of F_st values was assessed using a control chart. The LOWESS technique was efficient in removing noise from the raw data and was able to highlight selection signatures in chromosomes known to harbour genes affecting dairy and beef traits. Examples include the peaks detected for BTA2 in the region where the myostatin gene is located and for BTA6 in the region harbouring the ABCG2 locus. Moreover, several loci not previously reported in cattle studies were detected.