Chromosome identification and assignment of DNA clones in the dog using a red fox and dog comparative map

Chromosome-level genome assembly for giant panda provides novel insights into Carnivora chromosome evolution

BACKGROUND

Chromosome evolution is an important driver of speciation and species evolution. Previous studies have detected chromosome rearrangement events among different Carnivora species using chromosome painting strategies. However, few of these studies have focused on chromosome evolution at a nucleotide resolution due to the limited availability of chromosome-level Carnivora genomes. Although the de novo genome assembly of the giant panda is available, current short read-based assemblies are limited to moderately sized scaffolds, making the study of chromosome evolution difficult.

RESULTS

Here, we present a chromosome-level giant panda draft genome with a total size of 2.29 Gb. Based on the giant panda genome and published chromosome-level dog and cat genomes, we conduct six large-scale pairwise synteny alignments and identify evolutionary breakpoint regions. Interestingly, gene functional enrichment analysis shows that for all of the three Carnivora genomes, some genes located in evolutionary breakpoint regions are significantly enriched in pathways or terms related to sensory perception of smell. In addition, we find that the sweet receptor gene TAS1R2, which has been proven to be a pseudogene in the cat genome, is located in an evolutionary breakpoint region of the giant panda, suggesting that interchromosomal rearrangement may play a role in the cat TAS1R2 pseudogenization.

CONCLUSIONS

We show that the combined strategies employed in this study can be used to generate efficient chromosome-level genome assemblies. Moreover, our comparative genomics analyses provide novel insights into Carnivora chromosome evolution, linking chromosome evolution to functional gene evolution.

Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours

Strains of red fox (Vulpes vulpes) with markedly different behavioural phenotypes have been developed in the famous long-term selective breeding programme known as the Russian farm-fox experiment. Here we sequenced and assembled the red fox genome and re-sequenced a subset of foxes from the tame, aggressive and conventional farm-bred populations to identify genomic regions associated with the response to selection for behaviour. Analysis of the re-sequenced genomes identified 103 regions with either significantly decreased heterozygosity in one of the three populations or increased divergence between the populations. A strong positional candidate gene for tame behaviour was highlighted: SorCS1, which encodes the main trafficking protein for AMPA glutamate receptors and neurexins and suggests a role for synaptic plasticity in fox domestication. Other regions identified as likely to have been under selection in foxes include genes implicated in human neurological disorders, mouse behaviour and dog domestication. The fox represents a powerful model for the genetic analysis of affiliative and aggressive behaviours that can benefit genetic studies of behaviour in dogs and other mammals, including humans.

Genotyping-By-Sequencing (GBS) Detects Genetic Structure and Confirms Behavioral QTL in Tame and Aggressive Foxes (Vulpes vulpes)

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.

Sequence comparison of prefrontal cortical brain transcriptome from a tame and an aggressive silver fox (Vulpes vulpes)

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.

Chromosomal mapping of canine-derived BAC clones to the red fox and American mink genomes

High-quality sequencing of the dog (Canis lupus familiaris) genome has enabled enormous progress in genetic mapping of canine phenotypic variation. The red fox (Vulpes vulpes), another canid species, also exhibits a wide range of variation in coat color, morphology, and behavior. Although the fox genome has not yet been sequenced, canine genomic resources have been used to construct a meiotic linkage map of the red fox genome and begin genetic mapping in foxes. However, a more detailed gene-specific comparative map between the dog and fox genomes is required to establish gene order within homologous regions of dog and fox chromosomes and to refine breakpoints between homologous chromosomes of the 2 species. In the current study, we tested whether canine-derived gene-containing bacterial artificial chromosome (BAC) clones can be routinely used to build a gene-specific map of the red fox genome. Forty canine BAC clones were mapped to the red fox genome by fluorescence in situ hybridization (FISH). Each clone was uniquely assigned to a single fox chromosome, and the locations of 38 clones agreed with cytogenetic predictions. These results clearly demonstrate the utility of FISH mapping for construction of a whole-genome gene-specific map of the red fox. The further possibility of using canine BAC clones to map genes in the American mink (Mustela vison) genome was also explored. Much lower success was obtained for this more distantly related farm-bred species, although a few BAC clones were mapped to the predicted chromosomal locations.

A meiotic linkage map of the silver fox, aligned and compared to the canine genome

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.

FISH mapping of 10 canine BAC clones harbouring genes and microsatellites in the arctic fox and the Chinese raccoon dog genomes

