Y-chromosome microsatellite variation in Italian wolves: A contribution to the study of wolf-dog hybridization patterns

One major concern in wolf (Canis lupus) conservation is the risk of genetic contamination due to crossbreeding with domestic dogs. Although genetic monitoring of wolf populations has become widely used, the behavioural mechanisms involved in wolf-dog hybridization and the detrimental effects of genetic introgression are poorly known. In this study we analysed Y-chromosome microsatellite variation in the recovering Italian wolf population and detected strikingly different allele frequencies between wolves and dogs. Four Y haplotypes were found in 74 analysed male wolves, and all of them were present in a focus wolf population in the Apennines. On the other hand, only 1 haplotype was found in the recolonizing wolf population from the Western Alps. The most common haplotype in a sample of domestic dogs, was also found in 5 wolves, 2 of which revealing a signature of recent hybridization. Moreover, another suspect hybrid carried a private haplotype of possible canine origin. These results give support to the idea that female wolves can breed with male stray dogs in the wild. The Y-chromosome variation in Italian wolves contrasts with the previously observed lack of mitochondrial variation. Further investigations are needed to clarify at what extent historical or recent wolf-dog hybridization events may have contributed to the observed haplotype diversity. In conclusion, the two molecular markers employed in this study represent effective means to trace directional genetic introgression into the wolves male lineage and have the noteworthy advantage of being suitable for analyses on low-quality DNA samples.     

Characterization of 59 canine single nucleotide polymorphisms in the Italian wolf (Canis lupus) population

We characterized 59 canine single nucleotide polymorphisms (SNPs) in the endangered Italian wolf (Canis lupus) population, which were discovered by resequencing sequence‐tagged‐site (STS) DNA sequences that are known to contain SNPs in domestic dogs. Dog SNPs were usually found also in wolves. Additional SNPs unique in dogs or wolves were discovered, which is important for detecting hybrids between dogs and wolves. We developed new primer sets and analysed 15 SNPs by Pyrosequencing. The characterized SNPs will provide an important addition to the genetic markers that are currently available for studying wild populations of canids.     

Sequence analysis of three pigmentation genes in the Newfoundland population of Canis latrans links the Golden Retriever Mc1r variant to white coat color in coyotes

Three genes, Mc1r, Agouti, and CBD103, interact in a type-switching process that controls much of the pigmentation variation observed in mammals. A deletion in the CBD103 gene is responsible for dominant black color in dogs, while the white-phased black bear (“spirit bear”) of British Columbia, Canada, is the lightest documented color variant caused by a mutation in Mc1r. Rare all-white animals have recently been discovered in a new northeastern population of the coyote in insular Newfoundland and Labrador, Canada. To investigate the causative gene and mutation of white coat in coyotes, we sequenced the three type-switching genes in white and dark-phased animals from Newfoundland. The only sequence variants unambiguously associated with white color were in Mc1r, and one of these variants causes the amino acid variant R306Ter, a premature stop codon also linked to coat color in Golden Retrievers and other dogs with yellow/red coats. The allele carrying R306Ter in coyotes matches that in the Golden Retriever at other variable amino acid sites and hence may have originated in these dogs. Coyotes experienced introgression with wolves and dogs as they colonized northeastern North America, and coyote/Golden Retriever interactions have been observed in Newfoundland. We speculate that natural selection, with or without a founder effect, may contribute to the observed frequency of white coyotes in Newfoundland, as it has contributed to the high frequency of white bears, and of a domestic dog-derived CBD allele in gray wolves.     

A second KRT71 allele in curly coated dogs

