Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.     

Identification and characterization of a set of 100 tri- and dinucleotide microsatellites in the canine genome

A set of 100 canine microsatellite markers--83 dinucleotides and 17 trinucleotides--is reported. A study of their frequency in the dog genome showed that, while the frequency of the CA repeats is one (CA)n every 47 kb, the 10 trinucleotidic frequencies vary from one every 117 kb (AGG)n to one every 875 kb (AGT)n. Polymorphism analysis performed on 16 unrelated mongrel dogs showed that 80% of dinucleotides are polymorphic, while only 30% of the trinucleotides are so. Of this set of 100 markers, 56 have been mapped on the RHDF5000 dog/hamster whole genome radiation hybrid panel. Moreover, through systematic BLAST analogy searches of the microsatellite-containing clone sequence, three new dog genes could be identified, based on their human ortholog. All of the markers presented may prove useful in physical mapping methods, and polymorphic microsatellites in genetic linkage studies or parentage controls in dog.     

A comparison of complete genome sequences of a rabies virus chinese isolate SH06 with the vaccine strains

In this study, we determined the complete nucleotide and deduced amino acid sequence of a primary isolate of rabies virus (SH06) obtained from the brain of a rabid dog. The overall length of the genome was 11 924 nucleotides. Comparison of the genomic sequence showed the homology of SH06 at nucleotide level with full-length genomes of reference vaccine strains ranged from 82.2% with the PV strain to 86.9% with the CTN strain. A full-length genome-based phylogenetic analysis was performed with sequences available from GenBank. Phylogenetic analysis of the complete genome sequences indicated that the SH06 exhibited the highest homology with rabies street virus BD06 and CTN vaccine strain originated from China.     

The canine sarcoglycan delta gene: BAC clone contig assembly, chromosome assignment and interrogation as a candidate gene for dilated cardiomyopathy in Dobermann dogs

Dilated cardiomyopathy (DCM) is a common disease of the myocardium recognized in human, dog and experimental animals. Genetic factors are responsible for a large proportion of cases in humans, and 17 genes with DCM causing mutations have been identified. The genetic origin of DCM in the Dobermann dogs has been suggested, but no disease genes have been identified to date. In this paper, we describe the characterization and evaluation of the canine sarcoglycan delta (SGCD), a gene implicated in DCM in human and hamster. Bacterial artificial chromosomes (BACs) containing the canine SGCD gene were isolated with probes for exon 3 and exons 4-8 and were characterized by Southern blot analysis. BAC end sequences were obtained for four BACs. Three of the BACs overlapped and could be ordered relative to each other and the end sequences of all four BACs could be anchored on the preliminary assembly of the dog genome sequence (www. ensembl.org). One of the BACs of the partial contig was localized by fluorescent in situ hybridization to canine chromosome 4q22, in agreement with the dog genome sequence. Two highly informative polymorphic microsatellite markers in intron 7 of the SGCD gene were identified. In 25 DCM-affected and 13 non DCM-affected dogs seven different haplotypes could be distinguished. However, no association between any of the SGCD variants and the disease locus was apparent.     

Origin and evolution of a major feline satellite DNA

A major satellite DNA has been cloned from the domestic cat (Felis catus) and characterized. The satellite monomer, termed FA-SAT, is 483 base-pairs in size, 64% G + C, and represents about 1 to 2% of the cat genome. A consensus sequence based upon partial sequence data from 21 independently isolated clones demonstrates: (1) FA-SAT is not composed of a series of shorter repeats, although about 25 copies, primarily imperfect, of the hexanucleotide TAACCC appear in the sequence; (2) there are many more CpG dinucleotides present in FA-SAT than expected for a random sequence of its size; and (3) 61% of all base substitutions in FA-SAT involve the replacement of G and C residues by A and T residues, indicating that FA-SAT is rapidly becoming A + T-rich. FA-SAT-related sequences are found in many mammals, where they appear to be scattered throughout the genome and not tandemly arranged as in the cat. An FA-SAT-related sequence was cloned from the domestic dog genome and sequenced, and shown to contain multiple copies of the same TAACCC hexanucleotide found in the cat satellite.     

The dog genome map and its use in mammalian comparative genomics

