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

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

FIGURE 1. X-ray manifestations of the recruited CIPA patient. (A,B) Front and side images of the torso. (C) Hands and wrists. (D–F) Ankles in different angles.

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

FIGURE 2. Genetic findings. (A) Pedigree diagram and the two variants detected in this case. Dark blue and red blocks represent the carrying status of these variants, respectively. (B) The location of these two variant illustrated in gene and protein schematics. (C) The conservatism of the amino acid Arg748 (R748) affected by c. 2242C > T variant across species.

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

FIGURE 3. Results of structure modeling and molecular dynamic simulation. (A) The structures of domain containing Arg748 (R748) or the Arg748Trp (R748W) mutant. Hydrogen bonds formed between the Arg748/Arg748Trp and Cys752 residues are shown in stick representation. The dotted yellow lines represent the hydrogen bonds involving Arg748. (B) The trajectory of RMSD (Cα) (root mean square deviation) of the two proteins. (C) RMSF (root mean square fluctuation) of the two proteins calculated from each simulation. (D) The number of hydrogen bonds formed between Arg748 (WT)/Arg748Trp (R748W)and the rest residues in the two protein models, respectively. (E) Secondary structural components of corresponding region in the two models (WT and R748W mutant) as a function of time.

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

FIGURE 4. Results of in vitro study. (A) The GFP fluorescent signal at 48 h after transfection. NC, no vector; WT, with wild-type MFN2 cDNA plasmid; MUT, with MFN2: c.638T > C mutant cDNA plasmid. Scale bar, 1000 μm; magnification, ×100. (B) The relative NTRK1mRNA levels in three cell groups. (C) The volcano plot showing significantly different compounds. Red dots represent the compounds with difference >2.0 folds and p < 0.05, blue font represents the m/z molecular weight information. (D) Part of the compounds with vip (variable importance in projection) value > 1. The X- and Y-axes stand for vip score and m/z, respectively. (E) KEGG pathway enrichment result. Each point represents a pathway, with X- and Y-axes indicating the importance of a compound related to the pathway and-log10(P) value, respectively. (F) The relative mRNA levels of various enzymes that can catalyze the hydrolysis of the AMP in WT and Mutant transfected cell groups.

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).

TABLE 1

TABLE 1. Metabolic pathways that significantly affected by the NTRK1^Arg748Trp variant (in reverse order of comprehensive significance).

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.

Data Availability Statement

The sequence data presented in the study are deposited in the repository at: https://doi.org/10.6084/m9.figshare.16934713.v1.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of Hebei Medical University Third Affiliated Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin. Written informed consent was obtained from the minor(s)’ legal guardian/next of kin for the publication of any potentially identifiable images or data included in this article.

Author Contributions

C-HY and KY designed the whole study scheme. Y-SY, YL, and Y-QS recruited the case and performed the clinical evaluation. Y-CX and H-YH conducted the genetic detection and data analysis. QL performed the structural and MD analysis. Y-ZL, Y-YL, and Z-KF conducted the in vitro study. KY and H-YH wrote the original version of this manuscript, and all co-authors have read it and proved the submission.

Funding

The work was supported by the Beijing Hospitals Authority Youth Programme (no.QML20211401), the National Key Research and Development Program of China (no. 2016YFC1000307), and the Hebei Natural Science Foundation Precision Medicine Joint Fund Cultivation Project (no. H2021206242).

Conflict of Interest

The 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.

Acknowledgments

The authors acknowledge the enrolled family for their participation in this study, and the expertise of Dr. Yuqing You from Beijing MyGenostics Technology Co., Ltd.

Supplementary Material

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

References

Bertrand, T., Kothe, M., Liu, J., Dupuy, A., Rak, A., Berne, P. 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. doi:10.1016/j.jmb.2012.08.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Bibel, M., and Barde, Y.-A. (2000). Neurotrophins: Key Regulators of Cell Fate and Cell Shape in the Vertebrate Nervous System. GENES DEVELOPMENT 14, 2919–2937. doi:10.1101/gad.841400

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y.-C., Auer-Grumbach, M., Matsukawa, S., Zitzelsberger, M., Themistocleous, A. C., Strom, T. 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. doi:10.1038/ng.3308

PubMed Abstract | CrossRef Full Text | Google Scholar

