Domestication of the Dog from the Wolf Was Promoted by Enhanced Excitatory Synaptic Plasticity: A Hypothesis (original) (raw)

• Journal List
• Genome Biol Evol
• v.6(11); 2014 Nov
• PMC4255776

Genome Biol Evol. 2014 Nov; 6(11): 3115–3121.

Yan Li,1,2,† Guo-Dong Wang,1,2,† Ming-Shan Wang,1,2,3 David M. Irwin,1,4,5 Dong-Dong Wu,1,2,* and Ya-Ping Zhang1,2,*

Yan Li

1State Key Laboratory of Genetic Resources and Evolution, Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

2Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, China

Guo-Dong Wang

1State Key Laboratory of Genetic Resources and Evolution, Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

2Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, China

Ming-Shan Wang

1State Key Laboratory of Genetic Resources and Evolution, Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

2Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, China

3University of Chinese Academy of Sciences, Beijing, China

David M. Irwin

1State Key Laboratory of Genetic Resources and Evolution, Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

4Department of Laboratory Medicine and Pathobiology, University of Toronto, Ontario, Canada

5Banting and Best Diabetes Centre, University of Toronto, Ontario, Canada

Dong-Dong Wu

1State Key Laboratory of Genetic Resources and Evolution, Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

2Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, China

Ya-Ping Zhang

1State Key Laboratory of Genetic Resources and Evolution, Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

2Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, China

1State Key Laboratory of Genetic Resources and Evolution, Yunnan Laboratory of Molecular Biology of Domestic Animals, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, China

2Kunming College of Life Science, University of Chinese Academy of Sciences, Kunming, China

3University of Chinese Academy of Sciences, Beijing, China

4Department of Laboratory Medicine and Pathobiology, University of Toronto, Ontario, Canada

5Banting and Best Diabetes Centre, University of Toronto, Ontario, Canada

Corresponding author.

†These authors contributed equally to this work.

Associate editor: Greg Elgar

Copyright © The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact journals.permissions@oup.com

Supplementary Materials

Supplementary Data

GUID: CBC95463-E834-4AF8-96A7-C3E26CA11AEF

GUID: 8829C024-85A5-44FB-A97D-45D4D96F7315

Abstract

Dogs shared a much closer relationship with humans than any other domesticated animals, probably due to their unique social cognitive capabilities, which were hypothesized to be a by-product of selection for tameness toward humans. Here, we demonstrate that genes involved in glutamate metabolism, which account partially for fear response, indeed show the greatest population differentiation by whole-genome comparison of dogs and wolves. However, the changing direction of their expression supports a role in increasing excitatory synaptic plasticity in dogs rather than reducing fear response. Because synaptic plasticity are widely believed to be cellular correlates of learning and memory, this change may alter the learning and memory abilities of ancient scavenging wolves, weaken the fear reaction toward humans, and prompt the initial interspecific contact.

Keywords: gray wolf, self-domestication, fear response, learning, memory

Dogs have evolved unique social cognitive capabilities not found in their wolf progenitors (Hare et al. 2002; Miklósi et al. 2003; Topál et al. 2009). “Selection for communication” was proposed as the direct selective pressure that drove the evolution of these unusual abilities (Hare et al. 2002; Miklósi et al. 2003). Alternatively, the “correlated by-product” hypothesis proposed that these abilities was a by-product of selection for tameness toward humans, because tame foxes show greater skill in reading human gestures than control foxes (Hare et al. 2005), and hypothesized the reduced fearful-aggressive response, which largely shortened their distances from human presence, to be the prerequisite of dog domestication (Belyaev 1969). However, no genetic evidence has been reported that is directly associated with the precise aggressive-tame behavioral transformation, although several studies have identified genes that are involved in the neural system and are highly divergent from wolves (Saetre et al. 2004; Li et al. 2013; Wang et al. 2013).

Excess of Fixed Alleles within Stress-Related Genes in Dogs

