Investigation of a Novel NTRK1 Variation Causing Congenital Insensitivity to Pain With Anhidrosis

Abstract

Background: Congenital insensitivity to pain with anhidrosis (CIPA), a rare autosomal recessive sensory neuropathy, was caused mainly by biallelic mutations in the NTRK1 gene. The pathogenesis of CIPA still needs further elucidation.

Methods: Here, we recruited a CIPA case and introduced whole-exome sequencing (WES) to identify the causative variation. Subsequently, an in silico molecular dynamic (MD) analysis was performed to explore the intramolecular impact of the novel missense variant. Meanwhile, in vitro functional study on the novel variant from a metabolomic perspective was conducted via the liquid chromatography–mass spectrometry (LC-MS) approach, of which the result was verified by quantitative real-time PCR (qRT-PCR).

Results: A novel compound heterozygous variation in NTRK1 gene was detected, consisting of the c.851–33T > A and c.2242C > T (p.Arg748Trp) variants. MD result suggested that p.Arg748Trp could affect the intramolecular structure stability. The results of the LC-MS and metabolic pathway clustering indicated that the NTRK1^Arg748Trp variant would significantly affect the purine metabolism in vitro. Further analysis showed that it induced the elevation of NT5C2 mRNA level.

Conclusion: The findings in this study extended the variation spectrum of NTRK1, provided evidence for counseling to the affected family, and offered potential clues and biomarkers to the pathogenesis of CIPA.

1 Introduction

