Variation of BMP3 contributes to dog breed skull diversity

Copy-cat evolution: Divergence and convergence within and between cat and dog breeds

Many domesticated species exhibit remarkable phenotypic diversity. In nature, selection produces not only divergence but also convergence when organisms experience similar selective pressures. Whether artificial selection during domestication also produces convergence has received little attention. Three-dimensional shape analysis of domestic cat and dog skulls demonstrated convergence at multiple levels. Most broadly, cats and dogs have both diversified greatly: equaling or exceeding the morphological disparity among all modern-day species of their respective families. Moreover, as a result of artificial selection, some breeds of these two phenotypically distinct species, evolutionarily separated for 50 My, have converged to such an extreme extent that they are more similar to each other than they are to many members of their own species or their ancestors, a phenomenon never previously observed in domesticated species. Remarkably, this convergence evolved not only between dogs and cats but also multiple times within each taxon.

Genomic analyses in Cavalier King Charles spaniels identify loci associated with clinical signs of Chiari-like malformation and Syringomyelia

BACKGROUND

Chiari-like malformations (CM) and syringomyelia (SM) are common in Cavalier King Charles spaniels (CKCS) leading to variable manifestations of pain and scratch. Inheritance studies suggest a polygenic mode of inheritance and association studies have identified loci associated with the presence of SM on MRI. Given the poor correlation of clinical signs of CMSM with MRI findings, we hypothesized that an association study with clinical signs as the phenotype could reveal new loci of interest. The objectives of this study were to perform genome-wide association studies on CKCS using SM and clinical sign phenotypes of pain and scratch and to use whole genome sequencing (WGS) to identify variants in regions of interest. We collected DNA on 174 CKCS. Owners completed questionnaires to establish the clinical pain and scratch phenotype and magnetic resonance imaging (MRI) was used to identify CM and SM (linear T2 hyperintensity greater than 2 mm in height) in all dogs. Dogs were genotyped using the Axiom K9 HD (710,000 snps) array. GWAS analyses were performed using GEMMA and categorical and quantitative approaches were used to define clinical phenotypes. Whole genome sequencing (WGS) was performed on an Illumina HiSeq 4000 high-throughput sequencing system.

RESULTS

There were no regions associated with SM presence. The presence of signs of pain and scratch was associated with a region on Canis familiaris autosome (CFA) 26 downstream of ZWINT, previously associated with skull changes in CKCS with SM, although genome-wide significance was not reached. Loci were also associated with quantitative pain and scratch scores on CFA 13, 2 and 38. There were 66 variants that segregated with phenotype including 2 missense variants that were predicted to have moderate effects on ZWINT function.

CONCLUSIONS

The identification of a locus on CFA26 using the clinical phenotype of pain and scratch that coincided with a locus identified in a morphological study provides strong support for this as a region of interest.

Analysis of canine gene constraint identifies new variants for orofacial clefts and stature

Dog breeding promotes within-group homogeneity through conformation to strict breed standards, while simultaneously driving between-group heterogeneity. There are over 350 recognized dog breeds that provide the foundation for investigating the genetic basis of phenotypic diversity. Typically, breed standard phenotypes such as stature, pelage, and craniofacial structure are analyzed through genetic association studies. However, such analyses are limited to assayed phenotypes only, leaving difficult-to-measure phenotypic subtleties easily overlooked. We investigated coding variation from over 2000 dogs, leading to discoveries of variants related to craniofacial morphology and stature. Breed-enriched variants were prioritized according to gene constraint, which was calculated using a mutation model derived from trinucleotide substitution probabilities. Among the newly found variants is a splice-acceptor variant in PDGFRA associated with bifid nose, a characteristic trait of Çatalburun dogs, implicating the gene's role in midline closure. Two additional LCORL variants, both associated with canine body size are also discovered: a frameshift that causes a premature stop in large breeds (>25 kg) and an intronic substitution found in small breeds (<10 kg), thus highlighting the importance of allelic heterogeneity in selection for breed traits. Most variants prioritized in this analysis are not associated with genomic signatures for breed differentiation, as these regions are enriched for constrained genes intolerant to nonsynonymous variation. This indicates trait selection in dogs is likely a balancing act between preserving essential gene functions and maximizing regulatory variation to drive phenotypic extremes.

Whole exome sequencing as a screening tool in dogs: A pilot study

Background

Whole-exome sequencing (WES) is used to selectively sequence all exons of protein-coding genes. WES is considered as a cost-effective and direct approach for identifying phenotype-associated variants in protein-coding regions and is as such situated between the traditional Sanger sequencing and whole genome sequencing (WGS). While WES is already widely used as a clinical tool in human and medical genetics, its use in veterinary medicine is currently restricted to research purposes. In this article, we aimed to provide baseline performance characteristics of a WES design to assess its suitability with future applications in veterinary clinical genetics in mind.

Methods

To assess the potential of WES in a clinical setting for dogs, 49 canine samples underwent capture, sequencing and analysis for the presence of 352 known phenotype-associated variants. The sequencing performance was compared for three types of variants, based on their size and location: single nucleotide variants (SNVs) inside exons, larger indel variants (≤20 bp) inside exons and intronic variants.

Results

On average, 85 % and 82 % of the exonic SNPs and larger variants were sequenced at a sequencing depth of ≥ 10x in the 49 samples, respectively. In the best performing sample, 94 % of the exonic SNPs were covered at least 10x, whereas in the worst performing sample, still 71 % of the exonic SNPs had an average sequencing depth of more than 10x.

Conclusion

To our knowledge, this is the first report that describes the performance of a research-intended WES design if it would be used in clinical genetics. This study found that WES demonstrated high efficacy in detecting variants located within target regions, including those that were not initially included in the design. However, the performance varied across different variants. The next steps would be the development of improved designs and settings to ameliorate the results.

Genetic and Anatomical Determinants of Olfaction in Dogs and Wild Canids

Understanding the anatomical and genetic basis of complex phenotypic traits has long been a challenge for biological research. Domestic dogs offer a compelling model as they demonstrate more phenotypic variation than any other vertebrate species. Dogs have been intensely selected for specific traits and abilities, directly or indirectly, over the past 15,000 years since their initial domestication from the gray wolf. Because olfaction plays a central role in critical tasks, such as the detection of drugs, diseases, and explosives, as well as human rescue, we compared relative olfactory capacity across dog breeds and assessed changes to the canine olfactory system to their direct ancestors, wolves, and coyotes. We conducted a cross-disciplinary survey of olfactory anatomy, olfactory receptor (OR) gene variation, and OR gene expression in domestic dogs. Through comparisons to their closest wild canid relatives, the gray wolf and coyote, we show that domestic dogs might have lost functional OR genes commensurate with a documented reduction in nasal morphology as an outcome of the domestication process prior to breed formation. Critically, within domestic dogs alone, we found no genetic or morphological profile shared among functional or genealogical breed groupings, such as scent hounds, that might indicate evidence of any human-directed selection for enhanced olfaction. Instead, our results suggest that superior scent detection dogs likely owe their success to advantageous behavioral traits and training rather than an "olfactory edge" provided by morphology or genes.

GRAMMAR-Lambda Delivers Efficient Understanding of the Genetic Basis for Head Size in Catfish

The shape of the skull plays a crucial role in the evolution and adaptation of species to their environments. In the case of aquaculture fish, the size of the head is also an important economic trait, as it is linked to fillet yield and ornamental value. This study applies our GRAMMAR-Lambda method to perform a genome-wide association study analysis on loci related to head size in catfish. Compared with traditional GWAS methods, the GRAMMAR-Lambda method offers higher computational efficiency, statistical power, and stability, especially in complex population structures. This research identifies many candidate genes closely related to cranial morphology in terms of head length, width, and depth in catfish, including bmpr1bb, fgfrl1b, nipbl, foxp2, and pax5, etc. Based on the results of gene-gene interaction analysis, we speculate that there may be frequent genetic interactions between chromosome 19 and chromosome 29 in bone development. Additionally, many candidate genes, gene families, and mechanisms (such as SOCE mechanisms) affecting skeletal development and morphology have been identified. These findings contribute to our understanding of the genetic architecture of head size and will support marker-assisted breeding in aquaculture, also reflecting the potential application of the GRAMMAR-Lambda method in genetic studies of complex traits.

