Recent, full-length gene retrocopies are common in canids.

Gene retrocopies arise from the reverse transcription and insertion into the genome of processed mRNA transcripts. Although many retrocopies have acquired mutations that render them functionally inactive, most mammals retain active LINE-1 sequences capable of producing new retrocopies. New retrocopies, referred to as retro copy number variants (retroCNVs), may not be identified by standard variant calling techniques in high-throughput sequencing data. Although multiple functional FGF4 retroCNVs have been associated with skeletal dysplasias in dogs, the full landscape of canid retroCNVs has not been characterized. Here, retroCNV discovery was performed on a whole-genome sequencing data set of 293 canids from 76 breeds. We identified retroCNV parent genes via the presence of mRNA-specific 30-mers, and then identified retroCNV insertion sites through discordant read analysis. In total, we resolved insertion sites for 1911 retroCNVs from 1179 parent genes, 1236 of which appeared identical to their parent genes. Dogs had on average 54.1 total retroCNVs and 1.4 private retroCNVs. We found evidence of expression in testes for 12% (14/113) of the retroCNVs identified in six Golden Retrievers, including four chimeric transcripts, and 97 retroCNVs also had significantly elevated F ST across dog breeds, possibly indicating selection. We applied our approach to a subset of human genomes and detected an average of 4.2 retroCNVs per sample, highlighting a 13-fold relative increase of retroCNV frequency in dogs. Particularly in canids, retroCNVs are a largely unexplored source of genetic variation that can contribute to genome plasticity and that should be considered when investigating traits and diseases.

Genetic Variants Affecting Skeletal Morphology in Domestic Dogs.

Purebred dog breeds provide a powerful resource for the discovery of genetic variants affecting skeletal morphology. Domesticated and subsequently purebred dogs have undergone strong artificial selection for a broad range of skeletal variation, which include both the size and shapes of their bones. While the phenotypic variation between breeds is high, within-breed morphological variation is typically low. Approaches for defining genetic variants associated with canine morphology include quantitative within-breed analyses, as well as across-breed analyses, using breed standards as proxies for individual measurements. The ability to identify variants across the genomes of individual dogs can now be paired with precise measures of morphological variation to define the genetic interactions and the phenotypic effect of variants on skeletal morphology.

NME5 frameshift variant in Alaskan Malamutes with primary ciliary dyskinesia.

Primary ciliary dyskinesia (PCD) is a hereditary defect of motile cilia in humans and several domestic animal species. Typical clinical findings are chronic recurrent infections of the respiratory tract and fertility problems. We analyzed an Alaskan Malamute family, in which two out of six puppies were affected by PCD. The parents were unaffected suggesting autosomal recessive inheritance. Linkage and homozygosity mapping defined critical intervals comprising ~118 Mb. Whole genome sequencing of one case and comparison to 601 control genomes identified a disease associated frameshift variant, c.43delA, in the NME5 gene encoding a sparsely characterized protein associated with ciliary function. Nme5-/- knockout mice exhibit doming of the skull, hydrocephalus and sperm flagellar defects. The genotypes at NME5:c.43delA showed the expected co-segregation with the phenotype in the Alaskan Malamute family. An additional unrelated Alaskan Malamute with PCD and hydrocephalus that became available later in the study was also homozygous mutant at the NME5:c.43delA variant. The mutant allele was not present in more than 1000 control dogs from different breeds. Immunohistochemistry demonstrated absence of the NME5 protein from nasal epithelia of an affected dog. We therefore propose NME5:c.43delA as the most likely candidate causative variant for PCD in Alaskan Malamutes. These findings enable genetic testing to avoid the unintentional breeding of affected dogs in the future. Furthermore, the results of this study identify NME5 as a novel candidate gene for unsolved human PCD and/or hydrocephalus cases.

Canine DVL2 variant contributes to brachycephalic phenotype and caudal vertebral anomalies.

