Next Article in Journal
Overall Proportion of Total DNA Consistent with an Individual Briefly Handling a Firearm
Previous Article in Journal
Analysis of Nucleotide Variations in Human G-Quadruplex Forming Regions Associated with Disease States
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 12
10.3390/genes14122126
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single Nucleotide Polymorphisms Associated with AA-Amyloidosis in Siamese and Oriental Shorthair Cats

by
Stella L. Esders
1,
Kirsten Hülskötter
2,
Tom Schreiner
2,
Peter Wohlsein
2,
Jessica Schmitz
3,
Jan H. Bräsen
3 and
Ottmar Distl
1,*
1
Institute for Animal Breeding and Genetics, University of Veterinary Medicine Hannover (Foundation), 30559 Hannover, Germany
2
Department of Pathology, University of Veterinary Medicine Hannover (Foundation), 30559 Hannover, Germany
3
Nephropathology Unit, Institute of Pathology, Hannover Medical School, 30625 Hannover, Germany
*
Author to whom correspondence should be addressed.
Genes 2023, 14(12), 2126; https://doi.org/10.3390/genes14122126
Submission received: 13 September 2023 / Revised: 21 November 2023 / Accepted: 23 November 2023 / Published: 25 November 2023
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figure
Versions Notes

Abstract

:
AA-amyloidosis in Siamese and Oriental shorthair cats is a lethal condition in which amyloid deposits accumulate systemically, especially in the liver and the thyroid gland. The age at death of affected cats varies between one and seven years. A previous study indicated a complex mode of inheritance involving a major locus. In the present study, we performed a multi-locus genome-wide association study (GWAS) using five methods (mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB and ISIS EM-BLASSO) to identify variants associated with AA-amyloidosis in Siamese/Oriental cats. We genotyped 20 affected mixed Siamese/Oriental cats from a cattery and 48 healthy controls from the same breeds using the Illumina Infinium Feline 63 K iSelect DNA array. The multi-locus GWAS revealed eight significantly associated single nucleotide polymorphisms (SNPs) on FCA A1, D1, D2 and D3. The genomic regions harboring these SNPs contain 55 genes, of which 3 are associated with amyloidosis in humans or mice. One of these genes is SAA1, which encodes for a member of the Serum Amyloid A family, the precursor protein of Amyloid A, and a mutation in the promotor of this gene causes hereditary AA-amyloidosis in humans. These results provide novel knowledge regarding the complex genetic background of hereditary AA-amyloidosis in Siamese/Oriental cats and, therefore, contribute to future genomic studies of this disease in cats.

1. Introduction

Amyloidosis is a condition characterized by the deposition of amyloid, an insoluble, proteinaceous substance, that accumulates in various organs and tissues [1]. These deposits cause an impairment of organ architecture via pressure atrophy of adjacent cells, thus potentially causing loss of function of affected organs [1]. Amyloid possesses a fibrillary architecture and a ß-pleated sheet conformation, which is responsible for the affinity for Congo red stain and green birefringence in polarization microscopy [2]. In cats, a systemic form of amyloidosis is Amyloid A amyloidosis (AA-amyloidosis), with acute phase reactant Serum Amyloid A (SAA) as its precursor protein [3,4]. The predisposed breeds are Abyssinian and Somali along with Siamese and Oriental, in which hereditary forms of AA-amyloidosis are recognized [5,6,7,8]. In Abyssinian/Somali cats, the continuous failure of the kidneys until the age of five years due to a high degree of medullar and glomerular amyloid deposits is the main cause of death of affected cats [5,9,10]. In Siamese/Oriental cats, the organs mostly affected by the deposition of amyloid are the liver and the thyroid gland, but variable other organs and tissues can be affected as well [11]. The main cause of death of affected cats is rupture of the liver with subsequent blood loss in the abdominal cavity and hypovolemic shock [6,7,8,12]. The disease is lethal with a variable age at death between one and seven years [6]. Hereditary AA-amyloidosis in Oriental cats is autosomal, and neither the onset nor progression of the disease is influenced by the sex [11]. A segregation analysis of 17 clinically affected cats revealed a complex mode of inheritance, composed of a polygenic component with a major locus with Mendelian inheritance and heritability estimates based on the pedigree data of 0.56 ± 0.09 [11].
A study comparing Amyloid A protein sequences of one affected Abyssinian cat and two affected Siamese cats concluded that sequence differences may contribute to the different tissue affinities in these breeds, but this conclusion was not confirmed in past studies [13]. A recent study of two affected Abyssinian cats and 195 controls from different breeds concluded a polygenic background of AA-amyloidosis in this breed and found variants in Amyloid Beta Precursor Protein gene (APP) and Transthyretin (TTR) in affected cats [14]. Mutations in APP and TTR are associated with amyloidosis in humans [15,16,17,18].
The objectives of the present study were to identify variants associated with hereditary AA-amyloidosis in mixed Siamese/Oriental cats based on Illumina Infinium Feline 63 K iSelect DNA array data of cases and controls with multi-locus GWAS models. Within the genomic regions of these variants, we searched for putative functional candidate genes that have previously been reported to be associated with hereditary forms of amyloidosis in humans and mice [15,18,19,20].

2. Materials and Methods

2.1. Ethical Statement

All animal work was conducted following the national and international guidelines for animal welfare following the guidelines of the Declaration of Helsinki and was approved by the Institutional Review Board of the University of Veterinary Medicine Hannover (Foundation) and the responsible state veterinary office of Lower Saxony (registration number 33.9-42502-05-18A284). The sampling and handling of the animals followed European Union guidelines for animal care and handling and the Guidelines of Good Veterinary Practices. Cats were euthanized based on a thorough veterinary diagnosis and then submitted to the Department of Pathology, University of Veterinary Medicine Hannover (Foundation), by cat owners. EDTA-blood, hair root and buccal swab samples of controls were obtained from the biobank of the Institute for Animal Breeding and Genetics, University of Veterinary Medicine Hannover, and were taken for diagnostic purposes by accredited veterinarians.