Major characteristics of coat variation in dogs can be explained by variants in only a few genes. Until now, only one missense variant in the KRT71 gene, p.Arg151Trp, has been reported to cause curly hair in dogs. However, this variant does not explain the curly coat in all breeds as the mutant ¹⁵¹Trp allele, for example, is absent in Curly Coated Retrievers. We sequenced the genome of a Curly Coated Retriever at 22× coverage and searched for variants in the KRT71 gene. Only one protein‐changing variant was present in a homozygous state in the Curly Coated Retriever and absent or present in a heterozygous state in 221 control dogs from different dog breeds. This variant, NM_001197029.1:c.1266_1273delinsACA, was an indel variant in exon 7 that caused a frameshift and an altered and probably extended C‐terminus of the KRT71 protein NP_001183958.1:p.(Ser422ArgfsTer?). Using Sanger sequencing, we found that the variant was fixed in a cohort of 125 Curly Coated Retrievers and segregating in five of 14 additionally tested breeds with a curly or wavy coat. KRT71 variants cause curly hair in humans, mice, rats, cats and dogs. Specific KRT71 variants were further shown to cause alopecia. Based on this knowledge from other species and the predicted molecular consequence of the newly identified canine KRT71 variant, it is a compelling candidate causing a second curly hair allele in dogs. It might cause a slightly different coat phenotype than the previously published p.Arg151Trp variant and could potentially be associated with follicular dysplasia in dogs.     

Dietary Variation and Evolution of Gene Copy Number among Dog Breeds

Prolonged human interactions and artificial selection have influenced the genotypic and phenotypic diversity among dog breeds. Because humans and dogs occupy diverse habitats, ecological contexts have likely contributed to breed-specific positive selection. Prior to the advent of modern dog-feeding practices, there was likely substantial variation in dietary landscapes among disparate dog breeds. As such, we investigated one type of genetic variant, copy number variation, in three metabolic genes: glucokinase regulatory protein (GCKR), phytanol-CoA 2-hydroxylase (PHYH), and pancreatic α-amylase 2B (AMY2B). These genes code for proteins that are responsible for metabolizing dietary products that originate from distinctly different food types: sugar, meat, and starch, respectively. After surveying copy number variation among dogs with diverse dietary histories, we found no correlation between diet and positive selection in either GCKR or PHYH. Although it has been previously demonstrated that dogs experienced a copy number increase in AMY2B relative to wolves during or after the dog domestication process, we demonstrate that positive selection continued to act on amylase copy number in dog breeds that consumed starch-rich diets in time periods after domestication. Furthermore, we found that introgression with wolves is not responsible for deterioration of positive selection on AMY2B among diverse dog breeds. Together, this supports the hypothesis that the amylase copy number expansion is found universally in dogs.     

Patterns of Exposure of Iberian Wolves (<Emphasis Type="Italic">Canis lupus</Emphasis>) to Canine Viruses in Human-Dominated Landscapes

Wildlife inhabiting human-dominated landscapes is at risk of pathogen spill-over from domestic species. With the aim of gaining knowledge in the dynamics of viral infections in Iberian wolves (Canis lupus) living in anthropized landscapes of northern Spain, we analysed between 2010 and 2013 the samples of 54 wolves by serology and polymerase chain reaction (PCR) for exposure to four pathogenic canine viruses: canine distemper virus (CDV), canine parvovirus-2 (CPV), canine adenovirus 1 and 2 (CAV-1 and CAV-2) and canine herpesvirus. Overall, 76% of the studied wolves presented evidence of exposure to CPV (96% by HI, 66% by PCR) and 75% to CAV (75% by virus neutralization (VN), 76% by PCR, of which 70% CAV-1 and 6% CAV-2). This represents the first detection of CAV-2 infection in a wild carnivore. CPV/CAV-1 co-infection occurred in 51% of the wolves. The probability of wolf exposure to CPV was positively and significantly correlated with farm density in a buffer zone around the place where the wolf was found, indicating that rural dogs might be the origin of CPV infecting wolves. CPV and CAV-1 appear to be enzootic in the Iberian wolf population, which is supported by the absence of seasonal and inter-annual variations in the proportion of positive samples detected. However, while CPV may depend on periodical introductions by dogs, CAV-1 may be maintained within the wolf population. All wolves were negative for exposure to CDV (by VN and PCR) and CHV (by PCR). The absence of acquired immunity against CDV in this population may predispose it to an elevated rate of mortality in the event of a distemper spill-over via dogs.     

