Search

Individuals involved in a social interaction exhibit different behavioral traits that, in combination, form the individual’s behavioral responses. Selectively bred strains of silver foxes (Vulpes vulpes) demonstrate markedly different behaviors in their response to humans. To identify the genetic basis of these behavioral differences we constructed a large F2 population including 537 individuals by cross-breeding tame and aggressive fox strains. 98 fox behavioral traits were recorded during social interaction with a human experimenter in a standard four-step test. Patterns of fox behaviors during the test were evaluated using principal component (PC) analysis. Genetic mapping identified eight unique significant and suggestive QTL. Mapping results for the PC phenotypes from different test steps showed little overlap suggesting that different QTL are involved in regulation of behaviors exhibited in different behavioral contexts. Many individual behavioral traits mapped to the same genomic regions as PC phenotypes. This provides additional information about specific behaviors regulated by these loci. Further, three pairs of epistatic loci were also identified for PC phenotypes suggesting more complex genetic architecture of the behavioral differences between the two strains than what has previously been observed.

The foxes at Novosibirsk, Russia, are the only population of domesticated foxes in the world. These domesticated foxes originated from farm-bred silver foxes (Vulpes vulpes), whose genetic source is unknown. In this study we examined the origin of the domesticated strain of foxes and two other farm-bred fox populations (aggressive and unselected) maintained in Novosibirsk. To identify the phylogenetic origin of these populations we sequenced two regions of mtDNA, cytochrome b and D-loop, from 24 Novosibirsk foxes (8 foxes from each population) and compared them with corresponding sequences of native red foxes from Europe, Asia, Alaska and Western Canada, Eastern Canada, and the Western Mountains of the USA. We identified seven cytochrome b - D-loop haplotypes in Novosibirsk populations, four of which were previously observed in Eastern North America. The three remaining haplotypes differed by one or two base change from the most common haplotype in Eastern Canada. Φ(ST) analysis showed significant differentiation between Novosibirsk populations and red fox populations from all geographic regions except Eastern Canada. No haplotypes of Eurasian origin were identified in the Novosibirsk populations. These results are consistent with historical records indicating that the original breeding stock of farm-bred foxes originated from Prince Edward Island, Canada. Mitochondrial DNA data together with historical records indicate two stages in the selection of domesticated foxes: the first includes captive breeding for ~50 years with unconscious selection for behaviour; the second corresponds to over 50 further years of intensive selection for tame behaviour.

The silver fox (Vulpes vulpes) offers a novel model for studying the genetics of social behavior and animal domestication. Selection of foxes, separately, for tame and for aggressive behavior has yielded two strains with markedly different, genetically determined, behavioral phenotypes. Tame strain foxes are eager to establish human contact while foxes from the aggressive strain are aggressive and difficult to handle. These strains have been maintained as separate outbred lines for over 40 generations but their genetic structure has not been previously investigated. We applied a genotyping-by-sequencing (GBS) approach to provide insights into the genetic composition of these fox populations. Sequence analysis of EcoT22I genomic libraries of tame and aggressive foxes identified 48,294 high quality SNPs. Population structure analysis revealed genetic divergence between the two strains and more diversity in the aggressive strain than in the tame one. Significant differences in allele frequency between the strains were identified for 68 SNPs. Three of these SNPs were located on fox chromosome 14 within an interval of a previously identified behavioral QTL, further supporting the importance of this region for behavior. The GBS SNP data confirmed that significant genetic diversity has been preserved in both fox populations despite many years of selective breeding. Analysis of SNP allele frequencies in the two populations identified several regions of genetic divergence between the tame and aggressive foxes, some of which may represent targets of selection for behavior. The GBS protocol used in this study significantly expanded genomic resources for the fox, and can be adapted for SNP discovery and genotyping in other canid species.

The objective of this study was to develop a systematic and flexible method for assembling multiplex simple sequence repeat (SSR) marker panels for high throughput genome analysis in rice, Oryza sativa, and to test these panels on a set of cultivated rice germplasm. To do this, 159 microsatellite markers were fluorescently labeled and assembled into 21 multiplex panels for semiautomated genotyping, providing genome-wide coverage of the 12 rice chromosomes. Panels are comprised of an average of eight markers each, occurring at approximately 11-centimorgan (cM) intervals throughout the genome. On a standard set of 13 genetically diverse cultivars of Oryza, these markers detected an average of five alleles per locus and had a mean polymorphism information content (P.I.C. value) of 0.67. Polymerase chain reactions (PCR) were optimized on a per marker basis to generate a uniform amount of PCR product and each primer pair was assessed in replicated trials for reliability of allele size estimates. T4 DNA polymerase was used to treat PCR products where the standard deviation of allele molecular weight was greater than 0.5 base pairs (bp). This treatment minimized the variance so that, in the multiplex set reported here, the average std. devJmarker was 0.24 bp, allowing accurate discrimination of alleles that differed by a single nucleotide. The resulting data on allele sizes were then entered into GeneFlow analysis software for the evaluation of polymorphism patterns among diverse rice cultivars. The use of an automated software tool for designing multiplex panels on the basis of both highly polymorphic and more conservative SSR markers resulted in the development of a highly informative semiautomated genotyping system for applications in rice genetics and breeding.

