Characterization of a novel COL10A1 variant associated with Schmid-type metaphyseal chondrodysplasia and a literature review

Screening and Analysis of Core Genes for Osteoporosis Based on Bioinformatics Analysis and Machine Learning Algorithms

Objective

This study aimed to identify osteoporosis-related core genes using bioinformatics analysis and machine learning algorithms.

Methods

mRNA expression profiles of osteoporosis patients were obtained from the Gene Expression Profiles (GEO) database, with GEO35958 and GEO84500 used as training sets, and GEO35957 and GSE56116 as validation sets. Differential gene expression analysis was performed using the R software "limma" package. A weighted gene co-expression network analysis (WGCNA) was conducted to identify key modules and modular genes of osteoporosis. Kyoto Gene and Genome Encyclopedia (KEGG), Gene Ontology (GO), and gene set enrichment analysis (GSEA) were performed on the differentially expressed genes. LASSO, SVM-RFE, and RF machine learning algorithms were used to screen for core genes, which were subsequently validated in the validation set. Predicted microRNAs (miRNAs) from the core genes were also analyzed, and differential miRNAs were validated using quantitative real-time PCR (qPCR) experiments.

Results

A total of 1280 differentially expressed genes were identified. A disease key module and 215 module key genes were identified by WGCNA. Three core genes (ADAMTS5, COL10A1, KIAA0040) were screened by machine learning algorithms, and COL10A1 had high diagnostic value for osteoporosis. Four core miRNAs (has-miR-148a-3p, has-miR-195-3p, has-miR-148b-3p, has-miR-4531) were found by intersecting predicted miRNAs with differential miRNAs from the dataset (GSE64433, GSE74209). The qPCR experiments validated that the expression of has-miR-195-3p, has-miR-148b-3p, and has-miR-4531 was significantly increased in osteoporosis patients.

Conclusion

This study demonstrated the utility of bioinformatics analysis and machine learning algorithms in identifying core genes associated with osteoporosis.

Natural history and genetic spectrum of the Turkish metaphyseal dysplasia cohort, including rare types caused by biallelic COL10A1, COL2A1, and LBR variants

BACKGROUND

Metaphyseal chondrodysplasias are a heterogeneous group of diseases characterized by short and bowed long bones and metaphyseal abnormality. The aim of this study is to investigate the genetic etiology and prognostic findings in patients with metaphyseal dysplasia.

METHODS

Twenty-four Turkish patients were included in this study and 13 of them were followed for 2-21 years. COL10A1, RMRP sequencing and whole exome sequencing were performed.

RESULTS

Results: Seven heterozygous pathogenic variants in COL10A1 were detected in 17 patients with Schmid type metaphyseal chondrodysplasia(MCDS). The phenotype was more severe in patients with heterozygous missense variants (one in signal peptide domain at the N-terminus of the protein, the other, class-1 group mutation at NC1 domain) compared to the patients with truncating variants. Short stature and coxa vara deformity appeared after 3 and 5 years of age, respectively, while large femoral head resolved after the age of 13 years in MCDS group. Interestingly, one patient with severe phenotype also had a biallelic missense variant in NC1 domain of COL10A1. Three patients with biallelic mutations in RMRP had prenatal onset short stature with short limb, and typical findings of cartilage hair hypoplasia (CHH). While immunodeficiency or recurrent infections were not observed, resistant congenital anemia was detected in one. Biallelic mutation in LBR was described in a patient with prenatal onset short stature, short and curved limb and metaphyseal abnormalities. Unlike previously reported patients, this patient had ectodermal findings, similar to CHH. A biallelic COL2A1 mutation was also found in the patient with lower limb deformities and metaphyseal involvement without vertebral and epiphyseal changes.

CONCLUSION

Long-term clinical characteristics are presented in a metaphyseal dysplasia cohort, including rare types caused by biallelic COL10A1, COL2A1, and LBR variants. We also point out that the domains where mutations on COL10A1 take place are important in the genotype-phenotype relationship.

High Expression of COL10A1 Is an Independent Predictive Poor Prognostic Biomarker and Associated with Immune Infiltration in Advanced Gastric Cancer Microenvironment