A frameshift deletion variant in the Wnt pathway gene dishevelled 2 (DVL2) is associated with a truncated, kinked tail ("screw tail") in English Bulldogs, French Bulldogs and Boston Terriers. These breeds are also characterized by distinctive morphological traits, including a wide head, flat face and short-limbed dwarfism, which are characteristic of Robinow syndrome in humans, caused by defects in genes such as DVL1 and DVL3. Based on these phenotypic and genetic similarities, it has previously been hypothesized that the canine DVL2 variant results in a syndromic phenotype called the Robinow-like syndrome. In our study, we investigated the distribution of the DVL2 variant in 1954 dogs from 15 breeds, identifying breeds with allele variation and enabling the dissection of the genotype-phenotype correlation for the first time. With CT examinations in American Staffordshire Terriers, we confirmed that the DVL2 allele is associated with caudal vertebral malformations and a brachycephalic phenotype. We also hypothesize that the variant may be linked to additional health conditions, including brachycephalic obstructive airway syndrome and congenital heart defects. Altogether, our study strengthens the role of DVL2 as one of the contributors to the "bulldog type" morphology and features on the spectrum of human Robinow syndrome.

Variants at the ASIP locus contribute to coat color darkening in Nellore cattle.

BACKGROUND: Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls. RESULTS: We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons. CONCLUSIONS: Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin.

Whole genome variant association across 100 dogs identifies a frame shift mutation in DISHEVELLED 2 which contributes to Robinow-like syndrome in Bulldogs and related screw tail dog breeds.

Domestic dog breeds exhibit remarkable morphological variations that result from centuries of artificial selection and breeding. Identifying the genetic changes that contribute to these variations could provide critical insights into the molecular basis of tissue and organismal morphogenesis. Bulldogs, French Bulldogs and Boston Terriers share many morphological and disease-predisposition traits, including brachycephalic skull morphology, widely set eyes and short stature. Unlike other brachycephalic dogs, these breeds also exhibit vertebral malformations that result in a truncated, kinked tail (screw tail). Whole genome sequencing of 100 dogs from 21 breeds identified 12.4 million bi-allelic variants that met inclusion criteria. Whole Genome Association of these variants with the breed defining phenotype of screw tail was performed using 10 cases and 84 controls and identified a frameshift mutation in the WNT pathway gene DISHEVELLED 2 (DVL2) (Chr5: 32195043_32195044del, p = 4.37 X 10-37) as the most strongly associated variant in the canine genome. This DVL2 variant was fixed in Bulldogs and French Bulldogs and had a high allele frequency (0.94) in Boston Terriers. The DVL2 variant segregated with thoracic and caudal vertebral column malformations in a recessive manner with incomplete and variable penetrance for thoracic vertebral malformations between different breeds. Importantly, analogous frameshift mutations in the human DVL1 and DVL3 genes cause Robinow syndrome, a congenital disorder characterized by similar craniofacial, limb and vertebral malformations. Analysis of the canine DVL2 variant protein showed that its ability to undergo WNT-induced phosphorylation is reduced, suggesting that altered WNT signaling may contribute to the Robinow-like syndrome in the screwtail breeds.

FGF4 retrogene on CFA12 is responsible for chondrodystrophy and intervertebral disc disease in dogs.

Chondrodystrophy in dogs is defined by dysplastic, shortened long bones and premature degeneration and calcification of intervertebral discs. Independent genome-wide association analyses for skeletal dysplasia (short limbs) within a single breed (PBonferroni = 0.01) and intervertebral disc disease (IVDD) across breeds (PBonferroni = 4.0 × 10-10) both identified a significant association to the same region on CFA12. Whole genome sequencing identified a highly expressed FGF4 retrogene within this shared region. The FGF4 retrogene segregated with limb length and had an odds ratio of 51.23 (95% CI = 46.69, 56.20) for IVDD. Long bone length in dogs is a unique example of multiple disease-causing retrocopies of the same parental gene in a mammalian species. FGF signaling abnormalities have been associated with skeletal dysplasia in humans, and our findings present opportunities for both selective elimination of a medically and financially devastating disease in dogs and further understanding of the ever-growing complexity of retrogene biology.

