MFSD2A frameshift variant in Kerry Hill sheep with microcephaly.

Microcephaly is a rare neurodevelopmental disorder characterized by reduced skull circumference and brain volume that occurs sporadically in farm animals. We investigated an early-onset neurodegenerative disorder observed in seven lambs of purebred Kerry Hill sheep. Clinical signs included inability to stand or severe ataxia, convulsions, and early death. Diagnostic imaging and brain necropsy confirmed microcephaly. The pedigree of the lambs suggested monogenic autosomal recessive inheritance. We sequenced the genome of one affected lamb, and comparison with 115 control genomes revealed a single private protein-changing variant. This frameshift variant, MFSD2A: c.285dupA, p.(Asp96fs*9), represents a 1-bp duplication predicted to truncate 80% of the open reading frame. MFSD2A is a transmembrane protein that is essential for maintaining blood-brain barrier homeostasis and plays a key role in regulating brain lipogenesis. Human MFSD2A pathogenic variants are associated with a neurodevelopmental disorder with progressive microcephaly, spasticity, and brain imaging abnormalities (NEDMISBA, OMIM 616486). Here we present evidence for the occurrence of a recessively inherited form of microcephaly in sheep due to a loss-of-function variant in MFSD2A (OMIA 002371-9940). To the best of our knowledge, this is the first report of a spontaneous MFSD2A variant in domestic animals.

Loss-of-function variant in the ovine TMCO6 gene in North Country Cheviot sheep with motor neuron disease.

In North Country Cheviot lambs with early-onset progressive ataxia and motor neuron degeneration, whole-genome sequencing identified a homozygous loss-of-function variant in the ovine transmembrane and coiled-coil domains (TMCO6) gene. The familial recessive form of motor neuron disease in sheep is due to a pathogenic 4 bp deletion leading to a 50% protein truncation that is assumed to result in the absence of a functional TMCO6. This uncharacterised protein is proposed to interact with ubiquilin 1 which is associated with Alzheimer's disease, whereas sporadic forms of amyotrophic lateral sclerosis are caused by variants in UBQLN2. Our findings provide a first spontaneous animal model for TMCO6, which could have implications in the studies of other comparative neurodegenerative diseases. In addition, these results will allow the design of a genetic test to prevent the occurrence of this fatal disease in the affected sheep population.

Compound heterozygous PLA2G6 loss-of-function variants in Swaledale sheep with neuroaxonal dystrophy.

Sporadic occurrences of neurodegenerative disorders including neuroaxonal dystrophy (NAD) have been previously reported in sheep. However, so far no causative genetic variant has been found for ovine NAD. The aim of this study was to characterize the phenotype and the genetic aetiology of an early-onset neurodegenerative disorder observed in several lambs of purebred Swaledale sheep, a native English breed. Affected lambs showed progressive ataxia and stiff gait and subsequent histopathological analysis revealed the widespread presence of axonal spheroid indicating neuronal degeneration. Thus, the observed clinical phenotype could be explained by a novel form of NAD. After SNP genotyping and subsequent linkage mapping within a paternal half-sib pedigree with a total of five NAD-affected lambs, we identified two loss-of-function variants by whole-genome sequencing in the ovine PLA2G6 gene situated in a NAD-linked genome region on chromosome 3. All cases were carriers of a compound heterozygous splice site variant in intron 2 and a nonsense variant in exon 8. Herein we present evidence for the occurrence of a familial novel form of recessively inherited NAD in sheep due to allelic heterogeneity at PLA2G6. This study reports two pathogenic variants in PLA2G6 causing a novel form of NAD in Swaledale sheep which enables selection against this fatal disorder.

A frameshift insertion in FA2H causes a recessively inherited form of ichthyosis congenita in Chianina cattle.

