Spatial transcriptomics implicates impaired BMP signaling in NF1 fracture pseudarthrosis in murine and patient tissues.

The neurofibromatosis type 1 (NF1) RASopathy is associated with persistent fibrotic nonunions (pseudarthrosis) in human and mouse skeletal tissue. Here, we performed spatial transcriptomics to define the molecular signatures occurring during normal endochondral healing following fracture in mice. Within the control fracture callus, we observed spatially restricted activation of morphogenetic pathways, such as TGF-ß, WNT, and BMP. To investigate the molecular mechanisms contributing to Nf1-deficient delayed fracture healing, we performed spatial transcriptomic analysis on a Postn-cre;Nf1fl/- (Nf1Postn) fracture callus. Transcriptional analyses, subsequently confirmed through phospho-SMAD1/5/8 immunohistochemistry, demonstrated a lack of BMP pathway induction in Nf1Postn mice. To gain further insight into the human condition, we performed spatial transcriptomic analysis of fracture pseudarthrosis tissue from a patient with NF1. Analyses detected increased MAPK signaling at the fibrocartilaginous-osseus junction. Similar to that in the Nf1Postn fracture, BMP pathway activation was absent within the pseudarthrosis tissue. Our results demonstrate the feasibility of delineating the molecular and tissue-specific heterogeneity inherent in complex regenerative processes, such as fracture healing, and reconstructing phase transitions representing endochondral bone formation in vivo. Furthermore, our results provide in situ molecular evidence of impaired BMP signaling underlying NF1 pseudarthrosis, potentially informing the clinical relevance of off-label BMP2 as a therapeutic intervention.

RAB1A haploinsufficiency phenocopies the 2p14-p15 microdeletion and is associated with impaired neuronal differentiation.

Hereditary spastic parapareses (HSPs) are clinically heterogeneous motor neuron diseases with variable age of onset and severity. Although variants in dozens of genes are implicated in HSPs, much of the genetic basis for pediatric-onset HSP remains unexplained. Here, we re-analyzed clinical exome-sequencing data from siblings with HSP of unknown genetic etiology and identified an inherited nonsense mutation (c.523C>T [p.Arg175Ter]) in the highly conserved RAB1A. The mutation is predicted to produce a truncated protein with an intact RAB GTPase domain but without two C-terminal cysteine residues required for proper subcellular protein localization. Additional RAB1A mutations, including two frameshift mutations and a mosaic missense mutation (c.83T>C [p.Leu28Pro]), were identified in three individuals with similar neurodevelopmental presentations. In rescue experiments, production of the full-length, but not the truncated, RAB1a rescued Golgi structure and cell proliferation in Rab1-depleted cells. In contrast, the missense-variant RAB1a disrupted Golgi structure despite intact Rab1 expression, suggesting a dominant-negative function of the mosaic missense mutation. Knock-down of RAB1A in cultured human embryonic stem cell-derived neurons resulted in impaired neuronal arborization. Finally, RAB1A is located within the 2p14-p15 microdeletion syndrome locus. The similar clinical presentations of individuals with RAB1A loss-of-function mutations and the 2p14-p15 microdeletion syndrome implicate loss of RAB1A in the pathogenesis of neurodevelopmental manifestations of this microdeletion syndrome. Our study identifies a RAB1A-related neurocognitive disorder with speech and motor delay, demonstrates an essential role for RAB1a in neuronal differentiation, and implicates RAB1A in the etiology of the neurodevelopmental sequelae associated with the 2p14-p15 microdeletion syndrome.

Molecular Dissection of Somatic Skeletal Disease in Neurofibromatosis Type 1.

