Epidermolysis bullosa in Danish Hereford calves is caused by a deletion in LAMC2 gene.

BACKGROUND Heritable forms of epidermolysis bullosa (EB) constitute a heterogeneous group of skin disorders of genetic aetiology that are characterised by skin and mucous membrane blistering and ulceration in response to even minor trauma. Here we report the occurrence of EB in three Danish Hereford cattle from one herd. RESULTS Two of the animals were necropsied and showed oral mucosal blistering, skin ulcerations and partly loss of horn on the claws. Lesions were histologically characterized by subepidermal blisters and ulcers. Analysis of the family tree indicated that inbreeding and the transmission of a single recessive mutation from a common ancestor could be causative. We performed whole genome sequencing of one affected calf and searched all coding DNA variants. Thereby, we detected a homozygous 2.4 kb deletion encompassing the first exon of the LAMC2 gene, encoding for laminin gamma 2 protein. This loss of function mutation completely removes the start codon of this gene and is therefore predicted to be completely disruptive. The deletion co-segregates with the EB phenotype in the family and absent in normal cattle of various breeds. Verifying the homozygous private variants present in candidate genes allowed us to quickly identify the causative mutation and contribute to the final diagnosis of junctional EB in Hereford cattle. CONCLUSIONS Our investigation confirms the known role of laminin gamma 2 in EB aetiology and shows the importance of whole genome sequencing in the analysis of rare diseases in livestock.

DNA-based diagnosis of rare diseases in veterinary medicine: a 4.4 kb deletion of ITGB4 is associated with epidermolysis bullosa in Charolais cattle.

BACKGROUND Rare diseases in livestock animals are traditionally poorly diagnosed. Other than clinical description and pathological examination, the underlying causes have, for the most part, remained unknown. A single case of congenital skin fragility in cattle was observed, necropsy, histological and ultrastructural examinations were carried out and whole genome sequencing was utilized to identify the causative mutation. RESULTS A single purebred female Charolais calf with severe skin lesions was delivered full-term and died spontaneously after birth. The clinical and pathological findings exactly matched the gross description given by previous reports on epitheliogenesis imperfecta and epidermolysis bullosa (EB) in cattle. Histological and ultrastructural changes were consistent with EB junctionalis (EBJ). Genetic analysis revealed a previously unpublished ITGB4 loss-of-function mutation; the affected calf was homozygous for a 4.4 kb deletion involving exons 17 to 22, and the dam carried a single copy of the deletion indicating recessive inheritance. The homozygous mutant genotype did not occur in healthy controls of various breeds but some heterozygous carriers were found among Charolais cattle belonging to the affected herd. The mutant allele was absent in a representative sample of unrelated sires of the German Charolais population. CONCLUSION This is the first time in which a recessively inherited ITGB4 associated EBJ has been reported in cattle. The identification of heterozygous carriers is of importance in avoiding the transmission of this defect in future. Current DNA sequencing methods offer a powerful tool for understanding the genetic background of rare diseases in domestic animals having a reference genome sequence available.

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency.

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.

Molecular characterization and SNP development for the porcine IL6 and IL10 genes

Different cytokines are secreted in response to specific microbial molecules referred to as pathogen associated molecular patterns (PAMPs). Interleukin 6 (IL6) and interleukin 10 (IL10), both secreted by macrophages and lymphocytes, play a central role in the immunological response. In this work we obtained the genomic structure and complete DNA sequence of the porcine IL6 and IL10 genes and identified polymorphisms in the genomic sequences of these genes on a panel of ten different pig breeds. Comparative intra- and interbreed sequence analysis revealed a total of eight polymorphisms in the porcine IL6 gene and 21 in the porcine IL10 gene, which include single nucleotide polymorphisms (SNPs) and insertion deletion polymorphisms (indels). Additionally, the chromosomal localization of the IL10 gene was determined by FISH and RH mapping.

Allele-specific polymerase chain reaction diagnostic test for the functional MDR1 polymorphism in dogs

The major multidrug transporter P-glycoprotein (Pgp) contributes to the barrier function of several tissues and organs, including the brain. In a subpopulation of Collies and seven further dog breeds, a 4 base pair deletion has been described in the Pgp-encoding MDR1 gene. This deletion results in the absence of a functional form of Pgp and loss of its protective function. Severe intoxication with the Pgp substrate ivermectin has been attributed to the genetically determined lack of Pgp. An allele-specific polymerase chain reaction (PCR)-based screening method has been developed to detect the mutant allele and to determine if a dog is homozygous or heterozygous for the mutation. Based on this validation, the allele-specific PCR proved to be a robust, reproducible and specific tool, allowing rapid determination of the MDR1 genotype of dogs of at risk breeds using blood samples or buccal swabs.