The dog genome organization was extensively studied in the last ten years. The most important achievements are the well-developed marker genome maps, including over 3200 marker loci, and a survey of the DNA genome sequence. This knowledge, along with the most advanced map of the human genome, turned out to be very useful in comparative genomic studies. On the one hand, it has promoted the development of marker genome maps of other species of the family Canidae (red fox, arctic fox, Chinese raccoon dog) as well as studies on the evolution of their karyotype. But the most important approach is the comparative analysis of human and canine hereditary diseases. At present, causative gene mutations are known for 30 canine hereditary diseases. A majority of them have human counterparts with similar clinical and molecular features. Studies on identification of genes having a major impact on some multifactorial diseases (hip dysplasia, epilepsy) and cancers (multifocal renal cystadenocarcinoma and nodular dermatofibrosis) are advanced. Very promising are the results of gene therapy for certain canine monogenic diseases (haemophilia, hereditary retinal dystrophy, mucopolysaccharidosis), which have human equivalents. The above-mentioned examples prove a very important model role of the dog in studies of human genetic diseases. On the other hand, the identification of gene mutations responsible for hereditary diseases has a substantial impact on breeding strategy in the dog.     

Sleeping dogs

The dog genome sequence represents a landmark in the genomics era for a number of reasons.     

The complete mitochondrial genome and phylogenetic analysis of Nyctereutes procyonoides

The complete mitochondrial genome sequence of the raccoon dog (Nyctereutes procyonoides) was determined by using the long and accurate polymerase chain reaction. The entire mitochondrial genome sequence is 16,713 bp in length contains two ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA genes and 1 control region. Most mitochondrial genes are encoded on the H strand, except for the ND6 gene and 8 tRNA genes. The base compositions of mitochondrial genomes present clearly A–T skew. All the transfer RNA genes can be folded into the typical cloverleaf-shaped structure except tRNA-Ser (AGY), which lacks the dihydrouridine arm. Protein-coding genes mainly initiate with ATG and terminate with TAA. Some reading frame intervals and overlaps are found in the mitochondrial genome. The control region can be divided into three domains: the extended termination associated sequences (ETASs) domain, the central conserved domain and the conserved sequence blocks (CSBs) domain. Three conserved sequence blocks (CSBs) and one extended termination associated sequences (ETAS-1) is found in the control region. The phylogenetic analysis based on the concatenated data set of 14 genes in the mitochondrial genome of Canidae shows that the raccoon dog has close phylogenetic position with the red fox (Vulpes vulpes) and they constitute a clade which has an equil evolutionary position with the clade formed by the genera Canis and Cuon.     

Novel origins of copy number variation in the dog genome

ABSTRACT: BACKGROUND: Copy number variants (CNVs) account for substantial variation between genomes and are a major source of normal and pathogenic phenotypic differences. The dog is an ideal model to investigate mutational mechanisms that generate CNVs as its genome lacks a functional ortholog of the PRDM9 gene implicated in recombination and CNV formation in humans. Here we comprehensively assay CNVs using high-density array comparative genomic hybridization in 50 dogs from 17 dog breeds and 3 gray wolves. RESULTS: We use a stringent new method to identify a total of 430 high-confidence CNV loci, that range in size from 9 kb to 1.6 Mb and span 26.4 Mb, or 1.08%, of the assayed dog genome, overlapping 413 annotated genes. 98% of CNVs observed in each breed are also observed in multiple breeds. CNVs predicted to disrupt gene function are significantly less common than expected by chance. We identify a significant overrepresentation of peaks of GC content, previously shown to be enriched in dog recombination hotspots, in the vicinity of CNV breakpoints. CONCLUSIONS: A number of the CNVs identified by this study are candidates for generating breed-specific phenotypes. Purifying selection seems to be a major factor shaping structural variation in the dog genome, suggesting that many CNVs are deleterious. Localized peaks of GC content appear to be novel sites of CNV formation in the dog genome by non-allelic homologous recombination, potentially activated by the loss of PRDM9. These sequence features may have driven genome instability and chromosomal rearrangements throughout canid evolution.     

Lessons learned from the dog genome

Extensive genetic resources and a high-quality genome sequence position the dog as an important model species for understanding genome evolution, population genetics and genes underlying complex phenotypic traits. Newly developed genomic resources have expanded our understanding of canine evolutionary history and dog origins. Domestication involved genetic contributions from multiple populations of gray wolves probably through backcrossing. More recently, the advent of controlled breeding practices has segregated genetic variability into distinct dog breeds that possess specific phenotypic traits. Consequently, genome-wide association and selective sweep scans now allow the discovery of genes underlying breed-specific characteristics. The dog is finally emerging as a novel resource for studying the genetic basis of complex traits, including behavior.     

A cytogenetically characterized, genome-anchored 10-Mb BAC set and CGH array for the domestic dog

