Congenital hepatic fibrosis in the Franches-Montagnes horse is associated with the polycystic kidney and hepatic disease 1 (PKHD1) gene

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed

BACKGROUND

The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity.

RESULTS

To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders.

CONCLUSIONS

The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Does inbreeding contribute to pregnancy loss in Thoroughbred horses?

BACKGROUND

Excessive inbreeding increases the probability of uncovering homozygous recessive genotypes and has been associated with an increased risk of retained placenta and lower semen quality. No genomic analysis has investigated the association between inbreeding levels and pregnancy loss.

OBJECTIVES

To compare genetic inbreeding coefficients (F) of naturally occurring Thoroughbred Early Pregnancy Loss (EPLs), Mid and Late term Pregnancy Loss (MLPL) and Controls. The F value was hypothesised to be higher in cases of pregnancy loss (EPLs and MLPLs) than Controls.

STUDY DESIGN

Observational case-control study.

METHODS

Allantochorion and fetal DNA from EPL (n = 37, gestation age 14-65 days), MLPL (n = 94, gestational age 70 days-24 h post parturition) and Controls (n = 58) were genotyped on the Axiom Equine 670K SNP Genotyping Array. Inbreeding coefficients using Runs of Homozygosity (FROH) were calculated using PLINK software. ROHs were split into size categories to investigate the recency of inbreeding.

RESULTS

MLPLs had significantly higher median number of ROH (188 interquartile range [IQR], 180.8-197.3), length of ROH (3.10, IQR 2.93-3.33), and total number of ROH (590.8, IQR 537.3-632.3), and FROH (0.26, IQR 0.24-0.28) when compared with the Controls and the EPLs (p < 0.05). There was no significant difference in any of the inbreeding indices between the EPLs and Controls. The MLPLs had a significantly higher proportion of long (>10 Mb) ROH (2.5%, IQR 1.6-3.6) than the Controls (1.7%, IQR 0.6-2.5), p = 0.001. No unique ROHs were found in the EPL or MLPL populations.

MAIN LIMITATIONS

SNP-array data does not allow analysis of every base in the sequence.

CONCLUSIONS

This first study of the effect of genomic inbreeding levels on pregnancy loss showed that inbreeding is a contributor to MLPL, but not EPL in the UK Thoroughbred population. Mating choices remain critical, because inbreeding may predispose to MLPL by increasing the risk of homozygosity for specific lethal allele(s).

Predicted genetic burden and frequency of phenotype-associated variants in the horse

Disease-causing variants have been identified for less than 20% of suspected equine genetic diseases. Whole genome sequencing (WGS) allows rapid identification of rare disease causal variants. However, interpreting the clinical variant consequence is confounded by the number of predicted deleterious variants that healthy individuals carry (predicted genetic burden). Estimation of the predicted genetic burden and baseline frequencies of known deleterious or phenotype associated variants within and across the major horse breeds have not been performed. We used WGS of 605 horses across 48 breeds to identify 32,818,945 variants, demonstrate a high predicted genetic burden (median 730 variants/horse, interquartile range: 613-829), show breed differences in predicted genetic burden across 12 target breeds, and estimate the high frequencies of some previously reported disease variants. This large-scale variant catalog for a major and highly athletic domestic animal species will enhance its ability to serve as a model for human phenotypes and improves our ability to discover the bases for important equine phenotypes.

Relationships between congenital peritoneopericardial diaphragmatic hernia or congenital central diaphragmatic hernia and ductal plate malformations in dogs and cats

OBJECTIVE

To characterize the association between peritoneopericardial diaphragmatic hernia (PPDH) or congenital central diaphragmatic hernia (CCDH) and ductal plate malformations (DPMs) in dogs and cats.

ANIMALS

18 dogs and 18 cats with PPDH or CCDH and 19 dogs and 18 cats without PPDH or CCDH.

PROCEDURES

Evaluation of clinical details verified PPDH or CCDH and survival times. Histologic features of nonherniated liver samples were used to categorize DPM. Immunohistochemical staining for cytokeratin-19 distinguished bile duct profiles per portal tract and for Ki-67-assessed cholangiocyte proliferation. Histologic features of herniated liver samples from PPDH or CCDH were compared with those of pathological controls (traumatic diaphragmatic hernia, n = 6; liver lobe torsion, 6; ischemic hepatopathy, 2).

RESULTS