2.2. Animals

In the present study, we included 68 cats mainly collected from a German cattery consisting of mixes between Siamese and Oriental shorthair cats previously analyzed in the complex segregation analysis [11] (Supplementary Figure S1 and Table S1).
Cases include cats that were affected by clinical signs consistent with systemic AA-amyloidosis as rupture of the liver and died between the ages of one and seven years (n = 20). All cases belong to the aforementioned cattery. Histopathological confirmation of the diagnosis was available for 16/20 cats. Cats from the same breeds that never showed clinical signs consistent with systemic AA-amyloidosis and were at least ten years old at the time of the analyses were used as controls (n = 48). Relationships of the cats are shown in Supplementary Figure S1. Entry numbers 67 and 68 are not shown in Supplementary Figure S1 as they are not closely related to the other cats (1–66) in this study.

2.3. Histological and Immunohistochemical Examination

Histological and immunohistochemical examinations were performed as previously described [11]. Briefly, the tissue of the 16 autopsied cats was fixed, stained with Congo red and spectated via polarization light microscopy to confirm amyloid deposits. The localization of amyloid deposits in different organs and tissues was evaluated (Supplementary Table S2).
Additionally, in 9/16 cases, distributed across the different generations with affected individuals, amyloid classification was performed. In these nine animals, Amyloid A was classified with immunostaining, which was conducted with an automated platform (Ventana ULTRA; Ventana Medical Systems) with mouse monoclonal antibodies [11].

2.4. Genotyping on the Illumina Infinium Feline 63 K iSelect DNA Array

Genomic DNA was isolated from EDTA-blood samples, buccal swabs and hair roots of the 48 controls, and tissue samples of the 20 cases using phenol/chloroform extraction. DNA concentration per sample was adjusted to 50 ng/µL with Nanodrop ND-1000 (Peqlab Biotechnology, Erlangen, Germany). Genotyping for the 68 cats was performed using the Illumina Infinium Feline 63 K iSelect DNA array with 62,897 SNPs (Neogen, Lincoln, NE, USA). Chromosomal positions of the SNPs were identified using the current reference genome from the National Center for Biotechnology Information (NCBI) F. catus_Fca126_mat1.0 (https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_018350175.1/, accessed on 4 September 2023) [21] and BLAST [22].

2.5. Genome-Wide Association Studies

The genome-wide association analyses were performed using the R package mrMLM v4.0.2 [23], which implements five methods: the multi-locus random-SNP-effect mixed linear model (mrMLM) [24], the fast multi-locus random-SNP-effect mixed linear model (FASTmrMLM) [25], the fast multi-locus random-SNP-effect efficient mixed model association (FASTmrEMMA) [26], the polygenic-background-control-based least angle regression plus empirical Bayes (pLARmEB) [27] and the iterative sure independence screening expectation-maximization Bayesian least absolute shrinkage and selection operator (ISIS EM-BLASSO) [28].
In single-locus GWAS, usually the Bonferroni correction is applied. The Bonferroni correction incorporates the number of tests by dividing the significance threshold by the number of tests and is used in multiple testing scenarios to reduce the type I error since the likelihood of type I errors increases as the number of tests increases [29]. In complex traits, using single-locus methods bears the risk of missing associated SNPs since they do not surpass the significance threshold [30]. The implementation of multi-locus models in GWAS can help overcome this risk since there is no need for the often excessively conservative Bonferroni correction because of the multi-locus nature of these methods considering every SNP simultaneously [24]. In multi-locus GWAS, the analyses are conducted in two stages. In the first stage, an initial genome-wide single-locus screening selects all SNPs that surpass a moderate significance threshold of 0.01 or 0.005, depending on the multi-locus method. In the second stage, multi-locus models are implemented only for the subset of SNPs that passed the initial screening stage, and SNP effects are estimated simultaneously [26,28]. For the second stage, a significance threshold for the logarithm of the odds (LOD) score of 3 was applied, which is recommended based on the current literature [28]. The study design reached enough power for the detection of significant loci when loci are in high linkage disequilibrium with the potential causative mutation. The perfect association enables reaching p-values of 1.714 × 10−15, even when 4, 6 and 8 cases are misclassified by the markers p-values of 1.247 × 10−11, 6.49 × 10−10 and 2.50 × 10−8 are obtained, respectively (Supplementary Table S3).
Data quality control was performed using PLINK v.1.9 (www.cog-genomics.org/plink/1.9/, accessed on 1 August 2023) [31]. Applied quality criteria included minor allele frequency (MAF) of >0.05 and genotyping rate per SNP and animal of >0.9. Of the 39,246 SNPs that met these quality criteria, only SNPs identified by at least two of the five applied multi-locus methods were included in the putative candidate gene search. Manhattan and Q-Q plots were generated using the calculated −log10 p-values from the multi-locus GWAS using the mrMLM v4.0.2 package in R software, version 4.1.1.

2.6. Putative Candidate Gene Identification