Pathologic Features of the Intervertebral Disc in Young Nova Scotia Duck Tolling Retrievers Confirms Chondrodystrophy Degenerative Phenotype Associated With Genotype.

Chondrodystrophy results in predictable and progressive biochemical and structural changes to the intervertebral disc, resulting in early onset degeneration and dystrophic mineralization of the disc. Accelerated degeneration and mineralization of the intervertebral disc are common in multiple dog breeds and can result in compromised function, herniation, pain, and a variety of neurological sequelae. A mutation responsible for chondrodystrophy in dogs has been identified as an aberrant fibroblast growth factor 4 (FGF4) retrogene insertion on chromosome 12 (CFA12) and is associated with short stature of the Nova Scotia Duck Tolling Retriever. Segregation of the CFA12 FGF4 retrogene in this dog breed provides an opportunity to examine the effect of retrogene presence on radiographic and histologic appearance of chondrodystrophic disc degeneration within a single breed. Here we found that in the intervertebral discs isolated from 2 dogs with the CFA12 FGF4 genotype, the nucleus pulposus was largely replaced by cartilaginous tissue, and physaliferous notochordal cells were rarely if ever identified. These findings are in contrast to the normal histologic findings in 2 breed-matched dogs lacking the mutation. The findings are consistent with premature chondroid degeneration of the intervertebral disc and suggest that the presence of the CFA12 FGF4 retrogene is sufficient to cause the chondrodystrophic phenotype.

Utilizing the Dog Genome in the Search for Novel Candidate Genes Involved in Glioma Development-Genome Wide Association Mapping followed by Targeted Massive Parallel Sequencing Identifies a Strongly Associated Locus.

Gliomas are the most common form of malignant primary brain tumors in humans and second most common in dogs, occurring with similar frequencies in both species. Dogs are valuable spontaneous models of human complex diseases including cancers and may provide insight into disease susceptibility and oncogenesis. Several brachycephalic breeds such as Boxer, Bulldog and Boston Terrier have an elevated risk of developing glioma, but others, including Pug and Pekingese, are not at higher risk. To identify glioma-associated genetic susceptibility factors, an across-breed genome-wide association study (GWAS) was performed on 39 dog glioma cases and 141 controls from 25 dog breeds, identifying a genome-wide significant locus on canine chromosome (CFA) 26 (p = 2.8 x 10-8). Targeted re-sequencing of the 3.4 Mb candidate region was performed, followed by genotyping of the 56 SNVs that best fit the association pattern between the re-sequenced cases and controls. We identified three candidate genes that were highly associated with glioma susceptibility: CAMKK2, P2RX7 and DENR. CAMKK2 showed reduced expression in both canine and human brain tumors, and a non-synonymous variant in P2RX7, previously demonstrated to have a 50% decrease in receptor function, was also associated with disease. Thus, one or more of these genes appear to affect glioma susceptibility.

Evaluation of the major histocompatibility complex (MHC) class II as a candidate for sudden acquired retinal degeneration syndrome (SARDS) in Dachshunds.