Adaptation in the Alleyways: Candidate Genes Under Potential Selection in Urban Coyotes

In the context of evolutionary time, cities are an extremely recent development. Although our understanding of how urbanization alters ecosystems is well developed, empirical work examining the consequences of urbanization on adaptive evolution remains limited. To facilitate future work, we offer candidate genes for one of the most prominent urban carnivores across North America. The coyote (Canis latrans) is a highly adaptable carnivore distributed throughout urban and nonurban regions in North America. As such, the coyote can serve as a blueprint for understanding the various pathways by which urbanization can influence the genomes of wildlife via comparisons along urban-rural gradients, as well as between metropolitan areas. Given the close evolutionary relationship between coyotes and domestic dogs, we leverage the well-annotated dog genome and highly conserved mammalian genes from model species to outline how urbanization may alter coyote genotypes and shape coyote phenotypes. We identify variables that may alter selection pressure for urban coyotes and offer suggestions of candidate genes to explore. Specifically, we focus on pathways related to diet, health, behavior, cognition, and reproduction. In a rapidly urbanizing world, understanding how species cope and adapt to anthropogenic change can facilitate the persistence of, and coexistence with, these species.

Structure and scaling of the middle ear in domestic dog breeds

Although domestic dogs vary considerably in both body size and skull morphology, behavioural audiograms have previously been found to be similar in breeds as distinct as a Chihuahua and a St Bernard. In this study, we created micro-CT reconstructions of the middle ears and bony labyrinths from the skulls of 17 dog breeds, including both Chihuahua and St Bernard, plus a mongrel and a wolf. From these reconstructions, we measured middle ear cavity and ossicular volumes, eardrum and stapes footplate areas and bony labyrinth volumes. All of these ear structures scaled with skull size with negative allometry and generally correlated better with condylobasal length than with maximum or interaural skull widths. Larger dogs have larger ear structures in absolute terms: the volume of the St Bernard's middle ear cavity was 14 times that of the Chihuahua. The middle and inner ears are otherwise very similar in morphology, the ossicular structure being particularly well-conserved across breeds. The expectation that larger ear structures in larger dogs would translate into hearing ranges shifted towards lower frequencies is not consistent with the existing audiogram data. Assuming that the audiograms accurately reflect the hearing of the breeds in question, oversimplifications in existing models of middle ear function or limitations imposed by other parts of the auditory system may be responsible for this paradox.

Whole-Genome Sequencing Analyses Reveal the Evolution Mechanisms of Typical Biological Features of Decapterus maruadsi

Decapterus maruadsi is a typical representative of small pelagic fish characterized by fast growth rate, small body size, and high fecundity. It is a high-quality marine commercial fish with high nutritional value. However, the underlying genetics and genomics research focused on D. maruadsi is not comprehensive. Herein, a high-quality chromosome-level genome of a male D. maruadsi was assembled. The assembled genome length was 716.13 Mb with contig N50 of 19.70 Mb. Notably, we successfully anchored 95.73% contig sequences into 23 chromosomes with a total length of 685.54 Mb and a scaffold N50 of 30.77 Mb. A total of 22,716 protein-coding genes, 274.90 Mb repeat sequences, and 10,060 ncRNAs were predicted, among which 22,037 (97%) genes were successfully functionally annotated. The comparative genome analysis identified 459 unique, 73 expanded, and 52 contracted gene families. Moreover, 2804 genes were identified as candidates for positive selection, of which some that were related to the growth and development of bone, muscle, cardioid, and ovaries, such as some members of the TGF-β superfamily, were likely involved in the evolution of typical biological features in D. maruadsi. The study provides an accurate and complete chromosome-level reference genome for further genetic conservation, genomic-assisted breeding, and adaptive evolution research for D. maruadsi.

Exploring the origin of a unique mutant allele in twin-tail goldfish using CRISPR/Cas9 mutants

Artificial selection has been widely applied to genetically fix rare phenotypic features in ornamental domesticated animals. For many of these animals, the mutated loci and alleles underlying rare phenotypes are known. However, few studies have explored whether these rare genetic mutations might have been fixed due to competition among related mutated alleles or if the fixation occurred due to contingent stochastic events. Here, we performed genetic crossing with twin-tail ornamental goldfish and CRISPR/Cas9-mutated goldfish to investigate why only a single mutated allele-chdS with a E127X stop codon (also called chdAE127X)-gives rise to the twin-tail phenotype in the modern domesticated goldfish population. Two closely related chdS mutants were generated with CRISPR/Cas9 and compared with the E127X allele in F2 and F3 generations. Both of the CRISPR/Cas9-generated alleles were equivalent to the E127X allele in terms of penetrance/expressivity of the twin-tail phenotype and viability of carriers. These findings indicate that multiple truncating mutations could have produced viable twin-tail goldfish. Therefore, the absence of polymorphic alleles for the twin-tail phenotype in modern goldfish likely stems from stochastic elimination or a lack of competing alleles in the common ancestor. Our study is the first experimental comparison of a singular domestication-derived allele with CRISPR/Cas9-generated alleles to understand how genetic fixation of a unique genotype and phenotype may have occurred. Thus, our work may provide a conceptual framework for future investigations of rare evolutionary events in domesticated animals.

Genome-wide identification and structural analysis of the BMP gene family in Triplophysa dalaica

BACKGROUND

Bone morphogenetic proteins (BMPs) are part of the transforming growth factor beta (TGF-β) superfamily and play crucial roles in bone development, as well as in the formation and maintenance of various organs. Triplophysa dalaica, a small loach fish that primarily inhabits relatively high elevations and cooler water bodies, was the focus of this study. Understanding the function of BMP genes during the morphogenesis of T. dalaica helps to clarify the mechanisms of its evolution and serves as a reference for the study of BMP genes in other bony fishes. The data for the T. dalaica transcriptome and genome used in this investigation were derived from the outcomes of our laboratory sequencing.

RESULTS

This study identified a total of 26 BMP genes, all of which, except for BMP1, possess similar TGF-β structural domains. We conducted an analysis of these 26 BMP genes, examining their physicochemical properties, subcellular localization, phylogenetic relationships, covariance within and among species, chromosomal localization, gene structure, conserved motifs, conserved structural domains, and expression patterns. Our findings indicated that three BMP genes were associated with unstable proteins, while 11 BMP genes were located within the extracellular matrix. Furthermore, some BMP genes were duplicated, with the majority being enriched in the GO:0008083 pathway, which is related to growth factor activity. It was hypothesized that genes within the BMP1/3/11/15 subgroup (Group I) play a significant role in the growth and development of T. dalaica. By analyzing the expression patterns of proteins in nine tissues (gonad, kidney, gill, spleen, brain, liver, fin, heart, and muscle), we found that BMP genes play diverse regulatory roles during different stages of growth and development and exhibit characteristics of division of labor.

CONCLUSIONS

This study contributes to a deeper understanding of BMP gene family member expression patterns in high-altitude, high-salinity environments and provides valuable insights for future research on the BMP gene family in bony fishes.

Skeletal Class III phenotype: Link between animal models and human genetics: A scoping review

