Classification of feline hypertrophic cardiomyopathy-associated gene variants according to the American College of Medical Genetics and Genomics guidelines

Introduction: The correct labeling of a genetic variant as pathogenic is important as breeding decisions based on incorrect DNA tests can lead to the unwarranted exclusion of animals, potentially compromising the long-term health of a population. In human medicine, the American college of Medical Genetics (ACMG) guidelines provide a framework for variant classification. This study aims to apply these guidelines to six genetic variants associated with hypertrophic cardiomyopathy (HCM) in certain cat breeds and to propose a modified criterion for variant classification.

Methods: Genetic samples were sourced from five cat breeds: Maine Coon, Sphynx, Ragdoll, Devon Rex, and British Short- and Longhair. Allele frequencies were determined, and in the subset with phenotypes available, odds ratios to determine the association with HCM were calculated. In silico evaluation followed with joint evidence and data from other publications assisting in the classification of each variant.

Results: Two variants, MYBPC3:c.91G > C [A31P] and MYBPC3:c.2453C > T [R818W], were designated as pathogenic. One variant, MYH7:c.5647G > A [E1883K], was found likely pathogenic, while the remaining three were labeled as variants of unknown significance.

Discussion: Routine genetic testing is advised solely for the MYBPC3:c.91G > C [A31P] in the Maine Coon and MYBPC3:c.2453C > T [R818W] in the Ragdoll breed. The human ACMG guidelines serve as a suitable foundational tool to ascertain which variants to include; however, refining them for application in veterinary medicine might be beneficial.

1 Introduction

Hypertrophic cardiomyopathy (HCM) is the most common heart disease in cats. Prevalence estimates vary between 3% (Pawpeds database, purebred cats, n = 19.866), 8% (purebred cats, n = 144), and even 15% (shelter cats, majority of domestic shorthairs and a small number of purebred cats, n = 780) (1–3). Based upon a 2012 estimate of feline pet ownership and extrapolating prevalence of 0.5% (= conservative estimate based on human prevalence estimates) and 8%, upwards of 470.000 to 7.5 million cats, respectively, could be affected by HCM (4, 5). Clinical signs can present as early as 3 months; although the median age of diagnosis has been reported at approximately 6.5 years, it may vary when evaluated in different populations (4, 6). HCM occurs in both sexes, in all breeds (but in some more frequently), and is characterized by primary concentric left ventricular hypertrophy with focal to diffuse myocardial thickening, myocardial fibrosis, and myocyte fiber disarray (7, 8). These features may lead to impaired diastolic filling, atrial dilation, congestive heart failure, thromboembolism, or sudden cardiac death (4, 9). Unless one or more of these events occur, affected cats may not express overt clinical signs, and many remain pre-clinical (10).

HCM in humans and cats shares several morphological and clinical similarities (11–14). HCM was first described as an autosomal dominant disorder in humans in 1960 (15). The first sarcomeric gene variant for HCM was identified in 1989 (16). Currently, more than 1700 gene variants have been associated with HCM in humans. The myosin-binding protein C3 (MYBPC3) gene and myosin heavy chain 7 (MYH7) gene contain most of the disease-associated variants (17). In domestic cats, six variants in four different genes (MYBPC3:c.91G > A [A31P], MYBPC3:c.220G > A [A74T], MYBPC3:c.2453C > T [R818W (XP_019667956.1); originally known a R820W], ALMS1:c.7384G > C [G2462R], TNNT2:c.95-108G > A, and MYH7:c.5647G > A [E1883]) have been reported to be associated with feline HCM (18–22).

While the association with HCM has been independently confirmed for some of the variants, this is not the case for all six variants, for example, the MYPBC3:c.220G > A [A74T] variant. This may in part be due to the lack of objective guidelines to evaluate variants’ pathogenicity. In human medicine, the American College of Medical Genetics (ACMG) guidelines are used for the interpretation and classification of sequence variants (23). Variants are classified into five different categories (pathogenic, likely pathogenic, unknown significance, likely benign, and benign) based on a comprehensive analysis that includes criteria such as functional studies, computational predictive models, co-segregation with the disease in families, and previously published research. Only when sufficient evidence is available to support the classification of the variant into one of the first two categories can a variant be used for clinical decision-making and genetic testing. The reason to only include variants that meet those criteria is obvious: the benefit of genetic tests is that they are quick and allow the early identification of cats predisposed to develop certain diseases (10). This early detection can be crucial for making informed breeding decisions and/or to decide at an early age which individuals should be retained as breeding stock (24). DNA tests based on variants that are not correctly associated with disease can lead to unnecessary exclusion of animals from breeding, can reduce confidence in DNA testing, and can create reluctance to use other reliable tests. Ultimately, this can reduce the effective population size and negatively impact the long-term health and diversity of animal populations (25). The objective of this study was to perform this evaluation of pathogenicity for the first time on all six variants reported to be associated with HCM using the existing human ACMG guidelines and to report a modified criterion that might help distinguish future variants in veterinary medicine.

2 Materials and methods

2.1 Study design

This study was established as a multicenter four-phase study and was based on (1) the prospective genotyping of samples of five breeds (British Short- and Longhair, Ragdoll, Sphynx, Maine Coon, and Devon Rex) for which the sample size requirements were met (see below), based on the samples stored at universities and repositories of Cornell University (USA), Ghent University (Belgium), Swedish University of Agricultural Sciences (Sweden), VetAgro Sup (France), Utrecht University (the Netherlands), and University of Milan (Italy) (= batch one) and (2) retrospective collection of data from two commercial laboratories (Antagene and Wisdom Panel) for the MYBPC3:c.91G > C [A31P] and MYBPC3:c.2453C > T [R818W] variants (= batch two). A subset of the samples provided by the Wisdom panel has been published in an earlier study with a different scope (26). The number of overlapping samples is depicted in Supplementary Table S1. For every breed in this study, these samples only represent a minority of what was available in the current study.

In the first phase, allele frequencies of all six published variants (the three in MYBPC3 and one each in TNNT2, MYH7, and ALMS1) were determined. Inclusion criteria for samples used in this phase were (1) belonging to one of the five aforementioned breeds and (2) the necessary material had to be present to allow genotyping or genotypes had to be available. To be included in the second phase, an additional inclusion criterion had to be fulfilled: only samples for which a phenotype was available were used to calculate odds ratios (ORs) to evaluate the association between the variant and HCM. In this phase, three exclusion criteria were also used: one linked to the phenotype, one linked to the breed, and one linked to the genotype, as detailed in the statistical analysis section. In the third phase, all variants were assessed with in silico tools to predict their pathogenicity. In the final phase, the aforementioned evidence, together with additional data from other publications, was used to classify each variant using the ACMG guidelines.

2.2 Power analysis

A power analysis was conducted with the R package epiR to determine the minimal sample size necessary to detect a variant with a certain odds ratio (27). The odds ratio used for these calculations was based on the allelic odds ratio derived from the paper in which the variant was described for the first time. The lowest and highest disease prevalence reported in literature were used, leading to a HCM prevalence of 3% (= worst case scenario in terms of sample size) or 15% (= best case scenario in terms of sample size) in the group without the disease-causing allele, respectively (1, 3). The prevalence that resembled the prevalence in our sample most closely was retained. Furthermore, a balanced number of affected and healthy cats and a desired power of 0.80 were assumed. As the odds ratio differed for each variant, the required total sample sizes ranged between 4 and 32. An overview for each variant is presented in Supplementary Table S2.

2.3 Ethics

For batch one, all included cats had informed owner consent. For batch two, the owners provided written consent for data use in research upon submission of samples for commercial genetic testing. Samples were non-invasive buccal swabs or whole blood on EDTA. In some cats, blood was obtained as part of routine clinical procedures for diagnostic purposes, at the request and with the consent of the owner. As these samples were from client-owned cats for which no harmful invasive procedures were performed, there was no animal experimentation according to the legal definition in Europe (Subject 5f of Article1, Chapter I of the Directive 2010/63/UE of the European Parliament and of the Council). A subset of the blood samples was collected as part of a biobanking initiative, approved by the Ethics Committee of the Faculties of Veterinary Medicine and Bioscience Engineering at Ghent University (EC2017-86). Samples from Sweden were collected in agreement with the ethical committee of the Swedish Board of Agriculture [No. C2/12 (2012-02-24),C12/15 (2015-02-27), and Dnr 5.8.18–04682/2020]. Samples sourced from the Cornell Veterinary Biobank were collected under Institutional Animal Care and Use Committee (IACUC) approved animal protocol #2005–0151.

2.4 Genotyping

A detailed overview of the primers, product lengths and type of test, PCR temperatures, and other parameters used for genotyping of the samples of batch one can be found in Supplementary Table S3. For ALMS1:c.7384G > C, the recently published primers were used to avoid allelic drop-out (28). A subset of the samples of batch two (n = 10.851) were genotyped using a commercial custom genotyping microarray [i.e., MyCatDNA and the Optimal Selection Feline tests (26)].