A new and unusual family of LTR elements, Dasheng, has been discovered in the genome of Oryza sativa following database searches of ~100 Mb of rice genomic sequence and 78 Mb of BAC-end sequence information. With all of the cis-elements but none of the coding domains normally associated with retrotransposons (e.g., gag, pol), Dasheng is a novel nonautonomous LTR element with high copy number. Over half of the ~1000 Dasheng elements in the rice genome are full length (5.6–8.6 kb), and 60% are estimated to have amplified in the past 500,000 years. Using a modified AFLP technique called transposon display, 215 elements were mapped to all 12 rice chromosomes. Interestingly, more than half of the mapped elements are clustered in the heterochromatic regions around centromeres. The distribution pattern was further confirmed by FISH analysis. Despite clustering in heterochromatin, Dasheng elements are not nested, suggesting their potential value as molecular markers for these marker-poor regions. Taken together, Dasheng is one of the highest-copy-number LTR elements and one of the most recent elements to amplify in the rice genome.

A total of 57.8 Mb of publicly available rice (Oryza sativaL.) DNA sequence was searched to determine the frequency and distribution of different simple sequence repeats (SSRs) in the genome. SSR loci were categorized into two groups based on the length of the repeat motif. Class I, or hypervariable markers, consisted of SSRs ≥20 bp, and Class II, or potentially variable markers, consisted of SSRs ≥12 bp EcoRI- and HindIII-digested BAC clones was one SSR per 40 Kb, whereas in continuous genomic sequence (represented by 27 fully sequenced BAC and PAC clones), the frequency was one SSR every 16 kb. Class II SSRs were estimated to occur every 3.7 kb in BAC ends and every 1.9 kb in fully sequenced BAC and PAC clones. GC-rich trinucleotide repeats (TNRs) were most abundant in protein-coding portions of ESTs and in fully sequenced BACs and PACs, whereas AT-rich TNRs showed no such preference, and di- and tetranucleotide repeats were most frequently found in noncoding, intergenic regions of the rice genome. Microsatellites with poly(AT)n repeats represented the most abundant and polymorphic class of SSRs but were frequently associated with the Micropon family of miniature inverted-repeat transposable elements (MITEs) and were difficult to amplify. A set of 200 Class I SSR markers was developed and integrated into the existing microsatellite map of rice, providing immediate links between the genetic, physical, and sequence-based maps. This contribution brings the number of microsatellite markers that have been rigorously evaluated for amplification, map position, and allelic diversity in Oryza spp. to a total of 500.[Clone sequences for 199 markers (RM1–RM88, RM200–RM345) developed in this lab are available as GenBank accessions AF343840–AF343869 andAF344003–AF344169.]

The growing number of rice microsatellite markers warrants a comprehensive comparison of allelic variability between the markers developed using different methods, with various sequence repeat motifs, and from coding and non-coding portions of the genome. We have performed such a comparison over a set of 323 microsatellite markers; 194 were derived from genomic library screening and 129 were derived from the analysis of rice-expressed sequence tags (ESTs) available in public DNA databases. We have evaluated the frequency of polymorphism between parental pairs of six inter- subspecific crosses and one inter-specific cross widely used for mapping in rice. Microsatellites derived from genomic libraries detected a higher level of polymorphism than those derived from ESTs contained in the GenBank database (83.8% versus 54.0%). Similarly, the other measures of genetic variability [the number of alleles per locus, polymorphism information content (PIC), and allele size ranges] were all higher in genomic library-derived microsatellites than in their EST-database counterparts. The highest overall degree of genetic diversity was seen in GA-containing microsatellites of genomic library origin, while the most conserved markers contained CCG- or CAG-trinucleotide motifs and were developed from GenBank sequences. Preferential location of specific motifs in coding versus non-coding regions of known genes was related to observed levels of microsatellite diversity. A strong positive correlation was observed between the maximum length of a microsatellite motif and the standard deviation of the molecular-weight of amplified fragments. The reliability of molecular weight standard deviation (SDmw) as an indicator of genetic variability of microsatellite loci is discussed.