OBJECTIVE: Sudden acquired retinal degeneration syndrome (SARDS) is one of the leading causes of acute blindness in dogs, with an unknown etiology and no effective treatment. Certain breeds such as Dachshunds are overrepresented among SARDS patients, and therefore, the syndrome is suspected to have a genetic component. The objective of this study was to determine if a genetic locus associated with SARDS in Dachshunds could be identified using a genome-wide association study (GWAS). PROCEDURES: Genome-wide association mapping was performed in 15 SARDS-affected and 16 unaffected Dachshunds. Genotyping of three classical DLA class II genes (DLA-DRB1, DLA-DQA1, and DLA-DQB1) was performed in 34 SARDS-affected and 66 unaffected Dachshunds to evaluate for an association in this region. RESULTS: Although no single nucleotide polymorphisms (SNPs) were of genome-wide statistical significance (PBonferroni  < 0.05), 5 of the top 9 SNPs were in the major histocompatibility complex (MHC). Using DLA typing, the allele DLA-DRB1*09401 was identified as a risk factor for the development of SARDS (P = 0.0032, OR = 4.0). The alleles DLA-DQB1*00101 (P = 0.0050, OR = 0.31), DLA-DQA1*00901 (P = 0.0087, OR = 0.33), and a previously identified DLA-DRB1allele described as "DRB1-T" (P = 0.0284, OR = 0.37) were identified as protective factors. CONCLUSIONS: Although far from definitive, association of SARDS with alleles of immunologic importance further supports the hypothesis that autoimmunity may play a role in the pathogenesis of SARDS.

Identification of a Candidate Mutation in the COL1A2 Gene of a Chow Chow With Osteogenesis Imperfecta.

Osteogenesis imperfecta (OI) is a genetic disease that occurs in humans and animals. Individuals with OI exhibit signs of extreme bone fragility and osteopenia with frequent fractures and perinatal lethality in severe cases. In this study, we report the clinical diagnosis of OI in a dog and the use of targeted next-generation sequencing to identify a candidate autosomal dominant mutation in the COL1A2 gene. A 5-month-old male Chow Chow was examined with a fractured left humerus and resolving, bilateral femoral fractures. Radiographs revealed generalized osteopenia and bilateral humeral, radial, and femoral fractures. Targeted next-generation sequencing of genes associated with OI in humans (COL1A1, COL1A2, LEPRE1, SERPINH1, and CRTAP) revealed a G>A heterozygous mutation in the splice donor site of exon 18 of the COL1A2 gene (c.936 + 1G>A). The splice donor mutation was not detected among 91 control dogs representing 21 breeds. A comparative analysis of exon 18 and the exon-intron junction further showed that the mutated splice donor site is conserved among vertebrates. Altogether, these findings reveal a candidate autosomal splice donor site mutation causing OI in an individual Chow Chow.

SERPINB11 frameshift variant associated with novel hoof specific phenotype in Connemara ponies.

Horses belong to the order Perissodactyla and bear the majority of their weight on their third toe; therefore, tremendous force is applied to each hoof. An inherited disease characterized by a phenotype restricted to the dorsal hoof wall was identified in the Connemara pony. Hoof wall separation disease (HWSD) manifests clinically as separation of the dorsal hoof wall along the weight-bearing surface of the hoof during the first year of life. Parents of affected ponies appeared clinically normal, suggesting an autosomal recessive mode of inheritance. A case-control allelic genome wide association analysis was performed (ncases = 15, ncontrols = 24). Population stratification (λ = 1.48) was successfully improved by removing outliers (ncontrols = 7) identified on a multidimensional scaling plot. A genome-wide significant association was detected on chromosome 8 (praw = 1.37x10-10, pgenome = 1.92x10-5). A homozygous region identified in affected ponies spanned from 79,936,024-81,676,900 bp and contained a family of 13 annotated SERPINB genes. Whole genome next-generation sequencing at 6x coverage of two cases and two controls revealed 9,758 SNVs and 1,230 indels within the ~1.7-Mb haplotype, of which 17 and 5, respectively, segregated with the disease and were located within or adjacent to genes. Additional genotyping of these 22 putative functional variants in 369 Connemara ponies (ncases = 23, ncontrols = 346) and 169 horses of other breeds revealed segregation of three putative variants adjacent or within four SERPIN genes. Two of the variants were non-coding and one was an insertion within SERPINB11 that introduced a frameshift resulting in a premature stop codon. Evaluation of mRNA levels at the proximal hoof capsule (ncases = 4, ncontrols = 4) revealed that SERPINB11 expression was significantly reduced in affected ponies (p<0.001). Carrier frequency was estimated at 14.8%. This study describes the first genetic variant associated with a hoof wall specific phenotype and suggests a role of SERPINB11 in maintaining hoof wall structure.

