Osteochondrodysplasias are a heterogeneous group of genetic skeletal dysplasias. Mutations in the COL2A1 gene cause a spectrum of rare autosomal-dominant type II collagenopathies characterized by skeletal dysplasia, short stature, and with vision and auditory defects. In this study, we have investigated in more detail the phenotypic and genotypic characterization resulting from glycine to serine mutations in the COL2A1 gene in a 2-year-old boy.
Materials and methods
Detailed clinical and radiological phenotypic characterization was the baseline tool to guide the geneticists toward proper genotypic confirmation.
Results
Genetic analysis revealed a de novo mutation, c.1681G>A (p.Gly561Ser), in the collagen type II alpha-1 gene (COL2A1). The identified variant showed impaired protein stability, and lead to dysfunction of type II collagen. In addition to pre and postnatal growth retardation, remarkable retardation of gross motor development and intellectual disability were noted. The latter was connected to cerebral malformations. The overall clinical phenotype of our current patient resembles spondyloepiphyseal dysplasia congenita (SEDC), but with extra phenotypic criteria.
Conclusions
The aim of this paper is twofold; firstly, raising awareness among orthopaedic surgeons when dealing with children manifesting multiple deformities, and secondly to broaden the clinical phenotype in patients with COL2A1 mutations of amino acid substitution (glycine to serine).
Citation: Mohammad Shboul, Hela Sassi, Houweyda Jilani, Imen Rejeb, Yasmina Elaribi, Syrine Hizem, Lamia Ben Jemaa, Marwa Hilmi, Susanna Gerit Kircher, Ali Al Kaissi. The phenotypic spectrum in a patient with Glycine to Serine mutation in the COL2A1 gene: overview study[J]. AIMS Molecular Science, 2021, 8(1): 76-85. doi: 10.3934/molsci.2021006
Related Papers:
[1]
María Ángeles García-Ferrero, Angkana Rüland .
Strong unique continuation for the higher order fractional Laplacian. Mathematics in Engineering, 2019, 1(4): 715-774.
doi: 10.3934/mine.2019.4.715
[2]
Chiara Gavioli, Pavel Krejčí .
Deformable porous media with degenerate hysteresis in gravity field. Mathematics in Engineering, 2025, 7(1): 35-60.
doi: 10.3934/mine.2025003
[3]
Fernando Farroni, Gioconda Moscariello, Gabriella Zecca .
Lewy-Stampacchia inequality for noncoercive parabolic obstacle problems. Mathematics in Engineering, 2023, 5(4): 1-23.
doi: 10.3934/mine.2023071
[4]
Federico Cluni, Vittorio Gusella, Dimitri Mugnai, Edoardo Proietti Lippi, Patrizia Pucci .
A mixed operator approach to peridynamics. Mathematics in Engineering, 2023, 5(5): 1-22.
doi: 10.3934/mine.2023082
[5]
Francesco Maddalena, Danilo Percivale, Franco Tomarelli .
Signorini problem as a variational limit of obstacle problems in nonlinear elasticity. Mathematics in Engineering, 2024, 6(2): 261-304.
doi: 10.3934/mine.2024012
[6]
Xavier Fernández-Real, Alessio Figalli .
On the obstacle problem for the 1D wave equation. Mathematics in Engineering, 2020, 2(4): 584-597.
doi: 10.3934/mine.2020026
[7]
Patrizia Pucci, Letizia Temperini .
On the concentration–compactness principle for Folland–Stein spaces and for fractional horizontal Sobolev spaces. Mathematics in Engineering, 2023, 5(1): 1-21.
doi: 10.3934/mine.2023007
[8]
Petteri Harjulehto, Peter Hästö, Jonne Juusti .
Bloch estimates in non-doubling generalized Orlicz spaces. Mathematics in Engineering, 2023, 5(3): 1-21.
doi: 10.3934/mine.2023052
[9]
Luca Spolaor, Bozhidar Velichkov .
On the logarithmic epiperimetric inequality for the obstacle problem. Mathematics in Engineering, 2021, 3(1): 1-42.
doi: 10.3934/mine.2021004
[10]
Boubacar Fall, Filippo Santambrogio, Diaraf Seck .
Shape derivative for obstacles in crowd motion. Mathematics in Engineering, 2022, 4(2): 1-16.
doi: 10.3934/mine.2022012
Abstract
Objective
Osteochondrodysplasias are a heterogeneous group of genetic skeletal dysplasias. Mutations in the COL2A1 gene cause a spectrum of rare autosomal-dominant type II collagenopathies characterized by skeletal dysplasia, short stature, and with vision and auditory defects. In this study, we have investigated in more detail the phenotypic and genotypic characterization resulting from glycine to serine mutations in the COL2A1 gene in a 2-year-old boy.
Materials and methods
Detailed clinical and radiological phenotypic characterization was the baseline tool to guide the geneticists toward proper genotypic confirmation.
Results
Genetic analysis revealed a de novo mutation, c.1681G>A (p.Gly561Ser), in the collagen type II alpha-1 gene (COL2A1). The identified variant showed impaired protein stability, and lead to dysfunction of type II collagen. In addition to pre and postnatal growth retardation, remarkable retardation of gross motor development and intellectual disability were noted. The latter was connected to cerebral malformations. The overall clinical phenotype of our current patient resembles spondyloepiphyseal dysplasia congenita (SEDC), but with extra phenotypic criteria.
Conclusions
The aim of this paper is twofold; firstly, raising awareness among orthopaedic surgeons when dealing with children manifesting multiple deformities, and secondly to broaden the clinical phenotype in patients with COL2A1 mutations of amino acid substitution (glycine to serine).
1.
Introduction
Let n and k be two positive integers. Denote by p(n,k) the number of partitions of the positive number n on exactly k parts. Then the partition class k is the sequence p(1,k),p(2,k),…,p(n,k),… We already know, see [1], all these values can be divided into the highest d0=LCM(1,2,…,k) sub sequences, each of which is calculated by the same polynomial.
Choose a sequence of k natural numbers such that: the first member is arbitrary, and the rest form an arithmetic progression with a difference d=m⋅d0,m∈N, starting from the chosen first member. For example:
x1=j,x2=j+d,…,xk=j+(k−1)⋅d,j∈N.
(1.1)
The corresponding number of partitions of the class k for the elements of the previous arithmetic progression's values is:
p(x1,k),p(x2,k),…,p(xk,k).
(1.2)
If the values, which are calculated using the same polynomial, multiplied by the corresponding binomial coefficients, form the alternate sum, we notice that the sum always has a value which is independent of x1, no matter how we form the sequence (1.1).
For the partition function of classes we already know the following results, see [1,2] for some details:
ⅰ) The values of the partition function of classes is calculated with one quasi polynomial.
ⅱ) For each class k the quasi polynomial consists of at most LCM(1,2,…,k) different polynomials, each of them consists of a strictly positive and an alternating part.
ⅲ) All polynomials within one quasi polynomial p(n,k) are of degree k−1.
ⅳ) All the coefficients with the highest degrees down to [k2] are equal for all polynomials (all of strictly positive) and all polynomials differ only in lower coefficients (alternating part).
ⅴ) The form of any polynomial p(n,k) is:
p(n,k)=a1nk−1+a2nk−2+⋯+ak,
(1.3)
where the coefficients a1,a2,…,ak are calculated in the general form.
Let us forget for a moment that the coefficients a1,a2,… are known in general form. Knowing that all values for partitions class of the sequence (1.1) are obtained by one polynomial p(n,k), it is possible to determine all unknown coefficients in a completely different way from that given in papers [1,2]. To determine k unknowns, a k equation is required. For this purpose, it is sufficient to know all the values of the sequence (1.2). To this end, we must form the system (1.4) and solve it. (For k=10, see [3]).
The system (1.4) can be solved by Cramer's Rule. For further analysis, we need to find the following determinants. We will start with the known Vandermonde determinant, see [4].
When we remove the first column and an arbitrary row from the previous determinant we obtain the Vandermonde determinant of one order less. The following results are known, see [4] and are needed for further exposure. If we remove the second column and an arbitrary a-th row from the determinant (1.5) we get
The label Δ(a,b)m means that from Δm remove the a-th row and b-th column from the set of variables xa.
2.
Invariants of the partitions classes
2.1. The first partition invariant of classes
Theorem 1. Let m,j and k be three positive integers and
I1(k,j,d)=k−1∑i=0(−1)i(k−1i)p(j+i⋅d,k),
where d=m⋅LCM(1,2,3,…,k). Then I1(k,j,d)=(−1)k−1dk−1k! and is independent of j. (I1(k,j,d) is the first partition invariant which exists in all classes.)
Proof. Among the values of the class k we choose the ones corresponding to the sequence (1.1), and they are given with the sequence (1.2). According to [2], all the elements in (1.2) can be calculated using the same polynomial p(n,k) with degree k−1. Elements of the following sequence:
q,q+d,…,q+(k−1)⋅d,q≠j,
are calculated with not necessarily the same polynomial as the previous one. Let the polynomial p(n,k) have the form as in (1.3). To determine the coefficients a1,a2,…,ak it suffices to know the k values: p(x1,k),p(x2,k),…,p(xk,k) where x1=j,x2=j+d,…,xk=j+(k−1)d are different numbers. Since Δk≠0, system (1.4) always has a unique solution, because all the elements of the set {x1,x2,…,xk} are different from one another. According to Cramer's Rule, to determine the coefficient of the highest degree of the polynomial (1.3), which calculates the value of the number of partitions of class k, we have the following formula:
The coefficient a1 is already defined in [2] where it is shown that a1=1k!(k−1)!. Substituting into the previous equality and multiplying by (−1)k−1, we obtain
Multiplying the last equality with (k−1)!dk−1 we obtain
(−1)k−1dk−1k!=k−1∑i=0(−1)i(k−1i)p(j+i⋅d,k),
which was to be proved. As these values are equal to each observed number of objects (1.2) within a class, the sum is invariant for any observed class.
All classes of the partition do not contain all the invariants we will list. This primarily refers to the classes from the beginning. Only the first invariant appears in all classes. The second invariant holds starting from the third class. The third invariant holds starting from the fifth class. Fourth, from the seventh class, etc. This coincides with the appearance of the common coefficients {ak} in quasi polynomials p(n,k), k∈N.
Theorem 2. Let m, j and k be three positive integers, k≥3 and
where d=m⋅LCM(2,3,…,k). Then I2(k,j,d)=(−1)k(k−3)dk−14(k−2)! and is independent of j.
Remark. In the previous expression, we should not simplify as then the value for k=3 cannot be obtained. However, the value for k=3 exists and is equal to zero.
Proof. Analogously to Theorem 1, the fact that the sum does not depend on the parameter j is a consequence of the periodicity per modulo LCM(2,3,…,k) using the same polynomial to calculate the partition class values.
In [2] it is shown how the system of linear equations can determine the other unknown coefficient of the polynomials which are calculated values of the partition classes. This coefficient is obtained from Cramer's Rule on system (1.4) and a2 is given by
In every subsequent invariant, the proceedings become more complex. But, it is quite clear how further invariants can be calculated.
3.
Consideration of special cases
For each partitions class k, k∈N we determine d0=LCM(1,2,3,…,k), and then form d=m⋅d0, m∈N. In addition arbitrarily choose the natural number j and than form sequences (1.1) and (1.2). Finally, we form an appropriate sum which is for the first invariant:
Sum (3.1) has a constant value in each partitions class and can be nominated as the first partitions class invariant.
3.1. The first partitions class invariant
For k=1, sum (3.1) has a constant value of 1.
For k=2, d0=2. If we choose some m∈N and set d=2m, the sum (3.1) has the form: p(j,2)−p(j+d,2),j∈N. According to [1], it is known that p(n,2)=[n2]. Distinguishing between even and odd numbers of j (j and j+d have the same parity) and substituting into the sum, we obtain that the result, in both cases, is equal to −d2=−m.
For k=3, d0=6. If we choose some m∈N and set d=6m the sum (3.1) has the form:
Similar to case k=3, by distinguishing the even and odd j and replacing (3.5) in relation (3.4) we obtain that the corresponding sums in both cases are equal to: −72m3. (Note that: i1=jmod12, i2=(j+d)mod12, i3=(j+2d)mod12, i4=(j+3d)mod12 and wi1=wi2=wi3=wi4.)
The number of invariants increases, when the class number increases. Starting with class three, another invariant can be observed.
3.2. The second partitions class invariant
Form in the same way as in the previous section: d0, d and the sequences (1.1) and (1.2) as well as the sum:
k−1∑i=0(−1)i(j(k−1)+((k2)−i)d)(k−1i)p(j+i⋅d,k).
Previous sum has a constant value in each partitions class (starting from third class) and can be nominated as the second partitions class invariant.
For k=3, d0=6. If we choose some m∈N and set d=6m the general form of the second invariant in the third class can be written as
The last equations can be verified in an analogous manner, by using the same form of the known polynomial for the fourth class given in (3.5). Note that: i1=jmod12, i2=(j+d)mod12, i3=(j+2d)mod12, i4=(j+3d)mod12 and wi1=wi2=wi3=wi4. By distinguishing the even and odd j and replacing (3.5) in relation (3.6) we obtain that the corresponding sums in both cases are equal to: −216m3.
3.3. The third partition invariants
Form in the same way as in the previous two section: d0, d and the sequences (1.1) and (1.2) as well as the sum I3(k,j,d) (Theorem 3). For each class (starting from the fifth) I3(k,j,d) has constant values and can be nominated as the third partitions class invariant. It is known [1] that
Using formulas from (3.7), we find that: p(1,5)=0, p(61,5)=5608, p(121,5)=80631, p(181,5)=393369 and p(241,5)=1220122, and so by checking we are assured of the accuracy.
Remark 2. Obviously, p(n,k) define values only for n≥k. The invariants determine very precisely that values for n<k should be taken as zero.
4.
Conclusions
In this paper, authors have demonstrated a new approach to partitions class invariants, as a way of proving the relevance and accuracy of all formulas given in [1,2]. Also, it I can be considered to be another way to obtain some of the formulas in [2]. The quasi polynomials p(n,k) needed to calculate the number of partitions of a number n to exactly k parts consists of at most LCM(1,2,…,k) different polynomials. The invariants claim that the more different polynomials in one quasi polynomial, the more invariable sizes connect them.
Acknowledgments
The author thank to The Academy of Applied Technical Studies Belgrade for partial funding of this paper.
Conflict of interest
Authors declare no conflicts of interest in this paper.
Acknowledgments
We wish to thank Ms. Katharina Sigl, head of the musculo-skeletal group (ordensklinikum) Linz, Austria for her help in facilitating the required investigations.
Author contributions
AAK, LBJ prepare clinical and radiographic documentation, and write the MS. HS, HJ, YE, SH, and MH contributed in patient's and parents investigations. MS and SGK performed screening for lysosomal storage disorders and whole exome sequencing. All authors approved the final version.
Conflict of interest
The authors declare that they have no conflict of interest.
References
[1]
Kannu P, Bateman J, Savarirayan R (2012) Clinical phenotypes associated with type II collagen mutations. J Paediatr Child Health 48: E38-43. doi: 10.1111/j.1440-1754.2010.01979.x
[2]
Spranger J, Winterpacht A, Zabel B (1994) The type II collagenopathies: A spectrum of chondrodysplasias. Eur J Pediatr 153: 56-65.
[3]
Cao LH, Wang L, Ji CY, et al. (2012) Novel and recurrent COL2A1 mutations in Chinese patients with spondyloepiphyseal dysplasia. Genet Mol Res 11: 4130-4137. doi: 10.4238/2012.September.27.1
[4]
Deng H, Huang X, Yuan L (2016) Molecular genetics of the COL2A1-related disorders. Mutat Res Rev Mutat Res 768: 1-13. doi: 10.1016/j.mrrev.2016.02.003
[5]
Nenna R, Turchetti A, Mastrogiorgio G, et al. (2019) COL2A1 gene mutations: Mechanisms of spondyloepiphyseal dysplasia congenita. Appl Clin Genet 12: 235-238. doi: 10.2147/TACG.S197205
[6]
Nishimura G, Haga N, Kitoh H, et al. (2005) The phenotypic spectrum of COL2A1 mutations. Hum Mutat 26: 36-43. doi: 10.1002/humu.20179
[7]
Barat-Houari M, Sarrabay G, Gatinois V, et al. (2016) Mutation Update for COL2A1 Gene Variants Associated with Type II Collagenopathies. Hum Mutat 37: 7-15. doi: 10.1002/humu.22915
[8]
Rukavina I, Mortier G, Van Laer L, et al. (2014) Mutation in the type II collagen gene (COL2AI) as a cause of primary osteoarthritis associated with mild spondyloepiphyseal involvement. Semin Arthritis Rheum 44: 101-104. doi: 10.1016/j.semarthrit.2014.03.003
[9]
Nazme NI, Sultana J, Chowdhury RB (2016) Escobar Syndrome - A Case Report in A Newborn. Bangladesh J Child Heal 39: 50-53. doi: 10.3329/bjch.v39i1.28359
Tafakhori A, Yu Jin Ng A, Tohari S, et al. (2016) Mutation in twinkle in a large Iranian family with progressive external ophthalmoplegia, myopathy, dysphagia and dysphonia, and behavior change. Arch Iran Med 19: 87-91.
[12]
Zhang X, Minikel EV, O'Donnell-Luriaa AH, et al. (2017) ClinVar data parsing. Wellcome Open Res 2: 33. doi: 10.12688/wellcomeopenres.11640.1
[13]
Fokkema IFAC, Taschner PEM, Schaafsma GCP, et al. (2011) LOVD v.2.0: The next generation in gene variant databases. Hum Mutat 32: 557-563. doi: 10.1002/humu.21438
[14]
Karczewski KJ, Francioli LC, Tiao G, et al. (2020) The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581: 434-443. doi: 10.1038/s41586-020-2308-7
[15]
Lek M, Karczewski KJ, Minikel EV, et al. (2016) Analysis of protein-coding genetic variation in 60,706 humans. Nature 536: 285-291. doi: 10.1038/nature19057
[16]
Karczewski KJ, Weisburd B, Thomas B, et al. (2017) The ExAC browser: Displaying reference data information from over 60 000 exomes. Nucleic Acids Res 45: D840-845. doi: 10.1093/nar/gkw971
[17]
Richards S, Aziz N, Bale S, et al. (2015) 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 17: 405-424. doi: 10.1038/gim.2015.30
[18]
Sim NL, Kumar P, Hu J, et al. (2012) SIFT web server: Predicting effects of amino acid substitutions on proteins. Nucleic Acids Res 40: W452-457. doi: 10.1093/nar/gks539
[19]
Adzhubei IA, Schmidt S, Peshkin L, et al. (2010) A method and server for predicting damaging missense mutations. Nat Methods 7: 248-249. doi: 10.1038/nmeth0410-248
[20]
Grantham R (1974) Amino acid difference formula to help explain protein evolution. Science 185: 862-864. doi: 10.1126/science.185.4154.862
[21]
Jagadeesh KA, Wenger AM, Berger MJ, et al. (2016) M-CAP eliminates a majority of variants of uncertain significance in clinical exomes at high sensitivity. Nat Genet 48: 1581-1586. doi: 10.1038/ng.3703
[22]
Pollard KS, Hubisz MJ, Rosenbloom KR, et al. (2010) Detection of nonneutral substitution rates on mammalian phylogenies. Genome Res 20: 110-121. doi: 10.1101/gr.097857.109
[23]
Schwarz JM, Cooper DN, Schuelke M, et al. (2014) Mutationtaster2: Mutation prediction for the deep-sequencing age. Nat Methods 11: 361-362. doi: 10.1038/nmeth.2890
[24]
Kõressaar T, Lepamets M, Kaplinski L, et al. (2018) Primer3-masker: Integrating masking of template sequence with primer design software. Bioinformatics 34: 1937-1938. doi: 10.1093/bioinformatics/bty036
[25]
Corpet F (1988) Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res 16: 10881-10890. doi: 10.1093/nar/16.22.10881
[26]
Liu YF, Chen WM, Lin YF, et al. (2005) Type II Collagen Gene Variants and Inherited Osteonecrosis of the Femoral Head. N Engl J Med 352: 2294-2301. doi: 10.1056/NEJMoa042480
[27]
Barat-Houari M, Dumont B, Fabre A, et al. (2016) The expanding spectrum of COL2A1 gene variants in 136 patients with a skeletal dysplasia phenotype. Eur J Hum Genet 24: 992-1000. doi: 10.1038/ejhg.2015.250
[28]
Terhal PA, Nievelstein RJAJ, Verver EJJ, et al. (2015) A study of the clinical and radiological features in a cohort of 93 patients with a COL2A1 mutation causing spondyloepiphyseal dysplasia congenita or a related phenotype. Am J Med Genet Part A 167: 461-475. doi: 10.1002/ajmg.a.36922
[29]
Husar-Memmer E, Ekici A, Al Kaissi A, et al. (2013) Premature Osteoarthritis as Presenting Sign of Type II Collagenopathy: A Case Report and Literature Review. Semin Arthritis Rheum 42: 355-360. doi: 10.1016/j.semarthrit.2012.05.002
[30]
Al Kaissi A, Laccone F, Karner C, et al. (2013) Hip dysplasia and spinal osteochondritis (Scheuermann's disease) in a girl with type II manifesting collagenopathy. Orthopade 42: 963-968. doi: 10.1007/s00132-013-2182-1
[31]
Chen J, Ma X, Zhou Y, et al. (2017) Recurrent c.G1636A (p.G546S) mutation of COL2A1 in a Chinese family with skeletal dysplasia and different metaphyseal changes: a case report. BMC Pediatr 17: 175. doi: 10.1186/s12887-017-0930-9
[32]
Anderson CE, Sillence DO, Lachman RS, et al. (1982) Spondylometepiphyseal dysplasia, Strudwick type. Am J Med Genet 13: 243-256. doi: 10.1002/ajmg.1320130304
This article has been cited by:
1.
Ching-Lung Lin, Hongyu Liu, Catharine W. K. Lo,
Uniqueness principle for fractional (non)-coercive anisotropic polyharmonic operators and applications to inverse problems,
2024,
0,
1930-8337,
0,
10.3934/ipi.2024054
Mohammad Shboul, Hela Sassi, Houweyda Jilani, Imen Rejeb, Yasmina Elaribi, Syrine Hizem, Lamia Ben Jemaa, Marwa Hilmi, Susanna Gerit Kircher, Ali Al Kaissi. The phenotypic spectrum in a patient with Glycine to Serine mutation in the COL2A1 gene: overview study[J]. AIMS Molecular Science, 2021, 8(1): 76-85. doi: 10.3934/molsci.2021006
Mohammad Shboul, Hela Sassi, Houweyda Jilani, Imen Rejeb, Yasmina Elaribi, Syrine Hizem, Lamia Ben Jemaa, Marwa Hilmi, Susanna Gerit Kircher, Ali Al Kaissi. The phenotypic spectrum in a patient with Glycine to Serine mutation in the COL2A1 gene: overview study[J]. AIMS Molecular Science, 2021, 8(1): 76-85. doi: 10.3934/molsci.2021006
Figure 1. Clinical phenotype of the proband showed the clinical features. (A) anterior view showed dysmorphic craniofacial features of wide and bulging frontal area (tower like head), hypotelorism and upslanting palpebral fissures and small penis associated with unilateral cryptorchidism. Synophrys, long eyelashes, anteverted nares, full cheeks and a barrel-shaped chest with thoracic kyphosis are also observed; (B) Lateral view showed normal limbs length, exaggerated lumbar lordosis, short trunk dwarfism and protuberant abdomen; (C) Lateral upper limb showed hypertrichosis of the arm were; (D) Photos of the hands revealed small and stubby fingers with clinodactyly of the fifth fingers
Figure 2. Radiographs of the proband. (A) Lateral spine radiograph showed platyspondyly of the lumbar vertebrae with spondylolisthesis of the 5th lumbar vertebra; (B) Lateral cervical and upper thoracic showed odontoid hypoplasia associated with massive compression of the cervical vertebrae; (C) AP pelvis radiograph showed marked epiphyseal dysplasia of the capital femoral epiphyses associated with acetabular dysplasia; (D) and (E) AP and lateral knee radiographs showed epiphyseal dysplasia and bifidity of the inferior femoral epiphysis
Figure 3. Electropherograms and partial COL2A1 alignment. (A) Sanger sequencing of the proband and his parents showed de novo heterozygous mutation c.1681G>A in COL2A1 gene. The red arrow indicates the position of the nucleotide mutation. (B) Cross-species conservation of the protein COL2A1 in Gly561 amino acid residue. The partial amino acid sequences were retrieved from the NCBI database, using the following accession numbers: Homo sapiens (Human; NP_001835.3), Pan troglodytes (Chimpanzee; XP_509026.2), Bos taurus (Cattle; NP_001001135.2), Mus musculus (Mouse; NP_112440.2), Rattus norvegicus (Rat; NP_037061.1), Gallus gallus (Chicken; NP_989757.1), zebrafish (Fish; NP_571367.1), Xenopus tropicalis (Frog; NP_989220.1). Alignment performed using MultAlin software [25]. The red arrow shows the conserved mutated amino acid residue, Gly561. The letter N represents the normal allele