Darvish, H., Azcona, L. J., Tafakhori, A., Ahmadi, M., Ahmadifard, A., and Paisán-Ruiz, C. (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. doi:10.1038/s41525-017-0022-7

CrossRef Full Text | Google Scholar

Desiderio, S., Vermeiren, S., Van Campenhout, C., Kricha, S., Malki, E., Richts, S., 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. doi:10.1016/j.celrep.2019.02.097

PubMed Abstract | CrossRef Full Text | Google Scholar

Drissi, I., Woods, W. A., and Woods, C. G. (2020). Understanding the Genetic Basis of Congenital Insensitivity to Pain. Br. Med. Bull. 133, 65–78. PubMed PMID: 32219415. doi:10.1093/bmb/ldaa003

PubMed Abstract | CrossRef Full Text | Google Scholar

Duarte, R. R. R., Bachtel, N. D., Côtel, M.-C., Lee, S. H., Selvackadunco, S., Watson, I. 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. Psychiatry 86 (2), 120–130. doi:10.1016/j.biopsych.2019.03.977

PubMed Abstract | CrossRef Full Text | Google Scholar

Einarsdottir, E., Carlsson, A., Minde, J., Toolanen, G., Svensson, O., Solders, G., 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. doi:10.1093/hmg/ddh096

PubMed Abstract | CrossRef Full Text | Google Scholar

Esteban-Villarrubia, J., Soto-Castillo, J. J., Pozas, J., San Román-Gil, M., Orejana-Martín, I., Torres-Jiménez, J., et al. (2020). Tyrosine Kinase Receptors in Oncology. Ijms 21 (22), 8529. PubMed PMID: 33198314; PubMed Central PMCID: PMCPMC7696731. doi:10.3390/ijms21228529

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldberg, Y., MacFarlane, J., MacDonald, M., Thompson, J., Dube, M.-P., Mattice, M., 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. doi:10.1111/j.1399-0004.2007.00790.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, H. Y., Zhang, J., Qiu, W., Liang, C., Li, C. X., Wei, T. Y., et al. (2021). Comprehensive Strategy Improves the Genetic Diagnosis of Different Polycystic Kidney Diseases. J. Cel Mol Med 25 (13), 6318–6332. PubMed PMID: 34032358; PubMed Central PMCID: PMCPMC8256360. doi:10.1111/jcmm.16608

CrossRef Full Text | Google Scholar

Indo, Y. (2008). “NTRK1 Congenital Insensitivity to Pain with Anhidrosis,” in GeneReviews® [Internet]. M. P. Adam, H. H. Ardinger, R. A. Pagon, S. E. Wallace, L. J. H. Bean, G. Mirzaaet al. (Seattle (WA): University of Washington, Seattle)–2020.

Google Scholar

Indo, Y. (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. doi:10.1002/humu.1224

PubMed Abstract | CrossRef Full Text | Google Scholar

Indo, Y. (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. doi:10.1111/j.1399-0004.2012.01943.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Jia, C., Tang, Y., Kong, Y., Xia, Y., and Ma, L. (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. doi:10.3389/fped.2021.638190

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, N., Guo, S., Wang, Q., Duan, G., Sun, J., Liu, Y., et al. (2019). Heterogeneity of Clinical Features and Mutation Analysis of NTRK1 in Han Chinese Patients with Congenital Insensitivity to Pain with Anhidrosis. Jpr Vol. 12, 453–465. PubMed PMID: 30774415; PubMed Central PMCID: PMCPMC6348974. doi:10.2147/JPR.S188566

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Z., Liu, J., Liu, G., Cao, W., Liu, S., Chen, Y., 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. doi:10.1177/0300060517747164

CrossRef Full Text | Google Scholar

Naseer, M. I., Abdulkareem, A. A., Pushparaj, P. N., Bibi, F., and Chaudhary, A. 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. doi:10.3389/fgene.2020.00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Novarino, G., Fenstermaker, A. G., Zaki, M. S., Hofree, M., Silhavy, J. L., Heiberg, A. D., et al. (2014). Exome Sequencing Links Corticospinal Motor Neuron Disease to Common Neurodegenerative Disorders. Science 343 (6170), 506–511. PubMed PMID: 24482476; PubMed Central PMCID: PMCPMC4157572. doi:10.1126/science.1247363

PubMed Abstract | CrossRef Full Text | Google Scholar

Richards, S., Aziz, N., Aziz, N., Bale, S., Bick, D., Das, 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 (5), 405–423. PubMed PMID: 25741868; PubMed Central PMCID: PMCPMC4544753. doi:10.1038/gim.2015.30

PubMed Abstract | CrossRef Full Text | Google Scholar

Schon, K. R., Parker, A. P. J., and Woods, C. G. (2018). “Congenital Insensitivity to Pain Overview,” in GeneReviews® [Internet] Seattle (WA). M. P. Adam, H. H. Ardinger, R. A. Pagon, S. E. Wallace, L. J. H. Bean, G. Mirzaaet al. (Seattle: University of Washington), 1993–2020. [Updated 2020 Jun 11].

Google Scholar

Shatzky, S., Moses, S., Levy, J., Pinsk, V., Hershkovitz, E., Herzog, L., 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. doi:10.1002/1096-8628(20000619)92:5<353:aid-ajmg12>3.0.co;2-c

CrossRef Full Text | Google Scholar

Thomaz, A., Jaeger, M., Brunetto, A. L., Brunetto, A. T., Gregianin, L., de Farias, C. B., et al. (2020). Neurotrophin Signaling in Medulloblastoma. Cancers 12 (9), 2542. PubMed PMID: 32906676; PubMed Central PMCID: PMCPMC7564905. doi:10.3390/cancers12092542

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Li, M., and Hakonarson, H. (2010). ANNOVAR: Functional Annotation of Genetic Variants from High-Throughput Sequencing Data. Nucleic Acids Res. 38, e164. doi:10.1093/nar/gkq603

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, T., Li, H., Xiang, J., Wei, B., Zhang, Q., Zhu, Q., 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. doi:10.1177/0300060517691699

CrossRef Full Text | Google Scholar

Yang, K., Shen, M., Yan, Y., Tan, Y., Zhang, J., Wu, J., et al. (2019). Genetic Analysis in Fetal Skeletal Dysplasias by Trio Whole-Exome Sequencing. Biomed. Res. Int. 2019, 1–8. doi:10.1155/2019/2492590

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, F., Mao, B., Geng, X., Ren, X., Wang, Y., Guan, Y., 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. doi:10.1111/ene.14234

CrossRef Full Text | Google Scholar