A novel unstable duplication upstream of HAS2 predisposes to a breed-defining skin phenotype and a periodic fever syndrome in Chinese Shar-Pei dogs

Chromosome-scale assembly with improved annotation provides insights into breed-wide genomic structure and diversity in domestic cats

INTRODUCTION

Comprehensive genomic resources offer insights into biological features, including traits/disease-related genetic loci. The current reference genome assembly for the domestic cat (Felis catus), Felis_Catus_9.0 (felCat9), derived from sequences of the Abyssinian cat, may inadequately represent the general cat population, limiting the extent of deducible genetic variations.

OBJECTIVES

The goal was to develop Anicom American Shorthair 1.0 (AnAms1.0), a reference-grade chromosome-scale cat genome assembly.

METHODS

In contrast to prior assemblies relying on Abyssinian cat sequences, AnAms1.0 was constructed from the sequences of more popular American Shorthair breed, which is related to more breeds than the Abyssinian cat. By combining advanced genomics technologies, including PacBio long-read sequencing and Hi-C- and optical mapping data-based sequence scaffolding, we compared AnAms1.0 to existing Felidae genome assemblies (20 scaffolds, scaffolds N50 > 150 Mbp). Homology-based and ab initio gene annotation through Iso-Seq and RNA-Seq was used to identify new coding genes and splice variants.

RESULTS

AnAms1.0 demonstrated superior contiguity and accuracy than existing Felidae genome assemblies. Using AnAms1.0, we identified over 1.5 thousand structural variants and 29 million repetitions compared to felCat9. Additionally, we identified > 1,600 novel protein-coding genes. Notably, olfactory receptor structural variants and cardiomyopathy-related variants were identified.

CONCLUSION

