GWAS using low-pass whole genome sequence reveals a novel locus in canine congenital idiopathic megaesophagus.

Congenital idiopathic megaesophagus (CIM) is a gastrointestinal disorder of dogs wherein the esophagus is dilated and swallowing activity is reduced, causing regurgitation of ingesta. Affected individuals experience weight loss and malnourishment and are at risk for aspiration pneumonia, intussusception, and euthanasia. Great Danes have among the highest incidences of CIM across dog breeds, suggesting a genetic predisposition. We generated low-pass sequencing data for 83 Great Danes and used variant calls to impute missing whole genome single-nucleotide variants (SNVs) for each individual based on haplotypes phased from 624 high-coverage dog genomes, including 21 Great Danes. We validated the utility of our imputed data set for genome-wide association studies (GWASs) by mapping loci known to underlie coat phenotypes with simple and complex inheritance patterns. We conducted a GWAS for CIM with 2,010,300 SNVs, identifying a novel locus on canine chromosome 1 (P-val = 2.76 × 10-10). Associated SNVs are intergenic or intronic and are found in two clusters across a 1.7-Mb region. Inspection of coding regions in high-coverage genomes from affected Great Danes did not reveal candidate causal variants, suggesting that regulatory variants underlie CIM. Further studies are necessary to assess the role of these non-coding variants.

Congenital idiopathic megaesophagus in the German shepherd dog is a sex-differentiated trait and is associated with an intronic variable number tandem repeat in Melanin-Concentrating Hormone Receptor 2.

Congenital idiopathic megaesophagus (CIM) is a gastrointestinal (GI) motility disorder of dogs in which reduced peristaltic activity and dilation of the esophagus prevent the normal transport of food into the stomach. Affected puppies regurgitate meals and water, fail to thrive, and experience complications such as aspiration pneumonia that may necessitate euthanasia. The German shepherd dog (GSD) has the highest disease incidence, indicative of a genetic predisposition. Here, we discover that male GSDs are twice as likely to be affected as females and show that the sex bias is independent of body size. We propose that female endogenous factors (e.g., estrogen) are protective via their role in promoting relaxation of the sphincter between the esophagus and stomach, facilitating food passage. A genome-wide association study for CIM revealed an association on canine chromosome 12 (P-val = 3.12x10-13), with the lead SNPs located upstream or within Melanin-Concentrating Hormone Receptor 2 (MCHR2), a compelling positional candidate gene having a role in appetite, weight, and GI motility. Within the first intron of MCHR2, we identified a 33 bp variable number tandem repeat (VNTR) containing a consensus binding sequence for the T-box family of transcription factors. Across dogs and wolves, the major allele includes two copies of the repeat, whereas the predominant alleles in GSDs have one or three copies. The single-copy allele is strongly associated with CIM (P-val = 1.32x10-17), with homozygosity for this allele posing the most significant risk. Our findings suggest that the number of T-box protein binding motifs may correlate with MCHR2 expression and that an imbalance of melanin-concentrating hormone plays a role in CIM. We describe herein the first genetic factors identified in CIM: sex and a major locus on chromosome 12, which together predict disease state in the GSD with greater than 75% accuracy.

Congenital myasthenic syndrome in Golden Retrievers is associated with a novel COLQ mutation.

BACKGROUND: Congenital myasthenic syndromes (CMSs) are a group of inherited disorders of neuromuscular transmission that may be presynaptic, synaptic, or postsynaptic. Causative mutations have been identified in 4 breeds including the Labrador Retriever, Jack Russell Terrier, Heideterrier, and Danish Pointing Dog. HYPOTHESIS/OBJECTIVE: Clinical and genetic characterization of a neuromuscular disorder in Golden Retriever (GR) puppies. ANIMALS: Four GR puppies from California were evaluated for generalized muscle weakness beginning at weaning. Biological specimens were collected from the affected puppies, and familial information was obtained. Blood or buccal swabs were obtained from 63 unaffected GRs. METHODS: Complete physical, neurological, electrodiagnostic, and histological evaluations and biochemical quantification of muscle acetylcholine receptors were performed. Polymerase chain reaction was used to amplify the 17 exons of COLQ, and sequences were obtained by Sanger sequencing. Variant frequency was assessed in unrelated GRs and a public database. RESULTS: Clinical, neurological, and electrodiagnostic evaluations confirmed a disorder of neuromuscular transmission in a GR family. Sequencing of all exons and splice sites of a primary candidate gene, COLQ, identified a point mutation that predicts an amino acid substitution (G294R). The primary COLQ transcript was absent from affected muscle samples. All affected puppies were homozygous for the mutation, which was not detected outside this GR family or in other breeds. CONCLUSIONS AND CLINICAL IMPORTANCE: We confirmed the diagnosis of a CMS in GR puppies and identified a novel COLQ mutation. The COLQ gene encodes the collagenous tail of acetylcholinesterase, the enzyme responsible for termination of skeletal muscle contraction by clearing acetylcholine at the neuromuscular junction. Clinicians and breeders should be aware of this CMS in GR puppies with an early onset of weakness.