The generation of a 7.5x dog genome assembly provides exciting new opportunities to interpret tumor-associated chromosome aberrations at the biological level. We present a genomic microarray for array comparative genomic hybridization (aCGH) analysis in the dog, comprising 275 bacterial artificial chromosome (BAC) clones spaced at intervals of approximately 10 Mb. Each clone has been positioned accurately within the genome assembly and assigned to a unique chromosome location by fluorescence in situ hybridization (FISH) analysis, both individually and as chromosome-specific BAC pools. The microarray also contains clones representing the dog orthologues of 31 genes implicated in human cancers. FISH analysis of the 10-Mb BAC clone set indicated excellent coverage of each dog chromosome by the genome assembly. The order of clones was consistent with the assembly, but the cytogenetic intervals between clones were variable. We demonstrate the application of the BAC array for aCGH analysis to identify both whole and partial chromosome imbalances using a canine histiocytic sarcoma case. Using BAC clones selected from the array as probes, multicolor FISH analysis was used to further characterize these imbalances, revealing numerous structural chromosome rearrangements. We outline the value of a combined aCGH/FISH approach, together with a well-annotated dog genome assembly, in canine and comparative cancer studies.     

The dog: A powerful model for studying genotype–phenotype relationships

Within the last two years, series of studies have focused on the structure of the dog genome (Canis familiaris) and the characteristics of the dog population as it evolved since being domesticated from wolves about 14,000 years ago. In this review, we explain why the dog is a unique and promising model for determining genotype/phenotype relationships and why it should be easier with this model to identify the genes responsible for many genetic diseases. We also revisit the last ten years of developments in canine molecular genetics that culminated in the release of the entire genome sequence.     

貉源阿留申病毒全基因组测序及末端序列分子特征分析

貉源阿留申病毒(Raccoon dog and arctic fox amdoparvovirus,RFAV)是自然感染貉和蓝狐的新种阿留申病毒(Amdoparvovirus),为测序RFAV全基因组序列,预测分析RFAV末端发夹结构序列分子特征。本研究采用分段克隆成功获得3株长4832nt、4827nt、4830nt的RFAV全基因组序列,分别命名为RFAV-Y9J、RFAV-RD15、RFAV-HS-R,利用在线软件预测RFAV末端序列二级结构,并与水貂阿留申病毒(AMDV)末端序列进行同源性比对。结果显示阿留申病毒种间、种内3’末端基因组序列保守性强,均存在116nt的Y型发夹结构;RFAV-Y9J与RFAV-RD15毒株5′末端分别存在310nt、305nt的U型发夹结构,RFAV和AMDV种内5′末端基因组序列保守性强,而种间5′末端基因组序列有较大变异。本研究首次完整测序了RFAV的3′和5′末端序列,为其他种阿留申病毒的末端序列扩增提供一种有效方法,为构建RFAV的全基因组序列感染性克隆奠定了基础。     

Mapping and sequence of the gene for the pseudorabies virus glycoprotein which accumulates in the medium of infected cells.

RNA from pseudorabies virus (PRV)-infected cells was translated in a reticulocyte lysate with and without the addition of dog pancreas microsomes. Upon addition of the microsomes to the translation reaction, an additional prominent protein product was observed that was not present when microsomes were omitted. The gene coding for this processed protein and its lower-molecular-weight precursor was mapped within the small unique region of the genome by hybridization of mRNA to cloned fragments of PRV DNA and translation of the selected mRNAs. A fragment of the coding region of this gene was inserted into an open reading frame cloning vector to express part of this gene as a hybrid protein in Escherichia coli. This hybrid protein was injected into mice to raise an antiserum which was found to precipitate the glycoprotein which accumulates in the medium of PRV-infected cells. This allows us to conclude that the gene for the "excreted" glycoprotein (gX) maps to the small unique region of the genome, and that the precursor of this glycoprotein is readily processed by dog pancreas microsomes. The region of the PRV genome which codes for this glycoprotein was sequenced and found to include an open reading frame coding for 498 amino acids, flanked by sequences which contain features common to eucaryotic promoters and polyadenylation signals. The predicted protein sequence includes a hydrophobic sequence at the N-terminus which could be a signal sequence, and a hydrophobic sequence followed by a hydrophilic sequence at the C-terminus.     

Using comparative genomics to reorder the human genome sequence into a virtual sheep genome

Background

Is it possible to construct an accurate and detailed subgene-level map of a genome using bacterial artificial chromosome (BAC) end sequences, a sparse marker map, and the sequences of other genomes?

Results