AnAms1.0 facilitates the discovery of novel genes related to normal and disease phenotypes in domestic cats. The analyzed data are publicly accessible on Cats-I (https://cat.annotation.jp/), which we established as a platform for accumulating and sharing genomic resources to discover novel genetic traits and advance veterinary medicine.

Shar-Pei autoinflammatory disease with skin necrosis due to neutrophilic vasculitis: a case report

A 4-month-old female Shar-pei dog was admitted with apathy, anorexia, and vomiting. The patient had an appropriate vaccination protocol, with the last vaccine administered 2.5 weeks prior to the onset of clinical signs. Physical examination revealed tachycardia, fever and swelling of the tibiotarsal joints. Several diagnostic tests including complete blood cell count, biochemistry profile, urinalysis, thoracic radiographs, hind limbs radiographs, abdominal ultrasound, and infectious diseases tests, were conducted to determine the underlying cause. Shar-Pei Auto-inflammatory Disease (SPAID) was diagnosed. Additionally, the patient developed skin necrosis in the inner aspect of the tibiotarsal joints as a complication. A skin biopsy revealed cutaneous vasculopathy causing degeneration, abrupt ulceration, and ischemic necrosis with intense neutrophilic inflammation of the dermis and subcutis. Moreover, a hospital-acquired infection was identified by skin culture. Debridement of the necrotic skin was performed, and due to its' severity and extent, the wound was closed by secondary intention. A diagnostic protocol and the therapeutic dose of acetylsalicylic acid, which led to clinical improvement, are recommended here. The patient has continued to present episodic manifestations of SPAID mainly fever and swelling of the tibiotarsal joints, but there has been no recurrence of necrosis or other cutaneous lesion in the last two years.

Genome-wide identification of copy number variations in wrinkled skin cases of Xiang pigs

Copy number variation (CNV) tends to occur in genetically enriched regions and is likely associated with a number of complex diseases such as skin aging. In this study, we investigated the genome-wide CNVs in 20 wrinkled skin cases (WSC) of Xiang pigs and 63 controls, and identified 7893 copy number variable regions (CNVRs). We estimated the F-statistic (Fst) at each locus and identified that 93 case-controls stratified CNVRs (Fst > = 0.15) overlapped with 87 known genes. Functional enrichment analysis showed that most of these genes were predominantly enriched in pathways and terms related to the extracellular matrix. Finally, we found that some CNVs were predicted to have high effects on genes such as VCAN, TIMP1 and FOXO1 through transcriptional amplification, transcript ablation and so on. Most of the genes overlapped with those CNVRs have been reported to be related to aging in human or animals. The copy numbers presented the positive correlations with the transcript level of the genes in skins between the cases and controls. Our results suggested that those 22 CNVRs, including 19 CNV losses and 3 CNV gains, were putatively associated with the skin wrinkle of Xiang pigs.

Copy number variant scan in more than four thousand Holstein cows bred in Lombardy, Italy

Copy Number Variants (CNV) are modifications affecting the genome sequence of DNA, for instance, they can be duplications or deletions of a considerable number of base pairs (i.e., greater than 1000 bp and up to millions of bp). Their impact on the variation of the phenotypic traits has been widely demonstrated. In addition, CNVs are a class of markers useful to identify the genetic biodiversity among populations related to adaptation to the environment. The aim of this study was to detect CNVs in more than four thousand Holstein cows, using information derived by a genotyping done with the GGP (GeneSeek Genomic Profiler) bovine 100K SNP chip. To detect CNV the SVS 8.9 software was used, then CNV regions (CNVRs) were detected. A total of 123,814 CNVs (4,150 non redundant) were called and aggregated into 1,397 CNVRs. The PCA results obtained using the CNVs information, showed that there is some variability among animals. For many genes annotated within the CNVRs, the role in immune response is well known, as well as their association with important and economic traits object of selection in Holstein, such as milk production and quality, udder conformation and body morphology. Comparison with reference revealed unique CNVRs of the Holstein breed, and others in common with Jersey and Brown. The information regarding CNVs represents a valuable resource to understand how this class of markers may improve the accuracy in prediction of genomic value, nowadays solely based on SNPs markers.

Purebred dogs show higher levels of genomic damage compared to mixed breed dogs

Inbreeding is a common phenomenon in small, fragmented or isolated populations, typical conditions of many threatened species. In the present paper, we used a new non-invasive approach based on the buccal micronucleus assay to evaluate the possible relationships between inbreeding and genomic damage using the dog as model species. In particular, we assessed the frequencies of micronuclei and other nuclear aberrations in a group of purebred dogs (n = 77), comparing the obtained data with those from a control group represented by mixed breed dogs (n = 75). We found a significant increase of micronuclei, nuclear buds and total nuclear aberrations frequencies in purebred dogs compared to mixed-bred dogs. The absence of significant differences in the frequency of micronuclei and other nuclear aberrations amongst different breeds reinforces the hypothesis that the observed increased genomic damage amongst purebred dogs may not be due to a different genomic instability typical of a particular breed, but to inbreeding itself. This hypothesis is further confirmed by the fact that other endogen confounding factors, such as sex, age and weight, do not contribute significantly to the increase of genomic damage observed amongst purebred dogs. In conclusion, results presented in this study showed that, in purebred dogs, inbreeding may increase the levels of genomic damage. Considering that genomic damage is associated with increased physiological problems affecting animal health, the results we obtained may represent a stimulus to discourage the use of intensive inbreeding practices in captive populations and to reduce the fragmentation of wild populations.

Ethical Concerns about Fashionable Dog Breeding

The historical relationship between humans and dogs has involved selective breeding for various purposes, such as hunting, guarding, and service roles. However, over time, there has been a shift in preferences from functionality to aesthetics, which has influenced the diverse sizes, shapes, and coats of dog breeds. This review looks at fashionable dog breeding and questions the ethics of prioritising looks over health and behaviour. It aims to alert potential owners, breeders, and regulators to the importance of considering a dog's overall well-being, not just its appearance, which has resulted in fad breeding, leading to genetic disorders, health issues, and a loss of biodiversity. Ethical concerns arise from breeding brachycephalic breeds with respiratory conditions, inbreeding causing inherited disorders, and overbreeding popular breeds while shelter dogs remain unadopted. Additionally, the impact of cosmetic surgeries on popular dog breeds, as well as the neglect of behavioural traits in favour of physical characteristics and strict breeding practices are also considered. The current breeding model can have a negative impact on the emotional and cognitive well-being of dogs, resulting in issues such as aggression, anxiety, and other behavioural problems that can significantly reduce their overall quality of life. Unregulated breeding practices and the demand for rare breeds can lead to illegal breeding, compromising animal welfare. Prospective owners, veterinarians, kennel clubs, and legislators all need to play a responsible role in protecting animals.

Genetic Changes in Mastocytes and Their Significance in Mast Cell Tumor Prognosis and Treatment

Mast cell tumors are a large group of diseases occurring in dogs, cats, mice, as well as in humans. Systemic mastocytosis (SM) is a disease involving the accumulation of mast cells in organs. KIT gene mutations are very often seen in abnormal mast cells. In SM, high KIT/CD117 expression is observed; however, there are usually no KIT gene mutations present. Mastocytoma (MCT)-a form of cutaneous neoplasm-is common in animals but quite rare in humans. KIT/CD117 receptor mutations were studied as the typical changes for human mastocytosis. In 80% of human cases, the KIT gene substitution p.D816H was present. In about 25% of MCTs, metastasis was observed. Changes in the gene expression of certain genes, such as overexpression of the DNAJ3A3 gene, promote metastasis. In contrast, the SNORD93 gene blocks the expression of metastasis genes. The panel of miR-21-5p, miR-379, and miR-885 has a good efficiency in discriminating healthy and MCT-affected dogs, as well as MCT-affected dogs with and without nodal metastasis. Further studies on the pathobiology of mast cells can lead to clinical improvements, such as better MCT diagnosis and treatment. Our paper reviews studies on the topic of mast cells, which have been carried out over the past few years.

Is it now time to iron out the wrinkles? Health of Shar Pei dogs under primary veterinary care in the UK

BACKGROUND

The Shar Pei is a common dog breed with a distinctive appearance caused by hyaluronosis that has been linked with several health conditions. Anonymised primary-care veterinary clinical records were explored to extract data on the demography, common disorders and mortality of Shar Pei in the UK in 2013.

RESULTS

The study population of 455,557 dogs included 1913 (0.42%) Shar Pei. The mean adult bodyweight was 22.26 kg. The most prevalent fine-level precision disorders were entropion (prevalence 17.88%, 95% CI: 16.16-19.59), otitis externa (16.36%, 95% CI: 14.70-18.02), ear disorders (6.69%, 95% CI: 5.57-7.81), aggression (5.23%, 95% CI: 4.23-6.22), and pyoderma (4.29%, 95% CI: 3.38-5.19). The most prevalent disorder groups were ophthalmologic (prevalence = 22.27%, 95% CI: 20.40-24.13), dermatological (21.01%, 95% CI: 19.19-22.84), aural (18.66%, 95% CI: 16.92-20.41), traumatic injury (7.53%, 95% CI: 6.35-8.71) and behavioural (7.21%, 95% CI: 6.05-8.37). The median longevity of 190 Shar Pei that died during the study period was 7.28 years (IQR 5.04-10.05, range 0.04-15.04). Of 184 deaths with a recorded method of death, 157 (85.33%) deaths involved euthanasia and 27 (14.67%) deaths were unassisted. Among 136 (71.58%) deaths with a recorded biomedical cause of death, the most common causes of death at group level precision were neoplasia (15.44%, 95% CI: 9.37-21.51), renal disorders (13.24%, 95% CI: 7.54-18.93), and behavioural disorders (11.03%, 95% CI: 5.76-16.29).

CONCLUSIONS

Almost one fifth of Shar Pei receive veterinary care each year for entropion, a condition linked strongly with the extreme conformation of thickened and folded skin and bristly hair that characterises the Shar Pei breed. Several other common disorders are also linked to hyaluronosis. Current UK legislation can help support efforts to avoid breeding or acquiring animals with extreme conformations and to promote adequate veterinary care for already-owned animals with extreme conformations.

Hemangiosarcoma in dogs as a potential non-rodent animal model for drug discovery research of angiosarcoma in humans

Since the domestication of dogs 10,000 years ago, they have shared their living environment with humans and have co-evolved. The breeding process that dogs have undergone in only a few centuries has led to a significant accumulation of specific genetic alterations that could induce particular diseases in certain breeds. These canine diseases are similar to what is found in humans with several differences; therefore, comparing such diseases occurring in humans and dogs can help discover novel disease mechanisms, pathways, and causal genetic factors. Human angiosarcoma (AS) and canine hemangiosarcoma (HSA), which are sarcomas originating from endothelium, are examples of diseases shared between humans and dogs. They exhibit similar characteristics and clinical behaviors, although with some critical differences resulting from evolution. In this review, we will describe the similarities and differences in terms of clinical and molecular characteristics between human AS and canine HSA, and discuss how these similarities and differences can be applied to advance the treatment of these diseases.

A survey on Shar Pei autoinflammatory disease in the United Kingdom

OBJECTIVES

To characterise the fever episodes attributed to Shar Pei autoinflammatory disease and to identify common diagnostic and management strategies in the United Kingdom. A secondary objective was to determine risk factors associated with Shar Pei autoinflammatory disease fever episodes.

METHODS

A retrospective survey was performed to characterise episodes of Shar Pei autoinflammatory disease fever and to identify commonly used treatments in affected dogs. Clinical data were collected from owners and veterinarians. Frequencies of previously proposed risk factors (skin thickness and folding, muzzle conformation) and comorbid conditions were compared between dogs that had exhibited fever episodes consistent with Shar Pei autoinflammatory disease and those who had not.

RESULTS

At least one episode of fever attributed to Shar Pei autoinflammatory disease was reported in 52 of 106 (49%) Shar Pei. Nine other dogs had fever episodes consistent with Shar Pei autoinflammatory disease reported by their owners but not by veterinarians. Median rectal temperature at presentation for Shar Pei autoinflammatory disease fever was 40.1°C [104.2°F] (39.9 to 41.3°C [103.8 to 106.3°F]) and owners reported associated hyporexia (n=33, 63%) and vomiting (n=8, 15%) more frequently than veterinary records (n=22, 42% and n=0, 0%, respectively). The median number of veterinary appointments for Shar Pei autoinflammatory disease was two per dog (1 to 15) while owners reported a median of four episodes per dog per year. None of the assessed phenotypic variants or comorbidities were significantly associated with exhibiting Shar Pei autoinflammatory disease fever episodes.

CLINICAL SIGNIFICANCE

Episodes of Shar Pei autoinflammatory disease fever were reported approximately twice as frequently by owners compared to veterinary records, suggesting the burden of this condition may be underestimated by veterinarians. Specific risk factors for Shar Pei autoinflammatory disease fever were not identified.

Identification of Genetic Risk Factors for Monogenic and Complex Canine Diseases

Advances in DNA sequencing and other technologies have greatly facilitated the identification of genetic risk factors for inherited diseases in dogs. We review recent technological developments based on selected examples from canine disease genetics. The identification of disease-causing variants in dogs with monogenic diseases may become a widely employed diagnostic approach in clinical veterinary medicine in the not-too-distant future. Diseases with complex modes of inheritance continue to pose challenges to researchers but have also become much more tangible than in the past. In addition to strategies for identifying genetic risk factors, we provide some thoughts on the interpretation of sequence variants that are largely inspired by developments in human clinical genetics.

Evaluation of Genetic Diversity in Dog Breeds Using Pedigree and Molecular Analysis: A Review

Domestic dogs are important for many economic and social reasons, and they have become a well-known model species for human disease. According to research, dog breeds exhibit significant levels of inbreeding and genetic diversity loss, decreasing the population’s ability to adapt in certain conditions, and indicating the need of conservation strategies. Before the development of molecular markers, pedigree information was used for genetic diversity management. In recent years, genomic tools are frequently applied for accurate estimation of genetic diversity and improved genetic conservation due to incomplete pedigrees and pedigree errors. The most frequently used molecular markers include PCR-based microsatellite markers (STRs) and DNA sequencing-based single-nucleotide polymorphism markers (SNP). The aim of this review was to highlight genetic diversity studies on dog breeds conducted using pedigree and molecular markers, as well as the importance of genetic diversity conservation in increasing the adaptability and survival of dog breed populations.

Detection of genomic structure variants associated with wrinkled skin in Xiang pig by next generation sequencing

Wrinkling is prominent manifestation of aging skin. A mutant phenotype characterized by systemic wrinkles and thickened skin was discovered in Xiang pig populations with incidence about 1-3%. The feature in histological structure was epidermal hyperplasia and thickening, collagen fibers disorder. To uncover genetic mechanisms for the mutant phenotype of Xiang pigs with systemic wrinkle (WXP), a genome-wide of structural variations (SVs) in WXP was described by next generation resequencing, taking Xiang pigs (XP) and European pigs (EUP) as compares. Total of 32,308 SVs were detected from three pig groups and 965 SVs were identified specifically from WXP, involving 481 protein-coding genes. These genes were mainly enriched in nuclear structure, ECM components and immunomodulatory pathways. According to gene function and enrichment analysis, we found that 65 candidate SVs in 59 protein genes were probably related with the systemic wrinkle of WXP. Of these, several genes are reported to be associate with aging, such as EIF4G2, NOLC1, XYLT1, FUT8, MDM2 and so on. The insertion/deletion and duplication variations of SVs in these genes resulted in the loss of stop-codon or frameshift mutation, and aberrant alternative splicing of transcripts. These genes are involved in cell lamin filament, intermediate filament cytoskeleton, supramolecular complex, cell differentiation and regulation of macromolecule metabolic process etc. Our study suggested that the loss of function or aberrant expression of these genes lead to structural disorder of nuclear and the extracellular matrix (ECM) in skin cells, which probably was the genetic mechanisms for the mutant phenotype of systemic skin wrinkle of Xiang pig.

Genetics of inherited skin disorders in dogs

Canine genodermatoses represent a broad spectrum of diseases with diverse phenotypes. Modern genetic technology including whole genome sequencing has expedited the identification of novel genes and greatly simplified the establishment of genetic diagnoses in such heterogeneous disorders. The precise genetic diagnosis of a heritable skin disorder is essential for the appropriate counselling of owners regarding the course of the disease, prognosis and implications for breeding. Understanding the underlying pathophysiology is a prerequisite to developing specific, targeted or individualized therapeutic approaches. This review aims to create a comprehensive overview of canine genodermatoses and their respective genetic aetiology known to date. Raising awareness of genodermatoses in dogs is important and this review may help clinicians to apply modern genetics in daily clinical practice, so that a precise diagnoses can be established in suspected genodermatoses.

Canine coat pigmentation genetics: a review

Our understanding of canine coat colour genetics and the associated health implications is developing rapidly. To date, there are 15 genes with known roles in canine coat colour phenotypes. Many coat phenotypes result from complex and/or epistatic genetic interactions among variants within and between loci, some of which remain unidentified. Some genes involved in canine pigmentation have been linked to aural, visual and neurological impairments. Consequently, coat pigmentation in the domestic dog retains considerable ethical and economic interest. In this paper we discuss coat colour phenotypes in the domestic dog, the genes and variants responsible for these phenotypes and any proven coat colour-associated health effects.

A conditional multi-trait sequence GWAS discovers pleiotropic candidate genes and variants for sheep wool, skin wrinkle and breech cover traits

Background

Imputation to whole-genome sequence is now possible in large sheep populations. It is therefore of interest to use this data in genome-wide association studies (GWAS) to investigate putative causal variants and genes that underpin economically important traits. Merino wool is globally sought after for luxury fabrics, but some key wool quality attributes are unfavourably correlated with the characteristic skin wrinkle of Merinos. In turn, skin wrinkle is strongly linked to susceptibility to “fly strike” (Cutaneous myiasis), which is a major welfare issue. Here, we use whole-genome sequence data in a multi-trait GWAS to identify pleiotropic putative causal variants and genes associated with changes in key wool traits and skin wrinkle.

Results

A stepwise conditional multi-trait GWAS (CM-GWAS) identified putative causal variants and related genes from 178 independent quantitative trait loci (QTL) of 16 wool and skin wrinkle traits, measured on up to 7218 Merino sheep with 31 million imputed whole-genome sequence (WGS) genotypes. Novel candidate gene findings included the MAT1A gene that encodes an enzyme involved in the sulphur metabolism pathway critical to production of wool proteins, and the ESRP1 gene. We also discovered a significant wrinkle variant upstream of the HAS2 gene, which in dogs is associated with the exaggerated skin folds in the Shar-Pei breed.

Conclusions

The wool and skin wrinkle traits studied here appear to be highly polygenic with many putative candidate variants showing considerable pleiotropy. Our CM-GWAS identified many highly plausible candidate genes for wool traits as well as breech wrinkle and breech area wool cover.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00651-0.

Analysis of copy number variation in dogs implicates genomic structural variation in the development of anterior cruciate ligament rupture

Anterior cruciate ligament (ACL) rupture is an important condition of the human knee. Second ruptures are common and societal costs are substantial. Canine cranial cruciate ligament (CCL) rupture closely models the human disease. CCL rupture is common in the Labrador Retriever (5.79% prevalence), ~100-fold more prevalent than in humans. Labrador Retriever CCL rupture is a polygenic complex disease, based on genome-wide association study (GWAS) of single nucleotide polymorphism (SNP) markers. Dissection of genetic variation in complex traits can be enhanced by studying structural variation, including copy number variants (CNVs). Dogs are an ideal model for CNV research because of reduced genetic variability within breeds and extensive phenotypic diversity across breeds. We studied the genetic etiology of CCL rupture by association analysis of CNV regions (CNVRs) using 110 case and 164 control Labrador Retrievers. CNVs were called from SNPs using three different programs (PennCNV, CNVPartition, and QuantiSNP). After quality control, CNV calls were combined to create CNVRs using ParseCNV and an association analysis was performed. We found no strong effect CNVRs but found 46 small effect (max(T) permutation P<0.05) CCL rupture associated CNVRs in 22 autosomes; 25 were deletions and 21 were duplications. Of the 46 CCL rupture associated CNVRs, we identified 39 unique regions. Thirty four were identified by a single calling algorithm, 3 were identified by two calling algorithms, and 2 were identified by all three algorithms. For 42 of the associated CNVRs, frequency in the population was <10% while 4 occurred at a frequency in the population ranging from 10–25%. Average CNVR length was 198,872bp and CNVRs covered 0.11 to 0.15% of the genome. All CNVRs were associated with case status. CNVRs did not overlap previous canine CCL rupture risk loci identified by GWAS. Associated CNVRs contained 152 annotated genes; 12 CNVRs did not have genes mapped to CanFam3.1. Using pathway analysis, a cluster of 19 homeobox domain transcript regulator genes was associated with CCL rupture (P = 6.6E-13). This gene cluster influences cranial-caudal body pattern formation during embryonic limb development. Clustered genes were found in 3 CNVRs on chromosome 14 (HoxA), 28 (NKX6-2), and 36 (HoxD). When analysis was limited to deletion CNVRs, the association was strengthened (P = 8.7E-16). This study suggests a component of the polygenic risk of CCL rupture in Labrador Retrievers is associated with small effect CNVs and may include aspects of stifle morphology regulated by homeobox domain transcript regulator genes.

Atherosclerosis-associated endothelial cell apoptosis by miRNA let7-b-mediated downregulation of HAS-2

MicroRNAs (miRNAs) play essential roles in the regulation and pathophysiology of various types of human diseases including atherosclerosis. Increasing numbers of miRNAs have been identified to be important regulators in the progression of atherosclerosis by regulating gene expression. However, functional miRNAs and the underlying mechanisms involved in atherosclerosis need fully elucidation. In the present study, the function of miRNA let-7b was investigated in human aortic endothelial cells (HAECs). The results showed that downregulation of let-7b in the high-fat diet mice and HAECs was inversely correlated with the expression level of HAS-2. upregulation of let-7b significantly reduced apoptosis of HAECs. The results also revealed that HAS-2 was a target gene of let-7b and HAS-2 reduction reversed the antiapoptotic effect of let-7b through regulation of the P13K/Akt pathway. These results together suggest the potential of regulating the let-7b expression and endothelial apoptosis against development and progression of atherosclerosis.

Multiple FGF4 Retrocopies Recently Derived within Canids

Two transcribed retrocopies of the fibroblast growth factor 4 (FGF4) gene have previously been described in the domestic dog. An FGF4 retrocopy on chr18 is associated with disproportionate dwarfism, while an FGF4 retrocopy on chr12 is associated with both disproportionate dwarfism and intervertebral disc disease (IVDD). In this study, whole-genome sequencing data were queried to identify other FGF4 retrocopies that could be contributing to phenotypic diversity in canids. Additionally, dogs with surgically confirmed IVDD were assayed for novel FGF4 retrocopies. Five additional and distinct FGF4 retrocopies were identified in canids including a copy unique to red wolves (Canis rufus). The FGF4 retrocopies identified in domestic dogs were identical to domestic dog FGF4 haplotypes, which are distinct from modern wolf FGF4 haplotypes, indicating that these retrotransposition events likely occurred after domestication. The identification of multiple, full length FGF4 retrocopies with open reading frames in canids indicates that gene retrotransposition events occur much more frequently than previously thought and provide a mechanism for continued genetic and phenotypic diversity in canids.

Whole-genome genotyping and resequencing reveal the association of a deletion in the complex interferon alpha gene cluster with hypothyroidism in dogs

Background

Hypothyroidism is a common complex endocrinopathy that typically has an autoimmune etiology, and it affects both humans and dogs. Genetic and environmental factors are both known to play important roles in the disease development. In this study, we sought to identify the genetic risk factors potentially involved in the susceptibility to the disease in the high-risk Giant Schnauzer dog breed.

Results

By employing genome-wide association followed by fine-mapping (top variant p-value = 5.7 × 10^− 6), integrated with whole-genome resequencing and copy number variation analysis, we detected a ~ 8.9 kbp deletion strongly associated (p-value = 0.0001) with protection against development of hypothyroidism. The deletion is located between two predicted Interferon alpha (IFNA) genes and it may eliminate functional elements potentially involved in the transcriptional regulation of these genes. Remarkably, type I IFNs have been extensively associated to human autoimmune hypothyroidism and general autoimmunity. Nonetheless, the extreme genomic complexity of the associated region on CFA11 warrants further long-read sequencing and annotation efforts in order to ascribe functions to the identified deletion and to characterize the canine IFNA gene cluster in more detail.

Conclusions

Our results expand the current knowledge on genetic determinants of canine hypothyroidism by revealing a significant link with the human counterpart disease, potentially translating into better diagnostic tools across species, and may contribute to improved canine breeding strategies.

Body size, inbreeding, and lifespan in domestic dogs

Inbreeding poses a real or potential threat to nearly every species of conservation concern. Inbreeding leads to loss of diversity at the individual level, which can cause inbreeding depression, and at the population level, which can hinder ability to respond to a changing environment. In closed populations such as endangered species and ex situ breeding programs, some degree of inbreeding is inevitable. It is therefore vital to understand how different patterns of breeding and inbreeding can affect fitness in real animals. Domestic dogs provide an excellent model, showing dramatic variation in degree of inbreeding and in lifespan, an important aspect of fitness that is known to be impacted by inbreeding in other species. There is a strong negative correlation between body size and lifespan in dogs, but it is unknown whether the higher rate of aging in large dogs is due to body size per se or some other factor associated with large size. We used dense genome-wide SNP array data to calculate average inbreeding for over 100 dog breeds based on autozygous segment length and found that large breeds tend to have higher coefficients of inbreeding than small breeds. We then used data from the Veterinary medical Database and other published sources to estimate life expectancies for pure and mixed breed dogs. When controlling for size, variation in inbreeding was not associated with life expectancy across breeds. When comparing mixed versus purebred dogs, however, mixed breed dogs lived about 1.2 years longer on average than size-matched purebred dogs. Furthermore, individual pedigree coefficients of inbreeding and lifespans for over 9000 golden retrievers showed that inbreeding does negatively impact lifespan at the individual level. Registration data from the American Kennel Club suggest that the molecular inbreeding patterns observed in purebred dogs result from specific breeding practices and/or founder effects and not the current population size. Our results suggest that recent inbreeding, as reflected in variation within a breed, is more likely to affect fitness than historic inbreeding, as reflected in variation among breeds. Our results also indicate that occasional outcrosses, as in mixed breed dogs, can have a substantial positive effect on fitness.

Genome-wide association and evolutionary analyses reveal the formation of swine facial wrinkles in Chinese Erhualian pigs

Wrinkles are uneven concave-convex folds, ridges or creases in skin. Facial wrinkles appear in head, typically increasing along with aging. However in several Chinese indigenous pigs, such as Erhualian pigs, rich facial wrinkles have been generated during the growth stages as one of their breed characteristics. To investigate the genetic basis underlying the development of swine facial wrinkles, we estimated the folding extent of facial wrinkles in a herd of Erhualian pigs (n=332), and then conducted genome-wide association studies and multi-trait meta-analysis for facial wrinkles using 60K porcine chips. We found that facial wrinkles had high heritability estimates of ~0.7 in Erhualian pigs. Notably, only one genome-wide significant QTL was detected at 34.8 Mb on porcine chromosome 7. The most significant SNP rs80983858 located at the 3255-bp downstream of candidate gene GRM4, and the G allele was of benefit to increase facial wrinkles. Evolutionary and selection analyses suggested that the haplotypes containing G allele were under artificial selection, which was consistent with local animal sacrificial custom praying for longevity. Our findings made important clues for further deciphering the molecular mechanism of swine facial wrinkles formation, and shed light on the research of skin wrinkle development in human or other mammals.

Hyaluronan as tunable drug delivery system

The hyaluronan (HA) polymer is an important macromolecule of extracellular matrix with remarkable structure and functions: it is a linear and unbranched polymer without sulphate or phosphate groups and has key role in several biological processes in mammals. It is ubiquitous in mammalian tissues with several and specific functions, influencing cell proliferation and migration as well as angiogenesis and inflammation. To exert these important functions in tissues HA modifies the concentration and size. Considering this HA content in tissues is carefully controlled by different mechanisms including covalent modification of the synthetic enzymes and epigenetic control of their gene expression. The function of HA is also critical in several pathologies including cancer, diabetes and chronic inflammation. Among these biological roles, the structural properties of HA allow to use this polymer in regenerative medicine including cosmetics and drug delivery. HA takes advantage from its capacity to form gels even at concentration of 1% producing scaffolds with very intriguing mechanical properties. These hydrogels are useful in regenerative medicine as biocompatible material for advanced therapeutic uses. In this review we highlight the biological aspects of HA addressing the mechanisms controlling the HA content in tissues and its role as drug delivery system.

Natural models for retinitis pigmentosa: progressive retinal atrophy in dog breeds

Retinitis pigmentosa (RP) is a heterogeneous group of inherited retinal disorders eventually leading to blindness with different ages of onset, progression and severity. Human RP, first characterized by the progressive degeneration of rod photoreceptor cells, shows high genetic heterogeneity with more than 90 genes identified. However, about one-third of patients have no known genetic causes. Interestingly, dogs are also severely affected by similar diseases, called progressive retinal atrophy (PRA). Indeed, RP and PRA have comparable clinical signs, physiopathology and outcomes, similar diagnosis methods and most often, orthologous genes are involved. The many different dog PRAs often segregate in specific breeds. Indeed, undesired alleles have been selected and amplified through drastic selection and excessive use of inbreeding. Out of the 400 breeds, nearly 100 have an inherited form of PRA, which are natural animal models that can be used to investigate the genetics, disease progression and therapies in dogs for the benefit of both dogs and humans. Recent knowledge on the canine genome and access to new genotyping and sequencing technologies now efficiently allows the identification of mutations involved in canine genetic diseases. To date, PRA genes identified in dog breeds correspond to the same genes in humans and represent relevant RP models, and new genes found in dogs represent good candidate for still unknown human RP. We present here a review of the main advantages of the dog models for human RP with the genes already identified and an X-linked PRA in the Border collie as a model for orphan X-linked RPs in human.

Risk-modeling of dog osteosarcoma genome scans shows individuals with Mendelian-level polygenic risk are common

Background

Despite the tremendous therapeutic advances that have stemmed from somatic oncogenetics, survival of some cancers has not improved in 50 years. Osteosarcoma still has a 5-year survival rate of 66%. We propose the natural canine osteosarcoma model can change that: it is extremely similar to the human condition, except for being highly heritable and having a dramatically higher incidence. Here we reanalyze published genome scans of osteosarcoma in three frequently-affected dog breeds and report entirely new understandings with immediate translational indications.

Results

First, meta-analysis revealed association near FGF9, which has strong biological and therapeutic relevance. Secondly, risk-modeling by multiple logistic regression shows 22 of the 34 associated loci contribute to risk and eight have large effect sizes. We validated the Greyhound stepwise model in our own, independent, case-control cohort. Lastly, we updated the gene annotation from approximately 50 genes to 175, and prioritized those using cross-species genomics data. Mostly positional evidence suggests 13 genes are likely to be associated with mapped risk (including MTMR9, EWSR1 retrogene, TANGO2 and FGF9). Previous annotation included seven of those 13 and prioritized four by pathway enrichment. Ten of our 13 priority genes are in loci that contribute to risk modeling and thus can be studied epidemiologically and translationally in pet dogs. Other new candidates include MYCN, SVIL and MIR100HG.

Conclusions

Polygenic osteosarcoma-risk commonly rises to Mendelian-levels in some dog breeds. This justifies caninized animal models and targeted clinical trials in pet dogs (e.g., using CDK4/6 and FGFR1/2 inhibitors).

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5531-6) contains supplementary material, which is available to authorized users.

