Activation of the IGF1 pathway mediates changes in cellular contractility and motility in single-suture craniosynostosis

Integrated transcriptomic analysis of human induced pluripotent stem cell-derived osteogenic differentiation reveals a regulatory role of KLF16

Osteogenic differentiation is essential for bone development, metabolism, and repair; however, the underlying regulatory relationships among genes remain poorly understood. To elucidate the transcriptomic changes and identify novel regulatory genes involved in osteogenic differentiation, we differentiated mesenchymal stem cells (MSCs) derived from 20 human iPSC lines into preosteoblasts (preOBs) and osteoblasts (OBs). We then performed transcriptome profiling of MSCs, preOBs and OBs. The iPSC-derived MSCs and OBs showed similar transcriptome profiles to those of primary human MSCs and OBs, respectively. Differential gene expression analysis revealed global changes in the transcriptomes from MSCs to preOBs, and then to OBs, including the differential expression of 840 genes encoding transcription factors (TFs). TF regulatory network analysis uncovered a network comprising 451 TFs, organized into five interactive modules. Multiscale embedded gene co-expression network analysis (MEGENA) identified gene co-expression modules and key network regulators (KNRs). From these analyses, KLF16 emerged as an important TF in osteogenic differentiation. We demonstrate that overexpression of Klf16 in vitro inhibited osteogenic differentiation and mineralization, while Klf16 +/- mice exhibited increased bone mineral density, trabecular number, and cortical bone area. Our study underscores the complexity of osteogenic differentiation and identifies novel regulatory genes such as KLF16, which plays an inhibitory role in osteogenic differentiation both in vitro and in vivo.

Loss of MMP9 disturbs cranial suture fusion via suppressing cell proliferation, chondrogenesis and osteogenesis in mice

Cranial sutures function as growth centers for calvarial bones. Abnormal suture closure will cause permanent cranium deformities. MMP9 is a member of the gelatinases that degrades components of the extracellular matrix. MMP9 has been reported to regulate bone development and remodeling. However, the function of MMP9 in cranial suture development is still unknown. Here, we identified that the expression of Mmp9 was specifically elevated during fusion of posterior frontal (PF) suture compared with other patent sutures in mice. Interestingly, inhibition of MMP9 ex vivo or knockout of Mmp9 in mice (Mmp9-/-) disturbed the fusion of PF suture. Histological analysis showed that knockout of Mmp9 resulted in wider distance between osteogenic fronts, suppressed cell condensation and endocranial bone formation in PF suture. Proliferation, chondrogenesis and osteogenesis of suture cells were decreased in Mmp9-/- mice, leading to the PF suture defects. Moreover, transcriptome analysis of PF suture revealed upregulated ribosome biogenesis and downregulated IGF signaling associated with abnormal closure of PF suture in Mmp9-/- mice. Inhibition of the ribosome biogenesis partially rescued PF suture defects caused by Mmp9 knockout. Altogether, these results indicate that MMP9 is critical for the fusion of cranial sutures, thus suggesting MMP9 as a potential therapeutic target for cranial suture diseases.

The value of genome-wide analysis in craniosynostosis

Background: This study assessed the diagnostic yield of high-throughput sequencing methods in a cohort of craniosynostosis (CS) patients not presenting causal variants identified through previous targeted analysis. Methods: Whole-genome or whole-exome sequencing (WGS/WES) was performed in a cohort of 59 patients (from 57 families) assessed by retrospective phenotyping as having syndromic or nonsyndromic CS. Results: A syndromic form was identified in 51% of the unrelated cases. A genetic cause was identified in 38% of syndromic cases, with novel variants detected in FGFR2 (a rare Alu insertion), TWIST1, TCF12, KIAA0586, HDAC9, FOXP1, and NSD2. Additionally, we report two patients with rare recurrent variants in KAT6A and YY1 as well as two patients with structural genomic aberrations: one with a 22q13 duplication and one with a complex rearrangement involving chromosome 2 (2p25 duplication including SOX11 and deletion of 2q22). Moreover, we identified potentially relevant variants in 87% of the remaining families with no previously detected causal variants, including novel variants in ADAMTSL4, ASH1L, ATRX, C2CD3, CHD5, ERF, H4C5, IFT122, IFT140, KDM6B, KMT2D, LTBP1, MAP3K7, NOTCH2, NSD1, SOS1, SPRY1, POLR2A, PRRX1, RECQL4, TAB2, TAOK1, TET3, TGFBR1, TCF20, and ZBTB20. Conclusion: These results confirm WGS/WES as a powerful diagnostic tool capable of either targeted in silico or broad genomic analysis depending on phenotypic presentation (e.g., classical or unusual forms of syndromic CS).

