Search

ispartof: AMERICAN JOURNAL OF MEDICAL GENETICS PART A status: Published online

ObjectivesProgressive Pseudorheumatoid Arthropathy of Childhood (PPAC), caused by deficiency ofWNT1 inducible signaling pathway protein 3(WISP3), has been challenging to study because no animal model of the disease exists and cartilage recovered from affected patients is indistinguishable from common end-stage osteoarthritis. Therefore, to gain insights into why precocious articular cartilage failure occurs in this disease, we madein vitroderived articular cartilage using isogenicWISP3-deficient andWISP3-sufficient human pluripotent stem cells (hPSCs).MethodsWe generated articular cartilage-like tissues from induced-(i)PSCs from 2 patients with PPAC and 1 wild-type human embryonic stem cell line in which we knocked out WISP3. We compared these tissues toin vitro-derived articular cartilage tissues from 2 isogenicWISP3-sufficient control lines using histology, bulk RNA sequencing, single cell RNA sequencing, andin situhybridization.ResultsWISP3-deficient andWISP3-sufficient hPSCs both differentiated into articular cartilage-like tissues that appeared histologically similar. However, the transcriptomes ofWISP3-deficient tissues differed significantly fromWISP3-sufficient tissues and pointed to increased TGFβ, TNFα/NFkB, and IL-2/STAT5 signaling and decreased oxidative phosphorylation. Single cell sequencing andin situhybridization revealed thatWISP3-deficient cartilage contained a significantly higher fraction (∼ 4-fold increase,p< 0.001) of superficial zone chondrocytes compared to deeper zone chondrocytes than didWISP3-sufficient cartilage.ConclusionsWISP3-deficient andWISP3-sufficient hPSCs can be differentiated into articular cartilage-like tissues, but these tissues differ in their transcriptomes and in the relative abundances of chondrocyte sub-types they contain. These findings provide important starting points forin vivostudies when an animal model of PPAC or presymptomtic patient-derived articular cartilage becomes available.KEY MESSAGESWhat is already known on this topicLoss-of-function mutations inWISP3cause Progressive Pseudorheumatoid Arthropathy of Childhood (PPAC), yet the precise function ofWISP3in cartilage is unknown due to the absence of cartilage diseaseWisp3knockout mice and the lack of available PPAC patient cartilage that is not end-stage. Thus, most functional studies of WISP3 have been performedin vitrousing WISP3 over-expressing cell lines (i.e., not wild-type) andWISP3-deficient chondrocytes.What this study addsWe describe 3 newWISP3-deficient human pluripotent stem cell (hPSC) lines and show they can be differentiated into articular cartilage-like tissue.We comparein vitro-derived articular cartilage made fromWISP3-deficient and isogenicWISP3- sufficient hPSCs using bulk RNA sequencing, single cell RNA sequencing, andin situhybridization.We observe significant differences in the expression of genes previously associated with cartilage formation and homeostasis in the TGFβ, TNFα/NFkB, and IL-2/STAT5 signaling pathways. We also observe that WISP3-deficient cartilage-like tissues contain significantly higher fractions of chondrocytes that express superficial zone transcripts. These data suggest precocious cartilage failure in PPAC is the result of abnormal articular cartilage formation, dysregulated homeostatic signaling, or both.How this study might affect research, practice or policyThis study usesin vitro-derived articular cartilage to generate hypotheses for why cartilage fails in children with PPAC. This work prioritizes downstream studies to be performed when pre-symptomatic patient-derived cartilage samples or animal model of PPAC becomes available. It is essential to know how WISP3 functions in cartilage to develop therapies that benefit patients with PPAC and other degenerative joint diseases.

PURPOSE. In humans, loss-of-function mutations in the gene encoding Chordin-like 1 (CHRDL1) cause X-linked megalocornea (MGC1), characterized by bilateral corneal enlargement, decreased cornea thickness, and increased anterior chamber depth (ACD). We sought to determine whether Chrdl1 knockout (KO) mice would recapitulate the ocular findings found in patients with MGC1. METHODS. We generated mice with a Chrdl1 knockout allele and confirmed that male Chrdl1 hemizygous KO mice do not express Chrdl1 mRNA. We examined the eyes of male mice that were hemizygous for either the wild-type (WT) or KO allele and measured corneal diameter, corneal area, corneal thickness, endothelial cell density, ACD, tear volume, and intraocular pressure. We also harvested retinas and counted retinal ganglion cell numbers. Eye segregation pattern in the dorsal lateral geniculate nucleus were also compared between male Chrdl1 KO and WT mice. RESULTS. Male Chrdl1 KO mice do not have larger cornea diameters than WT mice. KO mice have significantly thicker central corneas (116.5 ± 3.9 vs 100.9 ± 4.2 μm, p = 0.013) and smaller ACD (325.7 ± 5.7 vs 405.6 ± 6.3 μm, p < 0.001) than WT mice, which is the converse of what occurs in patients who lack CHRDL1. Retinal- thalamic projections and other ocular measurements did not significantly differ between KO and WT mice. CONCLUSIONS. Male Chrdl1 KO mice do not have the same anterior chamber abnormalities seen in humans with CHRDL1 mutations. Therefore, Chrdl1 KO mice do not recapitulate the human MGC1 phenotype. Nevertheless, Chrdl1 plays a role during mouse ocular development, since corneas in KO mice differ from those in WT mice.

Serotonin (5HT) is a chemical messenger with biologic activities affecting multiple organs. Within the skeletal system, studies in mice and humans suggest blood 5HT levels affect bone, with elevations impairing and reductions enhancing bone accrual. Other studies, however, have not supported this hypothesis. Recently, administering 5HT to a Piezo1 mutant mouse strain with hyposerotonemia, intestinal dysmotility, and a doubling of femoral trabecular bone mass at 2 months of age, was reported to return the animals' intestinal motility and bone mass to normal. However, whether the 5HT directly affected bone metabolism or indirectly affected metabolism by improving intestinal function was not determined. Therefore, we administered 5HT to mice with normal intestinal function. We randomized female C57BL6/J mice and male and female mice that have increased bone mass due to a missense mutation in the WNT co-receptor LRP5 (Lrp5

