Aggressive behavior displays a high heritability in our study group of Golden Retriever dogs. Alterations in brain serotonin metabolism have been described in aggressive dogs before. Here, we evaluate whether four genes of the canine serotonergic system, coding for the serotonin receptors 1A, 1B, and 2A, and the serotonin transporter, could play a major role in aggression in Golden Retrievers. We performed mutation screens, linkage analysis, an association study, and a quantitative genetic analysis. There was no systematic difference between the coding DNA sequence of the candidate genes in aggressive and non-aggressive Golden Retrievers. An affecteds-only parametric linkage analysis revealed no strong major locus effect on human-directed aggression related to the candidate genes. An analysis of 41 single nucleotide polymorphisms (SNPs) in the 1 Mb regions flanking the genes in 49 unrelated human-directed aggressive and 49 unrelated non-aggressive dogs did not show association of SNP alleles, genotypes, or haplotypes with aggression at the candidate loci. We completed our analyses with a study of the effect of variation in the candidate genes on a collection of aggression-related phenotypic measures. The effects of the candidate gene haplotypes were estimated using the Restricted Maximum Likelihood method, with the haplotypes included as fixed effects in a linear animal model. We observed no effect of the candidate gene haplotypes on a range of aggression-related phenotypes, thus extending our conclusions to several types of aggressive behavior. We conclude that it is unlikely that these genes play a major role in the variation in aggression in the Golden Retrievers that we studied. Smaller phenotypic effects of these loci could not be ruled out with our sample size.