TSC2 loss in neural progenitor cells suppresses translation of ASD/NDD-associated transcripts in an mTORC1- and MNK1/2-reversible fashion.

Tuberous sclerosis complex (TSC) is an inherited neurodevelopmental disorder (NDD) with frequent manifestations of epilepsy and autism spectrum disorder (ASD). TSC is caused by inactivating mutations in TSC1 or TSC2 tumor suppressor genes, with encoded proteins hamartin (TSC1) and tuberin (TSC2) forming a functional complex inhibiting mechanistic target of rapamycin complex 1 (mTORC1) signaling. This has led to treatment with allosteric mTORC1 inhibitor rapamycin analogs ("rapalogs") for TSC tumors; however, rapalogs are ineffective for treating neurodevelopmental manifestations. mTORC1 signaling controls protein synthesis by regulating formation of the eIF4F complex, with further modulation by MNK1/2 kinases via phosphorylation of the eIF4F subunit eIF4E. While both these pathways modulate translation, comparing their impact on transcriptome-wide mRNA translation, as well as effects of inhibiting these pathways in TSC has not been explored. Here, employing CRISPR-modified, isogenic TSC2 patient-derived neural progenitor cells (NPCs), we have examined transcriptome-wide changes in mRNA translation upon TSC2 loss. Our results reveal dysregulated translation in TSC2 -Null NPCs, which significantly overlaps with the translatome from TSC1 -Null NPCs. Interestingly, numerous non-monogenic ASD-, NDD-and epilepsy-associated genes identified in patients harboring putative loss-of-function mutations, were translationally suppressed in TSC2 -Null NPCs. Importantly, translation of these ASD- and NDD-associated genes was reversed upon inhibition of either mTORC1 or MNK1/2 signaling using RMC-6272 or eFT-508, respectively. This study establishes the importance of mTORC1-eIF4F- and MNK-eIF4E-sensitive mRNA translation in TSC, ASD and other neurodevelopmental disorders laying the groundwork for evaluating drugs in clinical development that target these pathways as a treatment strategy for these disorders.

Proteasomal pathway inhibition as a potential therapy for NF2-associated meningioma and schwannoma.

BACKGROUND: Neurofibromatosis 2 (NF2) is an inherited disorder caused by bi-allelic inactivation of the NF2 tumor suppressor gene. NF2-associated tumors, including schwannoma and meningioma, are resistant to chemotherapy, often recurring despite surgery and/or radiation, and have generally shown cytostatic response to signal transduction pathway inhibitors, highlighting the need for improved cytotoxic therapies. METHODS: Leveraging data from our previous high-throughput drug screening in NF2 preclinical models, we identified a class of compounds targeting the ubiquitin-proteasome pathway (UPP), and undertook studies using candidate UPP inhibitors, ixazomib/MLN9708, pevonedistat/MLN4924, and TAK-243/MLN7243. Employing human primary and immortalized meningioma (MN) cell lines, CRISPR-modified Schwann cells (SCs), and mouse Nf2-/- SCs, we performed dose response testing, flow cytometry-based Annexin V and cell cycle analyses, and RNA-sequencing to identify potential underlying mechanisms of apoptosis. In vivo efficacy was also assessed in orthotopic NF2-deficient meningioma and schwannoma tumor models. RESULTS: Testing of three UPP inhibitors demonstrated potent reduction in cell viability and induction of apoptosis for ixazomib or TAK-243, but not pevonedistat. In vitro analyses revealed that ixazomib or TAK-243 downregulates expression of c-KIT and PDGFRα, as well as the E3 ubiquitin ligase SKP2 while upregulating genes associated with endoplasmic reticulum stress-mediated activation of the unfolded protein response (UPR). In vivo treatment of mouse models revealed delayed tumor growth, suggesting a therapeutic potential. CONCLUSIONS: This study demonstrates the efficacy of proteasomal pathway inhibitors in meningioma and schwannoma preclinical models and lays the groundwork for use of these drugs as a promising novel treatment strategy for NF2 patients.

