Open Access

Whole genome sequencing in cats, identifies new models for blindness in AIPL1 and somite segmentation in HES7

  • Leslie A. Lyons1Email author,
  • Erica K. Creighton1,
  • Hasan Alhaddad2,
  • Holly C. Beale3,
  • Robert A. Grahn4,
  • HyungChul Rah5,
  • David J. Maggs6,
  • Christopher R. Helps7 and
  • Barbara Gandolfi1
BMC Genomics201617:265

https://doi.org/10.1186/s12864-016-2595-4

Received: 5 November 2015

Accepted: 16 March 2016

Published: 31 March 2016

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Keywords

Aryl-hydrocarbon-interacting receptor protein-like 1 Domestic cat Felis silvestris catus Hairy and Enhancer of Split 7 LCA4 Progressive retinal atrophy

Background

Whole genome sequencing (WGS) is becoming the standard of health care in humans [19], and is commonly used to discover causal variants in children with otherwise undiagnosed congenital defects [1013]. The 50-h and now the 26-h genome efforts have demonstrated how genome medicine can be applied to health management for acute care patients with time-critical morbidity and mortalities [10, 11]. The improved efficiency and lower costs of WGS now apply to all animal species, including pet cats. Over 80 million cats are owned in the USA and their roles as family-members and their health care are increasingly priorities for owners [14, 15]. Clinical trials of client-owned animals will advance both veterinary and human medicine [16]. Research colonies of cats have demonstrated their role as biomedical models for the development of gene and enzyme therapies [17, 18]. Indeed, gene therapy for canine Leber’s congenital amaurosis (LCA) has restored vision to the blind [19] and a variety of LCA clinical trials are underway in human patients [2023].

Congenital amaurosis (MIM: 204000) was first described by Leber in the 1800’s [24, 25] and is now recognized as a heterogeneous group of early-onset childhood retinal dystrophies. LCA is classified as the most severe form of retinopathy [26]. LCA has a worldwide prevalence of 1 in 81,000 to 1 in 33,000 newborn babies, accounting for ≥ 5 % of all inherited retinopathies and approximately 20 % of children attending schools for the blind around the world [2729]. Nineteen different genes are currently associated with the various forms of LCA [30, 31] (http://omim.org/phenotypicSeries/PS204000). Variants at an autosomal recessive locus, LCA4 (MIM: 604393), are caused by DNA variants in AIPL1 (MIM: 604392) on human chromosome 17p13.1 [3234]. Aryl-hydrocarbon-interacting receptor protein-like 1 (AIPL1) variants cause approximately 7 % of LCA worldwide and may also cause dominant retinopathy [33, 34].

The Persian cat has been proposed as an animal model for retinal degeneration, specifically LCA [3537]. Previous studies have localized the Persian cat progressive retinal atrophy to a 1.3 Mb region on cat chromosome E1, which is homologous to human chromosome 17 [35]. At the time of analysis, the publically available reference sequence, an Abyssinian cat, had over 11 Mb of genomic sequence in unplaced scaffolds [38]. The incomplete reference obstructed a full examination of an LCA gene, AIPL1 within the haplotype block due to a gap in the cat genomic sequence.

WGS is feasible for domestic cats and is a highly efficient method to genetically examine several interesting phenotypes, in one cat, that do not have overlapping or compromising health effects. To maintain diversity in the Persian cat colony segregating for PRA, several innocuous phenotypes that are specific to pedigreed cats were introduced by outcrossing to cats of unrelated breeds. The bobbed tail of the Japanese bobtail cat breed results from abnormal caudal vertebrae morphology [39]. In addition, bobbed-tailed cats are typically lack a vertebra from the thoracic or lumbar regions, promoting this cat breed a model for spinal column development and somite segmentation [4044]. Genes involved with somite segmentation are responsible for a severe disease in humans and dogs, spondylocostal dysostosis [4547].

In this study, the Persian cat PRA is defined as a new model for LCA and the bobbed tail cat is a new model for spinal cord development. Using WGS of a trio of cats segregating for multiple traits, a gene within a ~1.3 Mb targeted region of interest from a previous genome-wide association study (GWAS) was implicated for Persian cat PRA. Without a priori knowledge of the GWAS targeted region, the LCA variant would have been identified based on proper segregation in the trio and the present of a stop gain variant in a known LCA gene. This feline PRA model is now available to the vision science community for gene therapies studies, and has potential to support both veterinary and human medicine. In addition, a gene associated with somite segmentation was implicated for the bobbed tail phenotype. The genomic resources for the domestic cat are now making WGS-based disease variant discovery feasible in Felis catus, further supporting the domestic cat as a comparative biomedical model for human studies.

Methods

Pedigree and clinical description

The breeding colony of Persian cats and their clinical diagnoses of PRA and bobbed tail are previously described [36, 39]. Cats were maintained under the Institutional Animal Care and Use protocols 11977, 15117 and 16691 at the University of California – Davis (UC Davis) and protocol 7808 at the University of Missouri. No additional ethical committee approval was required or requested over and above general adherence to protocols. PRA disease status was confirmed during a complete by ophthalmic examination performed by board-certified veterinary ophthalmologists from the Ophthalmology Service of the UC Davis Veterinary Medical Teaching Hospital as described [36]. All cats were over 16 weeks of age, therefore disease status was also confirmed by behavioral observations as described [36]. Tail phenotypes were confirmed by radiographs and or palpation.

Cats (N = 1740) from 40 different breeds, as well as random bred cats, were genotyped for the variants of interest (Additional file 1: Table S1). The majority of cats included were Persian or Persian-derived breeds (British shorthair, exotic, Himalayan, Scottish fold, Selkirk rex) [48, 49] (n = 1213), including 85 cats (59 sighted and 26 blind cats) from an established PRA colony [35, 36], 84 biased and 1044 unbiased Persian and Persian family samples, Table 1). A subset of these cats (n = 271) from 27 different breeds and random bred cats were genotyped for the HES7 variant, including 36 colony cats (24 normal tailed and 6 bobbed-tail cats). A cat trio was selected from the larger, extended pedigree segregating as for the maximum number of traits integrated into the colony, including two forms of PRA, bobbed-tail, dorsally curled ears, hair length and six coat colors (Fig. 1). The trio included a blind sire, a non-carrier dam with a bobbed tail and a carrier offspring with a bobbed tail, Fig. 1).
Table 1

AIPL1 c.577C > T Genotypes in Persian and Related Breeds

Typea

Population

No.

Wildtype (CC)

Carrier (CT)

Affected (TT)

Biased

Pedigreeb

85

19

40

26

 

Persian

51

43

8

0

 

British Shorthair

2

2

0

0

 

Exotic

29

26

3

0

 

Himalayan