A sheep BAC library, CHORI-243, was constructed and the BAC end sequences were determined and mapped with high sensitivity and low specificity onto the frameworks of the human, dog, and cow genomes. To maximize genome coverage, the coordinates of all BAC end sequence hits to the cow and dog genomes were also converted to the equivalent human genome coordinates. The 84,624 sheep BACs (about 5.4-fold genome coverage) with paired ends in the correct orientation (tail-to-tail) and spacing, combined with information from sheep BAC comparative genome contigs (CGCs) built separately on the dog and cow genomes, were used to construct 1,172 sheep BAC-CGCs, covering 91.2% of the human genome. Clustered non-tail-to-tail and outsize BACs located close to the ends of many BAC-CGCs linked BAC-CGCs covering about 70% of the genome to at least one other BAC-CGC on the same chromosome. Using the BAC-CGCs, the intrachromosomal and interchromosomal BAC-CGC linkage information, human/cow and vertebrate synteny, and the sheep marker map, a virtual sheep genome was constructed. To identify BACs potentially located in gaps between BAC-CGCs, an additional set of 55,668 sheep BACs were positioned on the sheep genome with lower confidence. A coordinate conversion process allowed us to transfer human genes and other genome features to the virtual sheep genome to display on a sheep genome browser.

Conclusion

We demonstrate that limited sequencing of BACs combined with positioning on a well assembled genome and integrating locations from other less well assembled genomes can yield extensive, detailed subgene-level maps of mammalian genomes, for which genomic resources are currently limited.     

The canine era: the rise of a biomedical model

Since the annotation of its genome a decade ago, the dog has proven to be an excellent model for the study of inherited diseases. A large variety of spontaneous simple and complex phenotypes occur in dogs, providing physiologically relevant models to corresponding human conditions. In addition, gene discovery is facilitated in clinically less heterogeneous purebred dogs with closed population structures because smaller study cohorts and fewer markers are often sufficient to expose causal variants. Here, we review the development of genomic resources from microsatellites to whole‐genome sequencing and give examples of successful findings that have followed the technological progress. The increasing amount of whole‐genome sequence data warrants better functional annotation of the canine genome to more effectively utilise this unique model to understand genetic contributions in morphological, behavioural and other complex traits.     

A comparative analysis of MC4R gene sequence, polymorphism, and chromosomal localization in Chinese raccoon dog and Arctic fox

Numerous mutations of the human melanocortin receptor type 4 (MC4R) gene are responsible for monogenic obesity, and some of them appear to be associated with predisposition or resistance to polygenic obesity. Thus, this gene is considered a functional candidate for fat tissue accumulation and body weight in domestic mammals. The aim of the study was comparative analysis of chromosome localization, nucleotide sequence, and polymorphism of the MC4R gene in two farmed species of the Canidae family, namely the Chinese raccoon dog (Nycterutes procyonoides procyonoides) and the arctic fox (Alopex lagopus). The whole coding sequence, including fragments of 3'UTR and 5'UTR, shows 89% similarity between the arctic fox (1276 bp) and Chinese raccoon dog (1213 bp). Altogether, 30 farmed Chinese raccoon dogs and 30 farmed arctic foxes were searched for polymorphisms. In the Chinese raccoon dog, only one silent substitution in the coding sequence was identified; whereas in the arctic fox, four InDels and two single-nucleotide polymorphisms (SNPs) in the 5'UTR and six silent SNPs in the exon were found. The studied gene was mapped by FISH to the Chinese raccoon dog chromosome 9 (NPP9q1.2) and arctic fox chromosome 24 (ALA24q1.2-1.3). The obtained results are discussed in terms of genome evolution of species belonging to the family Canidae and their potential use in animal breeding.     

应用生物信息学寻找山羊新的microRNA分子及其实验验证

摘要: microRNA(miRNA)是一类长约22个碱基的非编码RNA分子, 在转录后水平调节基因的表达及其在细胞的增殖、分化、凋亡等过程中起着重要的调控作用。根据miRNA分子具有一定的保守性, 文章将人、小鼠、牛、猪和狗5种哺乳动物已知的miRNA分子与NCBI公布的与山羊具有极高同源性的绵羊基因组序列对比, 获得11条miRNA候选分子, 然后通过逆转录聚合酶链反应(RT-PCR)验证, 发现在山羊脑组织中这11条分子均有表达, 肝脏组织中有5条分子表达, 初步确定为山羊新的miRNA分子, 为寻找山羊miRNA提供了新的 思路。     

The canine olfactory subgenome

We identified 971 olfactory receptor (OR) genes in the dog genome, estimated to constitute approximately 80% of the canine OR repertoire. This was achieved by directed genomic DNA cloning of olfactory sequence tags as well as by mining the Celera canine genome sequences. The dog OR subgenome is estimated to have 12% pseudogenes, suggesting a functional repertoire similar to that of mouse and considerably larger than for humans. No novel OR families were discovered, but as many as 34 gene subfamilies were unique to the dog. "Fish-like" Class I ancient ORs constituted 18% of the repertoire, significantly more than in human and mouse. A set of 122 dog-human-mouse ortholog triplets was identified, with a relatively high fraction of Class I ORs. The elucidation of a large portion of the canine olfactory receptor gene superfamily, with some dog-specific attributes, may help us understand the unique chemosensory capacities of this species.