Evidence of selection signatures that shape the Persian cat breed

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.

Chromosome-scale assembly with improved annotation provides insights into breed-wide genomic structure and diversity in domestic cats

INTRODUCTION

Comprehensive genomic resources offer insights into biological features, including traits/disease-related genetic loci. The current reference genome assembly for the domestic cat (Felis catus), Felis_Catus_9.0 (felCat9), derived from sequences of the Abyssinian cat, may inadequately represent the general cat population, limiting the extent of deducible genetic variations.

OBJECTIVES

The goal was to develop Anicom American Shorthair 1.0 (AnAms1.0), a reference-grade chromosome-scale cat genome assembly.

METHODS

In contrast to prior assemblies relying on Abyssinian cat sequences, AnAms1.0 was constructed from the sequences of more popular American Shorthair breed, which is related to more breeds than the Abyssinian cat. By combining advanced genomics technologies, including PacBio long-read sequencing and Hi-C- and optical mapping data-based sequence scaffolding, we compared AnAms1.0 to existing Felidae genome assemblies (20 scaffolds, scaffolds N50 > 150 Mbp). Homology-based and ab initio gene annotation through Iso-Seq and RNA-Seq was used to identify new coding genes and splice variants.

RESULTS

AnAms1.0 demonstrated superior contiguity and accuracy than existing Felidae genome assemblies. Using AnAms1.0, we identified over 1.5 thousand structural variants and 29 million repetitions compared to felCat9. Additionally, we identified > 1,600 novel protein-coding genes. Notably, olfactory receptor structural variants and cardiomyopathy-related variants were identified.

CONCLUSION

