Strong selection for behavioural resilience in Australian stock working dogs identified by selective sweep analysis
© Arnott et al.; licensee BioMed Central. 2015
Received: 6 January 2015
Accepted: 26 March 2015
Published: 7 May 2015
Working dog handlers and breeders have strong opinions on characteristics that are desirable in the breeds that they use to handle stock. Most of these characteristics are related to conformation or behaviour. This study explored whether the genetics underlying desirable working behaviour traits might be identified by selective sweep analysis; a method that identifies long regions of strong homozygosity combined with allelic divergence from a comparison group. For this analysis, we compared genomic haplotype architecture in two breeds derived from common founder stock but subjected to divergent selective pressures. The breeds studied were the Australian Kelpie, which is registered with the Australian National Kennel Council, and the Australian Working Kelpie, which is registered with the Working Kelpie Council.
A selective sweep spanning 3 megabases on chromosome 3 was identified in the Australian Working Kelpie. This region is the location of genes related to fear-memory formation and pain perception.
Selective sweep loci of similar magnitude were observed in the Australian Kelpie. On chromosome 8 is a locus which may be related to behavioural excitability and on chromosome 30 is a smaller locus which most likely is related to morphology.
Active working stock dogs of the Australian Working Kelpie breed have been bred primarily for gene loci influencing pain perception and fear memory formation. By contrast Australian Kelpies are commonly maintained in urban environments where these characteristics are not required and have been affected by selection for conformation and coat colour. The identified loci may aid in the identification of superior working dogs.
The term ‘selective sweep’ is used in the study of genetics to describe a reduction or loss of DNA sequence variation in regions of the genome of species, breeds or cultivars. This can occur as a consequence of strong selective pressure (positive evolutionary selection) due to a highly desirable DNA mutation in such regions which convey a survival advantage (positive evolutionary natural selection). Alternatively it can be artificially driven by man in other species through manipulating selection with eugenic breeding programmes. A selective sweep analysis can be very helpful in both determining the close genetic relationship of individuals and groups within a species, and also in some cases can be useful in identifying gene variants that cause disease.
In this study of Australian Kelpies we have applied a fresh approach to selective sweep analysis that interrogates a breed ‘split’ so as to learn more about the ‘external characteristics’ (often called the phenotype) that are regarded as desirable in two very different cohorts of dogs that share a common breed origin. One group of dogs is intensively selected for its ability to work with livestock, while the other group is bred for conformation and companionship – usually in an urban setting. This study suggests that “working Kelpies” have a strong genomic selection signal in a region of their genome that contains genes predominantly concerned with resilience traits such as the ability to sense pain and form memories associated with fear. The “nonworking” Kelpie group have selection signals that are most likely related to body shape and size (conformation traits).
A selective sweep is the reduction or elimination of variation among the nucleotides in neighboring DNA of a mutation as the result of recent and strong positive natural or artificial selection.
Selection, both natural and artificial, can be very powerful in shaping phenotype and, as widely described in the literature, this is broadly demonstrated in the morphologies and behaviours of domestic dogs. Previous studies have assessed selective sweep by comparison of disparate dog breeds [1-6]. In this work, we concentrate on comparing two closely related breeds that have been derived from a common ancestor but then bred for different purposes. By doing so, we expect to be able to more easily identify regions that underpin the observable behavioural and physical differences between the breeds.
Where selection is based upon a common breeding goal, selective progress in breed improvement is enhanced; a strategy widely used in the production animal industries. Selection in dog populations is assisted by the birth of progeny in litters that may number as high as 15 individuals but are more commonly in the range of 5–7 offspring. The prolificacy enables breeders to exert considerable selection pressure within a family, albeit with a demonstrably low correlation between pup behavioural assessment and adult working success . A relatively large amount of selection pressure in dog breeding occurs at a young age (approximately 8 weeks) when the breeder elects either to keep a pup potentially for breeding or to sell it. A second stage of selection occurs when the adult is identified as being suitable for breeding.
Methodologies for the assessment of selective sweep have been explored extensively in the literature [6,14,15]. For the current study, we elected to use a simple analysis of extended homozygosity over one megabase genomic windows combined with a long region of heterogeneous allele frequency (FST)  relative to the alternate population with the goal of identifying genomic regions that have responded best to the strong selective forces applied by breeders.
Results and discussion
Regions identified in the primary analysis for differential selective sweep including 12 Australian Kelpie and 12 Australian Working Kelpie dogs that were then assessed for validation in a further 10 AK and 28 AWK individuals
Position (1 Mb bin) canfam 3.0
Markers in window
maf £ AWK †
maf AK ‡
maf absolute difference
maf AWK †
maf AK ‡
maf absolute difference
Two loci were identified for the Australian Kelpie. The first is on chromosome 8 and spans 4 Mb from 64 to 68 Mb with the strongest signal at 67 Mb. The second is on chromosome 30 and is supported in the wider data for a single window at 23 Mb. Adjoining windows spanning an area of 21 Mb to 25 Mb maintained mean minor allele frequencies of lower than 0.07 in the larger data cohort.
Positional candidate genes for behaviour and morphology phenotypes according to the mouse genome browser
Positional candidate genes
Canine Locus (Canfam 3.0)
Mouse Locus (GRCm38/mm10)
Homer1, Arsb, Lhfpl2, Ap3b1,Crhbp, F2rl1
Msh3, Homer1, Ap3b1,F2rl1
Dicer1, Bdkrb2,Bdkrb1, Bcl11b
Dicer1, Bdkrb2,Bdkrb1, Bcl11b
Dicer1, Ak7, Bcl11b
Dicer1, Bdkrb2, Bcl11b
Both breeds appeared to have a single predominant common haplotype on the X chromosome (Additional file 1: Table S1). This is not surprising as both groups in our data were founded on a single female that had offered exceptional working ability . One additional locus of sweep for the AK was excluded from the reported results due to low density of markers in the window (at chr16:59000000) but is nonetheless observably supported in the dog phenotypes. The Australian Kelpie is frequently self-coloured brown or black (although the brown may be reported as either red or chocolate by breeders, depending on the hue). In contrast, Australian Working Kelpies may be the same colours as those reported for the AK, but most frequently exhibit tan markings (commonly referred to as black and tan or red and tan). One window of selective sweep for the AK (chr16:59000000) sits directly on the K-locus (DEFB103)  that affects the expression of pheomelanin versus eumelanin in the dog. No sweep was observed in either breed for the agouti locus that has been previously described as the major locus underlying the black and tan coat colour phenotype [20,21].
It is not the first time that the locus on chromosome 3 for the AWK has been identified in research on this breed. In an analysis to identify a putatively recessive locus underpinning an inherited disorder in Kelpies, Shearman et al.  identified this exact same region as a region of extended homozygosity in samples affected by cerebellar ataxia. We predict that the earlier disease association study had included a mixture of AK and AWK in the control cohort and that the signal identified is based on cohort stratification. While some dogs in the current population participated in the study of cerebellar ataxia, our analysis excluded all case samples from the data and doing so clearly reveals that homozygosity in this region is a characteristic of all AWK.
Presuming that the region in this sweep is affected by strong selection for working success, the most interesting of the positional candidate genes on chromosome 3 is the HOMER1 gene. This gene is associated with fear memory formation and pain perception (nociception) in the mouse [23-25]. The areas where working stock dogs are employed in Australia have large numbers of environmental hazards. In particular, a large proportion of the groundcovers are of a spiked nature and species such as Tribulus terrestris (cat-head burr), Xanthium spp (Bathurst burr, Noogoora burr), Onopordum acanthium (Scotch thistle), Alternanthera pungens (Khaki weed), Nassella neesiana (Chilean needle grass) and Cenchrus spp (Spiny burrgrass) are common (http://www.weeds.org.au/). In addition, traumatic injuries caused by livestock, fences and vehicles are well documented . A dog that can overcome pain to maintain sufficient focus and continue its work will be a strong asset to the working dog handler. Dogs that are resilient will also have an enhanced chance of working success in the challenging environment of Australian stockwork. Interestingly, an important founder sire for the AWK, a blue dog named “Coil”, is renowned for his exceptional pain tolerance and endurance. Coil won the 1898 Sydney trial achieving a perfect score despite competing with a fractured foreleg . Further research is required to explore whether pain thresholds truly differ between these cohorts of dogs.
For the Australian Kelpie, the major identified sweep locus on chromosome 8 contains genes that relate to both behaviour and morphology but none that stands out as an obvious candidate. The region contains a large gene poor section that has an enriched selection of highly conserved non-coding elements. The gene BCL11B (B-Cell Lymphoma/Leukemia 11B) is a transcriptional repressor and one gene, UNC79 (Protein Unc-79 Homolog) just outside of the swept interval but in the vicinity has been linked with hyperactivity in the mouse .
Given the activity of the genes in this region, the driver for the chromosome 30 sweep locus in the AK likely has a morphological, rather than behavioural, basis.
By focusing on breeds that are derived from a common foundation but then selected for different purposes we have been able to identify a major locus underlying ability of stock dogs, represented by the Australian Working Kelpie, to effectively work in harsh environmental conditions. While stock dog breeders may be selecting primarily for traits such as stock sense and boldness, we reveal that they are favouring breeding from dogs that can focus and continue working despite the presence of environmental hazards and discomforts. This requirement is not needed for dogs accustomed to an urban lifestyle. Australian Kelpies are not usually employed in stock work. They appear to have been subjected to selection that is predominantly based on morphology rather than behaviour.
Twelve Australian Working Kelpie (AWK) dogs (registered with the Working Kelpie Council) and twelve Australian Kelpie (AK) dogs (registered with the Australian National Kennel Council) were used in the primary analysis. Ten additional AK (representing data from an unrelated family collected for an unrelated disease study) and 28 additional AWK (collected for unrelated studies) were used to validate the loci identified in the primary analysis. Case samples from the other studies were excluded.
Peripheral blood samples were obtained from six dogs using EDTA blood collection tubes (BD Vacutainer, BD Franklin Lakes, NJ) and extracted using EZ1® DNA Blood Kits (Qiagen, Valencia, CA). The remaining two samples were collected using Oragene ANIMAL OA-400 saliva collection kits (DNA Genotek, Ontario Canada) and extracted following standard kit-issued protocol.
Both primary analysis and validation samples were collected with University of Sydney ethics clearance (N00/10-2012/3/5837 and N00/10-2012/3/5928).
Genotyping was conducted on the Illumina Canine High Density Genotyping array (170,000 markers) by Neogen/Geneseek Nebraska USA.
Results were assessed for minor allele frequency within dog registry cohorts using the --freq option in the Plink software package . Mean minor allele frequencies were binned in one megabase non-overlapping sliding windows across the genome within cohort.
For analysis of the Australian Working Kelpie selective sweep(s), one megabase windows with more than 40 single nucleotide polymorphism (SNP) observations and a mean minor allele frequency of lower than 0.05 were identified and compared with the same window for the Australian Kelpie. Where the absolute difference between minor allele frequencies in the compared groups was greater than the arbitrarily set value of 0.11 the region was selected for validation. Next, windows with a mean minor allele frequency of lower than 0.05 in the Australian Kelpie and an absolute difference in minor allele frequency (>0.11) relative to the Australian Working Kelpie were identified for validation as indicating sweep in the Australian Kelpie.
Finally, both data sets were expanded to include the additional 38 dogs (10 AK and 28 AWK) and the analysis was repeated.
The intervals identified as representing selective sweeps were compared with the conserved syntenic regions on the mouse genome using the mouse genome browser hosted by Mouse Genome Informatics . Genes with relevance for behaviour, integument, craniofacial and skeletal phenotypes were identified using the phenotype utility within the browser.
This work was funded by grants from the Rural Industries Research and Development Corporation, Meat and Livestock Australia, the Working Kelpie Council and the Manning Bequest administered by the University of Sydney Faculty of Veterinary Science.
- Vonholdt BM, Pollinger JP, Lohmueller KE, Han E, Parker HG, Quignon P, et al. Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication. Nature. 2010;464(7290):898–902. doi:nature08837 [pii]10.1038/nature08837.View ArticlePubMed CentralPubMedGoogle Scholar
- Vaysse A, Ratnakumar A, Derrien T, Axelsson E, Rosengren Pielberg G, Sigurdsson S, et al. Identification of genomic regions associated with phenotypic variation between dog breeds using selection mapping. PLoS Genet. 2011;7(10):e1002316. doi:10.1371/journal.pgen.1002316PGENETICS-D-11-00264 [pii].View ArticlePubMed CentralPubMedGoogle Scholar
- Boyko AR, Quignon P, Li L, Schoenebeck JJ, Degenhardt JD, Lohmueller KE, et al. A simple genetic architecture underlies morphological variation in dogs. PLoS Biol. 2010;8(8):e1000451. doi:10.1371/journal.pbio.1000451.View ArticlePubMed CentralPubMedGoogle Scholar
- Quilez J, Short AD, Martinez V, Kennedy LJ, Ollier W, Sanchez A, et al. A selective sweep of >8 Mb on chromosome 26 in the Boxer genome. BMC Genomics. 2011;12:339. doi:10.1186/1471-2164-12-339.View ArticlePubMed CentralPubMedGoogle Scholar
- Bannasch D, Young A, Myers J, Truve K, Dickinson P, Gregg J, et al. Localization of canine brachycephaly using an across breed mapping approach. PLoS One. 2010;5(3):e9632. doi:10.1371/journal.pone.0009632.View ArticlePubMed CentralPubMedGoogle Scholar
- Pollinger JP, Bustamante CD, Fledel-Alon A, Schmutz S, Gray MM, Wayne RK. Selective sweep mapping of genes with large phenotypic effects. Genome Res. 2005;15(12):1809–19. doi:10.1101/gr.4374505.View ArticlePubMed CentralPubMedGoogle Scholar
- Goddard ME, Beilharz RG. Genetics of traits which determine the suitability of dogs as guide-dogs for the blind. Appl Anim Ethol. 1983;9(3–4):299–315.View ArticleGoogle Scholar
- The_Working_Kelpie_Council_of_Australia. Origin of the Kelpie. 2014. http://www.wkc.org.au/About-Kelpies/Origin-of-the-Kelpie.php. Accessed 30/09/2014.
- Arnott ER, Early JB, Wade CM, McGreevy PD. Environmental Factors Associated with Success Rates of Australian Stock Herding Dogs. PloS One. 2014;9(8). doi:ARTN e104457DOI 10.1371/journal.pone.0104457.Google Scholar
- Arnott ER, Early JB, Wade CM, McGreevy PD. Estimating the economic value of Australian stock herding dogs. Anim Welf. 2014;23(2):189–97. doi:10.7120/096272220.127.116.11.View ArticleGoogle Scholar
- Early JB, Arnott ER, Wade CM, McGreevy PD. Manual muster: a critical analysis of the use of common terms in Australian working dog manuals. J Vet Behav. 2014;9(6):370–4.View ArticleGoogle Scholar
- van Rooy D, Arnott ER, Early JB, McGreevy PD, Wade CM. Holding back the genes: limitations of research into canine behavioural genetics. Canine Genet Epidemiol. 2014;1:1–7. doi:10.1186/2052-6687-1-7.View ArticleGoogle Scholar
- Australian_National_Kennel_Council. Australian Kelpie extended breed standard. 2015. http://ankc.org.au/media/pdf/635576344320930744_d0d9014f-a85b-407a-8dc0-d2010a9293e0.pdf. Accessed 18/03/2015.Google Scholar
- Sabeti PC, Reich DE, Higgins JM, Levine HZ, Richter DJ, Schaffner SF, et al. Detecting recent positive selection in the human genome from haplotype structure. Nature. 2002;419(6909):832–7. doi:10.1038/nature01140.View ArticlePubMedGoogle Scholar
- de Simoni Gouveia JJ, da Silva MV, Paiva SR, de Oliveira SM. Identification of selection signatures in livestock species. Genet Mol Biol. 2014;37(2):330–42.View ArticlePubMed CentralPubMedGoogle Scholar
- Felsenstein J. How Can We infer geography and history from gene-frequencies. J Theor Biol. 1982;96(1):9–20. doi:10.1016/0022-5193(82)90152-7.View ArticlePubMedGoogle Scholar
- Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005;438(7069):803–19. doi:10.1038/nature04338.View ArticlePubMedGoogle Scholar
- Hoeppner MP, Lundquist A, Pirun M, Meadows JRS, Zamani N, Johnson J, et al. An Improved Canine Genome and a Comprehensive Catalogue of Coding Genes and Non-Coding Transcripts. PloS One. 2014;9(3). doi:ARTN e91172DOI10.1371/journal.pone.0091172.Google Scholar
- Candille SI, Kaelin CB, Cattanach BM, Yu B, Thompson DA, Nix MA, et al. A -defensin mutation causes black coat color in domestic dogs. Science. 2007;318(5855):1418–23. doi:10.1126/science.1147880.View ArticlePubMed CentralPubMedGoogle Scholar
- Dreger DL, Parker HG, Ostrander EA, Schmutz SM. Identification of a mutation that is associated with the saddle Tan and black-and-Tan phenotypes in basset hounds and Pembroke welsh corgis. J Hered. 2013;104(3):399–406. doi:10.1093/jhered/est012.View ArticlePubMedGoogle Scholar
- Dreger DL, Schmutz SM. A SINE insertion causes the black-and-Tan and saddle Tan phenotypes in domestic dogs. J Hered. 2011;102:S11–8. doi:10.1093/jhered/esr042.View ArticlePubMedGoogle Scholar
- Shearman JR, Cook RW, McCowan C, Fletcher JL, Taylor RM, Wilton AN. Mapping cerebellar abiotrophy in Australian Kelpies. Anim Genet. 2011;42(6):675–8. doi:10.1111/j.1365-2052.2011.02199.x.View ArticlePubMedGoogle Scholar
- Inoue N, Nakao H, Migishima R, Hino T, Matsui M, Hayashi F, et al. Requirement of the immediate early gene vesl-1S/homer-1a for fear memory formation. Mol Brain. 2009;2:7. doi:10.1186/1756-6606-2-7.View ArticlePubMed CentralPubMedGoogle Scholar
- Klugmann M, Szumlinski KK. Targeting Homer genes using adeno-associated viral vector: lessons learned from behavioural and neurochemical studies. Behav Pharmacol. 2008;19(5–6):485–500. doi:10.1097/FBP.0b013e32830c369f.View ArticlePubMedGoogle Scholar
- Szumlinski KK, Lominac KD, Kleschen MJ, Oleson EB, Dehoff MH, Schwarz MK, et al. Behavioral and neurochemical phenotyping of Homer1 mutant mice: possible relevance to schizophrenia. Genes Brain Behav. 2005;4(5):273–88. doi:10.1111/j.1601-183X.2005.00120.x.View ArticlePubMedGoogle Scholar
- Singh I, Tucker LA, Gendall P, Rutherfurd-Markwick KJ, Cline J, Thomas DG. Age, breed, sex distribution and nutrition of a population of working farm dogs in New Zealand: results of a cross-sectional study of members of the New Zealand sheep dog trial association. N Z Vet J. 2011;59(3):133–8. doi:10.1080/00480169.2011.567967.View ArticlePubMedGoogle Scholar
- Lu B, Zhang Q, Wang H, Wang Y, Nakayama M, Ren D. Extracellular calcium controls background current and neuronal excitability via an UNC79-UNC80-NALCN cation channel complex. Neuron. 2010;68(3):488–99. doi:10.1016/j.neuron.2010.09.014.View ArticlePubMed CentralPubMedGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75. doi:10.1086/519795.View ArticlePubMed CentralPubMedGoogle Scholar
- Blake JA, Bult CJ, Eppig JT, Kadin JA, Richardson JE. The Mouse Genome Database: integration of and access to knowledge about the laboratory mouse. Nucleic Acids Res. 2014;42(Database issue):D810–7. doi:10.1093/nar/gkt1225.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.