This study aimed to identify evidence from animal studies examining genetic variants underlying maxillomandibular discrepancies resulting in a skeletal Class III (SCIII) malocclusion phenotype. Following the Manual for Evidence Synthesis of the JBI and the PRISMA extension for scoping reviews, a participant, concept, context question was formulated and systematic searches were executed in the PubMed, Scopus, WOS, Scielo, Open Gray, and Mednar databases. Of the 779 identified studies, 13 met the selection criteria and were included in the data extraction. The SCIII malocclusion phenotype was described as mandibular prognathism in the Danio rerio, Dicentrarchus labrax, and Equus africanus asinus models; and as maxillary deficiency in the Felis silvestris catus, Canis familiaris, Salmo trutta, and Mus musculus models. The identified genetic variants highlight the significance of BMP and TGF-β signaling. Their regulatory pathways and genetic interactions link them to cellular bone regulation events, particularly ossification regulation of postnatal cranial synchondroses. In conclusion, twenty genetic variants associated with the skeletal SCIII malocclusion phenotype were identified in animal models. Their interactions and regulatory pathways corroborate the role of these variants in bone growth, differentiation events, and ossification regulation of postnatal cranial synchondroses.

Important Role of the Ihh Signaling Pathway in Initiating Early Cranial Remodeling and Morphological Specialization in Cromileptes altivelis

In this study, we identified the important contribution of frontal bone remodeling in shaping the 'sunken head and humpback' appearance in C. altivelis. Our investigation identified a developmental milestone at a total length of 5-6 cm, making the onset of its morphologic specialization in this species. A comparative analysis with closely related species reveals heightened activity in the frontal osteoblasts of the humpback grouper, potentially providing a physiological basis for its remodeling. Furthermore, our findings highlight that a significant upregulation in the expression levels of Ihhb, Ptch1, and Gli2a genes was seen in C. altivelis within the specified developmental stage, indicating an important involvement of the Ihhb-Ptch1-Gli2a signaling pathway in initiating the morphological specialization. We hypothesized that Ihh signaling could be attributed to shifts in mechanical stress, resulting from muscle traction on the frontal bone due to changes in swimming patterns during development. This study not only offers significant insights into unraveling the molecular mechanisms that govern phenotypic specialization and ecological adaptations in the humpback grouper but also serves as a valuable reference for studies on fishes with a controversial morphology and molecular phylogeny.

Outcome of Surgical Extraction of Linguoverted Deciduous Mandibular Canine Teeth, Performed as an Interceptive Orthodontic Procedure in Puppies, on Permanent Occlusion

Linguoversion of deciduous mandibular canine teeth can be a painful condition, interferes with the development and growth of the jaws, and potentially leads to further malocclusions affecting permanent dentition. Extraction of linguoverted deciduous mandibular canines is considered an interceptive orthodontic procedure that would allow unimpeded development of the jaws and permanent teeth. This study assessed clinical records of 124 dogs that had linguoverted deciduous mandibular canine teeth surgically extracted between October 2010 and September 2019 in a veterinary dental referral clinic. Seventy-seven cases fulfilled the study criteria. Fifty-one percent of these patients required further orthodontic treatment of the permanent occlusion and forty-nine percent demonstrated atraumatic permanent occlusion. The study found no correlation of the outcome with age at the time of surgery. The class of malocclusion (class 1 or class 2) at the time of surgery was also not associated with the outcome.

Evolutionary mechanisms modulating the mammalian skull development

Mammals possess impressive craniofacial variation that mirrors their adaptation to diverse ecological niches, feeding behaviour, physiology and overall lifestyle. The spectrum of craniofacial geometries is established mainly during embryonic development. The formation of the head represents a sequence of events regulated on genomic, molecular, cellular and tissue level, with each step taking place under tight spatio-temporal control. Even minor variations in timing, position or concentration of the molecular drivers and the resulting events can affect the final shape, size and position of the skeletal elements and the geometry of the head. Our knowledge of craniofacial development increased substantially in the last decades, mainly due to research using conventional vertebrate model organisms. However, how developmental differences in head formation arise specifically within mammals remains largely unexplored. This review highlights three evolutionary mechanisms acknowledged to modify ontogenesis: heterochrony, heterotopy and heterometry. We present recent research that links changes in developmental timing, spatial organization or gene expression levels to the acquisition of species-specific skull morphologies. We highlight how these evolutionary modifications occur on the level of the genes, molecules and cellular processes, and alter conserved developmental programmes to generate a broad spectrum of skull shapes characteristic of the class Mammalia. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

The LASSIE MPS panel: Predicting externally visible traits in dogs for forensic purposes

Predicting the outward appearance of dogs via their DNA, also known as Canine DNA Phenotyping, is a young, emerging field of research in forensic genetics. The few previous studies published in this respect were restricted to the consecutive analysis of single DNA markers, a process that is time- and sample-consuming and therefore not a viable option for limited forensic specimens. Here, we report on the development and evaluation of a Massively Parallel Sequencing (MPS) based molecular genetic assay, the LASSIE MPS Panel. This panel aims to predict externally visible as well as skeletal traits, which include coat color, coat pattern, coat structure, tail morphology, skull shape, ear shape, eye color and body size from DNA using 44 genetic markers in a single molecular genetic assay. A biostatistical naïve Bayes classification approach was applied to identify the most informative marker combinations for predicting phenotypes. Overall, the predictive performance was characterized by a very high classification success for some of the trait categories, and high to moderate success for others. The performance of the developed predictive framework was further evaluated using blind samples from three randomly selected dog individuals, whose appearance was well predicted.

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.

That brachycephalic look: Infant-like facial appearance in short-muzzled dog breeds

Brachycephalic dog breeds are highly popular, yet their conformation-related disorders represent a major welfare concern. It has been suggested that the current popularity of such breeds can be explained by their cute, infant-like facial appearances. The concept of 'kindchenschema' refers to the observation that certain physical features of infant humans and other animals can automatically stimulate positive and nurturant feelings in adult observers. But the proposal that brachycephalic dogs possess heightened 'kindchenschema' facial features, even into adulthood, has never been formally investigated. Here, we hypothesised that relative muzzle shortening across a range of breeds would be associated with known 'kindchenschema' facial features, including a relatively larger forehead, larger eyes and smaller nose. Relative fronto-facial feature sizes in exemplar photographs of adult dogs from 42 popular breeds were measured and associated with existing data on the relative muzzle length and height-at-withers of the same breeds. Our results show that, in adulthood, shorter-muzzled breeds have relatively larger (taller) foreheads and relatively larger eyes (i.e. area of exposed eyeball relative to overall face area) than longer-muzzled breeds, and that this effect is independent of breed size. In sum, brachycephalic dog breeds do show exaggeration of some, but not all, known fronto-facial 'kindchenschema' features, and this may well contribute to their apparently cute appearance and to their current popularity as companion animals. We conclude that the challenge of addressing conformation-related disorders in companion dogs needs to take account of the cute, 'kindchenschema' looks that many owners are likely to be attracted to.

Across-breed genetic investigation of canine hip dysplasia, elbow dysplasia, and anterior cruciate ligament rupture using whole-genome sequencing

Here, we report the use of genome-wide association study (GWAS) for the analysis of canine whole-genome sequencing (WGS) repository data using breed phenotypes. Single-nucleotide polymorphisms (SNPs) were called from WGS data from 648 dogs that included 119 breeds from the Dog10K Genomes Project. Next, we assigned breed phenotypes for hip dysplasia (Orthopedic Foundation for Animals (OFA) HD, n = 230 dogs from 27 breeds; hospital HD, n = 279 dogs from 38 breeds), elbow dysplasia (ED, n = 230 dogs from 27 breeds), and anterior cruciate ligament rupture (ACL rupture, n = 279 dogs from 38 breeds), the three most important canine spontaneous complex orthopedic diseases. Substantial morbidity is common with these diseases. Previous within- and between-breed GWAS for HD, ED, and ACL rupture using array SNPs have identified disease-associated loci. Individual disease phenotypes are lacking in repository data. There is a critical knowledge gap regarding the optimal approach to undertake categorical GWAS without individual phenotypes. We considered four GWAS approaches: a classical linear mixed model, a haplotype-based model, a binary case-control model, and a weighted least squares model using SNP average allelic frequency. We found that categorical GWAS was able to validate HD candidate loci. Additionally, we discovered novel candidate loci and genes for all three diseases, including FBX025, IL1A, IL1B, COL27A1, SPRED2 (HD), UGDH, FAF1 (ED), TGIF2 (ED & ACL rupture), and IL22, IL26, CSMD1, LDHA, and TNS1 (ACL rupture). Therefore, categorical GWAS of ancestral dog populations may contribute to the understanding of any disease for which breed epidemiological risk data are available, including diseases for which GWAS has not been performed and candidate loci remain elusive.

