The genetic architecture of domestication in the chicken: effects of pleiotropy and linkage

Domestication affects sex-biased gene expression evolution in the duck

Genes with sex-biased expression are thought to underlie sexually dimorphic phenotypes and are therefore subject to different selection pressures in males and females. Many authors have proposed that sexual conflict leads to the evolution of sex-biased expression, which allows males and females to reach separate phenotypic and fitness optima. The selection pressures associated with domestication may cause changes in population architectures and mating systems, which in turn can alter their direction and strength. We compared sex-biased expression and genetic signatures in wild and domestic ducks (Anas platyrhynchos), and observed changes of sexual selection and identified the genomic divergence affected by selection forces. The extent of sex-biased expression in both sexes is positively correlated with the level of both d _N /d _S and nucleotide diversity. This observed changing pattern may mainly be owing to relaxed genetic constraints. We also demonstrate a clear link between domestication and sex-biased evolutionary rate in a comparative framework. Decreased polymorphism and evolutionary rate in domesticated populations generally matched life-history phenotypes known to experience artificial selection. Taken together, our work suggests the important implications of domestication in sex-biased evolution and the roles of artificial selection and sexual selection for shaping the diversity and evolutionary rate of the genome.

Genome-Wide Association Study Reveals the Genetic Basis of Duck Plumage Colors

Plumage color is an artificially and naturally selected trait in domestic ducks. Black, white, and spotty are the main feather colors in domestic ducks. Previous studies have shown that black plumage color is caused by MC1R, and white plumage color is caused by MITF. We performed a genome-wide association study (GWAS) to identify candidate genes associated with white, black, and spotty plumage in ducks. Two non-synonymous SNPs in MC1R (c.52G>A and c.376G>A) were significantly related to duck black plumage, and three SNPs in MITF (chr13:15411658A>G, chr13:15412570T>C and chr13:15412592C>G) were associated with white plumage. Additionally, we also identified the epistatic interactions between causing loci. Some ducks with white plumage carry the c.52G>A and c.376G>A in MC1R, which also compensated for black and spotty plumage color phenotypes, suggesting that MC1R and MITF have an epistatic effect. The MITF locus was supposed to be an upstream gene to MC1R underlying the white, black, and spotty colors. Although the specific mechanism remains to be further clarified, these findings support the importance of epistasis in plumage color variation in ducks.

Population structure and hybridisation in a population of Hawaiian feral chickens

Chickens are believed to have inhabited the Hawaiian island of Kauai since the first human migrations around 1200AD, but numbers have peaked since the tropical storms Iniki and Iwa in the 1980s and 1990s that destroyed almost all the chicken coops on the island and released large numbers of domestic chickens into the wild. Previous studies have shown these now feral chickens are an admixed population between Red Junglefowl (RJF) and domestic chickens. Here, using genetic haplotypic data, we estimate the time of the admixture event between the feral population on the island and the RJF to 1981 (1976-1995), coinciding with the timings of storm Iwa and Iniki. Analysis of genetic structure reveals a greater similarity between individuals inhabiting the northern and western part of the island to RJF than individuals from the eastern part of the island. These results point to the possibility of introgression events between feral chickens and the wild chickens in areas surrounding the Koke'e State Park and the Alaka'i plateau, posited as two of the major RJF reservoirs in the island. Furthermore, we have inferred haplotype blocks from pooled data to determine the most plausible source of the feral population. We identify a clear contribution from RJF and layer chickens of the White Leghorn (WL) breed. This work provides independent confirmation of the traditional hypothesis surrounding the origin of the feral populations and draws attention to the possibility of introgression of domestic alleles into the wild reservoir.

Effects of genotype, sex, and feed restriction on the biochemical composition of chicken preen gland secretions and their implications for commercial poultry production

Preen gland secretions spread on the feathers contain various chemical compounds dominated by fatty acids (FAs) and volatile organic compounds (VOCs). These chemicals may significantly affect plumage condition, microbial and ectoparasitic load on feathers, and chemical communication of birds. However, how chemical composition of preen secretions varies in commercially produced chickens with respect to their genotype, sex, and feeding regime remain largely unknown, as well as the welfare implications for farmed poultry. We found that while polyunsaturated fatty acids in chicken preen secretions differed significantly with genotype (P << 0.001), saturated fatty acids and monounsaturated fatty acids varied with genotype-dependent preen gland volume (P < 0.01). Chickens of meat-type fast-growing Ross 308 genotype had reduced preen gland volume and lower proportions of all FA categories in their preen secretions compared with dual-purpose slow-growing ISA Dual chickens. A total of 34 FAs and 77 VOCs with tens of unique FAs were detected in preen secretions of both genotypes. While differences in the relative proportion of 6 of the 10 most dominant VOCs in chicken preen gland secretions were related to genotype (P < 0.001), only 1 of the 10 most dominant VOCs showed a sex effect (P < 0.01), and only 2 of the 10 most dominant VOCs showed a genotype-dependent effect of feed restriction (P < 0.05). Feed restriction had no effect on the relative proportion of any of the FAs in chicken preen gland secretions. Moreover, we found that meat-type Ross 308 preen secretions were dominated by VOCs, which are proven attractants for poultry red mite and may also increase infestation with other ectoparasites and negatively influence overall odor-mediated intraspecific communication and welfare. This study shows that no feeding management, but long-term genetic selection in commercial breeding may be the main cause of the differences in the biochemistry and function of chicken preen secretions. This might have negative consequences for chemosignaling, antiparasitic, and antimicrobial potential of preen secretions and can lead to increased susceptibility to ectoparasites, plumage care disorders, and can affect the overall condition, welfare, and productivity of commercially bred chickens. Selection-induced preen gland impairments must therefore be considered and compensated by proper management of the chicken farm and increased care about animal well-being.

Genomic and gene expression associations to morphology of a sexual ornament in the chicken