The application of massively parallel sequencing technology to the field of skeletal disorders has boosted the discovery of the underlying genetic defect for many of these diseases. It has also resulted in the delineation of new clinical entities and the identification of genes and pathways that had not previously been associated with skeletal disorders. These rapid advances have prompted the Nosology Committee of the International Skeletal Dysplasia Society to revise and update the last (2015) version of the Nosology and Classification of Genetic Skeletal Disorders. This newest and tenth version of the Nosology comprises 461 different diseases that are classified into 42 groups based on their clinical, radiographic, and/or molecular phenotypes. Remarkably, pathogenic variants affecting 437 different genes have been found in 425/461 (92%) of these disorders. By providing a reference list of recognized entities and their causal genes, the Nosology should help clinicians achieve accurate diagnoses for their patients and help scientists advance research in skeletal biology.

Mechanical stimulation is a key regulator of bone mass, maintenance, and turnover. Wnt signaling is a key regulator of mechanotransduction in bone, but the role of β-catenin-an intracellular signaling node in the canonical Wnt pathway-in disuse mechanotransduction is not defined. Using the β-catenin exon 3 flox (constitutively active [CA]) mouse model, in conjunction with a tamoxifen-inducible, osteocyte-selective Cre driver, we evaluated the effects of degradation-resistant β-catenin on bone properties during disuse. We hypothesized that if β-catenin plays an important role in Wnt-mediated osteoprotection, then artificial stabilization of β-catenin in osteocytes would protect the limbs from disuse-induced bone wasting. Two disuse models were tested: tail suspension, which models fluid shift, and botulinum-toxin (botox)-induced muscle paralysis, which models loss of muscle force. Tail suspension was associated with a significant loss of tibial bone mass and density, reduced architectural properties, and decreased bone formation indices in uninduced (control) mice, as assessed by dual-energy X-ray absorptiometry (DXA), micro-computed tomography (µCT), and histomorphometry. Activation of the βcatCA allele in tail-suspended mice resulted in little to no change in those properties; ie, these mice were protected from bone loss. Similar protective effects were observed among botox-treated mice when the βcatCA was activated. RNAseq analysis of altered gene regulation in tail-suspended mice yielded 35 genes, including Wnt11, Gli1, Nell1, Gdf5, and Pgf, which were significantly differentially regulated between tail-suspended β-catenin stabilized mice and tail-suspended nonstabilized mice. Our findings indicate that selectively targeting/blocking of β-catenin degradation in bone cells could have therapeutic implications in mechanically induced bone disease. © 2019 American Society for Bone and Mineral Research.

MicroRNAs (miRNAs) are post-transcriptional regulators of gene expression. Heterozygous loss-of-function point mutations of miRNA genes are associated with several human congenital disorders1-5, but neomorphic (gain-of-new-function) mutations in miRNAs due to nucleotide substitutions have not been reported. Here we describe a neomorphic seed region mutation in the chondrocyte-specific, super-enhancer-associated MIR140 gene encoding microRNA-140 (miR-140) in a novel autosomal dominant human skeletal dysplasia. Mice with the corresponding single nucleotide substitution show skeletal abnormalities similar to those of the patients but distinct from those of miR-140-null mice6. This mutant miRNA gene yields abundant mutant miR-140-5p expression without miRNA-processing defects. In chondrocytes, the mutation causes widespread derepression of wild-type miR-140-5p targets and repression of mutant miR-140-5p targets, indicating that the mutation produces both loss-of-function and gain-of-function effects. Furthermore, the mutant miR-140-5p seed competes with the conserved RNA-binding protein Ybx1 for overlapping binding sites. This finding may explain the potent target repression and robust in vivo effect by this mutant miRNA even in the absence of evolutionary selection of miRNA-target RNA interactions, which contributes to the strong regulatory effects of conserved miRNAs7,8. Our study presents the first case of a pathogenic gain-of-function miRNA mutation and provides molecular insight into neomorphic actions of emerging and/or mutant miRNAs.

Summary Objective Growth of mandibular condylar cartilage (MCC) is associated with cell proliferation within the polymorphic cell layer and subsequent differentiation into chondrocytes that reside along the condylar surface and along the cartilage/subchondral bone interface. We examined whether cells in the polymorphic layer would proliferate and repopulate toxin-induced cell-depleted areas in MCCs of adult mice. Method We induced diphtheria toxin (DTA) expression (ROSA26l-s-lDTA) to cell-autonomously kill large fractions of MCC chondrocytes throughout the cartilage or along the articular cartilage surface with Aggrecan-CreERt2 (AcanCreERt2) or Lubricin-CreERt2 (Prg4CreERt2) Cre-recombinase-inducible mice, respectively. We examined MCCs from these mice shortly after cell killing or several months later with histology and confocal microscopy for evidence of chondrocyte proliferation and repopulation. Results AcanCreERt2-induced DTA expression killed an average of 53% MCC chondrocytes in adult mice after 1 week (39–66%, 95% confidence interval (CI)). Twelve weeks later, surviving chondrocytes had proliferated but not migrated to cell depleted areas. Prg4CreERt2-induced DTA expression killed an average of 24% surface chondrocytes in mice after 5 weeks (14–34% CI). After thirteen weeks there was 34% fewer surface chondrocytes (4–63% CI) in Prg4CreERt2 DTA-induced mice compared to controls. Conclusion In adult mice, after diphtheria toxin-mediated chondrocyte killing, cell depleted areas within MCC cartilage are not repopulated by new cells.

Evolution is replete with reuse of genes in different contexts, leading to multifunctional roles of signaling factors during development. Here, we explore osteoclast regulation during skeletal development through analysis of colony-stimulating factor 1 receptor (csf1r) function in the zebrafish. A primary role of Csf1r signaling is to regulate the proliferation, differentiation and function of myelomonocytic cells, including osteoclasts. We demonstrate the retention of two functional paralogues of csf1r in zebrafish. Mutant analysis indicates that the paralogues have shared, non-redundant roles in regulating osteoclast activity during the formation of the adult skeleton. csf1ra, however, has adopted unique roles in pigment cell patterning not seen in the second paralogue. We identify a unique noncoding element within csf1ra of fishes that is sufficient for controlling gene expression in pigment cells during development. As a role for Csf1r signaling in pigmentation is not observed in mammals or birds, it is likely that the overlapping roles of the two paralogues released functional constraints on csf1ra, allowing the signaling capacity of Csf1r to serve a novel function in the evolution of pigment pattern in fishes.