Linkage Disequilibrium (LD) for significantly associated SNPs was calculated using the −r2 command from PLINK 1.90 (www.cog-genomics.org/plink/1.9/, accessed on 1 August 2023) [31], and a threshold of >0.7 for squared correlation between a SNP and another SNP within the 1 Mb region was applied.
Genes within these genomic regions were screened using the current reference genome from NCBI F. catus_Fca126_mat1.0. (https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_018350175.1/, accessed on 4 September 2023) [21]. Putative candidate genes were functionally annotated according to the GeneCards database [32] and literature. The determined genes were further compared to a list of genes associated with AA-amyloidosis from NCBI (https://www.ncbi.nlm.nih.gov/gene/?term=AA-amyloidosis, accessed on 4 September 2023). Furthermore, we compared the results of the present study to the results of the multi-omic study in Abyssinian cats [14]. Potential networks between putative candidate genes and all genes in the region with LD to the associated SNP were analyzed using GeneMANIA [33].

3. Results

3.1. Multi-Locus GWAS

The multi-locus GWAS revealed eight SNPs with a significant association with hereditary AA-amyloidosis in mixed Siamese/Oriental cats with LOD scores > 3 that were detected using at least two multi-locus methods simultaneously (Figure 1 and Table 1). All SNPs that exceeded the LOD scores > 3 were only detected using one multi-locus method, which are provided in Supplementary Table S4.
The method FASTmrMLM detected five of the eight significant SNPs, whereas mrMLM, FASTmrEMMA, ISIS EM-BLASSO and pLARmEB detected four of the eight significant SNPs each (Supplementary Table S5). Of the eight identified SNPs, three are located on Felis catus chromosome (FCA) A1, three on FCA D1 and one each on FCA D2 and FCA D3. Two of the significantly associated SNPs on FCA A1 are located within an interval of 1.5 Mb. The most significantly associated SNP is CHRA1.33307498 on FCA A1 (LOD score 5.40–12.65), which also explains the highest proportion of the phenotypic variance in the mrMLM (26.61%), FASTmrMLM (28.50%) and FASTmrEMMA (16.75%) methods. For pLARmEB, the SNP CHRA1.35096216 explains the highest proportion of the phenotypic variance at 17.36% and the ISIS EM-BLASSO CHRD2.113529492 at 27.57%. The SNP identified by most multi-locus methods (n = 4; FASTmrMLM, FASTmrEMMA, pLARmEB and ISIS EM-BLASSO) is CHRD2.113529492 on FCA D2.
The phenotypic variance estimated for the trait is 0.2067 (Supplementary Table S5). The five multi-locus methods, mrMLM, FASTmrMLM, FASTmrEMMA, pLARmEB and ISIS EM-BLASSO explained 93.55%, 65.05%, 40.58%, 87.74% and 83.91% of the phenotypic variance, respectively, when considering every identified SNP (Supplementary Table S6).
A quantile–quantile plot (Q-Q-plot) shows the observed p-values plotted against the expected p-values under the null hypothesis (Supplementary Figure S2). The distribution of AA-amyloidosis-affected Siamese/Oriental cats by their genotypes of the significantly associated SNPs is shown in Supplementary Table S7. The susceptible genotype of four SNPs explained more than 60% of the cases when at least ten cases were present for the most associated genotype. By combining the four most strongly associated genotypes into a concatenated genotype, the phenotypes of all but three animals were correctly assigned.

3.2. Candidate Gene Search

The 1 Mb upstream and downstream of the identified genomic regions significantly associated with SNPs were screened and 55 genes in these regions were functionally annotated (Supplementary Table S8). Two of the significantly associated SNPs (CHRUN5.1424878 and CHRD1.112359805) are located within introns of genes (IGSSF9B and NAV2), whereas the other SNPs are located in intergenic regions. Functional annotation of these genes provided eight putative candidate genes in close vicinity to three of the significantly associated SNPs (Table 2).
Among the genes retrieved in LD with associated SNPs, SAA1 was identified, which encodes for the precursor of Amyloid A and thus may be a putative candidate gene. Two other putative candidate genes, namely SLC10A2 and ZDHHC13, have been reported to be associated with amyloidosis in humans and/or mice in previous studies (Table 2). ZDHHC13 may be a stronger putative candidate as it plays a role in thyroid gland disease in mice and the thyroid gland is one of the main target organs of AA-amyloidosis in Siamese and Oriental shorthair cats. Furthermore, a genomic region on FCA D1 harboring four genes (HPS5, GTF2H1, LDHA and LDHC) may be considered relevant as this region has a major impact on the precursor of Amyloid A, SAA, in humans (Table 2). The analysis of the genetic interactions of the genes within the 1 Mb genomic regions surrounding the significantly associated SNPs and the eight putative candidate genes revealed shared protein domains, illustrated with green lines, as well as physical interactions, illustrated with red lines, and co-expressions, illustrated with purple lines (Supplementary Figures S3 and S4).

4. Discussion

We identified eight significantly associated SNPs on four cat chromosomes (FCA A1, D1, D2, and D3), confirmed by at least two approaches with mrMLM. The most consistent genomic region among the mrMLM methods employed was identified on FCA A1. Within this 1.5 Mb region, each of the five multi-locus methods identified one significantly AA-amyloidosis-associated SNP (CHRA1.33307498 or CHRA1.35096216). These two SNPs on FCA A1 explained 16.75–28.50% of the phenotypic variance, depending on the multi-locus method applied. Each multi-locus method also identified at least one of the three significantly associated SNPs on FCA D1 (CHRUN5.1424878, CHRD1.83154354 and CHRD1.112359805). The simultaneous detection of the same SNP which was found to be significantly associated with AA-amyloidosis by at least two multi-locus methods corroborates the likelihood of the association of the SNP with the trait. The results of the present study propose the involvement of several loci in the development of AA-amyloidosis in the investigated breeding population of mixed Siamese/Oriental cats. These findings are consistent with the results of the complex segregation analysis, according to which the inheritance of systemic AA amyloidosis is most likely not monogenic but genetically complex [11].
A recent report based on whole genome sequencing (WGS) data of cats affected by AA-amyloidosis concludes that AA-amyloidosis in Abyssinian cats is also not a monogenic trait [14]. Genome variant detection revealed 101 missense variants and one frameshift deletion as disease-associated when only considering the coding gene variants in this WGS study [14]. Following candidate gene search within these variants revealed Amyloid Beta Precursor Protein (APP) and Transthyretin (TTR) genes that have previously been reported to be associated with the deposition of amyloid in humans [14,15,17]. APP and TTR are not located within the genomic regions surrounding the significantly AA-amyloidosis-associated SNPs in the present study. Nevertheless, the putative candidate gene PTPRE is located within the genomic region of CHRD2.113529492 on FCA D2, and a single nucleotide variant in this gene was solely identified in affected Abyssinian cats [14]. Protein Tyrosine Phosphatase Receptor Type E gene (PTPRE) encodes for a member of the protein tyrosine phosphatase family, which regulates various different cellular processes such as cell differentiation and oncogenic transformation [37]. This gene has further been reported to play a role in different forms of cancer including thyroid carcinoma [36]. The thyroid gland was affected by amyloid deposition in 100% of the cases in the present study and in one Abyssinian cat in the multi-omic study, thus a somehow compromised background of the thyroid gland might play a role in the affectedness of this organ. Nevertheless, PTPRE has not been reported in the context of amyloidosis and thus this assumption is speculative. Future analyses of the WGS data of identified genomic regions should nonetheless include this gene to explore a possible link between PTPRE and systemic AA-amyloidosis. Concordant with the multi-omic analyses performed in Abyssinian cats, we consider that a compromised organic background of the main target organ could facilitate the deposition of amyloid or might play a role in the differences in the mainly affected organs in the different cat breeds [14]. Even though there is a possibility of a common genetic background of AA-amyloidosis in Siamese/Oriental cats, Abyssinian cats and other feline species, it is also important to keep in mind that the genetic background could differ between these breeds and species.
It has previously been shown that human chromosome 11 and FCA D1 are syntenic for large segments [38]. A region within 1 Mb one of the significantly associated SNPs (CHRD1.112359805) on FCA D1 harbors the same genes, including HPS5, GTF2H1, LDHA and LDHC, as the region 11p15.5-p13 in humans. The region 11p15.5-p13 has a major impact on SAA, the precursor protein of Amyloid A, and thus represents a highly relevant genomic region for hereditary AA-amyloidosis in cats [35]. This region also harbors Serum Amyloid A1 (SAA1), an important putative candidate gene for AA-amyloidosis. A mutation in the SAA1 promoter causes a hereditary form of AA-amyloidosis in humans; thus, SAA1 could be of high importance in AA-amyloidosis in cats [19]. Nevertheless, a potential involvement of SAA1 in hereditary AA-amyloidosis in cats has not been reported yet. The WGS data of SAA1 need to be the focus of future investigations to review if SAA1 is directly involved in hereditary AA-amyloidosis in cats.
Approximately 340 kb downstream of CHRD1.112359805, the gene Zinc Finger DHHC-Type Palmitoyltransferase 13 (ZDHHC13) is located. The encoded enzyme modifies proteins via palmitoylation, which has been shown to be an important regulator of protein trafficking and protein stability [39]. Interestingly, a mutation in ZDHHC13 caused a severe phenotype in mice including systemic AA- and AL-amyloidosis. Therefore, a variant in ZDHHC13 might also play a role in systemic AA-amyloidosis in cats [20].
In other forms of amyloidosis, such as Aβ–amyloidosis in Alzheimer’s disease (AD), instead of a single gene, multiple genes mostly influence the pathogenesis.
Around 510 kb downstream of CHRA1.87876182 Solute Carrier Family 10 Member 2 (SLC10A2) is located. SLC10A2 encodes a sodium/bile acid cotransporter and is a crucial factor for cholesterol homeostasis [40]. A SNP in close vicinity to SLC10A2 is associated with late-onset Alzheimer’s disease (LOAD) in humans [34]. The potential role in the pathogenesis of LOAD has not been verified yet. Since this gene could be involved in the Aβ-amyloid accumulation in AD patients and not in AA-amyloidosis, drawing comparisons between these conditions is difficult. Nevertheless, the common occurrence of deposits of misfolded, insoluble proteins may be the result of some common features in the development of these conditions and thus justifies a closer look into genes possibly involved in Aβ-amyloid accumulation.
Histopathological confirmation of the clinical diagnosis of accredited veterinarians was available in 16/20 cases. The remaining four cases were closely related to the cases with histopathological confirmation and the similarities in age at death due to clinical signs of amyloidosis corroborate the assumption of the presence of hereditary AA-amyloidosis in these cats.
The results of the present study need to be assessed in a larger cohort of other Siamese/Oriental cats from different populations and cats from other breeds to evaluate if the results may be valid in other Siamese/Oriental breeding populations, or if the findings can solely be confirmed in individuals belonging to the investigated cattery. Inbreeding coefficients were tested for association with the amyloidosis status in the same population. However, we came to the conclusion that inbreeding effects were not significant [11]. In the present study, we calculated genomic inbreeding based on Runs of Homozygosity, but we were not able to show significant effects. Variants that explained a large amount of phenotypic variance should be validated to see if the genomic regions harboring these variants show a significant risk of developing AA-amyloidosis in other Siamese/Oriental cats. If this is the case, searching for further variants in these genomic regions seems advisable.
Depending on the multi-locus method, GWAS results from the present study were able to explain a large amount of the phenotypic variance, but not all, which might have occurred since we excluded SNPs from putative candidate gene searches that did not reach the significance threshold. A limitation of GWAS in complex polygenic traits is that SNPs that do not surpass the significance threshold could nevertheless be associated with the investigated trait if their effect is small. Such variants could play a role in deciphering the genetic background of the trait. We aimed at reducing this risk by choosing different multi-locus methods, which have previously been proven more suitable for polygenic traits [30,41].
Considering this, further work should explore genomic regions that were identified in the present study with WGS data and screen for variants involved in AA-amyloidosis in Siamese/Oriental cats. The regions explaining the highest amounts of phenotypic variance on FCA A1, D1 and D2 and the gene SAA1 should be the focus of these investigations. Increased knowledge of the underlying genetic background in Siamese/Oriental cats will help decipher the pathogenesis of the condition and enable genomic selection to reduce the prevalence of the lethal condition and further aid in identifying similarities and differences with other feline species. The limitations of this study are the low number of cases and controls as well as the origin of cases and controls from a smaller subpopulation of Siamese/Oriental cats. A sampling of cases and controls is difficult because cat breeders are not willing to provide euthanized animals for necropsy and donate samples from older and healthy cats. Due to these sampling limitations, the power to identify associated variants is reduced. Therefore, additional genetic and structural variants not detectable in these data may be relevant to AA amyloidosis. A further study with independent and more data is warranted to confirm these findings and to show whether a smaller number of genetic variants can capture more than 90% of the phenotypic trait variance. A similar approach to a previous report on Abyssinian cats [14], using cats of the genetically predisposed breed affected with AA amyloidosis and controls from a large sample of very different cat breeds not genetically disposed to AA amyloidosis, may also help to corroborate the present findings. The variants and genomic regions identified in this study may, therefore, be useful for further population studies to confirm possible associations. In addition, the identified putative candidate genes and variants may also be a source for screening whole-genome sequencing data and meta-analyses.

5. Conclusions

Hereditary AA-amyloidosis in mixed Siamese/Oriental cats is not a monogenic trait but a complex condition. The region harboring SAA1, HPS5, GTF2H1, LDHA and LDHC, which has a major impact on SAA in humans and the genomic regions that explain a large amount of the phenotypic variance, should be investigated in future studies based on WGS data of cases and controls to identify putatively causal variants for systemic AA-amyloidosis in Siamese/Oriental cats. Further studies with larger datasets are needed to elucidate the complexity of genetic and structural variants involved in AA-amyloidosis in Siamese/Oriental cats.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14122126/s1, Figure S1. Pedigree of the mixed Siamese/Oriental cats. Figure S2. Quantile-quantile (Q-Q) plot of expected p-values versus observed p-values of the GWAS for hereditary AA-amyloidosis in Siamese/Oriental cats. Figure S3. GeneMANIA network of 55 genes located within 1 Mb upstream and downstream of significantly associated SNPs. Figure S4. GeneMANIA network of eight putative candidate genes. Table S1. Animals included in the multi-locus GWAS. Table S2. Localization of amyloid deposits. Table S3. Considerations for identifying significant SNPs in this study depending on the correlation coefficient (r) between a SNP locus that perfectly matches the disease-causing SNP locus and correlated SNP loci. Table S4. Significant SNPs of the multi-locus GWAS. Table S5. Significant SNPs detected by multi-locus methods, phenotypic variance and residual variance of the GWAS-SNPs. Table S6. Phenotypic variance explained by identified SNPs using multi-locus methods. Table S7. Distribution (%) of AA-amyloidosis-affected Siamese/Oriental cats by their genotypes of significantly associated SNPs. Table S8. Putative candidate genes located 1 Mb upstream and downstream of significant SNPs.

Author Contributions

Conceptualization, O.D. and S.L.E.; methodology, O.D. and S.L.E.; software, O.D.; validation, O.D. and S.L.E.; formal analysis, O.D.; investigation, J.H.B., K.H., J.S., T.S., P.W. and O.D.; resources, O.D.; data curation, O.D.; writing—original draft preparation, S.L.E.; writing—review and editing, O.D., J.H.B., K.H., J.S., T.S. and P.W.; visualization, O.D. and S.L.E.; supervision, O.D.; project administration, O.D. All authors have read and agreed to the published version of the manuscript.

Funding

T.S. is the recipient of a PhD fellowship from the National Research Fund, Luxembourg (AFR15686728). This open-access publication was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)—491094227 “Open Access Publication Funding” and the University of Veterinary Medicine Hannover, Foundation.