How sexual selection affects the genome ultimately relies on the strength and type of selection, and the genetic architecture of the involved traits. While associating genotype with phenotype often utilizes standard trait morphology, trait representations in morphospace using geometric morphometric approaches receive less focus in this regard. Here, we identify genetic associations to a sexual ornament, the comb, in the chicken system (Gallus gallus). Our approach combined genome-wide genotype and gene expression data (>30k genes) with different aspects of comb morphology in an advanced intercross line (F8) generated by crossing a wild-type Red Junglefowl with a domestic breed of chicken (White Leghorn). In total, 10 quantitative trait loci were found associated to various aspects of comb shape and size, while 1,184 expression QTL were found associated to gene expression patterns, among which 98 had overlapping confidence intervals with those of quantitative trait loci. Our results highlight both known genomic regions confirming previous records of a large effect quantitative trait loci associated to comb size, and novel quantitative trait loci associated to comb shape. Genes were considered candidates affecting comb morphology if they were found within both confidence intervals of the underlying quantitative trait loci and eQTL. Overlaps between quantitative trait loci and genome-wide selective sweeps identified in a previous study revealed that only loci associated to comb size may be experiencing on-going selection under domestication.

A milk-sharing economy allows placental mammals to overcome their metabolic limits

Here, we demonstrate that a naturally evolving behavior (allonursing) has greater effect on reproductive power (mass per unit of time) and output (litter mass at birth) than does artificial selection (domestication). Additionally, we demonstrate the importance of resource optimization afforded by sociality (rather than resource abundance per se) in shaping a species’ life history profile and its ability to overcome its own physiological constraints.

Maternal resource availability and metabolism have a strong limiting effect on reproductive output. Allomaternal care and domestication increase the energy available to the mother and should correlate with an increase in reproductive output. Here, we take a comparative approach to understand how this increase is accomplished (e.g., litter mass, reproductive frequency, etc.) and the strength of the effect among different forms of external energetic supplementation. We find that domestication and all forms of allocare correlate with increased fertility. All forms of provisioning correlate with larger litters without compromising offspring size. The greatest increase we observe in reproductive power is in species that practice allonursing. Our results suggest that the ultimate factor limiting reproductive output in placental mammals is maternal metabolic power rather than resource availability.

A combination of genome-wide association study and selection signature analysis dissects the genetic architecture underlying bone traits in chickens

The bones of chicken play an important role in supporting and protecting the body. The growth and development of bones have a substantial influence on the health and production performance in chickens. However, genetic architecture underlying chicken bone traits is not well understood. The objectives of this study are to dissect the genetic basis of bone traits in chickens and to identify valuable genes and genetic markers for chicken breeding. We performed a combination of genome-wide association study (GWAS) and selection signature analysis (fixation index values and nucleotide diversity ratios) in an F₂ crossbred experimental population with different genetic backgrounds (broiler × layer) to identify candidate genes and significant variants related to femur, shank, keel length, chest width, metatarsal claw weight, metatarsal length, and metatarsal circumference. A total of 545 individuals were genotyped based on the whole genome re-sequencing method (26 F₀ individuals were re-sequenced at 10 × coverage; 519 F₂ individuals were re-sequenced at 3 × coverage). A total of 2 028 112 single-nucleotide polymorphisms (SNPs) remained to carry out analysis after quality control and imputation. The integration of GWAS and selection signature analysis indicated that all significant SNPs responsible for bone traits were mainly localized on chicken chromosomes 1, 4, and 27. Finally, we identified 21 positional candidate genes that might regulate chicken bone growth and development, including LRCH1, RB1, FNDC3A, MLNR, CAB39L, FOXO1, LHFP, TRPC4, POSTN, SMAD9, RBPJ, PPARGC1A, SLIT2, NCAPG, NKX3-2, CPZ, SPOP, NGFR, SOST, ZNF652, and HOXB3. Additionally, an array of uncharacterized genes was identified. The findings provide an in-depth understanding of the genetic architecture of chicken bone traits and offer a molecular basis for applying genomics in practical chicken breeding.

Genetic basis of chicken plumage color in artificial population of complex epistasis

Chicken plumage color, the genetic basis of which is often affected by epistasis, has long interested scientists. In the current study, a population of complex epistasis was constructed by crossing dominant White Leghorn chickens with recessive white feather chickens. Through a genome-wide association study, we identified single nucleotide polymorphisms and genes significantly associated with white and colored plumage in hens at different developmental stages. Interestingly, white plumage in adulthood was associated with the recessive white feather gene (TYR), whereas white feathers at birth stage were associated with the dominant white feather gene (PMEL), indicating age-related roles for these genes. TYR was shown to exert an epistatic effect on PMEL in adult hens. Additionally, TYR had an epistatic effect on barred plumage, while barred plumage had an epistatic effect on black plumage. TYR had no epistatic effect on the yellow plumage. We confirmed that the barred plumage gene is CDKN2A, as reported in previous studies. Golgb1 and REEP3, which play important roles in the Golgi network and affect the formation of feather pigments, are important candidate genes for yellow plumage. The candidate genes for black plumage are CAMKK1 and IFT22. Further research is warranted to elucidate the molecular mechanisms underlying these traits.

Genomic Analyses of Unveil Helmeted Guinea Fowl (Numida meleagris) Domestication in West Africa

Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300-5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.

CREBBP and WDR 24 Identified as Candidate Genes for Quantitative Variation in Red-Brown Plumage Colouration in the Chicken

Plumage colouration in birds is important for a plethora of reasons, ranging from camouflage, sexual signalling, and species recognition. The genes underlying colour variation have been vital in understanding how genes can affect a phenotype. Multiple genes have been identified that affect plumage variation, but research has principally focused on major-effect genes (such as those causing albinism, barring, and the like), rather than the smaller effect modifier loci that more subtly influence colour. By utilising a domestic × wild advanced intercross with a combination of classical QTL mapping of red colouration as a quantitative trait and a targeted genetical genomics approach, we have identified five separate candidate genes (CREBBP, WDR24, ARL8A, PHLDA3, LAD1) that putatively influence quantitative variation in red-brown colouration in chickens. By treating colour as a quantitative rather than qualitative trait, we have identified both QTL and genes of small effect. Such small effect loci are potentially far more prevalent in wild populations, and can therefore potentially be highly relevant to colour evolution.

Artificial selection for male winners in the Siamese fighting fish Betta splendens correlates with high female aggression