Zrsr2 Is Essential for the Embryonic Development and Splicing of Minor Introns in RNA and Protein Processing Genes in Zebrafish

ZRSR2 (zinc finger CCCH-type, RNA binding motif and serine/arginine rich 2) is an essential splicing factor involved in 3' splice-site recognition as a component of both the major and minor spliceosomes that mediate the splicing of U2-type (major) and U12-type (minor) introns, respectively. Studies of ZRSR2-depleted cell lines and ZRSR2-mutated patient samples revealed its essential role in the U12-dependent minor spliceosome. However, the role of ZRSR2 during embryonic development is not clear, as its function is compensated for by Zrsr1 in mice. Here, we utilized the zebrafish model to investigate the role of zrsr2 during embryonic development. Using CRISPR/Cas9 technology, we generated a zrsr2-knockout zebrafish line, termed zrsr2hg129/hg129 (p.Trp167Argfs*9) and examined embryo development in the homozygous mutant embryos. zrsr2hg129/hg129 embryos displayed multiple developmental defects starting at 4 days post fertilization (dpf) and died after 8 dpf, suggesting that proper Zrsr2 function is required during embryonic development. The global transcriptome analysis of 3 dpf zrsr2hg129/hg129 embryos revealed that the loss of Zrsr2 results in the downregulation of essential metabolic pathways and the aberrant retention of minor introns in about one-third of all minor intron-containing genes in zebrafish. Overall, our study has demonstrated that the role of Zrsr2 as a component of the minor spliceosome is conserved and critical for proper embryonic development in zebrafish.

Application of zebrafish in the study of craniomaxillofacial developmental anomalies

Craniomaxillofacial developmental anomalies are one of the most prevalent congenital defects worldwide and could result from any disruption of normal development processes, which is generally influenced by interactions between genes and the environment. Currently, with the advances in genetic screening strategies, an increasing number of novel variants and their roles in orofacial diseases have been explored. Zebrafish is recognized as a powerful animal model, and its homologous genes and similar oral structure and development process provide an ideal platform for studying the contributions of genetic and environmental factors to human craniofacial malformations. Here, we reviewed zebrafish models for the study of craniomaxillofacial developmental anomalies, such as human nonsyndromic cleft lip with or without an affected palate and jaw and tooth developmental anomalies. Due to its potential for gene expression and regulation research, zebrafish may provide new perspectives for understanding craniomaxillofacial diseaseand its treatment.

The Shape of the Nasal Cavity and Adaptations to Sniffing in the Dog (Canis familiaris) Compared to Other Domesticated Mammals: A Review Article

Dogs are a good starting point for the description and anatomical analysis of turbinates of the nose. This work aimed at summing up the state of knowledge on the shape of the nasal cavity and airflow in these domestic animals and dealt with the brachycephalic syndrome (BOAS) and anatomical changes in the initial airway area in dogs with a short and widened skull. As a result of artificial selection and breeding concepts, the dog population grew very quickly. Modern dog breeds are characterized by a great variety of their anatomical shape. Craniological changes also had a significant impact on the structure and physiology of the respiratory system in mammals. The shape of the nasal cavity is particularly distinctive in dogs. Numerous studies have established that dogs and their olfactory ability are of great importance in searching for lost people, detecting explosives or drugs as well as signaling disease in the human body. The manuscript describes the structure of the initial part of the respiratory system, including the nasal turbinates, and compares representatives of various animal species. It provides information on the anatomy of brachycephalic dogs and BOAS. The studies suggest that further characterization and studies of nasal turbinates and their hypertrophy are important.

BMP3 is a novel locus involved in the causality of ocular coloboma

Coloboma, a congenital disorder characterized by gaps in ocular tissues, is caused when the choroid fissure fails to close during embryonic development. Several loci have been associated with coloboma, but these represent less than 40% of those that are involved with this disease. Here, we describe a novel coloboma-causing locus, BMP3. Whole exome sequencing and Sanger sequencing of patients with coloboma identified three variants in BMP3, two of which are predicted to be disease causing. Consistent with this, bmp3 mutant zebrafish have aberrant fissure closure. bmp3 is expressed in the ventral head mesenchyme and regulates phosphorylated Smad3 in a population of cells adjacent to the choroid fissure. Furthermore, mutations in bmp3 sensitize embryos to Smad3 inhibitor treatment resulting in open choroid fissures. Micro CT scans and Alcian blue staining of zebrafish demonstrate that mutations in bmp3 cause midface hypoplasia, suggesting that bmp3 regulates cranial neural crest cells. Consistent with this, we see active Smad3 in a population of periocular neural crest cells, and bmp3 mutant zebrafish have reduced neural crest cells in the choroid fissure. Taken together, these data suggest that Bmp3 controls Smad3 phosphorylation in neural crest cells to regulate early craniofacial and ocular development.

Estimation of periodontal pocket surface area in small to medium dogs: a proof-of-concept study

Background

Periodontal disease is the most common dental disease in dogs. Although the systemic effects of periodontal disease have not been clarified in veterinary science, it is necessary to evaluate the effects of periodontal disease in clinical trials in the future. There have been a few clinical attempts made, however, to assess the severity of periodontal inflammation and its impact on the systemic health of dogs. Meanwhile, in the field of dentistry for humans, the periodontal inflamed surface area (PISA) and periodontal epithelial surface area (PESA) have been used to quantitatively assess the degree of periodontal disease affecting a single tooth as well as the overall extent of periodontitis. Recent studies have also suggested the use of these assessments to examine the relationship between periodontal inflammation and systemic health.

Results

The estimation formula for a dog’s periodontal pocket surface area (PPSA), an alternative to PISA and PESA in humans, was established using body weight and periodontal pocket depth. Actual values were measured using extracted teeth from various dog breeds and sizes (2.3–25.0 kg of body weight) to obtain universal regression equations for PPSA. Altogether, 625 teeth from 73 dogs of 16 breeds were extracted and subsequently analyzed for morphological information.

PPSA was measured in 61 dogs of 10 breeds with periodontal disease using the established estimation formulas, and the correlation between PPSA and preoperative blood chemistry data was analyzed accordingly. A strong correlation was found between PPSA and serum globulin (r = 0.71) while moderate correlations were found for C-reactive protein (r = 0.54) and serum albumin (r = -0.51).

Conclusions

Estimation formulas using body weight and the 6-point probing depth were established for determining PPSA. Direct correlations between PPSA and several blood test results were observed in the study sample. Taken together, these results suggest that PPSA could be useful for evaluating the effects of periodontitis on systemic conditions in dogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-03116-0.

Brachygnathia Inferior in Cloned Dogs Is Possibly Correlated with Variants of Wnt Signaling Pathway Initiators