Neurofibromatosis type 1 (NF1) is a tumor predisposition syndrome caused by heterozygous NF1 gene mutations. Patients with NF1 present with pleiotropic somatic secondary manifestations, including development of bone pseudarthrosis after fracture. Somatic NF1 gene mutations were reproducibly identified in patient-derived pseudarthrosis specimens, suggesting a local mosaic cell population including somatic pathologic cells. The somatic cellular pathogenesis of NF1 pseudarthroses remains unclear, though defects in osteogenesis have been posited. Here, we applied time-series single-cell RNA-sequencing (scRNA-seq) to patient-matched control and pseudarthrosis-derived primary bone stromal cells (BSCs). We show that osteogenic specification to an osteoblast progenitor cell population was evident for control bone-derived cells and haploinsufficient pseudarthrosis-derived cells. Similar results were observed for somatic patient fracture-derived NF1-/- cells; however, expression of genetic pathways associated with skeletal mineralization were significantly reduced in NF1-/- cells compared with fracture-derived NF1+/- cells. In mice, we show that Nf1 expressed in bone marrow osteoprogenitors is required for the maintenance of the adult skeleton. Results from our study implicate impaired Clec11a-Itga11-Wnt signaling in the pathogenesis of NF1-associated skeletal disease. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

DNA Polymerase Epsilon Deficiency Causes IMAGe Syndrome with Variable Immunodeficiency.

During genome replication, polymerase epsilon (Pol ε) acts as the major leading-strand DNA polymerase. Here we report the identification of biallelic mutations in POLE, encoding the Pol ε catalytic subunit POLE1, in 15 individuals from 12 families. Phenotypically, these individuals had clinical features closely resembling IMAGe syndrome (intrauterine growth restriction [IUGR], metaphyseal dysplasia, adrenal hypoplasia congenita, and genitourinary anomalies in males), a disorder previously associated with gain-of-function mutations in CDKN1C. POLE1-deficient individuals also exhibited distinctive facial features and variable immune dysfunction with evidence of lymphocyte deficiency. All subjects shared the same intronic variant (c.1686+32C>G) as part of a common haplotype, in combination with different loss-of-function variants in trans. The intronic variant alters splicing, and together the biallelic mutations lead to cellular deficiency of Pol ε and delayed S-phase progression. In summary, we establish POLE as a second gene in which mutations cause IMAGe syndrome. These findings add to a growing list of disorders due to mutations in DNA replication genes that manifest growth restriction alongside adrenal dysfunction and/or immunodeficiency, consolidating these as replisome phenotypes and highlighting a need for future studies to understand the tissue-specific development roles of the encoded proteins.

Mutations Preventing Regulated Exon Skipping in MET Cause Osteofibrous Dysplasia.

The periosteum contributes to bone repair and maintenance of cortical bone mass. In contrast to the understanding of bone development within the epiphyseal growth plate, factors that regulate periosteal osteogenesis have not been studied as intensively. Osteofibrous dysplasia (OFD) is a congenital disorder of osteogenesis and is typically sporadic and characterized by radiolucent lesions affecting the cortical bone immediately under the periosteum of the tibia and fibula. We identified germline mutations in MET, encoding a receptor tyrosine kinase, that segregate with an autosomal-dominant form of OFD in three families and a mutation in a fourth affected subject from a simplex family and with bilateral disease. Mutations identified in all families with dominant inheritance and in the one simplex subject with bilateral disease abolished the splice inclusion of exon 14 in MET transcripts, which resulted in a MET receptor (MET(Δ14)) lacking a cytoplasmic juxtamembrane domain. Splice exclusion of this domain occurs during normal embryonic development, and forced induction of this exon-exclusion event retarded osteoblastic differentiation in vitro and inhibited bone-matrix mineralization. In an additional subject with unilateral OFD, we identified a somatic MET mutation, also affecting exon 14, that substituted a tyrosine residue critical for MET receptor turnover and, as in the case of the MET(Δ14) mutations, had a stabilizing effect on the mature protein. Taken together, these data show that aberrant MET regulation via the juxtamembrane domain subverts core MET receptor functions that regulate osteogenesis within cortical diaphyseal bone.

Neurofibromin deficiency-associated transcriptional dysregulation suggests a novel therapy for tibial pseudoarthrosis in NF1.