Early Developmental and Evolutionary Origins of Gene Body DNA Methylation Patterns in Mammalian Placentas.

Over the last 20-80 million years the mammalian placenta has taken on a variety of morphologies through both divergent and convergent evolution. Recently we have shown that the human placenta genome has a unique epigenetic pattern of large partially methylated domains (PMDs) and highly methylated domains (HMDs) with gene body DNA methylation positively correlating with level of gene expression. In order to determine the evolutionary conservation of DNA methylation patterns and transcriptional regulatory programs in the placenta, we performed a genome-wide methylome (MethylC-seq) analysis of human, rhesus macaque, squirrel monkey, mouse, dog, horse, and cow placentas as well as opossum extraembryonic membrane. We found that, similar to human placenta, mammalian placentas and opossum extraembryonic membrane have globally lower levels of methylation compared to somatic tissues. Higher relative gene body methylation was the conserved feature across all mammalian placentas, despite differences in PMD/HMDs and absolute methylation levels. Specifically, higher methylation over the bodies of genes involved in mitosis, vesicle-mediated transport, protein phosphorylation, and chromatin modification was observed compared with the rest of the genome. As in human placenta, higher methylation is associated with higher gene expression and is predictive of genic location across species. Analysis of DNA methylation in oocytes and preimplantation embryos shows a conserved pattern of gene body methylation similar to the placenta. Intriguingly, mouse and cow oocytes and mouse early embryos have PMD/HMDs but their placentas do not, suggesting that PMD/HMDs are a feature of early preimplantation methylation patterns that become lost during placental development in some species and following implantation of the embryo.

Genome-wide association studies in dogs and humans identify ADAMTS20 as a risk variant for cleft lip and palate.

Cleft lip with or without cleft palate (CL/P) is the most commonly occurring craniofacial birth defect. We provide insight into the genetic etiology of this birth defect by performing genome-wide association studies in two species: dogs and humans. In the dog, a genome-wide association study of 7 CL/P cases and 112 controls from the Nova Scotia Duck Tolling Retriever (NSDTR) breed identified a significantly associated region on canine chromosome 27 (unadjusted p=1.1 x 10(-13); adjusted p= 2.2 x 10(-3)). Further analysis in NSDTR families and additional full sibling cases identified a 1.44 Mb homozygous haplotype (chromosome 27: 9.29 - 10.73 Mb) segregating with a more complex phenotype of cleft lip, cleft palate, and syndactyly (CLPS) in 13 cases. Whole-genome sequencing of 3 CLPS cases and 4 controls at 15X coverage led to the discovery of a frameshift mutation within ADAMTS20 (c.1360_1361delAA (p.Lys453Ilefs*3)), which segregated concordant with the phenotype. In a parallel study in humans, a family-based association analysis (DFAM) of 125 CL/P cases, 420 unaffected relatives, and 392 controls from a Guatemalan cohort, identified a suggestive association (rs10785430; p =2.67 x 10-6) with the same gene, ADAMTS20. Sequencing of cases from the Guatemalan cohort was unable to identify a causative mutation within the coding region of ADAMTS20, but four coding variants were found in additional cases of CL/P. In summary, this study provides genetic evidence for a role of ADAMTS20 in CL/P development in dogs and as a candidate gene for CL/P development in humans.

Genetic analysis of optic nerve head coloboma in the Nova Scotia Duck Tolling Retriever identifies discordance with the NHEJ1 intronic deletion (collie eye anomaly mutation).