Abnormalities in animals cloned via somatic cell nuclear transfer (SCNT) have been reported. In this study, to produce bomb-sniffing dogs, we successfully cloned four healthy dogs through SCNT using the same donor genome from the skin of a male German shepherd old dog. Veterinary diagnosis (X-ray/3D-CT imaging) revealed that two cloned dogs showed normal phenotypes, whereas the others showed abnormal shortening of the mandible (brachygnathia inferior) at 1 month after birth, even though they were cloned under the same conditions except for the oocyte source. Therefore, we aimed to determine the genetic cause of brachygnathia inferior in these cloned dogs. To determine the genetic defects related to brachygnathia inferior, we performed karyotyping and whole-genome sequencing (WGS) for identifying small genetic alterations in the genome, such as single-nucleotide variations or frameshifts. There were no chromosomal numerical abnormalities in all cloned dogs. However, WGS analysis revealed variants of Wnt signaling pathway initiators (WNT5B, DVL2, DACT1, ARRB2, FZD 4/8) and cadherin (CDH11, CDH1like) in cloned dogs with brachygnathia inferior. In conclusion, this study proposes that brachygnathia inferior in cloned dogs may be associated with variants in initiators and/or regulators of the Wnt/cadherin signaling pathway.

Forensic DNA phenotyping: Canis familiaris breed classification and skeletal phenotype prediction using functionally significant skeletal SNPs and indels

This review highlights a novel application of breed identification and prediction of skeletal traits in forensic investigations using canine DNA evidence. Currently, genotyping methods used for canine breed classification involve the application of highly polymorphic short tandem repeats in addition to larger commercially available SNP arrays. Both applications face technical challenges. An additional approach to breed identification could be through genotyping SNPs and indels that characterise the array of skeletal differences displayed across domestic dog populations. Research has shown that a small number of genetic variants of large effect drive differences in skeletal phenotypes among domestic dog breeds. This feature makes functionally significant canine skeletal variants a cost-effective target for forensic investigators to classify individuals according to their breed. Further analysis of these skeletal variants would enable the prediction of external appearance. To date, functionally significant genes with genetic variants associated with differences in size, bulk, skull shape, ear shape, limb length, digit type, and tail morphology have been uncovered. Recommendations of a cost-effective genotyping method that can be readily designed and applied by forensic investigators have been given. Further advances to improve the field of canine skeletal forensic DNA phenotyping include the refinement of phenotyping methods, further biological validation of the skeletal genetic variants and establishing a publicly available database for storage of allele frequencies of the skeletal genetic variants in the wider domestic dog population.

Nasal hyperkeratosis in Griffon breeds: Clinical, histopathological features and the prevalence in the Swedish population compared to a control group and other brachycephalic breeds

BACKGROUND

In the Griffon breeds (GB) nasal hyperkeratosis is common and develops already in early adulthood. Breed-related features and prevalence have not previously been documented.

HYPOTHESIS/OBJECTIVES

To describe clinical and histopathological features of nasal hyperkeratosis in GB and to document the prevalence.

MATERIALS AND METHODS

Seven GB dogs with nasal hyperkeratosis were examined. Three histopathological samples were analysed. Owners of 107 GB and 493 control dogs completed a questionnaire distributed via social media.

RESULTS

Typical features of nasal hyperkeratosis in GB included varying degrees of dry, firm, excessive proliferation of keratin, affecting the dorsal or dorsolateral aspect of the planum nasale. Histopathology was characterized by severe, lamellar orthokeratotic and focal parakeratotic hyperkeratosis and multiple small serum lakes. Thirty-four of 107 GB dogs (31.8%) and 65 of 493 (13.2%) control dogs had varying degree of nasal hyperkeratosis. No sex predisposition was noted. Median age of onset was 3 years for GB, similar to brachycephalic control dogs whereas non-brachycephalic control dogs had a significantly later age of onset (p = 0.0053).

CONCLUSIONS AND CLINICAL IMPORTANCE

Idiopathic nasal hyperkeratosis is very common in GB dogs and other brachycephalic breeds with nearly one third being affected, often already a young age.

Ophthalmic disorders in a referral population of seven breeds of brachycephalic dogs: 970 cases (2008-2017)

OBJECTIVE

To evaluate the frequency of ophthalmic disorders in 7 brachycephalic dog breeds referred to an academic veterinary ophthalmology service.

ANIMALS

970 client-owned dogs of 7 brachycephalic breeds that were evaluated by the ophthalmology service in a veterinary teaching hospital from January 2008 through December 2017.

PROCEDURES

Medical records of 7 brachycephalic breeds (ie, Boston Terriers, English Bulldogs, French Bulldogs, Lhasa Apsos, Pekingese, Pugs, and Shih Tzus) were reviewed to collect data regarding patient signalment, ophthalmic diagnoses, affected eyes, and number and dates of visits.

RESULTS

Median age at the first examination was 7 years (range, 23 days to 22 years). The number of dogs seen for a first examination increased with age. Corneal ulcers, keratoconjunctivitis sicca, corneal pigmentation, immature cataracts, and uveitis were each diagnosed in ≥ 100 dogs and represented 40.4% (1,161/2,873) of all diagnoses. On the basis of anatomic location, 66.3% (1,905/2,873) of all disorders were located in either the cornea (1,014/2,873 [35.2%]) or adnexa (891/2,873 [31%]). There was a significant difference in breed proportion in the study population; of the 7 breeds studied, Shih Tzus (34.3% [333/970]), Pugs (20.8% [202/970]), and Boston Terriers (16.6% [161/970]) were the most prevalent breeds. The frequency of some diseases within the referral population was associated with breed.

CONCLUSIONS AND CLINICAL RELEVANCE

Findings suggested that the most prevalent disorders for the brachycephalic breeds in this ophthalmic referral population were corneal ulcers, keratoconjunctivitis sicca, corneal pigmentation, immature cataracts, and uveitis. Although all dogs shared brachycephalic features, the frequency of specific ophthalmic diseases varied between breeds.

Complex genetic architecture of three-dimensional craniofacial shape variation in domestic pigeons

Deciphering the genetic basis of vertebrate craniofacial variation is a longstanding biological problem with broad implications in evolution, development, and human pathology. One of the most stunning examples of craniofacial diversification is the adaptive radiation of birds, in which the beak serves essential roles in virtually every aspect of their life histories. The domestic pigeon (Columba livia) provides an exceptional opportunity to study the genetic underpinnings of craniofacial variation because of its unique balance of experimental accessibility and extraordinary phenotypic diversity within a single species. We used traditional and geometric morphometrics to quantify craniofacial variation in an F2 laboratory cross derived from the straight-beaked Pomeranian Pouter and curved-beak Scandaroon pigeon breeds. Using a combination of genome-wide quantitative trait locus scans and multi-locus modeling, we identified a set of genetic loci associated with complex shape variation in the craniofacial skeleton, including beak shape, braincase shape, and mandible shape. Some of these loci control coordinated changes between different structures, while others explain variation in the size and shape of specific skull and jaw regions. We find that in domestic pigeons, a complex blend of both independent and coupled genetic effects underlie three-dimensional craniofacial morphology.

Exceptional Changes in Skeletal Anatomy under Domestication: The Case of Brachycephaly

"Brachycephaly" is generally considered a phenotype in which the facial part of the head is pronouncedly shortened. While brachycephaly is characteristic for some domestic varieties and breeds (e.g., Bulldog, Persian cat, Niata cattle, Anglo-Nubian goat, Middle White pig), this phenotype can also be considered pathological. Despite the superficially similar appearance of "brachycephaly" in such varieties and breeds, closer examination reveals that "brachycephaly" includes a variety of different cranial modifications with likely different genetic and developmental underpinnings and related with specific breed histories. We review the various definitions and characteristics associated with brachycephaly in different domesticated species. We discern different types of brachycephaly ("bulldog-type," "katantognathic," and "allometric" brachycephaly) and discuss morphological conditions related to brachycephaly, including diseases (e.g., brachycephalic airway obstructive syndrome). Further, we examine the complex underlying genetic and developmental processes and the culturally and developmentally related reasons why brachycephalic varieties may or may not be prevalent in certain domesticated species. Knowledge on patterns and mechanisms associated with brachycephaly is relevant for domestication research, veterinary and human medicine, as well as evolutionary biology, and highlights the profound influence of artificial selection by humans on animal morphology, evolution, and welfare.