We firstly compared published resequenced genomes of three wolves and ten dogs (including five ancient dogs and five modern dogs, supplementary material, Supplementary Material online) to identify the most significant genetic legacy in the dogs deviating from their progenitors. To avoid inaccurate estimation of population differentiation due to small sample size, we only count the single nucleotide polymorphisms (SNPs) that differentiate extremely between the wolves and the dogs (allele frequency is 1 in wolves but 0 in dogs, or vice versa), which were defined as fixed SNPs. We identified 204 genes that have at least six fixed SNPs (within the 95% percentile rank). These genes showed an extremely significant lower level of nucleotide diversity and Tajima’s D values (P = 5.22E-05 and 1.23E-30, respectively, Mann–Whitney U test) compared with other genes in the genome (fig. 1A), suggesting a potential selection effects on the divergence observed here. Because only a very small number of fixed SNPs (totally 26) were nonsynonymous substitutions, this may indicate that the positive selection operated mainly on expressional regulation. Actually, the 204 genes showed appreciable changes in expression patterns between dogs and wolves than others for two different measurements: Absolute expression change and fold change (P = 0.022 and P = 0.005, respectively) (fig. 1B), based on the transcriptome data for the frontal cortex (Albert et al. 2012). These results suggest that expressional variation rather than structural variation in protein sequence is the major contributor to the currently observed differentiation between dogs and wolves.

Analysis of selection in the dog genome. (A) Comparisons of the nucleotide diversity (left) and Tajima’s D values (right) between genes containing large numbers of fixed SNP differences and other genes ± S.D. were presented. (B) Comparisons of the difference in expression levels between wolves and dogs between genes containing large numbers of fixed SNP differences and other genes. The expression value for each gene was log2 transformed. Left: Expression difference of each gene between the wolf and the dog was calculated by the transformed value in the dog minus the transformed value in wolf. Right: Difference of each gene between the wolf and the dog was calculated by the transformed value in the dog divided by the transformed value in the wolf. (C) Left: Negative correlation between _F_ST values and recombination rates of genome wide SNPs. Right: Positive correlation between _F_ST values and recombination rates of SNPs at genes in GO categories: GO: 0001640 and GO: 0007216, both of which contain only one gene: GRIK3 in the Ensembl 72 dog annotation.

GO (gene ontology) analysis of the 204 genes revealed most overrepresentation in categories referring to “multicellular organismal response to stress” (P = 9.87E-4 adjusted by Benjamini–Hochberg FDR, False Discovery Rate), “behavioral fear response” (P = 1.41E-3) and “behavioral defense response” (P = 1.41E-3, table 1), thus supporting the hypothesis that positive selection caused a behavioral shift as dogs diverged from wolves. The first category, multicellular organismal response to stress, contained five genes: GRIK3, MECP2, BCL2, GRIK2, and GABRA5, whereas the other two categories each contained these same genes except GRIK3. All five of these genes are associated with the metabolism of glutamate (table 2), which is an important neurotransmitter in the brain (Purves et al. 2001). Because none of the fixed SNPs detected within these five genes were nonsynonymous, this suggests that shifted fear behavior that occurred during the initial domestication of the dog might be an outcome of a change in expression of the glutamate-related genes. In addition to the above genes, the gene HTR2C (5-hydroxytryptamine receptor 2C), which is involved in serotonin and dopamine pathway (Stam et al. 1994; Alex et al. 2005), has ten fixed SNPs differences between dogs and wolves, and also belongs to the behavioral fear response categories in the GO Annotation (www.geneontology.org). It shared interacting genes with its paralogue HTR2A, which has been suggested to modulate cognitive process by enhancing glutamate release (Feng et al. 2001).

Table 1

GO Analysis of Genes Containing Large Numbers of Fixed SNP Differences between Wolves and Dogs

Table 2

Description of the Functions of Five Genes Potentially under Selection in the Dog

Selective Signatures of the Dog in Glutamate-Related Signaling Pathway Genes