A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves

The domestic dog serves as an excellent model to investigate the genetic basis of disease. More than 400 heritable traits analogous to human diseases have been described in dogs. To further canine medical genetics research, we established the Dog Biomedical Variant Database Consortium (DBVDC) and present a comprehensive list of functionally annotated genome variants that were identified with whole genome sequencing of 582 dogs from 126 breeds and eight wolves. The genomes used in the study have a minimum coverage of 10× and an average coverage of ~24×. In total, we identified 23 133 692 single‐nucleotide variants (SNVs) and 10 048 038 short indels, including 93% undescribed variants. On average, each individual dog genome carried ～4.1 million single‐nucleotide and ~1.4 million short‐indel variants with respect to the reference genome assembly. About 2% of the variants were located in coding regions of annotated genes and loci. Variant effect classification showed 247 141 SNVs and 99 562 short indels having moderate or high impact on 11 267 protein‐coding genes. On average, each genome contained heterozygous loss‐of‐function variants in 30 potentially embryonic lethal genes and 97 genes associated with developmental disorders. More than 50 inherited disorders and traits have been unravelled using the DBVDC variant catalogue, enabling genetic testing for breeding and diagnostics. This resource of annotated variants and their corresponding genotype frequencies constitutes a highly useful tool for the identification of potential variants causative for rare inherited disorders in dogs.     

Long-distance dispersal of a wolf, <Emphasis Type="Italic">Canis lupus</Emphasis>, in northwestern Europe

Several mammal species have recolonized their historical ranges across Europe during the last decades. In November 2012, a wolf-looking canid was found dead in Thy National Park (56° 56′ N, 8° 25′ E) in Jutland, Denmark. DNA from this individual and nine German wolves were genotyped using a genome-wide panel of 22,163 canine single nucleotide polymorphism (SNP) markers and compared to existing profiles based on the same marker panel obtained from northeastern Polish (n?=?13) wolves, domestic dogs (n?=?13) and known wolf-dog hybrids (n?=?4). The Thy canid was confirmed to be a wolf from the German-western Polish population, approximately 800 km to the southeast. Access to the German reference database on DNA profiles based on 13 autosomal microsatellites of German wolves made it possible to pinpoint the exact pack origin of the Thy wolf in Saxony, Germany. This was the first documented observation of a wolf in Denmark in 200 years and another example of long-distance dispersal of a carnivore.     

Detection of an East European wolf haplotype puzzles mitochondrial DNA monomorphism of the Italian wolf population

Southern European wolves suffered from reiterated population declines during glacial periods and historically due to human persecution. Differently from other European wolf populations, a single mitochondrial DNA (mtDNA) control region haplotype (W14) has been so far described in the Italian wolves, although no intensive genetic sampling has ever been conducted in historical source populations from central and southern Italy. Using non-invasive genetic techniques, we report the occurrence of an unexpected mtDNA haplotype (W16) in the wolf population of the Abruzzo, Lazio and Molise National Park (PNALM), central Italy. This haplotype, detected in three out of 90 faecal samples from the PNALM, was previously reported in wolves from the North Carpathians, Slovakia and the Balkans only. Microsatellite analysis and molecular sex determination confirmed that the W16 samples belonged to three distinct wolves. Although alternative explanations can be formulated for the origin of this mtDNA haplotype in the otherwise monomorphic Italian wolf population, assignment procedures indicated the likely admixed ancestry of one W16 sample with East European wolves. Anthropogenic introgression with dogs has been detected in the Italian wolf population using nuclear DNA microsatellites, but no population-wide genetic survey had previously reported a mtDNA control region variant in Italian wolves. Our findings strongly suggest that, in addition to wolf × dog hybridization, captive-released wolves or wolf × dog hybrids may successfully interbreed with wolves in the wild, and that human-mediated introgression may occur even in well established protected areas.     