Novel Y Chromosome Retrocopies in Canids Revealed through a Genome-Wide Association Study for Sex.

The lack of an annotated reference sequence for the canine Y chromosome has limited evolutionary studies, as well as our understanding of the role of Y-linked sequences in phenotypes with a sex bias. In genome-wide association studies (GWASs), we observed spurious associations with autosomal SNPs when sex was unbalanced in case-control cohorts and hypothesized that a subset of SNPs mapped to autosomes are in fact sex-linked. Using the Illumina 230K CanineHD array in a GWAS for sex, we identified SNPs that amplify in both sexes but possess significant allele frequency differences between males and females. We found 48 SNPs mapping to 14 regions of eight autosomes and the X chromosome that are Y-linked, appearing heterozygous in males and monomorphic in females. Within these 14 regions are eight genes: three autosomal and five X-linked. We investigated the autosomal genes (MITF, PPP2CB, and WNK1) and determined that the SNPs are diverged nucleotides in retrocopies that have transposed to the Y chromosome. MITFY and WNK1Y are expressed and appeared recently in the Canidae lineage, whereas PPP2CBY represents a much older insertion with no evidence of expression in the dog. This work reveals novel canid Y chromosome sequences and provides evidence for gene transposition to the Y from autosomes and the X.

Length variations within the Merle retrotransposon of canine PMEL: correlating genotype with phenotype.

BACKGROUND: The antisense insertion of a canine short interspersed element (SINEC_Cf) in the pigmentation gene PMEL (or SILV) causes a coat pattern phenotype in dogs termed merle. Merle is a semi-dominant trait characterized by patches of full pigmentation on a diluted background. The oligo(dT) tract of the Merle retrotransposon is long and uninterrupted and is prone to dramatic truncation. Phenotypically wild-type individuals carrying shorter oligo(dT) lengths of the Merle allele have been previously described and termed cryptic merles. Two additional coat patterns, dilute merle (uniform, steely-grey coat) and harlequin merle (white background with black patches), also appear in breeds segregating the Merle allele. RESULTS: Sequencing of all PMEL exons in a dilute and a harlequin merle reveals that variation exists solely within the oligo(dT) tract of the SINEC_Cf insertion. In fragment analyses from 259 dogs heterozygous for Merle, we observed a spectrum of oligo(dT) lengths spanning 25 to 105 base pairs (bp), with ranges that correspond to the four varieties of the merle phenotype: cryptic (25-55 bp), dilute (66-74 bp), standard (78-86 bp), and harlequin (81-105 bp). Somatic contractions of the oligo(dT) were observed in 43% of standard and 51% of harlequin merle dogs. A small proportion (4.6%) of the study cohort inherited de novo contractions or expansions of the Merle allele that resulted in dilute or harlequin coat patterns, respectively. CONCLUSIONS: The phenotypic consequence of the Merle SINE insertion directly depends upon oligo(dT) length. In transcription, we propose that the use of an alternative splice site increases with oligo(dT) length, resulting in insufficient PMEL and a pigment dilution spectrum, from dark grey to complete hypopigmentation. We further propose that during replication, contractions and expansions increase in frequency with oligo(dT) length, causing coat variegation (somatic events in melanocytes) and the spontaneous appearance of varieties of the merle phenotype (germline events).

Beyond the MHC: A canine model of dermatomyositis shows a complex pattern of genetic risk involving novel loci.