2

2

0

0

 

Sub-total

84

73

11

0

Unbiased

Persian

707

685

22

0

 

British Shorthair

111

111

0

0

 

Exotic

103

102

1

0

 

Himalayan

42

42

0

0

 

Scottish Fold

68

67

1

0

 

Selkirk Rex

13

13

0

0

 

Sub-total

1044

1020

24

0

Total

 

1213

1112

75

0

aType implies if the cat samples were biased because an owner submitted a sample specifically for the PRA genotyping or unbiased because the laboratory did a population screen of samples submitted for other genetic testing. bPedigree cats are related and not considered for allele frequency calculation. British Shorthairs, exotic, Himalayan, Scottish fold and Selkirk rex are cat breeds derived from Persians and are members of the Persian-family of cats [48]

Fig. 1

First Trio of Cats WGS for Disease Variant Identification. (Left) Pedigree of cats segregating for the AIPL1 variant for Persian progressive retinal atrophy (PRA). Circles represent females, squares represent males, open symbols indicate phenotypically normal cats; solid symbols indicate Persian PRA-affected cats, half-filled symbols indicate obligate carrier cats. The arrows indicate the probands. AIPL1 c.577C > T genotypes are indicated above each cat. Cats with an asterisk “*” are the cats used for whole genome sequencing. (Right) The trio of cats segregating for the Persian PRA was whole genome sequenced to potentially identify the causal variant for this animal model of LCA. Cats included (top) a wildtype bobbed tail queen (S13230), (middle) an PRA-affected normal tail sire (S14056), and (bottom) an PRA-obligate carrier female offspring (S16628) which also displays the bobbed tail. The bobbed tail trait was introduced by female cat 11662, a pedigreed cat of the Japanese Bobtail breed

Cat whole genome sequencing