Towards Forensic DNA Phenotyping for Predicting Visible Traits in Dogs

The popularity of dogs as human companions explains why these pets regularly come into focus in forensic cases such as bite attacks or accidents. Canine evidence, e.g., dog hairs, can also act as a link between the victim and suspect in a crime case due to the close contact between dogs and their owners. In line with human DNA identification, dog individualization from crime scene evidence is mainly based on the analysis of short tandem repeat (STR) markers. However, when the DNA profile does not match a reference, additional information regarding the appearance of the dog may provide substantial intelligence value. Key features of the dog's appearance, such as the body size and coat colour are well-recognizable and easy to describe even to non-dog experts, including most investigating officers and eyewitnesses. Therefore, it is reasonable to complement eyewitnesses' testimonies with externally visible traits predicted from associated canine DNA samples. Here, the feasibility and suitability of canine DNA phenotyping is explored from scratch in the form of a proof of concept study. To predict the overall appearance of an unknown dog from its DNA as accurately as possible, the following six traits were chosen: (1) coat colour, (2) coat pattern, (3) coat structure, (4) body size, (5) ear shape, and (6) tail length. A total of 21 genetic markers known for high predicting values for these traits were selected from previously published datasets, comprising 15 SNPs and six INDELS. Three of them belonged to SINE insertions. The experiments were designed in three phases. In the first two stages, the performance of the markers was tested on DNA samples from dogs with well-documented physical characteristics from different breeds. The final blind test, including dogs with initially withheld appearance information, showed that the majority of the selected markers allowed to develop composite sketches, providing a realistic impression of the tested dogs. We regard this study as the first attempt to evaluate the possibilities and limitations of forensic canine DNA phenotyping.

Multi-omics approach identifies germline regulatory variants associated with hematopoietic malignancies in retriever dog breeds

Histiocytic sarcoma is an aggressive hematopoietic malignancy of mature tissue histiocytes with a poorly understood etiology in humans. A histologically and clinically similar counterpart affects flat-coated retrievers (FCRs) at unusually high frequency, with 20% developing the lethal disease. The similar clinical presentation combined with the closed population structure of dogs, leading to high genetic homogeneity, makes dogs an excellent model for genetic studies of cancer susceptibility. To determine the genetic risk factors underlying histiocytic sarcoma in FCRs, we conducted multiple genome-wide association studies (GWASs), identifying two loci that confer significant risk on canine chromosomes (CFA) 5 (P_wald = 4.83x10^-9) and 19 (P_wald = 2.25x10^-7). We subsequently undertook a multi-omics approach that has been largely unexplored in the canine model to interrogate these regions, generating whole genome, transcriptome, and chromatin immunoprecipitation sequencing. These data highlight the PI3K pathway gene PIK3R6 on CFA5, and proximal candidate regulatory variants that are strongly associated with histiocytic sarcoma and predicted to impact transcription factor binding. The CFA5 association colocalizes with susceptibility loci for two hematopoietic malignancies, hemangiosarcoma and B-cell lymphoma, in the closely related golden retriever breed, revealing the risk contribution this single locus makes to multiple hematological cancers. By comparison, the CFA19 locus is unique to the FCR and harbors risk alleles associated with upregulation of TNFAIP6, which itself affects cell migration and metastasis. Together, these loci explain ~35% of disease risk, an exceptionally high value that demonstrates the advantages of domestic dogs for complex trait mapping and genetic studies of cancer susceptibility.

We have identified two regions of the canine genome that explain a striking 35% of risk for developing histiocytic sarcoma in FCRs. The disease is uniformly lethal, affects 20% of FCRs, and parallels a cancer of the same name in humans. Both regions harbor genes involved in cell migration and cancer-related pathways. The first includes variants in regulatory regions at the tumor suppressor PIK3R6 locus that are strongly associated with histiocytic sarcoma and likely confer risk for other hematopoietic cancers. FCRs with risk alleles at the second locus demonstrate increased expression of TNFAIP6, which correlates with poor prognosis in multiple human cancers. In identifying genomic differences between affected and unaffected dogs, we advance our understanding of both canine and human health biology and set the stage for the development of diagnostic and therapeutic strategies.

Copy number variation underlies complex phenotypes in domestic dog breeds and other canids

Extreme phenotypic diversity, a history of artificial selection, and socioeconomic value make domestic dog breeds a compelling subject for genomic research. Copy number variation (CNV) is known to account for a significant part of inter-individual genomic diversity in other systems. However, a comprehensive genome-wide study of structural variation as it relates to breed-specific phenotypes is lacking. We have generated whole genome CNV maps for more than 300 canids. Our data set extends the canine structural variation landscape to more than 100 dog breeds, including novel variants that cannot be assessed using microarray technologies. We have taken advantage of this data set to perform the first CNV-based genome-wide association study (GWAS) in canids. We identify 96 loci that display copy number differences across breeds, which are statistically associated with a previously compiled set of breed-specific morphometrics and disease susceptibilities. Among these, we highlight the discovery of a long-range interaction involving a CNV near MED13L and TBX3, which could influence breed standard height. Integration of the CNVs with chromatin interactions, long noncoding RNA expression, and single nucleotide variation highlights a subset of specific loci and genes with potential functional relevance and the prospect to explain trait variation between dog breeds.

Dog breeds and body conformations with predisposition to osteosarcoma in the UK: a case-control study

Background

Osteosarcoma is an aggressive and painful bone neoplasm in dogs. Previous studies have reported epidemiological associations suggesting that large body mass, long bone length and the genetics of certain breeds including the Rottweiler are associated with elevated osteosarcoma risk. However, these studies were often limited by selection bias and confounding factors, and have rarely offered insights into breed-associated protection for osteosarcoma. The current study includes 1756 appendicular and axial osteosarcoma cases presenting to VPG Histology (Bristol, UK) compared against a control population of 905,211 dogs without osteosarcoma from primary care electronic patient records in the VetCompass™ dataset.

Methods and study design

Retrospective, case-control study. Multivariable logistic regression analysis explored associations between demographic risk factors (including breed, chondrodystrophy, age, sex/neuter status, skull-shape, and body mass) and osteosarcoma of all anatomical sites.

Results

We identified several breeds with increased and reduced odds of osteosarcoma. At highest risk were the Rottweiler and Great Dane, with > 10 times the odds of osteosarcoma compared with crossbreds, and the Rhodesian Ridgeback, which has not featured in previous lists of at-risk breeds for osteosarcoma, and had an odds ratio of 11.31 (95% confidence interval 7.37–17.35). Breeds at lowest risk of osteosarcoma (protected breeds) included the Bichon Frise, the French Bulldog and the Cavalier King Charles Spaniel, all with odd ratios of less than 0.30 compared with crossbreds. Body mass was strongly associated with osteosarcoma risk; dogs over 40 kg exhibited osteosarcoma odds of 45.44 (95% confidence interval 33.74–61.20) compared with dogs less than 10 kg. Chondrodystrophic breeds had an osteosarcoma odds ratio of 0.13 (95% confidence interval 0.11–0.16) compared with non-chondrodystrophic breeds.

Conclusions

This study provides evidence of strong breed-associated osteosarcoma risk and protection, suggesting a genetic basis for osteosarcoma pathogenesis. It highlights that breeds selected for long legs/large body mass are generally overrepresented amongst at-risk breeds, whilst those selected for short leg length/small body mass are generally protected. These findings could inform genetic studies to identify osteosarcoma risk alleles in canines and humans; as well as increasing awareness amongst veterinarians and owners, resulting in improved breeding practices and clinical management of osteosarcoma in dogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40575-021-00100-7.