Copy number variation and variant discovery in Bullmastiff dogs

Identification of genomic variants within dogs is important for understanding genetic factors contributing to breed diversity and phenotypic traits. This study aimed to identify sources of variation in the Bullmastiff using high-density signal intensity and whole-genome sequence data. Close to 3000 copy number variants (CNVs) were identified in Bullmastiff dogs using Canine HD BeadChip data. When CNVs were collated, 82 CNV regions (CNVRs) were detected, 50% in transcribed regions encompassing 432 genes. Fifty of the CNVRs detected have not been reported in other breeds and represent potential breed-specific variants. A proportion of the CNVR variants with predicted modifying effects on gene pathways may contribute to breed traits. Approximately 5 million putative variants per dog, inclusive of single nucleotide polymorphisms (SNPs), multi-nucleotide polymorphisms (MNPs) and insertion and deletions (INDELs), were identified from DNA sequence data on a small number of animals. Identification of genetic variants in the Bullmastiff highlights sources of variation in the breed and molecular markers that will assist in future trait and disease investigations in dogs.

First genome-wide CNV mapping in FELIS CATUS using next generation sequencing data

Background

Copy Number Variations (CNVs) have becoming very significant variants, representing a major source of genomic variation. CNVs involvement in phenotypic expression and different diseases has been widely demonstrated in humans as well as in many domestic animals. However, genome wide investigation on these structural variations is still missing in Felis catus. The present work is the first CNV mapping from a large data set of Next Generation Sequencing (NGS) data in the domestic cat, performed within the 99 Lives Consortium.