Evolution is replete with reuse of genes in different contexts, leading to multifunctional roles of signaling factors during development. Here, we explore osteoclast regulation during skeletal development through analysis of colony stimulating factor 1 receptor (csf1r) function in the zebrafish. A primary role of Csf1r signaling is to regulate the proliferation, differentiation and function of myelomonocytic cells, including osteoclasts. We demonstrate the retention of two functional paralogues of csf1r in zebrafish. Mutant analysis indicates that the paralogues have shared, non-redundant roles in regulating osteoclast activity during the formation of the adult skeleton. Csf1ra, however, has adopted unique roles in pigment cell patterning not seen in the second paralogue. We identify a unique noncoding element within csf1ra of fishes that is sufficient for controlling gene expression in pigment cells during development. As a role for Csf1r signaling in pigmentation is not observed in mammals or birds, it is likely that the overlapping roles of the two paralogues released functional constraints on csf1ra, allowing the signaling capacity of Csf1r to serve a novel function in the evolution of pigment pattern in fishes.

Congenital lipomatous overgrowth with vascular, epidermal, and skeletal (CLOVES) anomalies and Klippel-Trenaunay (KTS) syndromes are caused by somatic gain-of-function mutations in PIK3CA, encoding a catalytic subunit of phosphoinositide 3-kinase. Affected tissue is needed to find mutations, as mutant alleles are not detectable in blood. Because some patients with CLOVES develop Wilms tumor, we tested urine as a source of DNA for mutation detection. We extracted DNA from the urine of 17 and 24 individuals with CLOVES and KTS, respectively, and screened 5 common PIK3CA mutation hotspots using droplet digital polymerase chain reaction. Six of 17 CLOVES participants (35%) had mutant PIK3CA alleles in urine. Among 8 individuals in whom a mutation had been previously identified in affected tissue, 4 had the same mutant allele in the urine. One study participant with CLOVES had been treated for Wilms tumor. We detected the same PIK3CA mutation in her affected tissue, urine, and tumor, indicating Wilms tumors probably arise from PIK3CA mutant cells in patients with CLOVES. No urine sample from a participant with KTS had detectable PIK3CA mutations. We suggest that urine, which has the advantage of being collected non-invasively, is useful when searching for mutations in individuals with CLOVES syndrome.

Dominant negative mutations in CLCN7, which encodes a homodimeric chloride channel needed for matrix acidification by osteoclasts, cause Albers-Schonberg disease (also known as autosomal dominant osteopetrosis type 2). More than 25 different CLCN7 mutations have been identified in patients affected with Albers-Schonberg disease, but only one mutation (Clcn7G213R) has been introduced in mice to create an animal model of this disease. Here we describe a mouse with a different osteopetrosis-causing mutation (Clcn7F318L). Compared to Clcn7+/+ mice, 12-week-old Clcn7F318L/+ mice have significantly increased trabecular bone volume, consistent with Clcn7F318L acting as a dominant negative mutation. Clcn7F318L/F318L and Clcn7F318L/G213R mice die by 1month of age and resemble Clcn7 knockout mice, which indicate that p.F318L mutant protein is non-functional and p.F318L and p.G213R mutant proteins do not complement one another. Since it has been reported that treatment with interferon gamma (IFN-G) improves bone properties in Clcn7G213R/+ mice, we treated Clcn7F318L/+ mice with IFN-G and observed a decrease in osteoclast number and mineral apposition rate, but no overall improvement in bone properties. Our results suggest that the benefits of IFN-G therapy in patients with Albers-Schonberg disease may be mutation-specific.

Disuse-induced bone loss is seen following spinal cord injury, prolonged bed rest, and exposure to microgravity. We performed whole transcriptomic profiling of cortical bone using RNA sequencing (RNAseq) and RNA molecular barcoding (NanoString) on a hindlimb unloading (HLU) mouse model to identify genes whose mRNA transcript abundances change in response to disuse. Eleven-week old female C57BL/6 mice were exposed to ambulatory loading or HLU for 7 days (n = 8/group). Total RNA from marrow-flushed femoral cortical bone was analyzed on HiSeq and NanoString platforms. The expression of several previously reported genes associated with Wnt signaling and metabolism was altered by HLU. Furthermore, the increased abundance of transcripts, such as Pfkfb3 and Mss51, after HLU imply these genes also have roles in the cortical bone’s response to altered mechanical loading. Our study demonstrates that an unbiased approach to assess the whole transcriptomic profile of cortical bone can reveal previously unidentified mechanosensitive genes and may eventually lead to novel targets to prevent disuse-induced osteoporosis.

The murine Gnas (human GNAS) locus gives rise to Gαs and different splice variants thereof. The Gαs promoter is not methylated thus allowing biallelic expression in most tissues. In contrast, the alternative first Gnas/GNAS exons and their promoters undergo parent specific methylation, which limits transcription to the non-methylated allele. Pseudohypoparathyroidism type Ia (PHP1A) or type Ib (PHP1B) are caused by heterozygous maternal GNAS mutations suggesting that little or no Gαs is derived in some tissues from the non-mutated paternal GNAS thereby causing hormonal resistance. Previous data had indicated that Gαs is mainly derived from the maternal Gnas allele in brown adipose tissue (BAT) of newborn mice, yet it is biallelically expressed in adult BAT. This suggested that paternal Gαs expression is regulated by an unknown factor(s) that varies considerably with age. To extend these findings, we now used a strain-specific SNP in Gnas exon 11 (rs13460569) for evaluation of parent-specific Gαs expression through the densitometric quantification of BanII-digested RT-PCR products and digital droplet PCR (ddPCR). At all investigated ages, Gαs transcripts were derived in BAT predominantly from the maternal Gnas allele, while kidney and liver showed largely biallelic Gαs expression. Only low or undetectable levels of other paternally Gnas-derived transcripts were observed, making it unlikely that these are involved in regulating paternal Gαs expression. Our findings suggest that a cis-acting factor could be implicated in reducing paternal Gαs expression in BAT and presumably in proximal renal tubules, thereby causing PTH-resistance if the maternal GNAS/Gnas allele is mutated.