Canine DVL2 variant contributes to brachycephalic phenotype and caudal vertebral anomalies

A frameshift deletion variant in the Wnt pathway gene dishevelled 2 (DVL2) is associated with a truncated, kinked tail (“screw tail”) in English Bulldogs, French Bulldogs and Boston Terriers. These breeds are also characterized by distinctive morphological traits, including a wide head, flat face and short-limbed dwarfism, which are characteristic of Robinow syndrome in humans, caused by defects in genes such as DVL1 and DVL3. Based on these phenotypic and genetic similarities, it has previously been hypothesized that the canine DVL2 variant results in a syndromic phenotype called the Robinow-like syndrome. In our study, we investigated the distribution of the DVL2 variant in 1954 dogs from 15 breeds, identifying breeds with allele variation and enabling the dissection of the genotype–phenotype correlation for the first time. With CT examinations in American Staffordshire Terriers, we confirmed that the DVL2 allele is associated with caudal vertebral malformations and a brachycephalic phenotype. We also hypothesize that the variant may be linked to additional health conditions, including brachycephalic obstructive airway syndrome and congenital heart defects. Altogether, our study strengthens the role of DVL2 as one of the contributors to the “bulldog type” morphology and features on the spectrum of human Robinow syndrome.

Genetic factors contributing to skeletal class III malocclusion: a systematic review and meta-analysis

OBJECTIVES

The present systematic review aims to report and critically assess the findings of the available scientific evidence from genetic association studies examining the genetic variants underlying skeletal class III malocclusion and its sub-phenotypes.

MATERIAL AND METHODS

A pre-piloted protocol was registered and followed. The PubMed, Scopus, WOS, Cochrane Library, Gray Open literature, and CADTH databases were explored for genetic association studies following PICOS-based selection criteria. The research was reported in accordance with PRISMA statement and HuGE guidelines. The Q-genie tool was applied to assess the quality of genetic studies. Meta-analysis of genetic association studies was done by means of Meta-Genyo tool.

RESULTS

A total of 8258 articles were retrieved, of which 22 were selected for in-depth analysis. Most of the studies did not differentiate between sub-phenotypes, and the cohorts were heterogeneous regarding ethnicity. Four to five principal components of class III malocclusion explained the phenotypic variation, and gene variants at MYO1H(rs10850110), BMP3(rs1390319), GHR (rs2973015,rs6184, rs2973015), FGF7(rs372127537), FGF10(rs593307), and SNAI3(rs4287555) (p < .05) explained most of the variation across the studies, associated to vertical, horizontal, or combined skeletal discrepancies. Meta-analysis results identified a statistically significant association between risk of class III malocclusion of A allele of the FBN3 rs7351083 [OR 2.13; 95% CI 1.1-4.1; p 0.02; recessive model].

CONCLUSION

Skeletal class III is a polygenic trait substantially modulated by ethnicity. A multicentric approach should be considered in future studies to increase sample sizes, applying multivariate analysis such as PCA and cluster analysis to characterize existing sub-phenotypes warranting a deeper analysis of genetic variants contributing to skeletal class III craniofacial disharmony.

CLINICAL RELEVANCE

Grasping the underlying mechanisms of this pathology is critical for a fuller understanding of its etiology, allowing generation of preventive strategies, new individualized therapeutic approaches and more accurate treatment planification strategies.

Molecules, Mechanisms, and Disorders of Self-Domestication: Keys for Understanding Emotional and Social Communication from an Evolutionary Perspective

The neural crest hypothesis states that the phenotypic features of the domestication syndrome are due to a reduced number or disruption of neural crest cells (NCCs) migration, as these cells differentiate at their final destinations and proliferate into different tissues whose activity is reduced by domestication. Comparing the phenotypic characteristics of modern and prehistoric man, it is clear that during their recent evolutionary past, humans also went through a process of self-domestication with a simultaneous prolongation of the period of socialization. This has led to the development of social abilities and skills, especially language, as well as neoteny. Disorders of neural crest cell development and migration lead to many different conditions such as Waardenburg syndrome, Hirschsprung disease, fetal alcohol syndrome, DiGeorge and Treacher-Collins syndrome, for which the mechanisms are already relatively well-known. However, for others, such as Williams-Beuren syndrome and schizophrenia that have the characteristics of hyperdomestication, and autism spectrum disorders, and 7dupASD syndrome that have the characteristics of hypodomestication, much less is known. Thus, deciphering the biological determinants of disordered self-domestication has great potential for elucidating the normal and disturbed ontogenesis of humans, as well as for the understanding of evolution of mammals in general.

Unravelling selection signatures in a single dog breed suggests recent selection for morphological and behavioral traits

Strong selection has resulted in substantial morphological and behavioral diversity across modern dog breeds, which makes dogs interesting model animals to study the underlying genetic architecture of these traits. However, results from between-breed analyses may confound selection signatures for behavior and morphological features that were coselected during breed development. In this study, we assess population genetic differences in a unique resource of dogs of the same breed but with systematic behavioral selection in only one population. We exploit these different breeding backgrounds to identify signatures of recent selection. Selection signatures within populations were found on chromosomes 4 and 19, with the strongest signals in behavior-related genes. Regions showing strong signals of divergent selection were located on chromosomes 1, 24, and 32, and include candidate genes for both physical features and behavior. Some of the selection signatures appear to be driven by loci associated with coat color (Chr 24; ASIP) and length (Chr 32; FGF5), while others showed evidence of association with behavior. Our findings suggest that signatures of selection within dog breeds have been driven by selection for morphology and behavior. Furthermore, we demonstrate that combining selection scans with association analyses is effective for dissecting the traits under selection.

Domesticating olfaction: Dog breeds, including scent hounds, have reduced cribriform plate morphology relative to wolves

The domestic dog is assumed by nearly everyone to be the consummate smeller. Within the species Canis familiaris individual breeds, such as the bloodhound or beagle, are known as olfactory stars. These are "scent breeds," a grouping variably defined as a genetic clade or breed class commonly used for scent detection tasks. Previous work suggests that the dog has a more robust olfactory anatomy than many mammal species. Now we undertake a closer investigation of the dog's olfactory system, both in relationship to its closest wild relatives, the wolf and coyote, and across individual breeds. First, we seek to resolve whether the dog has lost olfactory capacity through its domestication from the wolf lineage. Second, we test the inertial lore that among dogs, "scent breeds," have a superior olfactory facility. As a measure of relative olfactory capacity, we look to the cribriform plate (CP), a bony cup in the posterior nasal cavity perforated by passageways for all olfactory nerve bundles streaming from the periphery to the brain. Using high-resolution computed tomography (CT) scans and digital quantification, we compare relative CP size in 46 dog breeds, the coyote and gray wolf. Results show the dog has a reduced CP surface area relative to the wolf and coyote. Moreover, we found no significant differences between CP size of "scent" and "non-scent" breeds. Our study suggests that the dog lost olfactory capacity as a result of domestication and this loss was not recovered in particular breed groupings through directed artificial selection for increased olfactory facility.

Variants in FtsJ RNA 2'-O-Methyltransferase 3 and Growth Hormone 1 are associated with small body size and a dental anomaly in dogs