Brigatinib causes tumor shrinkage in both NF2-deficient meningioma and schwannoma through inhibition of multiple tyrosine kinases but not ALK.

Neurofibromatosis Type 2 (NF2) is an autosomal dominant genetic syndrome caused by mutations in the NF2 tumor suppressor gene resulting in multiple schwannomas and meningiomas. There are no FDA approved therapies for these tumors and their relentless progression results in high rates of morbidity and mortality. Through a combination of high throughput screens, preclinical in vivo modeling, and evaluation of the kinome en masse, we identified actionable drug targets and efficacious experimental therapeutics for the treatment of NF2 related schwannomas and meningiomas. These efforts identified brigatinib (ALUNBRIG®), an FDA-approved inhibitor of multiple tyrosine kinases including ALK, to be a potent inhibitor of tumor growth in established NF2 deficient xenograft meningiomas and a genetically engineered murine model of spontaneous NF2 schwannomas. Surprisingly, neither meningioma nor schwannoma cells express ALK. Instead, we demonstrate that brigatinib inhibited multiple tyrosine kinases, including EphA2, Fer and focal adhesion kinase 1 (FAK1). These data demonstrate the power of the de novo unbiased approach for drug discovery and represents a major step forward in the advancement of therapeutics for the treatment of NF2 related malignancies.

A deep learning approach to identify gene targets of a therapeutic for human splicing disorders.

Pre-mRNA splicing is a key controller of human gene expression. Disturbances in splicing due to mutation lead to dysregulated protein expression and contribute to a substantial fraction of human disease. Several classes of splicing modulator compounds (SMCs) have been recently identified and establish that pre-mRNA splicing represents a target for therapy. We describe herein the identification of BPN-15477, a SMC that restores correct splicing of ELP1 exon 20. Using transcriptome sequencing from treated fibroblast cells and a machine learning approach, we identify BPN-15477 responsive sequence signatures. We then leverage this model to discover 155 human disease genes harboring ClinVar mutations predicted to alter pre-mRNA splicing as targets for BPN-15477. Splicing assays confirm successful correction of splicing defects caused by mutations in CFTR, LIPA, MLH1 and MAPT. Subsequent validations in two disease-relevant cellular models demonstrate that BPN-15477 increases functional protein, confirming the clinical potential of our predictions.

Physiological Characterization and Transcriptomic Properties of GnRH Neurons Derived From Human Stem Cells.

Gonadotropin-releasing hormone (GnRH) neurons in the hypothalamus play a key role in the regulation of reproductive function. In this study, we sought an efficient method for generating GnRH neurons from human embryonic and induced pluripotent stem cells (hESC and hiPSC, respectively). First, we found that exposure of primitive neuroepithelial cells, rather than neuroprogenitor cells, to fibroblast growth factor 8 (FGF8), was more effective in generating GnRH neurons. Second, addition of kisspeptin to FGF8 further increased the efficiency rates of GnRH neurogeneration. Third, we generated a fluorescent marker mCherry labeled human embryonic GnRH cell line (mCh-hESC) using a CRISPR-Cas9 targeting approach. Fourth, we examined physiological characteristics of GnRH (mCh-hESC) neurons: similar to GnRH neurons in vivo, they released the GnRH peptide in a pulsatile manner at ~60 min intervals; GnRH release increased in response to high potassium, kisspeptin, estradiol, and neurokinin B challenges; and injection of depolarizing current induced action potentials. Finally, we characterized developmental changes in transcriptomes of GnRH neurons using hESC, hiPSC, and mCh-hESC. The developmental pattern of transcriptomes was remarkably similar among the 3 cell lines. Collectively, human stem cell-derived GnRH neurons will be an important tool for establishing disease models to understand diseases, such as idiopathic hypothalamic hypogonadism, and testing contraceptive drugs.

De novo variants in SIAH1, encoding an E3 ubiquitin ligase, are associated with developmental delay, hypotonia and dysmorphic features.