BackgroundFacial infiltrating lipomatosis (FIL) is a congenital disorder that causes overgrowth of one side of the face. The purpose of this study was to determine whether PIK3CA mutations are present in tissues outside of the subcutaneous adipose.MethodsFIL tissues from three patients were dissected to enrich for cells from skin, subcutaneous tissue, orbicularis oris muscle, buccal fat, zygomatic bone, and mucosal neuroma. Endothelial cells within the affected tissue also were enriched using CD31 microbeads. Laser capture microdissection on formalin-fixed paraffin-embedded histologic sections was performed to collect specific cell types. DNA was extracted from each tissue and cell type, and measured for the abundance of mutant PIK3CA alleles using droplet digital PCR.ResultsWe detected mutant PIK3CA alleles in every tissue and cell type tested from each overgrown face; frequencies ranged from 1.5 to 53%. There were fewer mutant endothelial cells compared with nonendothelial cells, and the stromal cell compartment had the highest frequency of mutant cells in each tissue.ConclusionsPIK3CA mutations are not restricted to a single tissue or cell type in FIL. Overgrowth in this condition is likely due to the mutation arising in a cell that contributes to several different facial structures during embryogenesis.

Arteriovenous malformation (AVM) is a fast-flow, congenital vascular anomaly that may arise anywhere in the body. AVMs typically progress, causing destruction of surrounding tissue and, sometimes, cardiac overload. AVMs are difficult to control; they often re-expand after embolization or resection, and pharmacologic therapy is unavailable. We studied extracranial AVMs in order to identify their biological basis. We performed whole-exome sequencing (WES) and whole-genome sequencing (WGS) on AVM tissue from affected individuals. Endothelial cells were separated from non-endothelial cells by immune-affinity purification. We used droplet digital PCR (ddPCR) to confirm mutations found by WES and WGS, to determine whether mutant alleles were enriched in endothelial or non-endothelial cells, and to screen additional AVM specimens. In seven of ten specimens, WES and WGS detected and ddPCR confirmed somatic mutations in mitogen activated protein kinase kinase 1 (MAP2K1), the gene that encodes MAP-extracellular signal-regulated kinase 1 (MEK1). Mutant alleles were enriched in endothelial cells and were not present in blood or saliva. 9 of 15 additional AVM specimens contained mutant MAP2K1 alleles. Mutations were missense or small in-frame deletions that affect amino acid residues within or adjacent to the protein’s negative regulatory domain. Several of these mutations have been found in cancers and shown to increase MEK1 activity. In summary, somatic mutations in MAP2K1 are a common cause of extracranial AVM. The likely mechanism is endothelial cell dysfunction due to increased MEK1 activity. MEK1 inhibitors, which are approved to treat several forms of cancer, are potential therapeutic agents for individuals with extracranial AVM.

The Clustered Regularly Interspaced Palindromic Repeats (CRISPR)/CRISPR-associated protein (Cas) pathway is revolutionizing biological research. Modifications to this primitive prokaryotic immune system now enable scientists to efficiently edit DNA or modulate gene expression in living eukaryotic cells and organisms. Thus, many laboratories can now perform important experiments that previously were considered scientifically risky or too costly. Here, we describe the components of the CRISPR/Cas system that have been engineered for use in eukaryotes. We also explain how this system can be used to genetically modify cell lines and model organisms, or regulate gene expression in order to search for new participants in biological pathways. © 2017 American Society for Bone and Mineral Research.

Capillary malformation is a cutaneous vascular anomaly that is present at birth, darkens over time, and can cause overgrowth of tissues beneath the stain. The lesion is caused by a somatic activating mutation in GNAQ. In a previous study, we were unable to identify a GNAQ mutation in patients with a capillary malformation involving an overgrown lower extremity. We hypothesized that mutations in GNA11 or GNA14, genes closely related to GNAQ, also may cause capillary malformations. Human capillary malformation tissue obtained from 8 patients that had tested negative for GNAQ mutations were studied. Lesions involved an extremity (n = 7) or trunk (n = 1). Droplet digital PCR (ddPCR) was used to detect GNA11 or GNA14 mutant cells (p.Arg183) in the specimens. Single molecule molecular inversion probe sequencing (smMIP-seq) was performed to search for other mutations in GNA11. Mutations were validated by subcloning and sequencing amplimers. We found a somatic GNA11 missense mutation (c.547C > T; p.Arg183Cys) in 3 patients with a diffuse capillary malformation of an extremity. Mutant allelic frequencies ranged from 0.3 to 5.0%. GNA11 or GNA14 mutations were not found in 5 affected tissues or in unaffected tissues (white blood cell DNA). GNA11 mutations are associated with extremity capillary malformations causing overgrowth. Pharmacotherapy that affects GNA11 signaling may prevent the progression of capillary malformations.

Articular cartilage has little regenerative capacity. Recently, genetic lineage tracing experiments have revealed chondrocyte progenitors at the articular surface. We further characterized these progenitors by using in vivo genetic approaches. Histone H2B–green fluorescent protein retention revealed that superficial cells divide more slowly than underlying articular chondrocytes. Clonal genetic tracing combined with immunohistochemistry revealed that superficial cells renew their number by symmetric division, express mesenchymal stem cell markers, and generate chondrocytes via both asymmetric and symmetric differentiation. Quantitative analysis of cellular kinetics, in combination with phosphotungstic acid–enhanced micro–computed tomography, showed that superficial cells generate chondrocytes and contribute to the growth and reshaping of articular cartilage. Furthermore, we found that cartilage renewal occurs as the progeny of superficial cells fully replace fetal chondrocytes during early postnatal life. Thus, superficial cells are self-renewing progenitors that are capable of maintaining their own population and fulfilling criteria of unipotent adult stem cells. Furthermore, the progeny of these cells reconstitute adult articular cartilage de novo, entirely substituting fetal chondrocytes.—Li, L., Newton, P. T., Bouderlique, T., Sejnohova, M., Zikmund, T., Kozhemyakina, E., Xie, M., Krivanek, J., Kaiser, J., Qian, H., Dyachuk, V., Lassar, A. B., Warman, M. L., Barenius, B., Adameyko, I., Chagin, A. S. Superficial cells are self-renewing chondrocyte progenitors, which form the articular cartilage in juvenile mice.

