Prevalence of PKD1 gene mutation in cats in Turkey and pathogenesis of feline polycystic kidney disease

Abstract

Polycystic kidney disease (PKD) is one of the most common hereditary diseases in cats, with high prevalence in Persian and Persian-related cats. PKD is caused mainly by an inherited autosomal dominant (AD) mutation, and animals may be asymptomatic for years. We screened 16 cats from various breeds exhibiting a renal abnormality by ultrasound examination and genotyped them for the c.10063C>A transversion on exon 29 of the polycystin-1 (PKD1) gene, by PCR–restriction fragment length polymorphism (PCR-RFLP). Among these cats, a Siamese nuclear family of 4 cats with ancestral hereditary renal failure were screened by whole-genome sequencing (WGS) to determine novel variations in genes associated with both AD and autosomal recessive PKD in humans. During the study period, one cat died as a result of renal failure and was forwarded for autopsy. Additionally, we screened 294 cats asymptomatic for renal disease (Angora, Van, Persian, Siamese, Scottish Fold, Exotic Shorthair, British Shorthair, and mixed breeds) to determine the prevalence of the mutation in cats in Turkey. Ten of the symptomatic and 2 of the asymptomatic cats carried the heterozygous C → A transversion, indicating a prevalence of 62.5% and 0.68%, respectively. In the WGS analysis of 4 cats in the Siamese nuclear family, novel variations were determined in the fibrocystin gene (PKHD1), which was not compatible with dominant inheritance of PKD.

Keywords

ADPKD autosomal dominant cats hereditary disease kidney polycystic kidney disease

Introduction

Domestic cats are potential models for 218 genetic diseases and phenotypes, of which 72 are known to have causative mutations shared by cats and humans, including polycystic kidney disease (PKD; www.omia.org). PKD is characterized by uni- or bilateral renal cyst formation, and is a systemic, progressive hereditary disease with clinical signs that may develop at any age; cysts can also form in ductal organs, such as liver, seminal vesicle, pancreas, and arachnoid membrane.^19,27 In humans, PKD has genetic heterogeneity, with autosomal dominant (AD) and autosomal recessive (AR) inheritance patterns in various genes. Polycystin-1 (PKD1), polycystin-2 (PKD2), neutral alpha-glucosidase AB (GANAB, PKD3),²⁰ and DnaJ homolog subfamily B member 11(DNAJB11, PKD6),¹² are associated with ADPKD; fibrocystin (PKHD1, PKD4) and zinc finger protein DZIP1L (DZIP1L, PKD5)²⁰ are associated with the AR version of the disease (ARPKD). Products of these genes are involved in the formation or structure of cilia, which are microtubule-based organelles that function in fluid transport and chemo- and mechano-sensation.²⁰ Mutations in the PKD1 gene result in one of the AD forms of PKD, characterized by cyst formation in the kidney in utero or at later ages. Over 1,270 mutations causing end-stage renal disease have been determined in the PKD1 gene in humans⁹; only 1 causative mutation has been described in cats.²³

In cats, PKD was first described in a Persian cat family with renal cysts⁵ and an AD inheritance pattern^6,14 mediated by a point mutation of C → A at sequence position c.10063 in exon 29 of the PKD1 gene.²³ This mutation results in a stop codon (opal stop codon, UGA) at amino acid position 3,284 of the protein. Since then, several studies have reported high disease prevalence in mainly Persian and Persian-related cat breeds: 49% in the United Kingdom, 43% in Germany and Australia, 42% in France, 36% in Slovenia, 41% in Italy, and 40% in Japan.^{1
-3,7,8,13,26}

The aims of our study were to determine whether the PKD1 mutation was the cause of renal failure in cats referred to veterinary clinics in Turkey, and to measure the prevalence of this mutation in various breeds of cats in Turkey. We screened for the mutation in symptomatic and asymptomatic cats by PCR–restriction fragment length polymorphism (PCR-RFLP). Additionally, we performed whole-genome sequencing (WGS) screening of 4 cats, focusing on genes associated with PKD in humans.

Materials and methods

Sample collection