BACKGROUND: Ubiquitination has a central role in numerous biological processes, including cell development, stress responses and ageing. Perturbed ubiquitination has been implicated in human diseases ranging from cancer to neurodegenerative diseases. SIAH1 encodes a RING-type E3 ubiquitin ligase involved in protein ubiquitination. Among numerous other roles, SIAH1 regulates metabotropic glutamate receptor signalling and affects neural cell fate. Moreover, SIAH1 positively regulates Wnt signalling through ubiquitin-mediated degradation of Axin and accumulation of ß-catenin. METHODS: Trio exome sequencing followed by Sanger validation was undertaken in five individuals with syndromic developmental delay. Three-dimensional structural modelling was used to predict pathogenicity of affected residues. Wnt stimulatory activity was measured by luciferase reporter assays and Axin degradation assays in HEK293 cells transfected with wild-type and mutant SIAH1 expression plasmids. RESULTS: We report five unrelated individuals with shared features of developmental delay, infantile hypotonia, dysmorphic features and laryngomalacia, in whom exome sequencing identified de novo monoallelic variants in SIAH1. In silico protein modelling suggested alteration of conserved functional sites. In vitro experiments demonstrated loss of Wnt stimulatory activity with the SIAH1 mutants, suggesting variant pathogenicity. CONCLUSION: Our results lend support to SIAH1 as a candidate Mendelian disease gene for a recognisable syndrome, further strengthening the connection between SIAH1 and neurodevelopmental disorders. Furthermore, the results suggest that dysregulation of the Wnt/ß-catenin pathway may be involved in the pathogenesis.

mTOR kinase inhibition disrupts neuregulin 1-ERBB3 autocrine signaling and sensitizes NF2-deficient meningioma cellular models to IGF1R inhibition.

Meningiomas (MNs), arising from the arachnoid/meningeal layer, are nonresponsive to chemotherapies, with ∼50% showing loss of the Neurofibromatosis 2 (NF2) tumor suppressor gene. Previously, we established NF2 loss activates mechanistic target of rapamycin complex 1 (mTORC1) and mechanistic target of rapamycin complex 2 (mTORC2) signaling, leading to clinical trials for NF2 and MN. Recently our omics studies identified activated ephrin (EPH) receptor and Src family kinases upon NF2 loss. Here, we report increased expression of several ligands in NF2-null human arachnoidal cells (ACs) and the MN cell line Ben-Men-1, particularly neuregulin-1/heregulin (NRG1), and confirm increased NRG1 secretion and activation of V-ERB-B avian erythroblastic leukemia viral oncogene homolog 3 (ERBB3) receptor kinase. Conditioned-medium from NF2-null ACs or exogenous NRG1 stimulated ERBB3, EPHA2, and mTORC1/2 signaling, suggesting pathway crosstalk. NF2-null cells treated with an ERBB3-neutralizing antibody partially downregulated mTOR pathway activation but showed no effect on viability. mTORC1/2 inhibitor treatment decreased NRG1 expression and downregulated ERBB3 while re-activating pAkt T308, suggesting a mechanism independent of NRG1-ERBB3 but likely involving activation of another upstream receptor kinase. Transcriptomics after mTORC1/2 inhibition confirmed decreased ERBB3/ERBB4 while revealing increased expression of insulin-like growth factor receptor 1 (IGF1R). Drug treatment co-targeting mTORC1/2 and IGF1R/insulin receptor attenuated pAkt T308 and showed synergistic effects on viability. Our findings indicate potential autocrine signaling where NF2 loss leads to secretion/activation of NRG1-ERBB3 signaling. mTORC1/2 inhibition downregulates NRG1-ERBB3, while upregulating pAkt T308 through an adaptive response involving IGF1R/insulin receptor and co-targeting these pathways may prove effective for treatment of NF2-deficient MN.

De novo variant in AMOTL1 in infant with cleft lip and palate, imperforate anus and dysmorphic features.