Medical technology has become more and more sophisticated recently, which, however, fails to contribute to a better prognosis for patients suffering advanced gastric cancer (GC). Hence, new biomarkers specific to GC diagnosis and prognosis shall be identified urgently. This study screened differentially expressed genes (DEGs) between 375 GC samples and 32 paracancer tissue samples from TCGA datasets. The expression of Collagen type X alpha 1 (COL10A1) in GC was analyzed. The chi-square test assisted in analyzing the relevance of COL10A1 to the clinicopathologic characteristics. The Kaplan-Meier method helped to assess the survival curves and log-rank tests assisted in the examination of the differences. The Cox proportional hazard regression model served for analyzing the risk factors for GC. Then, we developed a nomogram that contained the COL10A1 expression and clinical information. Finally, how COL10A1 expression was associated with the immune infiltration was also evaluated. In this study, 7179 upregulated and 3771 downregulated genes were identified. Among them, COL10A1 expression was distinctly increased in GC specimens compared with nontumor specimens. High COL10A1 expression exhibited an obvious relation to tumor T and pathologic stage. ROC assays confirmed the diagnostic value of COL10A1 expression in screening GC samples from normal samples. Survival data displayed that patients with high COL10A1 expression exhibited a shorter OS and DSS than those with low COL10A1 expression. We obtained a predictive nomogram, which could better predict the COL10A1 expression by virtue of discrimination and calibration. The prognostic value of COL10A1 expression was further confirmed in GSE84426 datasets. Immune assays revealed that COL10A1 expression was associated with tumor-filtrating immune cells, like CD8 T cells, cytotoxic cells, DC, eosinophils, iDC, macrophages, mast cells, NK CD56dim cells, NK cells, pDC, T helper cells, Tem, Th1 cells, Th17 cells, and Treg. Overall, we firstly proved that COL10A1 may be a novel and valuable prognostic and diagnostic factor for GC patients. In addition, COL10A1 has potential to be an immune indicator in GC.

Cancer-Associated Fibroblasts Hinder Lung Squamous Cell Carcinoma Oxidative Stress-Induced Apoptosis via METTL3 Mediated m6A Methylation of COL10A1

Background. Cancer-associated fibroblasts (CAFs) within the tumor microenvironment are key players in tumorigenesis and tumor development. Nevertheless, the regulatory mechanisms of CAFs on lung squamous cell carcinoma- (LUSC-) associated remain poorly elucidated. Methods. The microarray dataset GSE22874, containing 30 specimens of primary culture of normal fibroblasts (NFs) and 8 specimens of cancer-associated fibroblasts (CAFs) samples derived from LUSC, was retrieved from the Gene Expression Omnibus (GEO) database and then calculated by using the R language (limma package) to identify differentially expressed genes (DEGs). CAF-conditioned medium (CAF-CM) was collected and used to culture LUSC cells, followed by assessment of cell proliferation, apoptosis, and oxidative stress levels by using CCK-8, annexin V-FITC/PI double staining and ELISA assays. Subsequently, COL10A1 was knocked down in CAFs to assess the role of COL10A1 in CAF regulation of LUSC behavior. Bioinformatics online analysis and MeRIP were applied to predict and test the m6A modification of COL10A1 mRNA and the regulatory relationship with METTL3. Rescue experiments were next performed to explore the effects of METTL3 and COL10A1 in CAFs on LUSC cell proliferation, apoptosis, and oxidative stress. LUSC tumor cells with or without (COL10A1-silenced) CAFs were subcutaneously inoculated in nude mice to evaluate the effect of COL10A1 in CAFs on LUSC tumor growth. Results. Elevated expression of COL10A1 was found in LUSC-derived CAFs by GSE22874 dataset analysis. We discovered that COL10A1 and METTL3 was expressed in both LUSC cells and matched CAFs, while COL10A1 expression was prominently higher in CAFs than in LUSC cells. CAF-CM memorably encouraged LUSC cell proliferation and suppressed apoptosis-induced oxidative stress, which was reversed by interfering with COL10A1 expression in CAFs, suggesting that COL10A1 might be secreted by CAFs into the culture medium to exert its effects inside LUSC cells. Global m6A modification was decreased in METTL3 knocked down CAFs. M6A modification, expression levels, and stability of COL10A1 mRNA were impaired upon METTL3 knockdown in CAFs. Overexpression of COL10A1 in CAFs partially reversed the effect of METTL3 knockdown on the malignant behavior of LUSC cells. In vivo studies confirmed that CAFs accelerated LUSC tumor growth, and this effect was counteracted by COL10A1 silencing. Conclusions. COL10A1 secreted by CAFs could facilitate LUSC cell proliferation and repress apoptosis-induced oxidative stress, and the mechanism was due to elevated expression mediated by METTL3 promoting its mRNA m6A modification, thereby accelerating tumor growth.