Results

Reads have been mapped on the reference assembly_6.2 by Maverix Biomics. CNV detection with cn.MOPS and CNVnator detected 592 CNVs. These CNVs were used to obtain 154 CNV Regions (CNVRs) with BedTools, including 62 singletons. CNVRs covered 0.26% of the total cat genome with 129 losses, 19 gains and 6 complexes. Cluster Analysis and Principal Component Analysis of the detected CNVRs showed that breeds tend to cluster together as well as cats sharing the same geographical origins. The 46 genes identified within the CNVRs were annotated.

Conclusion

This study has improved the genomic characterization of 14 cat breeds and has provided CNVs information that can be used for studies of traits in cats. It can be considered a sound starting point for genomic CNVs identification in this species.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5297-2) contains supplementary material, which is available to authorized users.

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.

Evidence of evolutionary history and selective sweeps in the genome of Meishan pig reveals its genetic and phenotypic characterization

Background

Meishan is a pig breed indigenous to China and famous for its high fecundity. The traits of Meishan are strongly associated with its distinct evolutionary history and domestication. However, the genomic evidence linking the domestication of Meishan pigs with its unique features is still poorly understood. The goal of this study is to investigate the genomic signatures and evolutionary evidence related to the phenotypic traits of Meishan via large-scale sequencing.