AMOTL1 belongs to the Motin family of proteins that are involved in organogenesis and tumorigenesis through regulation of cellular migration, tube formation, and angiogenesis. While involvement of all AMOTs in development or suppression of cancers is relatively well described, little is known about the congenital phenotype of pathogenic variants in these genes in humans. Recently, a heterozygous variant in AMOTL1 was published in association with orofacial clefts and cardiac abnormalities in an affected father and his daughter. However, studies in mice did not recapitulate the human phenotype and the case was summarized as inconclusive. We present a female infant with cleft lip and palate, imperforate anus and dysmorphic features, in whom trio exome sequencing revealed a de novo variant in AMOTL1 affecting a highly conserved amino acid (c.479C>T; p.[Pro160Leu]). Bioinformatic predictions and in silico modeling supported pathogenicity. This case reinforces the conjecture regarding the disruptive effect of pathogenic variants in AMOTL1 on organ formation in humans. Studies of additional families will reveal the full phenotypic spectrum associated with this multiple malformation syndrome.

Correlation between NF1 genotype and imaging phenotype on whole-body MRI: NF1 radiogenomics.

OBJECTIVE: To investigate the genotype-phenotype correlation between neurofibromatosis 1 (NF1) germline mutations and imaging features of neurofibromas on whole-body MRI (WBMRI) by using radiomics image analysis techniques. MATERIALS AND METHODS: Twenty-nine patients with NF1 who had known germline mutations determined by targeted next-generation sequencing were selected from a previous WBMRI study using coronal short tau inversion recovery sequence. Each tumor was segmented in WBMRI and a set of 59 imaging features was calculated using our in-house volumetric image analysis platform, 3DQI. A radiomics heatmap of 59 imaging features was analyzed to investigate the per-tumor and per-patient associations between the imaging features and mutation domains and mutation types. Linear mixed-effect models and one-way analysis of variance tests were performed to assess the similarity of tumor imaging features within mutation groups, between mutation groups, and between randomly selected groups. RESULTS: A total of 218 neurofibromas (97 discrete neurofibromas and 121 plexiform neurofibromas) were identified in 19 of the 29 patients. The unsupervised hierarchical clustering in heatmap analysis revealed 6 major image feature patterns that were significantly correlated with gene mutation domains and types with strong to very strong associations of genotype-phenotype correlations in both per-tumor and per-patient studies (p < 0.05, Cramer V > 0.5), whereas tumor size and locations showed no correlations with imaging features (p = 0.79 and p = 0.42, respectively). The statistical analyses revealed that the number of significantly different features (SDFs) within mutation groups were significantly lower than those between mutation groups (mutation domains: 10.9 ± 9.5% vs 31.9 ± 23.8% and mutation types: 31.8 ± 30.7% vs 52.6 ± 29.3%). The first and second quartile p values of within-patient groups were more than 2 times higher than those between-patient groups. However, the numbers of SDFs between randomly selected groups were much lower (approximately 5.2%). CONCLUSION: This preliminary study identified the NF1 radiogenomics linkage between NF1 causative mutations and MRI radiomic features, i.e., the correlation between NF1 genotype and imaging phenotype on WBMRI.

TSC patient-derived isogenic neural progenitor cells reveal altered early neurodevelopmental phenotypes and rapamycin-induced MNK-eIF4E signaling.