2.5 Phenotyping

HCM phenotyping was performed by echocardiography according to the American College of Veterinary Internal Medicine (ACVIM) consensus guidelines for classification, diagnosis, and management of cardiomyopathy in cats (29). Echocardiography was performed by board-certified veterinary cardiologists, residents under direct supervision of a cardiologist or experienced cardiology clinicians. An HCM phenotype was defined as the presence of diffuse or regional increased LV wall thickness with a non-dilated LV chamber. A combination of two-dimensional (2D)-guided M-mode (right parasternal short axis view) and 2D echocardiography (right parasternal long axis view and right parasternal short axis view) was used to measure LV wall thickness. For M-mode, and for 2D measurement of the LV free wall, end-diastolic wall thickness was measured using a leading edge-to-leading edge technique, whereas a leading edge-to-trailing edge technique was applied for 2D measurement of interventricular septum thickness. At least three cardiac cycles were measured and averaged. Left atrial size was measured via 2D echocardiography in the right parasternal heart base view (left atrium-to-aorta ratio) and right parasternal long axis view (left atrial long axis diameter). Additional features of HCM were also assessed, including papillary muscle hypertrophy, end-systolic chamber obliteration, septal or apical hypertrophy, mitral valve systolic anterior motion, left ventricular outflow tract obstruction, and spontaneous echocontrast. Diastolic function was investigated using pulsed-wave and tissue Doppler imaging. In the realm of feline cardiology, there is no universally accepted standard for determining the normal thickness of the left ventricular wall during diastole. It is overly simplistic to assume that a single numerical value can reliably distinguish between a healthy left ventricular wall and one affected by hypertrophy across all cats. The thickness of the left ventricular wall increases nonlinearly with body weight and is influenced by physiological factors such as hydration level and heart rate. Conditions like hyperthyroidism and systemic hypertension, which are common in older cats, can further impact left ventricular wall thickness, and the presence of these conditions were examined in cats, if suspected, before inclusion.

Nevertheless, for most cats that are not hyperthyroid or hypertensive, of average size and with a normal body condition (weighing between 3.5–5 kg), a left ventricular diastolic wall thickness of 5 mm or less is typically considered within the normal range, while a thickness of 6 mm or more suggests concentric hypertrophy. Cats with a left ventricular wall thickness falling between 5 and 6 mm were termed equivocal and were excluded from analysis. Cats outside this weight range were managed on a case to case basis with the aid of previously published 95% prediction intervals (1).

2.6 Statistical analysis

All analyses were performed in R version 4.1.1. and Excel version 16.54. Allele frequencies were calculated for each breed separately. In the second phase of our study, where ORs and their 95% confidence intervals were computed, only cats for which phenotypical data were available were included. Cats with ambiguous hypertrophic cardiomyopathy (equivocal cats), cats diagnosed with restrictive cardiomyopathy, and healthy cats under 6 years of age were excluded (phenotypical exclusion criteria). This age cutoff for the control group was based on the observation that median or mean age of diagnosis lies approximately at 6 years (4, 30). By setting this as the minimal age, the goal was to obtain a control group in which the overall age would actually be higher than the cutoff as an increased age reduces the probability of misclassifying cats further. The association of each individual variant with HCM was determined (1) in the breed in which the variant was originally described (breed exclusion criterion) and (2) if a variant was not significant, this analysis was repeated in all breeds in which the variant segregated. If several variants were originally described in a breed (i.e., the three variants reported in the Maine Coon), only homozygous wild-type individuals for those variants except the one that was assessed were included to avoid the potential influence of these other variants. For the analysis across breeds, cats had to be homozygous wild type for all variants other than the variant that was evaluated. The exclusion based on genotype thus represents the genotypic exclusion criterion. Both allelic ORs (i.e., the ratio between the odds of having the phenotype and a specific allele and the odds of the same phenotype without having the allele) and genotypic ORs (based on two inheritance patterns, i.e., autosomal dominant and autosomal recessive) were calculated for every variant, using unconditional maximum likelihood estimation and normal approximation [Wald] confidence intervals. The Haldane-Anscombe correction was applied when necessary (31).

2.7 In silico analysis

For each variant, in silico computational prediction methods were employed using three online tools. For every variant, except TNNT2:c.95-108G > A, the following missense prediction tools were used: Polyphen-2¹, Sift² (default settings), and Mutpred2³ (p-value threshold: 0.05). The MYBPC3 [uniprot: M3VYP3_FELCA, (XM_019812396.1 RefSeq)] protein sequence was used for the three MYBPC3 variants (MYBPC3:p.A31P, MYBPC3:p.A74T, MYBPC3:p.R818W). The MYH7 (XM_006932746.4 RefSeq) protein sequence was used for the MYH7:p.E1883K variant. The ALMS1 (A0A5F5XP77_FELCA) protein sequence was used for the ALSM1:c.7384G > C [G2462R] variant.

For the TNNT2:c.95-108G > A variant, the following tools were used to predict the effect on mRNA for splicing of the TNNT2 gene: ASSP⁴ (default parameters), GENSCAN⁵ (organism: vertebrate), Netgene2⁶ (organism: Human). The TNNT2 genomic sequence (ENSFCAG0000004613) was used both with and without the variant as input for each tool. Wild type and variant sequence predicted splice sites were compared.

2.8 The American College of Medical Genetics classification

The ACMG has developed a systematic approach for classifying genetic variants, based on a set of criteria. Each criterion has a specific weight, and the set of fulfilled criteria determines the final classification, which ranges from pathogenic to benign. A concise overview of the criteria that were relevant is available in Supplementary Table S4. The classification of the variants into five categories (pathogenic, likely pathogenic, variant of unknown significance, likely benign, and benign) was done in two stages. In the first stage, the ACMG guidelines were used except for the criteria related to allele frequency, in agreement with what has been done in earlier studies (21, 32, 33). In the second stage, the additional effect of using these allelic frequency criteria for classification was evaluated. In this second stage, an alternative cutoff, i.e., the alternative allele frequency threshold (Tv), was also calculated and used (Supplementary Table S5) (34). This alternative cutoff was used with the same weight as in the ACMG guidelines (i.e., a stand-alone benign (BA), a strong benign (BS), or a moderate pathogenic (PM) criterion). To allow distinction with the original criteria, they are denoted as “BA1*,” “BS1*,” and “PM2*” instead of the original “BA1,” “BS1,” and “PM2” allelic frequency cutoff related criteria, respectively.

3 Results

3.1 Allele frequencies

For the calculation of allelic frequencies (Table 1), all available samples regardless of phenotype were used. The calculation of allelic frequencies was based on the combination of samples from the universities/repositories (=batch one) and data from commercial laboratories (=batch two), regardless of whether a phenotype was available. From batch one, the number of samples per breed was 45 British Short- and Longhair, 20 Ragdoll, 107 Sphynx, 97 Maine Coon, and 33 Devon Rex cats. These were complemented by screening results from two commercial laboratories that have already tested for MYBPC3:c.91G > C [A31P] and MYBPC3:c.2453C > T [R818W] for several years, increasing the sample sizes for those two variants up to nearly a thousand to several thousand samples, depending on the breed. An overview of the sample sizes is provided in Table 1. Furthermore, the evolution over time for MYBPC3:c.91G > C [A31P] and MYBPC3:c.2453C > T [R818W] was available for a subset of the samples from one commercial laboratory and is depicted in Supplementary Table S6.

Table 1

Table 1. Overview of the allele frequencies of the six variants reported to be associated with hypertrophic cardiomyopathy in cats.

Overall, one variant (MYH7:c.5647G > A [E1883K]) was not observed in any cat, one variant (MYBPC3:c.2453C > T [R818W]) was found only in the Ragdoll cats, and the MYBPC3:c.91G > C [A31P] variant segregated in the Maine Coon and one Ragdoll. The remaining three variants segregated in all five breeds.

3.2 Association with HCM

Twenty-nine British Short- and Longhair (cases = 16, controls = 13), 16 Ragdoll (cases = 9, controls = 7), 84 Sphynx (cases = 53, controls = 31), 92 Maine Coon (cases = 47, controls = 45), and 25 Devon Rex (cases = 15, controls = 10) were retained after the application of the phenotypical inclusion and exclusion criteria. Subsequently, for the first analysis focusing on the individual breeds in which the variants were initially identified (=Maine Coon, Ragdoll, and Sphynx; the breed exclusion criterion), we restricted our sample to cats that were homozygous for the ‘wild-type’ variant of all other gene variants under investigation in that particular breed. This (genotype) exclusion criterion was put in place to prevent any confounding influence from other variants. After eliminating all non-homozygous wild-type cats for the MYBPC3:c.220G > A [A74T] and TNNT2:c.95-108G > A variants, the total number of Maine Coons available for the OR calculation of MYBPC3:c.91G > C [A31P] was 37. By employing the same methodology, 27 cats and 28 cats were retained for the MYBPC3:c.220G > A [A74T] and TNNT2:c.95-108G > A variant, respectively. The sample sizes for each variant, as well as the median age and the range for each of the three breeds in which the variants were originally found, are reported in Table 2. Overall, the median age for the cases was 4 years (range: 0.25–16.50), while the median age for the controls was 8 years (range: 6.00–23.00).