In order to enhance the resolution of an existing genetic map of rice, and to obtain a comprehensive picture of marker utility and genomic distribution of microsatellites in this important grain species, rice DNA sequences containing simple sequence repeats (SSRs) were extracted from several small-insert genomic libraries and from the database. One hundred and eighty eight new microsatellite markers were developed and evaluated for allelic diversity. The new simple sequence length polymorphisms (SSLPs) were incorporated into the existing map previously containing 124 SSR loci. The 312 microsatellite markers reported here provide whole-genome coverage with an average density of one SSLP per 6 cM. In this study, 26 SSLP markers were identified in published sequences of known genes, 65 were developed based on partial cDNA sequences available in GenBank, and 97 were isolated from genomic libraries. Microsatellite markers with different SSR motifs are relatively uniformly distributed along rice chromosomes regardless of whether they were derived from genomic clones or cDNA sequences. However, the distribution of polymorphism detected by these markers varies between different regions of the genome.

Ninety-four newly developed microsatellite markers were integrated into existing RFLP framework maps of four rice populations, including two doubled haploid, a recombinant inbred, and an interspecific backcross population. These simple sequence repeats (SSR) were predominantly poly(GA) motifs, targetted because of their abundance in rice. They were isolated from a previously described sheared library and a newly constructed enzyme-digested library. Differences in the average length of poly(GA) tracts were observed for clones isolated from the two libraries. The length of GA motifs averaged 21 repeat units for clones isolated from the Tsp-509-digested library, while motifs averaged 17 units for clones from the sheared library. There was no evidence of clustering of microsatellite markers near centromeres or telomeres. Mapping of the 94 newly developed markers as well as of 27 previously reported microsatellites provided genome-wide coverage of the 12 chromosomes, with an average distance of 1 SSLP (simple sequence repeat polymorphism) per 16–20 cM.

Background Two strains of the silver fox (Vulpes vulpes), with markedly different behavioral phenotypes, have been developed by long-term selection for behavior. Foxes from the tame strain exhibit friendly behavior towards humans, paralleling the sociability of canine puppies, whereas foxes from the aggressive strain are defensive and exhibit aggression to humans. To understand the genetic differences underlying these behavioral phenotypes fox-specific genomic resources are needed. Results cDNA from mRNA from pre-frontal cortex of a tame and an aggressive fox was sequenced using the Roche 454 FLX Titanium platform (> 2.5 million reads & 0.9 Gbase of tame fox sequence; >3.3 million reads & 1.2 Gbase of aggressive fox sequence). Over 80% of the fox reads were assembled into contigs. Mapping fox reads against the fox transcriptome assembly and the dog genome identified over 30,000 high confidence fox-specific SNPs. Fox transcripts for approximately 14,000 genes were identified using SwissProt and the dog RefSeq databases. An at least 2-fold expression difference between the two samples (p < 0.05) was observed for 335 genes, fewer than 3% of the total number of genes identified in the fox transcriptome. Conclusions Transcriptome sequencing significantly expanded genomic resources available for the fox, a species without a sequenced genome. In a very cost efficient manner this yielded a large number of fox-specific SNP markers for genetic studies and provided significant insights into the gene expression profile of the fox pre-frontal cortex; expression differences between the two fox samples; and a catalogue of potentially important gene-specific sequence variants. This result demonstrates the utility of this approach for developing genomic resources in species with limited genomic information.

The silver fox provides a rich resource for investigating the genetics of behavior, with strains developed by intensely selective breeding that display markedly different behavioral phenotypes. Until recently, however, the tools for conducting molecular genetic investigations in this species were very limited. In this review, the history of development of this resource and the tools to exploit it are described. Although the focus is on the genetics of domestication in the silver fox, there is a broader context. In particular, one expectation of the silver fox research is that it will be synergistic with studies in other species, including humans, to yield a more comprehensive understanding of the molecular mechanisms and evolution of a wider range of social cognitive behaviors.