AnAms1.0 facilitates the discovery of novel genes related to normal and disease phenotypes in domestic cats. The analyzed data are publicly accessible on Cats-I (https://cat.annotation.jp/), which we established as a platform for accumulating and sharing genomic resources to discover novel genetic traits and advance veterinary medicine.

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.

An ancient haplotype containing antimicrobial peptide gene variants is associated with severe fungal skin disease in Persian cats

Dermatophytosis, also known as ringworm, is a contagious fungal skin disease affecting humans and animals worldwide. Persian cats exhibit severe forms of the disease more commonly than other breeds of cat, including other long-haired breeds. Certain types of severe dermatophytosis in humans are reportedly caused by monogenic inborn errors of immunity. The goal of this study was to identify genetic variants in Persian cats contributing to the phenotype of severe dermatophytosis. Whole-genome sequencing of case and control Persian cats followed by a genome-wide association study identified a highly divergent, disease-associated haplotype on chromosome F1 containing the S100 family of genes. S100 calcium binding protein A9 (S100A9), which encodes a subunit of the antimicrobial heterodimer known as calprotectin, contained 13 nonsynonymous variants between cases and controls. Evolutionary analysis of S100A9 haplotypes comparing cases, controls, and wild felids suggested the divergent disease-associated haplotype was likely introgressed into the domestic cat lineage and maintained via balancing selection. We demonstrated marked upregulation of calprotectin expression in the feline epidermis during dermatophytosis, suggesting involvement in disease pathogenesis. Given this divergent allele has been maintained in domestic cat and wildcat populations, this haplotype may have beneficial effects against other pathogens. The pathogen specificity of this altered protein should be investigated before attempting to reduce the allele frequency in the Persian cat breed. Further work is needed to clarify if severe Persian dermatophytosis is a monogenic disease or if hidden disease-susceptibility loci remain to be discovered. Consideration should be given to engineering antimicrobial peptides such as calprotectin for topical treatment of dermatophytosis in humans and animals.

Fungal skin infections known as ringworm or dermatophytosis affect billions of humans and animals worldwide. Normally the disease is self-limiting in affected individuals. The Persian cat breed is a popular breed known for its long hair coat and short nose as well as its propensity to develop severe, chronic dermatophytosis. By examining the genomes of Persian cats, we discovered that a specific region of DNA is highly altered between cats with and without severe dermatophytosis. The DNA sequence in this region is particularly divergent within a cluster of genes involved in immune defense against pathogens. Notably, alterations to the DNA sequence cause several changes in the antimicrobial protein known as calprotectin, which defends against pathogens in the skin of cats. Persian cats with severe dermatophytosis have a version of calprotectin similar to a version maintained by certain desert-dwelling wild felids such as sand cats and Asiatic wildcats. Therefore, we think this version of the protein is beneficial in some environments or against certain pathogens but not against the fungus that causes ringworm in cats. Our findings suggest changes to calprotectin may affect pathogen specificity and engineered calprotectin could be considered as a novel therapy for dermatophytosis in humans and animals.

Facial conformation characteristics in Persian and Exotic Shorthair cats

OBJECTIVES

The primary objectives of the study were to examine the diversity in facial conformation characteristics within a group of Persian (PER) and Exotic Shorthair (EXO) show cats, and to contrast the results to findings within a group of non-purebred domestic shorthair (DSH) and domestic longhair (DLH) cats. The secondary objectives were to determine the PER/EXO show cat owners' perceptions of the breathing status of their cats, and to evaluate if remarks from the cat show judges concerning the cats' head and facial conformation were exclusively related to the aesthetic features of the cats.

METHODS

Sixty-four PER and 12 EXO show cats were prospectively examined at five international cat shows, and 20 DSH/DLH cats were examined at an animal hospital. Facial conformation characteristics were evaluated by examining photos of the cats. Owners of the PER/EXO show cats answered a questionnaire concerning their cats' health status, and they were encouraged to send in the judges' score sheets from the cat shows.

RESULTS

The PER/EXO show cats had higher diversity in facial conformation characteristics than the DSH/DLH cats, and high incidences of hypoplasia of the nose leather (95%), the nose leather top positioned above the level of the lower eyelid (93%), moderate-to-severe stenotic nares (86%), epiphora (83%) and entropion (32%). Owners of 6/76 PER/EXO show cats stated that their cat had increased respiratory sounds and/or trouble breathing at least once a week. The cat show judges' written comments were exclusively related to aesthetic features of the cats' head and facial conformation details.

CONCLUSIONS AND RELEVANCE

Hypoplasia of the nose leather, high position of the nose leather top, stenotic nares, epiphora and entropion were common findings in the PER/EXO show cats but not in the DSH/DLH cats. Few of the cat owners perceived that their cat had problems related to the airways.

More Than a Moggy; A Population Genetics Analysis of the United Kingdom's Non-Pedigree Cats

The domestic cat is one of the most popular pets in the world. It is estimated that 89–92% of domestic cats in the UK are non-pedigree Domestic shorthair (DSH), Domestic longhair (DLH), or Domestic semi-longhair cats (DSLH). Despite their popularity, little is known of the UK non-pedigree cats’ population structure and breeding dynamics. Using a custom designed single nucleotide variant (SNV) array, this study investigated the population genetics of 1344 UK cats. Principal components analysis (PCA) and fastSTRUCTURE analysis verified that the UK’s DSH, DLH, and DSLH cats are random-bred, rather than admixed, mix breed, or crossbred. In contrast to pedigree cats, the linkage disequilibrium of these random-bred cats was least extensive and decayed rapidly. Homozygosity by descent (HBD) analysis showed the majority of non-pedigree cats had proportionally less of their genome in HBD segments compared to pedigree cats, and that these segments were older. Together, these findings suggest that the DSH, DLH, and DSLH cats should be considered as a population of random-bred cats rather than a crossbred or pedigree-admixed cat. Unexpectedly, 19% of random-bred cat genomes displayed a higher proportion of HBD segments associated with more recent inbreeding events. Therefore, while non-pedigree cats as a whole are genetically diverse, they are not impervious to inbreeding and its health risks.

Comparison of intraocular pressure in healthy brachycephalic and nonbrachycephalic cats using the Icare® TONOVET Plus rebound tonometer

OBJECTIVE

To compare intraocular pressure using the Icare® TONOVET Plus rebound tonometer in healthy brachycephalic and nonbrachycephalic cats.

ANIMALS STUDIED

Both eyes of 78 healthy cats were investigated in this study. Cats were divided into two groups: brachycephalic (n = 39) and nonbrachycephalic (n = 39).

PROCEDURES

Nose position and muzzle ratio were photographically recorded and analyzed. Physical and ophthalmic examinations were performed. Intraocular pressure was measured using the Icare® TONOVET Plus rebound tonometry instrument. Quantitative mean values were statistically compared using an unpaired t-test at a significance level of p < .05.

RESULTS

Mean values of the nose position and muzzle ratio were significantly lower in the brachycephalic group (20.14 ± 5.43%, 9.61 ± 3.29%) compared with the nonbrachycephalic group (29.21 ± 4.30%, 13.97 ± 6.01%). The mean intraocular pressure for brachycephalic cats (15.76 ± 0.50 mmHg) was significantly lower (p < .001) than for nonbrachycephalic cats (18.77 ± 0.49 mmHg).

CONCLUSIONS

Intraocular pressure was significantly lower in brachycephalic cats using the Icare® TONOVET Plus rebound tonometer. Intraocular pressure values obtained in this study could be used as a guideline for measurements obtained using this tonometry device in healthy brachycephalic and nonbrachycephalic cats.

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.

Patterns of allele frequency differences among domestic cat breeds assessed by a 63K SNP array

Cats are ubiquitous companion animals that have been keenly associated with humans for thousands of years and only recently have been intentionally bred for aesthetically appealing coat looks and body forms. The intense selection on single gene phenotypes and the various breeding histories of cat breeds have left different marks on the genomes. Using a previously published 63K Feline SNP array dataset of twenty-six cat breeds, this study utilized a genetic differentiation-based method (di) to empirically identify candidate regions under selection. Defined as three or more overlapping (500Kb) windows of high levels of population differentiation, we identified a total of 205 candidate regions under selection across cat breeds with an average of 6 candidate regions per breed and an average size of 1.5 Mb per candidate region. Using the combined size of candidate regions of each breed, we conservatively estimate that a minimum of ~ 0.1-0.7% of the autosomal genome is potentially under selection in cats. As positive controls and tests of our methodology, we explored the candidate regions of known breed-defining genes (e.g., FGF5 for longhaired breeds) and we were able to detect the genes within candidate regions, each in its corresponding breed. For breed specific exploration of candidate regions under selection, eleven representative candidate regions were found to encompass potential candidate genes for several phenotypes such as brachycephaly of Persian (DLX6, DLX5, DLX2), curled ears of American Curl (MCRIP2, PBX1), and body-form of Siamese and Oriental (ADGRD1), which encourages further molecular investigations. The current assessment of the candidate regions under selection is empiric and detailed analyses are needed to rigorously disentangle effects of demography and population structure from artificial selection.

Keep Garfagnina alive. An integrated study on patterns of homozygosity, genomic inbreeding, admixture and breed traceability of the Italian Garfagnina goat breed

The objective of this study was to investigate the genetic diversity of the Garfagnina (GRF) goat, a breed that currently risks extinction. For this purpose, 48 goats were genotyped with the Illumina CaprineSNP50 BeadChip and analyzed together with 214 goats belonging to 9 other Italian breeds (~25 goats/breed), whose genotypes were available from the AdaptMap project [Argentata (ARG), Bionda dell'Adamello (BIO), Ciociara Grigia (CCG), Di Teramo (DIT), Garganica (GAR), Girgentana (GGT), Orobica (ORO), Valdostana (VAL) and Valpassiria (VSS)]. Comparative analyses were conducted on i) runs of homozygosity (ROH), ii) admixture ancestries and iii) the accuracy of breed traceability via discriminant analysis on principal components (DAPC) based on cross-validation. ROH analyses was used to assess the genetic diversity of GRF, while admixture and DAPC to evaluate its relationship to the other breeds. For GRF, common ROH (more than 45% in GRF samples) was detected on CHR 12 at, roughly 50.25-50.94Mbp (ARS1 assembly), which spans the CENPJ (centromere protein) and IL17D (interleukin 17D) genes. The same area of common ROH was also present in DIT, while a broader region (~49.25-51.94Mbp) was shared among the ARG, CCG, and GGT. Admixture analysis revealed a small region of common ancestry from GRF shared by BIO, VSS, ARG and CCG breeds. The DAPC model yielded 100% assignment success for GRF. Overall, our results support the identification of GRF as a distinct native Italian goat breed. This work can contribute to planning conservation programmes to save GRF from extinction and will improve the understanding of the socio-agro-economic factors related with the farming of GRF.

A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species

BACKGROUND

Runs of homozygosity (ROH) have become the state-of-the-art method for analysis of inbreeding in animal populations. Moreover, ROH are suited to detect signatures of selection via ROH islands and are used in other applications, such as genomic prediction and genome-wide association studies (GWAS). Currently, a vast amount of single nucleotide polymorphism (SNP) data is available online, but most of these data have never been used for ROH analysis. Therefore, we performed a ROH analysis on large medium-density SNP datasets in eight animal species (cat, cattle, dog, goat, horse, pig, sheep and water buffalo; 442 different populations) and make these results publicly available.

RESULTS

The results include an overview of ROH islands per population and a comparison of the incidence of these ROH islands among populations from the same species, which can assist researchers when studying other (livestock) populations or when looking for similar signatures of selection. We were able to confirm many known ROH islands, for example signatures of selection for the myostatin (MSTN) gene in sheep and horses. However, our results also included multiple other ROH islands, which are common to many populations and not identified to date (e.g. on chromosomes D4 and E2 in cats and on chromosome 6 in sheep).

CONCLUSIONS

We are confident that our repository of ROH islands is a valuable reference for future studies. The discovered ROH island regions represent a unique starting point for new studies or can be used as a reference for future studies. Furthermore, we encourage authors to add their population-specific ROH findings to our repository.

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.

Genomic Analysis Suggests KITLG is Responsible for a Roan Pattern in two Pakistani Goat Breeds

The roan coat color pattern is described as the presence of white hairs intermixed with pigmented hairs. This kind of pigmentation pattern has been observed in many domestic species, including the goat. The molecular mechanisms and inheritance that underlie this pattern are known for some species and the KITLG gene has been shown associated with this phenotype. To date, no research effort has been carried out to find the gene(s) that control(s) roan coat color pattern in goats. In the present study, after genotyping with the GoatSNP50 BeadChip, 35 goats that showed a roan pattern and that belonged to two Pakistan breeds (Group A) were analyzed and then compared to 740 goats of 39 Italian and Pakistan goat breeds that did not have the same coat color pattern (Group B). Runs of homozygosity-based and XP-EHH analyses were used to identify unique genomic regions potentially associated with the roan pattern. A total of 3 regions on chromosomes 5, 6, and 12 were considered unique among the group A versus group B comparisons. The A region > 1.7 Mb on chromosome 5 was the most divergent between the two groups. This region contains six genes, including the KITLG gene. Our findings support the hypothesis that the KITLG gene may be associated with the roan phenotype in goats.

Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity

Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations.

Persian cats under first opinion veterinary care in the UK: demography, mortality and disorders

Persian cats are a popular cat breed worldwide, and especially in the US, Europe and Asia. This study aimed to describe the demography, common disorders and mortality in Persians under general practice veterinary care in 2013 in the UK. The study population of 285,547 cats overall included 3235 (1.1%) Persians. Mean adult Persian bodyweight was 3.9 kg (SD 0.9) and median age was 7.0 years (IQR 3.3–11.6). At least one disorder was recorded in 2099 (64.9%) Persians. The most common specific disorders were haircoat disorders (411, 12.7%), periodontal disease (365, 11.3%), overgrown nails (234, 7.2%), and ocular discharge (188, 5.8%). The most common disorder groups were dermatological (578, 17.9%), ophthalmological (496, 15.3%) and dental (397, 12.3%). Median longevity was 13.5 years (IQR 9.9–16.0). The most common grouped causes of death were renal disease (102, 23.4%), neoplasia (37, 8.5%) and mass-associated disorder (35, 8.0%). This is the first study to use general practice data to examine the overall health of Persian cats. With haircoat, ocular and dental disorders being the predominant disorders identified, this study highlights the need for increased owner awareness to manage and prevent the typical health problems associated with this breed’s phenotype.

Genome-wide patterns of homozygosity provide clues about the population history and adaptation of goats

Background

Patterns of homozygosity can be influenced by several factors, such as demography, recombination, and selection. Using the goat SNP50 BeadChip, we genotyped 3171 goats belonging to 117 populations with a worldwide distribution. Our objectives were to characterize the number and length of runs of homozygosity (ROH) and to detect ROH hotspots in order to gain new insights into the consequences of neutral and selection processes on the genome-wide homozygosity patterns of goats.

Results

The proportion of the goat genome covered by ROH is, in general, less than 15% with an inverse relationship between ROH length and frequency i.e. short ROH (< 3 Mb) are the most frequent ones. Our data also indicate that ~ 60% of the breeds display low F_ROH coefficients (< 0.10), while ~ 30 and ~ 10% of the goat populations show moderate (0.10 < F_ROH < 0.20) or high (> 0.20) F_ROH values. For populations from Asia, the average number of ROH is smaller and their coverage is lower in goats from the Near East than in goats from Central Asia, which is consistent with the role of the Fertile Crescent as the primary centre of goat domestication. We also observed that local breeds with small population sizes tend to have a larger fraction of the genome covered by ROH compared to breeds with tens or hundreds of thousands of individuals. Five regions on three goat chromosomes i.e. 11, 12 and 18, contain ROH hotspots that overlap with signatures of selection.

Conclusions

Patterns of homozygosity (average number of ROH of 77 and genome coverage of 248 Mb; F_ROH < 0.15) are similar in goats from different geographic areas. The increased homozygosity in local breeds is the consequence of their small population size and geographic isolation as well as of founder effects and recent inbreeding. The existence of three ROH hotspots that co-localize with signatures of selection demonstrates that selection has also played an important role in increasing the homozygosity of specific regions in the goat genome. Finally, most of the goat breeds analysed in this work display low levels of homozygosity, which is favourable for their genetic management and viability.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0424-8) contains supplementary material, which is available to authorized users.

Genome-Wide Signatures of Selection Reveal Genes Associated With Performance in American Quarter Horse Subpopulations

Selective breeding for athletic performance in various disciplines has resulted in population stratification within the American Quarter Horse (QH) breed. The goals of this study were to utilize high density genotype data to: (1) identify genomic regions undergoing positive selection within and among QH subpopulations; (2) investigate haplotype structure within each QH subpopulation; and (3) identify candidate genes within genomic regions of interest (ROI), as well as biological pathways, predicted to play a role in elite performance in each group. For that, 65K SNP genotyping data on 143 elite individuals from 6 QH subpopulations (cutting, halter, racing, reining, western pleasure, and working cow) were imputed to 2M SNPs. Signatures of selection were identified using F_ST-based (d_i ) and haplotype-based (hapFLK) analyses, accompanied by identification of local haplotype structure and sharing within subpopulations (hapQTL). Regions undergoing positive selection were identified on all 31 autosomes, and ROI on 2 chromosomes were identified by all 3 methods combined. Genes within each ROI were retrieved and used to identify pathways and genes that might contribute to performance in each subpopulation. These included, among others, candidate genes associated with skeletal muscle development, metabolism, and central nervous system development. This work improves our understanding of equine breed development, and provides breeders with a better understanding of how selective breeding impacts the performance of QH populations.

Creating diversity in mammalian facial morphology: a review of potential developmental mechanisms

Mammals (class Mammalia) have evolved diverse craniofacial morphology to adapt to a wide range of ecological niches. However, the genetic and developmental mechanisms underlying the diversification of mammalian craniofacial morphology remain largely unknown. In this paper, we focus on the facial length and orofacial clefts of mammals and deduce potential mechanisms that produced diversity in mammalian facial morphology. Small-scale changes in facial morphology from the common ancestor, such as slight changes in facial length and the evolution of the midline cleft in some lineages of bats, could be attributed to heterochrony in facial bone ossification. In contrast, large-scale changes of facial morphology from the common ancestor, such as a truncated, widened face as well as the evolution of the bilateral cleft possessed by some bat species, could be brought about by changes in growth and patterning of the facial primordium (the facial processes) at the early stages of embryogenesis.

Applications and efficiencies of the first cat 63K DNA array

The development of high throughput SNP genotyping technologies has improved the genetic dissection of simple and complex traits in many species including cats. The properties of feline 62,897 SNPs Illumina Infinium iSelect DNA array are described using a dataset of over 2,000 feline samples, the most extensive to date, representing 41 cat breeds, a random bred population, and four wild felid species. Accuracy and efficiency of the array’s genotypes and its utility in performing population-based analyses were evaluated. Average marker distance across the array was 37,741 Kb, and across the dataset, only 1% (625) of the markers exhibited poor genotyping and only 0.35% (221) showed Mendelian errors. Marker polymorphism varied across cat breeds and the average minor allele frequency (MAF) of all markers across domestic cats was 0.21. Population structure analysis confirmed a Western to Eastern structural continuum of cat breeds. Genome-wide linkage disequilibrium ranged from 50–1,500 Kb for domestic cats and 750 Kb for European wildcats (Felis silvestris silvestris). Array use in trait association mapping was investigated under different modes of inheritance, selection and population sizes. The efficient array design and cat genotype dataset continues to advance the understanding of cat breeds and will support monogenic health studies across feline breeds and populations.

Oral, Maxillofacial and Dental Diseases in Captive Cheetahs (Acinonyx jubatus)

Descriptions of several oral, maxillofacial and dental conditions/diseases exist for a variety of captive large felids, but little is reported on the pathology of free roaming large felids. Apart from focal palatine erosions (FPEs) as initially described by Fitch and Fagan (1982) and some reference to absent incisor teeth, few data exist on diseases affecting the oral, maxillofacial and dental structures of cheetahs (Acinonyx jubatus), regardless of their captivity status. This study reports 18 different conditions affecting the teeth, bone and oral cavity soft tissue of cheetahs, based on initial assessment of 256 animals over 11 years (2002-2012) in South Africa and Namibia. This report excludes oral tumours or FPEs, but includes several acquired and developmental conditions never described before.

Weighted likelihood inference of genomic autozygosity patterns in dense genotype data

Background

Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns.

Methods

We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data.

Results

Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies.

Conclusions

This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4312-3) contains supplementary material, which is available to authorized users.