In Southeast Asia, males of the Siamese fighting fish Betta splendens have been selected across centuries for paired-staged fights. During the selection process, matched for size males fight in a small tank until the contest is resolved. Breeders discard losing batches and reproduce winner batches with the aim of increasing fight performance. We assessed the results of this long-term selection process by comparing under standard laboratory conditions male and female aggressive behaviour of one strain selected for staged fights ("fighters") and one strain of wild-types. The aggressive response of adult fish was tested against their mirror image or a size-matched conspecific. Fighter males were more aggressive than wild-type males for all measured behaviours. Differences were not only quantitative but the pattern of fight display was also divergent. Fighter males had an overall higher swimming activity, performing frequent fast strikes in the direction of the intruder and displaying from a distance. Wild-type males were less active and exhibited aggressive displays mostly in close proximity to the stimuli. Females of the fighter strain, which are not used for fights, were also more aggressive than wild-type females. Aggressive behaviours were correlated across male and female fighter siblings, suggesting common genetic and physiological mechanisms to male and female aggression in this species. The study further shows that results were largely independent of the stimulus type, with the mirror test inducing similar and less variable responses than the live conspecific presentation. These results suggest that selection for male winners co-selected for high-frequency and metabolic demanding aggressive display in males and also enhanced female aggression, opening a wide range of testable hypothesis about the ultimate and proximate mechanisms of male and female aggression in B. splendens.

Maladaptation in feral and domesticated animals

Selection regimes and population structures can be powerfully changed by domestication and feralization, and these changes can modulate animal fitness in both captive and natural environments. In this review, we synthesize recent studies of these two processes and consider their impacts on organismal and population fitness. Domestication and feralization offer multiple windows into the forms and mechanisms of maladaptation. Firstly, domestic and feral organisms that exhibit suboptimal traits or fitness allow us to identify their underlying causes within tractable research systems. This has facilitated significant progress in our general understandings of genotype-phenotype relationships, fitness trade-offs, and the roles of population structure and artificial selection in shaping domestic and formerly domestic organisms. Additionally, feralization of artificially selected gene variants and organisms can reveal or produce maladaptation in other inhabitants of an invaded biotic community. In these instances, feral animals often show similar fitness advantages to other invasive species, but they are also unique in their capacities to modify natural ecosystems through introductions of artificially selected traits. We conclude with a brief consideration of how emerging technologies such as genome editing could change the tempos, trajectories, and ecological consequences of both domestication and feralization. In addition to providing basic evolutionary insights, our growing understanding of mechanisms through which artificial selection can modulate fitness has diverse and important applications-from enhancing the welfare, sustainability, and efficiency of agroindustry, to mitigating biotic invasions.

Genomic Analysis Reveals Pleiotropic Alleles at EDN3 and BMP7 Involved in Chicken Comb Color and Egg Production

Artificial selection is often associated with numerous changes in seemingly unrelated phenotypic traits. The genetic mechanisms of correlated phenotypes probably involve pleiotropy or linkage of genes related to such phenotypes. Dongxiang blue-shelled chicken, an indigenous chicken breed of China, has segregated significantly for the dermal hyperpigmentation phenotype. Two lines of the chicken have been divergently selected with respect to comb color for over 20 generations. The red comb line chicken produces significantly higher number of eggs than the dark comb line chicken. The objective of this study was to explore potential mechanisms involved in the relationship between comb color and egg production among chickens. Based on the genome-wide association study results, we identified a genomic region on chromosome 20 involving EDN3 and BMP7, which is associated with hyperpigmentation of chicken comb. Further analyses by selection signatures in the two divergent lines revealed that several candidate genes, including EDN3, BMP7, BPIFB3, and PCK1, closely located on chromosome 20 are involved in the development of neural crest cell and reproductive system. The two genes EDN3 and BMP7 have known roles in regulating both ovarian function and melanogenesis, indicating the pleiotropic effect on hyperpigmentation and egg production in blue-shelled chickens. Association analysis for egg production confirmed the pleiotropic effect of selected loci identified by selection signatures. The study provides insights into phenotypic evolution due to genetic variation across the genome. The information might be useful in the current breeding efforts to develop improved breeds for egg production.

Selection for longer limbs in mice increases bone stiffness and brittleness, but does not alter bending strength

The ability of a bone to withstand loads depends on its structural and material properties. These tend to differ among species with different modes of locomotion, reflecting their unique loading patterns. The evolution of derived limb morphologies, such as the long limbs associated with jumping, may compromise overall bone strength. We evaluated bone mechanical properties in the Longshanks mouse, which was selectively bred for increased tibia length relative to body mass. We combined analyses of 3D shape and cross-sectional geometry of the tibia, with mechanical testing and bone composition assays, to compare bone strength, elastic properties and mineral composition in Longshanks mice and randomly bred controls. Our data show that, despite being more slender, cortical geometry and predicted bending strength of the Longshanks tibia were similar to controls. In whole bone bending tests, measures of bone bending strength were similar across groups; however, Longshanks tibiae were significantly more rigid, more brittle, and required less than half the energy to fracture. Tissue-level elastic properties were also altered in Longshanks mice, but the bones did not differ from the control in water content, ash content or density. These results indicate that while Longshanks bones are as strong as control tibiae, selection for increased tibia length has altered its elastic properties, possibly through changes in organic bony matrix composition. We conclude that selection for certain limb morphologies, and/or selection for rapid skeletal growth, can lead to tissue-level changes that can increase the risk of skeletal fracture, which in turn may favor the correlated evolution of compensatory mechanisms to mitigate increased fracture risk, such as delayed skeletal maturity.

Genetical Genomics of Tonic Immobility in the Chicken

Identifying the molecular mechanisms of animal behaviour is an enduring goal for researchers. Gaining insight into these mechanisms enables us to gain a greater understanding of behaviour and their genetic control. In this paper, we perform Quantitative Trait Loci (QTL) mapping of tonic immobility behaviour in an advanced intercross line between wild and domestic chickens. Genes located within the QTL interval were further investigated using global expression QTL (eQTL) mapping from hypothalamus tissue, as well as causality analysis. This identified five candidate genes, with the genes PRDX4 and ACOT9 emerging as the best supported candidates. In addition, we also investigated the connection between tonic immobility, meat pH and struggling behaviour, as the two candidate genes PRDX4 and ACOT9 have previously been implicated in controlling muscle pH at slaughter. We did not find any phenotypic correlations between tonic immobility, struggling behaviour and muscle pH in a smaller additional cohort, despite these behaviours being repeatable within-test.