Cytogenetic mapping of the arctic fox and the Chinese raccoon dog were performed using a set of canine probes derived from the Bacterial Artificial Chromosome (BAC) library. Altogether, 10 BAC clones containing sequences of selected genes (PAX3, HBB, ATP2A2, TECTA, PIT1, ABCA4, ESR2, TPH1, HTR2A, MAOA) and microsatellites were mapped by fluorescence in situ hybridization (FISH) experiments to chromosomes of the canids studied. At present, the cytogenetic map on the arctic fox and Chinese raccoon dog consists of 45 loci each. Chromosomal localization of the BAC clones was in agreement with data obtained by earlier independent comparative chromosome painting. However, two events of telomere-to-centromere inversions were tentatively identified while compared with assignments in the dog karyotype.

Cross-species chromosome painting unveils cytogenetic signatures for the Eulipotyphla and evidence for the polyphyly of Insectivora

Insectivore-like animals are traditionally believed among the first eutherian mammals that have appeared on the earth. The modern insectivores are thus crucial for understanding the systematics and phylogeny of eutherian mammals as a whole. Here cross-species chromosome painting, with probes derived from flow-sorted chromosomes of human, was used to delimit the homologous chromosomal segments in two Soricidae species, the common shrew (Sorex araneus, 2n = 20/21), and Asiatic short-tailed shrew (Blarinella griselda, 2n = 44), and one Erinaceidae species, the shrew-hedgehog (Neotetracus sinensis, 2n = 32), and human. We report herewith the first comparative maps for the Asiatic short-tailed shrew and the shrew-hedgehog, in addition to a refined comparative map for the common shrew. In total, the 22 human autosomal paints detected 40, 51 and 58 evolutionarily conserved segments in the genomes of common shrew, Asiatic short-tailed shrew, and shrew-hedgehog, respectively, demonstrating that the common shrew has retained a conserved genome organization while the Asiatic short-tailed shrew and shrew-hedgehog have relatively rearranged genomes. In addition to confirming the existence of such ancestral human segmental combinations as HSA 3/21, 12/22, 14/15 and 7/16 that are shared by most eutherian mammals, our study reveals a shared human segmental combination, HSA 4/20, that could phylogenetically unite the Eulipotyphlan (i.e., the core insectivores) species. Our results provide cytogenetic evidence for the polyphyly of the order Insectivora and additional data for the eventual reconstruction of the ancestral eutherian karyotype.

Karyotypic conservatism in the suborder Feliformia (Order Carnivora)

Multidirectional comparative chromosome painting was used to investigate the karyotypic relationships among representative species from three Feliformia families of the order Carnivora (Viverridae, Hyaenidae and Felidae). Complete sets of painting probes derived from flow-sorted chromosomes of the domestic dog, American mink, and human were hybridized onto metaphases of the spotted hyena (Crocuta crocuta, 2n = 40) and masked palm civet (Paguma larvata, 2n = 44). Extensive chromosomal conservation is evident in these two species when compared with the cat karyotype, and only a few events of chromosome fusion, fission and inversion differentiate the karyotypes of these Feliformia species. The comparative chromosome painting data have enabled the integration of the hyena and palm civet chromosomes into the previously established comparative map among the domestic cat, domestic dog, American mink and human and improved our understanding on the karyotype phylogeny of Feliformia species.

Comparative chromosome painting defines the karyotypic relationships among the domestic dog, Chinese raccoon dog and Japanese raccoon dog

The Chinese raccoon dog (Nyctereutes procyonoides procyonoides, 2n = 54 + 2-3 B) and Japanese raccoon dog (Nyctereutes p. viverrinus, 2n = 38 + 3-4 B) are two subspecies of the same species. The genome-wide comparative chromosome map between the Japanese raccoon dog and domestic dog (Canis familiaris) has been established by fluorescence in-situ hybridization with a set of domestic dog painting probes. In this study, we established the comparative chromosome map for the Chinese raccoon dog and domestic dog. In total, dog probes specific for the 38 autosomes delineated 41 conserved chromosomal segments in the Chinese raccoon dog. Probes from dog chromosomes 1, 13 and 19 each painted two Chinese raccoon dog chromosome segments. Fifteen dog autosomal probes each hybridized to one Chinese raccoon dog chromosome, while each of the other dog autosomal probes painted to a single Chinese raccoon dog chromosomal arm. Dog X chromosome probe delineated the entire X chromosome of the Chinese raccoon dog; the dog Y chromosome probe hybridized to the pseudoautosomal region at the Xpter as well as the entire Y chromosome of the Chinese raccoon dog. Comparative analysis of the distribution patterns of conserved segments defined by dog paints in the genomes of the Chinese and Japanese raccoon dogs demonstrates that their differences in the karyotypes of these two subspecies could have resulted from eight Robertsonian translocations. The large difference in chromosome number between the Chinese and Japanese raccoon dogs suggests that they should be considered as two distinct species.

Radiation hybrid map, physical map, and low-pass genomic sequence of the canine prcd region on CFA9 and comparative mapping with the syntenic region on human chromosome 17