If selection for stress response was an initial target during domestication, then these fixed alleles should keep in a near fixed state even with amplified sampling. To test this, we resequenced the genomes of an additional three wolves and three Chinese native dogs presenting very rich genetic diversity (see supplementary material, Supplementary Material online, for more details). The fixed SNPs in the five genes: MECP2, BCL2, GRIK2, GABRA5, and GRIK3 identified above were present as a single allele or singleton in dogs. Furthermore, we calculated _F_ST for each SNP between dogs and wolves to evaluate the population differentiation, and identified GO categories for genes containing SNPs with _F_ST values statistically significantly higher than the average for SNPs for genome-wide genes. The GO categories showing the greatest statistical significance were “adenylate cyclase inhibiting G-protein coupled glutamate receptor activity” (GO: 0001640) and “G-protein coupled glutamate receptor signaling pathway” (GO: 0007216). Similarly, two pathways involved in glutamate receptor activity, “glutamate receptor signaling pathway” (GO: 0007215) and adenylate cyclase-inhibiting G-protein coupled glutamate receptor signaling pathway (GO: 0007196), were also among the top ten categories with greatest significances.

Because the _F_ST parameter does not show the direction of selection, and cannot identify upon which lineage, dog or wolf, explains the divergence for these categories, we applied the parameter XP-EHH (Sabeti et al. 2007) to evaluate selection on the SNPs in the dog lineage after divergence form the wolf. XP-EHH values on the dog lineage for genes involved in the glutamate receptor pathway retained statistically significant high values, suggesting that positive selection on these glutamate receptor pathway genes potentially occurred during the domestication of dog from wolf, and account for their change in behavior.

A considerable proportion of selective signatures was due to hitchhiking accompanied with the high intensity of artificial selection on a selected few genes. Accordingly, we next examined whether the observed signature of selection seen in glutamate metabolism genes was due to selection or hitchhiking. The Hill–Robertson effect states that selection is most effective when variants freely recombine (Hill and Robertson 1966). Selective sweeps are expected to extend less far in regions of higher recombination rate, and thus allele frequency differentiation is expected to be negatively correlated with recombination rate under hitchhiking. At the genomic level in dogs, the evolutionary rate for a SNP correlates negatively with its recombination rate (fig. 1C, r = −6.096e-04, P < 2e-16), which is consistent with the overall pattern observed in rice (Lu et al. 2006) and humans (Keinan and Reich 2010). In contrast, SNPs within GO: 0001640 category (the most divergent category) showed positive correlation between evolutionary rate and recombination rate, although the correlation was not statistically significant as it referred to only one gene (fig. 1C, r = 0.102, P = 0.325). Furthermore, GO:0007216 category from Ensembl version 74 (containing four genes: GRIK3, GRM5, GRM6, and TRPM1) showed a statistically significant positive correlation (r = 0.014, P = 0.00133). Thus, our result indicated that positive selection occurred on glutamate metabolism genes during the domestication of the dog.

Potential Function of Candidate Genes with Changed Expression Direction

Glutamate is the major excitatory neurotransmitter in the brain that regulates many kinds of behaviors and emotions and plays a key role in cognitive ability, including learning and memory through influencing short- and/or long-term potentiation (LTP) (Purves et al. 2001). Both GRIK2 (glutamate receptor, ionotropic, kainate 2) and GRIK3 (glutamate receptor, ionotropic, kainate 3) are glutamate receptors. GRIK2 knock-out mice exhibit significant reduction in anxiety and fear memory (Ko et al. 2005). Although no clear function has been identified for GRIK3, it coassembles with GRIK2 to form the kainate glutamate receptor (Dingledine et al. 1999), and deficits in mossy fiber LTP were observed in GRIK2 and GRIK3 knock-out animals (Contractor et al. 2001; Schmitz et al. 2003; Breustedt and Schmitz 2004; Pinheiro et al. 2007). Our analysis of the frontal cortex transcriptome data showed that GRIK2 is expressed at a significantly higher level in the frontal cortex of the dog than in the wolf (P = 0.0006 by the Mann–Whitney U test). Intriguingly, we also found a consistent up-regulation of GRIK2 in other domesticated animals compared with wild counterpart (student’s _t_-test), including chicken (P = 0.249), rat (P = 0.068), guinea pig (P = 0.045, data from Albert et al. [2012]), and rabbit (P = 0.381, data from Albert et al. [2012]) (supplementary table S1, Supplementary Material online), which showed a convergent evolution among domesticated animals. Increased transcription of GRIK2 should increase anxiety and fear memory (table 2). Consistent with the changes in GRIK2, BCL2, and GABRA5 also present changes (but no statistical significance) in their levels of expression in dogs compared with wolves that should increase the fear response in dogs (table 2).

We note that the changes in expression levels for these divergent genes were moderate, but they presented changes that contradict the expected expression pattern by the correlated by-product hypothesis, which proposed the fear reduction in the primary dogs to explain the prerequisite of the domestication. These moderate changes may be attributed to the minor effects of many genes underlying the selective targets, which may often occurred during the initial phase of domestication. Actually, according to the weighted gene coexpression network analysis (WGCNA) analysis (Langfelder and Horvath 2008), GRIK2, GRIK3, GABRA5, and MECP2 showed coexpression pattern and belonged to the same gene coregulatory network (e.g., module) which presented special positive correlation with the frontal cortex of wolf and dog (P = 4e-05 and 6e-05, respectively) (see supplementary fig. S1 and material, Supplementary Material online, for details). Moreover, GRIK2, GRIK3, and GABRA5 all present to be hub genes in this module (MM = 0.943, 0.914, and 0.917, respectively), indicating their important functions within this module on nervous system.

A Hypothesis of “Enhanced Excitatory Synaptic Plasticity”

It should be noted that the roles predicted for these genes in the fear response research (table 2) were all tested under Pavlovian fear conditioning, from which fear (the conditioned response) was trained to accompany a noxious stimuli. These Pavlovian tests contrast with both the fox experiment (Trut 1999) and dog domestication, where punishments were not received when the animals became close to humans. Additional pleiotropic functions of glutamate may have also contributed to the successful domestication of the wolf. The direction of change in the expression of the five genes should tend to cause excitatory synaptic plasticity in neural cells and/or benefit memory ability (although gene MECP2 locates in X chromosome, the equal sex ratio for both the domesticated and wild groups should eliminate the sex-linked effects on dosage). Consistent with this suggestion, dogs exhibit more excitatory behaviors than wolves, which sometimes becomes an overreaction yielding anxiety, or even obsessive-compulsive disorder, which may be associated with glutamate-related genes (Sampaio et al. 2013). Changes in synaptic plasticity are thought to be associated with changes in learning and memory abilities, by affecting short- and LTP (Purves et al. 2001). Thus, our results partially support the selection for communication hypothesis, where a strengthened learning ability should help the skill of reading human communicative behaviors. However, interspecific communication would only begin after a long period of scavenging life that enhanced the interactions between humans and wolves. In the “self-domestication” model, wolves domesticated themselves into dogs overtime of scavenging lifestyle (Coppinger and Coppinger 2001). In such a wild environment, the reduced fear response proposed by the correlated by-product hypothesis may be hard for these dog progenitors to survive. It therefore could be reasoned that during the early stages, the wolves with better learning and memory abilities would come close to human settlements more frequently, acquire greater food resources, and thus had greater opportunities to survive (with little disadvantage). These individuals would perform nonaggressive response because they would understand that the presence of humans was harmless, and thus would have a weakened fear reaction. We therefore propose a “selection for excitatory synaptic plasticity” hypothesis to account for the successful domestication of dogs from wolves. Following this hypothesis, affected learning and memory abilities would facilitate the behavioral shift, prolonged exposure to humans, and helped the dogs to understand the meaning of our gestures. Comparison of the genome of experimental foxes that have been tamed, and the unselected controls, may be an approach to test this hypothesis.

Materials and Methods