Cell-Dependent Pathogenic Roles of Filamin B in Different Skeletal Malformations

Mutations of filamin B (FLNB) gene can lead to a spectrum of autosomal skeletal malformations including spondylocarpotarsal syndrome (SCT), Larsen syndrome (LRS), type I atelosteogenesis (AO1), type III atelosteogenesis (AO3), and boomerang dysplasia (BD). Among them, LRS is milder while BD causes a more severe phenotype. However, the molecular mechanism underlying the differences in clinical phenotypes of different FLNB variants has not been fully determined. Here, we presented two patients suffering from autosomal dominant LRS and autosomal recessive vitamin D-dependent rickets type IA (VDDR-IA). Whole-exome sequencing revealed two novel missense variants in FLNB, c.4846A>G (p.T1616A) and c.7022T>G (p.I2341R), which are located in repeat 15 and 22 of filamin B, respectively. The expression of FLNB^I2341R in the muscle tissue from our LRS patient was remarkably increased. And in vitro studies showed that both variants led to a lack of filopodia and accumulation of the mutants in the perinuclear region in HEK293 cells. We also found that c.4846A>G (p.T1616A) and c.7022T>G (p.I2341R) regulated endochondral osteogenesis in different ways. c.4846A>G (p.T1616A) activated AKT pathways through inhibiting SHIP2, suppressed the Smad3 pathway, and impaired the expression of Runx2 in both Saos-2 and ATDC5 cells. c.7022T>G (p.I2341R) activated both AKT and Smad3 pathways and increased the expression of Runx2 in Saos-2 cells, while in ATDC5 cells it activated AKT pathways through inhibiting SHIP2, suppressed the Smad3 pathway, and reduced the expression of Runx2. Our study demonstrated the pathogenic mechanisms of two novel FLNB variants in two different clinical settings and proved that FLNB variants could not only directly cause skeletal malformations but also worsen skeletal symptoms in the setting of other skeletal diseases. Besides, FLNB variants differentially affect skeletal development which contributes to clinical heterogeneity of FLNB-related disorders.

Natural Mutations Affect Structure and Function of gC1q Domain of Otolin-1

Otolin-1 is a scaffold protein of otoliths and otoconia, calcium carbonate biominerals from the inner ear. It contains a gC1q domain responsible for trimerization and binding of Ca²⁺. Knowledge of a structure-function relationship of gC1q domain of otolin-1 is crucial for understanding the biology of balance sensing. Here, we show how natural variants alter the structure of gC1q otolin-1 and how Ca²⁺ are able to revert some effects of the mutations. We discovered that natural substitutions: R339S, R342W and R402P negatively affect the stability of apo-gC1q otolin-1, and that Q426R has a stabilizing effect. In the presence of Ca²⁺, R342W and Q426R were stabilized at higher Ca²⁺ concentrations than the wild-type form, and R402P was completely insensitive to Ca²⁺. The mutations affected the self-association of gC1q otolin-1 by inducing detrimental aggregation (R342W) or disabling the trimerization (R402P) of the protein. Our results indicate that the natural variants of gC1q otolin-1 may have a potential to cause pathological changes in otoconia and otoconial membrane, which could affect sensing of balance and increase the probability of occurrence of benign paroxysmal positional vertigo (BPPV).