The aim of this study was to characterize the phenotype and to identify the genetic etiology of a syndromic form of ichthyosis congenita (IC) observed in Italian Chianina cattle and to estimate the prevalence of the deleterious allele in the population. Sporadic occurrence of different forms of ichthyosis including IC have been previously reported in cattle. However, so far, no causative genetic variant has been found for bovine IC. Nine affected cattle presenting congenital xerosis, hyperkeratosis and scaling of the skin as well as urolithiasis and cystitis associated with retarded growth were examined. Skin histopathology revealed a severe, diffuse orthokeratotic hyperkeratosis with mild to moderate epidermal hyperplasia. The pedigree records indicated a monogenic recessive trait. Homozygosity mapping and whole-genome sequencing allowed the identification of a homozygous frameshift 1 bp insertion in the FA2H gene (c.9dupC; p.Ala4ArgfsTer142) located in a 1.92 Mb shared identical-by-descent region on chromosome 18 present in all cases, while the parents were heterozygous as expected for obligate carriers. These findings enable the selection against this sub-lethal allele showing an estimated frequency of ~ 7.5% in Chianina top sires. A sporadic incidence of mild clinical signs in the skin of heterozygous carriers was observed. So far, pathogenic variants affecting the encoded fatty acid 2-hydroxylase catalyzing the synthesis of 2-hydroxysphingolipids have been associated with myelin disorders. In conclusion, this study represents the first report of an FA2H-related autosomal recessive inherited skin disorder in a mammalian species and adds FA2H to the list of candidate genes for ichthyosis in humans and animals. Furthermore, this study provides a DNA-based diagnostic test that enables selection against the identified pathogenic variant in the Chianina cattle population. However, functional studies are needed to better understand the expression of FA2H in IC-affected Chianina cattle.

Genomic diversity and population structure of the Leonberger dog breed.

BACKGROUND: Leonberger is a giant dog breed formed in the 1850s in Germany. Its post-World War II popularity has resulted in a current global population of ~ 30,000 dogs. The breed has predispositions to neurodegenerative disorders and cancer, which is likely due in large part to limited genetic diversity. However, to date there is no scientific literature on the overall demography and genomic architecture of this breed. RESULTS: We assessed extensive pedigree records, SNP array genotype data, and whole-genome sequences (WGS) on 142,072, 1203 and 39 Leonberger dogs, respectively. Pedigree analyses identified 22 founder animals and revealed an apparent popular sire effect. The average pedigree-based inbreeding coefficient of 0.29 and average kinship of 0.31 show a dramatic loss of genetic diversity. The observed average life span decreased over time from 9.4 years in 1989 to 7.7 years in 2004. A global health survey confirmed a high prevalence of cancer and neurological disorders. Analysis of SNP-based runs of homozygosity (ROH) identified 125,653 ROH with an average length of 5.88 Mb, and confirmed an average inbreeding coefficient of 0.28. Genome-wide filtering of the WGS data revealed 28 non-protein-changing variants that were present in all Leonberger individuals and a list of 22 potentially pathogenic variants for neurological disorders of which 50% occurred only in Leonbergers and 50% occurred rarely in other breeds. Furthermore, one of the two mtDNA haplogroups detected was present in one dog only. CONCLUSIONS: The increasing size of the Leonberger population has been accompanied by a considerable loss of genetic diversity after the bottleneck that occurred in the 1940s due to the intensive use of popular sires resulting in high levels of inbreeding. This might explain the high prevalence of certain disorders; however, genomic data provide no evidence for fixed coding variants that explain these predispositions. The list of candidate causative variants for polyneuropathy needs to be further evaluated. Preserving the current genetic diversity is possible by increasing the number of individuals for breeding while restricting the number of litters per sire/dam. In addition, outcrossing would help optimize long-term genetic diversity and contribute to the sustainability and health of the population.

Correction to: Genomic diversity and population structure of the Leonberger dog breed.

An amendment to this paper has been published and can be accessed via the original article.

Chromosomal imbalance in pigs showing a syndromic form of cleft palate.