Analysis of phenotypic- and Estimated Breeding Values (EBV) to dissect the genetic architecture of complex traits in a Scots pine three-generation pedigree design

In forest tree breeding, family-based Quantitative Trait Loci (QTL) studies are valuable as methods to dissect the complexity of a trait and as a source of candidate genes. In the field of conifer research, our study contributes to the evaluation of phenotypic and predicted breeding values for the identification of QTL linked to complex traits in a three-generation pedigree population in Scots pine (Pinus sylvestris L.). A total of 11 470 open pollinated F₂-progeny trees established at three different locations, were measured for growth and adaptive traits. Breeding values were predicted for their 360 mothers, originating from a single cross of two grand-parents. A multilevel LASSO association analysis was conducted to detect QTL using genotypes of the mothers with the corresponding phenotypes and Estimated Breeding Values (EBV). Different levels of genotype-by-environment (G × E) effects among sites at different years, were detected for survival and height. Moderate-to-low narrow sense heritabilities and EBV accuracies were found for all traits and all sites. We identified 18 AFLPs and 12 SNPs to be associated with QTL for one or more traits. 62 QTL were significant with percentages of variance explained ranging from 1.7 to 18.9%. In those cases where the same marker was associated to a phenotypic or an ebvQTL, the ebvQTL always explained higher proportion of the variance, maybe due to the more accurate nature of Estimated Breeding Values (EBV). Two SNP-QTL showed pleiotropic effects for traits related with hardiness, seed, cone and flower production. Furthermore, we detected several QTL with significant effects across multiple ages, which could be considered as strong candidate loci for early selection. The lack of reproducibility of some QTL detected across sites may be due to environmental heterogeneity reflected by the genotype- and QTL-by-environment effects.

Effective population sizes and adaptive genetic variation in a captive bird population

Captive populations are considered a key component of ex situ conservation programs. Research on multiple taxa has shown the differential success of maintaining demographic versus genetic stability and viability in captive populations. In typical captive populations, usually founded by few or related individuals, genetic diversity can be lost and inbreeding can accumulate rapidly, calling into question their ultimate utility for release into the wild. Furthermore, domestication selection for survival in captive conditions is another concern. Therefore, it is crucial to understand the dynamics of population sizes, particularly the effective population size, and genetic diversity at non-neutral and adaptive loci in captive populations. In this study, we assessed effective population sizes and genetic variation at both neutral microsatellite markers, as well as SNP variants from the MHC-B locus of a captive Red Junglefowl population. This population represents a rare instance of a population with a well-documented history in captivity, following a realistic scenario of chain-of-custody, unlike many captive lab populations. Our analyses, which included 27 individuals comprising the entirety of one captive population show very low neutral and adaptive genetic variation, as well as low effective sizes, which correspond with the known demographic history. Finally, our study also shows the divergent impacts of small effective size and inbreeding in captive populations on microsatellite versus adaptive genetic variation in the MHC-B locus. Our study provides insights into the difficulties of maintaining adaptive genetic variation in small captive populations.

Genetics and Genomics of Social Behavior in a Chicken Model

The identification of genes affecting sociality can give insights into the maintenance and development of sociality and personality. In this study, we used the combination of an advanced intercross between wild and domestic chickens with a combined QTL and eQTL genetical genomics approach to identify genes for social reinstatement, a social and anxiety-related behavior. A total of 24 social reinstatement QTL were identified and overlaid with over 600 eQTL obtained from the same birds using hypothalamic tissue. Correlations between overlapping QTL and eQTL indicated five strong candidate genes, with the gene TTRAP being strongly significantly correlated with multiple aspects of social reinstatement behavior, as well as possessing a highly significant eQTL.

QTL mapping of stress related gene expression in a cross between domesticated chickens and ancestral red junglefowl

Domestication of animals is associated with numerous alterations in physiology, morphology, and behavior. Lower reactivity of the hypothalamic-pituitary-adrenal (HPA) axis and reduced fearfulness is seen in most studied domesticates, including chickens. Previously we have shown that the physiological stress response as well as expression levels of hundreds of genes in the hypothalamus and adrenal glands are different between domesticated White Leghorn and the progenitor of modern chickens, the Red Junglefowl. To map genetic loci associated with the transcription levels of genes involved in the physiological stress response, we conducted an eQTL analysis in the F₁₂ generation of an inter-cross between White Leghorn and Red Junglefowl. We selected genes for further studies based on their known function in the regulation of the HPA axis or sympathoadrenal (SA) system, and measured their expression levels in the hypothalamus and the adrenal glands after a brief stress exposure (physical restraint). The expression values were treated as quantitative traits for the eQTL mapping. The plasma levels of corticosterone were also assessed. We analyzed the correlation between gene expression and corticosterone levels and mapped eQTL and their potential effects on corticosterone levels. The effects on gene transcription of a previously found QTL for corticosterone response were also investigated. The expression levels of the glucocorticoid receptor (GR) in the hypothalamus and several genes in the adrenal glands were correlated with the post-stress levels of corticosterone in plasma. We found several cis- and trans-acting eQTL for stress-related genes in both hypothalamus and adrenal. In the hypothalamus, one eQTL for c-FOS and one QTL for expression of GR were found. In the adrenal tissue, we identified eQTL for the genes NR0B1, RGS4, DBH, MAOA, GRIN1, GABRB2, GABRB3, and HSF1. None of the found eQTL were significant predictors of corticosterone levels. The previously found QTL for corticosterone was associated with GR expression in hypothalamus. Our data suggests that domestication related modification in the stress response is driven by changes in the transcription levels of several modulators of the HPA and SA systems in hypothalamus and adrenal glands and not by changes in the expression of the steroidogenic genes. The presence of eQTL for GR in hypothalamus combined with the negative correlation between GR expression and corticosterone response suggests GR as a candidate for further functional studies regarding modification of stress response during chicken domestication.

Genetic and Targeted eQTL Mapping Reveals Strong Candidate Genes Modulating the Stress Response During Chicken Domestication

The stress response has been largely modified in all domesticated animals, offering a strong tool for genetic mapping. In chickens, ancestral Red Junglefowl react stronger both in terms of physiology and behavior to a brief restraint stress than domesticated White Leghorn, demonstrating modified functions of the hypothalamic–pituitary–adrenal (HPA) axis. We mapped quantitative trait loci (QTL) underlying variations in stress-induced hormone levels using 232 birds from the 12th generation of an advanced intercross between White Leghorn and Red Junglefowl, genotyped for 739 genetic markers. Plasma levels of corticosterone, dehydroepiandrosterone (DHEA), and pregnenolone (PREG) were measured using LC-MS/MS in all genotyped birds. Transcription levels of the candidate genes were measured in the adrenal glands or hypothalamus of 88 out of the 232 birds used for hormone assessment. Genes were targeted for expression analysis when they were located in a hormone QTL region and were differentially expressed in the pure breed birds. One genome-wide significant QTL on chromosome 5 and two suggestive QTL together explained 20% of the variance in corticosterone response. Two significant QTL for aldosterone on chromosome 2 and 5 (explaining 19% of the variance), and one QTL for DHEA on chromosome 4 (explaining 5% of the variance), were detected. Orthologous DNA regions to the significant corticosterone QTL have been previously associated with the physiological stress response in other species but, to our knowledge, the underlying gene(s) have not been identified. SERPINA10 had an expression QTL (eQTL) colocalized with the corticosterone QTL on chromosome 5 and PDE1C had an eQTL colocalized with the aldosterone QTL on chromosome 2. Furthermore, in both cases, the expression levels of the genes were correlated with the plasma levels of the hormones. Hence, both these genes are strong putative candidates for the domestication-induced modifications of the stress response in chickens. Improved understanding of the genes associated with HPA-axis reactivity can provide insights into the pathways and mechanisms causing stress-related pathologies.

Feralisation targets different genomic loci to domestication in the chicken

Feralisation occurs when a domestic population recolonizes the wild, escaping its previous restricted environment, and has been considered as the reverse of domestication. We have previously shown that Kauai Island's feral chickens are a highly variable and admixed population. Here we map selective sweeps in feral Kauai chickens using whole-genome sequencing. The detected sweeps were mostly unique to feralisation and distinct to those selected for during domestication. To ascribe potential phenotypic functions to these genes we utilize a laboratory-controlled equivalent to the Kauai population-an advanced intercross between Red Junglefowl and domestic layer birds that has been used previously for both QTL and expression QTL studies. Certain sweep genes exhibit significant correlations with comb mass, maternal brooding behaviour and fecundity. Our analyses indicate that adaptations to feral and domestic environments involve different genomic regions and feral chickens show some evidence of adaptation at genes associated with sexual selection and reproduction.

Genetical Genomics of Behavior: A Novel Chicken Genomic Model for Anxiety Behavior

The identification of genetic variants responsible for behavioral variation is an enduring goal in biology, with wide-scale ramifications, ranging from medical research to evolutionary theory on personality syndromes. Here, we use for the first time a large-scale genetical genomics analysis in the brains of chickens to identify genes affecting anxiety as measured by an open field test. We combine quantitative trait locus (QTL) analysis in 572 individuals and expression QTL (eQTL) analysis in 129 individuals from an advanced intercross between domestic chickens and Red Junglefowl. We identify 10 putative quantitative trait genes affecting anxiety behavior. These genes were tested for an association in the mouse Heterogeneous Stock anxiety (open field) data set and human GWAS data sets for bipolar disorder, major depressive disorder, and schizophrenia. Although comparisons between species are complex, associations were observed for four of the candidate genes in mice and three of the candidate genes in humans. Using a multimodel approach we have therefore identified a number of putative quantitative trait genes affecting anxiety behavior, principally in chickens but also with some potentially translational effects as well. This study demonstrates that chickens are an excellent model organism for the genetic dissection of behavior.

Factors influencing the effect size distribution of adaptive substitutions

The distribution of effect sizes of adaptive substitutions has been central to evolutionary biology since the modern synthesis. Early theory proposed that because large-effect mutations have negative pleiotropic consequences, only small-effect mutations contribute to adaptation. More recent theory suggested instead that large-effect mutations could be favoured when populations are far from their adaptive peak. Here we suggest that the distributions of effect sizes are expected to differ among study systems, reflecting the wide variation in evolutionary forces and ecological conditions experienced in nature. These include selection, mutation, genetic drift, gene flow, and other factors such as the degree of pleiotropy, the distance to the phenotypic optimum, whether the optimum is stable or moving, and whether new mutation or standing genetic variation provides the source of adaptive alleles. Our goal is to review how these factors might affect the distribution of effect sizes and to identify new research directions. Until more theory and empirical work is available, we feel that it is premature to make broad generalizations about the effect size distribution of adaptive substitutions important in nature.

Quantitative Trait Locus and Genetical Genomics Analysis Identifies Putatively Causal Genes for Fecundity and Brooding in the Chicken

Life history traits such as fecundity are important to evolution because they make up components of lifetime fitness. Due to their polygenic architectures, such traits are difficult to investigate with genetic mapping. Therefore, little is known about their molecular basis. One possible way toward finding the underlying genes is to map intermediary molecular phenotypes, such as gene expression traits. We set out to map candidate quantitative trait genes for egg fecundity in the chicken by combining quantitative trait locus mapping in an advanced intercross of wild by domestic chickens with expression quantitative trait locus mapping in the same birds. We measured individual egg fecundity in 232 intercross chickens in two consecutive trials, the second one aimed at measuring brooding. We found 12 loci for different aspects of egg fecundity. We then combined the genomic confidence intervals of these loci with expression quantitative trait loci from bone and hypothalamus in the same intercross. Overlaps between egg loci and expression loci, and trait–gene expression correlations identify 29 candidates from bone and five from hypothalamus. The candidate quantitative trait genes include fibroblast growth factor 1, and mitochondrial ribosomal proteins L42 and L32. In summary, we found putative quantitative trait genes for egg traits in the chicken that may have been affected by regulatory variants under chicken domestication. These represent, to the best of our knowledge, some of the first candidate genes identified by genome-wide mapping for life history traits in an avian species.

Domestication Effects on Stress Induced Steroid Secretion and Adrenal Gene Expression in Chickens