Polymorphisms within the canine MLPH gene are associated with dilute coat color in dogs

BACKGROUND Pinschers and other dogs with coat color dilution show a characteristic pigmentation phenotype. The fur colors are a lighter shade, e.g. silvery grey (blue) instead of black and a sandy color (Isabella fawn) instead of red or brown. In some dogs the coat color dilution is sometimes accompanied by hair loss and recurrent skin inflammation, the so called color dilution alopecia (CDA) or black hair follicular dysplasia (BHFD). In humans and mice a comparable pigmentation phenotype without any documented hair loss is caused by mutations within the melanophilin gene (MLPH). RESULTS We sequenced the canine MLPH gene and performed a mutation analysis of the MLPH exons in 6 Doberman Pinschers and 5 German Pinschers. A total of 48 sequence variations was identified within and between the breeds. Three families of dogs showed co-segregation for at least one polymorphism in an MLPH exon and the dilute phenotype. No single polymorphism was identified in the coding sequences or at splice sites that is likely to be causative for the dilute phenotype of all dogs examined. In 18 German Pinschers a mutation in exon 7 (R199H) was consistently associated with the dilute phenotype. However, as this mutation was present in homozygous state in four dogs of other breeds with wildtype pigmentation, it seems unlikely that this mutation is truly causative for coat color dilution. In Doberman Pinschers as well as in Large Munsterlanders with BHFD, a set of single nucleotide polymorphisms (SNPs) around exon 2 was identified that show a highly significant association to the dilute phenotype. CONCLUSION This study provides evidence that coat color dilution is caused by one or more mutations within or near the MLPH gene in several dog breeds. The data on polymorphisms that are strongly associated with the dilute phenotype will allow the genetic testing of Pinschers to facilitate the breeding of dogs with defined coat colors and to select against Large Munsterlanders carrying BHFD.

Molecular characterization and chromosomal assignment of the bovine glycinamide ribonucleotide formyltransferase (GART) gene on cattle chromosome 1q12.1-q12.2

The mammalian glycinamide ribonucleotide formyltransferase (GART) genes encode a trifunctional polypeptide involved in the de novo purine biosynthesis. We isolated a bacterial artificial chromosome (BAC) clone containing the bovine GART gene and determined the complete DNA sequence of the BAC clone. Cloning and characterization of the bovine GART gene revealed that the bovine gene consists of 23 exons spanning approximately 27 kb. RT-PCR amplification of bovine GART in different organs showed the expression of two GART transcripts in cattle similar to human and mouse. The GART transcripts encode two proteins of 1010 and 433 amino acids, respectively. Eleven single nucleotide polymorphisms (SNPs) were detected in a mutation scan of 24 unrelated animals of three different cattle breeds, including one SNP that affects the amino acid sequence of GART. The chromosomal localization of the gene was determined by fluorescence in situ hybridization. Comparative genome analysis between cattle, human and mouse indicates that the chromosomal location of the bovine GART gene is in agreement with a previously published mapping report.

A 1.8-kb insertion in the 3'-UTR of RXFP2 is associated with polledness in sheep

Sheep breeds show a broad spectrum of different horn phenotypes. In most modern production breeds, sheep are polled (absence of horns), whereas horns occur mainly in indigenous breeds. Previous studies mapped the responsible locus to the region of the RXFP2 gene on ovine chromosome 10. A 4-kb region of the 3'-end of RXFP2 was amplified in horned and polled animals from seven Swiss sheep breeds. Sequence analysis identified a 1833-bp genomic insertion located in the 3'-UTR region of RXFP2 present in polled animals only. An efficient PCR-based genotyping method to determine the polled genotype of individual sheep is presented. Comparative sequence analyses revealed evidence that the polled-associated insertion adds a potential antisense RNA sequence of EEF1A1 to the 3'-end of RXFP2 transcripts.

Fine-mapping and mutation analysis of TRPM1: a candidate gene for leopard complex (LP) spotting and congenital stationary night blindness in horses.