Results

We found that the unique domestication of Meishan pigs occurred in the Taihu Basin area between the Majiabang and Liangzhu Cultures, during which 300 protein-coding genes have underwent positive selection. Notably, enrichment of the FoxO signaling pathway with significant enrichment signal and the harbored gene IGF1R were likely associated with the high fertility of Meishan pigs. Moreover, NFKB1 exhibited strong selective sweep signals and positively participated in hyaluronan biosynthesis as the key gene of NF-kB signaling, which may have resulted in the wrinkled skin and face of Meishan pigs. Particularly, three population-specific synonymous single-nucleotide variants occurred in PYROXD1, MC1R, and FAM83G genes; the T305C substitution in the MCIR gene explained the black coat of the Meishan pigs well. In addition, the shared haplotypes between Meishan and Duroc breeds confirmed the previous Asian-derived introgression and demonstrated the specific contribution of Meishan pigs.

Conclusions

These findings will help us explain the unique genetic and phenotypic characteristics of Meishan pigs and offer a plausible method for their utilization of Meishan pigs as valuable genetic resources in pig breeding and as an animal model for human wrinkled skin disease research.

Similar genomic proportions of copy number variation within gray wolves and modern dog breeds inferred from whole genome sequencing

BACKGROUND

Whole genome re-sequencing data from dogs and wolves are now commonly used to study how natural and artificial selection have shaped the patterns of genetic diversity. Single nucleotide polymorphisms, microsatellites and variants in mitochondrial DNA have been interrogated for links to specific phenotypes or signals of domestication. However, copy number variation (CNV), despite its increasingly recognized importance as a contributor to phenotypic diversity, has not been extensively explored in canids.