The Relationship between Brachycephalic Head Features in Modern Persian Cats and Dysmorphologies of the Skull and Internal Hydrocephalus

Background

Cat breeders observed a frequent occurrence of internal hydrocephalus in Persian cats with extreme brachycephalic head morphology.

Objective

To investigate a possible relationship among the grade of brachycephaly, ventricular dilatation, and skull dysmorphologies in Persian cats.

Animals

92 Persian‐, 10 Domestic shorthair cats.

Methods

The grade of brachycephaly was determined on skull models based on CT datasets. Cranial measurements were examined with regard to a possible correlation with relative ventricular volume, and cranial capacity. Persians with high (peke‐face Persians) and lower grades of brachycephaly (doll‐face Persians) were investigated for the presence of skull dysmorphologies.

Results

The mean cranial index of the peke‐face Persians (0.97 ± 0.14) was significantly higher than the mean cranial index of doll‐face Persians (0.66 ± 0.04; P < 0.001). Peke‐face Persians had a lower relative nasal bone length (0.15 ± 0.04) compared to doll‐face (0.29 ± 0.08; P < 0.001). The endocranial volume was significantly lower in doll‐face than peke‐face Persians (89.6 ± 1.27% versus 91.76 ± 2.07%; P < 0.001). The cranial index was significantly correlated with this variable (Spearman's r: 0.7; P < 0.0001). Mean ventricle: Brain ratio of the peke‐face group (0.159 ± 0.14) was significantly higher compared to doll‐face Persians (0.015 ± 0.01; P < 0.001).

Conclusion and Clinical Relevance

High grades of brachycephaly are also associated with malformations of the calvarial and facial bones as well as dental malformations. As these dysmorphologies can affect animal welfare, the selection for extreme forms of brachycephaly in Persian cats should be reconsidered.

The Valdostana goat: a genome-wide investigation of the distinctiveness of its selective sweep regions

The Valdostana goat is an alpine breed, raised only in the northern Italian region of the Aosta Valley. This breed's main purpose is to produce milk and meat, but is peculiar for its involvement in the "Batailles de Chèvres," a recent tradition of non-cruel fight tournaments. At both the genetic and genomic levels, only a very limited number of studies have been performed with this breed and there are no studies about the genomic signatures left by selection. In this work, 24 unrelated Valdostana animals were screened for runs of homozygosity to identify highly homozygous regions. Then, six different approaches (ROH comparison, Fst single SNPs and windows based, Bayesian, Rsb, and XP-EHH) were applied comparing the Valdostana dataset with 14 other Italian goat breeds to confirm regions that were different among the comparisons. A total of three regions of selection that were also unique among the Valdostana were identified and located on chromosomes 1, 7, and 12 and contained 144 genes. Enrichment analyses detected genes such as cytokines and lymphocyte/leukocyte proliferation genes involved in the regulation of the immune system. A genetic link between an aggressive challenge, cytokines, and immunity has been hypothesized in many studies both in humans and in other species. Possible hypotheses associated with the signals of selection detected could be therefore related to immune-related factors as well as with the peculiar battle competition, or other breed-specific traits, and provided insights for further investigation of these unique regions, for the understanding and safeguard of the Valdostana breed.