Leopard Complex spotting occurs in several breeds of horses and is caused by an incompletely dominant allele (LP). Homozygosity for LP is also associated with congenital stationary night blindness (CSNB) in Appaloosa horses. Previously, LP was mapped to a 6 cm region on ECA1 containing the candidate gene TRPM1 (Transient Receptor Potential Cation Channel, Subfamily M, Member 1) and decreased expression of this gene, measured by qRT-PCR, was identified as the likely cause of both spotting and ocular phenotypes. This study describes investigations for a mutation causing or associated with the Leopard Complex and CSNB phenotype in horses. Re-sequencing of the gene and associated splice sites within the 105 624 bp genomic region of TRPM1 led to the discovery of 18 SNPs. Most of the SNPs did not have a predictive value for the presence of LP. However, one SNP (ECA1:108,249,293 C>T) found within intron 11 had a strong (P < 0.0005), but not complete, association with LP and CSNB and thus is a good marker but unlikely to be causative. To further localize the association, 70 SNPs spanning over two Mb including the TRPM1 gene were genotyped in 192 horses from three different breeds segregating for LP. A single 173 kb haplotype associated with LP and CSNB (ECA1: 108,197,355- 108,370,150) was identified. Illumina sequencing of 300 kb surrounding this haplotype revealed 57 SNP variants. Based on their localization within expressed sequences or regions of high sequence conservation across mammals, six of these SNPs were considered to be the most likely candidate mutations. While the precise function of TRPM1 remains to be elucidated, this work solidifies its functional role in both pigmentation and night vision. Further, this work has identified several potential regulatory elements of the TRPM1 gene that should be investigated further in this and other species.

A Multi-Breed Genome-Wide Association Analysis for Canine Hypothyroidism Identifies a Shared Major Risk Locus on CFA12.

Hypothyroidism is a complex clinical condition found in both humans and dogs, thought to be caused by a combination of genetic and environmental factors. In this study we present a multi-breed analysis of predisposing genetic risk factors for hypothyroidism in dogs using three high-risk breeds-the Gordon Setter, Hovawart and the Rhodesian Ridgeback. Using a genome-wide association approach and meta-analysis, we identified a major hypothyroidism risk locus shared by these breeds on chromosome 12 (p = 2.1x10-11). Further characterisation of the candidate region revealed a shared ~167 kb risk haplotype (4,915,018-5,081,823 bp), tagged by two SNPs in almost complete linkage disequilibrium. This breed-shared risk haplotype includes three genes (LHFPL5, SRPK1 and SLC26A8) and does not extend to the dog leukocyte antigen (DLA) class II gene cluster located in the vicinity. These three genes have not been identified as candidate genes for hypothyroid disease previously, but have functions that could potentially contribute to the development of the disease. Our results implicate the potential involvement of novel genes and pathways for the development of canine hypothyroidism, raising new possibilities for screening, breeding programmes and treatments in dogs. This study may also contribute to our understanding of the genetic etiology of human hypothyroid disease, which is one of the most common endocrine disorders in humans.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication

The domestication of the horse revolutionized warfare, trade, and the exchange of people and ideas. This at least 5,500-y-long process, which ultimately transformed wild horses into the hundreds of breeds living today, is difficult to reconstruct from archeological data and modern genetics alone. We therefore sequenced two complete horse genomes, predating domestication by thousands of years, to characterize the genetic footprint of domestication. These ancient genomes reveal predomestic population structure and a significant fraction of genetic variation shared with the domestic breeds but absent from Przewalski’s horses. We find positive selection on genes involved in various aspects of locomotion, physiology, and cognition. Finally, we show that modern horse genomes contain an excess of deleterious mutations, likely representing the genetic cost of domestication.

A Nonsense Variant in COL6A1 in Landseer Dogs with Muscular Dystrophy.

A novel canine muscular dystrophy in Landseer dogs was observed. We had access to five affected dogs from two litters. The clinical signs started at a few weeks of age and the severe progressive muscle weakness led to euthanasia between 5 and 15 months of age. The pedigrees of the affected dogs suggested a monogenic autosomal recessive inheritance of the trait. Linkage and homozygosity mapping indicated two potential genome segments for the causative variant on chromosomes 10 and 31 harboring a total of 4.8 Mb of DNA or 0.2% of the canine genome. Using the illumina sequencing technology we obtained a whole genome sequence from one affected Landseer. Variants were called with respect to the dog reference genome and compared to the genetic variants of 170 control dogs from other breeds. The affected Landseer dog was homozygous for a single private non-synonymous variant in the critical intervals, a nonsense variant in the COL6A1 gene (Chr31:39,303,964G>T; COL6A1:c.289G>T; p.E97*). Genotypes at this variant showed perfect concordance with the muscular dystrophy phenotype in all five cases and more than one thousand control dogs. Variants in the human COL6A1 gene cause Bethlem myopathy or Ullrich congenital muscular dystrophy. We therefore conclude that the identified canine COL6A1 variant is most likely causative for the observed muscular dystrophy in Landseer dogs. Based on the nature of the genetic variant in Landseer dogs and their severe clinical phenotype these dogs represent a model for human Ullrich congenital muscular dystrophy.

A Non-Synonymous HMGA2 Variant Decreases Height in Shetland Ponies and Other Small Horses

The identification of quantitative trait loci (QTL) such as height and their underlying causative variants is still challenging and often requires large sample sizes. In humans hundreds of loci with small effects control the heritable portion of height variability. In domestic animals, typically only a few loci with comparatively large effects explain a major fraction of the heritability. We investigated height at withers in Shetland ponies and mapped a QTL to ECA 6 by genome-wide association (GWAS) using a small cohort of only 48 animals and the Illumina equine SNP70 BeadChip. Fine-mapping revealed a shared haplotype block of 793 kb in small Shetland ponies. The HMGA2 gene, known to be associated with height in horses and many other species, was located in the associated haplotype. After closing a gap in the equine reference genome we identified a non-synonymous variant in the first exon of HMGA2 in small Shetland ponies. The variant was predicted to affect the functionally important first AT-hook DNA binding domain of the HMGA2 protein (c.83G>A; p.G28E). We assessed the functional impact and found impaired DNA binding of a peptide with the mutant sequence in an electrophoretic mobility shift assay. This suggests that the HMGA2 variant also affects DNA binding in vivo and thus leads to reduced growth and a smaller stature in Shetland ponies. The identified HMGA2 variant also segregates in several other pony breeds but was not found in regular-sized horse breeds. We therefore conclude that we identified a quantitative trait nucleotide for height in horses.

Determination of vertebral range of motion using inertial measurement units in 27 Franches-Montagnes stallions and comparison between conditions and with a mixed population

REASONS FOR PERFORMING STUDY: The diagnosis of equine back disorders is challenging. Objectively determining movement of the vertebral column may therefore be of value in a clinical setting. OBJECTIVES: To establish whether surface-mounted inertial measurement units (IMUs) can be used to establish normal values for range of motion (ROM) of the vertebral column in a uniform population of horses trotting under different conditions. STUDY DESIGN: Vertebral ROM was established in Franches-Montagnes stallions and a general population of horses and the variability in measurements compared between the two groups. Repeatability and the influence of specific exercise condition (on ROM) were assessed. Finally, attempts were made to explain the findings of the study through the evaluation of factors that might influence ROM. METHODS: Dorsoventral (DV) and mediolateral (ML) vertebral ROM was measured at a trot under different exercise conditions in 27 Franches-Montagnes stallions and six general population horses using IMUs distributed over the vertebral column. RESULTS: Variability in the ROM measurements was significantly higher for general population horses than for Franches-Montagnes stallions (both DV and ML ROM). Repeatability was strong to very strong for DV measurements and moderate for ML measurements. Trotting under saddle significantly reduced the ROM, with sitting trot resulting in a significantly lower ROM than rising trot. Age is unlikely to explain the low variability in vertebral ROM recorded in the Franches-Montagnes horses, while this may be associated with conformational factors. CONCLUSIONS: It was possible to establish a normal vertebral ROM for a group of Franches-Montagnes stallions. While within-breed variation was low in this population, further studies are necessary to determine variation in vertebral ROM for other breeds and to assess their utility for diagnosis of equine back disorders.

A RAB3GAP1 SINE Insertion in Alaskan Huskies with Polyneuropathy, Ocular Abnormalities and Neuronal Vacuolation (POANV) Resembling Human Warburg Micro Syndrome 1 (WARBM1)

We observed a hereditary phenotype in Alaskan Huskies, which was characterized by polyneuropathy with ocular abnormalities and neuronal vacuolation (POANV). The affected dogs developed a progressive severe ataxia, which led to euthanasia between 8 and 16 months of age. The pedigrees were consistent with a monogenic autosomal recessive inheritance. We localized the causative genetic defect to a 4 Mb interval on chromosome 19 by a combined linkage and homozygosity mapping approach. Whole genome sequencing of one affected dog, an obligate carrier and an unrelated control revealed a 218 bp SINE insertion into exon 7 of the RAB3GAP1 gene. The SINE insertion was perfectly associated with the disease phenotype in a cohort of 43 Alaskan Huskies and it was absent from 541 control dogs of diverse other breeds. The SINE insertion induced aberrant splicing and led to a transcript with a greatly altered exon 7. RAB3GAP1 loss-of-function variants in humans cause Warburg Micro Syndrome 1 (WARBM1), which is characterized by additional developmental defects compared to canine POANV, whereas Rab3gap1 deficient mice have a much milder phenotype than either humans or dogs. Thus the RAB3GAP1 mutant Alaskan Huskies provide an interesting intermediate phenotype that may help to better understand the function of RAB3GAP1 in development. Furthermore, the identification of the presumed causative genetic variant will enable genetic testing to avoid the non-intentional breeding of affected dogs.