Canine models of human amelogenesis imperfecta: identification of novel recessive ENAM and ACP4 variants

Molecular Evidence for Relaxed Selection on the Enamel Genes of Toothed Whales (Odontoceti) with Degenerative Enamel Phenotypes

Different species of toothed whales (Odontoceti) exhibit a variety of tooth forms and enamel types. Some odontocetes have highly prismatic enamel with Hunter-Schreger bands, whereas enamel is vestigial or entirely lacking in other species. Different tooth forms and enamel types are associated with alternate feeding strategies that range from biting and grasping prey with teeth in most oceanic and river dolphins to the suction feeding of softer prey items without the use of teeth in many beaked whales. At the molecular level, previous studies have documented inactivating mutations in the enamel-specific genes of some odontocete species that lack complex enamel. At a broader scale, however, it is unclear whether enamel complexity across the full diversity of extant Odontoceti correlates with the relative strength of purifying selection on enamel-specific genes. Here, we employ sequence alignments for seven enamel-specific genes (ACP4, AMBN, AMELX, AMTN, ENAM, KLK4, MMP20) in 62 odontocete species that are representative of all extant families. The sequences for 33 odontocete species were obtained from databases, and sequences for the remaining 29 species were newly generated for this study. We screened these alignments for inactivating mutations (e.g., frameshift indels) and provide a comprehensive catalog of these mutations in species with one or more inactivated enamel genes. Inactivating mutations are rare in Delphinidae (oceanic dolphins) and Platanistidae/Inioidea (river dolphins) that have higher enamel complexity scores. By contrast, mutations are much more numerous in clades such as Monodontidae (narwhal, beluga), Ziphiidae (beaked whales), Physeteroidea (sperm whales), and Phocoenidae (porpoises) that are characterized by simpler enamel or even enamelless teeth. Further, several higher-level taxa (e.g., Hyperoodon, Kogiidae, Monodontidae) possess shared inactivating mutations in one or more enamel genes, which suggests loss of function of these genes in the common ancestor of each clade. We also performed selection (dN/dS) analyses on a concatenation of these genes and used linear regression and Spearman's rank-order correlation to test for correlations between enamel complexity and two different measures of selection intensity (# of inactivating mutations per million years, dN/dS values). Selection analyses revealed that relaxed purifying selection is especially prominent in physeteroids, monodontids, and phocoenids. Linear regressions and correlation analyses revealed a strong negative correlation between selective pressure (dN/dS values) and enamel complexity. Stronger purifying selection (low dN/dS) is found on branches with more complex enamel and weaker purifying selection (higher dN/dS) occurs on branches with less complex enamel or enamelless teeth. As odontocetes diversified into a variety of feeding modes, in particular, the suction capture of prey, a reduced reliance on the dentition for prey capture resulted in the relaxed selection of genes that are critical to enamel development.

Enamel Structure Defects in Kdf1 Missense Mutation Knock-in Mice

The Keratinocyte differentiation factor 1 (KDF1) is reported to take part in tooth formation in humans, but the dental phenotype of Kdf1 mutant mice has not been understood. Additionally, the role of the KDF1 gene in dental hard tissue development is rarely known. In this study, we constructed a Kdf1 missense mutation knock-in mouse model through CRISPR/Cas9 gene-editing technology. Enamel samples from wildtypes (WT) and Kdf1 homozygous mutants (HO) were examined using micro-computed tomography (micro-CT), scanning electron microscopy (SEM), an atomic force microscope (AFM) and Raman microspectroscopy. The results showed that a novel Kdf1 missense mutation (c. 908G>C, p.R303P) knock-in mice model was constructed successfully. The enamel of HO mice incisors appeared chalky and defective, exposing the rough interior of the inner enamel and dentin. Micro-CT showed that HO mice had lower volume and mineral density in their tooth enamel. In addition, declined thickness was found in the unerupted enamel layer of incisors in the HO mice. Using SEM and AFM, it was found that enamel prisms in HO mice enamel were abnormally and variously shaped with loose decussating crystal arrangement, meanwhile the enamel rods were partially fused and collapsed, accompanied by large gaps. Furthermore, misshapen nanofibrous apatites were disorderly combined with each other. Raman microspectroscopy revealed a compromised degree of order within the crystals in the enamel after the Kdf1 mutation. To conclude, we identified enamel structure defects in the Kdf1 missense mutation knock-in mice, which displayed fragmentary appearance, abnormally shaped prism structure, decreased mineral density, altered crystal ordering degree and chemical composition of the enamel layer. This may support the potential role of the KDF1 gene in the natural development of enamel.