Single cell RNA-seq (scRNA-seq) is emerging as a powerful technology to examine transcriptomes of individual cells. We determined whether scRNA-seq could be used to detect the effect of environmental and pharmacologic perturbations on osteoblasts. We began with a commonly used in vitro system in which freshly isolated neonatal mouse calvarial cells are expanded and induced to produce a mineralized matrix. We used scRNA-seq to compare the relative cell type abundances and the transcriptomes of freshly isolated cells to those that had been cultured for 12 days in vitro. We observed that the percentage of macrophage-like cells increased from 6% in freshly isolated calvarial cells to 34% in cultured cells. We also found that Bglap transcripts were abundant in freshly isolated osteoblasts but nearly undetectable in the cultured calvarial cells. Thus, scRNA-seq revealed significant differences between heterogeneity of cells in vivo and in vitro. We next performed scRNA-seq on freshly recovered long bone endocortical cells from mice that received either vehicle or Sclerostin-neutralizing antibody for 1 week. Bone anabolism-associated transcripts were also not significantly increased in immature and mature osteoblasts recovered from Sclerostin-neutralizing antibody treated mice; this is likely a consequence of being underpowered to detect modest changes in gene expression, since only 7% of the sequenced endocortical cells were osteoblasts, and a limited portion of their transcriptomes were sampled. We conclude that scRNA-seq can detect changes in cell abundance, identity, and gene expression in skeletally derived cells. In order to detect modest changes in osteoblast gene expression at the single cell level in the appendicular skeleton, larger numbers of osteoblasts from endocortical bone are required.

Single-cell RNA sequencing (scRNA-Seq) is emerging as a powerful technology to examine transcriptomes of individual cells. We determined whether scRNA-Seq could be used to detect the effect of environmental and pharmacologic perturbations on osteoblasts. We began with a commonly used in vitro system in which freshly isolated neonatal mouse calvarial cells are expanded and induced to produce a mineralized matrix. We used scRNA-Seq to compare the relative cell type abundances and the transcriptomes of freshly isolated cells to those that had been cultured for 12 days in vitro. We observed that the percentage of macrophage-like cells increased from 6% in freshly isolated calvarial cells to 34% in cultured cells. We also found that Bglap transcripts were abundant in freshly isolated osteoblasts but nearly undetectable in the cultured calvarial cells. Thus, scRNA-Seq revealed significant differences between heterogeneity of cells in vivo and in vitro. We next performed scRNA-Seq on freshly recovered long bone endocortical cells from mice that received either vehicle or sclerostin-neutralizing antibody for 1 week. We were unable to detect significant changes in bone anabolism-associated transcripts in immature and mature osteoblasts recovered from mice treated with sclerostin-neutralizing antibody; this might be a consequence of being underpowered to detect modest changes in gene expression, because only 7% of the sequenced endocortical cells were osteoblasts and a limited portion of their transcriptomes were sampled. We conclude that scRNA-Seq can detect changes in cell abundance, identity, and gene expression in skeletally derived cells. In order to detect modest changes in osteoblast gene expression at the single-cell level in the appendicular skeleton, larger numbers of osteoblasts from endocortical bone are required. © 2020 American Society for Bone and Mineral Research.

The majority of extracranial arteriovenous malformations (AVMs) are caused by somatic mutations in MAP2K1. We report a somatic HRAS mutation in a patient who has a facial AVM associated with subcutaneous adipose overgrowth. We performed whole exome sequencing on DNA from the affected tissue and found a HRAS mutation (p.Thr58_Ala59delinsValLeuAspVal). Mutant allelic frequency was 5% in whole tissue and 31% in isolated endothelial cells (ECs); the mutation was not present in blood DNA or non-ECs. Somatic mutations in HRAS can cause AVM.

Anterior cruciate ligament (ACL) transection surgery in the minipig induces post-traumatic osteoarthritis (PTOA) in a pattern similar to that seen in human patients after ACL injury. Prior studies have reported the presence of cartilage matrix-degrading proteases, such as Matrix metalloproteinase-1 (MMP-1) and A disintegrin and metalloproteinase with thrombospondin motifs 4 (ADAMTS-4), in the synovial fluid of injured or arthritic joints; however, the tissue origin of these proteases is unknown. The objective of this study was to identify transcriptional processes activated in the synovium after surgical induction of PTOA with ACL transection, and to determine if processes associated with proteolysis were enriched in the synovium after ACL transection. Unilateral ACL transection was performed in adolescent Yucatan minipigs and synovium samples were collected at 1, 5, 9, and 14 days post-injury. Transcriptome-wide gene expression levels were determined using bulk RNA-Sequencing in the surgical animals and control animals with healthy knees. The greatest number of transcripts with significant changes was observed 1 day after injury. These changes were primarily associated with cellular proliferation, consistent with measurements of increased cellularity of the synovium at the two-week time point. At five to 14 days, the expression of transcripts relating to proteolysis and cartilage development was significantly enriched. While protease inhibitor-encoding transcripts (TIMP2, TIMP3) represented the largest fraction of protease-associated transcripts in the uninjured synovium, protease-encoding transcripts (including MMP1, MMP2, ADAMTS4) predominated after surgery. Cartilage development-associated transcripts that are typically not expressed by synovial cells, such as ACAN and COMP, were enriched in the synovium following ACL-transection. The upregulation in both catabolic processes (proteolysis) and anabolic processes (cartilage development) suggests that the synovium plays a complex, balancing role in the early response to PTOA induction.