RESULTS

Here, we develop a new accurate probabilistic framework to create fine-scale genomic maps of segmental duplications (SDs), compare patterns of CNV across groups and investigate their role in the evolution of the domestic dog by using information from 34 canine genomes. Our analyses show that duplicated regions are enriched in genes and hence likely possess functional importance. We identify 86 loci with large CNV differences between dogs and wolves, enriched in genes responsible for sensory perception, immune response, metabolic processes, etc. In striking contrast to the observed loss of nucleotide diversity in domestic dogs following the population bottlenecks that occurred during domestication and breed creation, we find a similar proportion of CNV loci in dogs and wolves, suggesting that other dynamics are acting to particularly select for CNVs with potentially functional impacts.

CONCLUSIONS

This work is the first comparison of genome wide CNV patterns in domestic and wild canids using whole-genome sequencing data and our findings contribute to study the impact of novel kinds of genetic changes on the evolution of the domestic dog.

The deubiquitinating enzymes USP4 and USP17 target hyaluronan synthase 2 and differentially affect its function

The levels of hyaluronan, a ubiquitous glycosaminoglycan prominent in the extracellular matrix, is balanced through the actions of hyaluronan-synthesizing enzymes (HAS1, 2 and 3) and degrading hyaluronidases (Hyal 1, 2, 3 and PH20). Hyaluronan accumulates in rapidly remodeling tissues, such as breast cancer, due to deregulated expression of the HAS2 gene and/or alterations of HAS2 activity. The activity of HAS2 is regulated by post-translational modifications, including ubiquitination. In order to identify deubiquitinating enzymes (DUBs) that are involved in de-ubiquitination of HAS2, a complementary (cDNA) library of 69 Flag-HA-tagged human DUBs cloned into retroviral vectors was screened in human embryonic kidney (HEK) 293T cells for their ability to de-ubiquitinate myc-tagged HAS2. Several DUBs were found to decrease the ubiquitination of 6myc-HAS2, among which, the most effective were USP17 and USP4. USP17 efficiently removed polyubiquitination, whereas USP4 preferentially removed monoubiquitination of 6myc-HAS2. Co-immunoprecipitation studies revealed interactions between HAS2 and USP17, as well as between HAS2 and USP4, in membrane preparations of HEK293T cells. USP17 significantly stabilized 6myc-HAS2 protein levels, whereas USP4 did not. The silencing of USP17 led to decreased hyaluronan production, whereas the suppression of USP4 increased hyaluronan synthesis. Importantly, high levels of USP17 and HAS2 were detected in a panel of cancer cell lines compared to normal cells, and immunohistochemical stainings revealed higher expression of USP17 and HAS2 in tissues of lung cancer patients compared to normal tissue. In conclusion, USP17 and USP4 differently affect HAS2 ubiquitination, and the stability and function of HAS2.

Characterization of Copy Number Variation's Potential Role in Marek's Disease

Marek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.

Whole genome sequencing identifies missense mutation in MTBP in Shar-Pei affected with Autoinflammatory Disease (SPAID)

BACKGROUND