Background: Tuberous sclerosis complex (TSC) is a neurodevelopmental disorder with frequent occurrence of epilepsy, autism spectrum disorder (ASD), intellectual disability (ID), and tumors in multiple organs. The aberrant activation of mTORC1 in TSC has led to treatment with mTORC1 inhibitor rapamycin as a lifelong therapy for tumors, but TSC-associated neurocognitive manifestations remain unaffected by rapamycin. Methods: Here, we generated patient-specific, induced pluripotent stem cells (iPSCs) from a TSC patient with a heterozygous, germline, nonsense mutation in exon 15 of TSC1 and established an isogenic set of heterozygous (Het), null and corrected wildtype (Corr-WT) iPSCs using CRISPR/Cas9-mediated gene editing. We differentiated these iPSCs into neural progenitor cells (NPCs) and examined neurodevelopmental phenotypes, signaling and changes in gene expression by RNA-seq. Results: Differentiated NPCs revealed enlarged cell size in TSC1-Het and Null NPCs, consistent with mTORC1 activation. TSC1-Het and Null NPCs also revealed enhanced proliferation and altered neurite outgrowth in a genotype-dependent manner, which was not reversed by rapamycin. Transcriptome analyses of TSC1-NPCs revealed differentially expressed genes that display a genotype-dependent linear response, i.e., genes upregulated/downregulated in Het were further increased/decreased in Null. In particular, genes linked to ASD, epilepsy, and ID were significantly upregulated or downregulated warranting further investigation. In TSC1-Het and Null NPCs, we also observed basal activation of ERK1/2, which was further activated upon rapamycin treatment. Rapamycin also increased MNK1/2-eIF4E signaling in TSC1-deficient NPCs. Conclusion: MEK-ERK and MNK-eIF4E pathways regulate protein translation, and our results suggest that aberrant translation distinct in TSC1/2-deficient NPCs could play a role in neurodevelopmental defects. Our data showing upregulation of these signaling pathways by rapamycin support a strategy to combine a MEK or a MNK inhibitor with rapamycin that may be superior for TSC-associated CNS defects. Importantly, our generation of isogenic sets of NPCs from TSC patients provides a valuable platform for translatome and large-scale drug screening studies. Overall, our studies further support the notion that early developmental events such as NPC proliferation and initial process formation, such as neurite number and length that occur prior to neuronal differentiation, represent primary events in neurogenesis critical to disease pathogenesis of neurodevelopmental disorders such as ASD.

Traditional and systems biology based drug discovery for the rare tumor syndrome neurofibromatosis type 2.

Neurofibromatosis 2 (NF2) is a rare tumor suppressor syndrome that manifests with multiple schwannomas and meningiomas. There are no effective drug therapies for these benign tumors and conventional therapies have limited efficacy. Various model systems have been created and several drug targets have been implicated in NF2-driven tumorigenesis based on known effects of the absence of merlin, the product of the NF2 gene. We tested priority compounds based on known biology with traditional dose-concentration studies in meningioma and schwann cell systems. Concurrently, we studied functional kinome and gene expression in these cells pre- and post-treatment to determine merlin deficient molecular phenotypes. Cell viability results showed that three agents (GSK2126458, Panobinostat, CUDC-907) had the greatest activity across schwannoma and meningioma cell systems, but merlin status did not significantly influence response. In vivo, drug effect was tumor specific with meningioma, but not schwannoma, showing response to GSK2126458 and Panobinostat. In culture, changes in both the transcriptome and kinome in response to treatment clustered predominantly based on tumor type. However, there were differences in both gene expression and functional kinome at baseline between meningioma and schwannoma cell systems that may form the basis for future selective therapies. This work has created an openly accessible resource (www.synapse.org/SynodosNF2) of fully characterized isogenic schwannoma and meningioma cell systems as well as a rich data source of kinome and transcriptome data from these assay systems before and after treatment that enables single and combination drug discovery based on molecular phenotype.

Pain correlates with germline mutation in schwannomatosis.

Schwannomatosis has been linked to germline mutations in the SMARCB1 and LZTR1 genes, and is frequently associated with pain.In a cohort study, we assessed the mutation status of 37 patients with clinically diagnosed schwannomatosis and compared to clinical data, whole body MRI (WBMRI), visual analog pain scale, and Short Form 36 (SF-36) bodily pain subscale.We identified a germline mutation in LZTR1 in 5 patients (13.5%) and SMARCB1 in 15 patients (40.5%), but found no germline mutation in 17 patients (45.9%). Peripheral schwannomas were detected in 3 LZTR1-mutant (60%) and 10 SMARCB1-mutant subjects (66.7%). Among those with peripheral tumors, the median tumor number was 4 in the LZTR1 group (median total body tumor volume 30 cc) and 10 in the SMARCB1 group (median volume 85cc), (P=.2915 for tumor number and P = .2289 for volume). mutation was associated with an increased prevalence of spinal schwannomas (100% vs 41%, P = .0197). The median pain score was 3.9/10 in the LZTR1 group and 0.5/10 in the SMARCB1 group (P = .0414), and SF-36 pain-associated quality of life was significantly worse in the LZTR1 group (P = .0106). Pain scores correlated with total body tumor volume (rho = 0.32471, P = .0499), but not with number of tumors (rho = 0.23065, P = .1696).We found no significant difference in quantitative tumor burden between mutational groups, but spinal schwannomas were more common in LZTR1-mutant patients. Pain was significantly higher in LZTR1-mutant than in SMARCB1-mutant patients, though spinal tumor location did not significantly correlate with pain. This suggests a possible genetic association with schwannomatosis-associated pain.

Recurrent De Novo and Biallelic Variation of ATAD3A, Encoding a Mitochondrial Membrane Protein, Results in Distinct Neurological Syndromes.

ATPase family AAA-domain containing protein 3A (ATAD3A) is a nuclear-encoded mitochondrial membrane protein implicated in mitochondrial dynamics, nucleoid organization, protein translation, cell growth, and cholesterol metabolism. We identified a recurrent de novo ATAD3A c.1582C>T (p.Arg528Trp) variant by whole-exome sequencing (WES) in five unrelated individuals with a core phenotype of global developmental delay, hypotonia, optic atrophy, axonal neuropathy, and hypertrophic cardiomyopathy. We also describe two families with biallelic variants in ATAD3A, including a homozygous variant in two siblings, and biallelic ATAD3A deletions mediated by nonallelic homologous recombination (NAHR) between ATAD3A and gene family members ATAD3B and ATAD3C. Tissue-specific overexpression of borR534W, the Drosophila mutation homologous to the human c.1582C>T (p.Arg528Trp) variant, resulted in a dramatic decrease in mitochondrial content, aberrant mitochondrial morphology, and increased autophagy. Homozygous null bor larvae showed a significant decrease of mitochondria, while overexpression of borWT resulted in larger, elongated mitochondria. Finally, fibroblasts of an affected individual exhibited increased mitophagy. We conclude that the p.Arg528Trp variant functions through a dominant-negative mechanism that results in small mitochondria that trigger mitophagy, resulting in a reduction in mitochondrial content. ATAD3A variation represents an additional link between mitochondrial dynamics and recognizable neurological syndromes, as seen with MFN2, OPA1, DNM1L, and STAT2 mutations.

Whole-exome sequencing identifies homozygous GPR161 mutation in a family with pituitary stalk interruption syndrome.

CONTEXT: Pituitary stalk interruption syndrome (PSIS) is a rare, congenital anomaly of the pituitary gland characterized by pituitary gland insufficiency, thin or discontinuous pituitary stalk, anterior pituitary hypoplasia, and ectopic positioning of the posterior pituitary gland (neurohypophysis). The clinical presentation of patients with PSIS varies from isolated growth hormone (GH) deficiency to combined pituitary insufficiency and accompanying extrapituitary findings. Mutations in HESX1, LHX4, OTX2, SOX3, and PROKR2 have been associated with PSIS in less than 5% of cases; thus, the underlying genetic etiology for the vast majority of cases remains to be determined. OBJECTIVE: We applied whole-exome sequencing (WES) to a consanguineous family with two affected siblings who have pituitary gland insufficiency and radiographic findings of hypoplastic (thin) pituitary gland, empty sella, ectopic neurohypophysis, and interrupted pitiutary stalk-characteristic clinical diagnostic findings of PSIS. DESIGN AND PARTICIPANTS: WES was applied to two affected and one unaffected siblings. RESULTS: WES of two affected and one unaffected sibling revealed a unique homozygous missense mutation in GPR161, which encodes the orphan G protein-coupled receptor 161, a protein responsible for transducing extracellular signals across the plasma membrane into the cell. CONCLUSION: Mutations of GPR161 may be implicated as a potential novel cause of PSIS.

CHD8 regulates neurodevelopmental pathways associated with autism spectrum disorder in neural progenitors.

Truncating mutations of chromodomain helicase DNA-binding protein 8 (CHD8), and of many other genes with diverse functions, are strong-effect risk factors for autism spectrum disorder (ASD), suggesting multiple mechanisms of pathogenesis. We explored the transcriptional networks that CHD8 regulates in neural progenitor cells (NPCs) by reducing its expression and then integrating transcriptome sequencing (RNA sequencing) with genome-wide CHD8 binding (ChIP sequencing). Suppressing CHD8 to levels comparable with the loss of a single allele caused altered expression of 1,756 genes, 64.9% of which were up-regulated. CHD8 showed widespread binding to chromatin, with 7,324 replicated sites that marked 5,658 genes. Integration of these data suggests that a limited array of direct regulatory effects of CHD8 produced a much larger network of secondary expression changes. Genes indirectly down-regulated (i.e., without CHD8-binding sites) reflect pathways involved in brain development, including synapse formation, neuron differentiation, cell adhesion, and axon guidance, whereas CHD8-bound genes are strongly associated with chromatin modification and transcriptional regulation. Genes associated with ASD were strongly enriched among indirectly down-regulated loci (P < 10(-8)) and CHD8-bound genes (P = 0.0043), which align with previously identified coexpression modules during fetal development. We also find an intriguing enrichment of cancer-related gene sets among CHD8-bound genes (P < 10(-10)). In vivo suppression of chd8 in zebrafish produced macrocephaly comparable to that of humans with inactivating mutations. These data indicate that heterozygous disruption of CHD8 precipitates a network of gene-expression changes involved in neurodevelopmental pathways in which many ASD-associated genes may converge on shared mechanisms of pathogenesis.

Whole exome sequencing identifies three novel mutations in ANTXR1 in families with GAPO syndrome.

GAPO syndrome (OMIM#230740) is the acronym for growth retardation, alopecia, pseudoanodontia, and optic atrophy. About 35 cases have been reported, making it among one of the rarest recessive conditions. Distinctive craniofacial features including alopecia, rarefaction of eyebrows and eyelashes, frontal bossing, high forehead, mid-facial hypoplasia, hypertelorism, and thickened eyelids and lips make GAPO syndrome a clinically recognizable phenotype. While this genomic study was in progress mutations in ANTXR1 were reported to cause GAPO syndrome. In our study we performed whole exome sequencing (WES) for five affected individuals from three Turkish kindreds segregating the GAPO trait. Exome sequencing analysis identified three novel homozygous mutations including; one frame-shift (c.1220_1221insT; p.Ala408Cysfs*2), one splice site (c.411A>G; p.Gln137Gln), and one non-synonymous (c.1150G>A; p.Gly384Ser) mutation in the ANTXR1 gene. Our studies expand the allelic spectrum in this rare condition and potentially provide insight into the role of ANTXR1 in the regulation of the extracellular matrix.

NetWalker: a contextual network analysis tool for functional genomics.

BACKGROUND: Functional analyses of genomic data within the context of a priori biomolecular networks can give valuable mechanistic insights. However, such analyses are not a trivial task, owing to the complexity of biological networks and lack of computational methods for their effective integration with experimental data. RESULTS: We developed a software application suite, NetWalker, as a one-stop platform featuring a number of novel holistic (i.e. assesses the whole data distribution without requiring data cutoffs) data integration and analysis methods for network-based comparative interpretations of genome-scale data. The central analysis components, NetWalk and FunWalk, are novel random walk-based network analysis methods that provide unique analysis capabilities to assess the entire data distributions together with network connectivity to prioritize molecular and functional networks, respectively, most highlighted in the supplied data. Extensive inter-operability between the analysis components and with external applications, including R, adds to the flexibility of data analyses. Here, we present a detailed computational analysis of our microarray gene expression data from MCF7 cells treated with lethal and sublethal doses of doxorubicin. CONCLUSION: NetWalker, a detailed step-by-step tutorial containing the analyses presented in this paper and a manual are available at the web site http://netwalkersuite.org.