Table 2

Table 2. Overview of the allelic and genotypic odds ratios (ORs) of the six variants reported to be associated with HCM in cats.

A significant allelic OR was found for the MYBPC3:c.91G > C [A31P] variant and the MYBPC3:c.2453C > T [R818W] variant. A significant genotypic OR was also found for the MYBPC3:c.91G > C [A31P] (autosomal recessive mode of inheritance). None of the other variants were significantly associated with HCM. ORs could not be calculated for the MYH7:c.5647G > A [E1883K] variant due to its absence in the dataset. For the three variants that were not significant (MYBPC3:c.220G > A [A74T], TNNT2:c.95-108G > A and ALMS1:c.7384G > C [G2462R]), the analysis was repeated across all breeds in which the variant was found. However, none of the point estimates increased (Supplementary Table S7).

3.3 Prediction of harmful effects

Two sets of in silico computational prediction methods were used to investigate the impact of genetic variants on protein function and mRNA splicing, respectively. The algorithms were only consistent for MYBPC3:c.2453C > T [R818W] (predicted pathogenic), ALMS1:c.7384G > C [G2462R] (predicted pathogenic), and MYH7:c.5647G > A [E1883K] (predicted pathogenic). For the other three variants, the algorithms gave inconsistent results (Table 3).

Table 3

Table 3. Overview of the outcomes of six hypertrophic cardiomyopathy-associated variants analyzed with in silico prediction tools.

3.4 ACMG variant classification

The combination of the ORs, the results of the predictive algorithms, and publications regarding functional studies and evidence across species was used to classify the variants with the ACMG guidelines (Tables 4, 5). In addition, the allelic frequencies found here were compared with the calculated Tvs (called “PM2*,” “BA1*,” and “BS1*,” respectively) and the ACMG criteria related to allelic frequency cutoffs (called “PM2,” “BA1,” and “BS1,” respectively). As the assignment to one of the five categories (from pathogenic to benign) depends on the combined evidence of various criteria and as individual criteria have different weights, one additional criterion can sometimes be sufficient to alter the assigned category. As such, the impact of these comparisons on the final classification is illustrated in Table 6.

Table 4

Table 4. Summary of the relevant variant classification criteria and the final category for each of the variants located in the MYBPC3 gene.

Table 5

Table 5. Summary of relevant variant classification criteria for the variants located in the TNNT2, ALMS1, and MYH7 genes.

Table 6

Table 6. Summary of the classification of the variants reported to be associated with hypertrophic cardiomyopathy with and without using the allelic frequency (AF) criteria based on the American College of Medical Genetics and Genomics (ACMG) guidelines or the alternative allelic frequency threshold (Tv), respectively.

For MYBPC3:c.91G > C [A31P], the combination of three pathogenic criteria results in a final classification as pathogenic (Table 4). While the allele frequency of 5.6% would result in an immediate classification as benign based on the allelic frequency criteria from the ACMG (i.e., “BA1”), the frequency is lower than the alternative Tv. As such, the alternative Tv would have led to additional support (“PM2*”) for the classification of pathogenic, while the ACMG criteria conflict with this evaluation (Table 6).

For MYBPC3:c.220G > A [A74T], there was only one pathogenic criterion fulfilled (Table 4). As such, the variant is classified as a variant of unknown significance. While the allele frequency of 22.63% would result in an immediate classification as benign based on the allelic frequency criteria from the ACMG (i.e., “BA1”), the frequency is lower than the alternative Tv (leading to additional “PM2*” support). This would, however, not be sufficient evidence to change the classification to likely pathogenic (Table 6).

For MYBPC3:c.2453C > T [R818W], five pathogenic criteria were fulfilled, resulting in classification as pathogenic (Table 4). With an allele frequency of 2.4%, both the allelic frequency criteria (i.e., “PM2”) and the Tv (i.e., “PM2*”) support the classification as pathogenic (Table 6).

For TNNT2:c.95-108G > A, there were no pathogenic or benign criteria fulfilled, resulting in a classification as a variant of unknown significance (Table 5). With an allele frequency > 27% in the Maine Coon, the allelic frequency criteria of the ACMG alone (i.e., “BA1”) would immediately result in a classification as benign, while the Tv would add additional evidence for benign (i.e., “BS1*”) but this would be insufficient to change the classification to likely benign.

For ALMS1:c.7384G > C [G2462R], only one pathogenic criterion was fulfilled, resulting in a classification as a variant of unknown significance (Table 5). With an allele frequency of close to 52%, the allelic frequency criteria of the ACMG alone (i.e., “BA1”) would immediately result in a classification as benign, while the Tv would add additional evidence for benign (i.e., “BS1*”) but this would be insufficient to change the classification to likely benign.

For MYH7:c.5647G > A [E1883K], three pathogenic criteria were fulfilled, resulting in classification as likely pathogenic (Table 5). With an allele frequency of 0%, both the allelic frequency criteria of the ACMG (“PM2”) and the Tv (“PM2*”) would support the classification as likely pathogenic.

4 Discussion

As a general trend, tests based on genetic variants are often quickly commercially available after their discovery and used for breeding decisions. Compared to phenotypic testing, genetic tests have several advantages, but only when the association these tests are based on is valid. As several variants for HCM were only recently discovered, and as there is a history of spurious associations [i.e., for MYBPC3:c.220G > A [A74T] and TNNT2:c.95-108G > A, the initial associations have not been replicated in subsequent studies (33, 36)], our goal was to classify all six currently published variants based on the ACMG guidelines.

To do this, this study focused on cats of five breeds: the Maine Coon, Ragdoll, Sphynx, British Short- and Longhair, and the Devon Rex. The reasons were that (1) they represent the breed in which a variant was originally discovered (the first three) and/or, (2) because a breed predisposition for HCM has been reported (the first three and the British Short- and Longhair) and/or, (3) because it has been crossbred with the Sphynx (the Devon Rex, British Short- and Longhair). While there are other breeds predisposed to develop HCM, it was decided upfront to only include those breeds in which the sample size requirements were met to allow the association analysis to be performed with sufficient power. An evaluation of the various contributors of the samples shows that the origins of the samples are diverse: the largest number of the samples used in the allelic frequency part of the study originates from laboratories that commercially offer DNA testing, while others originate from a clinical work-up, which can involve screening, diagnostics, and so on. It is difficult to assess whether this might have biased the estimates: samples obtained in a clinical setting can lead to a higher prevalence of disease (43). However, also samples obtained from a genetic testing facility can be biased compared to samples obtained randomly from the population (44). From a practical point of view, it is difficult to truly randomly collect samples. Obtaining sample sizes similar to this study would have required an enormous effort to reach out to groups of cats throughout the world. Furthermore, even then, it cannot be excluded that owners interested in HCM are predominantly submitting samples. Also ethically, the approach is not ideal as it would likely have led to additional invasive sampling, at least if sampling would not be restricted to swabs. Overall, the current method of sampling seemed to be the most convenient and allowed us to collect a sample size of well over 10.000 cats. Furthermore, as discussed further, the allele frequency estimates are in agreement with other studies, which indicates that if there was a bias, it did not result in unexpected deviations.

The first three phases of this study dealt with determining allelic frequencies, ORs, and in silico predictions, respectively. As these results were used to subsequently classify the variants, it is important to assess whether our findings are in line with previous reports, if available.