Understanding the genetic basis of phenotypic diversity is a challenge in contemporary biology. Domestication provides a model for unravelling aspects of the genetic basis of stress sensitivity. The ancestral Red Junglefowl (RJF) exhibits greater fear-related behaviour and a more pronounced HPA-axis reactivity than its domesticated counterpart, the White Leghorn (WL). By comparing hormones (plasmatic) and adrenal global gene transcription profiles between WL and RJF in response to an acute stress event, we investigated the molecular basis for the altered physiological stress responsiveness in domesticated chickens. Basal levels of pregnenolone and dehydroepiandrosterone as well as corticosterone response were lower in WL. Microarray analysis of gene expression in adrenal glands showed a significant breed effect in a large number of transcripts with over-representation of genes in the channel activity pathway. The expression of the best-known steroidogenesis genes were similar across the breeds used. Transcription levels of acute stress response genes such as StAR, CH25 and POMC were upregulated in response to acute stress. Dampened HPA reactivity in domesticated chickens was associated with changes in the expression of several genes that presents potentially minor regulatory effects rather than by means of change in expression of critical steroidogenic genes in the adrenal.

The Genetic Architecture of Domestication in Animals

Domestication has been essential to the progress of human civilization, and the process itself has fascinated biologists for hundreds of years. Domestication has led to a series of remarkable changes in a variety of plants and animals, in what is termed the “domestication phenotype.” In domesticated animals, this general phenotype typically consists of similar changes in tameness, behavior, size/morphology, color, brain composition, and adrenal gland size. This domestication phenotype is seen in a range of different animals. However, the genetic basis of these associated changes is still puzzling. The genes for these different traits tend to be grouped together in clusters in the genome, though it is still not clear whether these clusters represent pleiotropic effects, or are in fact linked clusters. This review focuses on what is currently known about the genetic architecture of domesticated animal species, if genes of large effect (often referred to as major genes) are prevalent in driving the domestication phenotype, and whether pleiotropy can explain the loci underpinning these diverse traits being colocated.

Genetic regulation of bone metabolism in the chicken: similarities and differences to Mammalian systems

Birds have a unique bone physiology, due to the demands placed on them through egg production. In particular their medullary bone serves as a source of calcium for eggshell production during lay and undergoes continuous and rapid remodelling. We take advantage of the fact that bone traits have diverged massively during chicken domestication to map the genetic basis of bone metabolism in the chicken. We performed a quantitative trait locus (QTL) and expression QTL (eQTL) mapping study in an advanced intercross based on Red Junglefowl (the wild progenitor of the modern domestic chicken) and White Leghorn chickens. We measured femoral bone traits in 456 chickens by peripheral computerised tomography and femoral gene expression in a subset of 125 females from the cross with microarrays. This resulted in 25 loci for female bone traits, 26 loci for male bone traits and 6318 local eQTL loci. We then overlapped bone and gene expression loci, before checking for an association between gene expression and trait values to identify candidate quantitative trait genes for bone traits. A handful of our candidates have been previously associated with bone traits in mice, but our results also implicate unexpected and largely unknown genes in bone metabolism. In summary, by utilising the unique bone metabolism of an avian species, we have identified a number of candidate genes affecting bone allocation and metabolism. These findings can have ramifications not only for the understanding of bone metabolism genetics in general, but could also be used as a potential model for osteoporosis as well as revealing new aspects of vertebrate bone regulation or features that distinguish avian and mammalian bone.

In this work we seek to further the understanding of bone genetics by mapping bone traits and gene expression in the chicken. Bone in female birds is special due to egg production. In this study, we combine the genetic mapping of bone traits with bone gene expression to find candidate quantitative trait genes that explain the differences between wild and domestic chickens in terms of bone production. The concept of combining genetic mapping and gene expression mapping is not new, and has already been successful in isolating bone-related genes in mammals, however this is the first time it has been applied to an avian system with such unique bone modelling processes. We aim to reveal new molecular mechanisms of bone regulation, and many of the candidates we find are new, highlighting the potential this technique has to identify the potential differences between avian and mammalian bone biology.

A medium density genetic map and QTL for behavioral and production traits in Japanese quail

BACKGROUND

Behavioral traits such as sociability, emotional reactivity and aggressiveness are major factors in animal adaptation to breeding conditions. In order to investigate the genetic control of these traits as well as their relationships with production traits, a study was undertaken on a large second generation cross (F2) between two lines of Japanese Quail divergently selected on their social reinstatement behavior. All the birds were measured for several social behaviors (social reinstatement, response to social isolation, sexual motivation, aggression), behaviors measuring the emotional reactivity of the birds (reaction to an unknown object, tonic immobility reaction), and production traits (body weight and egg production).

RESULTS

We report the results of the first genome-wide QTL detection based on a medium density SNP panel obtained from whole genome sequencing of a pool of individuals from each divergent line. A genetic map was constructed using 2145 markers among which 1479 could be positioned on 28 different linkage groups. The sex-averaged linkage map spanned a total of 3057 cM with an average marker spacing of 2.1 cM. With the exception of a few regions, the marker order was the same in Japanese Quail and the chicken, which confirmed a well conserved synteny between the two species. The linkage analyses performed using QTLMAP software revealed a total of 45 QTLs related either to behavioral (23) or production (22) traits. The most numerous QTLs (15) concerned social motivation traits. Interestingly, our results pinpointed putative pleiotropic regions which controlled emotional reactivity and body-weight of birds (on CJA5 and CJA8) or their social motivation and the onset of egg laying (on CJA19).

CONCLUSION

This study identified several QTL regions for social and emotional behaviors in the Quail. Further research will be needed to refine the QTL and confirm or refute the role of candidate genes, which were suggested by bioinformatics analysis. It can be hoped that the identification of genes and polymorphisms related to behavioral traits in the quail will have further applications for other poultry species (especially the chicken) and will contribute to solving animal welfare issues in poultry production.

The strong selective sweep candidate gene ADRA2C does not explain domestication related changes in the stress response of chickens