Transcriptional Regulation of Structural and Functional Adaptations in a Developing Adulthood Myocardium

The development of the heart follows a synergic action of several signaling pathways during gestational, pre- & postnatal stages. The current study aimed to investigate whether the myocardium experiences transcriptional changes during the transition from post-natal to adult hood stages. Herein, we used C57/B16/J mice at 4 (28- days; post-natal/PN) and 20 weeks (adulthood/AH) of ages and employed the next generation RNAseq (NGS) to profile the transcriptome and echocardiography analysis to monitor the structural/functional changes in the heart. NGS-based RNA-seq revealed that 1215 genes were significantly upregulated and 2549 were down regulated in the AH versus PN hearts, indicating a significant transcriptional change during this transition.

A synchronized cardiac transcriptional regulation through cell cycle, growth hormones, redox homeostasis and metabolic pathways was noticed in both PN and AH hearts. Echocardiography reveals significant structural and functional (i.e. systolic/diastolic) changes during the transition of PN to adult stage. Particularly, a progressive decline in ejection fraction and cardiac output was observed in AH hearts. These structural adaptations are in line with critical signaling pathways that drive the maturation of heart during AH. Overall, we have presented a comprehensive transcriptomic analysis along with structural-functional relationship during the myocardial development in adult mice.

IGF1R controls mechanosignaling in myofibroblasts required for pulmonary alveologenesis

Ventilation throughout life is dependent on the formation of pulmonary alveoli, which create an extensive surface area in which the close apposition of respiratory epithelium and endothelial cells of the pulmonary microvascular enables efficient gas exchange. Morphogenesis of the alveoli initiates at late gestation in humans and the early postnatal period in the mouse. Alveolar septation is directed by complex signaling interactions among multiple cell types. Here, we demonstrate that IGF1 receptor gene (Igf1r) expression by a subset of pulmonary fibroblasts is required for normal alveologenesis in mice. Postnatal deletion of Igf1r caused alveolar simplification, disrupting alveolar elastin networks and extracellular matrix without altering myofibroblast differentiation or proliferation. Moreover, loss of Igf1r impaired contractile properties of lung myofibroblasts and inhibited myosin light chain (MLC) phosphorylation and mechanotransductive nuclear YAP activity. Activation of p-AKT, p-MLC, and nuclear YAP in myofibroblasts was dependent on Igf1r. Pharmacologic activation of AKT enhanced MLC phosphorylation, increased YAP activation, and ameliorated alveolar simplification in vivo. IGF1R controls mechanosignaling in myofibroblasts required for lung alveologenesis.

Bioinformatics analysis of differentially expressed genes in subchondral bone in early experimental osteoarthritis using microarray data

Background

Osteoarthritis (OA) is the most common arthritic disease in humans, affecting the majority of individuals over 65 years of age. The aim of this study is to identify the gene expression profile specific to subchondral bone in OA by comparing the different expression profiles in experimental and sham-operation groups.

Methods

Gene expression profile GSE30322 was downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were obtained by limma package. And Database for Annotation, Visualization and Integrated Discovery (DAVID) databases were further used to identify the potential gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Furthermore, a protein–protein interaction (PPI) network was constructed and significant modules were extracted.

Results

Totally, 588 DEGs were identified including 199 upregulated DEGs and 389 downregulated DEGs screened in OA and sham-operation. GO showed that DEGs were significantly enhanced for ribosomal subunit export from nucleus and molting cycle. KEGG pathway analysis revealed that target genes were enriched in thiamine metabolism.

Conclusion

These key candidate DEGs that affect the progression of OA, and these genes might serve as potential therapeutic targets for OA.

Integrated Transcriptome and Network Analysis Reveals Spatiotemporal Dynamics of Calvarial Suturogenesis

Craniofacial abnormalities often involve sutures, the growth centers of the skull. To characterize the organization and processes governing their development, we profile the murine frontal suture, a model for sutural growth and fusion, at the tissue- and single-cell level on embryonic days (E)16.5 and E18.5. For the wild-type suture, bulk RNA sequencing (RNA-seq) analysis identifies mesenchyme-, osteogenic front-, and stage-enriched genes and biological processes, as well as alternative splicing events modifying the extracellular matrix. Single-cell RNA-seq analysis distinguishes multiple subpopulations, of which five define a mesenchyme-osteoblast differentiation trajectory and show variation along the anteroposterior axis. Similar analyses of in vivo mouse models of impaired frontal suturogenesis in Saethre-Chotzen and Apert syndromes, Twist1+/- and Fgfr2+/S252W, demonstrate distinct transcriptional changes involving angiogenesis and ribogenesis, respectively. Co-expression network analysis reveals gene expression modules from which we validate key driver genes regulating osteoblast differentiation. Our study provides a global approach to gain insights into suturogenesis.