Juvenile dermatomyositis (JDM) is a chronic inflammatory myopathy and vasculopathy driven by genetic and environmental influences. Here, we investigated the genetic underpinnings of an analogous, spontaneous disease of dogs also termed dermatomyositis (DMS). As in JDM, we observed a significant association with a haplotype of the major histocompatibility complex (MHC) (DLA-DRB1*002:01/-DQA1*009:01/-DQB1*001:01), particularly in homozygosity (P-val = 0.0001). However, the high incidence of the haplotype among healthy dogs indicated that additional genetic risk factors are likely involved in disease progression. We conducted genome-wide association studies in two modern breeds having common ancestry and detected strong associations with novel loci on canine chromosomes 10 (P-val = 2.3X10-12) and 31 (P-val = 3.95X10-8). Through whole genome resequencing, we identified primary candidate polymorphisms in conserved regions of PAN2 (encoding p.Arg492Cys) and MAP3K7CL (c.383_392ACTCCACAAA>GACT) on chromosomes 10 and 31, respectively. Analyses of these polymorphisms and the MHC haplotypes revealed that nine of 27 genotypic combinations confer high or moderate probability of disease and explain 93% of cases studied. The pattern of disease risk across PAN2 and MAP3K7CL genotypes provided clear evidence for a significant epistatic foundation for this disease, a risk further impacted by MHC haplotypes. We also observed a genotype-phenotype correlation wherein an earlier age of onset is correlated with an increased number of risk alleles at PAN2 and MAP3K7CL. High frequencies of multiple genetic risk factors are unique to affected breeds and likely arose coincident with artificial selection for desirable phenotypes. Described herein is the first three-locus association with a complex canine disease and two novel loci that provide targets for exploration in JDM and related immunological dysfunction.

A CHRNE frameshift mutation causes congenital myasthenic syndrome in young Jack Russell Terriers.

Congenital myasthenic syndromes (CMSs) are a group of rare genetic disorders of the neuromuscular junction resulting in structural or functional causes of fatigable weakness that usually begins early in life. Mutations in pre-synaptic, synaptic and post-synaptic proteins have been demonstrated in human cases, with more than half involving aberrations in nicotinic acetylcholine receptor (AChR) subunits. CMS was first recognized in dogs in 1974 as an autosomal recessive trait in Jack Russell Terriers (JRTs). A deficiency of junctional AChRs was demonstrated. Here we characterize a CMS in 2 contemporary cases of JRT littermates with classic clinical and electromyographic findings, and immunochemical confirmation of an approximately 90% reduction in AChR protein content. Loci encoding the 5 AChR subunits were evaluated using microsatellite markers, and CHRNB1 and CHRNE were identified as candidate genes. Sequences of the splice sites and exons of both genes revealed a single base insertion in exon 7 of CHRNE that predicts a frameshift mutation and a premature stop codon. We further demonstrated this pathogenic mutation in CHRNE in archival tissues from unrelated JRTs studied 34 years ago.

Alleles of the major histocompatibility complex play a role in the pathogenesis of pancreatic acinar atrophy in dogs.

Exocrine pancreatic insufficiency (EPI) is a disease wherein pancreatic acinar cells fail to synthesize and secrete sufficient amounts of digestive enzymes for normal digestion of food. EPI affects many dog breeds, with a dramatically higher prevalence in the German shepherd dog (GSD) population. In this breed and perhaps others, EPI most often results from degeneration of the acinar cells of the pancreas, a hereditary disorder termed pancreatic acinar atrophy (PAA). Evidence of lymphocytic infiltration indicates that PAA is an autoimmune disease, but the genetic etiology remains unclear. Data from global gene expression and single nucleotide polymorphism profiles in the GSD suggest the involvement of the major histocompatibility complex [MHC; dog leukocyte antigen (DLA)]. To determine if alleles of the MHC influence development of EPI, genotyping of polymorphic class I (DLA-88) and II loci (DLA-DRB1, DLA-DQA1, and DLA-DQB1) was carried out for 70 affected and 63 control GSDs, and four-locus haplotypes were determined. One haplotype containing a novel allele of DLA-88 is very highly associated with EPI (OR > 17; P = 0.000125), while two haplotypes were found to confer protection from EPI (P = 0.00087 and 0.0115). Described herein is the genotyping of MHC class I and II loci in a GSD cohort, establishment of four-locus haplotypes, and association of alleles/haplotypes with EPI.