BACKGROUND: Palatoschisis or cleft palate is a known anomaly in pigs resulting in their death. However, little is known about its aetiology. A detailed description of the phenotype was derived from necropsy and by computed tomography revealing that all 20 cases also exhibited hypodontia and renal cysts. Furthermore, a genetic origin was assumed due to dominant inheritance as all 20 recorded cases were confirmed offspring of a single boar. RESULTS: Single nucleotide variant (SNV) genotyping data were used to map the defect in the porcine genome and led to the detection of a chromosomal imbalance in the affected offspring. Whole genome sequencing of an affected piglet and a normal full sib was used to identify a chromosomal translocation and to fine map the breakpoints in the genome. Finally, we proved that the boar, which sired the malformed piglets, carried a balanced translocation. The detected translocation of Mb-sized segments of chromosome 8 and 14 had not been previously observed during karyotyping. All affected offspring were shown to be carriers of a partial trisomy of chromosome 14 including the FGFR2 gene, which is associated with various dominant inherited craniofacial dysostosis syndromes in man, and partial monosomy of chromosome 8 containing MSX1 known to be associated with tooth agenesis and orofacial clefts in other species. CONCLUSIONS: This study illustrates the usefulness of recently established genomic resources in pigs. In this study, the application of genome-wide genotyping and sequencing methods allowed the identification of the responsible boar and the genetic cause of the observed defect. By implementing systematic surveillance, it is possible to identify genetic defects at an early stage and avoid further distribution of congenital disorders.

A GJA9 frameshift variant is associated with polyneuropathy in Leonberger dogs.

BACKGROUND: Many inherited polyneuropathies (PN) observed in dogs have clinical similarities to the genetically heterogeneous group of Charcot-Marie-Tooth (CMT) peripheral neuropathies in humans. The canine disorders collectively show a variable expression of progressive clinical signs and ages of onset, and different breed prevalences. Previously in the Leonberger breed, a variant highly associated with a juvenile-onset PN was identified in the canine orthologue of a CMT-associated gene. As this deletion in ARHGEF10 (termed LPN1) does not explain all cases, PN in this breed may encompass variants in several genes with similar clinical and histopathological features. RESULTS: A genome-wide comparison of 173 k SNP genotypes of 176 cases, excluding dogs homozygous for the ARHGEF10 variant, and 138 controls, was carried out to detect further PN-associated variants. A single suggestive significant association signal on CFA15 was found. The genome of a PN-affected Leonberger homozygous for the associated haplotype was sequenced and variants in the 7.7 Mb sized critical interval were identified. These variants were filtered against a database of variants observed in 202 genomes of various dog breeds and 3 wolves, and 6 private variants in protein-coding genes, all in complete linkage disequilibrium, plus 92 non-coding variants were revealed. Five of the coding variants were predicted to have low or moderate effect on the encoded protein, whereas a 2 bp deletion in GJA9 results in a frameshift of high impact. GJA9 encodes connexin 59, a connexin gap junction family protein, and belongs to a group of CMT-associated genes that have emerged as important components of peripheral myelinated nerve fibers. The association between the GJA9 variant and PN was confirmed in an independent cohort of 296 cases and 312 controls. Population studies showed a dominant mode of inheritance, an average age of onset of approximately 6 years, and incomplete penetrance. CONCLUSIONS: This GJA9 variant represents a highly probable candidate variant for another form of PN in Leonberger dogs, which we have designated LPN2, and a new candidate gene for CMT disease. To date, approximately every third PN-diagnosed Leonberger dog can be explained by the ARHGEF10 or GJA9 variants, and we assume that additional genetic heterogeneity in this condition exists in the breed.

Congenital syndromic Chiari-like malformation (CSCM) in Holstein cattle: towards unravelling of possible genetic causes.

BACKGROUND: Chiari malformation type II (CMII) was originally reported in humans as a rare disorder characterized by the downward herniation of the hindbrain and towering cerebellum. The congenital brain malformation is usually accompanied by spina bifida, a congenital spinal anomaly resulting from incomplete closure of the dorsal aspect of the spinal neural tube, and occasionally by other lesions. A similar disorder has been reported in several animal species, including cattle, particularly as a congenital syndrome. A cause of congenital syndromic Chiari-like malformation (CSCM) in cattle has not been reported to date. We collected a series of 14 CSCM-affected Holstein calves (13 purebred, one Red Danish Dairy F1 cross) and performed whole-genome sequencing (WGS). WGS was performed on 33 cattle, including eight cases with parents (trio-based; group 1), three cases with one parent (group 2), and three single cases (solo-based; group 3). RESULTS: Sequencing-based genome-wide association study of the 13 Holstein calves with CSCM and 166 controls revealed no significantly associated genome region. Assuming a single Holstein breed-specific recessive allele, no region of shared homozygosity was detected suggesting heterogeneity. Subsequent filtering for protein-changing variants that were only homozygous in the genomes of the individual cases allowed the identification of two missense variants affecting different genes, SHC4 in case 4 in group 1 and WDR45B in case 13 in group 3. Furthermore, these two variants were only observed in Holstein cattle when querying WGS data of > 5,100 animals. Alternatively, potential de novo mutational events were assessed in each case. Filtering for heterozygous private protein-changing variants identified one DYNC1H1 frameshift variant as a candidate causal dominant acting allele in case 12 in group 3. Finally, the presence of larger structural DNA variants and chromosomal abnormalities was investigated in all cases. Depth of coverage analysis revealed two different partial monosomies of chromosome 2 segments in cases 1 and 7 in group 1 and a trisomy of chromosome 12 in the WDR45B homozygous case 13 in group 3. CONCLUSIONS: This study presents for the first time a detailed genomic evaluation of CSCM in Holstein cattle and suggests an unexpected genetic and allelic heterogeneity considering the mode of inheritance, as well as the type of variant. For the first time, we propose candidate causal variants that may explain bovine CSCM in a certain proportion of affected calves. We present cattle as a large animal model for human CMII and propose new genes and genomic variants as possible causes for related diseases in both animals and humans.