In terms of allelic frequency, the MYBPC3:c.91G > C [A31P] variant was found in the Maine Coon and one Ragdoll. The presence of the variant in a Ragdoll was a new finding compared to what has been reported in the literature. Considering that the sample size was substantial in the current study (and to our knowledge, the highest reported yet) and that it was only observed in a heterozygous state in one Ragdoll, this seems to be an extremely rare case. Recently, this variant was also found in unique cases in additional breeds: one Pixie-Bob, one Siberian, one Scottish Fold, and one Munchkin (26, 45). These rare occurrences may not necessarily indicate a natural prevalence in these breeds. Instead, they could also result from (1) outcrossing with Maine Coon cats, as has been suggested in previous genetic studies (38), or (2) faulty description of the breed. A third option, a genotyping error, was deemed very unlikely as the raw data of the sample were manually checked and looked entirely normal and the genotype result for that variant is perfectly consistent across various replicate assays. Overall, these findings are notable, but the interpretation of their significance should be made cautiously. The allelic frequency in the Maine Coon was 6% in the present study. In the literature, allelic frequencies in the Maine Coon were slightly higher, with values ranging between 11 and 41%. The highest values were, however, found in the oldest studies, while the allelic frequency was lowest in the most recently published study (26, 36, 46–50). We hypothesized that this might have been the result of selection against this variant as it was discovered nearly two decades ago (18). The data from one of the commercial laboratories that contributed to this study were consistent with this hypothesis: their data showed that allelic frequencies decreased substantially from 2008 to 2021 from 19.54 to 0.95% and from 15.63 to 0% for the MYBPC3:c.91G > C [A31P] and the MYBPC3:c.2453C > T [R818W] variant, respectively. This aligns with the temporal trends observed in earlier research and is slightly lower in 2021 than our current estimate (Supplementary Table S6) (36, 46–50). Overall, in agreement with the literature, this is a variant that mainly occurs in the Maine Coon breed, accidentally in other breeds and selection may have reduced its frequency. The MYBPC3:c.220G > A [A74T] variant was found in all five breeds, with allelic frequencies ranging between 22 and 58%. In the literature, it has also been observed in Norwegian Forest, Persian, domestic shorthair, Siberian, and Bengal cats, with allelic frequencies ranging between 17 and 100% (36, 48). Overall, this variant seems widespread and occurs at high frequencies. The MYBPC3:c.2453C > T [R818W] variant was only found in Ragdoll cats in the current study, with an allelic frequency of 2%. While the allelic frequency in the literature also tended to be higher, ranging between 14 and 23% (36, 46), a more recent study showed an allele frequency of 3% in Ragdolls (26). Recently, this variant was also found in the American Bobtail (with an allele frequency of 11.4% in the longhaired variant), one Highlander, One Munchkin, and two RagaMuffins, all breeds in which crosses to Ragdolls possibly have been part of breed development (26). As allelic frequencies over nearly 15 years of selection have declined (Supplementary Table S6), a trend similar to that for MYBPC3:c.91G > C [A31P] is observed, which might be a consequence of breeder selection (33, 36, 46–48). The TNNT2:c.95-108G > A was found in all five breeds, with allelic frequencies ranging between 2 and 27%. In addition to the breeds studied here, the variant has also been found previously in Persian, Siamese, Tennessee Rex, and Thai breeds, as well as domestic shorthair cats and allelic frequencies were very similar to this study, ranging between 1 and 25% (33). Hence, this is a widespread variant, occurring at high allelic frequencies. MYH7:c.5647G > A [E1883K] was only found in one domestic shorthair cat in the initial report, and the variant has not been found here, nor in a different study of 103 HCM-affected cats (47). Here, no domestic shorthairs were included, so as a result, we can only state that the variant does not seem to segregate in the breeds included in this study. Overall, this seems to be an example of a rare variant. The ALMS1:c.7384G > C [G2462R] variant was observed here in all five breeds studied, at allelic frequencies ranging between 1 and 52%. In the original study, this variant was also reported to segregate in two non-Sphynx cats, which may have belonged to any of 14 breeds included in that study (20). Four of these 14 breeds are also present in this study. Recently, the variant has been reported in the Sphynx at allelic frequencies >50% and also reported in the Scottish Fold, the Exotic Shorthair, and the Minuet breeds, at allelic frequencies approximately 8–9% (28, 45). Hence, this is a widespread variant occurring at high allelic frequencies.

Overall, in terms of allelic frequencies and breeds in which the variant is reported, our data seem to be consistent with the literature and variant characteristics seem to range from rare in terms of number of breeds (MYBPC3:c.2453C > T [R818W], MYBPC3:c.91G > C [A31P]), allelic frequency, or both (MYH7:c.5647G > A [E1883K]) to extremely widespread and prevalent (MYBPC3:c.220G > A [A74T], TNNT2:c.95-108G > A and ALMS1:c.7384G > C [G2462R]).

In our study, a critical aspect was the evaluation of ORs. This is a fundamental tool for the quantification of the association between variants and phenotypes. Nevertheless, various challenges presented, stemming from methodological considerations, sample size constraints, and inconsistencies in the literature. First of all, various calculation methods exist for ORs: one can start from the alleles themselves, or start with the genotypes and group them according to a prespecified inheritance model (51). Instead of choosing one method, we opted for performing several. The reason is that inheritance modes were not always clear in the respective original study or were sometimes disputed in follow-up studies. The downside of this strategy is the multiple testing issue, i.e., the type I error rate or number of false positives increases. As correcting for this negatively affects the power, we chose not to perform a correction because one potential criticism might be that by deliberately increasing the number of tests, we missed out on correct associations. Second, due to the high allelic frequency of several variants (especially MYBPC3:c.220G > A [A74T], ALMS1:c.7384G > C [G2462R], and TNNT2:c.95-108G > A), the sample size quickly decreased when only individuals who were wild type for the other important variants in that breed were retained. Our initial sample size was sufficiently high, so this was not a pressing issue in the current study. However, as the number of potentially associated variants increases in future studies—especially if they are as common as the ones mentioned—securing an adequate sample size will likely become more challenging. A final remark is that a comparison with literature was difficult, as ORs were not always reported and if they were reported, it was not always specified how they were calculated. Nevertheless, for several variants, we were able to compare our results with previous reports in the literature.

For MYBPC3:c.91G > C [A31P], the initial allelic OR in one study was 62.14 (18). The values of subsequent studies tended to be lower. One study reported an allelic OR of 5.10 and a genotypic OR of 19.40 (based on an autosomal recessive mode of inheritance), and a more recent study reported an OR of 8.5 (Vt/Vt versus Wt/Wt) (36, 52). Our results are similar to the subsequent studies with allelic and genotypic ORs of 3.59 and 13.09, respectively, and an OR of 11.25, if calculated as described in the study of Stern et al. (52). For MYBPC3:c.220G > A [A74T], the initial study reported a genotypic OR (based on an autosomal recessive mode of inheritance) of 7.6 (53), while a follow-up study found a value of 0.89 and the latter also resembles our allelic OR result of 2.56 (36). For MYBPC3:c.2453C > T [R818W], it is more difficult to compare the results with other replication studies: sample sizes were sometimes too low to assess association or data presentation did not allow exact OR calculation. Therefore, we were only able to calculate the range in which the allelic OR would have fallen in the original study and that is between 2.56 and 3.3, which means that the allelic OR of 35.82 observed in the current study indicates a larger effect (19, 36, 54). For TNNT2:c.95-108G > A, no OR was reported in the original publication, but based on the data, an allelic OR of 5 was found (22). In a subsequent study, the allelic OR and genotypic OR (based on an autosomal recessive mode of inheritance) were 1.2 and 2.43, which is similar to our results of 1.84 and 2.17, respectively (33). For MYH7:c.5647G > A [E1883K], no ORs were calculated as the variant was never observed, which is also according to expectations as it was postulated to be an extremely rare variant (21). Finally, for ALMS1:c.7384G > C [G2462R], only a relative risk of 13.6 and not an OR was reported in the original paper. However, this relative risk could not be reproduced based on the data that were provided. One subsequent study was not able to find an association between the variant and HCM development (45). This is in line with our result of an allelic and genotypic OR of 1.02 and 1.24 (autosomal recessive mode of inheritance).

Based on these results, three main observations can be made. The first is that the magnitude of the effect often drops between the first “discovery” study and subsequent validation studies. This is a finding that is regularly reported, underscoring the importance of conducting replication studies over time and across different geographic populations (55). The second observation is that our results are generally similar to the results obtained in other validation studies. This is important because this also implies that the subsequent classification was based on replicable results. The third and probably most important observation is that significant associations are reported for two variants only (MYBPC3:c.91G > C [A31P], MYBPC3:c.2453C > T [R818W]), while no association with HCM was reported for three others (MYBPC3:c.220G > A [A74T], TNNT2:c.95-108G > A, and ALMS1:c.7384G > C [G2462R]).

As there are no standardized criteria in veterinary genetics to evaluate whether there is sufficient evidence to consider a variant disease-causing, a decision on when the threshold is reached remains subjective. In human medicine, the ACMG guidelines were developed for that purpose, i.e., to help with the interpretation of newly discovered variants and to standardize their classification. Based on those criteria, only DNA tests of variants that fall in the category “pathogenic” should be used for clinical decisions (23). Variants that are “likely pathogenic” are indicative but should be combined with other evidence (23). For variants of unknown significance, the consensus is that they should not influence clinical decision-making but that the focus should be on resolving their unclear status (23).