Osteocalcin (OCN), the most abundant non-collagenous protein in the bone matrix, is reported to be a bone-derived endocrine hormone with wide-ranging effects on many aspects of physiology, including glucose metabolism and male fertility. Many of these observations were made using an OCN-deficient mouse allele (Osc-) in which the 2 OCN-encoding genes in mice, Bglap and Bglap2, were deleted in ES cells by homologous recombination. Here we describe mice with a new Bglap and Bglap2 double knockout (dko) allele (Bglap/2p.Pro25fs17Ter) that was generated by CRISPR/Cas9-mediated gene editing. Mice homozygous for this new allele do not express full length Bglap or Bglap2 mRNA and have no immunodetectable OCN in their plasma. FTIR imaging of cortical and trabecular bone in these homozygous knockout animals finds alterations in the crystal size and maturity of the bone mineral, hydroxyapatite, compared to wild-type littermates; however, μCT and 3-point bending tests do not find differences from wild-type littermates with respect to bone mass and strength. In contrast to the previously reported OCN-deficient mice with the Osc- allele, blood glucose levels and male fertility in the OCN-deficient mice with Bglap/2pPro25fs17Ter allele did not have significant differences from wild-type littermates. We cannot explain the absence of endocrine effects in mice with this new knockout allele. Potential explanations include effects of each mutated allele on the transcription of neighboring genes, and differences in genetic background and environment. So that our findings can be confirmed and extended by other interested investigators, we are donating this new Bglap and Bglap2 double knockout strain to The Jackson Laboratory for academic distribution.Author SummaryCells that make and maintain bone express proteins that function locally or systemically. The former proteins, such as type 1 collagen, affect the material properties of the skeleton while the latter proteins, such as fibroblast growth factor 23, enable the skeleton to communicate with other organ systems. Mutations that affect the functions of most bone cell expressed proteins cause diseases that have similar features in humans and other mammals, such as mice; for example, brittle bone diseases for type 1 collagen mutations and hypophosphatemic rickets for fibroblast growth factor 23 mutations.Our study focuses on another bone cell expressed protein, osteocalcin, which has been suggested to function locally to affect bone strength and systemically as hormone. Studies using osteocalcin knockout mice led other investigators to suggest endocrine roles for osteocalcin in regulating blood glucose levels, male fertility, muscle mass, brain development, behavior and cognition. We therefore decided to generate a new strain of osteocalcin knockout mice that could also be used to investigate these non-skeletal effects.To our surprise the osteocalcin knockout mice we created do not significantly differ from wild-type mice for the 3 phenotypes we examined: bone strength, blood glucose levels, and male fertility. Our data are consistent with findings from osteocalcin knockout rats, but inconsistent with data from the original osteocalcin knockout mice. Because we do not know why our new strain of osteocalcin knockout mice fails to recapitulate phenotypes previously reported for another knockout mouse stain, we have donated our mice to a public repository so that they can be easily obtained and studied in other academic laboratories.

Osteocalcin (OCN), the most abundant noncollagenous protein in the bone matrix, is reported to be a bone-derived endocrine hormone with wide-ranging effects on many aspects of physiology, including glucose metabolism and male fertility. Many of these observations were made using an OCN-deficient mouse allele (Osc–) in which the 2 OCN-encoding genes in mice, Bglap and Bglap2, were deleted in ES cells by homologous recombination. Here we describe mice with a new Bglap and Bglap2 double-knockout (dko) allele (Bglap/2p.Pro25fs17Ter) that was generated by CRISPR/Cas9-mediated gene editing. Mice homozygous for this new allele do not express full-length Bglap or Bglap2 mRNA and have no immunodetectable OCN in their serum. FTIR imaging of cortical bone in these homozygous knockout animals finds alterations in the collagen maturity and carbonate to phosphate ratio in the cortical bone, compared with wild-type littermates. However, μCT and 3-point bending tests do not find differences from wild-type littermates with respect to bone mass and strength. In contrast to the previously reported OCN-deficient mice with the Osc−allele, serum glucose levels and male fertility in the OCN-deficient mice with the Bglap/2pPro25fs17Ter allele did not have significant differences from wild-type littermates. We cannot explain the absence of endocrine effects in mice with this new knockout allele. Possible explanations include the effects of each mutated allele on the transcription of neighboring genes, or differences in genetic background and environment. So that our findings can be confirmed and extended by other interested investigators, we are donating this new Bglap and Bglap2 double-knockout strain to the Jackson Laboratories for academic distribution., Author summary Cells that make and maintain bone express proteins that function either locally or systemically. The former proteins, such as type 1 collagen, affect the material properties of the skeleton, while the latter, such as fibroblast growth factor 23, enable the skeleton to communicate with other organ systems. Mutations that affect the functions of most bone-cell-expressed proteins cause diseases that have similar features in humans and other mammals such as mice, for example, brittle bone diseases for type 1 collagen mutations and hypophosphatemic rickets for mutations in fibroblast growth factor 23. Our study focuses on another bone-cell-expressed protein, osteocalcin, which has been suggested to function locally to affect bone strength and systemically as a hormone. Studies using osteocalcin knockout mice led other investigators to suggest endocrine roles for osteocalcin in regulating blood glucose, male fertility, muscle mass, brain development, behavior, and cognition. We therefore decided to generate a new strain of osteocalcin knockout mice that could also be used to investigate these nonskeletal effects. To our surprise, the osteocalcin knockout mice we created did not significantly differ from wild-type mice for the three phenotypes we examined: bone strength, blood glucose, and male fertility. Our data are consistent with findings from osteocalcin knockout rats but are inconsistent with data from the original osteocalcin knockout mice. Because we do not know why our new strain fails to recapitulate the phenotypes previously reported for another knockout mouse stain, we have donated our mice to a public repository so that they can be easily obtained and studied in other academic laboratories.