Domesticated dogs show unparalleled diversity in body size across breeds, but within breeds variation is limited by selective breeding. Many heritable diseases of dogs are found among breeds of similar sizes, suggesting that as in humans, alleles governing growth have pleiotropic effects. Here, we conducted independent genome-wide association studies in the small Shetland Sheepdog breed and discovered a locus on chromosome 9 that is associated with a dental abnormality called maxillary canine-tooth mesioversion (MCM) (P = 1.53 × 10^-7) as well as two body size traits: height (P = 1.67 × 10^-5) and weight (P = 1.16 × 10^-7). Using whole-genome resequencing data, we identified variants in two proximal genes: FTSJ3, encoding an RNA methyltransferase, and GH1, encoding growth hormone. A substitution in FTSJ3 and a splice donor insertion in GH1 are strongly associated with MCM and reduced body size in Shetland Sheepdogs. We demonstrated in vitro that the GH1 variant leads to exon 3 skipping, predicting a mutant protein known to cause human pituitary dwarfism. Statistical modeling, however, indicates that the FTSJ3 variant is the stronger predictor of MCM and that each derived allele reduces body size by about 1 inch and 5 pounds. In a survey of 224 breeds, both FTSJ3 and GH1 variants are frequent among very small "toy" breeds and absent from larger breeds. Our findings indicate that a chromosome 9 locus harboring tightly linked variants in FTSJ3 and GH1 reduces growth in the Shetland Sheepdog and toy breed dogs and confers risk for MCM through vertical pleiotropy.

Deciphering the puzzles of dog domestication

The domestic dog, as a highly successful domestication model, is well known as a favored human companion. Exploring its domestication history should provide great insight into our understanding of the prehistoric development of human culture and productivity. Furthermore, investigation on the mechanisms underpinning the morphological and behavioral traits associated with canid domestication syndrome is of significance not only for scientific study but also for human medical research. Current development of a multidisciplinary canine genome database, which includes enormous omics data, has substantially improved our understanding of the genetic makeup of dogs. Here, we reviewed recent advances associated with the original history and genetic basis underlying environmental adaptations and phenotypic diversities in domestic dogs, which should provide perspectives on improving the communicative relationship between dogs and humans.

Brain gyrification in wild and domestic canids: Has domestication changed the gyrification index in domestic dogs?

Over the last 15 years, research on canid cognition has revealed that domestic dogs possess a surprising array of complex sociocognitive skills pointing to the possibility that the domestication process might have uniquely altered their brains; however, we know very little about how evolutionary processes (natural or artificial) might have modified underlying neural structure to support species-specific behaviors. Evaluating the degree of cortical folding (i.e., gyrification) within canids may prove useful, as this parameter is linked to functional variation of the cerebral cortex. Using quantitative magnetic resonance imaging to investigate the impact of domestication on the canine cortical surface, we compared the gyrification index (GI) in 19 carnivore species, including six wild canid and 13 domestic dog individuals. We also explored correlations between global and local GI with brain mass, cortical thickness, white and gray matter volume and surface area. Our results indicated that GI values for domestic dogs are largely consistent with what would be expected for a canid of their given brain mass, although more variable than that observed in wild canids. We also found that GI in canids is positively correlated with cortical surface area, cortical thickness and total cortical gray matter volumes. While we found no evidence of global differences in GI between domestic and wild canids, certain regional differences in gyrification were observed.

From head to hind: Elucidating function through contrasting morphometrics of ancient and modern pedigree dogs

Used together, caliper- and geometric-based morphometric analyses provide complimentary approaches to classifying form and function of archaeozoological remains. Here we apply these analytical tools to the skeletal remains of an ancient male dog unearthed from a rural farm settlement of Roman date near present day Warmington, United Kingdom. Our comparisons of the Warmington Roman dog against the morphological characteristics of modern dog breeds enabled us to establish the former's size and shape. It was of medium stature. Analysis of viscerocrania and neurocrania indicate it falls within the meso- to dolichocephalic rankings of modern dogs. The neurocranium shape and the dimensions of its long bones strongly suggest that the Warmington dog shares similarities to modern sight hounds. Historically sight hounds were bred for speed, as necessitated of a hunter that runs down small prey. Our analysis suggests that the Warmington dog was likely bred for, or derived from, Roman hunting stock. By revealing the Warmington Roman dog's form from cranial and postcranial analyses, we shed light on Roman life in one of the furthest outposts of the Roman Empire.

Genetic Variants Affecting Skeletal Morphology in Domestic Dogs

Purebred dog breeds provide a powerful resource for the discovery of genetic variants affecting skeletal morphology. Domesticated and subsequently purebred dogs have undergone strong artificial selection for a broad range of skeletal variation, which include both the size and shapes of their bones. While the phenotypic variation between breeds is high, within-breed morphological variation is typically low. Approaches for defining genetic variants associated with canine morphology include quantitative within-breed analyses, as well as across-breed analyses, using breed standards as proxies for individual measurements. The ability to identify variants across the genomes of individual dogs can now be paired with precise measures of morphological variation to define the genetic interactions and the phenotypic effect of variants on skeletal morphology.

"Forever young"-Postnatal growth inhibition of the turbinal skeleton in brachycephalic dog breeds (Canis lupus familiaris)

In short snouted (brachycephalic) dogs (Canis lupus familiaris), several genetic mutations cause postnatal growth inhibition of the viscerocranium. Thus, for example, the pug keeps a snub nose like that observed in neonate dogs in general. However, little is known how far intranasal structures like the turbinal skeleton are also affected. In the present study, we provide the first detailed morphological and morphometric analyses on the turbinal skeleton of pug, Japanese chin, pekingese, King Charles spaniel, and Cavalier. In order to elucidate how a shortened snout affects turbinal shape, size, and density, our sample covers different degrees of brachycephaly. Macerated skulls of 1 juvenile and 17 adult individuals were investigated by μCT and virtual 3D reconstructions. In addition, histological serial sections of two prenatal and one neonate whippet were taken into account. All investigated postnatal stages show three frontoturbinals and three ethmoturbinals similar to longer snouted breeds, whereas the number of interturbinals is reduced. The shape of the entire turbinal skeleton simplifies with decreasing snout length, that is, within a minimized nasal cavity the turbinals decrease proportionally in surface area and surface density due to a looser arrangement. We interpret these apparent reductions as a result of spatial constraint which affects postnatal appositional bone growth and the position of the turbinals inside the nasal cavity. The turbinal skeleton of brachycephalic dogs arrests at an early ontogenetic stage, corresponding with previous studies on the dermal bones. Hence, we assume an association between the growth of intranasal structures and facial elongation.

Canine Brachycephaly: Anatomy, Pathology, Genetics and Welfare

Brachycephalic dog breeds have experienced a marked rise in popularity in recent years. While numerous people clearly desire this phenotype in their pets, many of these dogs unfortunately experience several concomitant sequelae, including major problems with respiration and thermoregulation, as well as gastrointestinal, ophthalmological, dermatological, reproductive and even dental problems. This mini review focuses on the anatomical and pathological changes associated with brachycephalic skull shape, including brachycephalic obstructive airway syndrome and other co-existent disorders. It then details the known genetic contributors to brachycephaly, and concludes with a brief discourse on the welfare of these animals.

Intrinsic Ball Retrieving in Wolf Puppies Suggests Standing Ancestral Variation for Human-Directed Play Behavior

Domestication dramatically alters phenotypes across animal species. Standing variation among ancestral populations often drives phenotypic change during domestication, but some changes are caused by novel mutations. In dogs (Canis familiaris) specifically, it has been suggested that the ability to interpret social-communicative behavior expressed by humans originated post-domestication and this behavior is thus not expected to occur in wolves (Canis lupus). Here we report the observation of three 8-week-old wolf puppies spontaneously responding to social-communicative behaviors from an unfamiliar person by retrieving a ball. This behavioral expression in wolves has significant implications for our understanding and expectations of the genetic foundations of dog behavior. Importantly, our observations indicate that behavioral responses to human social-communicative cues are not unique to dogs. This suggests that, although probably rare, standing variation in the expression of human-directed behavior in ancestral populations could have been an important target for early selective pressures exerted during dog domestication.

•

Dogs, unlike wolves, are suggested to have the ability to interpret human cues

•

Yet, untrained wolf puppies will spontaneously retrieve a ball based on human cues

•

Standing variation for human-directed behavior in ancestral wolf populations

•

Human-directed behavior important during early selection in dog domestication

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.