Considering the high prevalence of HCM and as several variants were recently reported to be associated with HCM, we aimed to objectively evaluate the current evidence for the six variants published, based on the human ACMG guidelines, with one modification. This modification is linked to the use of the allelic frequency criteria that are tailored to the human situation. Even in humans, these cutoffs are sometimes too low, especially for highly prevalent diseases (23). As the population history of cats (and dogs) is characterized by population bottlenecks, the popular sire effect, and inbreeding, the situation is even more extreme in those species. Associated with that, the allelic frequencies of the disease-causing variants linked to these diseases are also higher. As a consequence, the criteria in the ACMG guidelines have reportedly been deemed too strict (34). For that reason, these criteria were not used in the initial classification but separately considered and compared to an alternative approach that takes disease prevalence into account.

Practically, based on the (modified) ACMG guidelines, two variants were classified as pathogenic, namely, MYBPC3:c.91G > C [A31P] and MYBPC3:c.2453C > T [R818W]. MYH7:c.5647G > A [E1883K] was classified as likely pathogenic. The other three (MYBPC3:c.220G > A [A74T], TNNT2:c.95-108G > A, and ALMS1:c.7384G > C [G2462R]) were all classified as variants of unknown significance. For four (MYBPC3:c.91G > C [A31P], MYBPC3:c.2453C > T [R818W], MYBPC3:c.220G > A [A74T], TNNT2:c.95-108G > A) of these six variants, this finding is a confirmation of earlier studies. In more detail, both the MYBPC3:c.91G > C [A31P] variant in the Maine Coon and the MYBPC3:c.2453C > T [R818W] in the Ragdoll have already been studied intensively, and while there is some disagreement on their actual mode of inheritance, their association with HCM is generally consistent (10, 36, 46–48, 52). For the other two (MYBPC3:c.220G > A [A74T], TNNT2:c.95-108G > A), several publications have discussed (and disproved) the initial association, and the MYBPC3:c.220G > A [A74T] variant was never published by the group that originally reported it (33, 36, 48).

As MYH7:c.5647G > A [E1883K] and ALMS1:c.7384G > C [G2462R] were only recently published, independent studies re-evaluating their association are limited. For the latter of the two (ALMS1:c.7384G > C [G2462R]), this study and another recently published study could not confirm the association with HCM (45). As too few other criteria supporting the classification are fulfilled, there is insufficient evidence to consider this variant as pathogenic or likely pathogenic; hence, we conclude that it is a variant of unknown significance. For the other variant (MYH7:c.5647G > A [E1883K]), the ACMG guidelines were already used for classification in the original publication. Recently, an independent study confirmed this low allelic frequency (47), and also in our study, this variant was absent from the entire cat population studied, albeit no domestic shorthairs were included here. In humans, the same variant has been described (ClinVar accession ID: VCV000014121.13), and based on a set of criteria that partially overlaps, the variant is there considered a variant of unknown significance. An in-depth comparison showed an agreement on the fulfillment of PP3 and PM2 (i.e., the computational evidence and that it is a rare variant) and that PM1 was not fulfilled (i.e., the location of the mutational hotspot ranged between amino acids 181 and 937, 51). The PP2 criterion (i.e., missense variants are a common mechanism) is not present anymore in the guidelines that resulted in the ClinVar classification, although they do state that it is “the predominant class of pathogenic alleles in MYH7” (42). They have an additional criterion that is fulfilled, namely, PS4, but with the weight set to support based on the observation that the allele occurred in five individuals consistent with the expected phenotype (i.e., these five individuals had HCM) (17, 56). In the original ACMG guidelines that were used here, this option is not available. The only disagreement was on PS3, a criterion that we considered fulfilled, and they did not. The decision to consider it fulfilled here was based on four papers: three of which did similar analyses and all showed a structural effect (38, 40, 56). The fourth paper did a movement analysis based on C. elegans and did not find an effect for any of the variants scrutinized, which makes it difficult to assess the suitability of the method itself as the other variants that were tested in that same paper were still considered to be pathogenic even with that result (41). Based on the overall agreement on whether criteria are fulfilled or not (i.e., there is only one that is debated, three that are agreed on, and two that have no direct match), the classification as “likely pathogenic” or “variant of unknown significance” seems linked to the system that is used. Among the authors, it did, however, raise various interesting questions and sparked debate: should the data of humans and animals be jointly considered as evidence? If the data are combined, there are six cases (five in humans, one in a cat), which lead to an increased weight of the PS4 criterion in the ClinVar classification and the fulfillment of the PS4 criterion in the ACMG guidelines applied in our study. However, this raises a question: should the PS4 criterion be applied species-specifically? To what extent can data be extrapolated across species either way? How should extremely rare variants be dealt with? While these questions are beyond the scope of the current study, they present interesting avenues for future studies. Overall, in the context of screening programs, we do not consider it necessary to include the MYH7:c.5647G > A [E1883K] variant as it is extremely rare (10).

Aside from the evaluation of these six variants with the guidelines, an additional goal was to compare the allelic frequency criteria of the ACMG and the alternative Tv-method. This can be done at two levels. The first level is checking whether the fulfillment of the criterion agrees with the classification based on all other criteria. For the Tv-method, this means that if the observed sample allele frequency is lower or equal to the cutoff, the frequency is consistent with the disease characteristics. If it is higher, that implies an inconsistency with theoretical expectations and that is a warning sign and indicates a benign role. For the ACMG, the cutoffs do not change with disease prevalence, penetrance, or genetic heterogeneity: an allelic frequency of >5% indicates a benign role, and (extremely) low allelic frequencies support a classification as pathogenic. Based on the Tv threshold, the allelic frequency of two out of two pathogenic, one out of the one classified as likely pathogenic, and two of the three variants classified as variants of unknown significance is correct. Based on the allelic frequency criteria of the ACMG guidelines, the allelic frequency of one out of two pathogenic, one out of one likely pathogenic, and three out of three variants classified as variants of unknown significance is correct; hence, the total agreement of the two methods at this level is identical (i.e., an agreement for five out of six variants).

The second level involves an evaluation of the effect on the final classification if the criterion would have been used as it is normally used. This is linked to the assigned weight of the criterion in the decision-making process. For the Tv threshold, the classification never changes: the two pathogenic variants remain pathogenic, the one likely pathogenic remains likely pathogenic, and the three others remain variants of unknown significance. For the ACMG guidelines, one pathogenic variant becomes a benign variant, and all three variants of unknown significance also become benign variants. Practically, this means that in the ACMG guidelines, this one criterion is enough to overrule all other evidence for 4/6 variants and that simply incorporating this criterion would change the classification drastically. Based on that result, the weight of the ACMG allelic frequency criteria seems to be too dominant in the decision-making process. However, based on the observation that even for highly prevalent diseases in humans, these cutoffs are too strict and considering that high disease prevalence (and the associated high allelic frequencies) occurs even far more often in domestic animals due to the breeding practices often used in those species, we doubt that the fixed ACMG cutoffs themselves are actually that suitable, in agreement with an earlier report (34). Broader than for HCM alone, as a wide range of variants has high allelic frequencies in cats (and dogs), if the BA1 stand-alone criterion would be applied to all of them, a large proportion would be declassified from (likely) pathogenic, irrespective of all the other information.

Interesting and similar to previous reports, the autosomal dominant mode of inheritance that was proposed for MYBPC3:c.91G > C [A31P] for Maine Coon cats in the original publication is less of an optimal fit to the data (Table 2). The autosomal recessive mode of inheritance seems a better fit to the data and is thus more likely (10, 36, 52). While an autosomal dominant mode of inheritance might make sense biochemically, statistically reports do not seem to find an increased risk for heterozygous cats, which is something one would expect if that is the correct mode of inheritance (10, 36). The difference between autosomal recessive and autosomal dominant might at some point become a semantic discussion (36); however, when informing owners on the risk, they must realize that the risk of developing HCM does not seem to be increased in heterozygous cats based on our results. Taking that into account, altering breeding advice for that mutation to an approach similar to what would be done for autosomal recessive diseases, especially when genetic diversity is low(er), might also be considered (24, 25). Noteworthy, for MYBPC3:c.2453C > T [R818W], the distribution of genotypes across cases and controls made it impossible to discriminate an autosomal recessive or dominant mode of inheritance. Hence, there is no evidence against the originally proposed autosomal dominant mode of inheritance (19).

4.1 Limitations of the study