Genome-wide association studies for multiple diseases of the German Shepherd Dog.

The German Shepherd Dog (GSD) is a popular working and companion breed for which over 50 hereditary diseases have been documented. Herein, SNP profiles for 197 GSDs were generated using the Affymetrix v2 canine SNP array for a genome-wide association study to identify loci associated with four diseases: pituitary dwarfism, degenerative myelopathy (DM), congenital megaesophagus (ME), and pancreatic acinar atrophy (PAA). A locus on Chr 9 is strongly associated with pituitary dwarfism and is proximal to a plausible candidate gene, LHX3. Results for DM confirm a major locus encompassing SOD1, in which an associated point mutation was previously identified, but do not suggest modifier loci. Several SNPs on Chr 12 are associated with ME and a 4.7 Mb haplotype block is present in affected dogs. Analysis of additional ME cases for a SNP within the haplotype provides further support for this association. Results for PAA indicate more complex genetic underpinnings. Several regions on multiple chromosomes reach genome-wide significance. However, no major locus is apparent and only two associated haplotype blocks, on Chrs 7 and 12 are observed. These data suggest that PAA may be governed by multiple loci with small effects, or it may be a heterogeneous disorder.

A canine BCAN microdeletion associated with episodic falling syndrome.

Episodic falling syndrome (EFS) is a canine paroxysmal hypertonicity disorder found in Cavalier King Charles spaniels. Episodes are triggered by exercise, stress or excitement and characterized by progressive hypertonicity throughout the thoracic and pelvic limbs, resulting in a characteristic 'deer-stalking' position and/or collapse. We used a genome-wide association strategy to map the EFS locus to a 3.48 Mb critical interval on canine chromosome 7. By prioritizing candidate genes on the basis of biological plausibility, we found that a 15.7 kb deletion in BCAN, encoding the brain-specific extracellular matrix proteoglycan brevican, is associated with EFS. This represents a compelling causal mutation for EFS, since brevican has an essential role in the formation of perineuronal nets governing synapse stability and nerve conduction velocity. Mapping of the deletion breakpoint enabled the development of Multiplex PCR and Multiplex Ligation-dependent Probe Amplification (MLPA) genotyping tests that can accurately distinguish normal, carrier and affected animals. Wider testing of a larger population of CKCS dogs without a history of EFS from the USA revealed that carriers are extremely common (12.9%). The development of molecular genetic tests for the EFS microdeletion will allow the implementation of directed breeding programs aimed at minimizing the number of animals with EFS and enable confirmatory diagnosis and pharmacotherapy of affected dogs.

A missense mutation in the 20S proteasome ß2 subunit of Great Danes having harlequin coat patterning.

Harlequin is a pigmentary trait of the domestic dog that is controlled by two autosomal loci: the melanosomal gene, SILV, and a modifier gene, harlequin (H), previously localized to chromosome 9. Heterozygosity for a retrotransposon insertion in SILV and a mutation in H causes a pattern of black patches on a white background. Homozygosity for H is embryonic lethal. Fine mapping of the harlequin locus revealed a 25 kb interval wherein all harlequin Great Danes are heterozygous for a common haplotype. This region contains one gene, PSMB7, which encodes the ß2 catalytic subunit of the proteasome. Sequence analysis identified a coding variant in exon 2 that segregates with harlequin patterning. The substitution predicts the replacement of a highly conserved valine with a glycine. Described herein is the identification of a naturally-occurring mutation of the ubiquitin proteasome system that is associated with a discernable phenotype of dogs.

Evaluation of quantitative trait loci for hip dysplasia in Labrador Retrievers.

OBJECTIVE: To identify the quantitative trait loci (QTL) that contribute to hip dysplasia in dogs. ANIMALS: 192 Labrador Retrievers. PROCEDURES: Hip dysplasia was measured by use of the Norberg angle (NA), dorsolateral subluxation (DLS) score, and distraction index (DI). Genome-wide screening was conducted by use of 276 unique microsatellites. Linkage analysis was performed with a variance-based linear model. Logarithm of the odds (LOD) scores were reported when values were > 2.0. RESULTS: Canis familiaris autosomes (CFAs) 01, 02, 10, 20, 22, and 32 harbored significant QTL at LOD scores > 2.0. Among the 6 QTL, the QTL on CFA02 had not been reported to harbor QTL for hip dysplasia. The highest LOD score of 3.32 on CFA20 contributed to the second principal component of the DLS score and NA of the right hip joint. The QTL that was mapped on CFA01 (LOD score of 3.13 at 55 centimorgans) was located on the same chromosome reported to harbor a QTL for hip dysplasia in Portuguese Water Dogs and German Shepherd Dogs. In this study, CFAs 10, 20, 22, and 32 harbored QTL for hip dysplasia that have been identified in a Labrador Retriever-Greyhound pedigree and in German Shepherd Dogs. CONCLUSIONS AND CLINICAL RELEVANCE: Multiple QTL were clearly involved with hip dysplasia. Identification of these QTL will enable fine-resolution mapping and subsequent assessment of candidate genes within the refined intervals to enable researchers to develop genetic screening tests and preventative and novel therapeutic regimens.

Identification of quantitative trait loci for osteoarthritis of hip joints in dogs.

OBJECTIVE: To identify quantitative trait loci (QTL) associated with osteoarthritis (OA) of hip joints of dogs by use of a whole-genome microsatellite scan. ANIMALS: 116 founder, backcross, F1, and F2 dogs from a crossbred pedigree. PROCEDURES: Necropsy scores and an optimized set of 342 microsatellite markers were used for interval mapping by means of a combined backcross and F2 design module from an online statistical program. Breed and sex were included in the model as fixed effects. Age of dog at necropsy and body weight at 8 months of age were also included in the model as covariates. The chromosomal location at which the highest F score was obtained was considered the best estimate of a QTL position. Chromosome-wide significance thresholds were determined empirically from 10,000 permutations of marker genotypes. RESULTS: 4 chromosomes contained putative QTL for OA of hip joints in dogs at the 5% chromosome-wide significance threshold: chromosomes 5, 18, 23, and 31. CONCLUSIONS AND CLINICAL RELEVANCE: Osteoarthritis of canine hip joints is a complex disease to which many genes and environmental factors contribute. Identification of contributing QTL is a strategy to elucidate the genetic mechanisms that underlie this disease. Refinement of the putative QTL and subsequent candidate gene studies are needed to identify the genes involved in the disease process.

Genome-wide linkage scan localizes the harlequin locus in the Great Dane to chromosome 9.

Harlequin is a coat pattern of the Great Dane characterized by ragged patches of full color on a white background. Harlequin patterning is a bigenic trait, resulting from the interaction of the merle allele of SILV, and a dominant modifier locus, H. Breeding data suggest that H is embryonic recessive lethal and that all harlequins are Hh. To identify linkage with the harlequin phenotype, 46 Great Danes from 5 pedigrees were genotyped for 280 microsatellite markers in a whole genome screen. One marker on the telomeric end of chromosome 9 was suggestive of linkage. Fine mapping of this region using additional microsatellite markers and 10 Great Danes from a sixth pedigree resulted in significant LOD scores for 2 markers. Reported herein is linkage mapping of the H locus to a 3.27 Mb region of chromosome 9 containing approximately 20 genes.

Identification, amplification and characterization of miR-17-92 from canine tissue.

Recently, a novel group of genes encoding small RNA molecules, termed microRNAs (miRNAs), has been discovered to play a vital role in eukaryotic gene expression. Known to act in a post-transcriptional fashion, miRNAs can inhibit translation by binding to messenger RNA (mRNA) or by targeting mRNA for degradation. A search of genetic databases revealed significant conservation of miRNA genes between the domestic dog and the human. This finding suggests that expression patterns may also be conserved. Proof of principle experiments, including serial dilutions and sequencing, were performed to verify that primers made to amplify human mature miRNAs can be used to amplify canine miRNAs, providing that the mature sequences are conserved. TaqMan Real-time PCR techniques were used to isolate the first miRNA mature products from canine tissues. The expression levels of miR-17-3p, miR-17-5p, miR-18, miR-19a, miR-19b, miR-20, and miR-92 were evaluated in five canine tissues (heart, lung, brain, kidney, and liver) using the delta-delta Ct (critical threshold) method.

Understanding hereditary diseases using the dog and human as companion model systems.

Animal models are requisite for genetic dissection of, and improved treatment regimens for, human hereditary diseases. While several animals have been used in academic and industrial research, the primary model for dissection of hereditary diseases has been the many strains of the laboratory mouse. However, given its greater (than the mouse) genetic similarity to the human, high number of naturally occurring hereditary diseases, unique population structure, and the availability of the complete genome sequence, the purebred dog has emerged as a powerful model for study of diseases. The major advantage the dog provides is that it is afflicted with approximately 450 hereditary diseases, about half of which have remarkable clinical similarities to corresponding diseases of the human. In addition, humankind has a strong desire to cure diseases of the dog so these two facts make the dog an ideal clinical and genetic model. This review highlights several of these shared hereditary diseases. Specifically, the canine models discussed herein have played important roles in identification of causative genes and/or have been utilized in novel therapeutic approaches of interest to the dog and human.

Expression of multiple human endogenous retrovirus surface envelope proteins in ovarian cancer.

Individual classes of human endogenous retrovirus (HERV) genes and proteins are expressed in cancer, but expression of more than one type of HERV is rare. We report here the expression of multiple HERV genes and proteins in ovarian cell lines and tissues. Expression of HERV-K env mRNA was greater in ovarian epithelial tumors than in normal ovarian tissues (N = 254). The expression of this protein on the surface and in the cytoplasm of ovarian cancer cells was confirmed using anti-HERV-K specific antibody by flow cytometric analysis. The frequency of expression of HERV-K env protein in multitissue microarrays (N = 641) was determined by immunohistochemistry and a significant correlation with tumor histotype was found. A significantly increased expression of HERV-K was observed in tumors with low malignant potential and low grade, relative to expression in normal ovarian tissues. The increase in expression of HERV-K env protein took place in a stepwise fashion in serous papillary adenocarcinoma. Interestingly, we found that other classes of HERV env mRNAs, including ERV3 and HERV-E, are expressed in the same ovarian cancer tissues that expressed HERV-K. Furthermore, anti-HERV antibodies including anti-ERV3 (30%), anti-HERV-E (40%) and anti-HERV-K (55%) were detected in patients with ovarian cancer, but not in normal female controls. HERV env proteins are frequently transcribed and translated in ovarian epithelial tumors, and multiple HERV families are detectable in ovarian cancer. HERV env proteins, and especially those expressed on the cell surface, may serve as novel tumor targets for detection, diagnosis and immunotherapy of ovarian cancer.

Clinical and genetic assessments of hip joint laxity in the Boykin spaniel.

Canine hip dysplasia (CHD) is characterized by a malformation of the hip joint that leads to joint laxity and consequential degenerative joint disease. The most widely used method for diagnosis of CHD is the ventrodorsal hip-extended radiologic view, commonly referred to as the Orthopedic Foundation for Animals (OFA) method. The method of the University of Pennsylvania Hip Improvement Program (PennHIP), an alternative technique that is based on hip joint laxity, provides a quantitative assessment, the distraction index (DI), of the likelihood of the development of CHD because of increased laxity in the hip joint. Linear regression analysis showed that, across many breeds of dog, the incidence of CHD, as defined by the OFA, is positively correlated with the mean DI, the determination coefficient (r2) being 26%. We used families of Boykin spaniels (BSs) to determine the level of joint laxity in the breed and to conduct an initial whole-genome screening to identify markers that co-segregate with increased joint laxity. Although there was a positive correlation between the incidence of hip dysplasia and increased joint laxity, we did not find significant linkage in the 28 BSs that underwent genotyping, likely owing to the small size of the pedigree.

Chromosome-specific microsatellite multiplex sets for linkage studies in the domestic dog.

To expedite linkage studies and positional cloning efforts in the dog, Minimal Screening Set 2 (MSS-2) of 327 canine microsatellite markers has been multiplexed into chromosome-specific panels. MSS-2 provides 9 Mb coverage of the canine genome with no gaps larger than 17.1 Mb and is the most recent and comprehensive set of microsatellites available for whole-genome scans. Markers were labeled with fluorescent dyes based on locations and expected product sizes to facilitate the multiplexing of a maximum number of markers for each chromosome. All markers are amplified using a single thermal cycling program and PCR mix and are optimized for resolution on an ABI 3100 genetic analyzer. Sixty-nine chromosome-specific panels were created by coamplification of a maximum number of markers and subsequent coloading of the remaining markers.