Analysis of selective sweeps to pinpoint causative genomic regions involved in chicken domestication has revealed a strong selective sweep on chromosome 4 in layer chickens. The autoregulatory α-adrenergic receptor 2C (ADRA2C) gene is the closest to the selective sweep and was proposed as an important gene in the domestication of layer chickens. The ADRA2C promoter region was also hypermethylated in comparison to the non-selected ancestor of all domesticated chicken breeds, the Red Junglefowl, further supporting its relevance. In mice the receptor is involved in the fight-or-flight response as it modulates epinephrine release from the adrenals. To investigate the involvement of ADRA2C in chicken domestication, we measured gene expression in the adrenals and radiolabeled receptor ligand in three brain regions comparing the domestic White Leghorn strain with the wild ancestor Red Junglefowl. In adrenals ADRA2C was twofold greater expressed than the related receptor gene ADRA2A, indicating that ADRA2C is the predominant modulator of epinephrine release but no strain differences were measured. In hypothalamus and amygdala, regions associated with the stress response, and in striatum, receptor binding pIC50 values ranged between 8.1–8.4, and the level was not influenced by the genotyped allele. Because chicken strains differ in morphology, physiology and behavior, differences attributed to a single gene may be lost in the noise caused by the heterogeneous genetic background. Therefore an F₁₀ advanced intercross strain between White Leghorn and Red Junglefowl was used to investigate effects of ADRA2C alleles on fear related behaviors and fecundity. We did not find compelling genotype effects in open field, tonic immobility, aerial predator, associative learning or fecundity. Therefore we conclude that ADRA2C is probably not involved in the domestication of the stress response in chicken, and the strong selective sweep is probably caused by selection of some unknown genetic element in the vicinity of the gene.

Evidence of phenotypic and genetic relationships between sociality, emotional reactivity and production traits in Japanese quail

The social behavior of animals, which is partially controlled by genetics, is one of the factors involved in their adaptation to large breeding groups. To understand better the relationships between different social behaviors, fear behaviors and production traits, we analyzed the phenotypic and genetic correlations of these traits in Japanese quail by a second generation crossing of two lines divergently selected for their social reinstatement behavior. Analyses of results for 900 individuals showed that the phenotypic correlations between behavioral traits were low with the exception of significant correlations between sexual behavior and aggressive pecks both at phenotypic (0.51) and genetic (0.90) levels. Significant positive genetic correlations were observed between emotional reactivity toward a novel object and sexual (0.89) or aggressive (0.63) behaviors. The other genetic correlations were observed mainly between behavioral and production traits. Thus, the level of emotional reactivity, estimated by the duration of tonic immobility, was positively correlated with weight at 17 and 65 days of age (0.76 and 0.79, respectively) and with delayed egg laying onset (0.74). In contrast, a higher level of social reinstatement behavior was associated with an earlier egg laying onset (-0.71). In addition, a strong sexual motivation was correlated with an earlier laying onset (-0.68) and a higher number of eggs laid (0.82). A low level of emotional reactivity toward a novel object and also a higher aggressive behavior were genetically correlated with a higher number of eggs laid (0.61 and 0.58, respectively). These results bring new insights into the complex determinism of social and emotional reactivity behaviors in birds and their relationships with production traits. Furthermore, they highlight the need to combine animal welfare and production traits in selection programs by taking into account traits of sociability and emotional reactivity.

Detection of molecular signatures of selection at microsatellite loci in the South African abalone (Haliotis midae) using a population genomic approach

Identifying genomic regions that may be under selection is important for elucidating the genetic architecture of complex phenotypes underlying adaptation to heterogeneous environments. A population genomic approach, using a classical neutrality test and various Fst-outlier detection methods was employed to evaluate genome-wide polymorphism data in order to identify loci that may be candidates for selection amongst six populations (three cultured and three wild) of the South African abalone, Haliotis midae. Approximately 9% of the genome-wide microsatellite markers were putatively subject to directional selection, whilst 6-18% of the genome is thought to be influenced by balancing selection. Genetic diversity estimates for candidate loci under directional selection was significantly reduced in comparison to candidate neutral loci, whilst candidate balancing selection loci demonstrated significantly higher levels of genetic diversity (Kruskal-Wallis test, P<0.05). Pairwise Fst estimates based on candidate directional selection loci also demonstrated increased levels of differentiation between study populations. Various candidate loci under selection showed significant inter-chromosomal linkage disequilibrium, suggesting possible gene-networks underling adaptive phenotypes. Furthermore, several loci had significant hits to known genes when performing BLAST searches to NCBI's non-redundant databases, whilst others are known to be derived from expressed sequences even though homology to a known gene could not be established. A number of loci also demonstrated relatively high similarity to transposable elements. The association of these loci to functional and genomically active sequences could in part explain the observed signatures of selection.

QTL clustering as a mechanism for rapid multi-trait evolution

Cave-dwelling animals exhibit remarkable convergence in multiple cave-related traits, yet the genetic mechanisms responsible for the evolution and integration of many such traits remain unclear. Astyanax mexicanus is a model cave-dwelling fish with sighted surface-dwelling (surface fish) and blind cave-dwelling (cavefish) forms. Using a genetic cross between these morphs, we discovered significant correlations among several cave-related traits, including reduced eyes, increased superficial neuromast receptors located within the eye orbit (EO SN) and a vibration-attraction behavior (VAB) that facilitates foraging in darkness. Furthermore, we discovered that the quantitative trait loci (QTL) underlying these traits are clustered within the Astyanax genome. Following an ablation experiment that demonstrated that the EO SN contribute to VAB, we concluded that the adaptive evolution of VAB and EO SN has likely contributed to eye loss in cavefish. In this addendum, we further discuss the possible role of multi-trait QTL clustering in facilitating rapid adaptation.

Domestication-related variation in social preferences in chickens is affected by genotype on a growth QTL

A growth-related QTL on chicken chromosome 1 has previously been shown to influence domestication behaviour in chickens. In this study, we used Red Junglefowl (RJF) and White Leghorn (WL) as well as the intercross between them to investigate whether stress affects the way birds allocate their time between familiar and unfamiliar conspecifics in a social preference test ('social support seeking'), and how this is related to genotype at specific loci within the growth QTL. Red Junglefowl males spent more time with unfamiliar chickens before the stressful event compared to the other birds, whereas all birds except WL males tended to spend less time with unfamiliar ones after stress. A significant QTL locus was found to influence both social preference under undisturbed circumstances and social support seeking. The WL allele at this QTL was associated not only with a preference for unfamiliar individuals but also with a shift towards familiar ones in response to stress (social support seeking). A second, suggestive QTL also affected social support seeking, but in the opposite direction; the WL allele was associated with increased time spent with unfamiliar individuals. The region contains several possible candidate genes, and gene expression analysis of a number of them showed differential expression between RJF and WL of AVPR2 (receptor for vasotocin), and possibly AVPR1a (another vasotocin receptor) and NRCAM (involved in neural development) in the lower frontal lobes of the brains of RJF and WL animals. These three genes continue to be interesting candidates for the observed behavioural effects.