Genome Wide Analysis Indicates Genes for Basement Membrane and Cartilage Matrix Proteins as Candidates for Hip Dysplasia in Labrador Retrievers

Hip dysplasia, an abnormal laxity of the hip joint, is seen in humans as well as dogs and is one of the most common skeletal disorders in dogs. Canine hip dysplasia is considered multifactorial and polygenic, and a variety of chromosomal regions have been associated with the disorder. We performed a genome-wide association study in Dutch Labrador Retrievers, comparing data of nearly 18,000 single nucleotide polymorphisms (SNPs) in 48 cases and 30 controls using two different statistical methods. An individual SNP analysis based on comparison of allele frequencies with a χ² statistic was used, as well as a simultaneous SNP analysis based on Bayesian variable selection. Significant association with canine hip dysplasia was observed on chromosome 8, as well as suggestive association on chromosomes 1, 5, 15, 20, 25 and 32. Next-generation DNA sequencing of the exons of genes of seven regions identified multiple associated alleles on chromosome 1, 5, 8, 20, 25 and 32 (p<0.001). Candidate genes located in the associated regions on chromosomes 1, 8 and 25 included LAMA2, LRR1 and COL6A3, respectively. The associated region on CFA20 contained candidate genes GDF15, COMP and CILP2. In conclusion, our study identified candidate genes that might affect susceptibility to canine hip dysplasia. These genes are involved in hypertrophic differentiation of chondrocytes and extracellular matrix integrity of basement membrane and cartilage. The functions of the genes are in agreement with the notion that disruptions in endochondral bone formation in combination with soft tissue defects are involved in the etiology of hip dysplasia.     

Association of a P2RX7 gene missense variant with brachycephalic dog breeds

Missense variants are associated with various phenotypic traits and disorders in dogs. The canine P2RX7 gene, coding the ATP-gated P2X7 receptor ion channel, contains four known missense variants. The current study aimed to examine the presence of these variants in a random sample of pedigree and mixed-pedigree dogs. Exons 3, 8, 11 and 13 of the P2RX7 gene, encoding these four respective variants, in 65 dogs were assessed by Sanger sequencing and combined with existing sequencing data from another 69 dogs. The distribution of these variants was then evaluated in all 134 dogs combined and separately within individual breeds including 35 different pure breeds. The rs23314713 (p.Phe103Leu) and rs23315462 (p.Pro452Ser) variants were present in 47 and 40% of all dogs studied respectively, with the rs23314713 variant associated with brachycephalic breeds. Among pedigree dogs, the rs23314713 and rs23315462 variants were associated with brachycephalic and non-brachycephalic breeds respectively. The rs851148233 (p.Arg270Cys) and rs850760787 (p.Arg365Gln) variants were present only in dogs of Cocker Spaniel and Labrador Retriever pedigrees respectively. No other missense variants were found in exons 3, 8, 11 and 13 of the P2RX7 gene within the dogs. In conclusion, the rs23314713 and rs23315462 missense variants of the P2RX7 gene are present in a large proportion of dogs, with the rs23314713 variant associated with a number of brachycephalic breeds. However, the association of this variant with dogs of bulldog ancestry, not brachycephaly per se, cannot be excluded.     

Analysis of structural diversity in wolf-like canids reveals post-domestication variants

Background

Although a variety of genetic changes have been implicated in causing phenotypic differences among dogs, the role of copy number variants (CNVs) and their impact on phenotypic variation is still poorly understood. Further, very limited knowledge exists on structural variation in the gray wolf, the ancestor of the dog, or other closely related wild canids. Documenting CNVs variation in wild canids is essential to identify ancestral states and variation that may have appeared after domestication.

Results

In this work, we genotyped 1,611 dog CNVs in 23 wolf-like canids (4 purebred dogs, one dingo, 15 gray wolves, one red wolf, one coyote and one golden jackal) to identify CNVs that may have arisen after domestication. We have found an increase in GC-rich regions close to the breakpoints and around 1 kb away from them suggesting that some common motifs might be associated with the formation of CNVs. Among the CNV regions that showed the largest differentiation between dogs and wild canids we found 12 genes, nine of which are related to two known functions associated with dog domestication; growth (PDE4D, CRTC3 and NEB) and neurological function (PDE4D, EML5, ZNF500, SLC6A11, ELAVL2, RGS7 and CTSB).