Collie eye anomaly (CEA) encompasses a spectrum of different ophthalmic phenotypes from clinically inconsequential choroidal hypoplasia to blindness from coloboma of the optic nerve head (ONH). A previous study found a 7.8-kb deletion in intron 4 of the NHEJ1 gene to be associated with CEA. A genetic test based on this association is recommended for many breeds, including the Nova Scotia Duck Tolling Retriever (NSDTR). Collection of ONH coloboma-affected NSDTR showed lack of concordance of the NHEJ1 intronic deletion with ONH coloboma. Using genomewide single nucleotide polymorphism (SNP) genotyping in 7 ONH coloboma-affected NSDTR cases and 47 unaffected NSDTR controls with no ophthalmic signs, one SNP, located on chromosome 7, demonstrated genomewide significance. However, high genomic inflation may have confounded the results. Therefore, the genomewide association study was repeated using EMMAX to control for population structure in the cohort of 7 cases and 47 controls. However, no regions of the genome were significantly associated with ONH coloboma. These results failed to document significant association with the CEA locus. Due to the complex genetic etiology of ONH coloboma, the NHEJ1 intronic deletion test results should be carefully considered when making breeding decisions. If the goal is to select for visually competent dogs, our data suggest that eye examinations of puppies would be more effective as a guide in selection of breeding pairs than relying solely on currently available genetic tests.

Developing a 670k genotyping array to tag ~2M SNPs across 24 horse breeds.

BACKGROUND: To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array. RESULTS: Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation. CONCLUSIONS: Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

A LINE-1 insertion in DLX6 is responsible for cleft palate and mandibular abnormalities in a canine model of Pierre Robin sequence.

Cleft palate (CP) is one of the most commonly occurring craniofacial birth defects in humans. In order to study cleft palate in a naturally occurring model system, we utilized the Nova Scotia Duck Tolling Retriever (NSDTR) dog breed. Micro-computed tomography analysis of CP NSDTR craniofacial structures revealed that these dogs exhibit defects similar to those observed in a recognizable subgroup of humans with CP: Pierre Robin Sequence (PRS). We refer to this phenotype in NSDTRs as CP1. Individuals with PRS have a triad of birth defects: shortened mandible, posteriorly placed tongue, and cleft palate. A genome-wide association study in 14 CP NSDTRs and 72 unaffected NSDTRs identified a significantly associated region on canine chromosome 14 (24.2 Mb-29.3 Mb; p(raw )= 4.64 × 10(-15)). Sequencing of two regional candidate homeobox genes in NSDTRs, distal-less homeobox 5 (DLX5) and distal-less homeobox 6 (DLX6), identified a 2.1 kb LINE-1 insertion within DLX6 in CP1 NSDTRs. The LINE-1 insertion is predicted to insert a premature stop codon within the homeodomain of DLX6. This prompted the sequencing of DLX5 and DLX6 in a human cohort with CP, where a missense mutation within the highly conserved DLX5 homeobox of a patient with PRS was identified. This suggests the involvement of DLX5 in the development of PRS. These results demonstrate the power of the canine animal model as a genetically tractable approach to understanding naturally occurring craniofacial birth defects in humans.

Genome-wide SNP and haplotype analyses reveal a rich history underlying dog domestication.

Advances in genome technology have facilitated a new understanding of the historical and genetic processes crucial to rapid phenotypic evolution under domestication. To understand the process of dog diversification better, we conducted an extensive genome-wide survey of more than 48,000 single nucleotide polymorphisms in dogs and their wild progenitor, the grey wolf. Here we show that dog breeds share a higher proportion of multi-locus haplotypes unique to grey wolves from the Middle East, indicating that they are a dominant source of genetic diversity for dogs rather than wolves from east Asia, as suggested by mitochondrial DNA sequence data. Furthermore, we find a surprising correspondence between genetic and phenotypic/functional breed groupings but there are exceptions that suggest phenotypic diversification depended in part on the repeated crossing of individuals with novel phenotypes. Our results show that Middle Eastern wolves were a critical source of genome diversity, although interbreeding with local wolf populations clearly occurred elsewhere in the early history of specific lineages. More recently, the evolution of modern dog breeds seems to have been an iterative process that drew on a limited genetic toolkit to create remarkable phenotypic diversity.