We studied 310 domestic cats aged 2 mo to 6 y old, from Turkey. Buccal swab samples were collected for genetic analysis from various cats of various breeds: 35 Angora, 56 Van, 26 Siamese, 22 Scottish Fold, 40 Persian, 9 Exotic Shorthair, 5 British Shorthair, and 117 mixed-breed cats (Ankara University Local Ethical Board of Animal Experiments, Approvals 2017-18-149 and 2018-4-33). Ultrasonography was performed (linear probe 7.5 MHz, convex probe 3.5 MHz, AU5 Harmonic; Esaote Biomedica).

Pathology examinations

A postmortem examination was performed on a Persian cat that died as a result of renal failure. Kidney and liver sections were fixed in 4% paraformaldehyde solution for 48 h and processed routinely for histologic (hematoxylin and eosin, and Masson trichrome stains) and immunohistochemical examinations.²¹ Cysts were photographed, and their dimensions were measured by light microscopy (DFC450C digital color microscope camera; Leica).

Immunohistochemistry

Immunohistochemical assays of kidney and liver tissues were performed (UltraVision detection system; Thermo Fisher Scientific), according to the manufacturer’s directions, to demonstrate expression of caspase-3, B-cell lymphoma 2 (Bcl-2), transforming growth factor–beta (TGF-β), and tumor necrosis factor–alpha (TNF-α). The primary antibodies used were: anti–caspase-3 (Ready to use; NeoMarkers), to show the activated apoptotic pathway; anti–Bcl-2 (sc-7382, dilution 1:50; Santa Cruz Biotechnology) for the evaluation of apoptosis inhibition; anti–TGF-β (sc-130348, dilution 1:100; Santa Cruz Biotechnology) to evaluate the initiation of fibrosis; and anti–TNF-α (sc-52746, dilution 1:100; Santa Cruz Biotechnology) to show the degree of inflammation.

In brief, tissue sections were deparaffinized in xylene and hydrated through graded alcohols. To unmask antigens, tissue sections were placed in a microwave oven in citrate buffer, for 20 min at the high setting. Endogenous peroxidase activity was inhibited by 0.1% H₂O₂ in methanol for 10 min, and nonspecific labeling was blocked by pre-incubation with normal goat serum for 10 min. Thereafter, the sections were incubated with anti–caspase-3, anti–Bcl-2, anti–TGF-β, and anti–TNF-α for 60 min in a humidity chamber at room temperature. Sections were then incubated with biotinylated secondary antibody for 15 min, labeled with horseradish peroxidase for 15 min, and 3-amino-9-ethylcarbazole (AEC) chromogen substrate solution (Thermo Fisher Scientific) for ~ 10 min. Sections were counterstained with Mayer hematoxylin for 1 min and coverslipped with a water-based mounting medium (Thermo Fisher Scientific). The density of immunopositive staining was photographed by light microscopy (BX51 microscope with DP25 digital camera; Olympus).

Molecular genetic analysis

DNA was extracted from the buccal swabs (Quick-DNA universal kit; Zymo), following the manufacturer’s instruction. Briefly, swabs were shredded and placed into centrifuge tubes; 200 µL of genomic lysis buffer was added to the samples and let stand at room temperature for 10 min. The mixture was then transferred to a spin column; DNA was eluted with 50 µL of elution buffer. Obtained DNA was quantified spectrophotometrically (NanoDrop; Thermo Fisher Scientific) and visualized for integrity on dyed (SafeView Classic; ABM) 1% agarose gel; electrophoresis under 100 V. Quantified DNA was stored at –20°C until further analyzed.