The samples collected in this study were derived from two main sources. The added value of the samples of the commercial laboratories is clear: large sample sizes make by chance deviations less probable and increase the probability to detect rare events. There are, however, two downsides: breeds are reported by the owners themselves for these samples, which makes it of course impossible to check whether breed was correctly specified. Furthermore, while the potential presence of duplicate samples from batch one could be ruled out, an absence of overlap between the samples from the commercial laboratories could not entirely be guaranteed. Both might have had an influence on the allelic frequencies of the MYBPC3:c.91G > C [A31P] and the MYBPC3:c.2453C > T [R818W] variants. However, as none of the results related to allelic frequencies was unexpected, we deem it unlikely that this has influenced the results drastically. A second limitation is linked to the international collaboration that has increased the sample size in this study: this might have led to variation in phenotypical classification. In more detail, several veterinarians were involved, and it has been shown that the interpretation of cardiac ultrasound measurements can differ between observers (57). However, an attempt at standardization was made by using internationally accepted consensus guidelines for echocardiographic phenotypic classification (29). Furthermore, those cats for which the phenotype was not entirely clear (=the equivocal cats) were excluded from the association part of the study. As the median age of the healthy control group was very high, this makes it far less likely to have a substantial proportion of false negative cats (i.e., cats that still will develop HCM, but did not have it at the time of the screening) in the control group. Based on all these combinations, the probability of misclassification was reduced as much as possible. A third limitation is linked to the absence of domestic short- and longhair cats from the study population. It can be argued that this absence might have resulted in missing the MYH7:c.5647G > A [E1883K] variant. We deem this unlikely, however: both in the original publication and in a recently published one, a large set of domestic shorthairs were included, leading to a combined sample size of 196 cats (21, 47). As the variant was never found in either of the publications, it seems very rare and the current data only support this with the observed absence in a large dataset of several breeds, including one breed (=the Sphynx) that had not been tested yet. A final limitation is that while we considered several modes of inheritance, these are only a subset of all the modes that are possible, and this subset might not fully capture the intricacies of HCM’s genetic landscape.

5 Conclusion

Based on the separate phases in this study, the evaluation with the ACMG guidelines is in line with previous research, two variants were found to be pathogenic, and one variant was classified as likely pathogenic. However, practically we recommend limiting routine genetic testing to two variants in two breeds, i.e., the MYBPC3:c.91G > C [A31P] in Maine Coon and the MYBPC3:c.2453C > T [R818W] variant in Ragdolls (10). While their (homozygous) presence increases the risk of developing HCM in those two breeds, they are virtually absent in other major breeds: the recommendation to specifically screen for these variants in the breeds in which they were originally described, but not other breeds, thus remains (36, 58, 59). While the MYH7:c.5647G > A [E1883K] variant is categorized as likely pathogenic, its rare occurrence renders routine testing for this variant unnecessary (10). The other variants lack evidence for disease causation, leading to a consensus against screening for these variants in any breed (36, 58, 59). Overall, the human ACMG guidelines are a good starting point, but one modification has already been shown to perform at least as well as the traditional allelic frequency criteria. Furthermore, interesting questions were raised during the application of these guidelines. As such, a thorough and evidence-based evaluation (and potential modification) of the ACMG guidelines would be an interesting next step to optimize their applicability in veterinary medicine.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

All included cats had informed owner consent. Samples were non-invasive buccal swabs or whole blood on EDTA. In some cats, blood was obtained as part of routine clinical procedures for diagnostic purposes, at the request and with the consent of the owner. As these samples were from client-owned cats for which no harmful invasive procedures were performed, there was no animal experimentation according to the legal definition in Europe (Subject 5f of Article1, Chapter I of the Directive 2010/63/UE of the European Parliament and of the Council). A subset of the blood samples was collected as part of a biobanking initiative, approved by the Ethics Committee of the Faculties of Veterinary Medicine and Bioscience Engineering at Ghent University (EC2017-86). Samples from Sweden were collected in agreement with the ethical committee of the Swedish Board of Agriculture [No. C2/12 (2012-02-24),C12/15 (2015-02-27), and Dnr 5.8.18–04682/2020]. Samples sourced from the Cornell Veterinary Biobank were collected under Institutional Animal Care and Use Committee (IACUC) approved animal protocol #2005–0151. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions

FB: Conceptualization, Data curation, Formal analysis, Resources, Writing – original draft, Writing – review & editing. MA: Conceptualization, Resources, Writing – review & editing. HA: Conceptualization, Resources, Writing – review & editing. TD: Conceptualization, Resources, Writing – review & editing. PF: Conceptualization, Resources, Writing – review & editing. PRF: Conceptualization, Resources, Writing – review & editing. JJH: Conceptualization, Resources, Writing – review & editing. JH: Conceptualization, Resources, Writing – review & editing. SD: Conceptualization, Resources, Writing – review & editing. MK: Conceptualization, Resources, Writing – review & editing. FS: Conceptualization, Resources, Writing – review & editing. IL: Conceptualization, Resources, Writing – review & editing. LL: Conceptualization, Resources, Writing – review & editing. ML: Conceptualization, Resources, Writing – review & editing. ÅO: Conceptualization, Resources, Writing – review & editing. LP: Conceptualization, Resources, Writing – review & editing. CC: Conceptualization, Resources, Writing – review & editing. PS: Conceptualization, Resources, Writing – review & editing. MT: Conceptualization, Resources, Writing – review & editing. BB: Conceptualization, Data curation, Formal analysis, Resources, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was partially funded by the Bijzonder Onderzoeksfonds (BOF) starting Grant (01N04119), Maria Naundorf – van Gorkum Foundation, the LOOF (Livre Officiel des Origines Félines), and the French Sphynx Club. Feline DNA samples from VetAgro Sup are part of the Feli-DNA biobank, which is part of the CRB-Anim infrastructure (ANR-11-INBS-003) in the framework of the “Investing for the future” program (PIA). The Cornell Veterinary Biobank is a resource built with the support of NIH grant R24 GM082910 and the Cornell University College of Veterinary Medicine.

Acknowledgments

The authors wish to thank Linda Impe, Ruben Van Gansbeke, Dominique Vander Donckt, and Claudia Rozendom for their excellent technical assistance. The authors would also like to thank the cat owners and breeders for providing sample and clinical reports.

Conflict of interest

HA and SD are current and former employees of Wisdom Panel Mars Petcare Science & Diagnostics that offers canine and feline DNA testing as a commercial service. CC is a current employee of Antagene, a DNA testing and genetic analysis company for dogs, cats, horses, and wildlife. MT is employed by Genefast srl. LL has historically received funding from various commercial testing laboratories to support variant discovery but has no involvement with these HCM variants.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2024.1327081/full#supplementary-material

Footnotes

1. ^http://genetics.bwh.harvard.edu/pph2/

2. ^https://sift.bii.a-star.edu.sg/www/SIFT_seq_submit2.html

3. ^http://mutpred2.mutdb.org/#qform

4. ^http://wangcomputing.com/assp/

5. ^http://hollywood.mit.edu/GENSCAN.html

6. ^https://services.healthtech.dtu.dk/service.php?NetGene2-2.42

References

1. Häggström, J, Andersson, O, Falk, T, Nilsfors, L, Oisson, U, Kresken, JG, et al. Effect of body weight on echocardiographic measurements in 19, 866 pure-bred cats with or without heart disease. J Vet Intern Med. (2016) 30:1601–11. doi: 10.1111/jvim.14569

PubMed Abstract | Crossref Full Text | Google Scholar

2. Riesen, SC, Kovacevic, A, Lombard, CW, and Amberger, C. Echocardiographic screening of purebred cats: an overview from 2002 to 2005. Schweiz Arch Tierheilkd. (2007) 149:73–6. doi: 10.1024/0036-7281.149.2.73

PubMed Abstract | Crossref Full Text | Google Scholar

3. Payne, JR, Brodbelt, DC, and Luis, FV. Cardiomyopathy prevalence in 780 apparently healthy cats in rehoming centres (the CatScan study). J Vet Cardiol. (2015) 17:S244–57. doi: 10.1016/j.jvc.2015.03.008

Crossref Full Text | Google Scholar

4. Fox, PR, Keene, BW, Lamb, K, Schober, KA, Chetboul, V, Luis Fuentes, V, et al. International collaborative study to assess cardiovascular risk and evaluate long-term health in cats with preclinical hypertrophic cardiomyopathy and apparently healthy cats: the REVEAL study. J Vet Intern Med. (2018) 32:930–43. doi: 10.1111/jvim.15122

PubMed Abstract | Crossref Full Text | Google Scholar

5. Semsarian, C, Ingles, J, Maron, MS, and Maron, BJ. New perspectives on the prevalence of hypertrophic cardiomyopathy. J Am Coll Cardiol. (2015) 65:1249–54. doi: 10.1016/j.jacc.2015.01.019

Crossref Full Text | Google Scholar

6. Follby, A, Pettersson, A, Ljungvall, I, Ohlsson, Å, and Häggström, J. A questionnaire survey on long-term outcomes in cats breed-screened for feline cardiomyopathy. Animals. (2022) 12:1–18. doi: 10.3390/ani12202782

Crossref Full Text | Google Scholar

7. Kittleson, MD, and Côté, E. The feline cardiomyopathies: 1 General concepts. J Feline Med Surg. (2021) 23:1009–27. doi: 10.1177/1098612X211021819

Crossref Full Text | Google Scholar

8. Fox, PR, Liu, S-K, and Maron, BJ. Echocardiographic assessment of spontaneously occurring feline hypertrophic cardiomyopathy: an animal model of human disease. Circulation. (1995) 92:2645–51. doi: 10.1161/01.CIR.92.9.2645