Autoinflammatory diseases in dogs are characterized by complex disease processes with varying clinical signs. In Shar-Pei, signs of inflammation including fever and arthritis are known to be related with a breed-specific predisposition for Shar-Pei Autoinflammatory Disease (SPAID).

RESULTS

Clinical and histopathological examinations of two severely SPAID-affected Shar-Pei revealed signs of inflammation including fever, arthritis, and perivascular and diffuse dermatitis in both dogs. A multifocal accumulation of amyloid in different organs was found in one SPAID-affected case. Whole genome sequencing resulted in 37 variants, which were homozygous mutant private mutations in SPAID-affected Shar-Pei. Nine SNVs with predicted damaging effects and three INDELs were further investigated in 102 Shar-Pei affected with SPAID, 62 unaffected Shar-Pei and 162 controls from 11 different dog breeds. The results showed the missense variant MTBP:g.19383758G > A in MTBP to be highly associated with SPAID in Shar-Pei. In the region of this gene a large ROH (runs of homozygosity) region could be detected exclusively in the two investigated SPAID-affected Shar-Pei compared to control dog breeds. No further SPAID-associated variant with predicted high or moderate effects could be found in genes identified in ROH regions. This MTBP variant was predicted to affect the MDN2-binding protein domain and consequently promote proinflammatory reactions. In the investigated group of Shar-Pei older than six years all dogs with the mutant genotype A/A were SPAID-affected whereas SPAID-unaffected dogs harbored the homozygous wildtype (G/G). Shar-Pei with a heterozygous genotype (G/A) were shown to have a 2.13-fold higher risk for disease development, which gave evidence for an incomplete dominant mode of inheritance.

CONCLUSIONS

The results of this study give strong evidence for a variant in MTBP related with proinflammatory processes via MTBP-MDM2 pathway. Thus, these results enable a reliable detection of SPAID in Shar-Pei dogs.

Evaluation of hyaluronic acid, procollagen type III N-terminal peptide, and tissue inhibitor of matrix metalloproteinase-1 as serum markers of canine hepatic fibrosis

The only way to diagnose hepatic fibrosis in dogs is by histological assessment of a liver biopsy specimen. As this technique is invasive and susceptible to sampling variation, serum biomarkers are used to detect hepatic fibrosis in humans. The objective of this study was to assess the utility of hyaluronic acid (HA), procollagen type III N-terminal peptide (PIIINP), and tissue inhibitor of matrix metalloproteinase-1 (TIMP-1) as serum markers of canine hepatic fibrosis. Serum samples were collected from 47 dogs with histologically confirmed hepatobiliary disease and 24 healthy dogs in order to measure concentrations of HA, PIIINP, and TIMP-1. Hepatic fibrosis was staged using a 5-point scoring scheme. There was no correlation between serum concentrations of HA or PIIINP and the severity of hepatic fibrosis. There was a negative correlation between serum concentration of TIMP-1 and the severity of hepatic fibrosis (rs = -0.33; P = 0.036). It was not possible to use serum concentrations of HA, PIIINP, or TIMP-1 to discriminate between dogs with absent-to-moderate hepatic fibrosis and those with marked-to-very-marked fibrosis. The results of this study do not support the utility of measuring serum concentrations of HA, PIIINP, or TIMP-1 for diagnosing canine hepatic fibrosis. Further studies are needed to support this finding.

Adverse interaction between colchicine and ketoconazole in a Chinese shar pei

A Chinese shar pei with a 2 yr history of episodic fever, lethargy, and shifting lameness was presumptively diagnosed with familial shar pei fever but had never been treated for the syndrome. After being presented for a superficial pyoderma with possible dermatophyte coinfection, treatment with a cephalosporin and ketoconazole were prescribed. One wk later, colchicine was initiated for familial shar pei fever using cautious dose escalation. Nevertheless, gastrointestinal toxicity, skeletal muscle myopathy, and hepatotoxicity developed within 2 wk. Abrupt resolution of gastrointestinal toxicity and myopathy followed drug withdrawal. However, escalating liver enzyme activity and hyperbilirubinemia led to liver biopsy to rule out an antecedent hepatopathy. Biopsy characterized canalicular cholestasis and colchicine-associated metaphase arrest and ring mitoses reflecting repression of mitotic spindle formation. Signs of illness completely resolved 3 mo after drug discontinuation. Although avoidable adverse interactions between ketoconazole and drugs reliant on cytochrome oxidase biotransformation and/or drug efflux mediated by multiple drug-resistant transporters are well documented in humans, these are rarely reported in veterinary patients. This case exemplifies an important and avoidable ketoconazole/colchicine drug interaction from which the patient completely recovered. The dog tested negative for the canine MDR1 loss of function mutation that also might potentiate colchicine toxicity.

Genetic testing in veterinary dermatology

BACKGROUND

Molecular genetics has made significant advances in the analysis of hereditary dermatoses during the last several years.

OBJECTIVES

To provide an update on currently available genetic tests for skin diseases of dogs, cats and horses, and to aid the veterinary clinician in the appropriate selection and applications of genetic tests.

METHODS

The scientific literature on the topic was critically reviewed. The list of known causative variants for genodermatoses and hair morphology traits was compiled by searching the Online Mendelian Inheritance in Animals (OMIA) database.

RESULTS

Genetic testing has become an important diagnostic method in veterinary medicine. Genetic tests can help to establish the correct diagnosis in some diseases with relatively nonspecific signs. Genetic tests are also essential for sustainable breeding programmes and to help minimize the frequency of animals with hereditary diseases. Advances in genetic methodology and bioinformatics already allow genome-wide screening for potential disease causing mutations for research purposes. It is anticipated that this will become a routine process in clinical practice in the future.

CONCLUSION AND CLINICAL IMPORTANCE

As specific DNA tests and broad genome-wide analyses come into more common use, it is critical that clinicians understand the proper application and interpretation of these test results.

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.

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.

Genome-Wide Analysis Identifies Germ-Line Risk Factors Associated with Canine Mammary Tumours

Canine mammary tumours (CMT) are the most common neoplasia in unspayed female dogs. CMTs are suitable naturally occurring models for human breast cancer and share many characteristics, indicating that the genetic causes could also be shared. We have performed a genome-wide association study (GWAS) in English Springer Spaniel dogs and identified a genome-wide significant locus on chromosome 11 (p_raw = 5.6x10^-7, p_perm = 0.019). The most associated haplotype spans a 446 kb region overlapping the CDK5RAP2 gene. The CDK5RAP2 protein has a function in cell cycle regulation and could potentially have an impact on response to chemotherapy treatment. Two additional loci, both on chromosome 27, were nominally associated (p_raw = 1.97x10^-5 and p_raw = 8.30x10^-6). The three loci explain 28.1±10.0% of the phenotypic variation seen in the cohort, whereas the top ten associated regions account for 38.2±10.8% of the risk. Furthermore, the ten GWAS loci and regions with reduced genetic variability are significantly enriched for snoRNAs and tumour-associated antigen genes, suggesting a role for these genes in CMT development. We have identified several candidate genes associated with canine mammary tumours, including CDK5RAP2. Our findings enable further comparative studies to investigate the genes and pathways in human breast cancer patients.