Institutional Review Board Statement

The Institutional Review Board of the University of Veterinary Medicine Hannover (Foundation) and the responsible state veterinary office of Lower Saxony (registration number 33.9-42502-05-18A284) were obtained.

Informed Consent Statement

Not applicable.

Data Availability Statement

Restrictions apply to the availability of these data. Data were obtained from cat owners and are available upon reasonable request from the authors with the permission of the cat owners.

Conflicts of Interest

The authors declare no conflict of interest and there are no relevant financial or non-financial competing interests to report. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Kumar, V.; Abbas, A.; Fausto, N.; Aster, J. Diseases of the Immune System, Robbins and Cotran Pathologic Basis of Disease; Saunders Elsevier: Philadelphia, PA, USA, 2010; Volume 8, pp. 249–255. [Google Scholar]
  2. Glenner, G.G.; Eanes, E.D.; Page, D.L. The relation of the properties of Congo red-stained amyloid fibrils to the -conformation. J. Histochem. Cytochem. 1972, 20, 821–826. [Google Scholar] [CrossRef]
  3. Johnson, K.H.; Sletten, K.; Werdin, R.E.; Westermark, G.T.; O’Brien, T.D.; Westermark, P. Amino acid sequence variations in protein AA of cats with high and low incidences of AA amyloidosis. Comp. Biochem. Physiol. B 1989, 94, 765–768. [Google Scholar] [CrossRef] [PubMed]
  4. Husebekk, A.; Skogen, B.; Husby, G.; Marhaug, G. Transformation of amyloid precursor SAA to Protein AA and incorporation in amyloid fibrils in vivo. Scand. J. Immunol. 1985, 21, 283–287. [Google Scholar] [CrossRef] [PubMed]
  5. Boyce, J.T.; DiBartola, S.P.; Chew, D.J.; Gasper, P.W. Familial renal amyloidosis in Abyssinian cats. Vet. Pathol. 1984, 21, 33–38. [Google Scholar] [CrossRef] [PubMed]
  6. Van der Linde-Sipman, J.S.; Niewold, T.A.; Tooten, P.C.J.; de Neijs-Backer, M.; Gruys, E. Generalized AA-amyloidosis in Siamese and Oriental cats. Vet. Immunol. Immunopathol. 1997, 56, 1–10. [Google Scholar] [CrossRef]
  7. Beatty, J.A.; Barrs, V.R.; Martin, P.A.; Nicoll, R.G.; France, M.P.; Foster, S.F.; Lamb, W.A.; Malik, R. Spontaneous hepatic rupture in six cats with systemic amyloidosis. J. Small Anim. Pract. 2002, 43, 355–363. [Google Scholar] [CrossRef] [PubMed]
  8. Zuber, R.M. Systemic amyloidosis in Oriental and Siamese Cats. Aust. Vet. Practit. 1993, 23, 66–70. [Google Scholar]
  9. DiBartola, S.P.; Tarr, M.J.; Benson, M.D. Tissue distribution of amyloid deposits in Abyssinian cats with familial amyloidosis. J. Comp. Pathol. 1986, 96, 387–398. [Google Scholar] [CrossRef]
  10. Chew, D.J.; DiBartola, S.P.; Boyce, J.T.; Gasper, P.W. Renal amyloidosis in related Abyssinian cats. J. Am. Vet. Med. Assoc. 1982, 181, 139–142. [Google Scholar]
  11. Struck, A.K.; Klotz, D.; Hülskötter, K.; Wohlsein, P.; Schmitz, J.; Bräsen, J.H.; Distl, O. Complex segregation analysis of familial amyloidosis in Oriental shorthair cats. Vet. J. 2020, 265, 105552. [Google Scholar] [CrossRef]
  12. Godfrey, D.R.; Day, M.J. Generalised amyloidosis in two Siamese cats: Spontaneous liver haemorrhage and chronic renal failure. J. Small. Anim. Pract. 1998, 39, 442–447. [Google Scholar] [CrossRef]
  13. Niewold, T.A.; Van der Linde-Sipman, J.S.; Murphy, C.; Tooten, P.C.; Gruys, E. Familial amyloidosis in cats: Siamese and Abyssinian AA proteins differ in primary sequence and pattern of deposition. Amyloid 1999, 6, 205–209. [Google Scholar] [CrossRef] [PubMed]
  14. Genova, F.; Nonnis, S.; Maffioli, E.; Tedeschi, G.; Strillacci, M.G.; Carisetti, M.; Sironi, G.; Cupaioli, F.A.; Di Nanni, N.; Mezzelani, A.; et al. Multi-omic analyses in Abyssinian cats with primary renal amyloid deposits. Sci. Rep. 2021, 11, 8339. [Google Scholar] [CrossRef]
  15. Connors, L.H.; Lim, A.; Prokaeva, T.; Roskens, V.A.; Costello, C.E. Tabulation of human transthyretin (TTR) variants, 2003. Amyloid 2003, 10, 160–184. [Google Scholar] [CrossRef] [PubMed]
  16. Van Broeckhoven, C.; Haan, J.; Bakker, E.; Hardy, J.A.; Van Hul, W.; Wehnert, A.; Vegter-Van der Vlis, M.; Roos, R.A.C. Amyloid β; Protein precursor gene and hereditary cerebral hemorrhage with amyloidosis (dutch). Science 1990, 248, 1120–1122. [Google Scholar] [CrossRef] [PubMed]
  17. Hendriks, L.; van Duijn, C.M.; Cras, P.; Cruts, M.; Van Hul, W.; van Harskamp, F.; Warren, A.; McInnis, M.G.; Antonarakis, S.E.; Martin, J.J.; et al. Presenile dementia and cerebral haemorrhage linked to a mutation at codon 692 of the β–amyloid precursor protein gene. Nat. Genet. 1992, 1, 218–221. [Google Scholar] [CrossRef] [PubMed]
  18. Kamino, K.; Orr, H.T.; Payami, H.; Wijsman, E.M.; Alonso, M.E.; Pulst, S.M.; Anderson, L.; O’Dahl, S.; Nemens, E.; White, J.A.; et al. Linkage and mutational analysis of familial Alzheimer disease kindreds for the APP gene region. Am. J. Hum. Genet. 1992, 51, 998–1014. [Google Scholar]
  19. Sikora, J.; Kmochová, T.; Mušálková, D.; Pohludka, M.; Přikryl, P.; Hartmannová, H.; Hodaňová, K.; Trešlová, H.; Nosková, L.; Mrázová, L.; et al. A mutation in the SAA1 promoter causes hereditary amyloid A amyloidosis. Kidney Int. 2022, 101, 349–359. [Google Scholar] [CrossRef]
  20. Saleem, A.N.; Chen, Y.-H.; Baek, H.J.; Hsiao, Y.-W.; Huang, H.-W.; Kao, H.-J.; Liu, K.-M.; Shen, L.-F.; Song, I.w.; Tu, C.-P.D.; et al. Mice with alopecia, osteoporosis, and systemic amyloidosis due to mutation in Zdhhc13, a gene coding for palmitoyl acyltransferase. PLoS Genet. 2010, 6, e1000985. [Google Scholar] [CrossRef]
  21. O’Leary, N.A.; Wright, M.W.; Brister, J.R.; Ciufo, S.; Haddad, D.; McVeigh, R.; Rajput, B.; Robbertse, B.; Smith-White, B.; Ako-Adjei, D.; et al. Reference sequence (RefSeq) database at NCBI: Current status, taxonomic expansion, and functional annotation. Nucleic Acids Res. 2016, 44, D733–D745. [Google Scholar] [CrossRef]
  22. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  23. Zhang, Y.W.; Tamba, C.L.; Wen, Y.J.; Li, P.; Ren, W.L.; Ni, Y.L.; Gao, J.; Zhang, Y.M. mrMLM v4.0.2: An R Platform for multi-locus genome-wide association studies. Genom. Proteom. Bioinform. 2020, 18, 481–487. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, S.-B.; Feng, J.-Y.; Ren, W.-L.; Huang, B.; Zhou, L.; Wen, Y.-J.; Zhang, J.; Dunwell, J.M.; Xu, S.; Zhang, Y.-M. Improving power and accuracy of genome-wide association studies via a multi-locus mixed linear model methodology. Sci. Rep. 2016, 6, 19444. [Google Scholar] [CrossRef]
  25. Tamba, C.L.; Zhang, Y.M. A fast mrMLM algorithm for multi-locus genome-wide association studies. biorXiv 2018, 341784. [Google Scholar] [CrossRef]
  26. Wen, Y.J.; Zhang, H.; Ni, Y.L.; Huang, B.; Zhang, J.; Feng, J.Y.; Wang, S.B.; Dunwell, J.M.; Zhang, Y.M.; Wu, R. Methodological implementation of mixed linear models in multi-locus genome-wide association studies. Brief. Bioinform. 2018, 19, 700–712. [Google Scholar] [CrossRef] [PubMed]
  27. Zhang, J.; Feng, J.Y.; Ni, Y.L.; Wen, Y.J.; Niu, Y.; Tamba, C.L.; Yue, C.; Song, Q.; Zhang, Y.M. pLARmEB: Integration of least angle regression with empirical Bayes for multilocus genome-wide association studies. Heredity 2017, 118, 517–524. [Google Scholar] [CrossRef] [PubMed]
  28. Tamba, C.L.; Ni, Y.L.; Zhang, Y.M. Iterative sure independence screening EM-Bayesian LASSO algorithm for multi-locus genome-wide association studies. PLoS Comput. Biol. 2017, 13, e1005357. [Google Scholar] [CrossRef]
  29. Neyman, J.; Pearson, E.S. On the use and interpretation of certain test criteria for purposes of statistical inference: Part I. Biometrika 1928, 20A, 175–240. [Google Scholar] [CrossRef]
  30. Xu, Y.; Yang, T.; Zhou, Y.; Yin, S.; Li, P.; Liu, J.; Xu, S.; Yang, Z.; Xu, C. Genome-wide association mapping of starch pasting properties in maize using single-locus and multi-locus models. Front. Plant Sci. 2018, 9, 1311. [Google Scholar] [CrossRef]
  31. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-generation PLINK: Rising to the challenge of larger and richer datasets. GigaScience 2015, 4, 7. [Google Scholar] [CrossRef]
  32. Safran, M.; Dalah, I.; Alexander, J.; Rosen, N.; Iny Stein, T.; Shmoish, M.; Nativ, N.; Bahir, I.; Doniger, T.; Krug, H.; et al. GeneCards version 3: The human gene integrator. Database 2010, 2010, baq020. [Google Scholar] [CrossRef]
  33. Warde-Farley, D.; Donaldson, S.L.; Comes, O.; Zuberi, K.; Badrawi, R.; Chao, P.; Franz, M.; Grouios, C.; Kazi, F.; Lopes, C.T.; et al. The GeneMANIA prediction server: Biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010, 38, W214–W220. [Google Scholar] [CrossRef] [PubMed]
  34. Mez, J.; Chung, J.; Jun, G.; Kriegel, J.; Bourlas, A.P.; Sherva, R.; Logue, M.W.; Barnes, L.L.; Bennett, D.A.; Buxbaum, J.D.; et al. Two novel loci, COBL and SLC10A2, for Alzheimer’s disease in African Americans. Alzheimers. Dement. 2017, 13, 119–129. [Google Scholar] [CrossRef]
  35. Marzi, C.; Albrecht, E.; Hysi, P.G.; Lagou, V.; Waldenberger, M.; Tönjes, A.; Prokopenko, I.; Heim, K.; Blackburn, H.; Ried, J.S.; et al. Genome-Wide Association Study identifies two novel regions at 11p15.5-p13 and 1p31 with major impact on acute-phase Serum Amyloid A. PLoS Genet. 2010, 6, e1001213. [Google Scholar] [CrossRef] [PubMed]
  36. Peng, C.; Zhang, C.; Yu, W.; Li, L.; Zhang, Z.; Liu, T.; Zhang, Y.; Fan, G.; Huangfu, H. Receptor Type Protein Tyrosine Phosphatase Epsilon (PTPRE) plays an oncogenic role in thyroid carcinoma by activating the AKT and ERK1/2 signaling pathway. Curr. Cancer Drug Targets 2023, 23, 471–481. [Google Scholar] [CrossRef] [PubMed]
  37. Aschner, Y.; Downey, G.P. The importance of tyrosine phosphorylation control of cellular signaling pathways in respiratory disease: pY and pY not. Am. J. Respir. Cell Mol. Biol. 2018, 59, 535–547. [Google Scholar] [CrossRef] [PubMed]
  38. Rettenberger, G.; Klett, C.; Zechner, U.; Bruch, J.; Just, W.; Vogel, W.; Hameister, H. ZOO-FISH analysis: Cat and human karyotypes closely resemble the putative ancestral mammalian karyotype. Chromosome Res. 1995, 3, 479–486. [Google Scholar] [CrossRef] [PubMed]
  39. Abrami, L.; Leppla, S.H.; van der Goot, F.G. Receptor palmitoylation and ubiquitination regulate anthrax toxin endocytosis. J. Cell Biol. 2006, 172, 309–320. [Google Scholar] [CrossRef]
  40. Boyer, J.L. Bile formation and secretion. Compr. Physiol. 2013, 3, 1035–1078. [Google Scholar] [CrossRef]
  41. Li, C.; Fu, Y.; Sun, R.; Wang, Y.; Wang, Q. Single-locus and multi-locus genome-wide association studies in the genetic dissection of fiber quality traits in upland cotton (Gossypium hirsutum L.). Front. Plant Sci. 2018, 9, 1083. [Google Scholar] [CrossRef]