The 559-bp region harboring the causative mutation c.10063C>A was amplified by using specific oligonucleotides forward (5′-CAGGTAGACGGGATAGACGA-3′) and reverse (5′-TTCTTCCTGGTCAACGACTG-3′).²³ PCR was set up with 1 pmol of each oligonucleotide, 200 µM dNTP, 1× buffer, 1.5 µM MgCl₂, 2 IU of Taq polymerase (MBI Fermentas), and ~20 ng of genomic DNA (gDNA) added as the template in a 25-µL reaction. PCR conditions were initial denaturation at 94°C for 4 min, followed by 30 cycles of denaturation at 94°C for 30 s, annealing at 58°C for 30 s, elongation at 72°C for 30 s, and 10 min at 72°C for final elongation. Obtained PCR products were used for direct mutation analysis by restriction enzyme digestion with MLY1 (NEB) as described previously.²³ Enzyme digestion was done at 37°C for 45 min and enzyme deactivation at 67°C for 20 min. Products were visualized on dyed (SafeView Classic) 2% agarose electrophoresis gels. Determined mutations were also verified by bidirectional DNA sequencing analysis (ABI 310; Applied Biosystems) using the same oligonucleotides used in the PCR. Sequence alignment to the GenBank reference sequence AY612847 was performed using commercial software (BioEdit, https://bioedit.software.informer.com/7.2/).¹⁷

WGS analyses of the Siamese cat family were performed (HiSeq 2500; Illumina) by a private service provider (BGI). Briefly, gDNA of the 4 individuals was fragmented randomly in separate sample tubes. Sequencing libraries were constructed by adapter ligation (Nextera DNA flex library prep; Illumina), with amplification by PCR and gel purification of the expected bands. DNA cluster preparation was performed and sequenced by paired-end sequencing. Whole genomes were resequenced with an average of 20 coverage. Obtained fastq files were aligned to the latest Felis catus genome (FelCat9, University of California–Davis, Santa Cruz). Sequence variants on 6 genes causing ADPKD and ARPKD in humans [PKD1, PKD2, GANAB (PKD3), PKHD1 (PKD4), DZIP1L (PKD5), and DNAJB11 (PKD6)] were extracted by using the SnpSift tool (http://snpeff.sourceforge.net/SnpSift.html).¹⁰ To check the novelty of the obtained variants, Variant Effect Predictor tool²⁴ was used in the Ensembl database to provide the published cat genomes.

Results

Clinical examinations and ultrasonography

Among 310 cats, 16 had various signs including fatigue, anorexia, and vomiting, and were referred to veterinary clinics. Clinical examination by palpation revealed an increase in total kidney volume; ultrasonography was performed. Anechoic and/or cystic lesions were seen in the cortex of both kidneys in 12 cats; bilateral renomegaly was seen in 4 cats. Cystic lesions were not evident in any of the other parenchymal organs.

PCR-RFLP and WGS analysis

Extracted DNA amounts varied 20–150 ng, and PCR was used to amplify regions of interest from all of the samples. Ten of the 12 cats with cystic kidneys were found heterozygous for the PKD1 mutation by PCR-RFLP, with DNA sequencing confirming a C→A transversion. After the presence of the mutation had been confirmed, owners were informed of their cat’s disease status. To limit spread of the mutant allele in the cat population, sterilization was suggested. Interestingly, 6 of the 16 symptomatic cats (2 mixed-breed cats with bilateral kidney cysts, and 4 Siamese cats having renomegaly without cysts) were negative for C→A, and sequencing of the amplified region did not show any nucleotide base differentiation. None of the cats of the native breeds of Anatolia carried the mutation (35 Angora and 56 Van), nor did 9 Exotic Shorthair or 5 British Shorthair cats (Table 1).

Table 1.

Autosomal dominant polycystic kidney disease–related PKD1 gene mutation (c.10063C>A) in the screened cat population. Homozygous wild and heterozygous genotypes are presented as CC and CA, respectively.

Breed	n	Asymptomatic (n = 294)		Symptomatic (n = 16)		Percentage of mutation
Breed	n	CC	CA	CC	CA	Percentage of mutation
Angora	35	35	0	0	0	0
Van	56	56	0	0	0	0
Siamese	26	20	0	4	2	7.7
Scottish Fold	22	21	1	0	0	4.5
Persian	40	33	1	0	6	17.5
Exotic Shorthair	9	9	0	0	0	0
British Shorthair	5	5	0	0	0	0
Mixed	117	113	0	2	2	1.7
Total	310	292	2	6	10	3.9

WGS analysis of the causative dominant and recessive genes in the Siamese cat family revealed several novel variations in all genes as well as missense homozygous and heterozygous mutations in PKD2, DZIP1L, and PKHD1 genes. Missense homozygous variations in the PKHD1 gene had significant SIFT scores (Sorting Intolerant From Tolerant; SIFT ≤ 0.03²⁴; Table 2).

Table 2.

Novel missense homozygous and heterozygous variations in PKD2, DZIP1L, and PKHD1 genes, and SIFT scores of the variants.

Location	Allele	Gene symbol	cDNA position	Amino acids	SIFT score
chrB1:134986222	G	PKD2	2943	I/L	Tolerated low confidence (1)
chrB2:52046775	T	PKHD1	11772	V/I	Deleterious (0)
chrB2:52152625	C	PKHD1	9663	C/G	Tolerated (0.48)
chrB2:52152625	C	PKHD1	9424	C/G	Tolerated (0.46)
chrB2:52412222	T	PKHD1	5671	S/Y	Deleterious (0.03)
chrB2:52428830	T	PKHD1	3592	S/N	Tolerated (0.73)
chrB2:52458753	G	PKHD1	2272	E/A	Tolerated (1)
chrB2:52488628	C	PKHD1	365	Q/R	Tolerated (0.52)
chrB2:52492657	G	PKHD1	375	Y/H	Deleterious (0.03)
chrC2:127474195	A	DZIP1L	2025	V/M	Tolerated (0.07)

Histopathology and immunohistochemical analyses

Autopsy of the Persian cat that died in renal failure revealed multiple macroscopic cysts in the cortex of both kidneys (Fig. 1A) and liver (Fig. 1B). Histologic examination of the kidney revealed multiple cystic structures containing homogeneous pink material. Adjacent to cysts, normal, degenerate, cystic, and atrophic tubules and glomeruli were observed. In some tubules, severe epithelial degeneration and necrosis were observed, as well as debris and residues of 2 or 3 tubules, along with microcystic areas (Fig. 1C). In some of the proximal tubules, thyroidization and epithelial hyperplasia were observed, indicating regeneration. Interstitial fibrosis was common in the cystic-degenerate areas. The areas of fibrosis were confirmed with Masson trichrome staining (Fig. 1D).

Figure 1.

Polycystic kidney and liver in a Persian cat. A. Gross appearance of cysts in the cortex and corticomedullary border of the kidney (arrows). Bar = 0.5 cm. B. Two subserosal macroscopic cysts (arrows) in left lobe of the liver. Bar = 1 cm. C. Homogeneous eosinophilic material in the cyst lumen (asterisk) and widespread fibrosis (arrows). H&E. Bar = 150 µm. D. Cystic and dilated tubules (asterisk) and fibrosis (arrows) in the kidney. Masson trichrome. Bar = 50 µm.

Immunohistochemical analyses revealed Bcl-2 (Fig. 2A); degenerate tubule epithelium showed strong caspase-3 immunoreactivity (Fig. 2B), and TNF-α (Fig. 2C) immunoreactivity in degenerate tubule epithelium. TNF-α and caspase-3 immunopositive reactions had similar localization. TGF-β activity was detected in fibrotic areas and tubule epithelium (Fig. 2D).

Figure 2.

Immunoperoxidase staining of the polycystic kidney in a Persian cat. Indirect avidin–biotin complex immunoperoxidase test, AEC chromogen with Mayer hematoxylin counterstaining. Bars = 100 µm. A. Mild Bcl-2 immunopositivity (faded red color). B. Strong caspase-3 expressions (red staining) in tubular epithelia (arrows). C. Pronounced TNF-α expression in degenerate tubular epithelia (arrows). D. Strong TGF-β expression in the degenerate tubules (arrows) surrounding a cystic structure.

Discussion

We found PKD prevalence of 17.5%, 2.7%, and 3.9% in Persian, Persian-related cats (Scottish Fold, Exotic Shorthair, and British Shorthair), and the total cat sample, respectively. In previous studies, which were mainly of Persian and Persian-related cat breeds, the prevalence of PKD was 26% (lowest) to 49% (highest).^1,2,8,19 In our study, the mutation was not detected in Exotic and British Shorthair breeds, which might be the result of the small number of subjects. Moreover, molecular examination showed that none of the cats from Angora and Van breeds had the PKD1 mutation. Although Angora cats are known to be prone to chronic kidney disease (pers. comm., clinicians’ view), we suggest that it is not caused by the known PKD1 mutation. We conclude that native Anatolian cat breeds can be considered a low-risk group.

Previous studies also reported that several cats had renal cysts visible on ultrasonographic examination, but that their 29th exon of PKD1 was free of the ADPKD1 mutation,^16,19 suggesting that the disease has a multigenic background in cats. Mutations that cause ciliopathies (e.g., HNF1ß) or diseases such as nephronophthisis or Bardet–Biedl–Joubert syndrome can all mimic PKD and should be considered phenocopies.⁴ Several gene mutations, such as PKD1 and PKD3, affecting ciliary structures can result in different severities of disease, suggesting epistatic gene interaction.¹¹ Overall, there is a lack of studies on the genes associated with cilia formation in cats. Moreover, environmental factors such as industrial pollutants in the water, or chemicals or drugs such as diphenylamine, diphenylthiazole, and streptozotocin,¹⁵ can cause renal cysts, and should be kept in mind when evaluating cases of renal failure.⁴

Uncontrolled mating of cats harboring the mutation caused the spread of the mutant allele in the cat population.^22,23 Persian cat ancestry is a known major factor in the spread of the PKD mutation,^2,8,19 but other breeds, such as Siamese, should also be considered at high risk. These individuals should be screened before accepting them into the breeding programs. According to our WGS analysis, one of the Siamese cats had homozygous non-synonymous mutations in the PKHD1 gene (Y46H, S1811Y, and V3845I). These mutations should be further studied for their impact on protein function. Although this cat had no conclusive clinical findings, close monitoring of kidney function was recommended to the owner.

In humans with PKD, ~ 300 mutations have been identified in the PKHD1 gene, whereas in cats only 5 exonic variations have been reported in Ensembl, of which one missense change has been identified.¹⁸ Although associated with a severe PKD phenotype in humans, the PKHD1 gene has not been screened previously in cats, to our knowledge.

In PKD, the balance between tubular degeneration and activation of necrosis and apoptosis is an important factor in the occurrence of cysts. The affected polycystin-1 protein in PKD1-positive cats has an important role in proliferation and cell differentiation of tubule epithelium, and it also has apoptosis-inhibiting activity on the tubular epithelium.^23,25 From this perspective, a high incidence of caspase-3 expression, especially in tubule epithelia, may suggest that apoptosis is actively involved in the development of cysts in the kidney. Accordingly, Bcl-2 expression, an anti-apoptotic protein, in epithelia surrounding the microscopic cysts was negative, although morphologically normal tubule epithelia were immunopositive.

Progressive tubular degeneration and microscopic cysts seen in the immediate vicinity of the macroscopic cysts suggests that the pathogenesis of this genetic disease might be related to induction of programmed cell death in affected cells. In addition to cystic structures in the kidney tissue, fibrosis and interstitial tissue enlargement have been noticed, and increased TGF-β expression around these fibrotic areas is noteworthy. This suggests that, in animals with PKD, renal failure may not be only caused by cyst formation, but renal fibrosis might also be an important factor. There are numerous studies on the role of TGF-β in renal fibrosis and progressive renal failure.²⁸ Fibrosis might be noted as an important factor in renal failure in PKD, and apoptotic markers might be particularly useful in illumination of etiopathogenesis.

Genes associated with cilia formation and those causing ADPKD and ARPKD should be screened for mutations in cat populations to explore differential cyst formation and disease severity in PKD. Furthermore, WGS, clinical analysis, and histopathology should be considered in tandem to discover the candidate genes in cats.

Footnotes

Acknowledgements

We thank the veterinary practitioners who provided materials for our study. We also thank M. Krzywinski for editing the figures.

Declaration of conflicting interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The authors received no financial support for the research, authorship, and/or publication of this article.

ORCID iDs

Nüket Bilgen

Bengi Çınar Kul

Sevim Isparta

Özge Ş. Çıldır

Zeynep Bars

References

Barrs

, et al. Prevalence of autosomal dominant polycystic kidney disease in Persian cats and related breeds in Sydney and Brisbane. Aust Vet J 2001;79:257–259.

Barthez

, et al. Prevalence of polycystic kidney disease in Persian and Persian related cats in France. J Feline Med Surg 2016;5:345–347.

Beck

Lavelle

RB.

Feline polycystic kidney disease in Persian and other cats: a prospective study using ultrasonography. Aust Vet J 2001;79:181–184.

Bergmann

ARPKD and early manifestations of ADPKD: the original polycystic kidney disease and phenocopies. Pediatr Nephrol 2015;30:15–30.

Biller

, et al. Polycystic kidney disease in a family of Persian cats. J Am Vet Med Assoc 1990;196:1288–1290.

Biller

, et al. Inheritance of polycystic kidney disease in Persian cats. J Hered 1996;87:1–5.

Bonazzi

, et al. Comparison between ultrasound and genetic testing for the early diagnosis of polycystic kidney disease in Persian and Exotic Shorthair cats. J Feline Med Surg 2009;11:430–434.

Cannon

, et al. Prevalence of polycystic kidney disease in Persian cats in the United Kingdom. Vet Rec 2001;149:409–411.

Chebib

Torres

VE.

Autosomal dominant polycystic kidney disease: core curriculum 2016. Am J Kidney Dis 2016;67:792–810.

10.

Cingolani

, 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 (Austin) 2012;6:80–92.

11.

Cordido

, et al. The genetic and cellular basis of autosomal dominant polycystic kidney disease—a primer for clinicians. Front Pediatr 2017;5:279.

12.

Cornec-Le Gall

, et al. Monoallelic mutations to DNAJB11 cause atypical autosomal-dominant polycystic kidney disease. Am J Hum Genet 2018;102:832–844.

13.

Domanjko-Petric

, et al. Polycystic kidney disease: a review and occurrence in Slovenia with comparison between ultrasound and genetic testing. J Feline Med Surg 2008;10:115–119.

14.

Eaton

, et al. Autosomal dominant polycystic kidney disease in Persian and Persian-cross cats. Vet Pathol 1997;34:117–126.

15.

Gardner

Jr.

Cystic kidneys. Kidney Int 1988;33:610–621.

16.

Guerra

, et al. Congenital hepatic fibrosis and polycystic kidney disease not linked to C>A mutation in exon 29 of PKD1 in a Persian cat. JFMS Open Rep 2015;1:2055116915619191.

17.

Hall

TA.

BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 1999;41:95–98.

18.

Hunt

, et al. Ensembl variation resources. Database (Oxford) 2018;2018.

19.

Lee

, et al. Diagnosis of feline polycystic kidney disease by a combination of ultrasonographic examination and PKD1 gene analysis. Vet Rec 2010;167:614–618.

20.

, et al. Mutations in DZIP1L, which encodes a ciliary-transition-zone protein, cause autosomal recessive polycystic kidney disease. Nat Genet 2017;49:1025–1034.

21.

Luna

, ed. Manual of Histologic Staining Methods of Armed Forces Institute of Pathology. McGraw-Hill, 1968.

22.

Lyons

LA.

DNA mutations of the cat: the good, the bad and the ugly. J Feline Med Surg 2015;17:203–219.

23.

Lyons

, et al. Feline polycystic kidney disease mutation identified in PKD1. J Am Soc Nephrol 2004;15:2548–2555.

24.

McLaren

, et al. The Ensembl variant effect predictor. Genome Biol 2016;17:122.

25.

Nivy

, et al. Polycystic kidney disease in four British shorthair cats with successful treatment of bacterial cyst infection. J Small Anim Pract 2015;56:585–589.

26.

Sato

, et al. Epidemiological evaluation of cats associated with feline polycystic kidney disease caused by the feline PKD1 genetic mutation in Japan. J Vet Med Sci 2019;81:1006–1011.

27.

Sweeney

Avner

. Polycystic kidney disease, autosomal dominant. 2019 Feb 14. GeneReviews. University of Washington. Accessed 2020 May 27: https://www.ncbi.nlm.nih.gov/books/NBK1326/

28.

Yang

, et al. Caspase-3 and apoptosis in experimental chronic renal scarring. Kidney Int 2001;60:1765–1776.