(De)form and Function: Measuring Cellular Forces with Deformable Materials and Deformable Structures

The ability for biological cells to produce mechanical forces is important for the development, function, and homeostasis of tissue. The measurement of cellular forces is not a straightforward task because individual cells are microscopic in size and the forces they produce are at the nanonewton scale. Consequently, studies in cell mechanics rely on advanced biomaterials or flexible structures that permit one to infer these forces by the deformation they impart on the material or structure. Herein, the scientific progression on the use of deformable materials and deformable structures to measure cellular forces are reviewed. The findings and insights made possible with these approaches in the field of cell mechanics are summarized.

Transcriptome analysis reveals GPNMB as a potential therapeutic target for gastric cancer

Gastric cancer has the fifth highest incidence of disease and is the third leading cause of cancer-associated mortality in the world. The etiology of gastric cancer is complex and needs to be fully elucidated. Thus, it is necessary to explore potential pathogenic genes and pathways that contribute to gastric cancer. Gene expression profiles of the GSE33335 and GSE54129 datasets were downloaded from the Gene Expression Omnibus database. The differentially expressed genes (DEGs) were compared and identified using R software. The DEGs were then subjected to gene set enrichment analysis and Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. Survival analyses based on The Cancer Genome Atlas database were used to further screen the essential DEGs. A knockdown assay was performed to determine the function of the candidate gene in gastric cancer. Finally, the association between the candidate gene and immune-related genes was investigated. We found that GPNMB serves as an essential gene, with a high expression level, and predicts a worse outcome of gastric cancer. Knockdown of GPNMB inhibited gastric cancer cell proliferation and migration. In addition, GPNMB may augment the immunosuppressive ability of gastric cancer by recruiting immunosuppressive cells and promoting immune cell exhaustion through PI3K/AKT/CCL4 signaling axis. Collectively, these data suggest that GPNMB acts as an important positive mediator of tumor progression in gastric cancer, and GPNMB could exert multimodality modulation of gastric cancer-mediated immune suppression.

Calvarial osteoblast gene expression in patients with craniosynostosis leads to novel polygenic mouse model

Craniosynostosis is the premature fusion of the sutures of the calvaria and is principally designated as being either syndromic (demonstrating characteristic extracranial malformations) or non-syndromic. While many forms of syndromic craniosynostosis are known to be caused by specific mutations, the genetic etiology of non-syndromic, single-suture craniosynostosis (SSC) is poorly understood. Based on the low recurrence rate (4-7%) and the fact that recurrent mutations have not been identified for most cases of SSC, we propose that some cases of isolated, single suture craniosynostosis may be polygenic. Previous work in our lab identified a disproportionately high number of rare and novel gain-of-function IGF1R variants in patients with SSC as compared to controls. Building upon this result, we used expression array data from calvarial osteoblasts isolated from infants with and without SSC to ascertain correlations between high IGF1 expression and expression of other osteogenic genes of interest. We identified a positive correlation between increased expression of IGF1 and RUNX2, a gene known to cause SSC with increased gene dosage. Subsequent phosphorylation assays revealed that osteoblast cell lines from cases with high IGF1 expression demonstrated inhibition of GSK3β, a serine/threonine kinase known to inhibit RUNX2, thus activating osteogenesis through the IRS1-mediated Akt pathway. With these findings, we have utilized established mouse strains to examine a novel model of polygenic inheritance (a phenotype influenced by more than one gene) of SSC. Compound heterozygous mice with selective disinhibition of RUNX2 and either overexpression of IGF1 or loss of function of GSK3β demonstrated an increase in the frequency and severity of synostosis as compared to mice with the RUNX2 disinhibition alone. These polygenic mouse models reinforce, in-vivo, that the combination of activation of the IGF1 pathway and disinhibition of the RUNX2 pathway leads to an increased risk of developing craniosynostosis and serves as a model of human SSC.

Genetic Causes of Craniosynostosis: An Update

In 1993, Jabs et al. were the first to describe a genetic origin of craniosynostosis. Since this discovery, the genetic causes of the most common syndromes have been described. In 2015, a total of 57 human genes were reported for which there had been evidence that mutations were causally related to craniosynostosis. Facilitated by rapid technological developments, many others have been identified since then. Reviewing the literature, we characterize the most common craniosynostosis syndromes followed by a description of the novel causes that were identified between January 2015 and December 2017.

BBS9 gene in nonsyndromic craniosynostosis: Role of the primary cilium in the aberrant ossification of the suture osteogenic niche