Genes 14 02126 g001
Figure 1. Manhattan plot of p-values and LOD-scores for AA-amyloidosis in mixed Siamese/Oriental cats. Each dot represents one SNP of the multi-locus GWAS. Blue and green dots represent −log10 p-values of SNPs identified by one multi-locus method. The genome-wide −log10 p-values and LOD scores for each SNP effect are plotted against their chromosomal position. Pink dots represent SNPs identified by at least two of the applied methods with LOD scores > 3.
Figure 1. Manhattan plot of p-values and LOD-scores for AA-amyloidosis in mixed Siamese/Oriental cats. Each dot represents one SNP of the multi-locus GWAS. Blue and green dots represent −log10 p-values of SNPs identified by one multi-locus method. The genome-wide −log10 p-values and LOD scores for each SNP effect are plotted against their chromosomal position. Pink dots represent SNPs identified by at least two of the applied methods with LOD scores > 3.
Genes 14 02126 g001
Table 1. SNPs significantly associated with hereditary AA-amyloidosis in mixed Siamese/Oriental cats and detected by at least two methods simultaneously in the multi-locus GWAS.
Table 1. SNPs significantly associated with hereditary AA-amyloidosis in mixed Siamese/Oriental cats and detected by at least two methods simultaneously in the multi-locus GWAS.
SNP-IDFCAPositionMAFLOD ScorePVE (%)Method
CHRA1.33307498A129,635,1720.3612.6526.61mrMLM
0.355.4028.50FASTmrMLM
0.3510.2716.75FASTmrEMMA
CHRA1.35096216A131,080,8870.293.7817.36pLARmEB
0.296.0518.54ISIS EM-BLASSO
CHRA1.87876182A174,262,2750.397.2112.22mrMLM
0.404.304.21FASTmrMLM
0.403.883.74FASTmrEMMA
CHRUN5.1424878D130,498,4770.315.2015.65FASTmrMLM
0.316.8415.35ISIS EM-BLASSO
CHRD1.83154354D152,625,2460.453.073.23FASTmrEMMA
0.454.208.47pLARmEB
CHRD1.112359805D175,099,1820.236.3411.86mrMLM
0.233.491.96FASTmrMLM
0.235.595.01pLARmEB
CHRD2.113529492D284,125,5130.356.3214.73FASTmrMLM
0.356.7710.81FASTmrEMMA
0.356.2916.56pLARmEB
0.3513.3527.57ISIS EM-BLASSO
CHRA1.170535359D312,578,0420.443.325.39mrMLM
0.446.9110.35ISIS-EM-BLASSO
FCA: Felis catus chromosome according to F. catus_Fca126_mat1.0; Position: SNP position in base pairs according to F. catus_Fca126_mat1.0; MAF: minor allele frequency; LOD score: logarithm of the odds scores; PVE: percentage of phenotypic variance explained by multi-locus methods.
Table 2. Putative candidate genes for AA-amyloidosis in mixed Siamese/Oriental cats and their chromosomal position on Felis catus (FCA) chromosome according to F. catus_Fca126_mat1.0.
Table 2. Putative candidate genes for AA-amyloidosis in mixed Siamese/Oriental cats and their chromosomal position on Felis catus (FCA) chromosome according to F. catus_Fca126_mat1.0.
SNP-IDFCAGeneStartEndRelated Functions
CHRA1.87876182A1SLC10A273,755,03673,774,055Late onset Alzheimer’s disease [34]
CHRD1.112359805D1SAA174,193,15874,196,364AA-amyloidosis [19]
HPS5
GTF2H1
LDHA
LDHC		74,220,513
74,280,919 74,338,285
74,353,740		74,280,965 74,318,855
74,349,667 74,397,837		Major impact on SAA [35]
ZDHHC1374,760,62574,812,445AA-/AL-amyloidosis [20]
CHRD2.113529492D2PTPRE83,094,72283,258,667Thyroid carcinoma [36]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Esders, S.L.; Hülskötter, K.; Schreiner, T.; Wohlsein, P.; Schmitz, J.; Bräsen, J.H.; Distl, O. Single Nucleotide Polymorphisms Associated with AA-Amyloidosis in Siamese and Oriental Shorthair Cats. Genes 2023, 14, 2126. https://doi.org/10.3390/genes14122126

AMA Style

Esders SL, Hülskötter K, Schreiner T, Wohlsein P, Schmitz J, Bräsen JH, Distl O. Single Nucleotide Polymorphisms Associated with AA-Amyloidosis in Siamese and Oriental Shorthair Cats. Genes. 2023; 14(12):2126. https://doi.org/10.3390/genes14122126

Chicago/Turabian Style

Esders, Stella L., Kirsten Hülskötter, Tom Schreiner, Peter Wohlsein, Jessica Schmitz, Jan H. Bräsen, and Ottmar Distl. 2023. "Single Nucleotide Polymorphisms Associated with AA-Amyloidosis in Siamese and Oriental Shorthair Cats" Genes 14, no. 12: 2126. https://doi.org/10.3390/genes14122126

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop