Assessment of rare genetic variants to identify candidate modifier genes underlying neurological manifestations in neurofibromatosis 1 patients

doi:10.21203/rs.3.rs-2012218/v1

Download PDF

Research Article

Assessment of rare genetic variants to identify candidate modifier genes underlying neurological manifestations in neurofibromatosis 1 patients

https://doi.org/10.21203/rs.3.rs-2012218/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Neurological phenotypes such as intellectual disability occur in almost half of patients with neurofibromatosis 1 (NF1). Current genotype-phenotype studies have failed to reveal the mechanism underlying this clinical variability. Despite the presence of pathogenic variants of NF1, modifier genes likely determine the occurrence and severity of neurological phenotypes. Exome sequencing data were used to identify genetic variants in 13 NF1 patients and 457 healthy controls, and this information was used to identify candidate modifier genes underlying neurological phenotypes based on the Optimal Sequence Kernel Association Test. Thirty-six genes were identified as significant modifying factors in patients with neurological phenotypes, 35 of which were highly expressed in the nervous system. A review of the literature confirmed that 19 genes including CUL7, DPH1 and PRND were clearly associated with the alteration of neurological functioning and development. Our study revealed the enrichment of rare variants of 19 genes closely related to neurological development and functioning in patients with NF1 with neurological phenotypes, indicating possible modifier genes and variants affecting neurodevelopment. Further studies on rare genetic variants of candidate modifier genes may help explain the clinical heterogeneity of NF1.

NF1

rare variants

modifier gene

nervous system

genetic diagnosis

Neurofibromatosis type 1 (NF1; MIM#162200) is one of the most common inherited neurogenic disorders, with a worldwide birth incidence of 1 in 2,500 [1]. Dominant loss-of-function mutations involving NF1, located on chromosome 17q11.2, are responsible for NF1. Clinical characteristics of NF1 develop with age and include café-au-lait spots, axillary and inguinal freckling, cutaneous and subcutaneous neurofibromas, and Lisch nodules. Other less common features include plexiform neurofibromas, learning disabilities, and musculoskeletal features such as osteopenia, dysplasia of the long bones, and scoliosis [2]. Pharmacological therapy for plexiform neurofibromas can control the development of malignant peripheral nerve sheath tumors, which are thought to be the most serious disease feature. Neurological phenotypes tend to be of greatest consequence to patients, since there is no treatment that effectively improves quality of life. Although several studies have revealed genotype-phenotype correlations for specific NF1 mutations, the clinical expression of NF1 is highly variable and unpredictable, even among patients with the same mutation.

Currently, no evidence directly shows whether neurological abnormalities in children with NF1 including learning impairment, intellectual impairment, and autism are related to a specific mutation. Further, no correlation between the disease onset and neurological symptom severity has been shown, even in patients harboring the same genetic variants. Variations in the mutant allele of NF1 itself cannot fully explain disease variability, indicating that modifier genes, environmental factors, or a combination of both factors likely underly the phenotypic variability of NF1.

Neurological symptoms are likely to be affected by environmental factors. However, the incidence of neurological symptoms in children with NF1 is much higher than that which would be predicted due the impact of environmental factors. Evidence from twin studies suggests that variability may be caused by unidentified modifier genes [3, 4]. Genetic modifiers, defined as genetic variants that can modify the phenotypic outcome of a primary disease-causing variant, have been shown to affect phenotypes of some Mendelian disorders [5, 6]. Thus, the existence of rare variants of potential modifier genes and their impact on the phenotypic variability of neurological development in NF1 patients warrants further investigation.

In this study, we evaluated collective effects of rare variants on the modification of neurodevelopmental phenotypes of NF1. For this purpose, information regarding the presence of variants among 13 patients with NF1 who had various degrees of neurological involvement was analyzed. Subsequently, we performed an Optimal Sequence Kernel Association Test (SKAT-O) to identify genetic modifiers of neurological phenotypes. Our preliminary study provides new insights into possible modifier genes and variants in rare diseases. The list of novel candidate modifier genes and disease pathways may help researchers better understand the cause of phenotypic variability in patients with NF1.

2.1 Patient enrollment

Patients previously diagnosed with NF1 at the Department of Medical Genetics of Shanghai Children’s Medical Center, Shanghai Jiaotong University School of Medicine (Shanghai, China) were included in this study. Ethical approval for this study was obtained from the ethics committee of Shanghai Children’s Medical Center, Shanghai Jiaotong University School of Medicine (project identification code: SCMCIRB-Y2019021; date of approval: 2/25/2019). Written informed consent was obtained from each proband’s parents. Patients' medical charts were reviewed to identify information related to neurological involvement including learning disabilities, behavioral problems, autism spectrum disorder, attention deficit hyperactivity disorder, seizures, and sleep disturbance. Genomic data of healthy adults without any neurological involvement evaluated by physicians were selected as control.

2.2 Sample preparation

Peripheral blood samples were collected from patients and their parents. The SureSelect Human All Exon V6 enrichment kit (Agilent, Santa Clara, CA, USA) was used to prepare a sequencing library. The experimental process included the enzymatic digestion of DNA fragments, library hybridization, and library amplification and purification. The Hiseq2500 sequencing platform (Illumina, Inc., San Diego, CA, USA) was used for high-throughput sequencing. FastQC version 0.11.9 (Babraham Research Institute, Cambridge, UK) and Fastp version.0.20.1 (Visible Genetics, Inc., Toronto, Canada) were used for data quality control and to remove the adaptor sequence. Speedseq version 0.1.2 (Ira Hall Lab, St. Louis, MO, USA) was used to perform sequence alignment to reference genome. Further, bamdst (version 1.0.9) and mosdepth (version 0.3.1) were used to count sequencing indexes of BAM files after alignment including mapping rate, polymerase chain reaction duplication rate, average sequencing depth, and coverage rate. Genome Analysis Toolkit version 4.2.0.0 (Broad Institute, Cambridge, MA, USA) was used to detect variants in the BAM file passing a quality control test that was performed in accordance with the best practice guideline [7]. We began by calling variants per sample to produce a file in GVCF format. Next, we consolidated GVCFs from multiple samples within a GenomicsDB database. We then performed joint genotyping. Finally, VQSR filtering was applied to produce the final multi-sample call set with the desired balance of precision and sensitivity. Copy number variants (CNVs) were identified using the open-source CNVkit software, which is a tool kit that can infer and visualize copy numbers from targeted DNA sequencing data [8].

All VCF files of NF1 cases and controls were merged using bcftools (v.1.12). Variants were annotated using ANNOVAR (version Jun 8, 2020) with sequence ontologies, MAFs(Minor Allele Frequency) from the Genome Aggregation Database (gnomAD, v2.1.1), and in silico prediction scores from Combined Annotation Dependent Depletion [CADD], v1.3). Rare variants with a MAF < 0.01 in gnomAD were included in subsequent analyses. Based on annotations, we collected two groups of variants according to their variation type and possibility of pathogenicity, as follows: (1) putative loss-of-function (LOF) variants (including nonsynonymous, stop-gain, start-loss, and splicing) and (2) Likely pathogenic (LP) variants (including missense variants with a CADD phred ≥ 20) (Fig. 1).

2.3 Statistics

The SKAT-O was developed by estimating the correlation parameter in the kernel matrix, producing a value equal to the estimated weight of the linear combination of the burden and SKAT test statistics that maximizes power [9]. Thus, we were able to flexibly apply the test in the following two common analysis scenarios: (1) a region with many noncausal variants or if causal variants have different directions of association and (2) a region with a high proportion of causal variants with an association in the same direction.

The SKAT-O test was implemented using the R package “SKAT” method of Seunggeun Lee. The method aggregates individual SNP score statistics within an SNP set and efficiently computes SNP set p-values. SKAT-O statistics were calculated using putative LOF and LP variants. Only genes with p-values < 0.005 were reported (Fig. 1).

Follow-up analysis was performed on SKAT-O results with a detected signal using Firth logistic regression. The brglm2 R package was used to fit the regression model and apply a mean bias reduction accounting for the low variant-positive counts. Genes with total rare variant counts < 5 between the case and control cohorts were excluded from the analysis.

2.4 Pathway analysis

GO, KEGG, and Reactome pathway associations between genes flagged as significant via SKAT-O statistics were performed using the clusterProfiler R package [10–12]. For the gene-set enrichment analysis, all genes listed in the database were used as background. The main goal of the pathway analysis was to group genes of interest into pathways and explore the possibility that these pathways may be involved in the modulation of neurological phenotypes of NF1. An additional search of the PubMed database was performed to determine functional data of identified genes. Gene interaction networks were plotted using STRING and GeneMANIA [13, 14].

3.1 Patient Characteristics

Thirteen patients diagnosed with NF1 using NIH diagnostic criteria were included in the analysis. The average age of the patients was 7 years and 2 months, and the average age at initial diagnosis was 4 years and 9 months. Pathogenic variants in or encompassing NF1 were detected through whole exome sequencing and data analysis. Further validation within a pedigree confirmed that 10 patients carried de novo variants of NF1. Two patients (19654, 20586) inherited pathogenic NF1 variants from their parents who were also diagnosed with NF1. All patients were evaluated within the Department of Developmental Behavioral Pediatrics for neurological manifestations (Table 1). A total of 457 healthy adults (average age 32 years and 2 months) without any neuropsychiatric symptoms were selected as controls in this study.

Table 1

NF1 variants and related neurological involvement in patients
Patient ID	Gender	Age	Age at diagnosis	NF1 variants	inheritance	NF1 related neurological involvement
14646	male	10y	8y6m	c.2125T > C p.Cys709Arg	dn	ADHD,sleep disturbance
14811	male	12y6m	7y9m	c.3739_3742del p.Phe1247IlefsTer18	dn	Severe learning disabilities, irritablity
15098	female	6y	4y3m	c.6709C > T p.Arg2237Ter	dn	Autism, severe sleep disturbance
15966	male	8y	6y9m	c.5013_5041dup p.Gly1681ValfsTer6	dn	AHDH, autism, psychological problems
16169	male	14y	9y	c.2971_2972del p.Met991AspfsTer29	dn	Migraine headaches, sleep disturbance, mild intellectual disability
17296	female	9y	7y	c.2851-6_2851-3het_del	dn	ADHD, learning problem,
19654	female	7y7m	4y6m	c.3472G > C p.Asp1158His	paternal	Language developmental delay, learning disabilites
19917	female	7y	4y	c.2409 + 1G > C	N/A	autism, ADHD
19928	male	5y3m	3y6m	c.2180C > A p.Ser727Ter	dn	epilepsy, language development delay
20586	female	4y6m	3y	c.499_502het_del p.Cys167GlnfsTer10	maternal	autism, language development delay
21178	male	9m	6m	c.1019_1020del p.Ser340CysfsTer12	dn	epilepsy,Cerebral Hypoplasia,gross developmental delay
21265	male	8y	4y	c.1019_1020del p.Ser340CysfsTer12	dn	Abnormal brain MRI, intellectual disability
21432	female	15m	6m	c.6974_6977del p.Asp2325ValfsTer49	dn	epilepsy, gross developmental delay

3.2 Variant calling and filtering

GATK yielded an average of 325,476 variants among different samples that passed quality control measures, with a mean value of 1,056 rare variants (AF_eas < 0.01) per sample (Supplemental Fig. 1). After filtering using two criteria (Fig. 1), 439,148 rare LOF variants and 70,658 rare LP variants were considered in the SKAT-O analysis.

Whole exome sequencing was performed in 13 NF1 patients with neurological manifestation and 457 healthy controls. Rare variants were called using GATK, and filtered according to variant types or CADD phred. Finally, SKAT-O was performed leading to 130 and 126 potential modifier genes respectively, when applied to the LOF and LP variants. In total, 36 common potential modifier genes were identified (p-value < 0.005).

3.3 Enrichment of identified genes

SKAT-O statistics were calculated to determine whether enrichment of rare variants within putative LOF or LP variants was observed. A total of 130 and 127 genes that significantly modified neurological phenotype severity were detected, respectively (Supplemental Table 1). Overall, 36 genes (excluding NF1) predicted to significantly modify the neurological phenotype of NF1 were identified in both groups. All identified genes are expressed in the nervous system, and 19 of the 36 have previously been shown to be associated with nervous system development and functioning (Table 2). The association between rare LOF variants of potential modifier genes (≥ 5 unique variants) in case and control groups was clearly identified (Fig. 2). GO term analysis of identified potential modifier genes indicated that genes involved in neurological function and development were significantly enriched (Fig. 3, Supplemental Table 2). Interaction network evaluation using STRING and GeneMANIA showed close correlations and interactions between potential modifier genes.

Table 2

Identified candidate modifier genes of neurological involvement in NF1
Gene	Description	Associated monogenic phenotype (#OMIM)	expression in normal human nervous tissue	Association with neurological system	Pubmed ID
CUL7	Cullin 7	3-M syndrome 1	Y	Core component of a Cul7-RING ubiquitin-protein ligase, which mediates ubiquitination and consequent degradation of important proteins of nervous system.	21572988
DUSP26	Dual Specificity Phosphatase 26	None	Y	A member of the tyrosine phosphatase family of proteins, which regulates neuronal proliferation.	28701747
OR1L8	Olfactory Receptor Family 1 Subfamily L Member 8	None	Y	Olfactory receptors interact with odorant molecules in the nose, to initiate a neuronal response that triggers the perception of a smell.	14983052
HIGD1C	HIG1 Hypoxia Inducible Domain Family Member 1C	None	Y	None
SLC19A1	Solute Carrier Family 19 Member 1	?Megaloblastic anemia, folate-responsive	Y	A transporter of folate, which may play a role in the pathogenesis of pediatric brain tumor.	18406541
DPH1	Diphthamide Biosynthesis 1	Developmental delay with short stature, dysmorphic facial features, and sparse hair	Y	An enzyme involved in the biosynthesis of diphthamide, which is involved in the intellectual development.	25558065, 26220823
BCO1	Beta-Carotene Oxygenase 1	?Hypercarotenemia and vitamin A deficiency, autosomal dominant	Y	Play a role in regulating miRNAs associated with neuronal differentiation, cell-cell adhesion, and the Wnt signaling pathway	31048207
ZMIZ2	Zinc Finger MIZ-Type Containing 2	None	Y	None
GSTM3	Glutathione S-Transferase Mu 3	None	Y	Gene variations in GSTM3 are a risk factor for Alzheimer's disease.	18423940, 17904251
PRND	Prion Like Protein Doppel	None	Y	Polymorphisms within the PRND are associated with the Alzheimer's disease and Creutzfeldt-Jakob disease.	11702213, 14745079, 19422537
TWF2	Twinfilin Actin Binding Protein 2	None	Y	None
RASSF1	Ras Association Domain Family Member 1	None	Y	The inactivation of RASSF1A may be an important step in the tumorigenesis of the central nervous sytem.	14586413, 19165202, 17899687
ASCC1	Activating Signal Cointegrator 1 Complex Subunit 1	Barrett esophagus/esophageal adenocarcinoma Spinal muscular atrophy with congenital bone fractures 2	Y	May play a role in the development of neuromuscular junction.	28218388, 34204919
NFE2L2	NFE2 Like BZIP Transcription Factor 2	Immunodeficiency, developmental delay, and hypohomocysteinemia	Y	NFE2L2 promotes neuronal cell differentiation and controls neuroinflammation. NFE2L2 deregulation has been linked to neurodegenerative diseases.	19573594, 32545924, 29969760, 26541884
PARP6	Poly(ADP-Ribose) Polymerase Family Member 6	None	Y	A regulator of dendrite morphogenesis in hippocampal neurons	34067418, 26725726
LIPH	Lipase H	Hypotrichosis 7 Woolly hair, autosomal recessive 2 with or without hypotrichosis	Y	None
SF3A2	Splicing Factor 3a Subunit 2	None	Y	None
KIAA1522	KIAA1522	None	Y	None
LAMB1	Laminin Subunit Beta 1	Lissencephaly 5	Y	Involved in the organization of the laminar architecture of cerebral cortex.	23472759, 21370991
DYNC1LI1	Dynein Cytoplasmic 1 Light Intermediate Chain 1	None	Y	DYNC1I1 deficiency causes neuronal atrophy in primary hippocampal neurons.	27510948
MYOC	Myocilin	Glaucoma 1A, primary open angle	Y	Mediates myelination in the peripheral nervous system and plays a role in neurite outgrowth.	23897819, 12799138
U2SURP	U2 SnRNP Associated SURP Domain Containing	None	Y	None
DDX39B	DExD-Box Helicase 39B	None	Y	Association of DDX39B polymorphisms with Alzheimer's disease	18715507
ZMAT2	Zinc Finger Matrin-Type 2	None	Y	None
GFI1	Growth Factor Independent 1 Transcriptional Repressor	None	Y	Important for neuroendocrine differentiation and regulation of neurite outgrowth	15466176, 19026687
SNAPC5	Small Nuclear RNA Activating Complex Polypeptide 5	None	Y	None
RBPMS	RNA Binding Protein, MRNA Processing Factor	None	Y	A selective marker of ganglion cells in the mammalian retina	24318667
DBF4	DBF4 Zinc Finger	None	Y	None
DMAC2L	Distal Membrane Arm Assembly Component 2 Like	None	Y	None
HIST1H3E	H3 Clustered Histone 6	None	Y	None
FIS1	Fission, Mitochondrial 1	None	Y	Potential biomarkers for Alzheimer's disease	22340708
RIT1	Ras Like Without CAAX 1	Noonan syndrome 8	Y	RIT1 cooperates with nerve growth factor to promote neuronal development and regeneration.	15632082 12668729
TRAPPC5	Trafficking Protein Particle Complex Subunit 5	None	Y	None
CELF3	CUGBP Elav-Like Family Member 3	None	Y	None
LRRC3C	Leucine Rich Repeat Containing 3C	None	Y	None
GCLM	Glutamate-Cysteine Ligase Modifier Subunit	{Myocardial infarction, susceptibility to}	Y	None

Incomplete penetrance and phenotypic variability have hindered the clinical diagnosis and study of monogenic Mendelian disorders. Possible explanations for these phenomena are attributed to a combination of genetic, environmental, and lifestyle factors [15], which can make it challenging for genetics professionals to interpret a person’s family medical history and predict the risk of passing a genetic condition to future generations or developing specific phenotypes. At present, the early diagnosis of NF1 is very simple due to genetic testing; however, determining whether and when children will develop clinical phenotypes after early diagnosis requires additional research.

Most individuals with NF1 have normal intelligence, but learning disabilities or behavioral problems occur in 50–80% of patients [16]. Features of autism spectrum disorder occur in up to 30% of children with NF1 [17]. Various other learning and behavioral problems that persist into adulthood have been described [18]. The genetic burden related to neurological manifestations, such as intellectual disabilities, in patients with 22q11.2 deletion syndrome has been reported [19]. Evidence for the burden of de novo mutations and inherited rare single nucleotide variants in children with sensory processing dysfunction also indicates that genetic background likely influences neurological traits [20]. Hence, interpreting rare variants associated with neurological manifestations in NF1 patients with information generated from whole exome sequencing data for clinical diagnosis my help researchers identify mechanisms that underlie reduced penetrance and phenotypical variability in these patients.

Modifier genes and relevant variants are assumed to affect mechanisms other than those associated with loss or gain of function mutations, as observed in Mendelian genetic disease genes. These genes may affect disease progression and phenotype severity by influencing multiple processes such as gene expression, transcription, or translation. In this study, we explored the collective effects of rare variants as modifiers of neurological phenotypes of NF1. A total of 194 rare variants of 36 genes predicted to have damaging effects on phenotypes were enriched in NF1 patients with neurological phenotypes (Supplemental table). All identified genes are expressed in the nervous system, as indicated by various databases including BioGPS (http://biogps.org/) and Genotype-Tissue Expression (GTEx) (https://www.gtexportal.org/). Among the 36 genes, 19 have been associated with nervous system development and functioning in previous studies. Our analysis identified several interesting candidate modifier genes including CUL7 ( p = 8.91e-05 ), DPH1 ( p = 3.16e-04 ), and BCO1 ( p = 3.18e-04 ).

CUL7 functions as a core component of the Cul7-RING ubiquitin-protein ligase, which mediates ubiquitination and the subsequent degradation of important proteins of the nervous system [21]. Genes encoding the human Eag1 K + channel have been associated with congenital neurodevelopmental anomalies, and Cul7 is colocalized with Eag1 within synaptic regions of neurons. Cul7 also has been shown to facilitate protein degradation in a disease-associated rEag1 mutant, indicating its possible participation in neurological disorders [22]. DPH1 encodes an enzyme involved in the biosynthesis of diphthamide, which is involved in intellectual development [23]. Homozygous or compound heterozygous variants of DPH1 have been reported in patients with developmental delays, unusual skull shape with or without craniosynostosis, sparse hair, and facial dysmorphisms [24]. DPH1 is within the critical deletion region of Miller-Dieker syndrome (OMIM#247200) and DPH1 deficiency has been shown to contribute to craniofacial abnormalities [25]. BCO1 is the enzyme principally responsible for catalyzing the first step in vitamin A biosynthesis. Notably, hypercarotenemia and vitamin A deficiency are associated with mutations in this gene. It also plays a role in regulating miRNAs associated with neuronal differentiation, cell-cell adhesion, and the Wnt signaling pathway [26]. Among candidate modifier genes involved in neurological development and functioning, at least four genes are closely related to Alzheimer’s disease including FIS1, DDX39B, PRND, GSTM3. Polymorphisms in these genes are considered risk factors or potential biomarkers for Alzheimer disease [27, 28]. Enrichment of rare variants of these genes in NF1 patients with neurological manifestations suggests their potential function in neurodevelopmental aging in addition to neurodegeneration.

The current study had some limitations. The sample size used for assessing rare variants was limited, which is inevitable in similar studies focusing on rare diseases. To overcome this and increase statistical power from our relatively small cohort, we selected NF1 patients with similar phenotypes (neuropsychiatric symptoms) and screened putative LOF/LP rare variants using SKAT-O, rather than general rare variants. In the presence of positive correlations among the tests, the Bonferroni correction method is conservative and might lack power. For example, an identified disease-causing gene in NF1 patients reached a p-value of 5.2e-05 in our SKAT-O analysis which was above the Bonferroni-corrected significance threshold. Given the small sample size of our study, failing to pass the Bonferroni correction threshold could be expected, thus all genes with a p < 0.005 in the SKAT-O analysis were moved forward to functional follow-up. Setting the threshold at p < 0.005 for SKAT-O analysis could provide significant findings of rare variants, and many studies have already identified interesting signals using a higher p-value threshold in SKAT-O analysis [29–31]. Further, there was a lack of functional studies and modifier gene validations. Meanwhile, impacts of the accumulation of related rare genetic variants remain unknown. A detailed classification of neurological manifestations in a greater number of NF1 patients may help confirm the importance of identified modifier genes and rare variants. Further, long-term follow-up of diagnosed NF1 patients with and without neurological involvement will benefit in finding genes and variants related to late-onset symptoms.

In conclusion, our study revealed the possible role of rare variants in neurological manifestations of NF1. Analyses based on rare variant exome sequencing data of NF1 patients and normal controls identified 25 candidate modifier genes. More than half of these genes were found to be associated with normal nervous system development. Further validation of these candidate modifier genes in larger cohorts of NF1 patients is essential. The role of these genes in neurologic abnormalities of NF1 patients requires additional investigations. Interrogating possible modifier genes in genomic data of NF1 patients could provide auxiliary evidence for heterogeneity of the disease. Our study provides new insights into the pathophysiology of NF1 and may improve our ability to provide genetic counseling services to NF1 patients.

Funding: This work was supported by Shanghai "Rising Stars of Medical Talents" Youth Development Program and Clinical Research Plan of SHDC(SHDC2020CR6028).

Informed Consent Statement: The patients’ parents provided written informed consent.

Data Availability Statement: The data presented here are available on request from the corresponding author. The data are not publicly available due to privacy and ethical issues.

Acknowledgments: We would like to thank the patients’ family members for their participation.

Author Contributions: RY and JT contributed equally to the study. RY and TY contributed to study conception and design. NL and JW organized databases. JT and NL performed statistical analyses. RY and TY wrote the first draft of this manuscript. All authors revised, read, and approved the submitted version of the manuscript.

Conflicts of Interest: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.

Evans DG, Howard E, Giblin C et al (2010) Birth incidence and prevalence of tumor-prone syndromes: estimates from a UK family genetic register service. Am J Med Genet A 152A(2):327–332. doi: 10.1002/ajmg.a.33139
Kimberly J, Jan MF (2010) Clinical and genetic aspects of neurofibromatosis 1. Genet Med 12(1):1–11. doi: 10.1097/GIM.0b013e3181bf15e3
Margaret BR, David AS, David HV et al (2011) Variable expression of neurofibromatosis 1 in monozygotic twins. Am J Med Genet A 155A(3):478–485. doi: 10.1002/ajmg.a.33851
Emily RS, Teresa AS, Lisa JM et al (2017) Analysis of copy number variants in 11 pairs of monozygotic twins with neurofibroma-tosis type 1. Am J Med Genet A 173(3):647–653. doi: 10.1002/ajmg.a.38058
Dana MB, Sarah F, Stephan Z (2020) Genetic modifiers and non-Mendelian aspects of CMT. Brain Res 1726:146459. doi: 10.1016/j.brainres.2019.146459
Thomas WP, Adrian RK, Yimin H et al (2009) A positive modifier of spinal muscular atrophy in the SMN2 gene. Am J Hum Genet 85(3):408–413. doi: 10.1016/j.ajhg.2009.08.002
DePristo MA, Banks E, Poplin R et al (2011) A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet 43(5):491–498. doi:10.1038/ng.806
Talevich E, Shain AH, Botton T et al (2016) CNVkit: Genome-Wide Copy Number Detection and Visualization from Targeted DNA Sequencing. PLoS Comput Biol 12(4):e1004873 Published 2016 Apr 21. doi:10.1371/journal.pcbi.1004873
Lee S, Wu MC, Lin X (2012) Optimal tests for rare variant effects in sequencing association studies. Biostatistics 13(4):762–775. doi:10.1093/biostatistics/kxs014
Harris MA, Clark J, Ireland A et al (2004) The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res 32(Database issue):D258–D261. doi:10.1093/nar/gkh036
Kanehisa M, Goto S (2000) KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res 28(1):27–30. doi:10.1093/nar/28.1.27
Jassal B, Matthews L, Viteri G et al (2020) The reactome pathway knowledgebase. Nucleic Acids Res 48(D1):D498–D503. doi:10.1093/nar/gkz1031
Szklarczyk D, Gable AL, Lyon D et al STRING v11: protein-protein association networks with increased coverage, sup-porting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019 Jan 8;47(D1):D607-D613. doi: 10.1093/nar/gky1131. PMID: 30476243; PMCID: PMC6323986
Warde-Farley D, Donaldson SL, Comes O et al (2010) The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. Jul;38(Web Server issue):W214-20. doi: 10.1093/nar/gkq537. PMID: 20576703; PMCID: PMC2896186
Monica C, Alexandra PS, Jesus M et al (2017) Incomplete Penetrance and Variable Expressivity: Hallmarks in Channelopathies Associated with Sudden Cardiac Death. Biology (Basel) 7(1):3. doi: 10.3390/biology7010003
Annukka L, Emma H, Dorothy T et al (2013) Behaviour in children with neurofibromatosis type 1: cognition, executive function, attention, emotion, and social competence. Dev Med Child Neurol 55(2):111–125. doi: 10.1111/j.1469-8749.2012.04399.x
Stephanie MM, Maria TA, Shruti G et al (2016) Disease Burden and Symptom Structure of Autism in Neurofibromatosis Type 1: A Study of the International NF1-ASD Consortium Team (INFACT). AMA Psychiatry 73(12):1276–1284. doi: 10.1001/jamapsychiatry.2016.2600
Sofia G, Reinhard EF, Anna KL et al (2014) Influence of learning disabilities on the tumour predisposition syndrome NF1–survey from adult patients' perspective. Anticancer Res 34(7):3675–3681
Jensen M, Kooy RF, Simon TJ et al (2018) A higher rare CNV burden in the genetic background potentially contributes to intel-lectual disability phenotypes in 22q11.2 deletion syndrome. Eur J Med Genet 61(4):209–212. doi: 10.1016/j.ejmg.2017.11.016
Elysa JM, Anne BA, Vishnu PN et al (2018) Burden of de novo mutations and inherited rare single nucleotide variants in children with sensory processing dysfunction. BMC Med Genomics 11(1):50. doi: 10.1186/s12920-018-0362-x
Litterman N, Ikeuchi Y, Gallardo G et al (2011 May) An OBSL1-Cul7Fbxw8 ubiquitin ligase signaling mechanism regulates Golgi morphology and dendrite patterning. PLoS Biol 9(5):e1001060. doi: 10.1371/journal.pbio.1001060. Epub 2011 May 10. PMID: 21572988; PMCID: PMC3091842
Hsu PH, Ma YT, Fang YC et al Cullin 7 mediates proteasomal and lysosomal degradations of rat Eag1 potassium channels.Sci Rep. 2017 Jan18;7:40825. doi: 10.1038/srep40825. PMID: 28098200; PMCID: PMC5241692.
Alazami AM, Patel N, Shamseldin HE et al Accelerating novel candidate gene discovery in neurogenetic disorders via whole-exome sequencing of prescreened multiplex consanguineous families. Cell Rep. 2015 Jan 13;10(2):148 – 61. doi: 10.1016/j.celrep.2014.12.015. Epub 2014 Dec 31. PMID: 25558065
Cheng SSW, Luk HM, Lo IFM (2021) Jun;185(6):1925–1931 An adult Chinese patient with developmental delay with short stature, dysmorphic features, and sparse hair (Loucks-Innes syndrome). Am J Med Genet A. doi: 10.1002/ajmg.a.62164. Epub 2021 Mar 11. PMID: 33704902
Yu YR, You LR, Yan YT, Chen CM Role of OVCA1/DPH1 in craniofacial abnormalities of Miller-Dieker syndrome. Hum Mol Genet. 2014 Nov 1;23(21):5579-96. doi: 10.1093/hmg/ddu273. Epub 2014 Jun 3. PMID: 24895408
Kim YS, Gong X, Rubin LP et al (2019 Jul) β-Carotene 15,15'-oxygenase inhibits cancer cell stemness and metastasis by regulating differentiation-related miRNAs in human neuroblastoma. J Nutr Biochem 69:31–43. doi: 10.1016/j.jnutbio.2019.03.010. Epub 2019 Mar 28. PMID: 31048207
Wang S, Song J, Tan M et al (2012) Jul;19(7):1015-22 Mitochondrial fission proteins in peripheral blood lymphocytes are potential biomarkers for Alzheimer's disease. Eur J Neurol. doi: 10.1111/j.1468-1331.2012.03670.x. Epub 2012 Feb 16. PMID: 22340708
Gnjec A, D'Costa KJ, Laws SM et al Association of alleles carried at TNFA – 850 and BAT1 -22 with Alzheimer's disease.J Neuroinflammation. 2008 Aug20;5:36. doi: 10.1186/1742-2094-5-36. PMID: 18715507; PMCID: PMC2538517.
Cai C, Yin Z, Liu A et al (2022) Identifying Rare Genetic Variants of Immune Mediators as Risk Factors for Autism Spectrum Disorder. Genes (Basel). ;13(6):1098. Published 2022 Jun 20. doi:10.3390/genes13061098
De Vilder EYG, Martin L et al (2021) Rare Modifier Variants Alter the Severity of Cardiovascular Disease in Pseudoxanthoma Elasticum: Identification of Novel Candidate Modifier Genes and Disease Pathways Through Mixture of Effects Analysis. Front Cell Dev Biol 9:612581 Published 2021 Jun 8. doi:10.3389/fcell.2021.612581
Trevino CE, Rounds JC, Charen K et al (2021) Identifying susceptibility genes for primary ovarian insufficiency on the high-risk genetic background of a fragile X premutation. Fertil Steril 116(3):843–854. doi:10.1016/j.fertnstert.2021.04.021

Download PDF

Version 1

posted

You are reading this latest preprint version

Assessment of rare genetic variants to identify candidate modifier genes underlying neurological manifestations in neurofibromatosis 1 patients

Status:

Version 1

Abstract

Figures

1. Introduction

2. Methods

2.1 Patient enrollment

2.2 Sample preparation

2.3 Statistics

2.4 Pathway analysis

3. Results

3.1 Patient Characteristics

3.2 Variant calling and filtering

3.3 Enrichment of identified genes

4. Discussion

Declarations

References

Supplementary Files

Status:

Version 1