Genomic mapping of social behavior traits in a F2 cross derived from mice selectively bred for high aggression
© Nehrenberg et al; licensee BioMed Central Ltd. 2010
Received: 23 July 2010
Accepted: 31 December 2010
Published: 31 December 2010
Rapid response to selection was previously observed in mice selected for high levels of inter-male aggression based on number of attacks displayed in a novel social interaction test after isolation housing. Attack levels in this high aggression line (NC900) increased significantly within just four generations of selective breeding, suggesting the presence of a locus with large effect. We conducted an experiment using a small (n ≈ 100) F2 cross between the ICR-derived, non-inbred NC900 strain and the low aggression inbred strain C57BL/6J, genotyped for 154 fully informative SNPs, to determine if a locus with large effect controls the high-aggression selection trait. A second goal was to use high density SNP genotyping (n = 549,000) in the parental strains to characterize residual patterns of heterozygosity within NC900, and evaluate regions that are identical by descent (IBD) between NC900 and C57BL/6J, to determine what impacts these may have on accuracy and resolution of quantitative trait locus (QTL) mapping in the F2 cross.
No evidence for a locus with major effect on aggressive behavior in mice was identified. However, several QTL with genomewide significance were mapped for aggression on chromosomes 7 and 19 and other social behavior traits on chromosomes 4, 7, 14, and 19. High density genotyping revealed that 28% of the genome is still segregating among the six NC900 females used to originate the F2 cross, and that segregating regions are present on every chromosome but are of widely different sizes. Regions of IBD between NC900 and C57BL/6J are found on every chromosome but are most prominent on chromosomes 10, 16 and X. No significant differences were found for amounts of heterozygosity or prevalence of IBD in QTL regions relative to global analysis.
While no major gene was identified to explain the rapid selection response in the NC900 line, transgressive variation (i.e. where the allele from the C57BL/6J increased attack levels) and a significant role for dominant gene action were hallmarks of the genetic architecture for aggressive behavior uncovered in this study. The high levels of heterozygosity and the distribution of minor allele frequency observed in the NC900 population suggest that maintenance of heterozygosity may have been under selection in this line.
Fighting is a near universal survival trait expressed in animal species as varied as flies , mice  and humans. Its ubiquity suggests that it serves similar functions. Within a species, fighting may function to disperse its members in ways that reduce pressure on resources necessary for the species to survive . At an individual level, fighting is a strategy for winning competitions for territorial resources necessary for individuals and their relatives to survive [1, 2]. For example, the house mouse (Mus musculus) regularly patrols the borders of its territory and the highest levels of fighting are observed in areas containing vital resources .
A genetic basis for mouse aggression is clearly supported by the success of several different types of selective breeding programs. Divergent selection for attack latency in wild mice  and for aggression in Swiss albino  and Institute of Cancer Research (ICR; ) mice all produced significant differences in aggression levels within 5 generations of selective breeding. Because aggression selection effects occurred so rapidly across several different selective breeding programs, it is reasonable to postulate that the aggression selection response at least partly involves a genetic locus of major effect. Several selection studies have uncovered major genes, including the mini-muscle locus in mice selected for high-levels of voluntary wheel running (HR; ) and the high-growth locus . Several mouse aggression QTL have been reported [10, 11], and other available evidence to date argues against single locus control of aggression [12, 13]. However, it remains possible that QTL with large effects could account for the rapid selection effects of aggression in outbred lines.
The availability of a mouse line (NC900) selectively bred for high levels of attack behavior, created in a high and low aggression selection breeding program from an ICR base population , offers a unique opportunity to begin to dissect the genetic architecture of mouse aggression, because in the NC900 lines the number of attacks was the sole selection criterion, whereas in other mouse aggression selection programs attack latencies  or rated scores  were used. The NC900 selection criterion was attack counts displayed towards a group-reared unselected ICR mouse in a 10-min novel social interaction test following isolation housing at weaning . Levels of aggression in NC900 rapidly diverged from a contemporary low-aggression selection line within just four generations of selective breeding, and these high levels were maintained throughout the long-term selective breeding program . Therefore, we sought to determine whether a locus with large effects controls NC900 male aggression by performing QTL analyses of NC900 social phenotypes.
A straightforward index of attack behavior is to measure its frequency, duration, and latency, and these measures are highly correlated . But in most populations, there will be some mixture of attackers and non-attackers producing both continuous and categorical variables of attack and non-attack. Therefore, we designed a coding system based on the view that there is one continuum of social behavior  having appetitive (affiliative in the case of non-agonistic social behaviors) and aversive spectrums. The existence of an aversive, defensive spectrum of behavior is partly supported by evidence that intensity of mouse defensive responses increases in proportion to the intensity of threat [16, 17]. However, we recognize that it is not necessarily the case that affiliative and aversive social responses are inversely associated  and that agonistic mice are capable of expressing social affiliation at levels comparable to non-agonistic mice . The expression of agonism and non-agonistic sociability are likely a net result complex and perhaps conflicting motivations , stemming from somewhat overlapping and distinct biological mechanisms. Thus, assessing indices of affiliation and aversion taken from social interactions may improve our understanding underlying biological differences between purely agonistic, nonagonistic and mixed agonistic/nonagonistic phenotypes.
With this new social interaction coding system as well as the standard measures of attack, our experiments were designed to address questions based on results from a whole genome scan in a small F2 population created by crossing NC900 to the non-aggressive C57BL/6J (B6) inbred line. First, is NC900 aggression controlled by a locus with large effect? Second, can we detect QTL for social behavior measures other than the classic measures of aggression such as frequency and latency of fighting? Finally, we used high-density SNP genotyping in the parental strains to characterize residual patterns of heterozygosity within NC900, and IBD between NC900 and C57BL/6J, to evaluate their potential impacts on QTL mapping in the F2 cross. While no major gene was identified to explain the rapid selection response in the NC900 line, transgressive variation and a significant role for dominant gene action were hallmarks of the genetic architecture for aggressive behavior uncovered in this study. We found high levels of residual heterozygosity in the NC900 line, potentially posing limitations in QTL mapping using crosses of non-inbred lines, and suggesting that maintenance of heterozygosity may have been under selection for the high aggression phenotype.
Social behavior phenotype means by non-aggressive and aggressive subpopulations.
Social behavior phenotype
Non-aggression mean (± s.e.)
Aggression mean (± s.e.)
4.43 (± 0.39)
5.58 (± 0.48)
2.50 (± 0.49)***
3.01 (± 0.31)
3.97 (± 0.42)
1.75 (± 0.33)***
10.55 (± 0.74)
12.88 (± 0.91)
6.32 (± 0.96)***
13.86 (± 1.44)
17.27 (± 2.10)
9.38 (± 1.65)**
5.47 (± 0.53)
7.80 (± 0.73)
2.39 (± 0.41) ***
1.83 (± 0.21)
2.59 (± 0.31)
0.84 (± 0.14)***
19.93 (± 1.37)
26.56 (± 1.51)
11.21 (± 1.64)***
18.71 (± 1.59)
23.83 (± 2.14)
11.97 (± 1.91)***
3.24 (± 0.42)
4.72 (± 0.59)
1.29 (± 0.23)***
3.96 (± 0.64)
5.46 (± 0.81)
1.99 (± 0.94)**
1.97 (± 0.25)
2.24 (± 0.35)
1.61 (± 0.36)
1.12 (± 0.19)
1.28 (± 0.27)
0.91 (± 0.26)
0.32 (± 0.09)
0.44 (± 0.15)
0.16 (± 0.09)
0.32 (± 0.12)
0.46 (± 0.20)
0.13 (± 0.08)
5.42 (± 0.82)
0.00 (± 0.00)
12.55 (± 1.12)na
10.73 (± 1.56)
0.00 (± 0.00)
24.85 (± 1.96)na
10.94 (± 0.82)
7.40 (± 0.87)
15.61 (± 1.13)***
16.13 (± 1.51)
7.19 (± 0.99)
27.88 (± 2.02)***
30.87 (± 1.20)
26.82 (± 1.55)**
34.84 (± 1.53)
31.02 (± 1.95)
39.86 (± 2.23)**
17.07 (± 2.83)
0.00 (± 0.00)
39.53 (± 4.41)na
368.18 (± 26.10)
600.00 (± 0.00)
126.58 (± 18.98)na
Phenotypic Pearson partial correlations and p-values (with Bonferroni adjustments) for measured traits.
bidirectional count (1)
% bidirectional duration (2)
unidirectional count (3)
% unidirectional duration (4)
approach count (5)
% approach duration (6)
affiliative count (7)
% affiliative duration (8)
passive avoidance count (9)
% passive avoidance duration (10)
active avoidance count (11)
bidirectional count (1)
% bidirectional duration (2)
unidirectional count (3)
% unidirectional duration (4)
approach count (5)
approach duration (6)
affiliative count (7)
affiliative duration (8)
passive avoidance count (9)
passive avoidance duration (10)
active avoidance count (11)
active avoidance duration (12)
freezing count (13)
freezing duration (14)
aggression count (15)
aggression duration (16)
aversion count (17)
aversion count (18)
total social count (19)
% total social duration (20)
attack latency (22)
F2 social behavior QTL
QTL results of genome-wide scans for social interaction QTL in NC900 × B6 F2 mice.
36.95 - 59.25
-34.09 (± 43.15)
189.30 (± 53.13)
7.76 - 47.36
-12.86 (± 3.51)
-12.10 (± 5.29)
20.56 - 48.86
-3.58 (± 1.02)
-2.93 (± 1.56)
% total social duration
21.38 - 56.38
2.06 (± 1.90)
-10.00 (± 2.90)
% affilliative duration
37.80 - 56.60
3.81 (± 2.00)
-9.19 (± 3.01)
40.55 - 59.64
-0.06 (± 0.80)
3.74 (± 1.04)
% approach duration
37.80 - 56.60
0.62 (± 0.31)
1.21 (± 0.44)
0 - 36.18
4.31 (± 1.57)
-8.97 (± 2.64)
% bidirectional duration
0 - 15.26
0.87 (± 0.35)
-2.03 (± 0.59)
0 - 18.95
0.53 (± 0.43)
-2.52 (± 0.74)
% aggression duration
20.46 - 47.46
-7.15 (± 1.55)
-6.40 (± 2.99)
A significant QTL for number of attacks was detected on MMU19 with a peak location of 34.5 cM. The attack QTL exhibited both dominance and additive effects, with the B6 allele contributing to an increase in phenotypic value, and underdominance (Figure 3I). Two other QTL, for aggression count and percent aggression duration, also mapped to MMU19 very near the peak for attack (Figure 3J and 3K). Like the QTL for attack, the B6 allele contributed to the increase in the aggression phenotype value, and underdominance was present. A suggestive QTL for attack latency was detected on MMU7 displaying underdominant gene action (Figure 3C).
Two QTL with underdominant gene action were also detected on MMU4: one suggestive QTL for percent total social duration and one significant QTL for percent affiliative duration (Figure 3A and 3B). An overdominant QTL for approach count was detected on MMU7 (Figure 3D) with a QTL peak near the QTL peak for attack latency. A QTL with both additive and dominance effects for percent approach duration was detected on MMU14 (Figure 3E). Three QTL having a peak location at 7.76 cM were detected on MMU19: a suggestive QTL for affiliative count, a suggestive QTL for percent bidirectional duration and a significant QTL for bidirectional count, and all three exhibited underdominant effects (Figure 3F, G, and 3H) with the NC900 allele contributing to the increase in value for the additive component.
In summary, four QTL regions were detected controlling 12 social behavior traits. MMU4 controlled the percent duration of total social behaviors and affiliation. MMU7 controlled both attack latency and approach count. MMU14 controlled percent approach duration. Finally, MMU19 controlled a total of six social behavior traits, including affiliative count, bidirectional count and percent bidirectional duration, attack and aggression count, and percent aggression duration.
Effects of coat color on social behavior
Effects of coat color on social behavior means ± s.e in parentheses.
246 (± 55)a
3.11 (± 1.00) a
8.74 (± 2.00) a
395 (± 43)b
5.94 (± 0.86) a
5.16 (± 1.44) a
407 (± 39)b
6.65 (± 6.65) a
3.95 (± 1.02) a
Selection lines are maintained as outbred populations and thus typically segregate for significant but variable regions of the genome. Residual heterozygosity in the NC900 selection line and the amount and distribution of regions of IBD between NC900 and C57BL/6J may have important implications in the interpretation of QTL mapping in our experiment. Therefore, we characterized heterozygosity and haplotype similarity in the F0 parents of our experimental population.
NC900 female genomewide segregation analysis and C57BL/6J sequence similarity.
We found no evidence of regions with more than two haplotypes segregating in the NC900 population. Therefore, it is simple to identify regions of IBD between NC900 and C57BL/6J independently of the heterozygosity status. Regions of IBD are found on every chromosome but are most prominent on chromosomes X, 10, and 16 (54, 43 and 41%, respectively) (Table 5).
We have combined the heterozygosity and IBD analyses to generate a map of the expected inheritance patterns in our experimental F2 population. At one extreme are homozygous regions in the NC900 females that are also IBD to C57BL/6J. In these regions (22% of the genome) polymorphisms should be limited to de novo mutations accumulated since the divergence of the two strains. On the other extreme are regions of heterozygosity in the NC900 females with two haplotypes that differ from C57BL/6J. In these regions, the mapping population can have as many as six different genotypes and their frequency depends on the MAF in the six NC900 females (Figure 5 and Additional File 2, Figure S1).
Local analysis in the QTL confidence intervals.
Our first objective was to test whether the aggression phenotype of NC900 mice is controlled primarily by a locus with large effect using a small NC900 × B6 F2 cross well powered to detect QTL with large effect. While no such result was identified, we were able to map QTL for social behavior traits despite a relative lack of power to detect loci with modest effect. A single QTL for aggression (percent duration count and attack count) was detected on MMU19 but had a paradoxical genetic effect. The B6 parental line, which rarely exhibits aggression in our social interaction test, contributed the high aggression allele. While QTL describing transgressive variation (i.e. in this case where the allele from the C57BL/6J, rather than the NC900 increased attack levels) are relatively common (Rieseberg, 1999), they are usually detected within a broad architecture dominated by QTL allelic effects that are in the expected direction based on parental phenotypic divergence [21–25]. Our results suggest that the NC900 mouse inter-male aggression phenotype is likely controlled by many additional QTL, each having effects too small to detect in the present F2 population.
Consideration of NC900 aggression as a complex trait is consistent with the genetic architecture of aggression in Drosophila where a minimum of 5 aggression QTL and extensive epistasis were detected . A cross of NC900 and NC100, while technically more difficult to evaluate (e.g. allele sharing), would be required to rule out the possibility that a genetic locus with a large effect was segregating in the ICR base population and contributed to the NC900 selection response, but is fixed in the same direction in the B6 line and thus would not be detected in this cross. The aggression QTL we detected do not correspond to those previously reported in mice [10, 11], which is not surprising for several reasons: 1) different methods of aggression behavior measurement; 2) different pre-aggression test housing environments; 3) laboratory-to-laboratory variation, and 4) different genetic background. It is particularly important to note that the aggression QTL we detected are potentially unique to the type of open field, social novelty aggression testing we used vs. resident intruder test . As an example, Roubertoux and colleagues found that differences in rearing and testing conditions produced stark differences in the genetic correlates of aggressive behavior .
Extensive dominance gene action was evident for aggression traits, and significant over/under-dominance effects for other social behavior traits were also evident. Strong heterosis may imply that these traits are relevant to overall fitness . From an evolutionary perspective, social behaviors like intermale aggression are likely important components of fitness, because they affect the ability to produce offspring. In a laboratory setting, dominant males sire > 90% of the offspring , but it is not known to what degree social behavior like aggression impacts fitness in wild mice. We have observed that overdominant gene action is common for other traits that affect overall fitness in mice such as litter size .
We observed qualitative differences in individual attack styles. Highly aggressive mice were noted for higher attack speed and more attacks described as front attacks. We did not attempt to assess individual differences in counts of front, side, and rear attacks in this study, because these types of attacks were realized retrospectively. We also lacked means to measure attack speed in the present study. Given the degree of the striking qualitative differences in attack speed and style we observed, such measures are warranted in the future. Future studies should also consider the possibility that the genotypes of the standard opponents indirectly interact with the genotypes of the subjects . We selected B6 standard opponents to in an effort to avoid confounding aggressive subject vs. aggressive social partner interactions, but did not account for F2 phenotypic variation that may be attributable to indirect genetic interactions between B6 standard opponent genotypes and various mixtures of NC900 × B6 F2 genotypes.
The peaks of QTL for attack latency and approach on MMU7 fall ~10 Mb from the tyrosinase gene (tyr). This is noteworthy because it has long been known that the tyr locus is associated with behavioral differences . Since NC900 mice are albino, and albinism is caused by a recessive point mutation in tyr, resulting in the absence of melanin production in the hair, skin, and eyes, it is conceivable that coat color represents a marker for attack latency and approach. Even though attack latency was never a NC900 high-aggression selection criterion, mean attack latencies significantly decreased across generations of selection for high levels of attack (unpublished data). We have reported effects of coat color on wheel running speed and also detected a QTL controlling mouse voluntary wheel running speed linked to tyr locus . In addition, effects of coat color on open-field activity have been reported [32, 33] and a contextual fear conditioning QTL has also been associated with the tyr locus . It is tempting to speculate that attack latency, wheel running speed, and fear conditioning QTL linked to the tyr locus share common functional variation manifested by control over the pace with which animals perform motivated behaviors like attack, fear conditioning, and wheel running. Alternatively, the attack latency QTL could relate to visual defects caused by the absence of melanin. For example, it is feasible that a lack of melanin could augment sensitivity to light in albino mice  resulting in an enhanced light-induced stress response. Thus, the faster attack latencies observed by F2 albino mice may be attributable to indirect effects of the tyr locus on the stress response rather than direct effects of tyr attack latency and approach tendencies per se.
High-density genotyping can be used to identify regions that are identical by descent (IBD) between strains used in mapping populations to refine QTL candidate intervals. Regions of IBD within QTL confidence intervals should be excluded as candidate QTL regions while regions with more than one haplotype remain. Based on our analysis, we can exclude 22% of the QTL confidence intervals indentified in this study.
QTL mapping was performed using methods that assume that our mapping population is analogous to a standard F2 population derived from two inbred strains. In contrast with these assumptions, our high-density genotype analysis shows that one-third of the genome is, in fact, segregating in a more complex pattern. This complex pattern poses challenges in QTL mapping especially in crosses with widely variable MAF distributions across the genome (Additional File 2, Figure S1). We suggest that crosses involving selection lines use a combination of markers that distinguish genetic variation both between and within the selection lines. Ignoring segregating regions is likely to result in overlooking loci that have been otherwise identified. This should be particularly relevant to populations in which heterozygosity has been selected for a particular phenotype. The size of the segregating regions, the presence of a very narrow bottleneck of only 3 breeding pairs (see Materials and Methods), and the distribution of MAF in the six NC900 females used in this study suggest that maintenance of heterozygosity may have been under selection during creation of the high-aggression NC900 line. Furthermore, the size of the segregating regions is smaller than in other selection lines derived from similar ICR stock and genotyped with the same high-density platform (FPMV and DP unpublished).
If a locus with large effect does not control aggression, how can the genetic basis of NC900 aggression be determined? It is important to recognize that the only NC900 selection criterion is number of attacks displayed in a novel social interaction test after isolation housing. Group housed NC900 animals display significantly less aggression than isolated NC900 mice . This gene by environment interaction is not unique to NC900 mice. Aggression phenotypes are strongly influenced by many environmental factors including group dynamics , maternal social relationships  and childhood maltreatment . Considering that predisposition to violence can require genotype (e.g. low-activity MAOA) and environmental (severe maltreatment) interactions , we should expect that aggression phenotypes pivot on complex relationships between genetic, environmental, and ontogenetic sources of variance . The challenge is to unravel the complexity. Given the reliability with which an environmental factor (i.e. isolation housing) can induce aggression in mice , it is an ideal phenotypic target for dissecting how a complex behavioral trait develops. But if aggression is controlled by many QTL that each have a small effect and whose influence is only apparent in particular environmental conditions, and is characterized by extensive epitasis , then very large populations of segregating genotypes along with contrasting environmental conditions will be needed to understand its genetic architecture.
The size of the segregating regions, the presence of a very narrow bottleneck of only 3 breeding pairs (see Materials and Methods), and the distribution of MAF in the six NC900 females used in this study suggest that maintenance of heterozygosity may have been under selection during the creation of the high-aggression NC900 lines. While lack of complete pedigrees and the limited sample size in this study prevents us from reaching definitive conclusions, this speculative hypothesis bears further study as such a result would have important ramifications for the nature of selection response for behavioral phenotypes.
The NC900 mouse line was selectively bred for high-aggression using only one selection criterion - attack counts displayed by an isolation housed 45 day old male toward a group-housed partner mouse in a 10-min novel social interaction test . At generation 55 (in 2004), 3 families of NC900 mice were successfully rederived to produce specific pathogen-free mice. Selective breeding for social interaction attack count did not proceed for every subsequent generation due to space and logistical constraints, but retention of phenotype was confirmed by periodically assessing isolated male NC900 attack counts against a group-housed C57BL/6J (B6; Jackson Laboratory, Bar Harbor, ME) mouse. At the suggestion of the UNC-Chapel Hill Institutional Animal Care and Use Committee, the age for standard social interaction testing was increased from 45 d to a minimum of 8 weeks. B6 mice were selected as social partners and as F2 founders because they rarely attack or display any observable form of aggression in our social interaction test (data not shown).
All mice were housed in standard cages on a 10:14 hr light/dark cycle (18:00 - 08:00 dark) and provided ad libitum access to feed and water. Mice were fed Prolab Isopro RMH 3000 (Lab Diet: protein 26%, fat 14%, carbohydrates 60%) through the experimental period. All procedures were conducted in accordance with NIH guidelines for the care and use of experimental animals and based on protocols approved by the Institutional Animal Care and Use Committee of UNC-Chapel Hill.
Six NC900 females having brothers that previously tested positive for high aggression (≥ 20 attacks within a 10 minute social interaction test against a B6 social partner) were mated with six B6 males in single mating pairs. Fifty-nine F1 animals were produced (28 females, 31 males). At 21 days of age all F1 males were weaned into isolation housing and tested at 8-9 weeks of age in the social interaction test against 8 to 10 week old B6 males. Eighteen of the 31 F1 males attacked at least once, and all 6 families produced attacking males. The 18 attacking males, plus five other males that did not attack within 5 minutes, but did attack within 10 minutes (as per the NC900 aggression selection criteria), were bred to randomly selected, non-sibling F1 females in 23 single mating groups. Twenty F1 mating groups produced 100 F2 males, which were weaned into isolation housing at 21 days of age.
Social behavior phenotyping
At 8-9 weeks of age, F1 and F2 males were tested against group-housed 8-10 week old B6 mice in the social interaction test according to the standard aggression selection procedures (Cairns et al., 1983), except for the social test duration. Social interaction tests were conducted in a 20 × 21 × 31 cm Plexiglas open field arena containing bed-o-cob bedding and lit to 360 Lux during 1900 - 2300 hours of the dark cycle. During the initial 2 minutes, the partner and subject animals habituated to opposing sides of the field separated by an opaque sliding divider. Upon removing the divider, mice interacted freely for 5 minutes. Black F2 subject mice and black B6 partner mice were differentiated by marking B6 partners with black permanent marker at the tips of their tails. Tests were recorded using a digital video camera (Sony DCR-SR45). Digital video files were coded using Noldus Observer XT™ software (Leesburg, VA). Subject coat color and subject and social partner body weights were recorded immediately after the social interaction test.
Attack behavior was defined as a vigorous lunge accompanied by a bite on the social partner animal . Aggressive behavior included: bite, chase, lunge, tail-rattling, and feints. Attack behavior was coded as a point variable within states of aggression. The social interaction coding system was developed in three stages. First, we performed ad-hoc scoring of NC900 and B6 social interactions behaviors in an effort to identify all social behaviors exhibited by the parental strains. Next, we developed a pilot coding system that could capture the range of behaviors observed in the ad-hoc observations. This pilot coding system was then further refined and tested by two coders using a subset of the F2 social interaction tests in an iterative manner until a high-degree of inter-coder reliability was established for all behavior codes (Cohen's κ > 0.80; ) using reliability statistic analyses included in Observer XT. The entire F2 population was coded by a single individual who was periodically tested for maintenance of intra-coder reliability against coded video standards (Cohen's κ > 0.90) using Observer XT throughout the F2 social interaction-coding period.
Definition of social interaction phenotypes
a Affiliative - sum of bidirectional, unidirectional, and approach behavior
simultaneous display of affiliative behaviors by the subject and partner mouse (e.g. simultaneous facial sniffing, grooming, anogenital sniffing), irrespective of which mouse initiated the contact
affiliative behavior displayed by subject with no observable response from partner
subject walks towards or follows partner (not coded if attack behavior immediately follows)
no subject-partner interaction
a Aversive - sum of passive avoidance, active avoidance, freezing, and aggression
a passive avoidance
partner-initiated affiliative social contact passively ignored by subject without moving away
a active avoidance
partner-initiated affiliative social contact actively ignored by subject by moving away
crouched, prone, & immobile posture lasting ≥ 1.0 second
vigorous lunge and bite typically directed at the partners' flanks and back, but also including lunges, feints, chasing, tail-rattling, and bites without lunge.
vigorous lunge and bite
amount of time expired from beginning of the social interaction test till the first attack
Low-density Genotyping and Linkage Map
A total of 100 F2 mice and eight F0 parental mice (six NC900 dams and two B6 sires) were genotyped for 176 SNPs using matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS). SNPs were initially selected based on their relatively even spacing across the genome and their predicted complete informativeness between NC900 and B6 mice, using data from the Wellcome-CTC Mouse Strain SNP Genotype Set http://www.well.ox.ac.uk/mouse/INBREDS. Predicted informativeness was based on genotypes from 8 mice representing lines (M16, ICR; ) derived from the same general genetic background as the ICR population used as the base for selection of the NC900 lines, relative to the genotype of B6. After genotyping, we excluded SNPs showing allele sharing across NC900 and B6 parents, and SNPs whose F2 genotypic frequencies significantly departed from the χ2 distribution based on Mendelian expectation. The final set of 154 SNPs used for QTL analyses is provided in Additional File 3, Table S2.
Eighty-eight of the 100 F2 subjects produced were analyzed; six were excluded because B6 partners displayed aggression towards the subject and six additional F2 subjects were excluded due to incomplete genotype data. Histograms of social behavior were plotted to determine normality of traits using SPSS 16.0 (Chicago, IL). Social behavior trait means between aggressive and non-aggressive subpopulations of the F2 mice were compared by t-tests, and the effects of coat color on social behavior traits were assessed by ANOVA with Tukey-adjusted p-values using SPSS 16.0. Social behavior Pearson partial correlations with Bonferonni corrections were generated with SAS 9.1 (Cary, NC).
Genome-wide QTL scans and stepwise model selections were performed using the R/qtl package in the R 2.7.2 environment . We used the stepwise model selection procedure to determine the influences of fixed effects (experimental batch), random effects (dam) and covariates (litter size and number of brothers and sisters, subject weight and partner body weight). Since none of the social behaviors were significantly influenced by any of these effects or covariates, they were not included in the single QTL model genome-wide scans. All LOD significance thresholds were determined by permutation , and LOD scores exceeding the permuted 95th and 90th percentiles were deemed significant and suggestive, respectively. Additive and dominance effects were extracted using R/qtl.
The social behavior measures exhibited non-parametric distributions. Therefore we performed non-parametric QTL scans by specifying the non-parametric distribution of traits which is equivalent to Kruskal-Wallis test statistic . The ranks of the phenotypic values, rather than the phenotypic values themselves, were fitted into a standard linear regression model to extract the percent variation and p-values for the QTL.
High quality, high-molecular weight DNA was extracted from the six NC900 females used for generation of the F2 population used in this study using phenol-chlororform. Samples were normalized to 50 ng/ul, processed according to the Affymetrix Genome-Wide Human SNP Nsp/Sty Assay Kit 5.0/6.0 protocol and hybridized to the Affymetrix Mouse Diversity genotyping array at the Functional Genomics core at the University of North Carolina at Chapel Hill. A total of 549,000 SNPS were used in this study.
SNP genotype calling was performed as described previously . To identify regions segregating within the founders of the experimental population, we determined the frequency of H calls with 200 consecutive SNP windows in each one of the six NC900 females independently. Regions with > 2% of SNPs with heterozygous calls were deemed heterozygous based on the analysis of 101 fully and partially inbred strains from the Jackson Laboratory (FPMV unpublished). We used a similar approach to identify regions of homozygosity in the six NC900 females containing different haplotypes (i.e., regions segregating among NC900 females). We mapped the start and end positions of each segregating interval to the first heterozygous SNP for segregating regions within a mouse and to the first SNP with discordant genotypes in segregating regions among the six NC900 females. This approach provides a conservative estimate of the length and position of heterozygous regions. For each segregating region we also determined the MAF.
IBD between each of the six NC900 females and C57BL/6J was determined by analysis of SNP markers. Analysis was performed using a 100 SNP marker sliding window and a threshold of 98% identity to C57BL/6J. This threshold is based on the characterization of multiple sister strains and biological duplicates (FPMV unpublished). Boundaries of the regions of IBD were determined. QTL regions were converted from CM to bp using the Center for Genome Dynamics Mouse Map Converter .
We thank J. Steven Reznick for contributing to the cost of rederiving the NC900 mouse line selectively bred for high aggression. This work was partially supported by grants from the National Institutes of Health DK076050, MH090338 and GM076468.
- Dow MA, von Schilcher F: Aggression and mating success in Drosophila melanogaster. Nature. 1975, 254 (5500): 511-512. 10.1038/254511a0.View ArticlePubMedGoogle Scholar
- Anderson PK, Hill JL: Mus musculus: Experimental Induction of Territory Formation. Science. 1965, 148 (3678): 1753-1755. 10.1126/science.148.3678.1753.View ArticlePubMedGoogle Scholar
- Lorenz : On Aggression. 1966, New York: Harcourt, Brace and WorldGoogle Scholar
- Gray SJ, Jensen SP, Hurst JL: Effects of resource distribution on activity and territory defense in house mice, Mus domesticus. Animal Behaviour. 2002, 63 (3): 531-539. 10.1006/anbe.2001.1932.View ArticleGoogle Scholar
- van Oortmerssen GA, Bakker TC: Artificial selection for short and long attack latencies in wild Mus musculus domesticus. Behav Genet. 1981, 11 (2): 115-126. 10.1007/BF01065622.View ArticlePubMedGoogle Scholar
- Sandnabba NK: Selective breeding for isolation-induced intermale aggression in mice: associated responses and environmental influences. Behav Genet. 1996, 26 (5): 477-488. 10.1007/BF02359752.View ArticlePubMedGoogle Scholar
- Cairns RB, MacCombie DJ, Hood KE: A developmental-genetic analysis of aggressive behavior in mice: I. Behavioral outcomes. J Comp Psychol. 1983, 97 (1): 69-89. 10.1037/0735-7036.97.1.69.View ArticlePubMedGoogle Scholar
- Hannon RM, Kelly SA, Middleton KM, Kolb EM, Pomp D, Garland T: Phenotypic effects of the "mini-muscle" allele in a large HR × C57BL/6J mouse backcross. J Hered. 2008, 99 (4): 349-354. 10.1093/jhered/esn011.View ArticlePubMedGoogle Scholar
- Horvat S, Medrano JF: Interval mapping of high growth (hg), a major locus that increases weight gain in mice. Genetics. 1995, 139 (4): 1737-1748.PubMed CentralPubMedGoogle Scholar
- Brodkin ES, Goforth SA, Keene AH, Fossella JA, Silver LM: Identification of quantitative trait loci that affect aggressive behavior in mice. J Neurosci. 2002, 22 (3): 1165-1170.PubMedGoogle Scholar
- Roubertoux PL, Guillot PV, Mortaud S, Pratte M, Jamon M, Cohen-Salmon C, Tordjman S: Attack behaviors in mice: from factorial structure to quantitative trait loci mapping. Eur J Pharmacol. 2005, 526 (1-3): 172-185. 10.1016/j.ejphar.2005.09.026.View ArticlePubMedGoogle Scholar
- Schicknick H, Hoffmann HJ, Schneider R, Crusio WE: Genetic analysis of isolation-induced aggression in the mouse. III. Classical cross-breeding analysis of differences between two closely related inbred strains. Behav Neural Biol. 1993, 59 (3): 242-248. 10.1016/0163-1047(93)91010-K.View ArticlePubMedGoogle Scholar
- Edwards AC, Mackay TF: Quantitative trait loci for aggressive behavior in Drosophila melanogaster. Genetics. 2009, 182 (3): 889-897. 10.1534/genetics.109.101691.PubMed CentralView ArticlePubMedGoogle Scholar
- Catlett RH: An evaluation of methods for measuring fighting behaviour with special reference to Mus musculus. Animal Behaviour. 1961, 9: 98-104. 10.1016/0003-3472(61)90042-2.View ArticleGoogle Scholar
- Lang PJ, Bradley MM, Cuthbert BN: Emotion, motivation, and anxiety: brain mechanisms and psychophysiology. Biol Psychiatry. 1998, 44 (12): 1248-1263. 10.1016/S0006-3223(98)00275-3.View ArticlePubMedGoogle Scholar
- Bolles RC: Species-specific defensive reactions and avoidance learning. Psychological Review. 1970, 77 (1): 32-48. 10.1037/h0028589.View ArticleGoogle Scholar
- Blanchard RJ, Hebert MA, Ferrari PF, Palanza P, Figueira R, Blanchard DC, Parmigiani S: Defensive behaviors in wild and laboratory (Swiss) mice: the mouse defense test battery. Physiol Behav. 1998, 65 (2): 201-209. 10.1016/S0031-9384(98)00012-2.View ArticlePubMedGoogle Scholar
- Trainor BC, Workman JL, Jessen R, Nelson RJ: Impaired nitric oxide synthase signaling dissociates social investigation and aggression. Behav Neurosci. 2007, 121 (2): 362-369. 10.1037/0735-7044.121.2.362.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Tuathaigh CM, Babovic D, O'Sullivan GJ, Clifford JJ, Tighe O, Croke DT, Harvey R, Waddington JL: Phenotypic characterization of spatial cognition and social behavior in mice with 'knockout' of the schizophrenia risk gene neuregulin 1. Neuroscience. 2007, 147 (1): 18-27.View ArticlePubMedGoogle Scholar
- Brodkin ES, Hagemann A, Nemetski SM, Silver LM: Social approach-avoidance behavior of inbred mouse strains towards DBA/2 mice. Brain Res. 2004, 1002 (1-2): 151-157. 10.1016/j.brainres.2003.12.013.View ArticlePubMedGoogle Scholar
- Rieseberg LH, Archer MA, Wayne RK: Transgressive segregation, adaptation and speciation. Heredity. 1999, 83 (Pt 4): 363-372. 10.1038/sj.hdy.6886170.View ArticlePubMedGoogle Scholar
- Rocha JL, Eisen EJ, Siewerdt F, Van Vleck LD, Pomp D: A large-sample QTL study in mice: III. Reproduction. Mamm Genome. 2004, 15 (11): 878-886. 10.1007/s00335-004-2364-6.View ArticlePubMedGoogle Scholar
- Rocha JL, Eisen EJ, Van Vleck LD, Pomp D: A large-sample QTL study in mice: II. Body composition. Mamm Genome. 2004, 15 (2): 100-113. 10.1007/s00335-003-2308-6.View ArticlePubMedGoogle Scholar
- Rocha JL, Eisen EJ, Van Vleck LD, Pomp D: A large-sample QTL study in mice: I. Growth. Mamm Genome. 2004, 15 (2): 83-99. 10.1007/s00335-003-2312-x.View ArticlePubMedGoogle Scholar
- Allan MF, Eisen EJ, Pomp D: Genomic mapping of direct and correlated responses to long-term selection for rapid growth rate in mice. Genetics. 2005, 170 (4): 1863-1877. 10.1534/genetics.105.041319.PubMed CentralView ArticlePubMedGoogle Scholar
- Falconer DS, Mackay TF: Introduction to quantitative genetics. 1996, Harlow, UK: Longman Group LtdGoogle Scholar
- DeFries JC, McClearn GE: Social dominance and Darwinian fitness in the laboratory mouse. The Amerian Naturalist. 1970, 104: 408-411. 10.1086/282675.View ArticleGoogle Scholar
- Wilson AJ, Gelin U, Perron MC, Reale D: Indirect genetic effects and the evolution of aggression in a vertebrate system. Proc Biol Sci. 2009, 276 (1656): 533-541. 10.1098/rspb.2008.1193.PubMed CentralView ArticlePubMedGoogle Scholar
- DeFries JC: Pleiotropic effects of albinism on open field behaviour in mice. Nature. 1969, 221 (5175): 65-66. 10.1038/221065a0.View ArticlePubMedGoogle Scholar
- Yokoyama T, Silversides DW, Waymire KG, Kwon BS, Takeuchi T, Overbeek PA: Conserved cysteine to serine mutation in tyrosinase is responsible for the classical albino mutation in laboratory mice. Nucleic Acids Res. 1990, 18 (24): 7293-7298. 10.1093/nar/18.24.7293.PubMed CentralView ArticlePubMedGoogle Scholar
- Nehrenberg DL, Wang S, Hannon RM, Garland T, Pomp D: QTL underlying voluntary exercise in mice: interactions with the "mini muscle" locus and sex. J Hered. 2010, 101 (1): 42-53. 10.1093/jhered/esp066.PubMed CentralView ArticlePubMedGoogle Scholar
- Henderson ND, Turri MG, DeFries JC, Flint J: QTL analysis of multiple behavioral measures of anxiety in mice. Behav Genet. 2004, 34 (3): 267-293. 10.1023/B:BEGE.0000017872.25069.44.View ArticlePubMedGoogle Scholar
- Turri MG, Henderson ND, DeFries JC, Flint J: Quantitative trait locus mapping in laboratory mice derived from a replicated selection experiment for open-field activity. Genetics. 2001, 158 (3): 1217-1226.PubMed CentralPubMedGoogle Scholar
- Ponder CA, Huded CP, Munoz MB, Gulden FO, Gilliam TC, Palmer AA: Rapid selection response for contextual fear conditioning in a cross between C57BL/6J and A/J: behavioral, QTL and gene expression analysis. Behav Genet. 2008, 38 (3): 277-291. 10.1007/s10519-008-9203-6.View ArticlePubMedGoogle Scholar
- Thomas BB, Aramant RB, Sadda SR, Seiler MJ: Light response differences in the superior colliculus of albino and pigmented rats. Neurosci Lett. 2005, 385 (2): 143-147. 10.1016/j.neulet.2005.05.034.View ArticlePubMedGoogle Scholar
- Gariepy JL, Gendreau PL, Mailman RB, Tancer M, Lewis MH: Rearing conditions alter social reactivity and D1 dopamine receptors in high- and low-aggressive mice. Pharmacol Biochem Behav. 1995, 51 (4): 767-773. 10.1016/0091-3057(95)00028-U.View ArticlePubMedGoogle Scholar
- Higgins LA, Jones KM, Wayne ML: Quantitative genetics of natural variation of behavior in Drosophila melanogaster: the possible role of the social environment on creating persistent patterns of group activity. Evolution. 2005, 59 (7): 1529-1539.View ArticlePubMedGoogle Scholar
- Suomi SJ: Early determinants of behaviour: evidence from primate studies. Br Med Bull. 1997, 53 (1): 170-184.View ArticlePubMedGoogle Scholar
- Dodge KA, Bates JE, Pettit GS: Mechanisms in the cycle of violence. Science. 1990, 250 (4988): 1678-1683. 10.1126/science.2270481.View ArticlePubMedGoogle Scholar
- Caspi A, McClay J, Moffitt TE, Mill J, Martin J, Craig IW, Taylor A, Poulton R: Role of genotype in the cycle of violence in maltreated children. Science. 2002, 297 (5582): 851-854. 10.1126/science.1072290.View ArticlePubMedGoogle Scholar
- Cairns RB, Gariepy JL, Hood KE: Development, microevolution, and social behavior. Psychol Rev. 1990, 97 (1): 49-65. 10.1037/0033-295X.97.1.49.View ArticlePubMedGoogle Scholar
- Valzelli L: The "isolation syndrome" in mice. Psychopharmacologia. 1973, 31 (4): 305-320. 10.1007/BF00421275.View ArticlePubMedGoogle Scholar
- Cohen J: A coefficient of agreement for nominal scales. Educational and Psychological Measurement. 1960, 20: 2037-2046. 10.1177/001316446002000104.View ArticleGoogle Scholar
- Broman KW, Wu H, Sen S, Churchill GA: R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003, 19 (7): 889-890. 10.1093/bioinformatics/btg112.View ArticlePubMedGoogle Scholar
- Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138 (3): 963-971.PubMed CentralPubMedGoogle Scholar
- Kruglyak L, Lander ES: A nonparametric approach for mapping quantitative trait loci. Genetics. 1995, 139 (3): 1421-1428.PubMed CentralPubMedGoogle Scholar
- Yang H, Ding Y, Hutchins LN, Szatkiewicz J, Bell TA, Paigen BJ, Graber JH, de Villena FP, Churchill GA: A customized and versatile high-density genotyping array for the mouse. Nat Methods. 2009, 6 (9): 663-666. 10.1038/nmeth.1359.PubMed CentralView ArticlePubMedGoogle Scholar
- Shifman S, Bell JT, Copley RR, Taylor MS, Williams RW, Mott R, Flint J: A high-resolution single nucleotide polymorphism genetic map of the mouse genome. PLoS Biol. 2006, 4 (12): e395-10.1371/journal.pbio.0040395.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.