PAX3 haploinsufficiency in Maine Coon cats with dominant blue eyes and hearing loss resembling the human Waardenburg syndrome.

This study investigated the dominant blue eyes (DBE) trait linked to hearing impairment and variable white spotting in Maine Coon cats. Fifty-eight animals descending from two different DBE lineages, the Dutch and the Topaz lines, were sampled. They comprised 48 cats from the Dutch bloodline, including 9 green-eyed and 31 blue-eyed cats, with some individuals exhibiting signs of deafness, and 8 stillborn kittens. Samples from the Topaz lineage included ten blue-eyed animals. A brainstem auditory evoked potential test (BAER) revealed a reduced to absent response to auditory stimuli and absent physiological waveforms in all of the eight examined DBE animals. We sequenced the genome of two affected cats from the Dutch line and searched for variants in 19 candidate genes for the human Waardenburg syndrome and pigmentary disorders. This search yielded nine private protein-changing candidate variants in the genes PAX3, EDN3, KIT, OCA2, SLC24A5, HERC2 and TYRP1. The genotype-phenotype co-segregation was observed for the PAX3 variant within all animals from the Dutch lineage. The mutant allele was absent from 461 control genomes and 241 additionally genotyped green-eyed Maine Coons. We considered the PAX3 variant as the most plausible candidate -a heterozygous nonsense single basepair substitution in exon 6 of PAX3 (NC_051841.1: g.205,787,310G>A, XM_019838731.3:c.937C>T, XP_019694290.1:p.Gln313*), predicted to result in a premature stop codon. PAX3 variants cause auditory-pigmentary syndrome in humans, horses, and mice. Together with the comparative data from other species, our findings strongly suggest PAX3:c.937C>T (OMIA:001688-9685) as the most likely candidate variant for the DBE, deafness and minimal white spotting in the Maine Coon Dutch line. Finally, we propose the designation of DBERE (Rociri Elvis Dominant Blue Eyes) allele in the domestic cat.

Organoids of patient-derived medullary thyroid carcinoma: The first milestone towards a new in vitro model in dogs.

Organoid cultures could constitute a valuable in vitro model to explore new treatments for canine (c) medullary thyroid carcinoma (MTC). The study's objectives were to establish and characterize 3D organoid cultures of cMTC using histology and immunohistochemistry (IHC) and to evaluate the effect of antitumor drugs on organoids' viability. Five cMTC tissue samples were used to develop organoid cultures of which one organoid line, named cMTC N°2, could be passaged for an extended period. This cMTC N°2 organoid line was further compared to the primary tumour regarding morphology and IHC expression of thyroid transcription factor-1 (TTF-1), thyroglobulin, calcitonin, synaptophysin, vimentin, Ki-67, cyclooxygenase-2 (COX-2), P-glycoprotein and vascular endothelial growth factor (VEGF). Quality control of the cMTC N°2 organoid line was achieved by a single nucleotide polymorphism (SNP) array of the organoids, primary tumour and healthy blood cells of the same dog. The effect of carboplatin, meloxicam and toceranib phosphate (TOC) on cMTC N°2 organoids' viability was evaluated. The cMTC N°2 organoid line was cultured for 94 days and showed similar histological features with the primary tumour. Immunolabelling for TTF-1, thyroglobulin, calcitonin and VEGF was similar between the primary tumour and cMTC N°2 organoids. Compared to the primary tumour, organoids showed higher immunolabelling for vimentin and Ki-67, and lower immunolabelling for synaptophysin, COX-2 and P-glycoprotein. The SNP genotype was similar for each chromosome between healthy blood cells, primary tumour and cMTC N°2 organoids. Carboplatin, meloxicam and TOC had no effect on cMTC N°2 organoid cell viability within achievable in vivo concentration range. In conclusion, the cMTC N°2 organoid line is a promising first milestone towards an established in vitro organoid model to explore pathophysiology and new treatment modalities in cMTC.