Enamel defects in Acp4R110C/R110C mice and human ACP4 mutations

Human ACP4 (OMIM*606362) encodes a transmembrane protein that belongs to histidine acid phosphatase (ACP) family. Recessive mutations in ACP4 cause non-syndromic hypoplastic amelogenesis imperfecta (AI1J, OMIM#617297). While ACP activity has long been detected in developing teeth, its functions during tooth development and the pathogenesis of ACP4-associated AI remain largely unknown. Here, we characterized 2 AI1J families and identified a novel ACP4 disease-causing mutation: c.774_775del, p.Gly260Aspfs*29. To investigate the role of ACP4 during amelogenesis, we generated and characterized Acp4R110C mice that carry the p.(Arg110Cys) loss-of-function mutation. Mouse Acp4 expression was the strongest at secretory stage ameloblasts, and the protein localized primarily at Tomes' processes. While Acp4 heterozygous (Acp4+/R110C) mice showed no phenotypes, incisors and molars of homozygous (Acp4R110C/R110C) mice exhibited a thin layer of aplastic enamel with numerous ectopic mineralized nodules. Acp4R110C/R110C ameloblasts appeared normal initially but underwent pathology at mid-way of secretory stage. Ultrastructurally, sporadic enamel ribbons grew on mineralized dentin but failed to elongate, and aberrant needle-like crystals formed instead. Globs of organic matrix accumulated by the distal membranes of defective Tomes' processes. These results demonstrated a critical role for ACP4 in appositional growth of dental enamel probably by processing and regulating enamel matrix proteins around mineralization front apparatus.

Missense variant in LOXHD1 is associated with canine nonsyndromic hearing loss

Hearing loss is a common sensory deficit in both humans and dogs. In canines, the genetic basis is largely unknown, as genetic variants have only been identified for a syndromic form of hearing impairment. We observed a congenital or early-onset sensorineural hearing loss in a Rottweiler litter. Assuming an autosomal recessive inheritance, we used a combined approach of homozygosity mapping and genome sequencing to dissect the genetic background of the disorder. We identified a fully segregating missense variant in LOXHD1, a gene that is known to be essential for cochlear hair cell function and associated with nonsyndromic hearing loss in humans and mice. The canine LOXHD1 variant was specific to the Rottweiler breed in our study cohorts of pure-bred dogs. However, it also was present in some mixed-breed dogs, of which the majority showed Rottweiler ancestry. Low allele frequencies in these populations, 2.6% and 0.04%, indicate a rare variant. To summarize, our study describes the first genetic variant for canine nonsyndromic hearing loss, which is clinically and genetically similar to human LOXHD1-related hearing disorder, and therefore, provides a new large animal model for hearing loss. Equally important, the affected breed will benefit from a genetic test to eradicate this LOXHD1-related hearing disorder from the population.

In-frame deletion in canine PITRM1 is associated with a severe early-onset epilepsy, mitochondrial dysfunction and neurodegeneration

We investigated the clinical, genetic, and pathological characteristics of a previously unknown severe juvenile brain disorder in several litters of Parson Russel Terriers. The disease started with epileptic seizures at 6–12 weeks of age and progressed rapidly to status epilepticus and death or euthanasia. Histopathological changes at autopsy were restricted to the brain. There was severe acute neuronal degeneration and necrosis diffusely affecting the grey matter throughout the brain with extensive intraneuronal mitochondrial crowding and accumulation of amyloid-β (Aβ). Combined homozygosity mapping and genome sequencing revealed an in-frame 6-bp deletion in the nuclear-encoded pitrilysin metallopeptidase 1 (PITRM1) encoding for a mitochondrial protease involved in mitochondrial targeting sequence processing and degradation. The 6-bp deletion results in the loss of two amino acid residues in the N-terminal part of PITRM1, potentially affecting protein folding and function. Assessment of the mitochondrial function in the affected brain tissue showed a significant deficiency in respiratory chain function. The functional consequences of the mutation were modeled in yeast and showed impaired growth in permissive conditions and an impaired respiration capacity. Loss-of-function variants in human PITRM1 result in a childhood-onset progressive amyloidotic neurological syndrome characterized by spinocerebellar ataxia with behavioral, psychiatric and cognitive abnormalities. Homozygous Pitrm1-knockout mice are embryonic lethal, while heterozygotes show a progressive, neurodegenerative phenotype characterized by impairment in motor coordination and Aβ deposits. Our study describes a novel early-onset PITRM1-related neurodegenerative canine brain disorder with mitochondrial dysfunction, Aβ accumulation, and lethal epilepsy. The findings highlight the essential role of PITRM1 in neuronal survival and strengthen the connection between mitochondrial dysfunction and neurodegeneration.

A novel genomic region on chromosome 11 associated with fearfulness in dogs

The complex phenotypic and genetic nature of anxieties hampers progress in unravelling their molecular etiologies. Dogs present extensive natural variation in fear and anxiety behaviour and could advance the understanding of the molecular background of behaviour due to their unique breeding history and genetic architecture. As dogs live as part of human families under constant care and monitoring, information from their behaviour and experiences are easily available. Here we have studied the genetic background of fearfulness in the Great Dane breed. Dogs were scored and categorised into cases and controls based on the results of the validated owner-completed behavioural survey. A genome-wide association study in a cohort of 124 dogs with and without socialisation as a covariate revealed a genome-wide significant locus on chromosome 11. Whole exome sequencing and whole genome sequencing revealed extensive regions of opposite homozygosity in the same locus on chromosome 11 between the cases and controls with interesting neuronal candidate genes such as MAPK9/JNK2, a known hippocampal regulator of anxiety. Further characterisation of the identified locus will pave the way for molecular understanding of fear in dogs and may provide a natural animal model for human anxieties.

A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves

The domestic dog serves as an excellent model to investigate the genetic basis of disease. More than 400 heritable traits analogous to human diseases have been described in dogs. To further canine medical genetics research, we established the Dog Biomedical Variant Database Consortium (DBVDC) and present a comprehensive list of functionally annotated genome variants that were identified with whole genome sequencing of 582 dogs from 126 breeds and eight wolves. The genomes used in the study have a minimum coverage of 10× and an average coverage of ~24×. In total, we identified 23 133 692 single-nucleotide variants (SNVs) and 10 048 038 short indels, including 93% undescribed variants. On average, each individual dog genome carried ∼4.1 million single-nucleotide and ~1.4 million short-indel variants with respect to the reference genome assembly. About 2% of the variants were located in coding regions of annotated genes and loci. Variant effect classification showed 247 141 SNVs and 99 562 short indels having moderate or high impact on 11 267 protein-coding genes. On average, each genome contained heterozygous loss-of-function variants in 30 potentially embryonic lethal genes and 97 genes associated with developmental disorders. More than 50 inherited disorders and traits have been unravelled using the DBVDC variant catalogue, enabling genetic testing for breeding and diagnostics. This resource of annotated variants and their corresponding genotype frequencies constitutes a highly useful tool for the identification of potential variants causative for rare inherited disorders in dogs.