Early-Onset Progressive Retinal Atrophy Associated with an IQCB1 Variant in African Black-Footed Cats (Felis nigripes)

African black-footed cats (Felis nigripes) are endangered wild felids. One male and full-sibling female African black-footed cat developed vision deficits and mydriasis as early as 3 months of age. The diagnosis of early-onset progressive retinal atrophy (PRA) was supported by reduced direct and consensual pupillary light reflexes, phenotypic presence of retinal degeneration, and a non-recordable electroretinogram with negligible amplitudes in both eyes. Whole genome sequencing, conducted on two unaffected parents and one affected offspring was compared to a variant database from 51 domestic cats and a Pallas cat, revealed 50 candidate variants that segregated concordantly with the PRA phenotype. Testing in additional affected cats confirmed that cats homozygous for a 2 base pair (bp) deletion within IQ calmodulin-binding motif-containing protein-1 (IQCB1), the gene that encodes for nephrocystin-5 (NPHP5), had vision loss. The variant segregated concordantly in other related individuals within the pedigree supporting the identification of a recessively inherited early-onset feline PRA. Analysis of the black-footed cat studbook suggests additional captive cats are at risk. Genetic testing for IQCB1 and avoidance of matings between carriers should be added to the species survival plan for captive management.

genetic tests for diseases recognized and observed in zoo populations may prove to be critical for the long-term survival of species that are threatened with extinction and are managed in captive breeding programs.
The black-footed cat (Felis nigripes; Fnig; Burchell, 1824) ( Fig. 1a) is the smallest wild cat found in Africa, endemic to the arid steppe and savannah habitats of the southern African subregion 15 . Phylogenetic studies indicate that the black-footed cat is part of the domestic cat lineage along with several other wild felids, including the European wildcat (Felis silvestris silvestris), African wildcat (Felis silvestris lybica), Chinese desert cat (Felis bieti), Sand cat (Felis margarita), and Jungle cat (Felis chaus) (Fig. 1b) [16][17][18] . Unfortunately, the size of both the wild and captive populations (~50 individuals in zoos worldwide) are diminishing. The species is listed under Appendix I-threatened with extinction-by the Convention on the International Trade of Endangered Species (CITES) 19 , as "vulnerable" by the International Union for Conservation of Nature (IUCN) since 2002 20 , and as "Endangered" by the United States Fish and Wildlife Service 21 . Black-footed cats typically have shorter life spans in captivity than in the wild. The majority of captive adult black-footed cat deaths occur between 2 to 5 years, whereas free-ranging cats have been estimated to live approximately 5 to 6 years 22 . Recent reports have documented that renal amyloidosis plays an important role in the decreased survival of captive black-footed cats 23 . Information concerning other disease processes or conditions that occur in wild or captive black-footed cats is limited. Due to their endangered status and their relatively small captive breeding population, the presence of genetic defects, such as heritable blindness, could be detrimental for the long-term sustainability of captive black-footed cat populations.
Progressive retinal atrophy (PRA) refers to inherited retinal dystrophic or degenerative disorders present in animals. The age of onset and rate of progression of initial photoreceptor death varies with each species and each type of PRA. Initially, the disease usually affects the photoreceptor layer (specifically the rods) and proceeds to affect all retinal layers over time leading to diffuse retinal degeneration and ultimately blindness 24 . Four different types of PRA have been described in the domestic cat (Felis catus) [25][26][27][28] . In Abyssinians, an autosomal dominant early-onset cone-rod dystrophy (rdy) manifests from a mutation in CRX 27,[29][30][31][32][33] , while autosomal recessive late-onset rod-cone degeneration (rdAc) results from a mutation in CEP290 26,[34][35][36][37][38] . In Persians, an autosomal recessive early-onset rod-cone dystrophy has been described and a causal variant in AIPL1 detected 28,39,40 . In Bengal cats, an autosomal recessive PRA has been characterized 25 . Due to their similarity to human globe size, relatively accelerated lifespan and different represented modes of inheritance for PRA, domestic cats represent a useful animal model for retinal diseases in humans such as retinitis pigmentosa and Leber's congenital amaurosis 41 .
The present study describes an early-onset generalized PRA diagnosed in two full-sibling African black-footed cats (Felis nigripes). Clinical abnormalities in affected cats were present as early as 3 months of age. This is the first report of ocular disease in this species. A novel cross-species approach was used for the genetic analyses to rapidly identify a highly associated causal variant. Variants within the 99 Lives cat genome sequence database, which contains whole genome sequences of 51 domestic cats at the time of the study, were compared to the variants identified in the black-footed cats. Variants were prioritized as specific to the Black-footed cat and segregating within the trio with the PRA phenotype. A DNA-based test for management of the disease within the captive breeding program can now be implemented.