Progressive rod-cone degeneration (prcd) is a canine retinal disease that maps to the centromeric end of CFA9 in a region of synteny with the distal part of HSA17q. As such, prcd has been postulated as the only animal model of RP17, a human retinitis pigmentosa locus that maps to 17q22. In an effort to establish more detailed regions of synteny between dog CFA9 and the HSA17q-ter region, we created a robust gene-enriched CFA9-RH08(3000) map with 34 gene-based markers and 12 microsatellites, with the highest resolution and number of markers for the centromeric end of CFA9. Furthermore, we built an approximately 1.5-Mb physical map containing both GRB2 and GALK1, genes so far identified by meiotic linkage analysis as being closest to the prcd locus, and generated about 1.2 Mb low-pass (3.2x) canine sequence. Canine to human comparative sequence analysis identified 49 transcripts that had been previously mapped to the HSA17q25 region. The generated low-pass canine sequence was annotated with a working draft of human sequence from HSA17q25, and we used this scaffold to order and orient the canine sequence against human. This order and orientation are preliminary, as high-throughput genomic sequencing of HSA17q-ter has not been fully completed.

Reciprocal chromosome painting illuminates the history of genome evolution of the domestic cat, dog and human

Domestic cats and dogs are important companion animals and model animals in biomedical research. The cat has a highly conserved karyotype, closely resembling the ancestral karyotype of mammals, while the dog has one of the most extensively rearranged mammalian karyotypes investigated so far. We have constructed the first detailed comparative chromosome map of the domestic dog and cat by reciprocal chromosome painting. Dog paints specific for the 38 autosomes and the X chromosomes delineated 68 conserved chromosomal segments in the cat, while reverse painting of cat probes onto red fox and dog chromosomes revealed 65 conserved segments. Most conserved segments on cat chromosomes also show a high degree of conservation in G-banding patterns compared with their canine counterparts. At least 47 chromosomal fissions (breaks), 25 fusions and one inversion are needed to convert the cat karyotype to that of the dog, confirming that extensive chromosome rearrangements differentiate the karyotypes of the cat and dog. Comparative analysis of the distribution patterns of conserved segments defined by dog paints on cat and human chromosomes has refined the human/cat comparative genome map and, most importantly, has revealed 15 cryptic inversions in seven large chromosomal regions of conserved synteny between humans and cats.

A comparative chromosome map of the Arctic fox, red fox and dog defined by chromosome painting and high resolution G-banding

A complete set of paint probes, with each probe specific for a single type of dog chromosome, was generated by DOP-PCR amplification of flow-sorted chromosomes. These probes have been assigned to high-resolution G-banded chromosomes of the dog and Arctic fox by fluorescence in-situ hybridization. On the basis of these results we propose improved nomenclature for the G-banded karyotypes of the dog and Artic fox. A comparative map between the Arctic fox, red fox and dog has been established based on results from chromosome painting and high-resolution G-banding. This map demonstrates that the euchromatic complements of these three canid species consists of 42 conserved segments. Thirty-four of these 42 segments are each represented by a single dog chromosome with dog chromosomes 1, 13, 18 and 19 each retaining two segments, respectively. The autosomes of the Arctic fox and red fox could be reconstructed from these 42 blocks in different combinations through chromosomal fusions. Our findings suggest that chromosome fusion has been the principal mechanism of karyotype evolution occuring during speciation in canids.

Use of flow-sorted canine chromosomes in the assignment of canine linkage, radiation hybrid, and syntenic groups to chromosomes: refinement and verification of the comparative chromosome map for dog and human

The mapping of the canine genome has recently been accelerated by the availability of chromosome-specific reagents and publication of radiation hybrid (RH), genetic linkage, and dog/human comparative maps, but the assignment of mapping groups to chromosomes is incomplete. To assign published radiation hybrid, linkage, and "syntenic" groups to chromosomes, individual markers found within each group have been amplified from canine and vulpine flow-sorted, chromosome-specific DNAs as templates. Here a further 102 type I genetic markers (previously mapped in human) and 21 further type II markers are assigned to canine chromosomes using marker-specific PCR. We have assigned all linkage, RH, and syntenic groups in the two most recently published canine genome maps to chromosomes. This demonstrates directly that there is at least one published mapping group for each of the 38 canine autosomes and thus that the coverage of the canine chromosome map is approaching completion. The dog/human comparative map is one of the most complex so far described, with 90 separate segments of chromosomal homology previously seen in dog-on-human cross-species chromosome-painting studies. The total of 142 type I markers now placed on canine chromosomes using this method of marker mapping has allowed us to confirm the placement of the great majority (83) of the 90 homologous segments. The positions of the remaining homologous segments were confirmed in new cross-species chromosome-painting experiments (dog-on-human, fox-on-human).