Neurofibromatosis type 1 (NF1) is an autosomal dominant disease caused by mutations in NF1. Among the earliest manifestations is tibial pseudoarthrosis and persistent nonunion after fracture. To further understand the pathogenesis of pseudoarthrosis and the underlying bone remodeling defect, pseudoarthrosis tissue and cells cultured from surgically resected pseudoarthrosis tissue from NF1 individuals were analyzed using whole-exome and whole-transcriptome sequencing as well as genomewide microarray analysis. Genomewide analysis identified multiple genetic mechanisms resulting in somatic biallelic NF1 inactivation; no other genes with recurring somatic mutations were identified. Gene expression profiling identified dysregulated pathways associated with neurofibromin deficiency, including phosphoinositide 3-kinase (PI3K) and mitogen-activated protein kinase (MAPK) signaling pathways. Unlike aggressive NF1-associated malignancies, tibial pseudoarthrosis tissue does not harbor a high frequency of somatic mutations in oncogenes or other tumor-suppressor genes, such as p53. However, gene expression profiling indicates that pseudoarthrosis tissue has a tumor-promoting transcriptional pattern, despite lacking tumorigenic somatic mutations. Significant overexpression of specific cancer-associated genes in pseudoarthrosis highlights a potential for receptor tyrosine kinase inhibitors to target neurofibromin-deficient pseudoarthrosis and promote proper bone remodeling and fracture healing.

Somatic gain-of-function mutations in PIK3CA in patients with macrodactyly.

Macrodactyly is a discrete congenital anomaly consisting of enlargement of all tissues localized to the terminal portions of a limb, typically within a 'nerve territory'. The classic terminology for this condition is 'lipofibromatous hamartoma of nerve' or Type I macrodactyly. The peripheral nerve, itself, is enlarged both in circumference and in length. It is not related to neurofibromatosis (NF1), nor is it associated with vascular malformations, such as in the recently reported CLOVES syndrome. The specific nerve pathophysiology in this form of macrodactyly has not been well described and a genetic etiology for this specific form of enlargement is unknown. To identify the genetic cause of macrodactyly, we used whole-exome sequencing to identify somatic mutations present in the affected nerve of a single patient. We confirmed a novel mutation in PIK3CA (R115P) present in the patient's affected nerve tissue but not in blood DNA. Sequencing PIK3CA exons identified gain-of-function mutations (E542K, H1047L or H1047R) in the affected tissue of five additional unrelated patients; mutations were absent in blood DNA available from three patients. Immunocytochemistry confirmed AKT activation in cultured cells from the nerve of a macrodactyly patient. Additionally, we found that the most abnormal structure within the involved nerve in a macrodactylous digit is the perineurium, with additional secondary effects on the axon number and size. Thus, isolated congenital macrodactyly is caused by somatic activation of the PI3K/AKT cell-signaling pathway and is genetically and biochemically related to other overgrowth syndromes.

Total RNA isolation from stallion sperm and testis biopsies.

Sperm mRNA transcriptional profiles can be used to evaluate male fertility, yet differences between species in sperm attributes and packaging require adjustments in sperm RNA isolation protocols. The objective was to optimize RNA isolation methodology for fresh, frozen, and extended ejaculates, and epididymal sperm of stallions. Additionally, a protocol for RNA isolation from testis biopsies was established. Separation of sperm from somatic cells was critical for assuring the isolation of sperm-specific RNA. The highest purity was obtained by centrifuging ejaculates and epididymal sperm at 200 x g for 30 min through a 40% Equipure silanized silica particle solution. Sperm RNA isolation was more efficient with TRIzol reagent than with a spin-column based method; it resulted in 2 microg of total RNA per 100 x 10(6) sperm. To evaluate RNA quantity and quality, we used a NanoDrop spectrophotometer and Agilent Bioanalyzer. A protocol for reverse transcriptase PCR with equine primers for PRM2 and PTPRC genes was developed to determine sperm RNA contamination with genomic DNA or RNA from somatic cells. By these methods, hybridization- and sequencing-quality RNA was isolated from 11 samples of stallion sperm. Stallion testis biopsy with a 14 gauge 22 mm deep biopsy needle yielded approximately 12 microg of good quality total RNA, and could serve as an alternative to excision surgery for sample procurement. Compared to RNA isolation from testis, the sperm required advanced processing and RNA quality control. The described methodologies provided a foundation to establish functional genomic studies of stallion fertility.