Reads of genome sequences were mapped onto the reference genome by using BWA-MEM (bio-bwa.sourceforge.net), and SNPs were calling by Genome Analysis Toolkit (McKenna et al. 2010) (GenomeAnalysisTK-2.6-4-g3e5ff60). The RNA-seq data from the frontal cortex of the wolf and the dog were from Albert et al. (2012). Tophat (Trapnell et al. 2009) and Cufflinks (Trapnell et al. 2010) were used to assemble transcripts and calculate the expression value of genes. GO analysis was performed using g:profiler (http://biit.cs.ut.ee/gprofiler/). Weighted gene coexpression networks were performed by WGCNA package implemented in R (Langfelder and Horvath 2008).

The full experimental methods are provided in supplementary material S1, Supplementary Material online.

Acknowledgments

This work was supported by grants from the National Natural Science Foundation of China (31321002;, 91231108;, 31123005), and the 973 program (2013CB835200).

Literature Cited

• Albert FW, et al. A comparison of brain gene expression levels in domesticated and wild animals. PLoS Genet. 2012;8:e1002962. [PMC free article] [PubMed] [Google Scholar]
• Alex KD, Yavanian GJ, McFarlane HG, Pluto CP, Pehek EA. Modulation of dopamine release by striatal 5-HT2C receptors. Synapse. 2005;55:242–251. [PubMed] [Google Scholar]
• Belyaev DK. Domestication of animals. Science. 1969;5:47–52. [Google Scholar]
• Breustedt J, Schmitz D. Assessing the role of GLUK5 and GLUK6 at hippocampal mossy fiber synapses. J Neurosci. 2004;24:10093–10098. [PMC free article] [PubMed] [Google Scholar]
• Clement Y, et al. Gabra5-gene haplotype block associated with behavioral properties of the full agonist benzodiazepine chlordiazepoxide. Behav Brain Res. 2012;233:474–482. [PubMed] [Google Scholar]
• Contractor A, Swanson G, Heinemann SF. Kainate receptors are involved in short- and long-term plasticity at mossy fiber synapses in the hippocampus. Neuron. 2001;29:209–216. [PubMed] [Google Scholar]
• Coppinger R, Coppinger L. Dogs: a new understanding of canine origin, behavior, and evolution. Chicago: University of Chicago Press; 2001. [Google Scholar]
• Dingledine R, Borges K, Bowie D, Traynelis SF. The glutamate receptor ion channels. Pharmacol Rev. 1999;51:7–61. [PubMed] [Google Scholar]
• Einat H, Yuan P, Manji HK. Increased anxiety-like behaviors and mitochondrial dysfunction in mice with targeted mutation of the Bcl-2 gene: further support for the involvement of mitochondrial function in anxiety disorders. Behav Brain Res. 2005;165:172–180. [PubMed] [Google Scholar]
• Feng J, Cai X, Zhao J, Yan Z. Serotonin receptors modulate GABA(A) receptor channels through activation of anchored protein kinase C in prefrontal cortical neurons. J Neurosci. 2001;21:6502–6511. [PMC free article] [PubMed] [Google Scholar]
• Hare B, Brown M, Williamson C, Tomasello M. The domestication of social cognition in dogs. Science. 2002;298:1634–1636. [PubMed] [Google Scholar]
• Hare B, et al. Social cognitive evolution in captive foxes is a correlated by-product of experimental domestication. Curr Biol. 2005;15:226–230. [PubMed] [Google Scholar]
• Hill WG, Robertson A. The effect of linkage on limits to artificial selection. Genet Res. 1966;8:269–294. [PubMed] [Google Scholar]
• Keinan A, Reich D. Human population differentiation is strongly correlated with local recombination rate. PLoS Genet. 2010;6:e1000886. [PMC free article] [PubMed] [Google Scholar]
• Ko S, Zhao MG, Toyoda H, Qiu CS, Zhuo M. Altered behavioral responses to noxious stimuli and fear in glutamate receptor 5 (GluR5)- or GluR6-deficient mice. J Neurosci. 2005;25:977–984. [PMC free article] [PubMed] [Google Scholar]
• Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559. [PMC free article] [PubMed] [Google Scholar]
• Li Y, et al. Artificial selection on brain-expressed genes during the domestication of dog. Mol Biol Evol. 2013;30:1867–1876. [PubMed] [Google Scholar]
• Lu J, et al. The accumulation of deleterious mutations in rice genomes: a hypothesis on the cost of domestication. Trends Genet. 2006;22:126–131. [PubMed] [Google Scholar]
• McKenna A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–1303. [PMC free article] [PubMed] [Google Scholar]
• Miklósi Á, et al. A simple reason for a big difference: wolves do not look back at humans, but dogs do. Curr Biol. 2003;13:763–766. [PubMed] [Google Scholar]
• Na ES, Nelson ED, Kavalali ET, Monteggia LM. The impact of MeCP2 loss- or gain-of-function on synaptic plasticity. Neuropsychopharmacology. 2013;38:212–219. [PMC free article] [PubMed] [Google Scholar]
• O'Driscoll CM, Kaufmann WE, Bressler JP. MeCP2 deficiency enhances glutamate release through NF-kappaB signaling in myeloid derived cells. J Neuroimmunol. 2013;265:61–67. [PubMed] [Google Scholar]
• Pinheiro PS, et al. GluR7 is an essential subunit of presynaptic kainate autoreceptors at hippocampal mossy fiber synapses. Proc Natl Acad Sci U S A. 2007;104:12181–12186. [PMC free article] [PubMed] [Google Scholar]
• Ponder CA, et al. Selection for contextual fear conditioning affects anxiety-like behaviors and gene expression. Genes Brain Behav. 2007;6:736–749. [PubMed] [Google Scholar]
• 2001. Purves D, et al. editors. 2001. Neuroscience. Sunderland (MA): Sinauer Associates.
• Rondi-Reig L, et al. Fear decrease in transgenic mice overexpressing bcl-2 in neurons. Neuroreport. 1997;8:2429–2432. [PubMed] [Google Scholar]
• Rondi-Reig L, Mariani J. To die or not to die, does it change the function? Behavior of transgenic mice reveals a role for developmental cell death. Brain Res Bull. 2002;57:85–91. [PubMed] [Google Scholar]
• Sabeti PC, et al. Genome-wide detection and characterization of positive selection in human populations. Nature. 2007;449:913–918. [PMC free article] [PubMed] [Google Scholar]
• Saetre P, et al. From wild wolf to domestic dog: gene expression changes in the brain. Mol Brain Res. 2004;126:198–206. [PubMed] [Google Scholar]
• Sampaio AS, et al. Genetic association studies in obsessive-compulsive disorder. Rev Psiq Clin. 2013;40:177–190. [Google Scholar]
• Schelman WR, Andres RD, Sipe KJ, Kang E, Weyhenmeyer JA. Glutamate mediates cell death and increases the Bax to Bcl-2 ratio in a differentiated neuronal cell line. Brain Res Mol Brain Res. 2004;128:160–169. [PubMed] [Google Scholar]
• Schmitz D, Mellor J, Breustedt J, Nicoll RA. Presynaptic kainate receptors impart an associative property to hippocampal mossy fiber long-term potentiation. Nat Neurosci. 2003;6:1058–1063. [PubMed] [Google Scholar]
• Stam NJ, et al. Genomic organisation and functional expression of the gene encoding the human serotonin 5-HT2C receptor. Eur J Pharmacol. 1994;269:339–348. [PubMed] [Google Scholar]
• Topál J, Gergely G, Erd hegyi Á, Csibra G, Miklósi Á. Differential sensitivity to human communication in dogs, wolves, and human infants. Science. 2009;325:1269–1272. [PubMed] [Google Scholar]
• Trapnell C, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28:511–515. [PMC free article] [PubMed] [Google Scholar]
• Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105–1111. [PMC free article] [PubMed] [Google Scholar]
• Trut LN. Early canid domestication: the farm-fox experiment. Am Sci. 1999;87:160–169. [Google Scholar]
• Wang DS, et al. Memory deficits induced by inflammation are regulated by alpha5-subunit-containing GABAA receptors. Cell Rep. 2012;2:488–496. [PMC free article] [PubMed] [Google Scholar]
• Wang GD, et al. The genomics of selection in dogs and the parallel evolution between dogs and humans. Nat Commun. 2013;4:1860. [PubMed] [Google Scholar]
• Wei M, Caballero A, Hill WG. Selection response in finite populations. Genetics. 1996;144:1961–1974. [PMC free article] [PubMed] [Google Scholar]

Articles from Genome Biology and Evolution are provided here courtesy of Oxford University Press