Independent COL17A1 Variants in Cats with Junctional Epidermolysis Bullosa.

Epidermolysis bullosa (EB), characterized by defective adhesion of the epidermis to the dermis, is a heterogeneous disease with many subtypes in human patients and domestic animals. We investigated two unrelated cats with recurring erosions and ulcers on ear pinnae, oral mucosa, and paw pads that were suggestive of EB. Histopathology confirmed the diagnosis of EB in both cats. Case 1 was severe and had to be euthanized at 5 months of age. Case 2 had a milder course and was alive at 11 years of age at the time of writing. Whole genome sequencing of both affected cats revealed independent homozygous variants in COL17A1 encoding the collagen type XVII alpha 1 chain. Loss of function variants in COL17A1 lead to junctional epidermolysis bullosa (JEB) in human patients. The identified splice site variant in case 1, c.3019+1del, was predicted to lead to a complete deficiency in collagen type XVII. Case 2 had a splice region variant, c.769+5G>A. Assessment of the functional impact of this variant on the transcript level demonstrated partial aberrant splicing with residual expression of wildtype transcript. Thus, the molecular analyses provided a plausible explanation of the difference in clinical severity between the two cases and allowed the refinement of the diagnosis in the affected cats to JEB. This study highlights the complexity of EB in animals and contributes to a better understanding of the genotype-phenotype correlation in COL17A1-related JEB.

Genomic Diversity and Runs of Homozygosity in Bernese Mountain Dogs.

Bernese mountain dogs are a large dog breed formed in the early 1900s in Switzerland. While originally farm dogs that were used for pulling carts, guarding, and driving cattle, today they are considered multi-purpose companion and family dogs. The breed is predisposed to several complex diseases, such as histiocytic sarcoma, degenerative myelopathy, or hip dysplasia. Using whole-genome sequencing (WGS) data, we assessed the genomic architecture of 33 unrelated dogs from four countries: France, Sweden, Switzerland, and the United States. Analysis of runs of homozygosity (ROH) identified 12,643 ROH with an average length of 2.29 Mb and an average inbreeding coefficient of 0.395. Multidimensional scaling analysis of the genetic relatedness revealed limited clustering of European versus USA dogs, suggesting exchanges of breeding stock between continents. Furthermore, only two mtDNA haplotypes were detected in the 33 studied dogs, both of which are widespread throughout multiple dog breeds. WGS-based ROH analyses revealed several fixed or nearly fixed regions harboring discreet morphological trait-associated as well as disease-associated genetic variants. Several genes involved in the regulation of immune cells were found in the ROH shared by all dogs, which is notable in the context of the breed's strong predisposition to hematopoietic cancers. High levels of inbreeding and relatedness, strongly exaggerated in the last 30 years, have likely led to the high prevalence of specific genetic disorders in this breed.

SGCD Missense Variant in a Lagotto Romagnolo Dog with Autosomal Recessively Inherited Limb-Girdle Muscular Dystrophy.