DNA for WGS was isolated by organic extraction from ~ 3 ml of EDTA anti-coagulated whole blood that was collected by jugular venipuncture of the trio [50]. DNA quality and quantity was visualized using ethidium bromide staining after agarose gel electrophoresis. Approximately 4 μg of high molecular weight DNA from each cat was submitted to the University of Missouri DNA Core for sequencing library preparation and NextGeneration sequencing. Sequencing libraries were constructed following the manufacturer’s protocol with reagents supplied in Illumina’s TruSeq DNA PCR-Free sample preparation kit (#FC-121-3001) (Illumina, San Diego, CA). Briefly, 1–2 μg of genomic DNA was sheared using standard Covaris (Woburn, MA) methods to prepare a 350 bp and a 550 bp library for each individual cat of the trio. The resulting 3′ and 5′ overhangs were converted to blunt ends by an end repair reaction which uses 3′ to 5′ exonuclease activity and polymerase activity. The desired size of fragment (~350 or 550 bp) was selected by sample purification beads (AMPure XP). A single adenosine nucleotide was added to the 3′ ends of the blunt fragment followed by the ligation of Illumina indexed paired-end adapters. The adaptor ligated libraries were purified twice with sample purification beads. The purified libraries then were quantified with a Qubit assay (Life Technologies, Carlsbad, CA) and library fragment size confirmed by Fragment Analyzer (Advanced Analytical Technologies, Inc., Ames, IA). Finally, the libraries were diluted and sequenced on a paired-end, 100 bp read length run according to Illumina’s standard sequencing protocol for the HiSeq 2000 (Illumina). The 350 bp and 550 bp libraries from all three cats were pooled and analyzed across nine lanes of a HiSeq 2000 (Illumina, Inc.). An average of ~30X sequencing coverage per cat was expected based on typical HiSeq output per lane and the use of PCR-free libraries.

Genome alignment and variant calling

The demultiplexed 100 bp paired-end reads generated from the two sequencing libraries for each cat were transferred to Maverix Biomics (San Mateo, CA) for variant detection. Raw sequencing reads were checked for potential sequencing issues and contaminants using FastQC [51]. Adapter sequences, primers, Ns, and reads with quality scores below 28 were trimmed using Trimmomatic [52]. Reads with a length < 20 bp after trimming were discarded. After reassessment for quality improvement using FastQC [53], trimmed reads were mapped to the domestic cat reference genome, Felis catus 6.2 (http://www.ncbi.nlm.nih.gov/assembly/320798), using BWA-MEM [54]. Duplicated reads were identified and removed using SAMBLASTER [55]. Sequencing coverage across the genome for each cat was measured using DepthOfCoverage from GATK [56]. Only bases with a Phred quality score above 28 and reads with a mapping quality score above 20 were included in the coverage analysis. FreeBayes [57] was used to detect variants from the read alignments. Variants with quality scores below 20 and homopolymers were filtered from the data. The remaining variants were annotated with dbSNP ID [58] and effect predictions using SnpEff [59] based on the Ensembl gene model [60] of the reference genome. The reported variant effects and severities (e.g. high or moderate) are generated by SnpEff using variant calls, the reference genome and transcript definitions as input. Transcript definitions were obtained from the publicly available Ensembl website (www.ensembl.org). Variant effect predictions were summarized and reformatted into tab-delimited annotation files for easy interpretation and filtering. Data tracks that contained genotypes of the cat trio were prepared for visualization and data comparison in the UCSC Genome Browser [61]. Annotated variants were loaded into the variant database of the Maverix Analytic Platform that included an interactive variant exploration tool for dynamic filtering of variants based on selected effects, quality scores, gene loci, and phenotypes. To verify dynamic filtering accuracy, known phenotypic coat color variants within the trio, reviewed in Lyons [62], were verified and visually validated, including the loci for Agouti (ASIP), Brown (TYRP1), Dilution (MLPH), Long fur (FGF5) and ventral white Spotting (KIT).

Genome-wide and haplotype variants analysis

A previous GWAS localized the PRA phenotype from position 713,552 - 2,076,816 on cat chromosome E1, therefore the region was scanned for disease-causing genes and variants [35]. In addition, a list of genes causing retinal diseases from the RetNet database (www.RetNet.org) was obtained and variants were visually investigated within those genes. Various studies on somite segmentation suggested HES7 as a candidate for cat bobbed tail [40, 41, 43]. The gene is also located on cat chromosome E1, from position 2,816,901 to position 2,819,548. For the tabulations in Table 2, vcftools v0.2.12a were used to identify variants in the GWAS haplotype or in Retnet genes [63]. Segregating variants were defined as those in which the dam was homozygous for the reference allele, the sire was homozygous for the alternate allele and the offspring was heterozygote for PRA and for bobbed-tail, the dam and offspring were heterozygote for the alternative allele and the sire was homozygous for the reference allele.
Table 2

WGS Variants Identified in 30x Coverage of a Persian Cat PRA Trio

 

Analysis type

Variant*

WGS trio

Trio Segregation

Retnet genes

GWAS haplotype

SNP

14,215,442

487,324

74,852

9607

MNP

822,020

11,856

4214

644

Insertion

1,390,990

16,656

7439

636

Deletion

1,279,028

22,134

6760

663

Mixed

304,341

3314

1588

198

Total

18,011,821

541,284

94,853

11,748

*Each alternate allele identified during variant calling is considered separately for tabulation of variant types, using the definitions at http://snpeff.sourceforge.net/SnpEff_manual.html#eff. The analyses are based on feline reference genome sequence (V6.2). The GWAS identified an associated haplotype on cat chromosome E1. For Trio Segregation, the variants segregated with the disease phenotype as well, not just parent to offspring

AIPL1 RNA analysis and variant genotyping

Total RNA was isolated from the retina of a sighted cat using the RNA mini kit (Invitrogen, Carlsband, CA). Complementary cDNA was synthesized as previously described (Gandolfi et al. [64]). The cDNA template was subject to PCR using 10 μM of the PCR RNA forward primer (Additional file 1: Table S2) combined with the Invitrogen polyT primer provided in the kit. The PCR conditions were: 94 °C x 4 min followed by 10 cycles conducted at 94 °C for 10 s, 66 °C for 30 s, and 68 °C for 2 min, and 20 additional cycles at 94 °C for 10 s, 62 °C for 30 s, and 68 °C for 2 min with an additional 20 s per cycle. The PCR products with appropriate lengths were purified and sequenced using four internal primers (Additional file 1: Table S2) as previously described [64].

To further correlate the suspected causative variants with the phenotypes, archival samples representing a subset of the colony cats and unrelated individuals from different breeds and random bred populations (Additional file 1: Table S1) were genotyped by direct Sanger sequencing of PCR-generated amplicons for AIPL1 exon 4 and the HES7 exon 1 confirmed the variants. Primer sequences are presented in Additional file 1: Table S2. PCR and thermocycling conditions were conducted as previously described with annealing at 58 °C [65]. PCR products were purified and sequenced as previously described [66]. Genotypes were also determined in broader cat populations (Table 1, Additional file 1: Table S1) by the commercial genetic typing services of the UC Davis Veterinary Genetics Laboratory using an allele-specific oligo assay and the University of Bristol, Langford Veterinary Services using a pyrosequencing-based assay (Additional file 1: Table S2). For the allele-specific assay, specificity of each primer was increased by incorporating a sequence mismatch and the addition of four mismatched bases that created a detectable product size variant (Additional file 1: Table S2). Amplicons were visualized on an ABI3730 DNA analyzer (Applied Biosystems, Foster City, CA) and analyzed using STRand software [67]. Primers for pyrosequencing were designed using PyroMark Assay Design Ver 2.0 (Qiagen, UK) (Additional file 1: Table S2). Pyrosequencing was undertaken after PCR amplification using GoTaq Master Mix (Promega, UK) of genomic DNA isolated from mouth swabs using the Nucleospin Blood kit (Macherey-Nagel, Germany) according to the manufacturer’s instructions (PyroGold, Qiagen) on a PyroMark Q24 (Qiagen). Pyrosequencing PCR was conducted using 95 °C for 2 min denaturation, followed by 38 cycles of 95 °C for 20 s, 58 °C for 40 s. For the HES7 variant, pedigreed cats were genotyped using MALDI-TOF by submission to GeneSeek (Neogene, Inc., Lincoln, NE).

Results

Availability of data and materials

Sequences were submitted to the NCBI Sequence Read Archive (http://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?cmd=show&f=sra_sub_expl) under BioProject PRJNA288177 and included all three cats: PRA – affected, normal - tailed sire (S14056: SRS1050848), a bobbed tail dam (S13230: SRS1053073), and an obligate PRA carrier bobbed tail offspring (S16628: SRS1051523). The complete AIPL1 CDS was obtained by RT-PCR of mRNA from the retina of a sighted cat and the sequence was submitted to the NCBI (accession #KP682504.1).

Whole genome sequence alignment, read depth and dataset statistics

A trio of cats segregating for PRA and bobbed tail, including a PRA – affected, normal - tailed sire (S14056: SRS1050848), a bobbed tail dam (S13230: SRS1053073), and an obligate PRA carrier bobbed tail offspring (S16628: SRS1051523), was selected from the colony pedigree (Alhaddad et al. 2014) for WGS (Fig. 1). Over one billion 100 bp paired - end raw reads were obtained for each cat (1,027,975,352 in S13230, 1,100,610,902 in S14056, and 1,246,386,562 in S16628). The sequencing results were of relatively high quality with approximately 19 % of the reads removed after filtering out adapter sequences and those with low quality scores (Additional files 2 and 3: Figures S1-S2). Because no amplification was required to prepare the PCR-free libraries, an average of only 0.017 % duplicated reads among the six libraries was observed, which minimized the removal of duplicated read alignments used for variant calling and retained the maximum of read coverage. The mapping rates of the libraries to the cat reference genome were between 86.6 and 86.8 %. Average depth of coverage for each cat in the trio was 28.3X, 30X, and 34.1X with 94.3, 93.6 and 96.1 % of the genome having above 15X - coverage (Additional file 1: Table S3, Additional file 4: Figure S3).

Variants analysis - genome-wide, candidate gene and associated haplotype

Compared to the reference, over 18,000,000 variants were detected by genome-wide sequencing the cat trio. Details for variant counts and their effects are presented in Tables 2 and 3. The number of variants that segregated concordantly with the PRA phenotype reduced to approximately 541,000 when accounting for segregation – homozygous affected sire, obligate carrier offspring and absent in the dam and reference sequence. A list of over 220 genes associated with blindness is available in the RetNet database (www.RetNet.org). Nearly 95,000 variants were found within these candidate genes and nearly 12,000 variants within the previously associated haplotype [35] (Table 3).
Table 3

WGS variants Identified in 30x Coverage of a Persian Cat PRA Trio

   

Analysis typeb

  

Variant impacta

Functional Class

WGS trio

Trio Segregation

RetNet genes

GWAS haplotype

High

Stop gain

379

11

6

1

 

Start/Stop loss

85

5

2

-

 

Splice donor/acceptor

1262

9

9

1

 

Exon deletion

2

-

-

-

 

Frameshift

2254

18

20

4

 

Rare amino acid

-

-

-

-

Moderate

Codon alteration

971

15

11

4

 

Missense

35,444

1014

544

41

 

Splice branch

-

-

-

-

 

5′ or 3′ UTR Deletion

1

-

-

-

Low

 

64,791

2072

1155

105

Modifier

 

16,876,221

538,148

87,816

11,018

aImpact determined by snpEff as defined by http://snpeff.sourceforge.net/SnpEff_manual.html#eff. b The analyses is based on feline reference genome sequence (V6.2). Effect counts are higher than variant counts because they include the effects of each alternate allele on each nearby gene isoform. GWAS haplotype are the variants identified within the haplotype region. For Trio Segregation, the variants segregated with the disease phenotype as well, not just parent to offspring

Considering SNP effects and their impacts (Table 3), 43 high impact variants segregated concordantly with the PRA disease within the trio. Thirty-seven highly damaging variants were within the RetNet candidate genes (Table 3). The 1.3 Mb haplotype block [35], harbored six highly damaging variants, including the AIPL1 variant, and 45 genic moderate variants (Table 3). From the phenotypes and specific genotypes of the cats, known mutations for Agouti, Brown, Color, Dilution, Long and Spotting were correctly identified by the variant calling (data not shown) [62]. The stop gain identified in AIPL1 was considered the highest priority candidate gene mutation and was further genotyped. Approximately 101 high and moderate impact variants that segregated properly with the bobbed-tail phenotype in the trio were identified in 35 genes. The HES7 gene was examined as an obvious candidate and a missense variant was identified in the coding region.

Gene analysis, experimental validation and mutation genotyping

A stop gain was identified in AIPL1, a gene that causes LCA, at position c.577C > T of the cat coding DNA sequence. The complete AIPL1 CDS was obtained by RT-PCR of mRNA from the retina of a sighted cat and the sequence was submitted to the NCBI (accession #KP682504.1). The complete gene coding sequence was obtained, including the complete sequence of exons 4, 5 and 6, which are not present in the current feline genome assembly V6.2. No polymorphisms were identified in the obtained sequence. The cDNA amplification generated a 987 bp long CDS and in silico sequence translation of feline wild-type AIPL1 generated from RNA (GenBank: KP682504) predicted a length of 329 amino acids, matching Mus musculus predicted AIPL1 (Additional file 5: Figure S4). The mutated AIPL1 transcript encoded a predicted p.Arg193* protein truncation, preventing translation of approximately 40 % of the protein when compared to the wild-type sequence. When compared to human, the feline and mouse predicted AIPL1 amino acid sequences were 56 amino acids shorter, confirming that the proline-rich protein C-terminal is absent in these species. Feline DNA identity with human coding sequence (GenBank: NM_014336.4) was 76 and 80 % when compared to mouse (GenBank: XM_053245.2) while feline protein identity was 77 % when compared with both species. Feline protein sequence alignment suggested the presence of highly conserved domains when compared to the human segment: an N-terminal FK506 binding protein (FKBP-like domain) and three C-terminal tetratricopeptide repeats (TPR1, TPR2, TPR3). Homology between the feline and human FKBP, TPR1, TPR2 and TPR3 domains was 91, 94, 97 and 97 % respectively. The truncated cat protein lacked part of the TPR1 domain and the entire TPR2 and TPR3 domains.

A c.5A > G causing a p.V2A amino acid substitution was the only coding region alteration identified in HES7, which also segregated properly in the trio and was not present in the reference sequence. Sequence homology of cat (XM_003996191.2) and human (NM_001165967.1) is 89.9 % and protein identity is 91.4 %. The p.V2A amino acid change is prior to the defined basic helix-loop-helix and orange domains. The cat protein is also proline rich from p.128 to the C – terminus.

Cats (N = 1740, distributed as follows: unbiased individuals = 1558, Persian colony pedigree = 85, biased individuals = 97) from 40 different breeds, including 36 random bred cats and 61 cats with unknown origins were genotyped for the AIPL1 variant (Table 1, Additional file 1: Table S1). A sub-set of cats (n = 271) was genotyped for the c.5A > G HES7 variant. Twenty-six of 85 cats from the multi-generational pedigree and verified by ophthalmic examination to be phenotypically affected were confirmed as homozygous for the c.577C > T variant, including all three pedigree founders (Fig. 1). For 59 cats verified by ophthalmic assessment to be unaffected, 19 were homozygous wild-type and 40 heterozygous (Table 1, Additional file 1: Table S1). A majority of sighted offspring in the pedigree were expected to be heterozygous because most of the breedings were performed as test crosses of an obligate carrier to an affected mate (Fig. 1). Cats from breeds related to the Persian breed and at risk for the AIPL1 variant included Persian, Scottish fold, Selkirk rex, British shorthair, burmilla, exotic shorthair, and Himalayan cats. No cats from a breed genetically unrelated to Persian cats had the AIPL1 variant. No cats except those in the colony with bobbed tails (n = 6) and cats of the Japanese bobtail breed (n = 14) had the HES7 variant. All colony cats were heterozygous for the HES7 variant, while pedigreed Japanese bobtails were homozygotes. To estimate the variant frequency within the Persian and Persian family of breeds, 1044 cats were genotyped by two commercial service laboratories (Table 1, Additional file 1: Table S1). These samples were originally submitted for other genetic testing services. In this unbiased sampling, 24 cats were heterozygous, including 22 Persians, one Scottish fold and one exotic, suggesting an allele frequency of 1.15 %. Additionally, 84 cats were specifically submitted for PRA genotyping by cat owners, eleven were heterozygous and none were homozygous for the variant, suggesting a 6.55 % allele frequency. Because these samples were submitted by owners for PRA genotyping, the sampling is biased and the allele frequency is likely high, especially because some cats may be siblings in a litter.

Discussion

Besides traditional rodent models, WGS has been used in cattle [68] and dogs [6971] to identify disease mutations, however, the datasets have included hundreds of individuals for comparison. More recently, a WGS of a single cat was used to confirm a previously identified phenotypic variant for White [72] and a single cat compared to a variant database including 18 cats identified a feline model for Congenital Myasthenic Syndrome [47, 73]. However, the cat variant database needs to be improved for continued success, particularly for genetic studies in non-pedigreed, random bred cats [38, 74, 75]. WGS of more than one individual with an inherited disease that is identical by descent can overcome the lack of a deep variant database.

The present study represents the first comprehensive WGS effort to identify causal mutations in the domestic cat, allowing assessment of the efficiency of the WGS and variant calling process using the current cat reference assembly. In addition, other traits also segregating in this cat trio, including the bobbed tail of the Japanese bobtail breed, dorsally curled ears of the American curl, hair length, and at least six color variants were available for analysis. A targeted sequencing depth 30X genomic coverage was selected because of the number of variants that can be identified and the fraction of the genome that is callable plateaus above this depth of coverage [76]. Deep WGS of the trio of cats, including one non - PRA affected bobbed tail parent, one PRA - affected non – bobbed tail parent, and one obligate PRA carrier, bobbed tail offspring was conducted to fill sequence gaps in the targeted GWAS region for PRA in Persian cats and to support accurate variant calling. Variation was high in this trio because the family included outcrosses to several different cat breeds, including Persian, Oriental shorthair, American curl, Somali, and Bengal. Of the approximately 18 million variants identified, only 541,284 (3.0 %) had genotypes that segregated with the Persian cat PRA phenotype, 94,853 were in genes associated with vision loss (RetNet), and 11,748 that were within the GWAS haplotype region. By examining high-impact variants, stop gains were reduced from 379 overall to 11 segregating within the trio and with the phenotypes. Six stop gains were in RetNet genes and only one was in the haplotype region within an obvious LCA candidate gene, AIPL1. Thus, without the GWAS and a list of candidate genes, eleven stop gains would likely have been prioritized and examined over the other 32 high-impact variants segregating in the trio. Without the trio and other cat genomes for comparison, many thousands of variants could have been considered of high impact and would have required prioritization.

The DNA variant identified in the present study is a few bases upstream of a gap in the reference sequence. RNA sequencing provided the complete coding DNA sequence of AIPL1. AIPL1 mutations affecting vision are grouped into three classes: the first class of missense mutations is located in the N-terminus, the second class of missense and stop mutations is located in the TPR motifs, and the third class of mutations includes small in-frame deletions located in the C-terminus. The first and second classes of mutations are associated with autosomal recessive LCA, while mutations in the third class appear to be associated with autosomal dominant cone-rod dystrophy and juvenile RP [34].

The abnormal vertebral presentation in the Japanese bobtail cats was an unexpected finding [39] that could be clearly associated with genes involved with somite segmentation [4044]. The cat trio was purposely selected to segregate for a variety of traits to assist interpretation of the efficiency of the cat reference assembly and WGS efforts. The identified c.5A > G missense variant predicts a p.V2A alteration in HES7 that appears to act dominantly to disrupt feline vertebral column development, leading to absence of either one thoracic or lumbar vertebrae and misplacement of ribs [39]. Homozygous cats have no additional abnormalities, but do have more predictable “show quality” tails. In humans, HES7 missense changes p.R25W, p.I58V and p.D186Y cause autosomal recessive spondylocostal dysostosis (SCDO4, OMIM:613686) [45, 46].

Segregation of the AIPL1 and HES7 variants was confirmed in an extended pedigree of Persian cats with PRA and the introduced bobbed-tail trait. The PRA founders of the pedigree were from three geographically wide-spread locations in the USA, and from distinct lines of Persian cats. These cats were ascertained as adults in 1999 for the breeding colony, thus this form of PRA has been in the Persian population for at least 15 years. The Persian breed is the largest by population in the Cat Fanciers’ Association with over 60 % of registered cats being represented by Persians and related breeds, such as exotic shorthairs, Himalayans, Selkirk rex, Scottish fold, and British shorthairs [77]. With an allele frequency of 1.15 %, and considering an average lifespan of 12 years and a population size of at least 500,000 pedigree cats from Persian-family breeds, approximately 66 may have vision loss due to the AIPL1 variant with approximately 11,367 carriers in the population. The HES7 variant was private to the Japanese bobtail breed and not identified in any normal tailed breeds, however, additional bobbed tailed breeds should be examined.

The cat eye has many structural and functional similarities with the human eye, including a dual retina with rod photoreceptors that greatly outnumber cones [78]. Due to similarities between the feline and human eye, the cat has been used extensively for research involving retinal structure and visual function. Two inherited feline retinal degenerations with known causal variants are defined in Abyssinian cats. The centrosomal protein 290 (CEP290) cat model [79] harbors variants also seen in approximately 30 % of humans with LCA10 [80] and the cone rod homeobox gene (CRX) is a LCA7 model [81]. The present study describes a third model, characterized by a rapidly progressive retinal dystrophy due likely to AIPL1 deficiency. This new model should advance understanding of LCA4 pathology.

Conclusions

The three cats reported here represent the first individuals to be sequenced as part of the 99 Lives Cat Genome Sequencing Initiative (www.felinegenetics.missouri.edu/99lives). Identified is a stop mutation in AIPL1 of the Persian cat at position c.577C > T of the CDS producing a p.Arg193*, as predicted by the human sequence, knocking out approximately 40 % of the normal protein. Additionally, a HES7 c.5A > G variant is associated with the bobbed tail of the Japanese bobtail breed. Cats have a longer life span than mice, affording researchers the opportunity for repeated and extended trials of gene and stem cell administration. Due to their increased longevity and larger globe relative to mice, and physiology more similar to humans, cats represent an efficient and effective model for LCA gene - and stem cell therapies.

Abbreviations

LCA: 

Leber’s congenital amaurosis

AIPL1

aryl-hydrocarbon-interacting receptor protein-like 1

WGS: 

whole genome sequencing

UC Davis: 

University of California – Davis

GWAS: 

genome-wide association study

PRA: 

progressive retinal atrophy

HES7

Hairy And Enhancer of Split 7

FKBP: 

FK506 binding protein

TPR: 

tetratricopeptide repeat

SNP: 

single nucleotide polymorphism

Declarations

Acknowledgements

Funding was provided by the National Center for Research Resources R24 RR016094 and is currently supported by the Office of Research Infrastructure Programs/OD R24OD010928, the University of Missouri – Columbia Gilbreath-McLorn Endowment, the Winn Feline Foundation (W06-020, W10 -014), the Phyllis and George Miller Trust (MT08-015), the Cat Health Network (D12FE-509) and the University of California – Davis, Center for Companion Animal Health (2008-36-F) (LAL). We acknowledge the assistance of the veterinary ophthalmologists at the UC Davis Veterinary Medical Teaching Hospital who examined some cats. We appreciate the assistance of Nicholas Gustafson and the dedication of the Persian cat breeders, particularly Kerrie Meeks.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Department of Veterinary Medicine and Surgery, College of Veterinary Medicine, University of Missouri - Columbia
(2)
College of Science, Kuwait University
(3)
Maverix Biomics, Inc.
(4)
Veterinary Genetics Laboratory, School of Veterinary Medicine, University of California - Davis
(5)
Graduate School of Health Science Business Convergence, College of Medicine, Chungbuk National University
(6)
Department of Surgical and Radiological Sciences, School of Veterinary Medicine, University of California - Davis
(7)
Langford Veterinary Services, University of Bristol

References

  1. McCarthy JJ, McLeod HL, Ginsburg GS. Genomic medicine: a decade of successes, challenges, and opportunities. Sci Transl Med. 2013;5(189):189sr4.PubMedGoogle Scholar
  2. Biesecker LG, Green RC. Diagnostic clinical genome and exome sequencing. N Engl J Med. 2014;371(12):1170.View ArticlePubMedGoogle Scholar
  3. Choi M, Scholl UI, Ji W, Liu T, Tikhonova IR, Zumbo P, et al. Genetic diagnosis by whole exome capture and massively parallel DNA sequencing. Proc Natl Acad Sci U S A. 2009;106(45):19096–101.View ArticlePubMedPubMed CentralGoogle Scholar
  4. Gomez CM, Das S. Clinical exome sequencing: the new standard in genetic diagnosis. JAMA Neurol. 2014;71(10):1215–6.View ArticlePubMedGoogle Scholar
  5. Lee H, Deignan JL, Dorrani N, Strom SP, Kantarci S, Quintero-Rivera F, et al. Clinical exome sequencing for genetic identification of rare Mendelian disorders. JAMA. 2014;312(18):1880–7.View ArticlePubMedPubMed CentralGoogle Scholar
  6. Sawyer SL, Hartley T, Dyment DA, Beaulieu CL, Schwartzentruber J, Smith A, et al. Utility of whole-exome sequencing for those near the end of the diagnostic odyssey: time to address gaps in care. Clin Genet. 2015. doi:10.1111/cge.12654.Google Scholar
  7. Valencia CA, Husami A, Holle J, Johnson JA, Qian Y, Mathur A, et al. Clinical impact and cost-effectiveness of whole exome sequencing as a diagnostic tool: a pediatric center’s experience. Front Pediatr. 2015;3:67.View ArticlePubMedPubMed CentralGoogle Scholar
  8. Virani A, Austin J. Diagnostic clinical genome and exome sequencing. N Engl J Med. 2014;371(12):1169–70.View ArticlePubMedGoogle Scholar
  9. Westerink J, Visseren FL, Spiering W. Diagnostic clinical genome and exome sequencing. N Engl J Med. 2014;371(12):1169.View ArticlePubMedGoogle Scholar
  10. Saunders CJ, Miller NA, Soden SE, Dinwiddie DL, Noll A, Alnadi NA, et al. Rapid whole-genome sequencing for genetic disease diagnosis in neonatal intensive care units. Sci Transl Med. 2012;4(154):135–54.Google Scholar
  11. Miller NA, Farrow EG, Gibson M, Willig LK, Twist G, Yoo B, et al. A 26-hour system of highly sensitive whole genome sequencing for emergency management of genetic diseases. Genome Med. 2015;7(1):100.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Soden SE, Saunders CJ, Willig LK, Farrow EG, Smith LD, Petrikin JE, et al. Effectiveness of exome and genome sequencing guided by acuity of illness for diagnosis of neurodevelopmental disorders. Sci Transl Med. 2014;6(265):265ra168.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Willig LK, Petrikin JE, Smith LD, Saunders CJ, Thiffault I, Miller NA, et al. Whole-genome sequencing for identification of Mendelian disorders in critically ill infants: a retrospective analysis of diagnostic and clinical findings. Lancet Respir Med. 2015;3(5):377–87.View ArticlePubMedGoogle Scholar
  14. APPMA. National Pet Owner’s Survey. Greenwich: American Pet Product Manufacturing Association; 2008.Google Scholar
  15. AVMA. US Pet Ownership and Demographics Sourcebook. Schaumburg: American Veterinary Medical Association; 2007.Google Scholar
  16. Kol A, Boaz A, Athanasiou KA, Farmer DL, Nolta JA, Rebhun RB, et al. Companion animals: Translational scientist’s new best friends. Sci Transl Med. 2015;7:308ps21.View ArticlePubMedGoogle Scholar
  17. Cotugno G, Annunziata P, Tessitore A, O’Malley T, Capalbo A, Faella A, et al. Long-term amelioration of feline mucopolysaccharidosis VI after AAV-mediated liver gene transfer. Mol Ther. 2011;19(3):461–9.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Ponder KP, O’Malley TM, Wang P, O’Donnell PA, Traas AM, Knox VW, et al. Neonatal gene therapy with a gamma retroviral vector in mucopolysaccharidosis VI cats. Mol Ther. 2012;20(5):898–907.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Bennicelli J, Wright JF, Komaromy A, Jacobs JB, Hauck B, Zelenaia O, et al. Reversal of blindness in animal models of Leber congenital amaurosis using optimized AAV2-mediated gene transfer. Mol Ther. 2008;16(3):458–65.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Bainbridge JW, Mehat MS, Sundaram V, Robbie SJ, Barker SE, Ripamonti C, et al. Long-term effect of gene therapy on Leber’s congenital amaurosis. N Engl J Med. 2015;372(20):1887–97.View ArticlePubMedPubMed CentralGoogle Scholar
  21. Gueven N, Faldu D. Therapeutic strategies for Leber’s hereditary optic neuropathy: a current update. Intractable Rare Dis Res. 2013;2(4):130–5.PubMedPubMed CentralGoogle Scholar
  22. Jacobson SG, Cideciyan AV, Roman AJ, Sumaroka A, Schwartz SB, Heon E, et al. Improvement and decline in vision with gene therapy in childhood blindness. N Engl J Med. 2015;372(20):1920–6.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Rakoczy EP, Narfstrom K. Gene therapy for eye as regenerative medicine? Lessons from RPE65 gene therapy for Leber’s Congenital Amaurosis. Int J Biochem Cell Biol. 2014;56:153–7.View ArticlePubMedGoogle Scholar
  24. Leber T. Ueber retinitis pigmentosa und angeborene amaurose. Graefes Arch Clin Exp Ophthalmol. 1869;15(3):1–25.View ArticleGoogle Scholar
  25. Leber T. Ueber anomale formen der retinitis pigmentosa. Graefes Arch Clin Exp Ophthalmol. 1871;17(1):314–41.View ArticleGoogle Scholar
  26. Foxman SG, Heckenlively JR, Bateman JB, Wirtschafter JD. Classification of congenital and early onset retinitis pigmentosa. Arch Ophthalmol. 1985;103(10):1502–6.View ArticlePubMedGoogle Scholar
  27. Alström CH, Olson O. Heredo-retinopathia Congenitalis: Monohybrida Recessiva Autosomalis. Hereditas. 1957;43;1–178.Google Scholar
  28. Schappert-Kimmijser J, Henkes HE, Van Den Bosch J. Amaurosis congenita (Leber). Arch Ophthalmol. 1959;61(2):211–8.View ArticleGoogle Scholar
  29. Stone EM. Leber congenital amaurosis–a model for efficient genetic testing of heterogeneous disorders: LXIV Edward Jackson Memorial Lecture. Am J Ophthalmol 2. 2007;144(6):791–811. e6.View ArticleGoogle Scholar
  30. Pattnaik BR, Shahi PK, Marino MJ, Liu X, York N, Brar S, et al. A novel KCNJ13 nonsense mutation and loss of Kir7.1 channel function causes Leber congenital amaurosis (LCA16). Hum Mutat. 2015;36(7):720–7. doi:10.1002/humu.22807. Epub 2015 May 20.View ArticlePubMedGoogle Scholar
  31. Tan MH, Mackay DS, Cowing J, Tran HV, Smith AJ, Wright GA, et al. Leber congenital amaurosis associated with AIPL1: challenges in ascribing disease causation, clinical findings, and implications for gene therapy. PLoS One. 2012;7(3):e32330.View ArticlePubMedPubMed CentralGoogle Scholar
  32. Perrault I, Rozet J-M, Gerber S, Ghazi I, Leowski C, Ducroq D, et al. Leber congenital amaurosis. Mol Gen Metab. 1999;68(2):200–8.View ArticleGoogle Scholar
  33. Sohocki MM, Bowne SJ, Sullivan LS, Blackshaw S, Cepko CL, Payne AM, et al. Mutations in a new photoreceptor-pineal gene on 17p cause Leber congenital amaurosis. Nature Genet. 2000;24(1):79–83.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Sohocki MM, Perrault I, Leroy BP, Payne AM, Dharmaraj S, Bhattacharya SS, et al. Prevalence of AIPL1 mutations in inherited retinal degenerative disease. Mol Gen Metab. 2000;70(2):142–50.View ArticleGoogle Scholar
  35. Alhaddad H, Gandolfi B, Grahn RA, Rah HC, Peterson CB, Maggs DJ, et al. Genome-wide association and linkage analyses localize a progressive retinal atrophy locus in Persian cats. Mamm Genome. 2014;25(7–8):354–62.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Rah H, Maggs DJ, Blankenship TN, Narfstrom K, Lyons LA. Early-onset, autosomal recessive, progressive retinal atrophy in Persian cats. Invest Ophthalmol Vis Sci. 2005;46(5):1742–7.View ArticlePubMedGoogle Scholar
  37. Rah H, Maggs DJ, Lyons LA. Lack of genetic association among coat colors, progressive retinal atrophy and polycystic kidney disease in Persian cats. J Feline Med Surg. 2006;8(5):357–60.View ArticlePubMedGoogle Scholar
  38. Montague MJ, Li G, Gandolfi B, Khan R, Aken BL, Searle SM, et al. Comparative analysis of the domestic cat genome reveals genetic signatures underlying feline biology and domestication. PNAS. 2014;111(48):17230–5.View ArticlePubMedPubMed CentralGoogle Scholar
  39. Pollard RE, Koehne AL, Peterson CB, Lyons LA. Japanese Bobtail: vertebral morphology and genetic characterization of an established cat breed. J Feline Med Surg. 2015;17(8):719–26.View ArticlePubMedGoogle Scholar
  40. Bessho Y, Sakata R, Komatsu S, Shiota K, Yamada S, Kageyama R. Dynamic expression and essential functions of Hes7 in somite segmentation. Genes Dev. 2001;15(20):2642–7.View ArticlePubMedPubMed CentralGoogle Scholar
  41. Fujimuro T, Matsui T, Nitanda Y, Matta T, Sakumura Y, Saito M, et al. Hes7 3’UTR is required for somite segmentation function. Sci Rep. 2014;4:6462.View ArticlePubMedPubMed CentralGoogle Scholar
  42. Harima Y, Kageyama R. Oscillatory links of Fgf signaling and Hes7 in the segmentation clock. Curr Opin Genet Dev. 2013;23(4):484–90.View ArticlePubMedGoogle Scholar
  43. Hirata H, Bessho Y, Kokubu H, Masamizu Y, Yamada S, Lewis J, et al. Instability of Hes7 protein is crucial for the somite segmentation clock. Nat Genet. 2004;36(7):750–4.View ArticlePubMedGoogle Scholar
  44. Niwa Y, Masamizu Y, Liu T, Nakayama R, Deng CX, Kageyama R. The initiation and propagation of Hes7 oscillation are cooperatively regulated by Fgf and notch signaling in the somite segmentation clock. Dev Cell. 2007;13(2):298–304.View ArticlePubMedGoogle Scholar
  45. Sparrow DB, Guillen-Navarro E, Fatkin D, Dunwoodie SL. Mutation of Hairy-and-Enhancer-of-Split-7 in humans causes spondylocostal dysostosis. Hum Mol Genet. 2008;17(23):3761–6.View ArticlePubMedGoogle Scholar
  46. Sparrow DB, Sillence D, Wouters MA, Turnpenny PD, Dunwoodie SL. Two novel missense mutations in Hairy-and-Enhancer-of-Split-7 in a family with spondylocostal dysostosis. Eur J Hum Genet. 2010;18(6):674–9.View ArticlePubMedPubMed CentralGoogle Scholar
  47. Willet CE, Makara M, Reppas G, Tsoukalas G, Malik R, Haase B, et al. Canine disorder mirrors human disease: exonic deletion in HES7 causes autosomal recessive spondylocostal dysostosis in miniature Schnauzer dogs. PLoS One. 2015;10(2):e0117055.View ArticlePubMedPubMed CentralGoogle Scholar
  48. Filler S, Alhaddad H, Gandolfi B, Kurushima JD, Cortes A, Veit C, et al. Selkirk Rex: morphological and genetic characterization of a new cat breed. J Hered. 2012;103(5):727–33.View ArticlePubMedPubMed CentralGoogle Scholar
  49. Gandolfi B, Alhaddad H, Joslin SE, Khan R, Filler S, Brem G, et al. A splice variant in KRT71 is associated with curly coat phenotype of Selkirk rex cats. Sci Rep. 2013;3:2000. doi:10.1038/srep02000
  50. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual. 2nd ed. New York: Cold Spring Laboratory Press; 1989.Google Scholar
  51. Andrews KR, Luikart G. Recent novel approaches for population genomics data analysis. Mol Ecol. 2014;23(7):1661–7.View ArticlePubMedGoogle Scholar
  52. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. doi:10.1093/bioinformatics/btu170.View ArticlePubMedPubMed CentralGoogle Scholar
  53. FastQC: a quality control tool for high throughput sequence data [http://www.bioinformatics.bbsrc.ac.uk/projects/fastqc/]
  54. Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25(14):1754–60.View ArticlePubMedPubMed CentralGoogle Scholar
  55. Faust GG, Hall IM. SAMBLASTER: fast duplicate marking and structural variant read extraction. Bioinformatics. 2014;30(17):2503–5. doi:10.1093/bioinformatics/btu314.View ArticlePubMedPubMed CentralGoogle Scholar
  56. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.View ArticlePubMedPubMed CentralGoogle Scholar
  57. Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. arXiv. 2012.Google Scholar
  58. Sherry ST, Ward M-H, Kholodov M, Baker J, Phan L, Smigielski EM, et al. dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001;29(1):308–11.View ArticlePubMedPubMed CentralGoogle Scholar
  59. Cingolani P, Platts A, Wang LL, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. 2012;6(2):80–92.View ArticlePubMedPubMed CentralGoogle Scholar
  60. Cunningham F, Amode MR, Barrell D, Beal K, Billis K, Brent S, et al. Ensembl 2015. Nucleic Acids Res. 2015;43(D1):D662–9.View ArticlePubMedPubMed CentralGoogle Scholar
  61. Rosenbloom KR, Armstrong J, Barber GP, Casper J, Clawson H, Diekhans M, et al. The UCSC Genome Browser database: 2015 update. Nucleic Acids Res. 2015;43(D1):D670–81.View ArticlePubMedPubMed CentralGoogle Scholar
  62. Lyons LA. DNA mutations of the cat: the good, the bad and the ugly. J Feline Med Surg. 2015;17(3):203–19.View ArticlePubMedGoogle Scholar
  63. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8.View ArticlePubMedPubMed CentralGoogle Scholar
  64. Gandolfi B, Daniel RJ, O’Brien DP, Guo LT, Youngs MD, Leach SB, et al. A novel mutation in CLCN1 associated with feline myotonia congenita. PLoS One. 2014;9(10):e109926.View ArticlePubMedPubMed CentralGoogle Scholar
  65. Gandolfi B, Gruffydd-Jones TJ, Malik R, Cortes A, Jones BR, Helps CR, et al. First WNK4-hypokalemia animal model identified by genome-wide association in Burmese cats. PLoS One. 2012;7(12):e53173.View ArticlePubMedPubMed CentralGoogle Scholar
  66. Gandolfi B, Alhaddad H, Affolter VK, Brockman J, Haggstrom J, Joslin SE, et al. To the root of the curl: a signature of a recent selective sweep identifies a mutation that defines the Cornish rex cat breed. PLoS One. 2013;8(6):e67105.View ArticlePubMedPubMed CentralGoogle Scholar
  67. Toonen RJ, Hughes S. Increased throughput for fragment analysis on an ABI PRISM 377 automated sequencer using a membrane comb and STRand software. Bio Techniques. 2001;31(6):1320–4.Google Scholar
  68. Daetwyler HD, Capitan A, Pausch H, Stothard P, Van Binsbergen R, Brøndum RF, et al. Whole-genome sequencing of 234 bulls facilitates mapping of monogenic and complex traits in cattle. Nature Genet. 2014;46(8):858–65.View ArticlePubMedGoogle Scholar
  69. Gilliam D, Kolicheski A, Johnson GS, Mhlanga-Mutangadura T, Taylor JF, Schnabel RD, et al. Golden Retriever dogs with neuronal ceroid lipofuscinosis have a two-base-pair deletion and frameshift in CLN5. Mol Genet Metab. 2015;115(2–3):101–9.View ArticlePubMedGoogle Scholar
  70. Guo J, Johnson GS, Brown HA, Provencher ML, da Costa RC, Mhlanga-Mutangadura T, et al. A CLN8 nonsense mutation in the whole genome sequence of a mixed breed dog with neuronal ceroid lipofuscinosis and Australian Shepherd ancestry. Mol Genet Metab. 2014;112(4):302–9.View ArticlePubMedGoogle Scholar
  71. Guo J, O’Brien DP, Mhlanga-Mutangadura T, Olby NJ, Taylor JF, Schnabel RD, et al. A rare homozygous MFSD8 single-base-pair deletion and frameshift in the whole genome sequence of a Chinese Crested dog with neuronal ceroid lipofuscinosis. BMC Vet Res. 2014;10:960.View ArticlePubMed CentralGoogle Scholar
  72. Frischknecht M, Jagannathan V, Leeb T. Whole genome sequencing confirms KIT insertions in a white cat. Anim Genet. 2014;46(1):98.View ArticlePubMedGoogle Scholar
  73. Gandolfi B, Grahn RA, Creighton EK, Williams DC, Dickinson PJ, Sturges BK, et al. COLQ variant associated with Devon Rex and Sphynx feline hereditary myopathy. Anim Genet. 2015;46(6):711–5.View ArticlePubMedGoogle Scholar
  74. Mullikin JC, Hansen NF, Shen L, Ebling H, Donahue WF, Tao W, et al. Light whole genome sequence for SNP discovery across domestic cat breeds. BMC Genomics. 2010;11(1):406.View ArticlePubMedPubMed CentralGoogle Scholar
  75. Pontius J, Mullikin J, Smith D, Lindblad-Toh K, Gnerre S, Clamp M, et al. Agencourt Sequencing Team; NISC Comparative Sequencing Program (2007) Initial sequence and comparative analysis of the cat genome. Genome Res. 2007;17(11):1675–89.View ArticlePubMedPubMed CentralGoogle Scholar
  76. Ajay SS, Parker SC, Abaan HO, Fajardo KV, Margulies EH. Accurate and comprehensive sequencing of personal genomes. Genome Res. 2011;21(9):1498–505.View ArticlePubMedPubMed CentralGoogle Scholar
  77. CFA: Registration Statistics 2010. In: Cat Fanciers’ Association Almanac. vol. 2011. Trenton, New Jersey: Cat Fanciers’ Association, Inc.; 2010.Google Scholar
  78. Steinberg RH, Reid M, Lacy PL. The distribution of rods and cones in the retina of the cat (Felis domesticus). J Comp Neurol. 1973;148(2):229–48.View ArticlePubMedGoogle Scholar
  79. Menotti-Raymond M, David VA, Schäffer AA, Stephens R, Wells D, Kumar-Singh R, et al. Mutation in CEP290 discovered for cat model of human retinal degeneration. J Hered. 2007;98(3):211–20.View ArticlePubMedGoogle Scholar
  80. den Hollander AI, Koenekoop RK, Yzer S, Lopez I, Arends ML, Voesenek KE, et al. Mutations in the CEP290 (NPHP6) gene are a frequent cause of Leber congenital amaurosis. Am J Hum Genet. 2006;79(3):556–61.View ArticleGoogle Scholar
  81. Menotti-Raymond M, Deckman KH, David V, Myrkalo J, O’Brien SJ, Narfström K. Mutation discovered in a feline model of human congenital retinal blinding disease. Invest Ophthalmol Vis Sci. 2010;51(6):2852–9.View ArticlePubMedPubMed CentralGoogle Scholar

Copyright

© Lyons et al. 2016