Nonsyndromic craniosynostosis (NCS) is the premature ossification of skull sutures, without associated clinical features. Mutations in several genes account for a small number of NCS patients; thus, the molecular etiopathogenesis of NCS remains largely unclear. Our study aimed at characterizing the molecular signaling implicated in the aberrant ossification of sutures in NCS patients. Comparative gene expression profiling of NCS patient sutures identified a fused suture-specific signature, including 17 genes involved in primary cilium signaling and assembly. Cells from fused sutures displayed a reduced potential to form primary cilia compared to cells from control patent sutures of the same patient. We identified specific upregulated splice variants of the Bardet Biedl syndrome-associated gene 9 (BBS9), which encodes a structural component of the ciliary BBSome complex. BBS9 expression increased during in vitro osteogenic differentiation of suture-derived mesenchymal cells of NCS patients. Also, Bbs9 expression increased during in vivo ossification of rat sutures. BBS9 functional knockdown affected the expression of primary cilia on patient suture cells and their osteogenic potential. Computational modeling of the upregulated protein isoforms (observed in patients) predicted that their binding affinity within the BBSome may be affected, providing a possible explanation for the aberrant suture ossification in NCS.

Cell Mechanics of Craniosynostosis

Craniosynostosis is the premature fusion of the calvarial sutures that is associated with a number of physical and intellectual disabilities spanning from pediatric to adult years. Over the past two decades, techniques in molecular genetics and more recently, advances in high-throughput DNA sequencing have been used to examine the underlying pathogenesis of this disease. To date, mutations in 57 genes have been identified as causing craniosynostosis and the number of newly discovered genes is growing rapidly as a result of the advances in genomic technologies. While contributions from both genetic and environmental factors in this disease are increasingly apparent, there remains a gap in knowledge that bridges the clinical characteristics and genetic markers of craniosynostosis with their signaling pathways and mechanotransduction processes. By linking genotype to phenotype, outlining the role of cell mechanics may further uncover the specific mechanotransduction pathways underlying craniosynostosis. Here, we present a brief overview of the recent findings in craniofacial genetics and cell mechanics, discussing how this information together with animal models is advancing our understanding of craniofacial development.

Bioinformatics analysis of fibroblasts exposed to TGF‑β at the early proliferation phase of wound repair

The aim of the current study was to identify gene signatures during the early proliferation stage of wound repair and the effect of TGF-β on fibroblasts and reveal their potential mechanisms. The gene expression profiles of GSE79621 and GSE27165 were obtained from GEO database. Differentially expressed genes (DEGs) were identified using Morpheus and co-expressed DEGs were selected using Venn Diagram. Gene ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs were performed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool. Protein-protein interaction (PPI) networks of the DEGs were constructed using Cytoscape software. PPI interaction network was divided into subnetworks using the MCODE algorithm and the function of the top one module was analyzed using DAVID. The results revealed that upregulated DEGs were significantly enriched in biological process, including the Arp2/3 complex-mediated actin nucleation, positive regulation of hyaluronan cable assembly, purine nucleobase biosynthetic process, de novo inosine monophosphate biosynthetic process, positive regulation of epithelial cell proliferation, whereas the downregulated DEGs were enriched in the regulation of blood pressure, negative regulation of cell proliferation, ossification, negative regulation of gene expression and type I interferon signaling pathway. KEGG pathway analysis showed that the upregulated DEGs were enriched in shigellosis, pathogenic Escherichia coli infection, the mitogen-activated protein kinase signaling pathway, Ras signaling pathway and bacterial invasion of epithelial cells. The downregulated DEGs were enriched in systemic lupus erythematosus, lysosome, arachidonic acid metabolism, thyroid cancer and allograft rejection. The top 10 hub genes were identified from the PPI network. The top module analysis revealed that the included genes were involved in ion channel, neuroactive ligand-receptor interaction pathway, purine metabolism and intestinal immune network for IgA production pathway. The functional analysis revealed that TGF-β may promote fibroblast migration and proliferation and defend against microorganisms at the early proliferation stage of wound repair. Furthermore, these results may provide references for chronic wound repair.

Genetic advances in craniosynostosis

Craniosynostosis, the premature ossification of one or more skull sutures, is a clinically and genetically heterogeneous congenital anomaly affecting approximately one in 2,500 live births. In most cases, it occurs as an isolated congenital anomaly, that is, nonsyndromic craniosynostosis (NCS), the genetic, and environmental causes of which remain largely unknown. Recent data suggest that, at least some of the midline NCS cases may be explained by two loci inheritance. In approximately 25-30% of patients, craniosynostosis presents as a feature of a genetic syndrome due to chromosomal defects or mutations in genes within interconnected signaling pathways. The aim of this review is to provide a detailed and comprehensive update on the genetic and environmental factors associated with NCS, integrating the scientific findings achieved during the last decade. Focus on the neurodevelopmental, imaging, and treatment aspects of NCS is also provided.