Congenital insensitivity to pain with anhidrosis (CIPA, MIM #256800), also known as hereditary sensory and autonomic neuropathy IV (HSAN IV), is a rare autosomal recessive mendelian disorder mainly characterized by absence of nociception, anhidrosis, and increased S. aureus infection risk (Indo et al., 2008). Patients with CIPA do not feel pain from any noxious stimuli, including inflammation and heat (Goldberg et al., 2007), so they often develop finger-tip and oral cavity wounds due to self-mutilation in childhood and bony deformities due to recurrent fractures (Schon et al., 2018).

The prevalence of CIPA is currently unknown, whereas biallelic mutations in the NTRK1 gene account for the vast majority of all cases (3). Still, there were a few CIPA cases found to be caused by mutations in NGF and PRDM12 genes (Einarsdottir et al., 2004; Chen et al., 2015), suggesting the importance of differential diagnosis at molecular level. All these genes are involved in the differentiation and maturation of nociceptor, the specialized neuron for pain sensing (Drissi et al., 2020), suggesting their synergistic interaction in this process.

The NTRK1 (MIM #191315) gene, encoding the neurotrophic tyrosine kinase receptor type 1 (TrkA), belongs to a family of nerve growth factor receptors whose ligands include neurotrophins (Indo, 2012), which plays a pivotal role in regulating development of both the central and peripheral neurons (Bibel and Barde, 2000). NTRK1 is located in chromosome 1q23.1 region, spans over 25-kb (kilobase) region, and contains 17 exons (https://www.omim.org/). To date, more than 130 pathogenic variants for CIPA in NTRK1 have been indexed (HGMD professional v2019.4, The Human Gene Mutation Database, http://www.hgmd.cf.ac.uk/ac/index.php). Among different ethnic groups, there are different hotspot variants with high incidence (Shatzky et al., 2000; Indo, 2001; Zhao et al., 2020); yet, the effects of specific NTRK1 variants, especially missense variants, on molecular structure and protein function still deserve further investigation.

In this study, a family with a boy affected by CIPA was recruited and underwent a comprehensive clinical evaluation and genetic detection. Further investigation was conducted in two aspects. On one hand, molecular dynamic (MD) analysis was performed to predict the intramolecular impact of the novel variant and its potential biological effects; on the other hand, expression vector carrying the mutant NTRK1 cDNA was transfected into tool cell lines to discover its impact on cellular metabolic profile.

2 Materials and Methods

2.1 Subjects

A family with a 15-year-old boy exhibiting “intellectual disability, multiple fractures” was recruited in our center. We performed a comprehensive clinical evaluation on the proband.

2.2 Genetic Analysis

For all participants, genomic DNA was extracted from 200 μl of peripheral blood using the DNA Blood Mini Kit (Qiagen GmbH) according to the protocol of manufacturer.

Whole-exome sequencing (WES) was performed as described in previous study (Yang et al., 2019). Briefly, DNA fragments were hybridized and captured by xGenExome Research Panel of IDT (Integrated DNA Technologies, San Diego, United States). Libraries were tested for enrichment by quantitative fluorescence PCR, and the size distribution and concentration were determined using Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, United States). The NovaSeq 6000 platform (Illumina, San Diego, United States), along with 150-bp pair-end reads, was used for the sequencing of DNA with ∼300pM per sample using the NovaSeq Reagent Kit (Illumina). Sequencing raw reads (with the quality level Q30% > 90%) were aligned to the human reference genome (accession no. hg19/GRCh37, ftp://hgdownload.cse.ucsc.edu/goldenPath/hg19/chromosomes/) using the Burrows–Wheeler Aligner tool, and PCR duplicates were removed by using Picard v1.57 (http://picard.sourceforge.net/). Variant calling was performed with the Verita Trekker^® Variants Detection system (v 2.0, Berry Genomics, Inc., China) and Genome Analysis Toolkit (https://software.broadinstitute.org/gatk/) and then annotated and interpreted using in silico tools of ANNOVAR (v 2.0) and Enliven^® Variants Annotation Interpretation systems (Berry Genomics, Inc.) (Wang et al., 2010) on the basis of the common guideline by ACMG (American College of Medical Genetics and Genomics) (Richards et al., 2015).

Sanger sequencing was performed as validation method for suspected variations. The conservatism of amino acid (AA) affected by missense variant was analyzed using MEGA7 (http://www.megasoftware.net) with default parameters. In addition, the pathgenecity of missense variant was predicted with Revel (an ensemble method for predicting the pathogenicity of missense variants using the following tools individually: MutPred, FATHMM, VEST, PolyPhen, SIFT, PROVEAN, MutationAssessor, MutationTaster, LRT, GERP, SiPhy, phyloP, and phastCons; with the cutoff value of >0.7) as in our previous study (Hu et al., 2021).

2.3 Structural Modeling and Molecular Dynamic Simulation

SWISS-MODEL tool (https://swissmodel.expasy.org/) was used to model the wild-type (WT) and mutant (Mut) domains using 4F0I template covering the NTRK1:p.Arg748 site (http://www1.rcsb.org/).

We then ran the MD simulation for both models. The program CHARMM22 was used to add hydrogen atoms, N- and C-terminal patches to the models (Hu et al., 2021). The generated models were solvated and neutralized in a box with TIP3P water at a minimum of 13 Å between the model and the wall of the box. All simulations were run using NAMD 2.9 with periodic boundary conditions applied (Desiderio et al., 2019). The temperature was held at 300 K, whereas the pressure was controlled at 1 atm. The time step was set to 2 fs, the particle mesh Ewald method was applied to model the electrostatics, and the van der Waals interactions cutoff was set at 12 Å. Both models followed a three-step pre-equilibration totaling 600 ps, the last snapshots of which were chosen as the starting structures for 20-ns productive simulations without constraints.

2.4 Plasmids Construction, Cell Culture, and Transfection

To explore the functional impact on cellular metabolism by the novel missense variant, expression plasmid vectors were constructed using a lentiviral backbone (pLV-hef1a-mNeongreen-P2A-Puro-WPRE-CMV-MCS-3×flag) containing the WT or Mut (with variant c.2242C > T) NTRK1 cDNA (see detailed data in Supplementary Material S1).

HEK-293 (human embryonic kidney) cell lines were purchased commercially (SyngenTech Inc., Beijing, China) and cultured by regular means with DMEM medium +10% BSA (Invitrogen, United States). The plasmid vectors were transfected using Lipofectamine 3000 (Thermo Fisher Scientific, United States). Forty-eight hours later, the intensity of green fluorescent protein (GFP) was observed by fluorescence microscope and the expression of NTRK1 was detected by quantitative real-time PCR (qRT-PCR) with ABI 7500 system (Thermo Fisher Scientific, United States) to determine the transfection effect (detailed methods in Supplementary Material S1).

2.5 Metabonomic Analysis

The cells were submitted to metabonomic analysis 48 h after transfection. The ultra-performance liquid chromatography combined with quadrupole time-of-flight tandem mass spectrometry analysis was performed on Nexera X2 system (Shimadzu, Japan) coupled with a TripleTOF 5600 quadrupole-time-of-flight mass spectrometer (AB SCIEX, United States). Liquid chromatography separation was performed on a ZORBAX Eclipse Plus C18 column (2.1 × 100 mm, 3.5 μm, Agilent, United States) maintained at 45°C. Independent reference lock mass ions via Analyst TF 1.6 and MarkerView 1.2.1 were used to ensure mass accuracy during data acquisition. The assigned modified metabolite ions were identified by database searches in the HMDB (http://www.hmdb.ca/spectra/ms/search) databases. The mass tolerance for the HMDB database search was set at 0.05 Da.

2.5.1 Statistical Analysis

Mann–Whitney U test were first performed to compare the WT (no. pHS-AVC-LY059) group with Mut (no. pHS-AVC-LY 060) groups. Orthogonal projection to latent structure–discriminant analysis (PLS-DA) was used to determine the distributions and find the metabolic differences between the two groups using MetaboAnalyst 5.0 (http://www.metaboanalyst.ca/MetaboAnalyst/). The PLS-DA models were cross-validated using a 10-fold method with unit variance scaling. The parameter R2 was used to evaluate the fitting condition for the PLS-DA models, and Q2 was used to assess predictive ability. Negative or very low Q2 values indicated that the differences between groups were not statistically significant. The PLS-DA model removes variation in the X matrix that is not correlated with the Y matrix. Thus, normally, only one predictive component is used for the discrimination between two classes.

2.5.2 Pathway Analysis

Chemical metabolites were analyzed using the MetaboAnalyst web portal for pathway analysis and visualization (http://www.metaboanalyst.ca/). Pearson correlations were used to evaluate relationships between metabolites (p < 0.05 and impact > 0.01 were considered significant).

2.5.3 qRT-PCR

It was used to detect the mRNA relative expression level of those metabolic enzymes that were suspectedly affected (detailed data shown in Supplementary Material S1).

3 Results

3.1 Clinical Manifestations

According to medical history survey, the patient has developed bilateral ankle swelling and fractures since 18 months. His intellectual development was significantly behind peers in motor and language development. Multiple fractures to his right ankle and swelling of his ankle and knee limited his movement. X-ray result showed scoliosis, osteomyelitis-like changes in the right ankle, bone deformation, bone destruction on multiple joints, and osteomyelitis-like changes after partial vertebral deformation (Figure 1).

FIGURE 1

3.2 Variation Analysis of the NTRK1 Gene

WES identified a compound heterozygous variation in the NTRK1 gene, consisting of two variant, namely, (NM_002529.3) c.851–33T > A and c.2242C > T (p. Arg748Trp) (Figure 2A). Figure 2B showed the location of the two variants in the NTRK1 gene and peptide chain schematics. Among these two, the former is prevalent in East-Asian populations, whereas the latter is an unreported novel missense variant. The revel score was 0.650 indicating that this variant was uncertain but likely to be pathogenic (Supplementary Table S1 in Supplementary Material S1).

FIGURE 2

The Arg748 residue maintained evolutionary conserved among species (Figure 2C). To elucidate the effect of p. Arg748Trp on molecular structure and protein function, we conducted MD simulation and in vitro experiments. The results are as follows.

3.3 Intramolecular Impact of the NTRK1: p.Arg748Trp Variant

The WT and Mut models were shown in Figure 3A. It was indicated that the p.Arg748Trp variant replacing the strongly basic arginine by a large AA with benzene rings broke the hydrogen bonds formed by the side chain of Arg748 and expectantly changed its potential distribution.

FIGURE 3

As shown by MD, it indicated that the Arg748Trp Mut was more flexible than the WT according to the trajectory of root mean square deviation and root mean square fluctuation (Figures 3B,C). Besides, the hydrogen bonds formed between Trp748, and the rest residues were more than those between Arg748 and the rest (Figure 3D), which might cause the corresponding loop to be less flexible. Moreover, Arg748Trp could influence the secondary structure of the connection region between two β sheets (Figure 3E).

3.4 Metabolic Impact of the NTRK1^Arg748Trp Mutant

According to the intracellular fluorescence signal of GFP (Figure 4A) and the relative expression level of NTRK1 mRNA (Figure 4B), the transfection efficiency was ideal.

FIGURE 4

Metabolites with significant difference between the two groups (WT, Lab No. pHS-AVC-LY061; Mut, pHS-AVC-LY062) after transfection were demonstrated in Figures 4C,D, and also listed in Supplementary Table S1 (the screening criteria for compounds with difference were FC > 2.0 or < 0.5, p value <0.05, vip >1. FC, fold count; vip, variable importance in projection). Pathway clustering result indicated that the metabolic pathway most significantly affected by NTRK1^Arg748Trp Mut was the purine metabolism (Figure 4E; Table 1). To be specific, the volume of AMP declined, whereas that of adenosine increased (Supplementary Table S1). Subsequently, among all enzymes catalyzing the hydrolysis of the AMP, qRT-PCR revealed a significant up-regulation of NT5C2 (5-Ppime-Nucleotidase, Cytosolic II, MIM #600417) mRNA expression in the Mut-transfected cells compared to WT cells (Figure 4F).

4 Discussion

Pain is essential in teaching us how to use our bodies optimally and avoid or respond to injuries while being permanently painless results in a significant morbidity and mortality (Drissi et al., 2020). The normal development and function of nociceptors is essential for pain sensation, and at least, nine genes have been recognized to play crucial roles in these progresses (Drissi et al., 2020). Among them, NGF, TrkA (NTRK1), and Prdm12 (PRDM12) collaboratively participate in the promotion of nociceptor fate (Desiderio et al., 2019); correspondingly, pathogenic variations in these genes result in autosomal recessive conditions with similar symptoms, namely, HSAN V (MIM #608654), HSAN IV (CIPA), and HSAN VIII (MIM #616488).

In this study, the patient displayed typical CIPA manifestations such as multiple fractures, ankle swelling, intellectual delay, rachioscoliosis, and osteomyelitis. The subsequent WES detection identified the causative compound heterozygous NTRK1variation consisting of c.851–33T > A and c.2242C > T (p. Arg748Trp) (NM_001012331). As we stated previously, the former variant is a hotspot mutation in East Asian people (Indo, 2001; Wang et al., 2017; Liu et al., 2018; Li et al., 2019; Zhao et al., 2020; Li et al., 2021); whereas the latter one is a novel missense variant. According to the common interpretation guidelines provided by ACMG (14), it was recognized as “uncertain significance” (with evidences of PM1+PM2+PP2+BP4). The BP4 evidence was due to the “uncertain” result by Revel and benign computational verdict by several tools in the Varsome predicting website (https://varsome.com/). Therefore, further investigation was necessary to confirm the pathogenicity of this variant.

The conservatism of Arg748 AA residue across species indirectly suggested the pathogenicity of p. Arg748Trp variant. Besides, according to the structural modeling and MD simulation result, the p.Arg748Trp variant might change the conformation of the loop which the Arg748 belongs to (Bertrand et al., 2012). Thus, the variant could probably cause the structural instability of TrkA protein and further broke its binding to NGF. Yet, more mechanistic experiments are needed to validate whether this variant could cause loss of protein function.

With regard to the in vitro study, we intended to analyze the cellular effect of NTRK1^Arg748Trp variant from the perspective of metabolic changes. It was demonstrated that the purine metabolism pathway was significantly affected. Furthermore, it seemed that the up-regulation of NT5C2 mRNA level should be a major contribution to this effect. Duarte et al. demonstrated that NT5C2 could participate in the neural development and motility of Drosophila by regulating the AMPK signaling pathway (Duarte et al., 2019). Singgih et al. showed that impairment in cNT5-II (NT5C2 orthologs in Drosophila) resulted in disease-relevant behaviors in Drosophila. Novarino et al. revealed that NT5C2 mutations could lead to a severe neurodegenerative motor neuron disease, hereditary spastic paraplegia (Spastic paraplegia 45, MIM #613162) (Novarino et al., 2014), which was validated by other studies (Darvish et al., 2017; Naseer et al., 2020). Collectively, it can be deduced that NT5C2 plays a crucial role in the development and functional maintenance of neurosystem. Up to our knowledge, the NGF/TrkA signaling could participate in the process of cell survival, proliferation, and migration through at least three pathways, namely, RAS/RAF/MEK/ERK, PI3K/Akt, and PLC/PCK (Thomaz et al., 2020; Esteban-Villarrubia et al., 2020). Through in silico analysis, we found that there were multiple binding sites of transcriptional factors such as STAT1β, ELK1, and NF-κB at the upstream promoter region of NT5C2 gene (Supplementary Material 1), and all these factors are associated with at least one of the above pathways (27). That is to say, NGF/TrkA could probably regulate NT5C2 transcription via one of these signaling pathways and jointly involved in the development of neurons; yet, detailed functional experiments are also needed to verify this hypothesis.

The future pregnancy of this couple will have a 25% chance of getting affected, so potential interventions such as prenatal diagnosis or pre-implantation diagnosis need to be recommended. The major limitation of our study was that the evidences were still not solid to elucidate the functional impact of the novel NTRK1 variant. As a necessary supplement, the alteration of interaction between NGF and TrkA by it and how TrkA regulate NT5C2 should be clarified in our further investigation.

In summary, we present with a novel compound heterozygous NTRK1 variation including one novel variant, c.2242C > T (p. Arg748Trp). The findings in this study extend the NTRK1 mutation spectrum and provide critical information for genetic counseling to guide the future pregnancy of the family. In addition, a potential mechanism of how TrkA relates NT5C2 was also suggested, which might be essential in the development of neurons.

References

• 1

  BertrandT.KotheM.LiuJ.DupuyA.RakA.BerneP. F.et al (2012). The crystal Structures of TrkA and TrkB Suggest Key Regions for Achieving Selective Inhibition. J. Mol. Biol.423 (3), 439–453. PubMed PMID: 22902478. 10.1016/j.jmb.2012.08.002

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BibelM.BardeY.-A. (2000). Neurotrophins: Key Regulators of Cell Fate and Cell Shape in the Vertebrate Nervous System. GENES DEVELOPMENT14, 2919–2937. 10.1101/gad.841400

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  ChenY.-C.Auer-GrumbachM.MatsukawaS.ZitzelsbergerM.ThemistocleousA. C.StromT. M.et al (2015). Transcriptional Regulator PRDM12 Is Essential for Human Pain Perception. Nat. Genet.47 (7), 803–808. PubMed PMID: 26005867; PubMed Central PMCID: PMCPMC7212047. 10.1038/ng.3308

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  DarvishH.AzconaL. J.TafakhoriA.AhmadiM.AhmadifardA.Paisán-RuizC. (2017). Whole Genome Sequencing Identifies a Novel Homozygous Exon Deletion in the NT5C2 Gene in a Family with Intellectual Disability and Spastic Paraplegia. Npj Genomic Med.2. PubMed PMID: 29123918; PubMed Central PMCID: PMCPMC5675118. 10.1038/s41525-017-0022-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  DesiderioS.VermeirenS.Van CampenhoutC.KrichaS.MalkiE.RichtsS.et al (2019). Prdm12 Directs Nociceptive Sensory Neuron Development by Regulating the Expression of the NGF Receptor TrkA. Cel Rep.26 (13), 3522–3536e5. PubMed PMID: 30917309. 10.1016/j.celrep.2019.02.097

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  DrissiI.WoodsW. A.WoodsC. G. (2020). Understanding the Genetic Basis of Congenital Insensitivity to Pain. Br. Med. Bull.133, 65–78. PubMed PMID: 32219415. 10.1093/bmb/ldaa003

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  DuarteR. R. R.BachtelN. D.CôtelM.-C.LeeS. H.SelvackaduncoS.WatsonI. A.et al (2019). The Psychiatric Risk Gene NT5C2 Regulates Adenosine Monophosphate-Activated Protein Kinase Signaling and Protein Translation in Human Neural Progenitor Cells. Biol. Psychiatry86 (2), 120–130. 10.1016/j.biopsych.2019.03.977

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  EinarsdottirE.CarlssonA.MindeJ.ToolanenG.SvenssonO.SoldersG.et al (2004). A Mutation in the Nerve Growth Factor Beta Gene (NGFB) Causes Loss of Pain Perception. Hum. Mol. Genet.13 (8), 799–805. PubMed PMID: 14976160. 10.1093/hmg/ddh096

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  Esteban-VillarrubiaJ.Soto-CastilloJ. J.PozasJ.San Román-GilM.Orejana-MartínI.Torres-JiménezJ.et al (2020). Tyrosine Kinase Receptors in Oncology. Ijms21 (22), 8529. PubMed PMID: 33198314; PubMed Central PMCID: PMCPMC7696731. 10.3390/ijms21228529

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  GoldbergY.MacFarlaneJ.MacDonaldM.ThompsonJ.DubeM.-P.MatticeM.et al (2007). Loss-of-function Mutations in the Nav1.7 Gene Underlie Congenital Indifference to Pain in Multiple Human Populations. Clin. Genet.71 (4), 311–319. PubMed PMID: 17470132. 10.1111/j.1399-0004.2007.00790.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 11

  HuH. Y.ZhangJ.QiuW.LiangC.LiC. X.WeiT. Y.et al (2021). Comprehensive Strategy Improves the Genetic Diagnosis of Different Polycystic Kidney Diseases. J. Cel Mol Med25 (13), 6318–6332. PubMed PMID: 34032358; PubMed Central PMCID: PMCPMC8256360. 10.1111/jcmm.16608

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  IndoY. (2008). “NTRK1 Congenital Insensitivity to Pain with Anhidrosis,” in GeneReviews® [Internet]. AdamM. P.ArdingerH. H.PagonR. A.WallaceS. E.BeanL. J. H.MirzaaG.et al (Seattle (WA): University of Washington, Seattle)–2020.

  • Google Scholar
  • View reference in article

• 13

  IndoY. (2001). Molecular Basis of Congenital Insensitivity to Pain with Anhidrosis (CIPA): Mutations and Polymorphisms inTRKA (NTRK1) Gene Encoding the Receptor Tyrosine Kinase for Nerve Growth Factor. Hum. Mutat.18, 462–471. 10.1002/humu.1224

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  IndoY. (2012). Nerve Growth Factor and the Physiology of Pain: Lessons from Congenital Insensitivity to Pain with Anhidrosis. Clin. Genet.82 (4), 341–350. PubMed PMID: 22882139. 10.1111/j.1399-0004.2012.01943.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  LiL.JiaC.TangY.KongY.XiaY.MaL. (2021). Novel Gross Deletion Mutations in NTRK1 Gene Associated with Congenital Insensitivity to Pain with Anhidrosis. Front. Pediatr.9, 638190. PubMed PMID: 33748046; PubMed Central PMCID: PMCPMC7969531. 10.3389/fped.2021.638190

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  LiN.GuoS.WangQ.DuanG.SunJ.LiuY.et al (2019). Heterogeneity of Clinical Features and Mutation Analysis of NTRK1 in Han Chinese Patients with Congenital Insensitivity to Pain with Anhidrosis. JprVol. 12, 453–465. PubMed PMID: 30774415; PubMed Central PMCID: PMCPMC6348974. 10.2147/JPR.S188566

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  LiuZ.LiuJ.LiuG.CaoW.LiuS.ChenY.et al (2018). Phenotypic Heterogeneity of Intellectual Disability in Patients with Congenital Insensitivity to Pain with Anhidrosis: A Case Report and Literature Review. J. Int. Med. Res.46 (6), 2445–2457. PubMed PMID: 29619836; PubMed Central PMCID: PMCPMC6023048. 10.1177/0300060517747164

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  NaseerM. I.AbdulkareemA. A.PushparajP. N.BibiF.ChaudharyA. G. (2020). Exome Analysis Identified Novel Homozygous Splice Site Donor Alteration in NT5C2 Gene in a Saudi Family Associated with Spastic Diplegia Cerebral Palsy, Developmental Delay, and Intellectual Disability. Front. Genet.11, 14. PubMed PMID: 32153630; PubMed Central PMCID: PMCPMC7050623. 10.3389/fgene.2020.00014

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  NovarinoG.FenstermakerA. G.ZakiM. S.HofreeM.SilhavyJ. L.HeibergA. D.et al (2014). Exome Sequencing Links Corticospinal Motor Neuron Disease to Common Neurodegenerative Disorders. Science343 (6170), 506–511. PubMed PMID: 24482476; PubMed Central PMCID: PMCPMC4157572. 10.1126/science.1247363

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  RichardsS.AzizN.AzizN.BaleS.BickD.DasS.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 (5), 405–423. PubMed PMID: 25741868; PubMed Central PMCID: PMCPMC4544753. 10.1038/gim.2015.30

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  SchonK. R.ParkerA. P. J.WoodsC. G. (2018). “Congenital Insensitivity to Pain Overview,” in GeneReviews® [Internet] Seattle (WA). AdamM. P.ArdingerH. H.PagonR. A.WallaceS. E.BeanL. J. H.MirzaaG.et al (Seattle: University of Washington), 1993–2020. [Updated 2020 Jun 11].

  • Google Scholar
  • View reference in article

• 22

  ShatzkyS.MosesS.LevyJ.PinskV.HershkovitzE.HerzogL.et al (2000). Congenital Insensitivity to Pain with Anhidrosis (CIPA) in Israeli-Bedouins: Genetic Heterogeneity, Novel Mutations in theTRKA/NGF Receptor Gene, Clinical Findings, and Results of Nerve Conduction Studies. Am. J. Med. Genet.92, 353–360. 10.1002/1096-8628(20000619)92:5<353:aid-ajmg12>3.0.co;2-c

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  ThomazA.JaegerM.BrunettoA. L.BrunettoA. T.GregianinL.de FariasC. B.et al (2020). Neurotrophin Signaling in Medulloblastoma. Cancers12 (9), 2542. PubMed PMID: 32906676; PubMed Central PMCID: PMCPMC7564905. 10.3390/cancers12092542

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  WangK.LiM.HakonarsonH. (2010). ANNOVAR: Functional Annotation of Genetic Variants from High-Throughput Sequencing Data. Nucleic Acids Res.38, e164. 10.1093/nar/gkq603

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  WangT.LiH.XiangJ.WeiB.ZhangQ.ZhuQ.et al (2017). Identification of a Novel Nonsense Mutation of the Neurotrophic Tyrosine Kinase Receptor Type 1 Gene in Two Siblings with Congenital Insensitivity to Pain with Anhidrosis. J. Int. Med. Res.45 (2), 549–555. PubMed PMID: 28345382; PubMed Central PMCID: PMCPMC5536645. 10.1177/0300060517691699

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  YangK.ShenM.YanY.TanY.ZhangJ.WuJ.et al (2019). Genetic Analysis in Fetal Skeletal Dysplasias by Trio Whole-Exome Sequencing. Biomed. Res. Int.2019, 1–8. 10.1155/2019/2492590

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  ZhaoF.MaoB.GengX.RenX.WangY.GuanY.et al (2020). Molecular Genetic Analysis in 21 Chinese Families with Congenital Insensitivity to Pain with or without Anhidrosis. Eur. J. Neurol.27, 1697–1705. PubMed PMID: 32219930. 10.1111/ene.14234

  • CrossRef
  • Google Scholar
  • View reference in article