Localizing F(ST) outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation-with-gene-flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence-with-gene-flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F(ST) outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between-population recombination and gene exchange are reduced over large 'divergence hitchhiking' (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F(ST) outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation-with-gene-flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.

A sexual ornament in chickens is affected by pleiotropic alleles at HAO1 and BMP2, selected during domestication

Domestication is one of the strongest forms of short-term, directional selection. Although selection is typically only exerted on one or a few target traits, domestication can lead to numerous changes in many seemingly unrelated phenotypes. It is unknown whether such correlated responses are due to pleiotropy or linkage between separate genetic architectures. Using three separate intercrosses between wild and domestic chickens, a locus affecting comb mass (a sexual ornament in the chicken) and several fitness traits (primarily medullary bone allocation and fecundity) was identified. This locus contains two tightly-linked genes, BMP2 and HAO1, which together produce the range of pleiotropic effects seen. This study demonstrates the importance of pleiotropy (or extremely close linkage) in domestication. The nature of this pleiotropy also provides insights into how this sexual ornament could be maintained in wild populations.

The genetic analysis of phenotypes and the identification of the causative underlying genes remain central to molecular and evolutionary biology. By utilizing the domestication process, it is possible to exploit the large differences between domesticated animals and their wild counterparts to study both this and the mechanism of domestication itself. Domestication has been central to the advent of modern civilization; and yet, despite domesticated animals displaying similar adaptations in morphology, physiology, and behaviour, the genetic basis of these changes are unknown. In addition, though sexual selection theory has been the subject of a vast amount of study, very little is known about which genes are underpinning such traits. We have generated multiple intercrosses and advanced intercrosses based on wild-derived and domestic chickens to fine-map genomic regions affecting a sexual ornament. These regions have been over-laid with putative selective sweeps identified in domestic chickens and found to be significantly associated with them. By using expression QTL analysis, we show that two genes in one region, HAO1 and BMP2, are controlling multiple aspects of the domestication phenotype, from a sexual ornament to multiple life history traits. This demonstrates the importance of pleiotropy (or extremely close linkage) in controlling these genetic changes.

Genetic divergence and the genetic architecture of complex traits in chromosome substitution strains of mice

Background

The genetic architecture of complex traits strongly influences the consequences of inherited mutations, genetic engineering, environmental and genetic perturbations, and natural and artificial selection. But because most studies are under-powered, the picture of complex traits is often incomplete. Chromosome substitution strains (CSSs) are a unique paradigm for these genome surveys because they enable statistically independent, powerful tests for the phenotypic effects of each chromosome on a uniform inbred genetic background. A previous CSS survey in mice and rats revealed many complex trait genes (QTLs), large phenotypic effects, extensive epistasis, as well as systems properties such as strongly directional phenotypic changes and genetically-determined limits on the range of phenotypic variation. However, the unusually close genetic relation between the CSS progenitor strains in that study raised questions about the impact of genetic divergence: would greater divergence between progenitor strains, with the corresponding changes in gene regulation and protein function, lead to significantly more distinctive phenotypic features, or alternatively would epistasis and systems constraints, which are pervasive in CSSs, limit the range of phenotypic variation regardless of the extent of DNA sequence variation?

Results

We analyzed results for an extensive survey of traits in two new panels of CSSs where the donor strains were derived from inbred strains with more distant origins and discovered a strong similarity in genetic and systems properties among the three CSS panels, regardless of divergence time.

Conclusion

Our results argue that DNA sequence differences between host and donor strains did not substantially affect the architecture of complex traits, and suggest instead that strong epistasis buffered the phenotypic effects of genetic divergence, thereby constraining the range of phenotypic variation.

Divergence hitchhiking and the spread of genomic isolation during ecological speciation-with-gene-flow

In allopatric populations, geographical separation simultaneously isolates the entire genome, allowing genetic divergence to accumulate virtually anywhere in the genome. In sympatric populations, however, the strong divergent selection required to overcome migration produces a genetic mosaic of divergent and non-divergent genomic regions. In some recent genome scans, each divergent genomic region has been interpreted as an independent incidence of migration/selection balance, such that the reduction of gene exchange is restricted to a few kilobases around each divergently selected gene. I propose an alternative mechanism, 'divergence hitchhiking' (DH), in which divergent selection can reduce gene exchange for several megabases around a gene under strong divergent selection. Not all genes/markers within a DH region are divergently selected, yet the entire region is protected to some degree from gene exchange, permitting genetic divergence from mechanisms other than divergent selection to accumulate secondarily. After contrasting DH and multilocus migration/selection balance (MM/SB), I outline a model in which genomic isolation at a given genomic location is jointly determined by DH and genome-wide effects of the progressive reduction in realized migration, then illustrate DH using data from several pairs of incipient species in the wild.

Heritability and genetic correlations of fear-related behaviour in Red Junglefowl--possible implications for early domestication

Domesticated species differ from their wild ancestors in a number of traits, generally referred to as the domesticated phenotype. Reduced fear of humans is assumed to have been an early prerequisite for the successful domestication of virtually all species. We hypothesized that fear of humans is linked to other domestication related traits. For three generations, we selected Red Junglefowl (ancestors of domestic chickens) solely on the reaction in a standardized Fear of Human-test. In this, the birds were exposed for a gradually approaching human, and their behaviour was continuously scored. This generated three groups of animals, high (H), low (L) and intermediate (I) fearful birds. The birds in each generation were additionally tested in a battery of behaviour tests, measuring aspects of fearfulness, exploration, and sociality. The results demonstrate that the variation in fear response of Red Junglefowl towards humans has a significant genetic component and is genetically correlated to behavioural responses in other contexts, of which some are associated with fearfulness and others with exploration. Hence, selection of Red Junglefowl on low fear for humans can be expected to lead to a correlated change of other behavioural traits over generations. It is therefore likely that domestication may have caused an initial suite of behavioural modifications, even without selection on anything besides tameness.