An 8-month-old female Lagotto Romagnolo dog was presented for a 1-month history of an initial severe reluctance to move, rapidly progressing to a marked stiff gait and progressive muscular weakness and evolving to tetraparesis, which persuaded the owner to request euthanasia. A primary muscle pathology was supported by necropsy and histopathological findings. Macroscopically, the muscles were moderately atrophic, except for the diaphragm and the neck muscles, which were markedly thickened. Histologically, all the skeletal muscles examined showed atrophy, hypertrophy, necrosis with calcification of the fibers, and mild fibrosis and inflammation. On immunohistochemistry, all three dystrophin domains and sarcoglycan proteins were absent. On Western blot analysis, no band was present for delta sarcoglycan. We sequenced the genome of the affected dog and compared the data to more than 900 control genomes of different dog breeds. Genetic analysis revealed a homozygous private protein-changing variant in the SGCD gene encoding delta- sarcoglycan in the affected dog. The variant was predicted to induce a SGCD:p.(Leu242Pro) change in the protein. In silico tools predicted the change to be deleterious. Other 770 Lagotto Romagnolo dogs were genotyped for the variant and all found to be homozygous wild type. Based on current knowledge of gene function in other mammalian species, including humans, hamsters, and dogs, we propose the SGCD missense variant as the causative variant of the observed form of muscular dystrophy in the index case. The absence of the variant allele in the Lagotto Romagnolo breeding population indicates a rare allele that has appeared recently.

A living biobank of canine mammary tumor organoids as a comparative model for human breast cancer.

Mammary tumors in dogs hold great potential as naturally occurring breast cancer models in translational oncology, as they share the same environmental risk factors, key histological features, hormone receptor expression patterns, prognostic factors, and genetic characteristics as their human counterparts. We aimed to develop in vitro tools that allow functional analysis of canine mammary tumors (CMT), as we have a poor understanding of the underlying biology that drives the growth of these heterogeneous tumors. We established the long-term culture of 24 organoid lines from 16 dogs, including organoids derived from normal mammary epithelium or benign lesions. CMT organoids recapitulated key morphological and immunohistological features of the primary tissue from which they were derived, including hormone receptor status. Furthermore, genetic characteristics (driver gene mutations, DNA copy number variations, and single-nucleotide variants) were conserved within tumor-organoid pairs. We show how CMT organoids are a suitable model for in vitro drug assays and can be used to investigate whether specific mutations predict therapy outcomes. Specifically, certain CMT subtypes, such as PIK3CA mutated, estrogen receptor-positive simple carcinomas, can be valuable in setting up a preclinical model highly relevant to human breast cancer research. In addition, we could genetically modify the CMT organoids and use them to perform pooled CRISPR/Cas9 screening, where library representation was accurately maintained. In summary, we present a robust 3D in vitro preclinical model that can be used in translational research, where organoids from normal, benign as well as malignant mammary tissues can be propagated from the same animal to study tumorigenesis.

Phenotypic and Genomic Analysis of Cystic Hygroma in Pigs.

Cystic hygroma is a malformation of the lymphatic and vascular system and is recognized as a benign congenital tumor that affects humans and animals in the perinatal period. This congeni-tal disorder is rarely described in animals, and until today, cystic hygroma in pigs has not been described in the literature. In a purebred Piètrain litter with twelve live-born piglets, cystic hy-groma was noticed on the rump of two male pigs within the first week of life. In addition, a third case of a crossbred weaner (Large White × Landrace) was detected during a herd examina-tion. To rule out common differential diagnoses, e.g., abscess or hematoma, further clinical and pathological investigations were conducted. During clinical examination, a painless and soft mass, which was compressible, was detected on the rump of all affected animals. The ultra-sonographic examination revealed a fluid-filled and cavernous subcutaneous structure. In addi-tion, a puncture of the cyst was conducted, revealing a serosanguinous fluid with negative bacte-riological culture. In all cases, a necropsy was performed, showing that the animals had fluid-filled cysts lined by well-differentiated lymphatic endothelium. Based on the clinicopathological examination, cystic hygroma was diagnosed. Furthermore, SNP array genotyping and whole-genome sequencing was performed and provided no evidence for a chromosomal disorder. In the Piètrain family, several genome regions were homozygous in both affected piglets. None-theless, a dominant acting de novo germline variant could not be ruled out, and therefore differ-ent filtering strategies were used to find pathogenic variants. The herein presented lists of pri-vate variants after filtering against hundreds of control genomes provide no plausible candidate and no shared variants among the two sequenced cases. Therefore, further studies are needed to evaluate possible genetic etiology. In general, systematic surveillance is needed to identify ge-netic defects as early as possible and to avoid the occurrence of losses in the pig population.