DPM occurred in 13 of 18 dogs with the proliferative-like phenotype predominating and in 15 of 18 cats with evenly distributed proliferative-like and Caroli phenotypes. Congenital hepatic fibrosis DPM was noted in 3 dogs and 2 cats and renal DPM in 3 dogs and 3 cats. No signalment, clinical signs, or clinicopathologic features discriminated DPM. Kaplan Meier survival curves were similar in dogs and cats. Bile duct profiles per portal tract in dogs (median, 5.0; range, 1.4 to 100.8) and cats (6.6; 1.9 to 11.0) with congenital diaphragmatic hernias significantly exceeded those in healthy dogs (1.4; 1.2 to 1.6) and cats (2.3; 1.7 to 2.6). Animals with DPM lacked active cholangiocyte proliferation. Histologic features characterizing malformative bile duct profiles yet without biliary proliferation were preserved in herniated liver lobes in animals with DPM.

CONCLUSIONS AND CLINICAL RELEVANCE

DPM was strongly associated with PPDH and CCDH. Because DPM can impact health, awareness of its coexistence with PPDH or CCDH should prompt biopsy of nonherniated liver tissue during surgical correction of PPDH and CCDH.

Ductal Plate Malformations in Captive Snakes

Ductal plate malformations are abnormalities in the liver that arise from inappropriate or incomplete remodeling of the embryologic ductal plate. Various types of ductal plate malformations are reported in the human and veterinary literature, most commonly affecting domestic mammalian species but also fish. We investigated the occurrence and described the histopathologic features of ductal plate malformations in captive snakes. Malformations were identified in 18 snakes: 10 colubrids, 6 vipers, and 2 boids. There was no sex predilection, and the mean age was 17 years. The majority of lesions were incidental with most snakes having one or more comorbidities, most commonly neoplasia or systemic inflammation, that resulted in natural death or euthanasia. Ductal plate malformations in all livers were broadly characterized by a well-demarcated nodule of irregular bile ducts embedded within a varying amount of fibrous stroma. Malformations were further categorized based on the amount of fibrous stroma and dilation of the bile ducts as von Meyenburg complexes, cystic liver disease, and/or an intermediate hybrid subtype representative of cysts arising within von Meyenburg complexes. Histochemical and immunohistochemical staining, including Gomori's trichome and pan-cytokeratin, respectively, were applied on select cases to confirm histologic features. Malignant transformation was not identified within this population.

Quality of DNA extracted from formalin-fixed, paraffin-embedded canine tissues

Veterinary pathology tissue banks are valuable resources for genetic studies. However, limited data exist as to whether quality DNA can be extracted from these tissues for use in canine genotyping studies. We extracted DNA from 44 formalin-fixed, paraffin-embedded (FFPE) tissue blocks from dogs; 9 of these dogs had DNA available from whole blood samples that had been banked. We genotyped DNA from 30 of 44 tissue blocks and 9 whole blood samples on the Illumina CanineHD BeadChip; DNA quality was insufficient in 14 of 44 samples from tissue blocks. There was significant correlation between the 260/280 ratio and single-nucleotide variation (SNV) call rate (p = 0.0276; r² = 0.162); 23 of 30 samples from FFPE were genotyped with > 65% call rates. Median pairwise identical-by-state (IBS) analysis was 0.99 in 8 pairs of dogs with call rates > 65%. Neither age of tissue block nor specific tissue types were associated with significant differences in DNA concentration, 260/280 ratio, or SNV call rate. DNA extracted from tissue blocks can have variable quality, although comparable levels of homozygosity suggest that extracts from FFPE with call rates > 65% might provide similar results to samples from whole blood when analyzed on the Illumina CanineHD BeadChip.

Functionally Annotating Regulatory Elements in the Equine Genome Using Histone Mark ChIP-Seq

One of the primary aims of the Functional Annotation of ANimal Genomes (FAANG) initiative is to characterize tissue-specific regulation within animal genomes. To this end, we used chromatin immunoprecipitation followed by sequencing (ChIP-Seq) to map four histone modifications (H3K4me1, H3K4me3, H3K27ac, and H3K27me3) in eight prioritized tissues collected as part of the FAANG equine biobank from two thoroughbred mares. Data were generated according to optimized experimental parameters developed during quality control testing. To ensure that we obtained sufficient ChIP and successful peak-calling, data and peak-calls were assessed using six quality metrics, replicate comparisons, and site-specific evaluations. Tissue specificity was explored by identifying binding motifs within unique active regions, and motifs were further characterized by gene ontology (GO) and protein–protein interaction analyses. The histone marks identified in this study represent some of the first resources for tissue-specific regulation within the equine genome. As such, these publicly available annotation data can be used to advance equine studies investigating health, performance, reproduction, and other traits of economic interest in the horse.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Origin and Evolution of Deleterious Mutations in Horses

Domestication has changed the natural evolutionary trajectory of horses by favoring the reproduction of a limited number of animals showing traits of interest. Reduced breeding stocks hampered the elimination of deleterious variants by means of negative selection, ultimately inflating mutational loads. However, ancient genomics revealed that mutational loads remained steady during most of the domestication history until a sudden burst took place some 250 years ago. To identify the factors underlying this trajectory, we gather an extensive dataset consisting of 175 modern and 153 ancient genomes previously published, and carry out the most comprehensive characterization of deleterious mutations in horses. We confirm that deleterious variants segregated at low frequencies during the last 3500 years, and only spread and incremented their occurrence in the homozygous state during modern times, owing to inbreeding. This independently happened in multiple breeds, following both the development of closed studs and purebred lines, and the deprecation of horsepower in the 20th century, which brought many draft breeds close to extinction. Our work illustrates the paradoxical effect of some conservation and improvement programs, which reduced the overall genomic fitness and viability.

Caroli Syndrome in a 6-Year-Old Rottweiler Dog

A 6-year-old entire female rottweiler dog with a recent history of ascites and respiratory disease was submitted for necropsy examination. The dog had been diagnosed ultrasonographically with biliary cysts as a puppy. Grossly, the liver was smaller than expected with an irregular surface. Islands of hepatocytes were separated by bands of fibrosis and many bile ducts were markedly dilated. Histologically, extensive fibrosis extended beyond the limiting plate and into the surrounding hepatic parenchyma and was associated with abundant small bile ducts throughout. In conjunction with the detection of biliary cysts early in life, the gross and histological findings were consistent with a diagnosis of Caroli syndrome. In man, Caroli syndrome is frequently associated with renal and pancreatic cysts; a single renal cyst was identified in this case.

Congenital Hepatic Fibrosis in a Purebred Spanish Horse Foal: Pathology and Genetic Studies on PKHD1 Gene Mutations

A 1-month-old Purebred Spanish Horse (PSH) foal presented with progressive hepatic failure culminating in death. Hepatic lesions were consistent with congenital hepatic fibrosis (CHF). Genetic studies in the PKHD1 gene in the affected foal revealed that it was heterozygous for the 2 previously described single-nucleotide polymorphisms (SNPs) linked to CHF in Swiss Franches-Montagnes (SFM) horses. In addition, 2 novel mutations were detected, the foal being homozygous for one of them and heterozygous for the other. Genetic studies in a healthy PSH population ( n = 35) showed a 3-fold higher genotypic frequency for PKHD1 SNP g.49,630,834G>A and a 5-fold higher genotypic frequency for PKHD1 SNP g.49,597,760A>T compared with those reported for SFM horses. SNPs in the PKHD1 gene in CHF-affected SFM horses might not fully explain the CHF observed in the PSH. Other mutations in the PKHD1 gene could play a more important role in the PSH.

Ductal Plate Malformation in the Liver of Boxer Dogs: Clinical and Histological Features

Ductal plate malformations (DPMs) represent developmental biliary disorders with a wide phenotypic spectrum. This study characterizes DPM in 30 Boxer dogs. Median age was 1.5 (range, 0.3-10.0) years, with 12 dogs <1 year. Clinical features included increased serum levels of liver enzymes (28), gastrointestinal signs (16), poor body condition (14), abdominal effusion (9), and hepatic encephalopathy (2). Additional malformations included gallbladder atresia (8), atrophied left liver (2), absent quadrate lobe with left-displaced gallbladder (1), portal vasculature atresia (left liver, 1), intrahepatic portosystemic shunt (1), and complex intrahepatic arteriovenous malformation (1). All dogs had portal tracts dimensionally expanded by a moderate-to-severe multiple small bile duct phenotype embedded in abundant extracellular matrix; 80% displayed variable portal-to-portal bridging. Quantitative analysis confirmed significantly increased fibrillar collagen and a 3-fold increased portal tract area relative to 6 Boxer and 10 non-Boxer controls. Biliary phenotype was dominated by tightly formed CK19-positive ductules, typically 10 to 15 μm in diameter, with 3 to >30 profiles per portal tract, reduced luminal apertures, and negative Ki-67 immunoreactivity. CK19-positive biliary epithelium intersected directly with zone 1 hepatocytes as a signature feature when considered with other DPM characteristics. Phenotypic variation included a multiple small bile duct phenotype (all dogs), predominantly thin-walled sacculated ducts (4), well-formed saccular ducts (4), and sacculated segmental, interlobular, and intralobular ducts (Caroli malformation, 2 dogs, one with bridging portal fibrosis). Histologic evidence of portal venous hypoperfusion accompanied increased biliary profiles in every case. We propose that this spectrum of disorders be referred to as DPM with appropriate modifiers to characterize the unique phenotypes.