Results
Ophthalmic examination and diagnostic results. No significant abnormal findings were observed on physical assessment of the male proband cat. However, both pupils were maximally dilated in well-lit conditions, circular in shape, and unresponsive to light (Fig. 2a). Menace responses and dazzle reflexes were absent, and significantly reduced direct and indirect PLRs were observed bilaterally. The cat did not show appropriate tracking responses. Palpebral reflexes, corneal sensitivity and oculocephalic reflexes were normal. On slit-lamp biomicroscopy, the eyelids, conjunctiva, and anterior chamber were normal. Other abnormalities noted included faint linear superficial dendritic corneal facets and a relucency of the embryonic lens nucleus in the left eye. Indirect ophthalmoscopy identified generalized tapetal hyperreflectivity and severe vascular attenuation bilaterally (Fig. 2b) consistent with end-stage PRA. Schirmer I tear test values (mean ± standard deviation) were 5 mm/min in both eyes, which are decreased compared to the normal reference range for domestic felids (14.3 ± 5.73 mm/min) 42,43 . Bilaterally decreased Schimer I tear test values was most likely secondary to sedation 44,45 . Intraocular pressure (IOP) data (mean ± standard deviation) collected with rebound tonometry revealed 17 mmHg bilaterally, which was within the normal reference interval for domestic felids (22.6 ± 4.0 mmHg) 46 . The cornea was fluorescein negative in both eyes. Physical, ophthalmic, and neuro-ophthalmic examination of the female sibling was similar to the male proband. Schirmer I tear test values were decreased in the right (3 mm/min) and left (5 mm/min) eye, presumably due to the influence of general anesthesia. The IOP measured 21 mmHg and 24 mmHg in the right and left eye, respectively.
In both the male proband and its female full-sibling, the electroretinogram (ERG) was non-recordable with a "flat line" in both eyes, consistent with end-stage PRA. Specific photoreceptor function could not be evaluated due to the flat (non-recordable) waveform in both cats. Complete blood hematology, chemistry, and urinalysis were unremarkable, and plasma taurine levels were within normal reference intervals in both cats. Indirect blood pressure measurements obtained under general anesthesia in both cats were found to be within normal reference ranges for domestic felids 47-49 . CEP290, CRX, and AIPL1 genetic testing. Genetic testing of the black-footed female proband and its two unaffected parents was performed at the commercial service at the Veterinary Genetics Laboratory (VGL) at the University of California-Davis. The cats were confirmed to be wild type for the variants known to cause PRA in Abyssinian type domestic cats, specifically the variants within CRX 27 Fig. 3). Average WGS coverage for each individual was 30X. The available database used for variant filtering included variants from 52 individuals, 51 domestic cats and one opportunistic Pallas cat, each with an average 30X genome coverage. Over 22 million variants were identified when the WGS of the trio of cats was aligned to the available annotated cat genome reference (Felis_catus 6.2 genome assembly) 50 . After excluding variants in common between the three black-footed cats and one Pallas cat (Otocolobus manul), over 13 million were identified. Thus, the one Pallas cat eliminated 41% of identified variants that were common between the two wild felid species and found in domestic cats. Over 234,000 variants segregated concordantly when considering the phenotype within the trio and uniqueness to black-footed cats. Only fifty variants had a predicted high impact on a protein product, including, eight predicted stop gains, 20 predicted splice donor-acceptors and predicted 22 frameshift variants. The prioritization of the variants by filtering and impact is presented in Table 2.
A list of over 240 genes associated with blindness is available on the RetNet database (www.RetNet.org). Of 50 genes with high impact variants, only three genes (RASGRF1, IQCB1 and CCDC114) had been previously associated with blindness (Supplementary Table 1). The RASGFR1 and CCDC114 variants were associated with several complex clinical features including retinal atrophy, whereas the phenotype associated with the IQCB1 variant was always PRA with extrarenal manifestations in 10-15% of cases. The high impact variant detected within the coding region of IQCB1 is a 2 base pair (bp) deletion in exon 13 (c.1282delCT) and was considered the highest priority. The identified 2 bp variant was homozygous for the deletion in the affected cat and heterozygous in both phenotypically normal parental obligate carriers ( Table 3). The polymorphism is predicted to cause a protein truncation at position 428 of the IQCB1 amino acid chain (p.L428*), truncating the last 170 amino acids at the C-terminal of nephrocystin-5 (NPHP5). A pedigree for all living cats (< 50) in the North American black-footed cat studbook, cats with available DNA samples, and related individuals is shown (Fig. 3). The identified c.1282delCT was screened via direct sequencing in 16 available samples; blind individuals (n = 2) tested homozygous for the mutated allele while sighted individuals (n = 13) were heterozygous (n = 5) or homozygous (n = 8) wild type. One sample failed to amplify. In addition, seven moderate impact missense variants were identified in RetNet genes that were unique to black-footed cats and segregated concordantly with the disease phenotypes of the WGS trio ( seven genes were not on the same chromosome as the putative causal IQCB1. The polymorphisms within CHD23, CDH3, CNGB1, DTHD1 (2 variants), RBP3 and USH1C were investigated in the additional black-footed cats using mass spectrometry ( Table 3). Five of seven variants (CHD3, CNGB1, DTHD1_942, RBP3, and USH1C) were excluded because the cats' genotype was discordant with PRA phenotype, assuming a recessive mode of inheritance model as suggested by pedigree analysis. Specifically, these variants were either not homozygous in the second PRA-affected cat or homozygous in other sighted black-footed cats. The remaining two missense variants (CDH23 and DTHD1_144) could not be excluded based on the genotype and phenotype discordance but were heterozygous each in six sighted cats from the pedigree. These variants were considered less likely associated with the phenotype because the variant alleles in CDH23 and DTHD1_144 were detected in unrelated black-footed cats that have not been interbred with the lineage that segregates for vision impairment (Fig. 3). For example, a visually sighted cat (223) is a DTHD1_144-carrier and this cat has no relationship to the affected lineage. WGS also confirmed that the cats of the trio are wild type for the known mutations within AIPL1, CRX and CEP290, confirming the commercial service genetic testing.   Table 3. Comparison of genotype versus phenotype between seven moderate impact missense variants. Gender: male (M), female (F); genotype: deleted 2 base pairs cytosine and thymine (delCT), wild type (wt), IQCB1-mutant (delCT/delCT), IQCB1-carrier (wt/delCT), wild type (wt/wt); nucleotides: guanine (G), adenine (A), thymine (T), cytosine (C); "/" no data due to sample error; a whole genome sequenced trio of cats; b missense variants excluded based on genotype versus phenotype mismatch; c missense variants excluded when compared to the pedigree (Fig. 3). The moderate impact variant alleles are detected in unrelated black-footed cats that never interbred with the line of cats with vision loss.

Discussion
The genetics of the patient is a critical aspect of an individual's health care for some diseases. The reduced cost of WGS and WES has made feasible the examination of entire genomes for individualized healthcare. Animals now have the same benefit since lower sequencing costs have allowed the development of DNA variant databases from normal control animals for comparison to individuals with diseases. The study of rare, simply inherited recessive diseases particularly benefit from the WGS approach. Although a large database of normal DNA variants within the genome of black-footed cats is not available, because these wild felids are phylogenetically closely related to domestic cats, the domestic cat genome variant database was used for comparison of polymorphisms to exclude variants in common between the species. Herein is the first cross-species comparison of WGS data to support the captive management of an endangered species. Vision loss is a significant concern in the vulnerable African black-footed cat as poor eyesight could result in a poor quality of life and failure to thrive for individuals, impacting healthy captive breeding practices for the species as a whole. Clinical examination of the affected male proband and female sibling cat revealed findings consistent with PRA. The ERG confirmed that both cats retained minimal or no retinal function, consistent with end-stage PRA. Clinical abnormalities in affected cats were present as early as 3 months of age and rapidly progressed to end-stage PRA, supporting an early-onset retinal disorder affecting both rod and cone photoreceptors. The male proband and female cat are full-siblings, but from different litters. Pedigree analysis supported an autosomal recessive mode of inheritance since visually normal parents produced both affected (n = 2) and unaffected kittens (n = 3), the parents of affected kittens had consanguinity and the affected kittens were either male (n = 3) and female (n = 2) (Fig. 3).
Results of ophthalmic examination and testing, pedigree analysis, and lack of evidence for nutritional, toxic, and hypertensive etiology for the retinal degeneration in the affected African black-footed cats strongly suggested a genetic basis for early-onset PRA. The male proband and female sibling cat were identified at 7 and 3 months of age, respectively, however initial signs of compromised photoreceptor integrity and function may have occurred at an earlier time point. The phenotype and electrodiagnostic findings in PRA-affected black-footed cats are chronologically most consistent with early-onset PRAs in Persians and Bengals, and inconsistent with autosomal dominant early-onset PRA of CRX (rdy) and autosomal recessive late-onset PRA of CEP290 (rdAc)-mutant Abyssinians. Although the presence of a domestic cat variant for vision loss would be highly unlikely in these black-footed cats, genetic testing for cat PRAs is rapid and low cost. Genotyping excluded the Persian autosomal recessive early-onset AIPL1 variant as well as Abyssinian rdAc and rdy PRA variants. The WGS would allow more exhaustive consideration of these same genes for a species-specific variant.
To identify the causal gene, WGS of a trio of cats, comprising two normal obligate carrier parents, and one affected offspring, was conducted. This study of the single black-footed cat was assisted by comparison variants from the closely related domestic cat database, which included variants from 51 domestic cats and one Pallas cat (Otocolobus manul) at the time of the study. The Pallas cat was an opportunistic inclusion in the 99 Lives dataset, included as part of a separate study focusing on polycystic kidney disease. The Pallas cat is the next nearest feline lineage and did support the filtering of shared wild felid variants and shared variants with domestics that would not be causal. Forty-one percent of variants were eliminated by the cross-species comparison. One of the fifty high impact variants resided in a gene known to cause retinal degeneration (IQCB1) and listed on RetNet. A 2 bp deletion within the coding region of the black-footed cats IQCB1 was detected and segregated concordantly with the disease phenotype (Fig. 3). The mutation is predicted to cause a protein truncation at position 428 of the amino acid chain (p.L428*), deleting an IQ calmodulin-binding protein domain and the complete protein C-terminal 51 . Seven moderate impact missense variants were also detected in RetNet genes that were unique to the three sequenced black-footed cats and segregated with the disease. Mismatch genotype and phenotype in addition to the probability of an unrelated cat being affected excluded five of the missense variants as potential candidate genes. Additionally, the remaining two variants were most likely not associated with the retinal degeneration since they were present in individuals unrelated to the affected cats. Although more likely not associated, the variants within CDH23 and DTHD1_144 will be tested along with the IQCB1 2 bp deletion in all samples that will be received in the future, until when excluded or confirmed by the phenotype.
Variants in IQCB1 (IQ calmodulin-binding motif-containing protein-1) are known to historically cause Senior-Loken syndrome (SLSN) in humans, a rare autosomal recessive disease characterized by both nephronophthisis (NPHP) and retinal degeneration 51,52 . All affected human patients with SLSN exhibit Leber congenital amauorsis (LCA). Leber congenital amauorsis is generally inherited in an autosomal recessive manner and associated with over 21 gene variants in animals and humans related to retinal function, including RPE65, CRX, CEP290, AIPL1, and IQCB1 53 . The IQCB1 gene encodes nephrocystin-5 (NPHP5), one of thirteen ciliary proteins collectively known as nephronophthisis (NPHP) proteins 54 . The NPHP5 protein is present in the photoreceptor connecting cilium, a zone that connects the photoreceptor cell body to the outer segments, where IQCB1 Interacts with CEP290 51,55 . The ciliary protein contains two calmodulin binding sites and co-localizes with nephrocystin-1, nephrocystin-4 51 . The NPHP5 protein also interacts with nephrocystin-6 55 and retinal ciliopathy gene retinitis pigmentosa GTPase regulator (RPGR) 51 . Premature terminating mutations in NPHP5 lead to a truncated protein, which disrupts peptide transportation, compromising the structural and functional integrity of the photoreceptor, leading to cone and rod death and subsequent retinal degeneration 51,55 . Variants in NPHP proteins also cause a broad-spectrum of autosomal recessive kidney diseases that are the most frequent genetic cause of end-stage kidney disease in children and young adults 56 . Extra-renal manifestations occur in 10-15% of cases and retinal degeneration are the most common abnormality. Thus, kidney function of PRA-affected black-footed cats should be monitored.
IQCB1 variants have not been recognized in the domestic cat but an insertional mutation in NPHP5 (IQCB1) has been recently documented in American Pit Bull Terrier dogs with early-onset recessively inherited cone-rod dystrophy (crd2) 57,58 . In mutant dogs, the resulting ciliopathy causes early loss of rod photoreceptors and relative retention of central retinal cone photoreceptors that lack function 58 . In the affected dogs, non-ophthalmic problems were neither observed nor reported. Similarly, a variant in NPHP4 is associated with recessive cone-rod dystrophy in Standard Wire-haired Dachshunds. No kidney involvement has been reported in affected dachshunds 59 .
The pedigree reveals three unique lineages in the African black-footed cat captive population (Fig. 3). The IQCB1-mutant gene was isolated within one lineage, but has crossed to another lineage with a recent breeding of cats. Ideally, the IQCB1 variant should be eliminated within the breeding population of captive black-footed cats. This would require all living cats (< 50) to be genotyped. Specifically cats with an increased risk for affected or carrier status based on the pedigree analysis and known genotypes ( Fig. 3; 209, 260, 261, 262, 264, 265, and  266). IQCB1-carriers (or CDH23 and DTHD1 carriers) would need to be mated with wild type cats. If significant population changes occurred in African black-footed cats that necessitated breeding of IQCB1-mutant cats from a conservation standpoint, they should be bred to wild type individuals to maintain maximum genetic diversity within this species. Additionally, the black-footed cats related to the founder cats of the PRA lineage could be genotyped to prevent further influx of the variant into captive breeding programs worldwide. This is the first report of early-onset recessively inherited PRA in the African black-footed cat. Genetic analysis revealed the etiology to likely be a frameshift mutation in the IQCB1 gene, or a missense mutation within CDH23 or DTHD1. IQCB1 represents the most likely variant associated with the phenotype, and the identified polymorphism is predicted to code for a truncated NPHP5 protein in the connecting cilium which is present in both rod and cone photoreceptors. Genetic testing all the candidate variants in more affected and sighted individuals will help confirm the causal variant associated with the PRA. The identification of the putative causal gene provides a DNA-based test for appropriate captive breeding practices and conservation of this vulnerable wild felid. Moreover, the WGS variant database comprised of 51 domestic cats and one Pallas cat at the time of the study has proven useful for cross-species comparisons for species that have recent common ancestors.

Methods
Ethical statement. All procedures used to evaluate the animals' health were performed during the animals' routine health screening at each participating zoo. Direct ophthalmoscopy and phlebotomy via jugular venipuncture were performed as part of the complete physical exam on each individual. The main purpose of the blood specimen collection was to evaluate overall systemic health, and surplus EDTA anti-coagulated whole blood was submitted for genetic analyses. This study was completed in strict accordance with the recommendations from the Kansas City Zoo and Cleveland Metroparks Zoo Scientific Committees. Experimental protocols and applications for use of biomaterials were approved by the Kansas City, Cleveland Metroparks and Omaha's Henry Doorly Zoo Research Review Committees. Examinations, blood or hair specimen collection, and advanced diagnostics were performed under sedation or general anesthesia. Samples were received for the genetic analyses under the University of Missouri Institutional Animal Care and Use Committee protocol 8313. Blood or hair samples submitted by other participating zoological parks were collected opportunistically during the cat's animal health exam and were submitted for diagnostic purposes.
Animals and Clinical descriptions. The 7 month-old neutered male proband African black-footed cat from the Kansas City Zoo presented to the University of Missouri, College of Veterinary Medicine (MU CVM), Columbia, MO for suspected blindness, dilated pupils, absent dazzle reflex in both eyes, and inappropriate tracking of moving objects during quarantine examination. Mild horizontal nystagmus was also present in the right eye. Before presentation, the cat had been screened and treated for parasites, vaccinated for feline viral rhinotracheitis, calicivirus, and panleukopenia (FVRCP) (Boehringer Ingelheim, St. Joseph, MO) in addition to Purevax Ferret Distemper (Merial Inc., Athens, GA) and was considered otherwise healthy. The cat was not receiving any medications and had never been administered fluoroquinolones. Daily diet provided once a day included Evo turkey and chicken formula (15 grams), Nebraska Canine Premium (40 grams) and a small mouse.
From a different litter, the 3 month-old intact female sibling African black-footed cat from the Cleveland Metroparks Zoo presented to the Metropolitan Veterinary Hospital, Akron, OH for fixed and dilated pupils and vision deficits bilaterally. The female sibling had a history of poor weight gain that resolved with treatment and increased caloric intake. On presentation, the cat had received its first vaccination for feline viral rhinotracheitis, calicivirus, panleukopenia, and chlamydia (FVRCP-C) (Boehringer Ingelheim, St. Joseph, MO), was not on any medications, and had never been administered fluroquinolone antibiotics. Daily diet included Nebraska Brand Premium Feline Diet (75 grams) and a mouse (25 grams). Both cats were housed indoors in an approximately 10 × 10 foot enclosure under a 12 hour light -dark cycle.
Ophthalmic examination. Both affected African black-footed cats underwent complete physical, ophthalmic and neuro-ophthalmic examinations. Assessment of the male proband cat was completed under sedation and physical restraint by a board-certified veterinary ophthalmologist (JWP) and board-certified zoo veterinarian (WKS) at the MU CVM. Sedation was achieved with an intramuscular injection of ketamine (Pfizer, New York, NY, USA; 7 mg) and dexmedetomidine (Pfizer, New York, NY, USA; 0.55 mg), and 100% oxygen was provided via mask. A 24-gauge catheter was placed in the saphenous vein and Lactate Ringer's Solution (LRS) was administered at a rate of 10 ml/kg/hour. The animals' vital signs (heart rate, respiratory rate, and temperature) were evaluated every 5 minutes. Following the procedure, the sedation was reversed with an intramuscular injection of atipamazole (Pfizer, New York, NY, USA; 0.55 mg). Evaluation of the female sibling was performed under general anesthesia by a board-certified veterinary ophthalmologist (EBB) and zoo veterinarian (MS) at the Metropolitan Veterinary Hospital, Akron, OH. Anesthesia was 'chamber induced' to effect and maintained with sevoflurane delivered in 100% oxygen via an endotracheal tube. Ophthalmic diagnostics included Schirmer tear test I (Schirmer Tear Test, Schering-Plough Animal Health, Union, NJ, USA), rebound tonometry (TonoVet ® ; Icare, Espoo, Finland) and fluorescein stain (Flu-Glo ® Fluorescein sodium 1 mg stain; Akorn Inc., Buffalo Grove, IL, USA). Neuro-ophthalmic examinations included evaluation of direct and consensual pupillary light reflexes (PLRs) using a 3.5 V Finnoff transilluminator (Welch Allyn, Skaneateles Falls, NY, USA), assessment of dazzle reflex using white light directly from the fiberoptic cable of a slit lamp biomicroscope (SL-2; Kowa Optimed Inc, Torrance, CA, USA), behavioral testing of vision using menace response, and measurement of horizontal pupil diameter with Jameson caliper. Diffuse and focal illumination with a portable hand-held slit-lamp biomicroscope (SL-14; Kowa Col Ltd., Tokyo, Japan) was utilized to examine the anterior segment of both animals. Indirect ophthalmoscopy was performed with a wireless binocular indirect headset (12,500; Welch Allyn Inc., Skaneateles Falls, NY, USA) and a condensing lens (Pan Retinal 30D; Volk Optical Inc., Mentor, OH, USA). Mydriasis was not induced pharmacologically for examination of the fundus due to marked pupillary dilation present in both patients. The fundus of the male proband was photographed with a NIDEK NM-100 digital fundus camera (NIDEK CO. LTD., Aichi, Japan).

Additional diagnostic procedures.
A complete blood hematology, chemistry panel and urinalysis were completed in both affected African black-footed cats. Plasma taurine levels were also determined. Both affected cats completed electroretinography examinations (ERG) at their respective clinics. The ERG was completed in both cats under general anesthesia based on previous studies 60 . Both cats were dark adapted for at least 20 minutes prior to recording. Briefly, the ERG was performed utilizing a monopolar electrode-contact lens (ERG-jet; Nicolet Instruments, Madison, WI, USA) applied to the cornea with hypromellose 2.5% gel as a coupling agent (Gonak; Akorn, Inc., Buffalo Grove, IL, USA) and male subdermal platinum needle electrodes (FD-E2-24; Astro-Medical, Inc. Warwick, RI, USA). The ERG for the male proband was performed with a Handheld Multispecies ERG (HMsERG) (HMsERG Model 1000; RetVetCorp, Columbia, MO, USA) using the Quick Ret Check Protocol. The ERG for the female sibling was performed with a Retinographics unit (Retinographics Inc., Norwalk, CT). For both cats, mixed rod and cone responses were assessed with a standard flash intensity of 3 cd.s/m 2 . Mixed rod and cone responses of the male cat were also evaluated with a high intensity flash (10 cd.s/m 2 ). Indirect blood pressure measurement was performed under general anesthesia in both cats.
All additional cats in the study with available DNA samples (n = 14), including the parents of the siblings, were reported to behave in a visually appropriate manner and failed to demonstrate any visual deficits.
Genetic analyses. African black-footed cats (n = 16; eight males and eight female cats) from the captive populations at seven zoological parks were evaluated (Table 1). An aliquot of the DNA was submitted for commercial genotyping to the Veterinary Genetics Laboratory at the University of California, Davis, to screen for the three known variants associated with PRA in cats: AIPL1, CRX and CEP290. DNA for WGS was isolated by organic extraction from EDTA anti-coagulated whole blood collected from a trio of African black-footed cats, including the obligate carrier sire (205) and dam (208), and affected female offspring (247) ( Table 1). The WGS libraries and sequencing were performed at the University of Missouri DNA Core on an Illumina HiSeq 2000 as previously described 11 . Cat WGS filtering and genome-wide variant calling were performed as previously described by Maverix Biomics, Inc. (San Mateo, CA) 11 and as part of the 99 Lives whole genome sequencing initiative (http://felinegenetics.missouri.edu/99lives). Briefly, the WGS from the trio was aligned to the available annotated cat genome reference (Felis_catus 6.2 genome assembly) to identify variants 50 . All variants detected in the black-footed cat trio were filtered based on the known individuals' phenotype and using the available variant dataset comprised of 52 whole genome sequenced cats, including one Pallas cat and 51 domestic cats. Variants were filtered as follows (1) considering variants only segregating in the black-footed cat trio, (2) excluding variants in common with black-footed cat and Pallas cat, (3) excluding variants shared in domestic and black-footed cats, (4) variants in concordant segregation with the autosomal recessive disease phenotype within the black-footed cat trio, and, (5) identifying variants in known retinal genes. A list of genes causing retinal diseases from the RetNet database (www.RetNet.org) was obtained and variants were visually investigated within those genes. Variants were prioritized by their effect on the protein using SNPeff (http://snpeff.sourceforge.net/ SnpEff_manual.html#eff).

IQ calmodulin-binding motif-containing protein-1 (IQCB1) variant genotyping. A direct
sequencing assay was performed to screen a possible causative variant within the candidate gene IQCB1 in all available black-footed cats with known phenotype (n = 16, Table 1). Two primers were designed on the variant flanking regions (IQCB1-FCCTCTCCTCATACACACTCAGAGTTA, IQCB1-R CGGTGTCTTTGGAGATAGTCATC) and PCR was performed at 58 °C annealing temperature. Cycles, reagents final concentrations, PCR purification and direct Sanger sequencing were performed as previously described 61 . Sanger sequencing was completed at the DNA Core at the University of Missouri on an ABI 3730 DNA Analyzer (Applied Biosystems, Foster City, CA).

Mass spectroscopy variant genotyping. Candidate variants associated with a missense change in a
RetNet gene were selected for genotyping. Assays were designed using the MassARRAY Assay Design included in the package Typer 4.0 (Agena Bioscience, San Diego, CA) according to manufacturer's instructions. A multiplex consisting of seven SNPs in six different genes (CHD3, CNGB1, DTHD1, RBP3, USH1C and CDH23) was selected (Supplementary Table 1), genotyping was performed on all available samples using the Agena Biosciences iPLEX Gold Genotyping reagent set and products were typed with the MassARRAY System with Nanodispenser RS1000 (Agena Bioscience).
Scientific RepoRts | 7:43918 | DOI: 10.1038/srep43918 Assay performance and validation. All data was inspected with Typer Analyzer to assess assay quality and assays with weak or ambiguous signals were discarded. Next, the Autocluster algorithm of Typer Analyzer was selected and only assays with "Conservative" or "Moderate" calls were considered. Moreover, the three WGS black-footed cats functioned as positive controls for assay validation.