Fam20c regulates the calpain proteolysis system through phosphorylating Calpasatatin to maintain cell homeostasis
Journal of Translational Medicine volume 21, Article number: 417 (2023)
The family with sequence similarity 20-member C (FAM20C) kinase, a Golgi casein kinase, which is responsible for phosphorylating the majority of the extracellular phosphoproteins within S-x-E/pS motifs, and is fundamentally associated with multiple biological processes to maintain cell proliferation, biomineralization, migration, adhesion, and phosphate homeostasis. In dissecting how FAM20C regulates downstream molecules and potential mechanisms, however, there are multiple target molecules of FAM20C, particularly many phenomena remain elusive, such as changes in cell-autonomous behaviors, incompatibility in genotypes and phenotypes, and others.
Here, assay for transposase-accessible chromatin using sequencing (ATAC-seq), RNA sequencing (RNA-seq), proteomics, and phosphoproteomics were performed in Fam20c-dificient osteoblasts and to facilitate an integrated analysis and determine the impact of chromatin accessibility, genomic expression, protein alterations, signaling pathway, and post translational modifcations.
By combining ATAC-seq and RNA-seq, we identified TCF4 and Wnt signaling pathway as the key regulators in Fam20c-dificient cells. Further, we showed Calpastatin/Calpain proteolysis system as a novel target axis for FAM20C to regulate cell migration and F-actin cytoskeleton by integrated analysis of proteomics and phosphoproteomics. Furthermore, Calpastatin/Calpain proteolysis system could negatively regulate the Wnt signaling pathway.
These observations implied that Fam20c knockout osteoblasts would cause cell homeostatic imbalance, involving changes in multiple signaling pathways in the conduction system.
The family with sequence similarity 20C (FAM20C), the bona fide “Golgi casein kinase”, answers for the phosphorylation of over 100 secreted phosphoproteins with a consensus motif Ser-x-Glu/pSer and is ubiquitously expressed in multiple tissues and bodily fluids [1,2,3,4]. Moreover, FAM20C figures prominently in a wide range of cellular processes, including phosphate metabolism, lipid homeostasis, wound healing, biomineralization, cardiac function, cell adhesion, and migration [1, 5,6,7,8]. These abundant and useful functions make FAM20C an important protein kinase. At the physiological and pathological level, FAM20C gene mutations or aberrant function of kinase lead to many diseases, such as Rains Syndrome, cancers, and other diseases. Notwithstanding, the exploration of the pivotal role of FAM20C in health and disease remains in an infant stage.
Aberrant FAM20C mutation causes lethal or non-lethal Raine syndrome in humans, which is clinically characterized by generalized osteosclerotic dysplasia and hypophosphatemia [9, 10]. To better understand the pathogenesis of Raine syndrome in humans, from work in animal and cellular models established by ablating Fam20c gene, diseases and phenotypes associated with Raine syndrome have been studied. Mice lacking Fam20c present the phenotype similar to human non-lethal Raine syndrome, including growth retardation, rickets/osteomalacia, defects in the growth plate, defective osteocytes differentiation, reduction of serum phosphate, and elevation of serum fibroblast growth factor (FGF23) . This has been partially attributed to elevated FGF23, which regulates serum phosphate levels, and its activity is regulated by FAM20C-mediated phosphorylation at Ser180 and O-glycosylation [6, 12]. Furthermore, Fam20c−/− mice develop abnormal dental morphological changes and differentiation defects, which can be caused by the suppressing of BMP signaling pathways [13,14,15]. Thus, considering Fam20c-deficient mice develop hard tissue dysplasia and FAM20C is highly expressed in bone and tooth tissues [2, 16, 17], studies have sought to identify the cellular changes in FAM20C ablation cells. Intriguingly, in Fam20c-deficient dental mesenchymal cells and osteoblasts, changes in gene expression of key factors related to osteogenesis and mineralization are consistent with Fam20c−/− mice, and impaired cellular mineralization could not be rescued by the addition of normal bone-derived extracellular matrix proteins [18, 19]. These findings suggest that FAM20C regulates osteoblast behaviors in a cell-autonomous manner. However, little is known regarding the mechanisms by which FAM20C influences cell autonomous behaviors.
A comprehensive understanding of human diseases requires exploration of their complexity and heterogeneity of the pathogenesis at the epigenetic, genomic, and proteomic levels [20, 21]. In the study of Raine syndrome regard, only 70 cases have been reported so far due to its extremely rare and low incidence rate , making it difficult to apply high-throughput sequencing technology in the research field of Raine syndrome. To investigate the mechanisms underlying the FAM20C-associated phenotype in Raine syndrome, some studies have been conducted to study the biological functions of FAM20C by using high-throughput sequencing technology. Zhang et al.  analyzed the proteins that interact with FAM20C in the endoplasmic reticulum (ER) and Golgi apparatus through immunoprecipitation and mass spectrometry, revealing that FAM20C could phosphorylate ER oxidoreductin 1α (Ero1α) to maintain the ER redox homeostasis. Besides, our previous research conducted a transcriptional analysis of Fam20c-deficient osteoblasts, indicating the cells undergo mesenchymal to epithelial transformation . Besides a single sequencing method, the combination of multi-omics analysis provides a more comprehensive picture linking “molecular information” to “phenotype” through functional and signaling networks [25,26,27]. Little is known, however, about multi-omics analysis of FAM20C remains to be determined, especially chromatin sequencing.
Here, we use established Fam20c-deficient cells in our previously published work  to conduct a comprehensive omics analysis. To probe the relationship between chromatin, transcriptional and protein changes response to Fam20c ablation in osteoblasts, we performed integrated omics analyses using an assay for transposase-accessible chromatin using sequencing (ATAC-seq), RNA sequencing (RNA-seq), proteomics, and phosphoproteomics. Fam20c ablation specifically regulates a subset of genes comprising Wnt signaling pathway molecules and TCF4 was highly enriched in the chromatin opened of Fam20c deficient osteoblasts, as revealed by ATAC-seq. We identified Calpastatin/Calpain proteolysis system as a novel target axis for FAM20C to regulate cell migration and F-actin cytoskeleton. Subsequently, we demonstrated that Calpain agonist treatment rescued Fam20c-deficient cells migration in vitro, and negatively manage the Wnt signaling pathway. Further, our results suggested that the homeostatic imbalance of Fam20c knockout osteoblasts may involve changes in multiple signaling pathways in the conduction system.
Materials and methods
Cell lines and cell culture
Immortalized mouse Fam20cf/f osteoblast cells (referred as “OB Fam20cf/f”) were purchased from Applied Biological Materials Inc. (Lot# T0038). Fam20c-deficient osteoblasts (referred as “OB Fam20cKO”) were established with infection of Cre recombinase lentivirus as previously published work described . OB Fam20cf/f and OB Fam20cKO were cultured at 37 °C with 5% CO2. All cells were cultured in α-Minimum Essential Medium (α-MEM) (Biosharp, BL306A) containing 10% fetal bovine serum (Gibco, #10091148) and 1% penicillin–streptomycin (Gibco, #15140122). Additionally, cells were passaged with 0.25% trypsin (Gibco, #12604013).
Cell samples for ATAC-seq were prepared as described in Current protocols in molecular biology by Buenrostro et al. . In brief, 50,000 cells were obtained from each sample, resuspended in cold PBS, and centrifugated at 500 g for 5 min at 4 °C. Cell pellets were washed once with 50 μl cold PBS on ice, followed by a 5 min centrifugation at 4 °C. Pelleted nuclei were taken for transposition reaction with Tn5 enzyme, and the DNA fragments after transposition were recovered with MinElute PCR Purification Kit (Qiagen). Accessible DNA was amplified by PCR with 1 X NEBNext High-Fidelity Pcr Master Mix (New England Biolabs, MA). Libraries were sequenced on a Novaseq 6000 (Illumina).
ATAC raw reads data from each sample was filtered using Trimmomatic software (Version 0.36) , after data filtering, the overall quality of the resulting clean reads was assessed. The clean reads were then aligned to reference genome mm10_gencode using the BWA program (Version 0.7.13-r1126). Three repeated samples for each group were used for callpeak by MACS2 (Version 2.1.2)  with the parameters qvalue < 0.05. Reads distributions across genes were presented using deeptools (Version 3.4.3), and genes were represented as lines sorted in descending order based on the signal strength in the heatmap. Motif analysis was performed using the HOMER's findMotifsGenome.pl tool (v4.11). Differential accessible peak analysis was performed by using bedtools software  and DESeq2 (Version 1.16.0), log2 fold change ≥ 1 or ≤ -1, and p-value < 0.05 were considered as the cutoff values. The raw data and processed data were uploaded into the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) with an accession number GSE233445.
Gene Ontology (GO) and Kyoto gene and genome encyclopedia (KEGG) signal pathway analysis of differential peaks nearest genes were performed using ClusterProfiler with a False Discovery Rate (FDR)-adjusted p-value cutoff of 0.05. Thus the significant GO categories and pathways were identified.
Integration Analysis of ATAC-seq and RNA-seq
RNA-seq analysis was as previously reported . The associated genes in open chromatin regions with enhanced and attenuated ATAC-seq signal overlapped with the up and down differentially expression genes (DEGs) in the transcriptomes, respectively. Further, GO and KEGG pathways analyses were performed as described above.
Proteomic and phosphoproteomic analysis
Cells were washed three times with cold PBS and lysed in lysis buffer containing 50 μl 50 mM Tris–HCl (pH 8.5), 57 μl ddH2O, 10 μl 400 mM 2-chloroacetamide (CAA) (Sigma, 22,790), 2 μl tris (2-carboxyethyl) phosphine (Sigma, 646547) and 1 μl phosphatase inhibitor cocktail. Samples were then boiled in a 95 °C water bath for 10 min in darkness. Subsequently, cooled at room temperature and diluted five-fold. The protein concentration was measured by using BCA protein assay, and the rest part of protein samples was digested with Lys-C (Wako) at 37 °C for 3 h. Protein samples underwent trypsin digestion (1:50) at 37 °C overnight. Afterward, samples were acidified by 1% trifluoroacetic acid (TFA) and ethyl acetate. The mixture of lysates was centrifuged at 15000g for 3 min, the bottom aqueous phase was collected followed by vacuum-dried centrifugation and desalting. 95% and 5% portions of each sample were used for phosphoproteomics and proteomics experiments, respectively.
Peptide samples were analyzed on the nanoElute coupled online with a timsTOF Pro mass spectrometer (Bruker). Peptides were re-dissolved in 0.1% formic acid (FA) and loaded onto a 25 cm in-house column n (360 μm OD × 75 μm inner diameter) packed with C18 resin (particle size 2.2 μm, pore size 100 Å, Michrom Bioresources). All peptides were separated onto the analytical column with a 120 min gradient (buffer A: 0.1% FA in ultrapure water; buffer B: 0.1% FA in acetonitrile) at a constant flow rate of 300 nL/min (90 min, 0 to 37% buffer B using a linear AB gradient of 2 to 22% of buffer B; 10 min, 22 to 37% of buffer B, 10 min, 37 to 80% of buffer B; 10 min, 80% of buffer B). Mass spectrometry was set under a data-dependent acquisition mode. Mass range was 100–1700 m/z at a resolution of 40,000.
MS raw files produced by LC–MS/MS were searched against the Uniprot human proteome database using PEAKS Studio X + software (Bioinformatics Solutions Inc.). Mass tolerances were 15 ppm for initial precursor mass and 6 ppm for final tolerance. Carbamidomethyl of cysteine (+ 57.0214 Da) was considered a fixed modification. Variable modifications were acetylation (+ 42.011 Da) at the N terminus of proteins, oxidation (+ 15.9949 Da) on methionine residues, and phosphorylation (+ 79.996 Da) on serine, threonine, or tyrosine residues. Up to 2 missed cleavages were allowed. The cutoff of FDR was 1% for all proteins, peptides, and phosphosites. Differential expression analyses were determined based on p-value < 0.05 and |log2(Fold change)|> 1.
Cell migration assays.
To measure the relative migration ability of cells, scratch-wound healing assays were performed. Cells were plated in a 6-well plate at a density of 50*104/well. After cells reached confluency, manually scraped the monolayer cell with a 10 μl pipette tip to create the wounds. The cells were then washed with phosphate saline buffer (PBS), replenished with serum-free α-MEM, photographed with the inverted phase-contrast microscope, and recorded the exact location of each wound. Afterward, the cells were placed back in the cell-culture incubator, and the recorded scratch regions were photographed 24 h after wounding. For statistical analysis, using Image J to outline the pictures, measure the area of the scratch, analysis and plot were performed using GraphPad Prism 8.
Cells were seeded in a 24-well plate at approximately 100% confluence and subsequently scraped with the pipette tip to create wounds. For scratch-wound healing experiments, cells were treated as described above. At 24 h after wounding, cells were fixed by 4% formaldehyde solution for 10 min, and washed three times with PBS, each time for 5 min. Cells were then permeabilized with 0.1% Triton X-100 for 5 min, and washed three times with PBS. To reduce nonspecific background staining, cells were blocked in 1% bovine serum albumin (BSA) for 30 min at room temperature. The rhodamine phalloidin staining solution (MESGEN, MF8204) was diluted in 1% BSA and incubated cells for 20 min at room temperature. The cells were rinsed by PBS three times, each time for 5 min. 200 μl per well of Hoechst 33,258 (Beyotime, C1018) was used to label the nuclear of cells. Cells were observed on fluorescence microscope (Nikon), and photographed at least three different representative regions at 4× magnification or 10× magnification.
Quantitative real-time PCR analysis
RNA was extracted from cells by using RNAiso Plus reagent (Takara, 9108) following the manufacturer’s protocols. RNA quality and concentration were determined by using a Nanovue spectrophotometer (GE Healthcare Life Sciences, Marlborough, MA, USA). cDNA was prepared by RNA reverse transcripted with PrimeScript RT reagent kit (Takara, RR047A). Subsequently, SYBR TB Green Premix Ex Taq™ kit (Takara, RR820A) was used for one-step real-time RT-PCR analysis on the MxPro-Mx3000P real-time PCR System. The experiment was repeated three times.
Normalize the expression value of the target gene in a given sample to the corresponding GAPDH expression. The relative expression value of the targeted genes was calculated by the 2 − △ △Ct method. The primers were: RhoA-F, 5’- GAAACTGGTGATTGTTGGTGATG-3’, RhoA-R, 5’- ACCGTGGGCACATAGACCT-3’, Rac1-F, 5’- ACGGAGCTGTTGGTAAAACCT-3’, Rac1-R, 5’- AGACGGTGGGGATGTACTCTC-3’, Cdc42-F, 5’- CCCATCGGAATATGTACCAACTG-3’, Cdc42-R, 5’- CCAAGAGTGTATGGCTCTCCAC-3’, Tcf4-F, 5’-GATGGGACTCCCTATGACCAC-3’, Tcf4-R, 5’- GAAAGGGTTCCTGGATTGCCC-3’, β-catenin-F, 5’- ATGGAGCCGGACAGAAAAGC-3’, β-catenin-R, 5’- TGGGAGGTGTCAACATCTTCTT-3’, Calpain1-F, 5’- ATGACAGAGGAGTTAATCACCCC-3’, Calpain1-R, 5’- GCCCGAAGCGTTTCATAATCC-3’, Calpain2-F, 5’- GGTCGCATGAGAGAGCCATC-3’, Calpain2-R, 5’- ATGCCCCGAGTTTTGCTGG-3’, Calpastatin-F, 5’- GGAAGGACAAACCAGAGAAGC-3’, Calpastatin -R, 5’- AGGGGCAGCTATCCAAATCTT-3’, GSK-3β-F, 5’- TTGGACAAAGGTCTTCCGGC-3’, GSK-3β-R, 5’- AAGAGTGCAGGTGTGTCTCG-3’, GAPDH-F, 5’- AACTCCCACTCTTCCACCTTC-3’, GAPDH-R, 5’- CCTGTTGCTGTAGCCGTATTC-3’.
Western blot and antibodies
Cells were collected and lysed in radio immunoprecipitation assay (RIPA) lysis buffer (Beyotime, P0013B) containing phenylmethanesulfonyl fluoride (PMSF) (Beyotime, ST506) and phosphatase inhibitor cocktail (Sigma, P5726). Total proteins were extracted in the supernatant. To prepare the nuclear and cytoplasmic extracts, the nuclear and cytoplasmic protein extraction kit (Beyotime, P0028) was applied to lyse the cells following the manufacturer’s instructions. After vigorous vortexing, the lysates were centrifuged at 12000g for 5 min at 4 °C. The supernatant was collected as the cytoplasmic protein. Add the nuclear protein extraction reagent to the pellet, vortexing and ice bathing alternately for 30 min, the lysates were centrifuged at 12000g for 10 min at 4 °C, and the supernatant was collected as the nuclear fraction. All extracted proteins were quantified using the Enhanced BCA Protein Assay Kit (Beyotime, P0010). The protein concentration between different groups was adjusted to be consistent before boiling.
For western blot assay, 40 µg of protein lysates were separated by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) (Beyotime, P0012AC) and transferred to polyvinylidene difluoride (PVDF) membranes (Biosharp, BS-PVDF-45/BS-PVDF-22). After blocking in 5% non-fat milk (Biosharp, BS102), PVDF membranes were then incubated overnight at 4 °C with primary antibodies and probed for Calpastatin (Cell Signaling Technology, #4146), Calpain 1 (Proteintech, 10538-1-AP), Calpain 2 (Proteintech, 11472-1-AP), β-catenin (Proteintech, 51067-2-AP), Tcf4 (Proteintech, 22337-1-AP), GSK3-β (Proteintech, 22104-1-AP), Lamin B (Wanleibio, WL01775), Rho A (Wanleibio, WL02853), Cdc42 (Proteintech, 10155-1-AP), β-Tubulin (Abmart, M20005), and GAPDH (Proteintech, 10494-1-AP). Peroxidase-conjugated IgG (Biosharp, BL003A) was used as the secondary antibody, and blots were developed using the enhanced chemiluminescence kit (Meilunbio, MA0186) detection.
Cells were seeded in a 24-well plate at approximately 60% confluence. Cells were fixed by 4% paraformaldehyde (PFA) for 30 min at room temperature, and washed three times with PBS, each time for 5 min. Cells were then blocked in immunol staining blocking buffer (Beyotime, P0102) for 1 h at room temperature. The primary antibodies were diluted in primary antibody dilution buffer for Immunol staining (Beyotime, P0262) and incubated the cells for overnight at 4 °C. PBS rinsed the cells five times, each time for 5 min. Secondary antibodies (Cell Signaling Technology, #8890S) were diluted in antibody dilution buffer (1:400) and incubated for 2 h at room temperature. Hoechst 33258 was used to tag the nuclear of cells. Cells were observed on fluorescence microscope, and photographed at least three different representative regions at 4× magnification or 10× magnification.
To determine calpain activity, casein zymography analysis was performed. The main principle is that the enzymes in the lysates are separated by the polyacrylamide gel containing casein. Casein molecules in zymography gels act as substrates for calpain in solutions containing calcium ions, and their degradation reflects enzymatic activity [35, 36]. In resolving gel, 0.21 mg casein was copolymerized with 4.9 ml ddH2O, 2.5 ml 1.5 mol/L Tris–HCl (pH 8.8), 2.5 ml polyacrylamide solution (4:0.16), 40 μl ammonium persulfate solution (APS), and 28 μl N,N,N’,N’-Tetramethylethylenediamine (TEMED). For stacking gel, it contained 6.5 ml ddH2O, 2.5 ml 1.32 mol/L Tris–HCl (pH 6.8), 1 ml polyacrylamide solution, 50 μl APS, and 10 μl TEMED. The casein gel was pre-run at 160 V for 30 min at 4 °C with electrode running buffer. The running buffer (10x) contained 25 mM Tris–HCl (pH 8.3), 192 mM glycine, 1 mM ethylene glycol-bis (β-aminoethyl ether)-N’,N’,N’,N’-tetraacetic acid (EGTA), 0.05% (v/v) 2-mercaptoethanol (2-MCE), and 1 mM ethylenediaminetetraacetic acid (EDTA). Subsequently, lysates were loaded with sample buffer (1:3, v/v), which included 150 mM Tris–HCl (pH 6.8), 0.04% (w/v) bromophenol blue, 20% (v/v) glycerol, and 0.75% (v/v) 2-MCE. The gel with samples loaded was run for electrophoresis at 125 V for 3 h at 4 °C. At the end of the electrophoresis, the gel was incubated in proteolysis buffer for 24 h at 20 °C. The proteolysis buffer (10x) was comprised of 20 mM Tris–HCl (pH 7.4), 10 mM dithiothreitol (DTT), and 4 mM CaCl2. The gel was stained with Coomassie Blue Fast Staining Solution (Beyotime, P0017) for 60 min, and the thermal cycler imaging system was applied to photograph the bands of the gel.
All mouse models are on a C57BL6/J background and breading in a SPF grade facility. The generation of conditional knockout mice was as follows: (1) Fam20cflox/flox mice (Texas A&M University College of Dentistry, America) mated with Osx-Cre mice (Biocytogen, China) from the progeny of which the Fam20cflox/+; Osx-Cre mice were selected. (2) By mating these mice with Fam20cflox/flox mice, the Osx-Cre; Fam20cflox/flox mice were obtained from the progeny, whose gene Fam20c was deleted in the pre-osteoblast. In this research, Osx-Cre; Fam20cflox/flox (conditional knockout, cKO) mice were used as the experimental group, and their littermates with Fam20cflox/flox genotype were put into the control group. The genotype identification was performed on the tails of postnatal mice by polymerase chain reaction, the genotyping of each mouse was repeated three times.
The bone tissue samples of 4-week-old mice were collected and fixed for 48 h, then placed in the decalcified fluid (15% EDTA). After the end of decalcification, the mouse tissues were placed in the automatic tissue dehydrator. The wax-soaked tissues were placed into the embedding machine. Conventional sections of tissue wax blocks were cut with a thickness of 3 µm, and the position and depth of the tissue sections in the cKO group and the control group were basically the same as far as possible. Immunohistochemical staining was performed on 3 µm thickness paraffin sections, the sections were sequentially deparaffinized, washed with pure water, and soaked in 3% H2O2 to block endogenous peroxidases. Antigen retrieval was performed using citrate buffer under high pressure, and then the slides were incubated in 10% normal goat serum and 5% bovine serum albumin (Solarbio, China) for 60 min, and incubated with the antibody for 2 h at room temperature. The sections were washed in PBST (PBS plus 0.1% Tween) and incubated with the secondary antibody at room temperature for 60 min. Then the slides were washed in PBST, and immunopositive reactions were visualized using a 3,3’-diaminobenzidine tetrahydrochloride solution. The nuclei were counterstained with hematoxylin. The slides were washed in distilled water and dehydrated in graded alcohol, cleared in xylene. Images were taken with a biological microscope (Nikon, Japan).
Each experiment was repeated at least three times. Statistical analysis was performed with GraphPad Prism8.0 software. Student’s t test was used for univariate comparison between two groups. P < 0.05 was considered to be different and statistically significant and presented as * P < 0.05, ** P < 0.01, or *** P < 0.001.
Mapping the ATAC-seq landscape of Fam20c-deficient osteoblasts
We first sought to map the chromatin landscape underlying the observed FAM20C-related transcriptional programs in the osteoblasts, and thus ATAC-seq was employed. This method for capturing chromatin accessibility regions by Tn5 transposases, which will integrate adapters payload into accessible chromatin regions such as enhancers, promoters, and insulators . We isolated nuclei from OB Fam20cf/f and OB Fam20cKO cells. Samples were massively sequenced in parallel, and the high-confidence regions of chromatin accessible regions were identified via overlapping peaks of the replicate samples. For data quality control, we examined whether our data between replicates were of good sequencing quality (Additional file 8: Table S1) and alignment of reads on the reference genome sequence (Additional file 9: Table S2). We then evaluated the length distribution of insert size, the chromatin was fragmented into nucleosome free, mono nucleosome, and di-nucleosome patterns (Additional file 1: Figure S1A). This analysis indicated that all samples present similar distribution of fragment sizes and there was a similar degree for Tn5 transposase access to chromatin, also revealing all sample reads were mainly enriched near the transcriptional start site (TSS) (Additional file 1: Figure S1B, C). Furthermore, examination of the Spearman correlation between samples based on the signals on merged ATAC-seq peaks, showed that the three samples were very close to each other in each group (Additional file 1: Figure S1D, E).
Enriched analysis of 9,028 filtered differential accessible peaks created a detailed global overview map between OB Fam20cf/f and OB Fam20cKO (Fig. 1A). These peaks were screened based on log2 (Fold Change) > 1 or < -1 and P value < 0.05, allowing identification the differentially accessible regions (DARs) between the two groups. Specifically, as shown in Fig. 1B, 5527 (62%) DARs were significantly intensified, and 3501 (38%) DARs were attenuated in Fam20c-deficient osteoblasts (Fig. 1B). Having demonstrated that peak annotation could link chromatin accessibility to gene regulation, and classify genome-wide functional regions into promoter, 5’ untranslated region (UTR), 3’ UTR, exon, intron and intergenic . We thus proceeded to test the distribution of DARs on the above genomic features. The results indicated that the genomic distribution of differentially peaks in Fam20c-deficient osteoblasts showed a slight shift toward intergenic region and a slight decrease in promoter compared with OB Fam20cf/f (Fig. 1C). Sequence motif analysis using HOMER revealed a prominent enrichment in OB Fam20cKO for the basic helix-loop-helix (bHLH) transcription factor family, including transcription 21 (TCF21) and atonal bHLH transcription factor 1 (ATOH1), and the basic leucine zipper (bZIP) transcription factor family including fos-like antigen 2 (FRA2), FRA1, junB proto-oncogene (JUNB), activating transcription factor 3 (ATF3), basic leucine zipper ATF-like transcription factor (BATF), TCF4, and others (Fig. 1E and Table 1). By contrast, some members of bZIP and zinc finger (ZC) family were highly downregulated in Fam20c-deficient osteoblasts, including FRA1, ATF3, BATF, activator protein-1 (AP-1), FRA2, JUNB, fos-like antigen 2 (FOSL2), JUN-AP1, BTB and CNC homology 2 (BACH2), and CCCTC-binding factor (CTCF) (Fig. 1F and Table 2).
Joint profiling of chromatin accessibility and gene expression
To create a chromatin accessibility and mRNA expression co-profiling analysis that is responsible for transcription factors regulation of downstream genes, we combined our ATAC-seq data with RNA-seq data. RNA-seq data was performed in our previous publication . Detailed gene expression and Q value were summarized in Additional file 10: Table S3. The filter criteria for differentially expressed genes (DEGs) were as follows: log2 (Fold Change) > 1 or < −1 and Q value < 0.05. We first intersected related genes of DARs in ATAC-seq with DEGs of RNA-seq to obtain overlapping genes. Overall, OB Fam20cKO versus OB Fam20cf/f identified 3875 up-regulated genes nearest to DARs, and 127 overlapped with up-regulated DEGs (Fig. 2A and Additional file 11: Table S4). While we obtained 2174 down-regulated genes associated with DARs, 321 down-regulated DEGs, and the overlap part was 109 genes (Fig. 2B and Additional file 12: Table S5). We further assessed gene functions that may be regulated by transcription factors by performing GO and KEGG pathway analysis on the corresponding overlapping gene between ATAC-seq and RNA-seq. GO categories associated with Fam20c-deficient genes included synapse organization, membrane depolarization during action potential, integrin-mediated signaling pathway, regulation of synapse structure or activity, and regulation of cell shape, linking the observed chromatin alternated events to cellular shape and synapse changes (Fig. 2C and Additional file 13: Table S6). KEGG signal pathway analysis revealed enrichment for Wnt signaling pathway, glycosaminoglycan biosynthesis, and extracellular matrix (ECM)-receptor interaction (Fig. 2D and Additional file 13: Table S7).
Proteomics characterization of Fam20c-deficient osteoblasts
Next, we extended our analyses to identify changes that span transcription factors to RNA to protein in Fam20c-deficient osteoblasts. Proteomic data can better reflect gene function than transcriptomic data . Using label-free quantitative proteomics, we quantified 8085 and 8146 peptides representing 4080 and 3874 proteins from OB Fam20cf/f and OB Fam20cKO with a false discovery rate (FDR) less than 1% at both peptide (Additional file 14: Table S8). We performed a heatmap analysis based on proteins differentially expressed, together, and identified 61 upregulated proteins and 51 downregulated proteins in the OB Fam20cKO group (p-value < 0.05, log2Fold change > 1). The upregulated and downregulated top 20 proteins were shown in Fig. 3A. For the phosphoproteomic data, a total of 8453 and 8503 phosphopeptides representing 4504 and 4444 phosphoproteins were identified from OB Fam20cf/f and OB Fam20cKO (Additional file 15: Table S9). Difference analysis-based phosphoproteomic data revealed 64 upregulated and 80 downregulated phosphoproteins in OB Fam20cKO, with normalization by the z-score (p-value < 0.05, log2Fold change > 1). Figure 3B showed the top 20 phosphorylated proteins up- and down-regulated.
Our following step was to learn more about the signaling pathways that these proteins and phosphoproteins regulated the biological behaviors of OB Fam20cKO. GSEA on the expression profiles of proteomics and phosphoproteomics was performed. The proteins associated with lamellipodium assembly and apoptotic mitochondrial were strongly enriched in OB Fam20cKO (Fig. 3C). OB Fam20cf/f showed high amounts enriched in cell growth, division, spreading, and protein catabolic process (Fig. 3C). In addition, for phosphorylated proteins, depletion of Fam20c resulted in increasing enriched in epithelial to mesenchymal transition (EMT) (Fig. 3D), which was consistent with our previous bioinformatic analysis of Fam20C in pan-cancer . Likewise, the phosphorylated proteins of OB Fam20cf/f showed high amounts enriched in cell division, cell cycle, cysteine-type endopeptidase activity, and protein kinase activity (Fig. 3D). Altogether, the deficiency of cell proliferation, migration, and protein catabolic was possibly a feature of OB Fam20cKO cell.
We attempted to identify key molecules from the proteome and phosphoproteome subgroups by combining the multi-omics data. Since there was a similar GSEA trend between proteomics and phosphorylated proteomics, we intersected the differential proteins in this two proteomics. As shown in Fig. 3E, only one protein expressed level and phosphorylated level were upregulated, while the expression and phosphorylated levels of 8 proteins were down-regulated in OB Fam20cKO. The volcano plots visualized these 9 intersected proteins (Fig. 3F, G). Noteworthy, these 9 intersected proteins contained Calpastatin (its protein entry was ICAL). This protein is an endogenous inhibitor of Calpain , and its phosphorylation state mediated the enzyme activity of Calpain has been reported . Calpain is a calcium-dependent cysteine protease, and, importantly, it plays a key role in regulating cellular physiological processes, including proteolysis, cytoskeleton remodeling, cell cycle regulation, apoptosis, cell signal transduction, and others [43, 44]. Moreover, the biological processes regulated by Calpain was highly related to the results of the above GSEA.
Furthermore, the activation of Calpain could promote proteolysis, thus realizing negative regulation of the Wnt pathway , which was highly enriched in joint profiling between ATAC-seq and RNA-seq (Fig. 2D). We next examined Calpastatin expression and phosphorylation levels, as well as Calpain levels.
Fam20c knockout influences Calpastatin/Calpain proteolysis system in vitro/vivo
Calpastatin/Calpain proteolysis system emerged in our research as a potential major protein, thus, we examined and verified the key molecular expression and activity of proteolysis system. Among the known 15 kinds of Calpain, Calpain 1 and Calpain 2 are the most widely distributed and studied. qPCR, Western blot, and Immunofluorescence analyses showed that Fam20c knockout resulted in a decreased expression of Calpastatin in OB Fam20cKO cells (Fig. 4A, B, E), whereas accompanied by a decreased phosphorylated expression level of Calpastatin (Fig. 4C). These were consistent with the results of proteomics and phosphoproteomics. Furthermore, Fam20c-deficient osteoblasts displayed a diminished expression of Calpain 1 and Calpain 2 (Fig. 4A, B, E) and attenuation of Calpain activity (Fig. 4D).
To validate our in vitro results, we used a conditional Fam20c knockout model (Osx-Cre; Fam20cflox/flox, abbreviated as “cKO mice” in this report) and matching control (Fam20cflox/flox, referred to as “Ctrl”). We isolated femurs of 4-week-old mice, removed bone marrow, and exacted mRNA from bone tissue. Analysis of the mRNA expression of Calpastatin in Ctrl versus cKO mice identified no significant change, whereas consistent with in vitro results, the gene expression of Calpain 1 and Calpain 2 decreased significantly (Fig. 4F). An immunohistochemistry (IHC) staining analysis displayed that Calpastatin immunoreactivity was low in preosteoblasts on trabecular bone surfaces and cartilage calcification area in cKO mice (Fig. 4G). No significant difference was found in the expression of cartilage proliferation area in cKO mice compared with Ctrl mice. Specifically, although there was less expression of Calpastatin in the hypertrophic cartilage area compared with the proliferation area in Ctrl mice, there was almost no expression in cKO mice. For the expression of Calpain 1 and Calpain 2, as anticipated, IHC staining showed a significant reduction in their immunoreactivity (Fig. 4H). Calpain 1 was mainly expressed in the articular surface cartilage proliferation layer (arrows A1 and B1), while the intensity of Calpain 2 expression was mainly in the cartilage proliferation layer (arrows A2 and B2). Hence, these data demonstrate that Calpain 1 and Calpain 2 might play a role in different locations.
Fam20c knockout-related changes in migration velocity correlate with changes in the organization of the F-actin cytoskeleton
To determine if the weakened Calpastatin/Calpain proteolysis system seen in Fam20c-deficient cells and cKO mice would cause dysfunction. We measured cell proliferation and migration, based on the results of GSEA analysis of proteomics and phosphoproteomics and the biological process regulated by Calpastatin/Calpain proteolysis system as noted in previous studies (Fig. 3C, D) [43, 44]. OB Fam20cKO showed less potent proliferation than control cells (Additional file 2: Figure S2). To dissect the role of FAM20C in osteoblast migration, we analyzed migration ability of OB Fam20cf/f and OB Fam20cKO by applying a monolayer scratch-wound assay in vitro. Fam20c-deficient cells displayed a relatively reduced cell migration (Fig. 5A). To confirm this, we quantified the scratch area at 0 h and 24 h applying Image J (MRI wound healing tool plugin), respectively, and through the formusla: (0 h scratch area -24 h scratch area)/ 0 h scratch area, to obtain migration rate. The results suggested that Fam20c-deficient cells had a lower migration rate, which was around 3 times slower than the OB Fam20cf/f (Fig. 5B).
Coordinated actin morphology and cytoskeleton arrangement are key steps in migration [46, 47]. In this regard, Calpain proteases also participate mediate cytoskeleton events [48, 49]. To this end, we assessed the indicators of actin cytoskeleton arrangement in OB Fam20cf/f and OB Fam20cKO. Using phalloidin staining for F-actin in a scratch wound assay, we observed the changes in the cytoskeleton of cells moving towards the site of wounding. Similarly, the loss of Fam20c was less adept in wound closure in cell migration (Fig. 5C). In addition, the differences in the F-actin filaments at the leading edge of the wound between OB Fam20cf/f and OB Fam20cKO are significant (Fig. 5C). At the leading edge of the scratch wound, OB Fam20cf/f cells exhibited a highly elongated shape through the extension of the pseudopod-like structure, while OB Fam20cKO cells adopt a square round in a relatively quiescent state. Furthermore, we found that the deletion of Fam20c contributed to decreased F-actin formation as indicated by the arrows in Fig. 5C, which was not observed in OB Fam20cf/f. These findings suggested that Fam20c may participate in the regulation of appropriate cellular actin morphology. Expression of the actin cytoskeleton organization regulators Ras homolog gene family member A (RhoA), Ras-related C3 botulinum toxin substrate 1 (Rac1), and Cell division cycle 42 (Cdc42) as measured by PCR and Western blot were decreased in OB Fam20cKO compared with OB Fam20cf/f (Fig. 5D, E).
Impaired Calpastatin/Calpain proteolysis system in Fam20c-deficient cells is associated with weak mobility and F-actin formation disorder signature
To further link the weak mobility phenotype in OB Fam20cKO to the impaired Calpastatin/Calpain proteolysis system, we treated OB Fam20cKO with low doses of Calcium chloride (CaCl2) (600 μmol/L) and repeated the above experiments. The concentration of CaCl2 was determined by the CCK8 experiment, the detailed information was seen in Additional file 3: Figure S3. In comparison with OB Fam20cf/f, OB Fam20cKO treated with CaCl2 showed significantly increased expression of Calpastatin, Calpain 1, and Calpain 2 (Fig. 6A, B). Of note, the gene expression of Calpain 1 and Calpain 2 exhibited a concentration-dependent elevation with the addition of CaCl2 in OB Fam20cKO (Fig. 6C). No significant phosphorylated expression level of Calpastatin difference was observed in OB Fam20cKO with or without CaCl2 (Fig. 6D), whereas OB Fam20cKO with CaCl2 showed higher Calpain activity (Fig. 6E). This treatment rescued almost entirely the deficit of OB Fam20cKO in migrating (Fig. 6F). In view of the results presented above, we further used PCR and Western blot assays to detect the expression of Rac1, RhoA, and Cdc42. Corroborating our predictions, there was an elevated expression in OB Fam20cKO treated with or without CaCl2 (Fig. 6G, H).
Changes of Wnt signal pathway in Fam20c-deficient cells
The chromatin accessibility signature of OB Fam20cKO suggested an enrichment of TCF4 motif (Fig. 1E), and Fam20c-deficient associated genes enrichment of Wnt signaling pathway was specifically observed in the joint profiling of ATAC-seq and RNA-seq data (Fig. 2D). Additionally, the depletion of Fam20c significantly changed Calpastatin/Calpain proteolysis system, as indicated by proteomics and phosphorylated proteomics analysis (Fig. 4). Importantly, previous studies have implicated Wnt signaling pathway could be regulated by Calpain activation to promote β-catenin degradation [45, 50]. We, therefore, set out to verify the expression of key molecules of Wnt signaling pathway in the subsequent experiments. RT-PCR revealed that β-catenin and Tcf4 mRNA becomes prominently expressed in OB Fam20cKO, and glycogen synthase kinase-3 (GSK-3β) was expressed at low levels (Fig. 7A). Consistent with the PCR results, protein levels displayed a similar trend (Fig. 7B). To confirm whether these changes mediated canonical Wnt signaling regulates nuclear translocation, we observed in OB Fam20cKO that Tcf4 and β-catenin were mainly localized in the nuclear, and a decreased expression of GSK-3β was detected (Fig. 7C–F). Among the Wnt family, RNA-seq analysis revealed the mRNA level of canonical Wnt ligands, Wnt family member 7a (Wnt7a), was suppressed (Fig. 7G). And other canonical Wnt ligands were not detected in the RNA-seq data. These data suggest that Fam20c activates Wnt/β-catenin signaling through β-catenin and Tcf4 transcription and translation. To extend this observation in vivo, we analyzed these protein expressions from Ctrl and cKO mice. During bone formation, GSK-3β was strongly expressed in the area where cartilage calcified to form bone trabecular in Ctrl mice, and a decrease was observed in the cKO mice (Fig. 7H). Furthermore, compared with Ctrl mice, Tcf4 was highly accumulated in the nucleus of chondrocytes with an obvious elevated expression level in cKO mice (Fig. 7I). The role of β-catenin was also measured in vivo, ~ 90% of the nuclei displayed increased expression of Tcf4 in the proliferating layer of articular cartilage of the femur and hypertrophic cartilaginous layer of the articular surface (Fig. 7J).
We then investigated whether the Wnt signaling pathway is a target of Calpastatin/Calpain proteolysis system. After the treatment of CaCl2, we observed reduced Tcf4 and β-catenin and elevated GSK-3β in OB Fam20cKO, immunofluorescence further documented a decreased location of Tcf4 and β-catenin in OB Fam20cKO nuclei (Fig. 7K–N).
Increasing evidence has pointed to essential roles for FAM20C phosphorylated substrates within S-x-E/pS motifs in regulating many physiological processes [1, 7, 51]. Despite this intrigue, the knowledge about FAM20C has remained in the infant stage resulting from its complex regulatory mechanisms. The widespread multiple FAM20C substrates, ranging from the secreted proteins to intracellular proteins, of which phosphorylated levels are altered [52, 53]. Accordingly, the complex regulatory network by FAM20C, including epigenetics, gene expression, and protein interactions, is attracting broad attention. Comprehensive multi-omics analysis could broaden our knowledge of the molecular events relevant to biological processes [54, 55]. Here, we hypothesized that the joint profiling analysis of chromatin accessibility, gene expression, proteomics, and phosphorylated proteomics data between OB Fam20cf/f and OB Fam20cKO, may capture some new regulatory mechanism of FAM20C.
ATAC-seq allows the identification of genomic regions associated with gene-regulatory activity, thus providing a method to infer transcription factor (TF) activity . Our data demonstrated that lack of Fam20c generally intensified 12% of the accessible regions in osteoblasts, leading to a global activation of downstream TFs (Fig. 1B). Motif analysis illuminated that TCF4 was highly enriched in the chromatin opened of Fam20c deficient osteoblasts (Fig. 1E). TCF4, a member belonging to the bHLH TF family, was regarded to play a vital role in many developmental processes [56, 57]. Furthermore, when Wnt ligands bind to the Frizzled-Lrp5/6 receptor complex, β-catenin degradation is blocked and translocated to the nucleus, allowing TCF to bind and to activate Wnt target genes [58, 59]. Consistently, we also found a significant enrichment of the KEGG pathway in the Wnt signaling pathway through joint profiling ATAC-seq and RNA-seq (Fig. 2D). Previous studies have shown that tissue homeostasis and cell maintenance could be impacted by Wnt signaling pathway [60, 61]. In this research, PCR, Western blot, and Immunofluorescence both elucidated that Fam20c knockout significantly activated the transcriptional activity of β-catenin and Tcf4 signaling, as well as elevated the protein expression levels of them, especially, the increased expression levels in the nucleus. While GSK-3β displayed a contrary trend, that is, the total expression level and the nuclear expression level decreased (Fig. 7A–F).
Moreover, Qin’s research found that TCF1 and LEF1 were down-regulated in the vertebrae of Sox2-Cre; Fam20cfoxl/flox mice . Further, TCF4 is combined with β-catenin upstream of BMP. In our research, knockout of Fam20c in osteoblasts in vivo, changes were only found in the hypertrophic cartilage layer, indicating that there were differences in vivo and in vitro, which might be caused by the paracrine and juxtacrine functions of cells. In addition, cartilage thickening has been observed in several condition knockout animal models, whether this phenomenon was related to the Wnt signaling pathway remains to be further investigated. Furthermore, studies have shown that Wnt signal can inhibit the expression of activated Col2a1, and the overexpression of β-catenin in chondrocytes can cause serious congenital chondrosis, indicating that its overexpression affects the bone development process, supporting the results of our study . Therefore, we believe that after Osx-Cre knockout of Fam20c, the expression of downstream GSK-3β in the cytoplasm is inhibited, thus increasing the entry of β-catenin into the nucleus, and the overexpression of β-catenin in the nucleus overactivated TCF4. Finally, the proliferation of chondrocytes is induced through the Wnt signaling pathway, thus affecting the process of bone development.
In the current study, extensive Mass Spectrometry (MS)-based proteomics and phosphoproteomics analysis permitted in-depth screening of the changes in protein expression, post-translational phosphorylated modification, and biological processes alteration in Fam20c deficient osteoblasts. These included but were not limited to: (I) weakened changes in the cell growth, where the predominant changes were the regulation of mitotic cell cycle, cell division, and regulation of cell cycle G2 M phase transition; (II) a decrease in alteration in regulation of proteolysis involved in protein catabolic process; (III) alteration in cell migration (Fig. 3C, D). The deficiency of cell growth in OB Fam20cKO was proved in previous studies [18, 24]. Of note, in this research, we verified the major driver protein was Calpastatin, which was less expressed and decreased phosphorylated levels in OB Fam20cKO than in OB Fam20cf/f (Fig. 4). We also performed experiments to knock out FAM20C in human 293 T cells and obtained a similar trend (Additional file 6: Figure S6). Calpastatin belongs to the calpain system; mediates cell motility, cell cycle, signal pathways transduction, apoptosis, and regulates gene expression via proteolytic degradation [43, 64, 65]. The other two members of the calpain system, Calpain 1 and Calpain 2, activity was mediated by Calpastatin phosphorylation state [42, 66, 67]. Specifically, the phosphorylated Calpastatin presented relatively static, located near the cell nucleus, and presented an aggregation state, that is, the activity of Calpastatin was inhibited, thus enhancing the activity of Calpain. On the contrary, the dephosphorylation of Calpastatin made its distribution in cells relatively favorable, and the activity of Calpastatin increased, which then weakens the activity of Calpain. Consistently, we observed Fam20c-deficient osteoblasts displayed a diminished expression of Calpain 1 and Calpain 2, and attenuation of Calpain activity (Fig. 4D). OB Fam20cKO cells under the influence of impaired Calpastatin/Calpain proteolysis system exhibited weak mobility and guide differential remodeling of cytoskeletons at the leading edge versus the OB Fam20cf/f, with the inhibition of RhoA, Rac1, and Cdc42 (Fig. 5). The Rho family of small guanosine triphosphatases (GTPases), particularly RhoA, Rac1, and Cdc42, play essential roles in regulating the cell cycle and actin cytoskeleton, affecting cell adhesion and migration [68,69,70]. For cell migration, specifically, RhoA participated in the assembly of stress fibers, Rac1 was responsible for the formation of lamellar pseudopods, and Cdc42 regulated the shape of filamentous pseudopods [71, 72]. Moreover, we observed the weakened mobility was rescued by CaCl2 in OB Fam20cKO, accompanied by a concentration-dependent elevation of Calpain 1 and Calpain 2 expression and activity but no significant changes in Calpastatin phosphorylation level (Fig. 6), further indicating that cells migration barriers were associated with Calpastatin/Calpain proteolysis system. Calpain could regulate a variety of signaling pathways through proteolysis of target proteins, including cell cycle (cyclin D1 and cyclin E), cell survival (nuclear factor- κ B), and apoptosis (Bcl2 family and caspases) . It was also possible that Calpain mediated cell migration through multiple mechanisms, as noted in previous studies, including the cleavage of integrin-related complex, adhesive plaque, Focal adhesion kinase (FAK), Paxillin, Talin, and other components, resulting in reduced cell adhesion and increased mobility [74, 75]. Similar to our in vitro results, we observed the expression of Calpastatin, Calpain 1, and Calpain 2 in the growth plate region in animal models at the time when articular cartilage changes were most visible (i.e. 4 weeks after birth), and the results were consistent with cytological trends. In the knockout group, the expression of Calpastatin in the calcified area of growth plate cartilage, the pre-osteoblasts on the surface of bone trabeculae and hypertrophic cartilage decreased significantly. Calpain 1 and Calpain 2 were expressed in different positions in the cartilage, among which Calpain 1 was mainly expressed in the surface and middle layers of articular cartilage, and Calpain 2 was mainly expressed in the deep layer. The expressions of both decreased significantly in the knockout group. Real-time quantitative PCR detection showed that the expression levels of Calpain 1 and Calpain 2 genes in bone tissue decreased. These results further prove that the Calpastatin/Calpain proteolysis system is inhibited overall after the conditional knockout of Fam20c, and the dysfunction of this system can cause abnormalities in its downstream signaling pathway, inhibit the functions of pre-osteoblasts and hypertrophic chondrocytes, and ultimately affect the process of cartilage osteogenesis and cancellous bone formation.
Regarding Calpastatin/Calpain proteolysis system’s influence on Wnt signaling pathway, we took note of some previous studies showing more than 100 substrates were downstream of Calpain, including β-catenin and GSK-3β of Wnt signaling pathway . Briefly, Calpain activity could regulate β-catenin expression level, intracellular localization, and function .
The down-regulated of Calpain 2 expression, resulted in elevated β-catenin expression and accumulate in the nucleus. Conversely, Calpain activation could promote β-catenin degradation to achieve negative regulation of Wnt signaling pathway . Lade et al.  reported that Calpain 1 induced N-terminal truncation of β-catenin in the mouse liver development process, this might play a vital role in regulating the differentiation of hepatoblasts. In addition, Goni-Oliver and Ma et al. [78, 79] reported that Calpain truncated C- and N-terminal self-suppressing domains of GSK-3β to activate GSK-3β. In this regard, we added CaCl2 to OB Fam20cKO with weak calpain activity to activate calpain activity, and then the cell migration ability was enhanced, β-catenin and Tcf4 expression levels decreased, a decreased location of Tcf4 and β-catenin in OB Fam20cKO nuclei (Fig. 7K–N). Another key to note in our study is that previous studies have indicated that activation of Wnt signaling pathway could promote cell proliferation and migration [80, 81]. However, activating Wnt signaling pathway had less pronounced benefits, we do not detect accelerated cell migration in OB Fam20cKO. Given that we found Fam20c depletion to cause a weakened Calpastatin/Calpain proteolysis system, which could promote the hydrolysis of migration-related proteins to accelerate cell migration and negatively manage the Wnt signaling pathway. Thus, in the absence of Fam20c, this phenomenon could be explained by the possible explanation that Calpastatin/Calpain proteolysis system acts predominantly in regulating cell migration.
Additional analysis of the Wnt ligands, as shown here, only Wnt7a was observed significant changes in canonical Wnt ligands (Fig. 7G). The overlapped genes obtained from ATAC-seq and RNA-seq were projected into the KEGG pathway and visualized, which also showed the enrichment of Wnt7a (Additional file 4: Figure S4). Several studies have highlighted the importance of Wnt7a in regulating cell migration and invasion [82, 83], however, its role was cell-specific. Xie et al. reported that Wnt7a promoted oral squamous cell carcinoma cell migration induced by epidermal growth factor via the activation of the β-catenin/Matrix metalloproteinase-9 (MMP-9) signal pathway . Nevertheless, Lan et al. proposed the overexpression of Wnt7a inhibited the growth and migration of hepatocellular carcinoma in β-catenin-independent manner . Taken together, our findings supported Fam20c depletion led to Wnt ligands down-regulated in osteoblasts, however, β-catenin increasingly translocated to the nucleus, triggering Wnt signaling pathway ligands-independent manner.
Furthermore, a series of elegant papers have reported that FAM20C could phosphorylate bone morphogenetic protein 4 (BMP4), while BMP2 and BMP7 might be possible substrates of FAM20C due to possessing S-x-E/pS motifs [7, 85]. Intriguingly, we have previously observed that in Fam20c salivary gland conditional knockout mice, the transcriptional and protein levels of BMP2 and BMP7 were up-regulated, while protein expression of BMP4 was reduced but did not change its transcript level . The distribution pattern of BMP4 expression in salivary glands was almost not observed in the extracellular matrix, but concentrated in the cytoplasm, indicating that knockout of Fam20c caused abnormal secretion of BMP4. Also, Liu et al. showed mRNA levels of BMP2 and BMP7 were similarly increased in OB Fam20cKO, however, their downstream remains in a repressed state, implying that phosphorylation may affect the activity of BMP ligands . In addition, as illustrated by our analyses, Fam20c depletion induced Calpastatin dephosphorylation, decreased Calpain activity, elevated β-catenin expression levels, and translocated into nucleus to combine with Tcf4, leading to the activation of Wnt signaling pathway in a ligands-independent manner. Since Wnt signaling pathway regulated the transcription of BMP [87, 88], thus, our findings provided a possibility to explain the enhanced transcription of BMP in Fam20c-deficient animal and cell models. On the other hand, increased transcription of BMP, but not activation of BMP signaling pathway, indicating that conduction inhibition occurred in BMP signaling pathway.
Finally, several prior studies have shown the loss of bone maturation phenotype in Fam20c cKO mice and loss-of-function of the in vivo gained osteoblasts [11, 15, 89]. Also, the absence of Fam20c distorted cell biological processes. Osteogenic differentiation blockade occurs in OB Fam20cKO, which could not be rescued by the addition of extracellular matrix proteins extracted from normal bone tissue, suggesting that FAM20C might regulate cell behaviors in a cell-autonomous manner . Another characteristic of OB Fam20cKO, as above elucidated, was their paucity of BMP signaling pathway activation even with the elevated BMP ligands transcription, indicating that there was conduction abnormality in BMP signal pathway. Our integrated analysis revealed alterations of Calpastatin/Calpain proteolysis system and Wnt signaling pathway among the most dramatic differences between OB Fam20cf/f and OB Fam20cKO. The aberrant Calpastatin/Calpain proteolysis system could alter the cleavage of specific proteins, regulating pathological and normal physiological processes that are thought to play critical roles in governing homeostasis [90, 91]. Also, the coordination of cell fates and tissue homeostasis is controlled by Wnt signaling pathway that regulates cell-to-cell communication [92, 93]. Then again, we found the signaling pathways related to cell signal transduction were significantly influenced in OB Fam20cKO, including mitogen-activated protein kinases (MAPK) signaling pathway, Ras signaling pathway, mammalian target of rapamycin (mTOR) signaling pathway, cAMP signaling pathway, etc.  (Additional file 5: Figure S5). These data provide experimental support for a paradigm in which OB Fam20cKO appears the phenomenon of cellular homeostasis imbalance. Homeostasis is the basis of maintaining all kinds of important functions of the body, which involves tissue, cell, signal pathway, molecule, and others, breaking cell homeostasis will cause various pathological processes [94, 95]. These results comprise that the homeostatic imbalance of FAM20C knockout osteoblasts may involve changes in multiple signaling pathways in the conduction system.
Taken together, our current work provides an integrated and comprehensive analysis of FAM20C using multiple omics, including ATAC-seq, RNA-seq, proteomics, and phosphoproteomics. We revealed that Calpastatin/Calpain proteolysis system and Wnt signaling pathway alterations were possibly the key factors between OB Fam20cf/f and OB Fam20cKO. Targeting Calpastatin phosphorylation may provide a promising research direction to study the role of FAM20C in development. However, this research still existed restrictions. Initially, the protein–protein interaction between FAM20C and Calpastatin needs to be examined in follow-up studies. In the future experiment, we will add FAM20C recombinant protein to OB Fam20cKO to observe the change of phosphorylation of Calpastatin. Applying the radiolabeled technique to test the phosphorylation site of Calpastatin, constructing a mutant vector of the Calpastatin site to transfect OB Fam20cKO to activate the phosphorylation of Calpastatin, and observing the expression and activity of Calpastatin/Calpain proteolytic system and its effect on osteoblast function. Second, we performed phosphorylated proteomics, and detected the difference in Calpastatin phosphorylation and the changes in the phosphorylated site after Fam20c gene knockout. Nevertheless, the exact phosphorylated site of Calpastatin has not been verified. Third, in the present study, we found the Calpastatin/Calpain proteolysis system could negatively regulate the Wnt signaling pathway. Further studies are needed to conduct experiments associated with Calpain activators/ inhibitors as well as activators/antagonists of the Wnt pathway to further define this regulatory relationship.
Availability of data and materials
Any data and R script in this study can be obtained from the corresponding author upon reasonable request. The final manuscript was read and approved by all authors.
Tagliabracci VS, Engel JL, Wen J, Wiley SE, Worby CA, Kinch LN, et al. Secreted kinase phosphorylates extracellular proteins that regulate biomineralization. Science. 2012;336:1150–3. https://doi.org/10.1126/science.1217817.
Wang X, Hao J, Xie Y, Sun Y, Hernandez B, Yamoah AK, et al. Expression of FAM20C in the osteogenesis and odontogenesis of mouse. J Histochem Cytochem. 2010;58:957–67. https://doi.org/10.1369/jhc.2010.956565.
Du EX, Wang XF, Yang WC, Kaback D, Yee SP, Qin CL, et al. Characterization of Fam20C expression in odontogenesis and osteogenesis using transgenic mice. Int J Oral Sci. 2015;7:89–94. https://doi.org/10.1038/ijos.2014.67.
Nalbant D, Youn H, Nalbant SI, Sharma S, Cobos E, Beale EG, et al. FAM20: an evolutionarily conserved family of secreted proteins expressed in hematopoietic cells. BMC Genomics. 2005;6:11. https://doi.org/10.1186/1471-2164-6-11.
Ishikawa HO, Xu A, Ogura E, Manning G, Irvine KD. The Raine syndrome protein FAM20C is a Golgi kinase that phosphorylates bio-mineralization proteins. PLoS One. 2012;7:e42988. https://doi.org/10.1371/journal.pone.0042988.
Tagliabracci VS, Engel JL, Wiley SE, Xiao J, Gonzalez DJ, Nidumanda Appaiah H, et al. Dynamic regulation of FGF23 by Fam20C phosphorylation, GalNAc-T3 glycosylation, and furin proteolysis. Proc Natl Acad Sci USA. 2014;111:5520–5. https://doi.org/10.1073/pnas.1402218111.
Tagliabracci VS, Wiley SE, Guo X, Kinch LN, Durrant E, Wen J, et al. A single kinase generates the majority of the secreted phosphoproteome. Cell. 2015;161:1619–32. https://doi.org/10.1016/j.cell.2015.05.028.
Pollak AJ, Haghighi K, Kunduri S, Arvanitis DA, Bidwell PA, Liu GS, et al. Phosphorylation of serine96 of histidine-rich calcium-binding protein by the Fam20C kinase functions to prevent cardiac arrhythmia. Proc Natl Acad Sci USA. 2017;114:9098–103. https://doi.org/10.1073/pnas.1706441114.
Simpson MA, Hsu R, Keir LS, Hao J, Sivapalan G, Ernst LM, et al. Mutations in FAM20C are associated with lethal osteosclerotic bone dysplasia (Raine syndrome), highlighting a crucial molecule in bone development. Am J Hum Genet. 2007;81:906–12. https://doi.org/10.1086/522240.
Kinoshita Y, Hori M, Taguchi M, Fukumoto S. Functional analysis of mutant FAM20C in Raine syndrome with FGF23-related hypophosphatemia. Bone. 2014;67:145–51. https://doi.org/10.1016/j.bone.2014.07.009.
Wang X, Wang S, Li C, Gao T, Liu Y, Rangiani A, et al. Inactivation of a novel FGF23 regulator, FAM20C, leads to hypophosphatemic rickets in mice. PLoS Genet. 2012;8:e1002708. https://doi.org/10.1371/journal.pgen.1002708.
Bhattacharyya N, Chong WH, Gafni RI, Collins MT. Fibroblast growth factor 23: state of the field and future directions. Trends Endocrinol Metab. 2012;23:610–8. https://doi.org/10.1016/j.tem.2012.07.002.
Li L, Liu P, Lv X, Yu T, Jin X, Wang R, et al. Ablation of FAM20C caused short root defects via suppressing the BMP signaling pathway in mice. J Orofac Orthop. 2022. https://doi.org/10.1007/s00056-022-00386-7.
Wang X, Wang J, Liu Y, Yuan B, Ruest LB, Feng JQ, et al. The specific role of FAM20C in dentinogenesis. J Dent Res. 2015;94:330–6. https://doi.org/10.1177/0022034514563334.
Vogel P, Hansen GM, Read RW, Vance RB, Thiel M, Liu J, et al. Amelogenesis imperfecta and other biomineralization defects in Fam20a and Fam20c null mice. Vet Pathol. 2012;49:998–1017. https://doi.org/10.1177/0300985812453177.
Hao J, Narayanan K, Muni T, Ramachandran A, George A. Dentin matrix protein 4, a novel secretory calcium-binding protein that modulates odontoblast differentiation. J Biol Chem. 2007;282:15357–65. https://doi.org/10.1074/jbc.M701547200.
Oya K, Ishida K, Nishida T, Sato S, Kishino M, Hirose K, et al. Immunohistochemical analysis of dentin matrix protein 1 (Dmp1) phosphorylation by Fam20C in bone: implications for the induction of biomineralization. Histochem Cell Biol. 2017;147:341–51. https://doi.org/10.1007/s00418-016-1490-z.
Liu C, Zhang H, Jani P, Wang X, Lu Y, Li N, et al. FAM20C regulates osteoblast behaviors and intracellular signaling pathways in a cell-autonomous manner. J Cell Physiol. 2018;233:3476–86. https://doi.org/10.1002/jcp.26200.
Liu C, Zhou N, Wang Y, Zhang H, Jani P, Wang X, et al. Abrogation of Fam20c altered cell behaviors and BMP signaling of immortalized dental mesenchymal cells. Exp Cell Res. 2018;363:188–95. https://doi.org/10.1016/j.yexcr.2018.01.004.
Chen D, Zhao X, Sui Z, Niu H, Chen L, Hu C, et al. A multi-omics investigation of the molecular characteristics and classification of six metabolic syndrome relevant diseases. Theranostics. 2020;10:2029–46. https://doi.org/10.7150/thno.41106.
Sun YV, Hu YJ. Integrative analysis of multi-omics data for discovery and functional studies of complex human diseases. Adv Genet. 2016;93:147–90. https://doi.org/10.1016/bs.adgen.2015.11.004.
Palma-Lara I, Perez-Ramirez M, Garcia Alonso-Themann P, Espinosa-Garcia AM, Godinez-Aguilar R, Bonilla-Delgado J, et al. FAM20C overview classic and novel targets, pathogenic variants and raine syndrome phenotypes. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22158039.
Zhang J, Zhu Q, Wang X, Yu J, Chen X, Wang J, et al. Secretory kinase Fam20C tunes endoplasmic reticulum redox state via phosphorylation of Ero1α. Embo J. 2018. https://doi.org/10.15252/embj.201798699.
Geng YW, Zhang Z, Jin H, Da JL, Zhang K, Wang JQ, et al. Mesenchymal-to-epithelial transition of osteoblasts induced by Fam20c knockout. Genes Genomics. 2022;44:155–64. https://doi.org/10.1007/s13258-021-01170-4.
Yan J, Risacher SL, Shen L, Saykin AJ. Network approaches to systems biology analysis of complex disease: integrative methods for multi-omics data. Brief Bioinform. 2018;19:1370–81. https://doi.org/10.1093/bib/bbx066.
Subramanian I, Verma S, Kumar S, Jere A, Anamika K. Multi-omics data integration, interpretation, and its application. Bioinform Biol Insights. 2020;14:1177932219899051. https://doi.org/10.1177/1177932219899051.
Pinu FR, Beale DJ, Paten AM, Kouremenos K, Swarup S, Schirra HJ, et al. Systems biology and multi-omics integration: viewpoints from the metabolomics research community. Metabolites. 2019. https://doi.org/10.3390/metabo9040076.
Geng YW, Zhang Z, Jin H, Da JL, Zhang K, Wang JQ, et al. Mesenchymal-to-epithelial transition of osteoblasts induced by Fam20c knockout. Genes Genomics. 2022. https://doi.org/10.1007/s13258-021-01170-4.
Buenrostro JD, Wu B, Chang HY, Greenleaf WJ. ATAC-seq: a method for assaying chromatin accessibility genome-wide. Curr Protoc Mol Biol. 2015. https://doi.org/10.1002/0471142727.mb2129s109.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20. https://doi.org/10.1093/bioinformatics/btu170.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137. https://doi.org/10.1186/gb-2008-9-9-r137.
Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. deepTools: a flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42:W187-191. https://doi.org/10.1093/nar/gku365.
Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2. https://doi.org/10.1093/bioinformatics/btq033.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. https://doi.org/10.1186/s13059-014-0550-8.
Biswas AK, Tandon S. Casein zymography for analysis of Calpain-1 and Calpain-2 activity. Methods Mol Biol. 2019;1915:31–8. https://doi.org/10.1007/978-1-4939-8988-1_3.
Biswas AK, Tandon S, Beura CK. Identification of different domains of calpain and calpastatin from chicken blood and their role in post-mortem aging of meat during holding at refrigeration temperatures. Food Chem. 2016;200:315–21. https://doi.org/10.1016/j.foodchem.2016.01.031.
Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods. 2013;10:1213–8. https://doi.org/10.1038/nmeth.2688.
Yan F, Powell DR, Curtis DJ, Wong NC. From reads to insight: a hitchhiker’s guide to ATAC-seq data analysis. Genome Biol. 2020;21:22. https://doi.org/10.1186/s13059-020-1929-3.
Wang J, Ma Z, Carr SA, Mertins P, Zhang H, Zhang Z, et al. Proteome profiling outperforms transcriptome profiling for coexpression based gene function prediction. Mol Cell Proteomics. 2017;16:121–34. https://doi.org/10.1074/mcp.M116.060301.
Liu X, Zhan Y, Xu W, Liu X, Geng Y, Liu L, et al. Prognostic and immunological role of Fam20C in pan-cancer. 2021. Biosci Rep. https://doi.org/10.1042/BSR20201920.
Flores-Toro J, Chun SK, Shin JK, Campbell J, Lichtenberger M, Chapman W, et al. Critical roles of Calpastatin in Ischemia/Reperfusion injury in aged livers. Cells. 2021. https://doi.org/10.3390/cells10081863.
Bassett EA, Palanichamy K, Pearson M, McElroy JP, Haque SJ, Bell EH, et al. Calpastatin phosphorylation regulates radiation-induced calpain activity in glioblastoma. Oncotarget. 2018;9:14597–607. https://doi.org/10.18632/oncotarget.24523.
Nian H, Ma B. Calpain-calpastatin system and cancer progression. Biol Rev Camb Philos Soc. 2021;96:961–75. https://doi.org/10.1111/brv.12686.
Chen Y, Su Z, Liu F. Effects of functionally diverse calpain system on immune cells. Immunol Res. 2021;69:8–17. https://doi.org/10.1007/s12026-021-09177-5.
Benetti R, Copetti T, Dell’Orso S, Melloni E, Brancolini C, Monte M, et al. The calpain system is involved in the constitutive regulation of beta-catenin signaling functions. J Biol Chem. 2005;280:22070–80. https://doi.org/10.1074/jbc.M501810200.
Ridley AJ, Schwartz MA, Burridge K, Firtel RA, Ginsberg MH, Borisy G, et al. Cell migration: integrating signals from front to back. Science. 2003;302:1704–9. https://doi.org/10.1126/science.1092053.
Chodniewicz D, Klemke RL. Guiding cell migration through directed extension and stabilization of pseudopodia. Exp Cell Res. 2004;301:31–7. https://doi.org/10.1016/j.yexcr.2004.08.006.
de Lourdes O-G, Bastian Y, Rivera-Carrera V, Garcia-Tiscareno IA, Zapata-Zuniga M, Lopez JA, et al. Calpain participates in cortical cytoskeleton modification and DNA release during neutrophil extracellular trap formation. Int Arch Allergy Immunol. 2021;182:877–87. https://doi.org/10.1159/000515201.
Fettucciari K, Quotadamo F, Noce R, Palumbo C, Modesti A, Rosati E, et al. Group B Streptococcus (GBS) disrupts by calpain activation the actin and microtubule cytoskeleton of macrophages. Cell Microbiol. 2011;13:859–84. https://doi.org/10.1111/j.1462-5822.2011.01584.x.
Lade A, Ranganathan S, Luo J, Monga SP. Calpain induces N-terminal truncation of β-catenin in normal murine liver development: diagnostic implications in hepatoblastomas. J Biol Chem. 2012;287:22789–98. https://doi.org/10.1074/jbc.M112.378224.
Wu Y, Wang H, Liu C. From biomineralization to tumorogenesis-the expanding insight of the physiological and pathological roles of Fam20C. 2021. Biosci Rep. https://doi.org/10.1042/BSR20210040.
Gerson-Gurwitz A, Worby CA, Lee KY, Khaliullin R, Bouffard J, Cheerambathur D, et al. Ancestral roles of the Fam20C family of secreted protein kinases revealed in C. elegans. J Cell Biol. 2019;218:3795–811. https://doi.org/10.1083/jcb.201807041.
Lin JH, Lin IP, Ohyama Y, Mochida H, Kudo A, Kaku M, et al. FAM20C directly binds to and phosphorylates Periostin. Sci Rep. 2020;10:17155. https://doi.org/10.1038/s41598-020-74400-6.
Guo L, Li S, Yan X, Shen L, Xia D, Xiong Y, et al. A comprehensive multi-omics analysis reveals molecular features associated with cancer via RNA cross-talks in the Notch signaling pathway. Comput Struct Biotechnol J. 2022;20:3972–85. https://doi.org/10.1016/j.csbj.2022.07.036.
da Silveira WA, Fazelinia H, Rosenthal SB, Laiakis EC, Kim MS, Meydan C, et al. Comprehensive multi-omics analysis reveals mitochondrial stress as a central biological hub for spaceflight impact. Cell. 2020;183:1185-1201.e1120. https://doi.org/10.1016/j.cell.2020.11.002.
Forrest MP, Hill MJ, Quantock AJ, Martin-Rendon E, Blake DJ. The emerging roles of TCF4 in disease and development. Trends Mol Med. 2014;20:322–31. https://doi.org/10.1016/j.molmed.2014.01.010.
Hennig KM, Fass DM, Zhao WN, Sheridan SD, Fu T, Erdin S, et al. WNT/β-catenin pathway and epigenetic mechanisms regulate the Pitt-Hopkins syndrome and schizophrenia risk gene TCF4. Mol Neuropsychiatry. 2017;3:53–71. https://doi.org/10.1159/000475666.
Behrens J, Jerchow BA, Wurtele M, Grimm J, Asbrand C, Wirtz R, et al. Functional interaction of an axin homolog, conductin, with beta-catenin, APC, and GSK3beta. Science. 1998;280:596–9. https://doi.org/10.1126/science.280.5363.596.
MacDonald BT, Tamai K, He X. Wnt/beta-catenin signaling: components, mechanisms, and diseases. Dev Cell. 2009;17:9–26. https://doi.org/10.1016/j.devcel.2009.06.016.
Baron R, Kneissel M. WNT signaling in bone homeostasis and disease: from human mutations to treatments. Nat Med. 2013;19:179–92. https://doi.org/10.1038/nm.3074.
Augustin I. Wnt signaling in skin homeostasis and pathology. J Dtsch Dermatol Ges. 2015;13:302–6. https://doi.org/10.1111/ddg.12620.
Huang Y, Chen H, Zhang H, Lu Y, Qin C. FAM20C plays a critical role in the development of mouse vertebra. Spine J. 2022;22:337–48. https://doi.org/10.1016/j.spinee.2021.07.022.
Akiyama H, Lyons JP, Mori-Akiyama Y, Yang X, Zhang R, Zhang Z, et al. Interactions between Sox9 and beta-catenin control chondrocyte differentiation. Genes Dev. 2004;18:1072–87. https://doi.org/10.1101/gad.1171104.
Goll DE, Thompson VF, Li H, Wei W, Cong J. The calpain system. Physiol Rev. 2003;83:731–801. https://doi.org/10.1152/physrev.00029.2002.
Franco SJ, Huttenlocher A. Regulating cell migration: calpains make the cut. J Cell Sci. 2005;118:3829–38. https://doi.org/10.1242/jcs.02562.
Averna M, de Tullio R, Passalacqua M, Salamino F, Pontremoli S, Melloni E. Changes in intracellular calpastatin localization are mediated by reversible phosphorylation. Biochem J. 2001;354:25–30. https://doi.org/10.1042/0264-6021:3540025.
Loonat AA, Martin ED, Sarafraz-Shekary N, Tilgner K, Hertz NT, Levin R, et al. p38γ MAPK contributes to left ventricular remodeling after pathologic stress and disinhibits calpain through phosphorylation of calpastatin. Faseb J. 2019;33:13131–44. https://doi.org/10.1096/fj.201701545R.
Lawson CD, Ridley AJ. Rho GTPase signaling complexes in cell migration and invasion. J Cell Biol. 2018;217:447–57. https://doi.org/10.1083/jcb.201612069.
Tang Y, He Y, Zhang P, Wang J, Fan C, Yang L, et al. LncRNAs regulate the cytoskeleton and related Rho/ROCK signaling in cancer metastasis. Mol Cancer. 2018;17:77. https://doi.org/10.1186/s12943-018-0825-x.
Villalonga P, Ridley AJ. Rho GTPases and cell cycle control. Growth Factors. 2006;24:159–64. https://doi.org/10.1080/08977190600560651.
Aspenström P. The intrinsic GDP/GTP exchange activities of Cdc42 and Rac1 are critical determinants for their specific effects on mobilization of the actin filament system. Cells. 2019. https://doi.org/10.3390/cells8070759.
Millar FR, Janes SM, Giangreco A. Epithelial cell migration as a potential therapeutic target in early lung cancer. Eur Respir Rev. 2017. https://doi.org/10.1183/16000617.0069-2016.
Storr SJ, Carragher NO, Frame MC, Parr T, Martin SG. The calpain system and cancer. Nat Rev Cancer. 2011;11:364–74. https://doi.org/10.1038/nrc3050.
Carragher NO, Levkau B, Ross R, Raines EW. Degraded collagen fragments promote rapid disassembly of smooth muscle focal adhesions that correlates with cleavage of pp125(FAK), paxillin, and talin. J Cell Biol. 1999;147:619–30. https://doi.org/10.1083/jcb.147.3.619.
Carragher NO, Fincham VJ, Riley D, Frame MC. Cleavage of focal adhesion kinase by different proteases during SRC-regulated transformation and apoptosis. Distinct roles for calpain and caspases. J Biol Chem. 2001;276:4270–5. https://doi.org/10.1074/jbc.M008972200.
Konze SA, van Diepen L, Schroder A, Olmer R, Moller H, Pich A, et al. Cleavage of E-cadherin and beta-catenin by calpain affects Wnt signaling and spheroid formation in suspension cultures of human pluripotent stem cells. Mol Cell Proteomics. 2014;13:990–1007. https://doi.org/10.1074/mcp.M113.033423.
Li G, Iyengar R. Calpain as an effector of the Gq signaling pathway for inhibition of Wnt/beta -catenin-regulated cell proliferation. Proc Natl Acad Sci U S A. 2002;99:13254–9. https://doi.org/10.1073/pnas.202355799.
Goni-Oliver P, Lucas JJ, Avila J, Hernandez F. N-terminal cleavage of GSK-3 by calpain: a new form of GSK-3 regulation. J Biol Chem. 2007;282:22406–13. https://doi.org/10.1074/jbc.M702793200.
Ma S, Liu S, Huang Q, Xie B, Lai B, Wang C, et al. Site-specific phosphorylation protects glycogen synthase kinase-3β from calpain-mediated truncation of its N and C termini. J Biol Chem. 2012;287:22521–32. https://doi.org/10.1074/jbc.M111.321349.
VanderVorst K, Dreyer CA, Konopelski SE, Lee H, Ho HH, Carraway KL 3rd. Wnt/PCP signaling contribution to carcinoma collective cell migration and metastasis. Cancer Res. 2019;79:1719–29. https://doi.org/10.1158/0008-5472.CAN-18-2757.
Chen S, Liu Y, Zhang K, Chen L. CST1 promoted gastric cancer migration and invasion through activating Wnt pathway. Cancer Manag Res. 2021;13:1901–7. https://doi.org/10.2147/cmar.S277770.
Polakis P. Wnt signaling in cancer. Cold Spring Harb Perspect Biol. 2012. https://doi.org/10.1101/cshperspect.a008052.
Xie H, Ma Y, Li J, Chen H, Xie Y, Chen M, et al. WNT7A promotes EGF-induced migration of oral squamous cell carcinoma cells by activating β-Catenin/MMP9-mediated signaling. Front Pharmacol. 2020;11:98. https://doi.org/10.3389/fphar.2020.00098.
Lan L, Wang W, Huang Y, Zhao C, Bu X. WNT7A Overexpression Inhibits Growth and Migration of Hepatocellular Carcinoma via the β-Catenin Independent Pathway. Biomed Res Int. 2019;2019:3605950. https://doi.org/10.1155/2019/3605950.
Zuo H, Yang D, Wan Y. Fam20C regulates bone resorption and breast cancer bone metastasis through osteopontin and BMP4. Cancer Res. 2021;81:5242–54. https://doi.org/10.1158/0008-5472.Can-20-3328.
Miao N, Zhan Y, Xu Y, Yuan H, Qin C, Lin F, et al. Loss of Fam20c causes defects in the acinar and duct structure of salivary glands in mice. Int J Mol Med. 2019;43:2103–17. https://doi.org/10.3892/ijmm.2019.4126.
Cho YD, Kim WJ, Yoon WJ, Woo KM, Baek JH, Lee G, et al. Wnt3a stimulates Mepe, matrix extracellular phosphoglycoprotein, expression directly by the activation of the canonical Wnt signaling pathway and indirectly through the stimulation of autocrine Bmp-2 expression. J Cell Physiol. 2012;227:2287–96. https://doi.org/10.1002/jcp.24038.
Zhu M, Tang D, Wu Q, Hao S, Chen M, Xie C, et al. Activation of beta-catenin signaling in articular chondrocytes leads to osteoarthritis-like phenotype in adult beta-catenin conditional activation mice. J Bone Miner Res. 2009;24:12–21. https://doi.org/10.1359/jbmr.080901.
Wang X, Wang J, Yuan B, Lu Y, Feng JQ, Qin C. Overexpression of Dmp1 fails to rescue the bone and dentin defects in Fam20C knockout mice. Connect Tissue Res. 2014;55:299–303. https://doi.org/10.3109/03008207.2014.923414.
Metwally E, Zhao G, Wang Q, Zhang YQ. Ttm50 facilitates calpain activation by anchoring it to calcium stores and increasing its sensitivity to calcium. Cell Res. 2021;31:433–49. https://doi.org/10.1038/s41422-020-0388-4.
Vandenberghe M, Raphael M, Lehen’kyi V, Gordienko D, Hastie R, Oddos T, et al. ORAI1 calcium channel orchestrates skin homeostasis. Proc Natl Acad Sci U S A. 2013;110:E4839-4848. https://doi.org/10.1073/pnas.1310394110.
Banziger C, Soldini D, Schutt C, Zipperlen P, Hausmann G, Basler K. Wntless, a conserved membrane protein dedicated to the secretion of Wnt proteins from signaling cells. Cell. 2006;125:509–22. https://doi.org/10.1016/j.cell.2006.02.049.
Zhang X, Abreu JG, Yokota C, MacDonald BT, Singh S, Coburn KL, et al. Tiki1 is required for head formation via Wnt cleavage-oxidation and inactivation. Cell. 2012;149:1565–77. https://doi.org/10.1016/j.cell.2012.04.039.
Stanger BZ. Cellular homeostasis and repair in the mammalian liver. Annu Rev Physiol. 2015;77:179–200. https://doi.org/10.1146/annurev-physiol-021113-170255.
Yamaguchi M. Role of regucalcin in maintaining cell homeostasis and function (review). Int J Mol Med. 2005;15:371–89.
We thank all of the authors listed in this manuscript.
This work was supported by National Natural Science Foundation of China (Grant Nos. 81801040, 82001488, 81870736).
Ethics approval and consent to participate
All animal procedures in this study were reviewed and approved by the Ethics Committee at The Second Affiliated Hospital of Harbin Medical University (Harbin, China; approved protocol nos. SYDW2021-079) and performed in strict accordance with the National Institute of Health Guide for the Care and Use of Laboratory Animals.
Consent for publication
There was no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Fig. S1. Quality control for ATAC-seq samples generated in this study.
Fig. S2. Cell proliferation assay of OB Fam20cf/f and OB Fam20cKO vitro.
Fig. S3 Cell proliferation assay and cytotoxicity assay by cell counting kit-8 in Fam20cKO treated with CaCl2.
Fig. S4 Pathway map of overlapped genes from ATAC-seq and RNA-seq.
Fig. S5 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis based on ATAC-seq.
Fig. S6 The changes of Calpastatin/Calpain proteolysis system in human 293T cells after Fam20c knock out
Table S1 ATAC-seq sequencing data quality control statistics.
Table S2 Statistics of the alignment results of Reads on the reference genome sequence.
Table S3 Gene expression and Q value of RNA-seq.
Table S4 127 overlapped up-regulated differentially expressed genes between ATAC-seq and RNA-seq.
Table S5 109 overlapped down-regulated differentially expressed genes between ATAC-seq and RNA-seq.
Table S6 Biological Processes of Gene Ontology (GO) enrichment analysis based on the corresponding overlapping gene between ATAC-seq and RNA-seq.
Table S7 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis based on the corresponding overlapping gene between ATAC-seq and RNA-seq.
Table S8 The detailed information of differentially expressed peptides and proteins between OB Fam20cf/f and OB Fam20cKO in proteomics.
Table S9 The detailed information of differentially expressed phosphorylated peptides and proteins between OB Fam20cf/f and OB Fam20cKO in phosphoproteomics.
About this article
Cite this article
Liu, X., Jiang, L., Zhang, W. et al. Fam20c regulates the calpain proteolysis system through phosphorylating Calpasatatin to maintain cell homeostasis. J Transl Med 21, 417 (2023). https://doi.org/10.1186/s12967-023-04275-4