Dogs provide an excellent model system for several human diseases, including cancer. Heavy breeding for certain behavioural or phenotypic traits has created genetic isolates–breeds–characterised by low levels of genetic variation and a limited number of genetic disease variants within each breed. Cancer is the most common cause of death in dogs today, and canine mammary tumours (CMT) are the most prevalent tumour type in unspayed female dogs. These tumours are very similar to human breast cancer and could therefore be used as a naturally occurring model for the human disease. We have investigated genetic variants associated with CMT in English Springer Spaniels pointing to a gene involved in cell cycle regulation (CDK5RAP2). The CDK5RAP2 could therefore have a key role in the development of mammary tumours and we suggest that further studies should be performed in both dogs and women to investigate CDK5RAP2 and its possible effect on disease and treatment response.

An Improved microRNA Annotation of the Canine Genome

The domestic dog, Canis familiaris, is a valuable model for studying human diseases. The publication of the latest Canine genome build and annotation, CanFam3.1 provides an opportunity to enhance our understanding of gene regulation across tissues in the dog model system. In this study, we used the latest dog genome assembly and small RNA sequencing data from 9 different dog tissues to predict novel miRNAs in the dog genome, as well as to annotate conserved miRNAs from the miRBase database that were missing from the current dog annotation. We used both miRCat and miRDeep2 algorithms to computationally predict miRNA loci. The resulting, putative hairpin sequences were analysed in order to discard false positives, based on predicted secondary structures and patterns of small RNA read alignments. Results were further divided into high and low confidence miRNAs, using the same criteria. We generated tissue specific expression profiles for the resulting set of 811 loci: 720 conserved miRNAs, (207 of which had not been previously annotated in the dog genome) and 91 novel miRNA loci. Comparative analyses revealed 8 putative homologues of some novel miRNA in ferret, and one in microbat. All miRNAs were also classified into the genic and intergenic categories, based on the Ensembl RefSeq gene annotation for CanFam3.1. This additionally allowed us to identify four previously undescribed MiRtrons among our total set of miRNAs. We additionally annotated piRNAs, using proTRAC on the same input data. We thus identified 263 putative clusters, most of which (211 clusters) were found to be expressed in testis. Our results represent an important improvement of the dog genome annotation, paving the way to further research on the evolution of gene regulation, as well as on the contribution of post-transcriptional regulation to pathological conditions.

Absolute quantification reveals the stable transmission of a high copy number variant linked to autoinflammatory disease

BACKGROUND

Dissecting the role copy number variants (CNVs) play in disease pathogenesis is directly reliant on accurate methods for quantification. The Shar-Pei dog breed is predisposed to a complex autoinflammatory disease with numerous clinical manifestations. One such sign, recurrent fever, was previously shown to be significantly associated with a novel, but unstable CNV (CNV_16.1). Droplet digital PCR (ddPCR) offers a new mechanism for CNV detection via absolute quantification with the promise of added precision and reliability. The aim of this study was to evaluate ddPCR in relation to quantitative PCR (qPCR) and to assess the suitability of the favoured method as a genetic test for Shar-Pei Autoinflammatory Disease (SPAID).

RESULTS

One hundred and ninety-six individuals were assayed using both PCR methods at two CNV positions (CNV_14.3 and CNV_16.1). The digital method revealed a striking result. The CNVs did not follow a continuum of alleles as previously reported, rather the alleles were stable and pedigree analysis showed they adhered to Mendelian segregation. Subsequent analysis of ddPCR case/control data confirmed that both CNVs remained significantly associated with the subphenotype of fever, but also to the encompassing SPAID complex (p < 0.001). In addition, harbouring CNV_16.1 allele five (CNV_16.1|5) resulted in a four-fold increase in the odds for SPAID (p < 0.001). The inclusion of a genetic marker for CNV_16.1 in a genome-wide association test revealed that this variant explained 9.7 % of genetic variance and 25.8 % of the additive genetic heritability of this autoinflammatory disease.

CONCLUSIONS

This data shows the utility of the ddPCR method to resolve cryptic copy number inheritance patterns and so open avenues of genetic testing. In its current form, the ddPCR test presented here could be used in canine breeding to reduce the number of homozygote CNV_16.1|5 individuals and thereby to reduce the prevalence of disease in this breed.

Th17 Pathway As a Target for Multipotent Stromal Cell Therapy in Dogs: Implications for Translational Research

Detrimental Th17 driven inflammatory and autoimmune disease such as Crohn's disease, graft versus host disease and multiple sclerosis remain a significant cause of morbidity and mortality worldwide. Multipotent stromal/stem cell (MSC) inhibit Th17 polarization and activation in vitro and in rodent models. As such, MSC based therapeutic approaches are being investigated as novel therapeutic approaches to treat Th17 driven diseases in humans. The significance of naturally occurring diseases in dogs is increasingly recognized as a realistic platform to conduct pre-clinical testing of novel therapeutics. Full characterization of Th17 cells in dogs has not been completed. We have developed and validated a flow-cytometric method to detect Th17 cells in canine blood. We further demonstrate that Th17 and other IL17 producing cells are present in tissues of dogs with naturally occurring chronic inflammatory diseases. Finally, we have determined the kinetics of a canine specific Th17 polarization in vitro and demonstrate that canine MSC inhibit Th17 polarization in vitro, in a PGE2 independent mechanism. Our findings provide fundamental research tools and suggest that naturally occurring diseases in dogs, such as inflammatory bowel disease, may be harnessed to translate novel MSC based therapeutic strategies that target the Th17 pathway.

Genome-Wide Association Study of Golden Retrievers Identifies Germ-Line Risk Factors Predisposing to Mast Cell Tumours

Canine mast cell tumours (CMCT) are one of the most common skin tumours in dogs with a major impact on canine health. Certain breeds have a higher risk of developing mast cell tumours, suggesting that underlying predisposing germ-line genetic factors play a role in the development of this disease. The genetic risk factors are largely unknown, although somatic mutations in the oncogene C-KIT have been detected in a proportion of CMCT, making CMCT a comparative model for mastocytosis in humans where C-KIT mutations are frequent. We have performed a genome wide association study in golden retrievers from two continents and identified separate regions in the genome associated with risk of CMCT in the two populations. Sequence capture of associated regions and subsequent fine mapping in a larger cohort of dogs identified a SNP associated with development of CMCT in the GNAI2 gene (p = 2.2x10-16), introducing an alternative splice form of this gene resulting in a truncated protein. In addition, disease associated haplotypes harbouring the hyaluronidase genes HYAL1, HYAL2 and HYAL3 on cfa20 and HYAL4, SPAM1 and HYALP1 on cfa14 were identified as separate risk factors in European and US golden retrievers, respectively, suggesting that turnover of hyaluronan plays an important role in the development of CMCT.