Evolution is replete with reuse of genes in different contexts, leading to multifunctional roles of signaling factors during development. Here, we explore osteoclast regulation during skeletal development through analysis of colony-stimulating factor 1 receptor (

BACKGROUND: The term “intramuscular hemangioma capillary type” (IHCT) refers to a fast-flow vascular lesion that is classified as a tumor, although its phenotype overlaps with arteriovenous malformation (AVM). The purpose of this study was to identify somatic mutations in IHCT. METHODS: Affected tissue specimens were obtained during a clinically indicated procedure. The diagnosis of IHCT was based on history, physical examination, imaging, and histopathology. Because somatic mutations in cancer-associated genes can cause vascular malformations, we sequenced exons from 446 cancer-related genes in DNA from 7 IHCT specimens. We then performed mutation-specific droplet digital PCR (ddPCR) to independently test for the presence of a somatic mutation found by sequencing and to screen one additional IHCT sample. RESULTS: We detected somatic mutations in 6 of 8 IHCT specimens. Four specimens had a mutation in MAP2K1 (p.Q58_E62del, p.P105_I107delinsL, p.Q56P) and 2 specimens had mutations in KRAS (p.K5E and p.G12D, p.G12D and p.Q22R). Mutant allele frequencies detected by sequencing and confirmed by ddPCR ranged from 2 to 15%. CONCLUSIONS: IHCT lesions are phenotypically similar to AVMs and contain the same somatic MAP2K1 or KRAS mutations, suggesting that IHCT is on the AVM spectrum. We propose calling this lesion “intramuscular fast-flow vascular anomaly”.

Enhancing LRP5 signaling and inhibiting TGFβ signaling have each been reported to increase bone mass and improve bone strength in wild-type mice. Monotherapy targeting LRP5 signaling, or TGFβ signaling, also improved bone properties in mouse models of Osteogenesis Imperfecta (OI). We investigated whether additive or synergistic increases in bone properties would be attained if enhanced LRP5 signaling was combined with TGFβ inhibition. We crossed an Lrp5 high bone mass (HBM) allele (Lrp5(A214V)) into the Col1a2(G610C/+) mouse model of OI. At 6-weeks-of-age we began treating mice with an antibody that inhibits TGFβ1, β2, and β3 (mAb 1D11), or with an isotype-matched control antibody (mAb 13C4). At 12-weeks-old, we observed that combining enhanced LRP5 signaling with inhibited TGFβ signaling produced an additive effect on femoral and vertebral trabecular bone volumes, but not on cortical bone volumes. Although enhanced LRP5 signaling increased femur strength in a 3-point bending assay in Col1a2(G610C/+) mice, femur strength did not improve further with TGFβ inhibition. Neither enhanced LRP5 signaling nor TGFβ inhibition, alone or in combination, improved femur 3-point-bending post-yield displacement in Col1a2(G610C/+) mice. These pre-clinical studies indicate combination therapies that target LRP5 and TGFβ signaling should increase trabecular bone mass in patients with OI more than targeting either signaling pathway alone. Whether additive increases in trabecular bone mass will occur in, and clinically benefit, patients with OI needs to be determined.

Osteogenesis imperfecta (OI) comprises a genetically heterogeneous group of skeletal fragility diseases. Here, we report on five independent families with a progressively deforming type of OI, in whom we identified four homozygous truncation or frameshift mutations in MESD. Affected individuals had recurrent fractures and at least one had oligodontia. MESD encodes an endoplasmic reticulum (ER) chaperone protein for the canonical Wingless-related integration site (WNT) signaling receptors LRP5 and LRP6. Because complete absence of MESD causes embryonic lethality in mice, we hypothesized that the OI-associated mutations are hypomorphic alleles since these mutations occur downstream of the chaperone activity domain but upstream of ER-retention domain. This would be consistent with the clinical phenotypes of skeletal fragility and oligodontia in persons deficient for LRP5 and LRP6, respectively. When we expressed wild-type (WT) and mutant MESD in HEK293T cells, we detected WT MESD in cell lysate but not in conditioned medium, whereas the converse was true for mutant MESD. We observed that both WT and mutant MESD retained the ability to chaperone LRP5. Thus, OI-associated MESD mutations produce hypomorphic alleles whose failure to remain within the ER significantly reduces but does not completely eliminate LRP5 and LRP6 trafficking. Since these individuals have no eye abnormalities (which occur in individuals completely lacking LRP5) and have neither limb nor brain patterning defects (both of which occur in mice completely lacking LRP6), we infer that bone mass accrual and dental patterning are more sensitive to reduced canonical WNT signaling than are other developmental processes. Biologic agents that can increase LRP5 and LRP6-mediated WNT signaling could benefit individuals with MESD-associated OI.

The cysteine knot protein sclerostin is an osteocyte-derived secreted inhibitor of the Wnt co-receptors LRP5 and LRP6. LRP5 plays a dominant role in bone homeostasis, but we previously reported that Sost/sclerostin suppression significantly increased osteogenesis regardless of Lrp5 presence or absence. Those observations suggested that the bone forming effects of sclerostin inhibition can occur through Lrp6 (when Lrp5 is suppressed), or through other yet undiscovered mechanisms independent of Lrp5/6. To distinguish between these two possibilities, we generated mice with compound deletion of Lrp5 and Lrp6 selectively in bone, and treated them with sclerostin monoclonal antibody (Scl-mAb). All mice were homozygous flox for both Lrp5 and Lrp6 (Lrp5f/f; Lrp6f/f), and varied only in whether or not they carried the Dmp1-Cre transgene. Positive (Cre+) and negative (Cre−) mice were injected with Scl-mAb or vehicle from 4.5 to 14 weeks of age. Vehicle-treated Cre+ mice exhibited significantly reduced skeletal properties compared to vehicle-treated Cre− mice, as assessed by DXA, μCT, pQCT, and histology, indicating that Lrp5/6 deletions were effective and efficient. Scl-mAb treatment improved nearly every bone-related parameter among Cre− mice, but the same treatment in Cre+ mice resulted in little to no improvement in skeletal properties. For the few endpoints where Cre+ mice responded to Scl-mAb, it is likely that antibody-induced promotion of Wnt signaling occurred in cell types earlier in the mesenchymal/osteoblast differentiation pathway than the Dmp1-expressing stage. This latter conclusion was supported by changes in some histomorphometric parameters. In conclusion, unlike with the deletion of Lrp5 alone, the bone-selective late-stage co-deletion of Lrp5 and Lrp6 significantly impairs or completely nullifies the osteogenic action of Scl-mAb, and highlights a major role for both Lrp5 and Lrp6 in the mechanism of action for the bone-building effects of sclerostin antibody.

In humans, somatic activating mutations in PIK3CA are associated with skeletal overgrowth. In order to determine if activated PI3K signaling in bone cells causes overgrowth, we used Tg(BGLAP-Cre) and Tg(DMP1-Cre) mouse strains to somatically activate a disease-causing conditional Pik3ca allele (Pik3caH1047R) in osteoblasts and osteocytes. We observed Tg(BGLAP-Cre);Pik3caH1047R/+ offspring were born at the expected Mendelian frequency. However, these mice developed cutaneous lymphatic malformations and died before 7 weeks of age. In contrast, Tg(DMP1-Cre);Pik3caH1047R/+ offspring survived and had no cutaneous lymphatic malformations. Assuming that Cre-activity outside of the skeletal system accounted for the difference in phenotype between Tg(BGLAP-Cre);Pik3caH1047R/+ and Tg(DMP1-Cre);Pik3caH1047R/+ mice, we developed sensitive and specific droplet digital PCR (ddPCR) assays to search for and quantify rates of Tg(BGLAP-Cre)- and Tg(DMP1-Cre)-mediated recombination in non-skeletal tissues. We observed Tg(BGLAP-Cre)-mediated recombination in several tissues including skin, muscle, artery, and brain; two CNS locations, hippocampus and cerebellum, exhibited Cre-mediated recombination in >5% of cells. Tg(DMP1-Cre)-mediated recombination was also observed in muscle, artery, and brain. Although we cannot preclude that differences in phenotype between mice with Tg(BGLAP-Cre)- and Tg(DMP1-Cre)-mediated PIK3CA activation are due to Cre-recombination being induced at different stages of osteoblast differentiation, differences in recombination at non-skeletal sites are the more likely explanation. Since unanticipated sites of recombination can affect the interpretation of data from experiments involving conditional alleles, we recommend ddPCR as a good first step for assessing efficiency, leakiness, and off-targeting in experiments that employ Cre-mediated or Flp-mediated recombination.

The Col1a2(+/G610C) knock-in mouse, models osteogenesis imperfecta in a large old order Amish family (OOA) with type IV OI, caused by a G-to-T transversion at nucleotide 2098, which alters the gly-610 codon in the triple-helical domain of the α2(I) chain of type I collagen. Mineral and matrix properties of the long bones and vertebrae of male Col1a2(+/G610C) and their wild-type controls (Col1a2(+/+)), were characterized to gain insight into the role of α2-chain collagen mutations in mineralization. Additionally, we examined the rescuability of the composition by sclerostin inhibition initiated by crossing Col1a2(+/G610C) with an LRP(+/A214V) high bone mass allele. At age 10-days, vertebrae and tibia showed few alterations by micro-CT or Fourier transform infrared imaging (FTIRI). At 2-months-of-age, Col1a2(+/G610C) tibias had 13% fewer secondary trabeculae than Col1a2(+/+), these were thinner (11%) and more widely spaced (20%) than those of Col1a2(+/+) mice. Vertebrae of Col1a2(+/G610C) mice at 2-months also had lower bone volume fraction (38%), trabecular number (13%), thickness (13%) and connectivity density (32%) compared to Col1(a2+/+). The cortical bone of Col1a2(+/G610C) tibias at 2-months had 3% higher tissue mineral density compared to Col1a2(+/+); Col1a2(+/G610C) vertebrae had lower cortical thickness (29%), bone area (37%) and polar moment of inertia (38%) relative to Col1a2(+/+). FTIRI analysis, which provides information on bone chemical composition at ~7μm-spatial resolution, showed tibias at 10-days did not differ between genotypes. Comparing identical bone types in Col1a2(+/G610C) to Col1a2(+/+) at 2-months-of-age, tibias showed higher mineral-to-matrix ratio in trabeculae (17%) and cortices (31%). and in vertebral cortices (28%). Collagen maturity was 42% higher at 10-days-of-age in Col1a2(+/G610C) vertebral trabeculae and in 2-month tibial cortices (12%), vertebral trabeculae (42%) and vertebral cortices (12%). Higher acid-phosphate substitution was noted in 10-day-old trabecular bone in vertebrae (31%) and in 2-month old trabecular bone in both tibia (31%) and vertebrae (4%). There was also a 16% lower carbonate-to-phosphate ratio in vertebral trabeculae and a correspondingly higher (22%) carbonate-to-phosphate ratio in 2month-old vertebral cortices. At age 3-months-of-age, male femurs with both a Col1a2(+/G610C) allele and a Lrp5 high bone mass allele (Lrp5+/A214V) showed an improvement in bone composition, presenting higher trabecular carbonate-to-phosphate ratio (18%) and lower trabecular and cortical acid-phosphate substitutions (8% and 18%, respectively). Together, these results indicate that mutant collagen α2(I) chain affects both bone quantity and composition, and the usefulness of this model for studies of potential OI therapies such as anti-sclerostin treatments.

Background A somatic mutation in GNAQ (c.548G>A; p.R183Q), encoding Gαq, has been found in syndromic and sporadic capillary malformation tissue. However, the specific cell type containing the mutation is unknown. The purpose of this study was to determine which cells in capillary malformations have the GNAQ mutation. Methods Human capillary malformation tissue was obtained from 13 patients during a clinically indicated procedure. Droplet digital polymerase chain reaction, capable of detecting mutant allelic frequencies as low as 0.1 percent, was used to quantify the abundance of GNAQ mutant cells in capillary malformation tissue. Six specimens were fractionated by fluorescence-activated cell sorting into hematopoietic, endothelial, perivascular, and stromal cells. The frequency of GNAQ mutant cells in these populations was quantified by droplet digital polymerase chain reaction. Results Eight capillary malformations contained GNAQ p.R183Q mutant cells, two lesions had novel GNAQ mutations (p.R183L and p.R183G), and three capillary malformations did not have a detectable GNAQ p.R183 mutation. Mutant allelic frequencies ranged from 2 to 11 percent. Following fluorescence-activated cell sorting, the GNAQ mutation was found in the endothelial but not the platelet-derived growth factor receptor-β-positive cell population; mutant allelic frequencies were 3 to 43 percent. Conclusion Endothelial cells in capillary malformations are enriched for GNAQ mutations and are likely responsible for the pathophysiology underlying capillary malformation.

Enhancing LRP5 signaling and inhibiting TGFβ signaling have each been reported to increase bone mass and improve bone strength in wild-type mice. Monotherapy targeting LRP5 signaling, or TGFβ signaling, also improved bone properties in mouse models of Osteogenesis Imperfecta (OI). We investigated whether additive or synergistic increases in bone properties would be attained if enhanced LRP5 signaling was combined with TGFβ inhibition. We crossed an Lrp5 high bone mass (HBM) allele (Lrp5A214V) into the Col1a2G610C/+ mouse model of OI. At 6-weeks-of-age we began treating mice with an antibody that inhibits TGFβ1, β2, and β3 (mAb 1D11), or with an isotype-matched control antibody (mAb 13C4). At 12-weeks-old, we observed that combining enhanced LRP5 signaling with inhibited TGFβ signaling produced an additive effect on femoral and vertebral trabecular bone volumes, but not on cortical bone volumes. Although enhanced LRP5 signaling increased femur strength in a 3-point bending assay in Col1a2G610C/+ mice, femur strength did not improve further with TGFβ inhibition. Neither enhanced LRP5 signaling nor TGFβ inhibition, alone or in combination, improved femur 3-point-bending post-yield displacement in Col1a2G610C/+ mice. These pre-clinical studies indicate combination therapies that target LRP5 and TGFβ signaling should increase trabecular bone mass in patients with OI more than targeting either signaling pathway alone. Whether additive increases in trabecular bone mass will occur in, and clinically benefit, patients with OI needs to be determined.