Conclusions

Our results provide insight into the evolution of structural variation in canines, where recombination is not regulated by PRDM9 due to the inactivation of this gene. We also identified genes within the most differentiated CNV regions between dogs and wolves, which could reflect selection during the domestication process.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-465) contains supplementary material, which is available to authorized users.     

Molecular Characterization of Three Canine Models of Human Rare Bone Diseases: Caffey,van den Ende-Gupta,and Raine Syndromes

One to two percent of all children are born with a developmental disorder requiring pediatric hospital admissions. For many such syndromes, the molecular pathogenesis remains poorly characterized. Parallel developmental disorders in other species could provide complementary models for human rare diseases by uncovering new candidate genes, improving the understanding of the molecular mechanisms and opening possibilities for therapeutic trials. We performed various experiments, e.g. combined genome-wide association and next generation sequencing, to investigate the clinico-pathological features and genetic causes of three developmental syndromes in dogs, including craniomandibular osteopathy (CMO), a previously undescribed skeletal syndrome, and dental hypomineralization, for which we identified pathogenic variants in the canine SLC37A2 (truncating splicing enhancer variant), SCARF2 (truncating 2-bp deletion) and FAM20C (missense variant) genes, respectively. CMO is a clinical equivalent to an infantile cortical hyperostosis (Caffey disease), for which SLC37A2 is a new candidate gene. SLC37A2 is a poorly characterized member of a glucose-phosphate transporter family without previous disease associations. It is expressed in many tissues, including cells of the macrophage lineage, e.g. osteoclasts, and suggests a disease mechanism, in which an impaired glucose homeostasis in osteoclasts compromises their function in the developing bone, leading to hyperostosis. Mutations in SCARF2 and FAM20C have been associated with the human van den Ende-Gupta and Raine syndromes that include numerous features similar to the affected dogs. Given the growing interest in the molecular characterization and treatment of human rare diseases, our study presents three novel physiologically relevant models for further research and therapy approaches, while providing the molecular identity for the canine conditions.     

Single Nucleotide Variant rs2232710 in the Protein Z-Dependent Protease Inhibitor (ZPI,SERPINA10) Gene Is Not Associated with Deep Vein Thrombosis

Rare mutations in PROC, PROS1 or SERPINC1 as well as common variants in F5, F2, F11 and SERPINC1 have been identified as risk factors for deep vein thrombosis (DVT). To identify novel genetic risk factors for DVT, we have developed and applied next-generation DNA sequencing (NGS) of the coding area of hemostatic and proinflammatory genes. Using this strategy, we previously identified a single nucleotide variant (SNV) rs6050 in the FGA gene and novel, rare SNVs in the ADAMTS13 gene associated with DVT. To identify novel coding variants in the genetic predisposition to DVT, we applied NGS analysis of the coding area of 186 hemostatic and proinflammatory genes in 94 DVT cases and 98 controls and we identified 18 variants with putative role in DVT. A group of 585 Italian idiopathic DVT patients and 550 healthy controls was used to genotype all the 18 risk-associated variants identified by NGS. Replication study in the Italian population identified the rs2232710 variant in the protein Z-dependent protease inhibitor (ZPI) gene to be associated with an increased risk of DVT (OR 2.74; 95% CI 1.33–5.65; P = 0.0045; Bonferroni P = 0.081). However, the rs2232710 SNV showed no association with DVT in two Dutch replication cohorts the LETS study (454 patients and 451 controls) and the MEGA study (3799 patients and 4399 controls), indicating that the rs2232710 variant is not a risk factor for DVT.     

Association between the g.296596G > A genetic variant of RELN gene and susceptibility to autism in a Chinese Han population