PubMed Abstract | Crossref Full Text | Google Scholar

9. Côté, E, MacDonald, KA, Meurs, KM, and Sleeper, MM. Hypertrophic cardiomyopathy In: Feline Cardiology. Chichester: John Wiley & Sons (2011). 103–75.

Google Scholar

10. Schipper, T, Peelman, LJ, Smets, P, and Broeckx, BJG. An overview of the current genetic and phenotypical selection strategies to reduce the prevalence of feline hypertrophic cardiomyopathy. Vlaams Diergeneeskd Tijdschr. (2020) 89:69–80. doi: 10.21825/vdt.v89i2.16355

Crossref Full Text | Google Scholar

11. Maron, BJ, and Fox, PR. Hypertrophic cardiomyopathy in man and cats. J Vet Cardiol. (2015) 17:S6–9. doi: 10.1016/j.jvc.2015.03.007

PubMed Abstract | Crossref Full Text | Google Scholar

12. Hogan, DF, Fox, PR, Jacob, K, Keene, B, Laste, NJ, Rosenthal, S, et al. Secondary prevention of cardiogenic arterial thromboembolism in the CAT: the double-blind, randomized, positive-controlled feline arterial thromboembolism; Clopidogrel vs. aspirin trial (FAT CAT). J Vet Cardiol. (2015) 17:S306–17. doi: 10.1016/j.jvc.2015.10.004

Crossref Full Text | Google Scholar

13. Ueda, Y, and Stern, JA. A one health approach to hypertrophic cardiomyopathy. Yale J Biol Med. (2017) 90:433–48.

PubMed Abstract | Google Scholar

14. Freeman, LM, Rush, JE, Stern, JA, Huggins, GS, and Maron, MS. Feline hypertrophic cardiomyopathy: a spontaneous large animal model of human HCM. Cardiol Res. (2017) 8:139–42. doi: 10.14740/cr578w

PubMed Abstract | Crossref Full Text | Google Scholar

15. Hollman, A, Goodwin, JF, Teare, D, and Renwick, JW. A family with obstructive cardiomyopathy (asymmetrical hypertrophy). Br Heart J. (1960) 22:449–56. doi: 10.1136/hrt.22.4.449

Crossref Full Text | Google Scholar

16. Geisterfer-Lowrance, AAT, Kass, S, Tanigawa, G, Vosberg, HP, McKenna, W, Seidman, CE, et al. A molecular basis for familial hypertrophic cardiomyopathy: a β cardiac myosin heavy chain gene missense mutation. Cell. (1990) 62:999–1006. doi: 10.1016/0092-8674(90)90274-I

PubMed Abstract | Crossref Full Text | Google Scholar

17. Walsh, R, Thomson, KL, Ware, JS, Funke, BH, Woodley, J, McGuire, KJ, et al. Reassessment of Mendelian gene pathogenicity using 7,855 cardiomyopathy cases and 60,706 reference samples. Genet Med. (2017) 19:192–203. doi: 10.1038/gim.2016.90

PubMed Abstract | Crossref Full Text | Google Scholar

18. Meurs, KM, Sanchez, X, David, RM, Bowles, NE, Towbin, JA, Reiser, PJ, et al. A cardiac myosin binding protein C mutation in the Maine coon cat with familial hypertrophic cardiomyopathy. Hum Mol Genet. (2005) 14:3587–93. doi: 10.1093/hmg/ddi386

PubMed Abstract | Crossref Full Text | Google Scholar

19. Meurs, KM, Norgard, MM, Ederer, MM, Hendrix, KP, and Kittleson, MD. A substitution mutation in the myosin binding protein C gene in ragdoll hypertrophic cardiomyopathy. Genomics. (2007) 90:261–4. doi: 10.1016/j.ygeno.2007.04.007

PubMed Abstract | Crossref Full Text | Google Scholar

20. Meurs, KM, Williams, BG, DeProspero, D, Friedenberg, SG, Malarkey, DE, Ezzell, JA, et al. A deleterious mutation in the ALMS1 gene in a naturally occurring model of hypertrophic cardiomyopathy in the sphynx cat. Orphanet J Rare Dis. (2021) 16:108–7. doi: 10.1186/s13023-021-01740-5

PubMed Abstract | Crossref Full Text | Google Scholar

21. Schipper, T, Van Poucke, M, Sonck, L, Smets, P, Ducatelle, R, Broeckx, BJG, et al. A feline orthologue of the human MYH7 c. 5647G > a (p. (Glu1883Lys)) variant causes hypertrophic cardiomyopathy in a domestic shorthair cat. Eur J Hum Genet [Internet]. (2019) 27:1724–1730. doi: 10.1038/s41431-019-0431-4

Crossref Full Text | Google Scholar

22. McNamara, JW, Schuckman, M, Becker, RC, and Sadayappan, S. A novel homozygous Intronic variant in TNNT2 associates with feline cardiomyopathy. Front Physiol. (2020) 11:1–7. doi: 10.3389/fphys.2020.608473

Crossref Full Text | Google Scholar

23. Richards, S, Aziz, N, Bale, S, Bick, D, Das, S, Gastier-Foster, J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. (2015) 17:405–24. doi: 10.1038/gim.2015.30

PubMed Abstract | Crossref Full Text | Google Scholar

24. Broeckx, BJG. The dog 2.0: lessons learned from the past. Theriogenology. (2020) 150:20–6. doi: 10.1016/j.theriogenology.2020.01.043

PubMed Abstract | Crossref Full Text | Google Scholar

25. Moses, L, Niemi, S, and Karlsson, E. Pet genomics medicine runs wild. Nature. (2018) 559:470–2. doi: 10.1038/d41586-018-05771-0

PubMed Abstract | Crossref Full Text | Google Scholar

26. Anderson, H, Davison, S, Lytle, KM, Honkanen, L, Freyer, J, Mathlin, J, et al. Genetic epidemiology of blood type, disease and trait variants, and genome-wide genetic diversity in over 11,000 domestic cats. PLoS Genet. (2022) 18:1–30. doi: 10.1371/journal.pgen.1009804

Crossref Full Text | Google Scholar

27. Stevenson, M, Nunes, T, Heuer, C, Marshall, J, Sanchez, J, Thomton, R, et al. (eds.) ‘epiR’: tools for the analysis of epidemiological data, R package version 0.9-79. Available at: https://cran.r-project.org/web/packages/epiR/index.html and https://cran.r-project.org/web/packages/epiR/epiR.pdf (2016)

Google Scholar

28. Turba, ME, Ferrari, P, Milanesi, R, Gentilini, F, and Longeri, M. HCM-associated ALMS1 variant: allele drop-out and frequency in Italian sphynx cats. Anim Genet. (2023)54:1–4. doi: 10.1111/age.13340

Crossref Full Text | Google Scholar

29. Luis Fuentes, V, Abbott, J, Chetboul, V, Côté, E, Fox, PR, Häggström, J, et al. ACVIM consensus statement guidelines for the classification, diagnosis, and management of cardiomyopathies in cats. J Vet Intern Med. (2020) 34:1062–77. doi: 10.1111/jvim.15745

Crossref Full Text | Google Scholar

30. Ferasin, L, Sturgess, CP, Cannon, MJ, Caney, SMA, Gruffydd-Jones, TJ, and Wotton, PR. Feline idiopathic cardiomyopathy: a retrospective study of 106 cats (1994-2001). J Feline Med Surg. (2003) 5:151–9. doi: 10.1016/S1098-612X(02)00133-X

Crossref Full Text | Google Scholar

31. Ruxton, GD, and Neuhäuser, M. Review of alternative approaches to calculation of a confidence interval for the odds ratio of a 2 × 2 contingency table. Methods Ecol Evol. (2013) 4:9–13. doi: 10.1111/j.2041-210x.2012.00250.x

Crossref Full Text | Google Scholar

32. Beckers, E, Cornelis, I, Bhatti, SFM, Smets, P, Shelton, GD, Guo, LT, et al. A nonsense variant in the DMD gene causes X-linked muscular dystrophy in the Maine coon cat. Animals. (2022) 12:1–12. doi: 10.3390/ani12212928

Crossref Full Text | Google Scholar

33. Schipper, T, Ohlsson, Å, Longeri, M, Hayward, JJ, Mouttham, L, Ferrari, P, et al. TheTNNT2:c.95‐108G>a variant is common in Maine coons and shows no association with hypertrophic cardiomyopathy. Anim Genet. (2022) 53:526–9. doi: 10.1111/age.13223

Crossref Full Text | Google Scholar

34. Broeckx, BJG, Peelman, L, Saunders, JH, Deforce, D, and Clement, L. Using variant databases for variant prioritization and to detect erroneous genotype-phenotype associations. BMC Bioinformatics. (2017) 18:1–10. doi: 10.1186/s12859-017-1951-y