The foxes at Novosibirsk, Russia, are the only population of domesticated foxes in the world. These domesticated foxes originated from farm-bred silver foxes (Vulpes vulpes), whose genetic source is unknown. In this study we examined the origin of the domesticated strain of foxes and two other farm-bred fox populations (aggressive and unselected) maintained in Novosibirsk. To identify the phylogenetic origin of these populations we sequenced two regions of mtDNA, cytochrome b and D-loop, from 24 Novosibirsk foxes (8 foxes from each population) and compared them with corresponding sequences of native red foxes from Europe, Asia, Alaska and Western Canada, Eastern Canada, and the Western Mountains of the USA. We identified seven cytochrome b - D-loop haplotypes in Novosibirsk populations, four of which were previously observed in Eastern North America. The three remaining haplotypes differed by one or two base change from the most common haplotype in Eastern Canada. Φ(ST) analysis showed significant differentiation between Novosibirsk populations and red fox populations from all geographic regions except Eastern Canada. No haplotypes of Eurasian origin were identified in the Novosibirsk populations. These results are consistent with historical records indicating that the original breeding stock of farm-bred foxes originated from Prince Edward Island, Canada. Mitochondrial DNA data together with historical records indicate two stages in the selection of domesticated foxes: the first includes captive breeding for ~50 years with unconscious selection for behaviour; the second corresponds to over 50 further years of intensive selection for tame behaviour.

During the second part of the twentieth century, Belyaev selected tame and aggressive foxes (Vulpes vulpes), in an effort known as the “farm-fox experiment”, to recapitulate the process of animal domestication. Using these tame and aggressive foxes as founders of segregant backcross and intercross populations we have employed interval mapping to identify a locus for tame behavior on fox chromosome VVU12. This locus is orthologous to, and therefore validates, a genomic region recently implicated in canine domestication. The tame versus aggressive behavioral phenotype was characterized as the first principal component (PC) of a PC matrix made up of many distinct behavioral traits (e.g. wags tail; comes to the front of the cage; allows head to be touched; holds observer’s hand with its mouth; etc.). Mean values of this PC for F1, backcross and intercross populations defined a linear gradient of heritable behavior ranging from tame to aggressive. The second PC did not follow such a gradient, but also mapped to VVU12, and distinguished between active and passive behaviors. These data suggest that (1) there are at least two VVU12 loci associated with behavior; (2) expression of these loci is dependent on interactions with other parts of the genome (the genome context) and therefore varies from one crossbred population to another depending on the individual parents that participated in the cross.

Microsatellites are simple, tandemly repeated di- to tetra-nucleotide sequence motifs flanked by unique sequences. They are valuable as genetic markers because they are co-dominant, detect high levels of allelic diversity, and are easily and economically assayed by the polymerase chain reaction (PCR). Results from screening a rice genomic library suggest that there are an estimated 5700-10 000 microsatellites in rice, with the relative frequency of different repeats decreasing with increasing size of the motif. A map consisting of 120 microsatellite markers demonstrates that they are well distributed throughout the 12 chromosomes of rice. Five multiple copy primer sequences have been identified that could be mapped to independent chromosomal locations. The current level of genome coverage provided by these simple sequence length polymorphisms (SSLPs) in rice is sufficient to be useful for genotype identification, gene and quantitative trait locus (QTL) analysis, screening of large insert libraries, and marker-assisted selection in breeding. Studies of allelic diversity have documented up to 25 alleles at a single locus in cultivated rice germplasm and provide evidence that amplification in wild relatives of Oryza sativa is generally reliable. The availability of increasing numbers of mapped SSLP markers can be expected to complement existing RFLP and AFLP maps, increasing the power and resolution of genome analysis in rice.

A meiotic linkage map is essential for mapping traits of interest and is often the first step toward understanding a cryptic genome. Specific strains of silver fox (a variant of the red fox, Vulpes vulpes), which segregate behavioral and morphological phenotypes, create a need for such a map. One such strain, selected for docility, exhibits friendly dog-like responses to humans, in contrast to another strain selected for aggression. Development of a fox map is facilitated by the known cytogenetic homologies between the dog and fox, and by the availability of high resolution canine genome maps and sequence data. Furthermore, the high genomic sequence identity between dog and fox allows adaptation of canine microsatellites for genotyping and meiotic mapping in foxes. Using 320 such markers, we have constructed the first meiotic linkage map of the fox genome. The resulting sex-averaged map covers 16 fox autosomes and the X chromosome with an average inter-marker distance of 7.5 cM. The total map length corresponds to 1480.2 cM. From comparison of sex-averaged meiotic linkage maps of the fox and dog genomes, suppression of recombination in pericentromeric regions of the metacentric fox chromosomes was apparent, relative to the corresponding segments of acrocentric dog chromosomes. Alignment of the fox meiotic map against the 7.6x canine genome sequence revealed high conservation of marker order between homologous regions of the two species. The fox meiotic map provides a critical tool for genetic studies in foxes and identification of genetic loci and genes implicated in fox domestication.