Autism is a childhood neuro-developmental disorder, and Reelin (RELN) is an important candidate gene for influencing autism. This study aimed at investigating the influence of genetic variants of the RELN gene on autism susceptibility. In this study, 205 autism patients and 210 healthy controls were recruited and the genetic variants of the RELN gene were genotyped by the created restriction site-polymerase chain reaction (CRS-PCR) method. The influence of genetic variants on autism susceptibility was analyzed by association analysis, and the g.296596G > A genetic variant in exon10 of the RELN gene was detected. The frequencies of allele/genotype in autistic patients were significantly different from those in healthy controls, and a statistically significant association was detected between this genetic variant and autism susceptibility. Our data lead to the inference that the g.296596G > A genetic variant in the RELN gene has a potential influence on autism susceptibility in the Chinese Han population.     

Genetic Mapping of Novel Loci Affecting Canine Blood Phenotypes

Since the publication of the dog genome and the construction of high-quality genome-wide SNP arrays, thousands of dogs have been genotyped for disease studies. For many of these dogs, additional clinical phenotypes are available, such as hematological and clinical chemistry results collected during routine veterinary care. Little is known about the genetic basis of variation in blood phenotypes, but this variation may play an important role in the etiology and progression of many diseases. From a cohort of dogs that had been previously genotyped on a semi-custom Illumina CanineHD array for various genome-wide association studies (GWAS) at Cornell University Hospital for Animals, we chose 353 clinically healthy, adult dogs for our analysis of clinical pathologic test results (14 hematological tests and 25 clinical chemistry tests). After correcting for age, body weight and sex, genetic associations were identified for amylase, segmented neutrophils, urea nitrogen, glucose, and mean corpuscular hemoglobin. Additionally, a strong genetic association (P = 8.1×10⁻¹³) was evident between a region of canine chromosome 13 (CFA13) and alanine aminotransferase (ALT), explaining 23% of the variation in ALT levels. This region of CFA13 encompasses the GPT gene that encodes the transferase. Dogs homozygous for the derived allele exhibit lower ALT activity, making increased ALT activity a less useful marker of hepatic injury in these individuals. Overall, these associations provide a roadmap for identifying causal variants that could improve interpretation of clinical blood tests and understanding of genetic risk factors associated with diseases such as canine diabetes and anemia, and demonstrate the utility of holistic phenotyping of dogs genotyped for disease mapping studies.     

A Deletion in the VLDLR Gene in Eurasier Dogs with Cerebellar Hypoplasia Resembling a Dandy-Walker-Like Malformation (DWLM)

Dandy-Walker-like malformation (DWLM) is the result of aberrant brain development and mainly characterized by cerebellar hypoplasia. DWLM affected dogs display a non-progressive cerebellar ataxia. Several DWLM cases were recently observed in the Eurasier dog breed, which strongly suggested a monogenic autosomal recessive inheritance in this breed. We performed a genome-wide association study (GWAS) with 9 cases and 11 controls and found the best association of DWLM with markers on chromosome 1. Subsequent homozygosity mapping confirmed that all 9 cases were homozygous for a shared haplotype in this region, which delineated a critical interval of 3.35 Mb. We sequenced the genome of an affected Eurasier and compared it with the Boxer reference genome and 47 control genomes of dogs from other breeds. This analysis revealed 4 private non-synonymous variants in the critical interval of the affected Eurasier. We genotyped these variants in additional dogs and found perfect association for only one of these variants, a single base deletion in the VLDLR gene encoding the very low density lipoprotein receptor. This variant, VLDLR:c.1713delC is predicted to cause a frameshift and premature stop codon (p.W572Gfs*10). Variants in the VLDLR gene have been shown to cause congenital cerebellar ataxia and mental retardation in human patients and Vldlr knockout mice also display an ataxia phenotype. Our combined genetic data together with the functional knowledge on the VLDLR gene from other species thus strongly suggest that VLDLR:c.1713delC is indeed causing DWLM in Eurasier dogs.