Crossref Full Text | Google Scholar

35. Tudurachi, B, Alexandra, Z, Leonte, A, Laura, T, Ureche, C, Bîrgoan, SG, et al. An update on MYBPC3 gene mutation in hypertrophic cardiomyopathy. Int J Mol Sci. (2023) 24:10510. doi: 10.3390/ijms241310510

Crossref Full Text | Google Scholar

36. Longeri, M, Ferrari, P, Knafelz, P, Mezzelani, A, Marabotti, A, Milanesi, L, et al. Myosin-binding protein C DNA variants in domestic cats (A31P, A74T, R820W) and their association with hypertrophic cardiomyopathy. J Vet Intern Med. (2013) 27:275–85. doi: 10.1111/jvim.12031

PubMed Abstract | Crossref Full Text | Google Scholar

37. Vera, TR, Iglesias, LM, Prieto, MH, Ortiz, M, Garcia, IR, Gallizo, NG, et al. The R820W mutation in the MYBPC3 gene, associated with hypertrophic cardiomyopathy in cats, causes hypertrophic cardiomyopathy and left ventricular non-compaction in humans. Int J Cardiol. (2010) 145:405–7. doi: 10.1016/j.ijcard.2010.04.032

Crossref Full Text | Google Scholar

38. Sohn, RL, Vikstrom, KL, Strauss, M, Cohen, C, Szent-Gyorgyi, A, and Leinwand, LA. A 29 residue region of the Sarcomeric Myosin rod is necessary for filament formation. J Mol Biol. (1997) 266:317–30. doi: 10.1006/jmbi.1996.0790

Crossref Full Text | Google Scholar

39. Armel, TZ, and Leinwand, LA. Mutations in the β-myosin rod cause myosin storage myopathy via multiple mechanisms. Proc Natl Acad Sci U S A. (2009) 106:6291–6. doi: 10.1073/pnas.0900107106

Crossref Full Text | Google Scholar

40. Viswanathan, MC, Tham, RC, Kronert, WA, Sarsoza, F, Trujillo, AS, Cammarato, A, et al. Myosin storage myopathy mutations yield defective myosin filament assembly in vitro and disrupted myofibrillar structure and function in vivo. Hum Mol Genet. (2017) 26:4799–813. doi: 10.1093/hmg/ddx359

PubMed Abstract | Crossref Full Text | Google Scholar

41. Dahl-Halvarsson, M, Pokrzywa, M, Rauthan, M, Pilon, M, and Tajsharghi, H. Myosin storage myopathy in C. elegans and human cultured muscle cells. PLoS One. (2017) 12:1–14. doi: 10.1371/journal.pone.0170613

Crossref Full Text | Google Scholar

42. Kelly, MA, Caleshu, C, Morales, A, Buchan, J, Wolf, Z, Harrison, SM, et al. Adaptation and validation of the ACMG/AMP variant classification framework for MYH7-associated inherited cardiomyopathies: recommendations by ClinGen’s inherited cardiomyopathy expert panel. Genet Med. (2018) 20:351–9. doi: 10.1038/gim.2017.218

PubMed Abstract | Crossref Full Text | Google Scholar

43. Broeckx, BJG, Verhoeven, G, Coopman, F, Van Haeringen, W, Bosmans, T, Gielen, I, et al. The effects of positioning, reason for screening and the referring veterinarian on prevalence estimates of canine hip dysplasia. Vet J. (2014) 201:378–84. doi: 10.1016/j.tvjl.2014.05.023

PubMed Abstract | Crossref Full Text | Google Scholar

44. Beckers, E, Casselman, I, Soudant, E, Daminet, S, Paepe, D, Peelman, L, et al. The prevalence of the ABCB1-1Δ variant in a clinical veterinary setting: the risk of not genotyping. PLoS One. (2022) 17:e0273706. doi: 10.1371/journal.pone.0273706

PubMed Abstract | Crossref Full Text | Google Scholar

45. Akiyama, N, Suzuki, R, Saito, T, Yuchi, Y, Ukawa, H, and Matsumoto, Y. Presence of known feline ALMS1 and MYBPC3 variants in a diverse cohort of cats with hypertrophic cardiomyopathy in Japan. PLoS One. (2023) 18:e0283433. doi: 10.1371/journal.pone.0283433

PubMed Abstract | Crossref Full Text | Google Scholar

46. Casamian-Sorrosal, D, Chong, SK, Fonfara, S, and Helps, C. Prevalence and demographics of the MYBPC3-mutations in ragdolls and Maine coons in the British Isles. J Small Anim Pract. (2014) 55:269–73. doi: 10.1111/jsap.12201

PubMed Abstract | Crossref Full Text | Google Scholar

47. O’Donnell, K, Adin, D, Atkins, CE, DeFrancesco, T, Keene, BW, Tou, S, et al. Absence of known feline MYH7 and MYBPC3 variants in a diverse cohort of cats with hypertrophic cardiomyopathy. Anim Genet. (2021) 52:542–4. doi: 10.1111/age.13074

PubMed Abstract | Crossref Full Text | Google Scholar

48. Wess, G, Schinner, C, Weber, K, Küchenhoff, H, and Hartmann, K. Association of A31P and A74T polymorphisms in the myosin binding protein C3 gene and hypertrophic cardiomyopahty in Maine Coon and other breed cats. J Vet Intern Med. (2010) 24:527–32. doi: 10.1111/j.1939-1676.2010.0514.x

PubMed Abstract | Crossref Full Text | Google Scholar

49. Mary, J, Chetboul, V, Sampedrano, CC, Abitbol, M, Gouni, V, Trehiou-Sechi, E, et al. Prevalence of the MYBPC3-A31P mutation in a large European feline population and association with hypertrophic cardiomyopathy in the Maine Coon breed. J Vet Cardiol. (2010) 12:155–61. doi: 10.1016/j.jvc.2010.06.004

PubMed Abstract | Crossref Full Text | Google Scholar

50. Sukumolanan, P, and Petchdee, S. Prevalence of cardiac myosin-binding protein C3 mutations in Maine coon cats with hypertrophic cardiomyopathy. Vet World. (2022) 15:502–8. doi: 10.14202/vetworld.2022.502-508

PubMed Abstract | Crossref Full Text | Google Scholar

51. Bush, WS, and Moore, JH. Genome-wide association studies. PLoS Comput Biol. (2012) 8:e1002822. doi: 10.1371/journal.pcbi.1002822

Crossref Full Text | Google Scholar

52. Stern, JA, Rivas, VN, Kaplan, JL, Ueda, Y, Oldach, MS, Ontiveros, ES, et al. Hypertrophic cardiomyopathy in purpose - bred cats with the A31P mutation in cardiac myosin binding protein - C. Sci Rep. (2023) 13:1–13. doi: 10.1038/s41598-023-36932-5

PubMed Abstract | Crossref Full Text | Google Scholar

53. Nyberg, M, Koch, J, and Christiansen, M. Intra-allelic genetic heterogeneity of hypertrophic cardiomyopathy in the Maine Coon cat. In Proceedings of the 12th human genome meeting. Montreal, Canada (2007)

Google Scholar

54. Borgeat, K, Stern, J, Meurs, KM, Fuentes, VL, and Connolly, DJ. The influence of clinical and genetic factors on left ventricular wall thickness in ragdoll cats. J Vet Cardiol. (2015) 17:S258–67. doi: 10.1016/j.jvc.2015.06.005

PubMed Abstract | Crossref Full Text | Google Scholar

55. König, IR. Validation in genetic association studies. Brief Bioinform. (2011) 12:253–8. doi: 10.1093/bib/bbq074

Crossref Full Text | Google Scholar

56. Tajsharghi, H, Oldfors, A, Macleod, DP, and Swash, M. Homozygous mutation in MYH7 in myosin storage myopathy and cardiomyopathy. Neurology. (2007) 68:962–3. doi: 10.1212/01.wnl.0000257131.13438.2c

PubMed Abstract | Crossref Full Text | Google Scholar

57. van Hoek, I, Payne, JR, Feugier, A, and Connolly, DJ. Inter-observer variability for cardiac ultrasound measurements in cats repeated at different time points in early adult life. Vet Anim Sci. (2017) 5:44–6. doi: 10.1016/j.vas.2018.03.002

PubMed Abstract | Crossref Full Text | Google Scholar

58. Côté, E, MacDonald, KA, Meurs, KM, and Sleeper, MM. (Eds.). Restrictive/unclassified cardiomyopathy In: Feline Cardiology. Chichester: John Wiley & Sons (2012). 177–81.

Google Scholar

59. Lyons, LA. DNA mutations of the cat: the good, the bad and the ugly. J Feline Med Surg. (2015) 17:203–19. doi: 10.1177/1098612X15571878

Crossref Full Text | Google Scholar