Hair is a unique structure, characteristic of mammals, controlling body homeostasis, as well as cell and tissue integration. Previous studies in dog, mouse, and rat have identified polymorphisms in Keratin 71 (KRT71) as responsible for the curly/wavy phenotypes. The coding sequence and the 3′ UTR of KRT71 were directly sequenced in randomly bred and pedigreed domestic cats with different pelage mutations, including hairless varieties. A SNP altering a splice site was identified in the Sphynx breed and suggested to be the hairless (hr) allele, and a complex sequence alteration, also causing a splice variation, was identified in the Devon Rex breed and suggested to be the curly (re) allele. The polymorphisms were genotyped in approximately 200 cats. All the Devon Rex were homozygous for the complex alterations and most of the Sphynx were either homozygous for the hr allele or compound heterozygotes with the Devon-associated re allele, suggesting that the phenotypes are a result of the identified SNPs. Two Sphynx carrying the proposed hr mutation did not carry the Devon-associated alteration. No other causative mutations for eight different rexoid and hairless cat phenotypes were identified. The allelic series KRT71 ⁺ > KRT71 ^hr > KRT71 ^re is suggested.

Domestic cats have a unique breeding history and can be used as models for human hereditary and infectious diseases. In the current era of genome-wide association studies, insights regarding linkage disequilibrium (LD) are essential for efficient association studies. The objective of this study is to investigate the extent of LD in the domestic cat, Felis silvestris catus, particularly within its breeds. A custom illumina GoldenGate Assay consisting of 1536 single nucleotide polymorphisms (SNPs) equally divided over ten 1 Mb chromosomal regions was developed, and genotyped across 18 globally recognized cat breeds and two distinct random bred populations. The pair-wise LD descriptive measure (r²) was calculated between the SNPs in each region and within each population independently. LD decay was estimated by determining the non-linear least-squares of all pair-wise estimates as a function of distance using established models. The point of 50% decay of r² was used to compare the extent of LD between breeds. The longest extent of LD was observed in the Burmese breed, where the distance at which r² ≈ 0.25 was ∼380 kb, comparable to several horse and dog breeds. The shortest extent of LD was found in the Siberian breed, with an r² ≈ 0.25 at approximately 17 kb, comparable to random bred cats and human populations. A comprehensive haplotype analysis was also conducted. The haplotype structure of each region within each breed mirrored the LD estimates. The LD of cat breeds largely reflects the breeds’ population history and breeding strategies. Understanding LD in diverse populations will contribute to an efficient use of the newly developed SNP array for the cat in the design of genome-wide association studies, as well as to the interpretation of results for the fine mapping of disease and phenotypic traits.

The reduced cost and improved efficiency of whole genome sequencing (WGS) is drastically improving the development of cats as biomedical models. Persian cats are models for Leber’s congenital amaurosis (LCA), the most severe and earliest onset form of visual impairment in humans. Cats with innocuous breed-defining traits, such as a bobbed tail, can also be models for somite segmentation and vertebral column development. The first WGS in cats was conducted on a trio segregating for LCA and the bobbed tail abnormality. Variants were identified using FreeBayes and effects predicted using SnpEff. Variants within a known haplotype block for cat LCA and specific candidate genes for both phenotypes were prioritized by the predicted variant effect on the proteins and concordant segregation within the trio. The efficiency of WGS of a single trio of domestic cats was evaluated. A stop gain was identified at position c.577C > T in cat AIPL1, a predicted p.Arg193*. A c.5A > G variant causing a p.V2A was identified in HES7. The variants segregated concordantly in a Persian – Japanese bobtail pedigree. Over 1700 cats from 40 different breeds and populations were genotyped for the AIPL1 variant, defining an allelic frequency in only Persian –related breeds of 1.15 %. A sub-set of cats was genotyped for the HES7 variant, supporting the variant as private to the Japanese bobtail breed. Approximately 18 million SNPs were identified for application in cat research. The cat AIPL1 variant would have been considered a high priority variant for evaluation, regardless of a priori knowledge from previous genetic studies. This study represents the first effort of the 99 Lives Cat Genome Sequencing Initiative to identify disease - causing variants in the domestic cat using WGS. The current cat reference assembly is efficient for gene and variant identification. However, as the feline variant database improves, development of cats as biomedical models for human disease will be more efficient, providing an alternative, large animal model for drug and gene therapy trials. Undiagnosed human patients with early-onset blindness should be screened for this AIPL1 variant. The HES7 variant should further calibrate the somite segmentation clock.

One of the salient features of the domestic cat is the aesthetics of its fur. The Selkirk Rex breed is defined by an autosomal dominant woolly rexoid hair (ADWH) abnormality that is characterized by tightly curled hair shafts. A genome-wide case – control association study was conducted using 9 curly coated Selkirk Rex and 29 controls, including straight-coated Selkirk Rex, British Shorthair and Persian, to localize the Selkirk autosomal dominant rexoid locus (SADRE). Although the control cats were from different breed lineages, they share recent breeding histories and were validated as controls by Bayesian clustering, multi-dimensional scaling and genomic inflation. A significant association was found on cat chromosome B4 (P_raw = 2.87 × 10⁻¹¹), and a unique haplotype spanning ~600 Kb was found in all the curly coated cats. Direct sequencing of four candidate genes revealed a splice site variant within the KRT71 gene associated with the hair abnormality in Selkirk Rex.