Gonadectomy effects on the risk of immune disorders in the dog: a retrospective study.

BACKGROUND: Gonadectomy is one of the most common procedures performed on dogs in the United States. Neutering has been shown to reduce the risk for some diseases although recent reports suggest increased prevalence for structural disorders and some neoplasias. The relation between neuter status and autoimmune diseases has not been explored. This study evaluated the prevalence and risk of atopic dermatitis (ATOP), autoimmune hemolytic anemia (AIHA), canine myasthenia gravis (CMG), colitis (COL), hypoadrenocorticism (ADD), hypothyroidism (HYPO), immune-mediated polyarthritis (IMPA), immune-mediated thrombocytopenia (ITP), inflammatory bowel disease (IBD), lupus erythematosus (LUP), and pemphigus complex (PEMC), for intact females, intact males, neutered females, and neutered males. Pyometra (PYO) was evaluated as a control condition. RESULTS: Patient records (90,090) from the William R. Pritchard Veterinary Medical Teaching Hospital at the University of California, Davis from 1995 to 2010 were analyzed in order to determine the risk of immune-mediated disease relative to neuter status in dogs. Neutered dogs had a significantly greater risk of ATOP, AIHA, ADD, HYPO, ITP, and IBD than intact dogs with neutered females being at greater risk than neutered males for all but AIHA and ADD. Neutered females, but not males, had a significantly greater risk of LUP than intact females. Pyometra was a greater risk for intact females. CONCLUSIONS: The data underscore the importance of sex steroids on immune function emphasizing a role of these hormones on tissue self-recognition. Neutering is critically important for population control, reduction of reproductive disorders, and offers convenience for owners. Despite these advantages, the analyses of the present study suggest that neutering is associated with increased risk for certain autoimmune disorders and underscore the need for owners to consult with their veterinary practitioner prior to neutering to evaluate possible benefits and risks associated with such a procedure.

Genome-wide association mapping in dogs enables identification of the homeobox gene, NKX2-8, as a genetic component of neural tube defects in humans.

Neural tube defects (NTDs) is a general term for central nervous system malformations secondary to a failure of closure or development of the neural tube. The resulting pathologies may involve the brain, spinal cord and/or vertebral column, in addition to associated structures such as soft tissue or skin. The condition is reported among the more common birth defects in humans, leading to significant infant morbidity and mortality. The etiology remains poorly understood but genetic, nutritional, environmental factors, or a combination of these, are known to play a role in the development of NTDs. The variable conditions associated with NTDs occur naturally in dogs, and have been previously reported in the Weimaraner breed. Taking advantage of the strong linkage-disequilibrium within dog breeds we performed genome-wide association analysis and mapped a genomic region for spinal dysraphism, a presumed NTD, using 4 affected and 96 unaffected Weimaraners. The associated region on canine chromosome 8 (pgenome  =3.0 × 10(-5)), after 100,000 permutations, encodes 18 genes, including NKX2-8, a homeobox gene which is expressed in the developing neural tube. Sequencing NKX2-8 in affected Weimaraners revealed a G to AA frameshift mutation within exon 2 of the gene, resulting in a premature stop codon that is predicted to produce a truncated protein. The exons of NKX2-8 were sequenced in human patients with spina bifida and rare variants (rs61755040 and rs10135525) were found to be significantly over-represented (p=0.036). This is the first documentation of a potential role for NKX2-8 in the etiology of NTDs, made possible by investigating the molecular basis of naturally occurring mutations in dogs.