Skip to main content
Wiley Open Access Collection logoLink to Wiley Open Access Collection
. 2024 Jul 5;38(13):e23779. doi: 10.1096/fj.202400449RR

The G9a histone methyltransferase represses osteoclastogenesis and bone resorption by regulating NFATc1 function

Koichiro Komatsu 1, Hisashi Ideno 1, Kazuhisa Nakashima 1, Nobuyuki Udagawa 2, Yasuhiro Kobayashi 3, Hiroshi Kimura 4, Makoto Tachibana 5, Teruhito Yamashita 3,✉, Akira Nifuji 1,✉
PMCID: PMC13687911  PMID: 38967255

Abstract

Epigenetic modifications affect cell differentiation via transcriptional regulation. G9a/EHMT2 is an important epigenetic modifier that catalyzes the methylation of histone 3 lysine 9 (H3K9) and interacts with various nuclear proteins. In this study, we investigated the role of G9a in osteoclast differentiation. When we deleted G9a by infection of Cre‐expressing adenovirus into bone marrow macrophages (BMMs) from G9a fl/fl (Ehmt2 fl/fl ) and induced osteoclastic differentiation by the addition of macrophage colony‐stimulating factor (M‐CSF) and receptor activator of NF‐κB ligand (RANKL), the number of TRAP‐positive multinucleated osteoclasts significantly increased compared with control. Furthermore, the mRNA expression of osteoclast markers, TRAP, and cathepsin K, and to a lesser extent, NFATc1, a critical transcription factor, increased in G9a KO cells. Infection of wild‐type (WT) G9a‐expressing adenovirus in G9a KO cells restored the number of TRAP‐positive multinucleated cells. In G9a KO cells, increased nuclear accumulation of NFATc1 protein and decreased H3K9me2 accumulation were observed. Furthermore, ChIP experiments revealed that NFATc1 binding to its target, Ctsk promoter, was enhanced by G9a deletion. For in vivo experiments, we created G9a conditional knock‐out (cKO) mice by crossing G9a fl/fl mice with Rank Cre/+ (Tnfrsf11a Cre/+ ) mice, in which G9a is deleted in osteoclast lineage cells. The trabecular bone volume was significantly reduced in female G9a cKO mice. The serum concentration of the C‐terminal telopeptide of type I collagen (CTX), a bone‐resorbing indicator, was higher in G9a cKO mice. In addition, osteoclasts differentiated from G9a cKO BMMs exhibited greater bone‐resorbing activity. Our findings suggest that G9a plays a repressive role in osteoclastogenesis by modulating NFATc1 function.

Keywords: bone development, bone resorption, cell differentiation, methyltransferases, osteoclast


G9a/EHMT2, an important epigenetic modifier that catalyzes the methylation of histone 3 lysine 9 (H3K9), plays a repressive role in osteoclast formation and subsequent bone resorption by regulating NFATc1 function in osteoclast progenitor cells.

graphic file with name FSB2-38-e23779-g001.webp


Abbreviations

Ad

adenovirus

BMMs

bone marrow macrophages

BS/BV

bone surface per bone volume

BS/TV

bone surface per tissue volume

BV/TV

bone volume per tissue volume

CAG

CMV early enhancer/chicken β actin

ChIP

chromatin immunoprecipitation

cKO

conditional knockout

Ct.Ar

cortical bone area

Ct.Pm

cortical bone perimeter

CTX

cross‐linked C‐telopeptide of type I collagen

DAPI

4',6‐diamidino‐2‐phenylindole

Ec.Pm

endocortical perimeter

EGFP

enhanced green fluorescent protein

FBS

fetal bovine serum

GFP

green fluorescent protein

Il2ra

interleukin 2 receptor subunit alpha

IRES

internal ribosome entry site

Ma.Ar

bone marrow area

M‐CSF

macrophage colony‐stimulating factor

MNCs

multinucleated cells

N.Oc/BS

number of osteoclasts per bone surface

Oc.S/BS

surface area of osteoclasts per bone surface

RANK

receptor activator of NF‐κB

RANKL

receptor activator of NF‐κB ligand

ROI

region of interest

sRANKL

soluble RANK ligand

Tb.N

number of trabecular bones

Tb.Sp

space between trabecular bones

Tb.Th

trabecular bone thickness

TBS

Tris‐buffered saline

TRAP

tartrate‐resistant acid phosphatase

μCT

micro‐computed tomography

1. INTRODUCTION

Cellular differentiation is regulated by a combination of transcription factors and positively and/or negatively interacting cofactors. This interaction is affected by posttranslational and conformational modifications of the interacting chromatin components. Posttranslational modification of histones is critical for establishing and maintaining the chromatin structure. Histones are posttranscriptionally modified through modifications including the methylation of lysine residues; methylation status is governed by two kinds of modifiers: methyltransferases and demethylases. Histone H3 lysine 9 dimethylation, H3K9me2, often accumulates in euchromatic regions of the genome as well as in heterochromatin regions. Genome structure analysis has shown that H3K9me2 in euchromatic compartments is regulated by three methyltransferases: G9a/EHMT2, GLP, and SETDB1. 1 G9a and H3K9me2 maintain topologically associated domain (TAD) boundaries and chromatin loop structures that affect chromatin accessibility and transcriptional silencing. 2

Recent studies have shown that the interactions between G9a and other transcription factors play critical roles in cell differentiation. In hematopoietic stem cells, G9a promotes lineage commitment via H3K9me2 accumulation and suppresses non‐hematopoietic gene expression. 3 In T cell lineages, G9a suppresses the expression of CD4+ genes during switching from CD4+CD8+ to CD4−CD8+ T cells, which requires the formation of the G9a‐RUNX3 complex. 4 G9a is recruited by PRMD1 (Blimp‐1) to silence Il2ra and Cd27 expression in CD8+ T cells. 5 During tumorigenesis, MYC interacts with G9a to repress target gene transcription. 6 The G9a‐Runx2 interaction is required for proper bone development and osteoblastic differentiation. 7 Thus, G9a interacts with various nuclear proteins via H3K9me2‐dependent and/or ‐independent mechanisms, regulating cellular differentiation.

Bone is maintained by the balance between bone resorption and formation by osteoclasts and osteoblasts. Excessive bone resorption leads to pathological conditions such as osteoporosis and tooth loss. One cause for increased bone resorption is activation of osteoclasts, which are generated from hemopoietic progenitors of the monocyte– macrophage lineage. The progenitor cells differentiate multinucleated osteoclasts by addition of macrophage colony‐stimulating factor (M‐CSF) and receptor activator of NF‐κB ligand (RANKL). NFATc1, which is upregulated during osteoclast differentiation, has been identified as a critical transcription factor and shown to be essential for osteoclast development in vitro and in vivo. 8 , 9 NFATc1 is translocated to the nucleus and is assembled on multiple NFAT‐binding motifs in the promoter regions of osteoclast‐specific genes such as TRAP (Acp5) and cathepsin K (Ctsk). 10 The function of NFATc1 is determined by several mechanisms such as binding to its partner and nuclear translocation. Epigenetic regulation of NFATc1 may be also important, however, its association with epigenetic modifiers in osteoclastic differentiation has not been elucidated.

Histone modifications and their associated enzymes regulate osteoclastic differentiation. G9a catalyzes monomethylation of H3K27 and facilitates MMP‐9‐dependent H3N‐terminal tail (H3NT) proteolysis, leading to osteoclastic differentiation. 11 , 12 In another study, one of G9a inhibitors, BIX01294, inhibits osteoclastic differentiation in RAW264.7 cells. 13 The accumulation of trimethylated H3K27 (H3K27me3) around the transcription start site of the Nfatc1 gene decreases with RANKL‐induced osteoclast differentiation. Concomitantly, the expression of the H3K27me3 demethylase Jumonji domain‐containing 3 (Jmjd3) is induced during osteoclast differentiation. 14 Enhancer of zeste homolog 2 (EZH2), a histone methyltransferase that catalyzes H3K27me3, represses the expression of negative regulators MafB during osteoclastic differentiation. 15 EZH2 deletion in mature macrophages increases bone volume in vivo. 16

As G9a is involved in the differentiation of the hematopoietic lineage including osteoclast progenitor cells, we hypothesized that G9a plays a pivotal role during osteoclastic differentiation. To explore its role, we first examined the effects of G9a deletion on the in vitro osteoclastic differentiation of bone marrow macrophages (BMMs) isolated from G9a fl/fl mice. We also investigated its effects on bone resorption in vivo in G9a cKO mice by crossing G9a fl/fl mice with RANK‐Cre mice.

2. MATERIALS AND METHODS

2.1. Reagents

Human M‐CSF (Leukoprol) and soluble human recombinant RANKL (sRANKL) were purchased from Kyowa Kirin (Tokyo, Japan) and Oriental Yeast (Tokyo, Japan), respectively.

The following antibodies were used: rabbit polyclonal anti‐human NFATc1 (Abcam, 25961) and mouse monoclonal anti‐G9a (PP‐A8602A‐00; Perseus Proteomics, Tokyo, Japan). Mouse monoclonal anti‐H3K9me2 antibodies were prepared by H. Kimura. 17 Following secondary antibodies were obtained from Invitrogen, Carlsbad, CA, USA: Goat Alexa 488‐labeled anti‐mouse IgG(H+L) (A11001), goat Alexa 405‐labeled anti‐rabbit IgG(H+L) (A31556), and goat Alexa 594‐labeled anti‐mouse IgG(H+L) (A11032). All chemical‐grade reagents were obtained from WAKO Chemicals (Osaka, Japan).

2.2. Mice

To generate G9a fl/fl (Ehmt2 fl/fl ) mice, G9a deleted in ES cells (TT2 cell line) were injected into ICR mouse embryos. 18 G9a fl/+ mice were sequentially backcrossed with the C57BL/6 background mice, and offspring later than F15 generation were used for further studies. To generate Rank Cre/+ (Tnfrsf11a Cre/+ ) mice, an eGFP‐fused Cre recombinase cassette was inserted after the STOP codon of Tnfrsf11a (RANK‐encoding gene) followed by an IRES stuffer. 19 This Rank Cre/+ line was backcrossed to C57BL/6 background mice for more than 10 generations. G9a cKO (osteoclast lineage‐specific G9a gene deletion: G9a fl/fl ;RANK Cre/+ ) and control (G9a +/+ ;RANK Cre/+ ) mice were generated by crossing G9a fl/+ ;RANK Cre/+ parents. In vivo experiments were conducted on 10‐week‐old mice. All animal experiments were approved by the Institutional Animal Care Committee, the Recombination Experiment and Biosafety Committee of the Tsurumi University School of Dental Medicine, and the Animal Management Committee of Matsumoto Dental University.

2.3. In vitro osteoclastogenesis

Mouse bone marrow macrophages (BMMs) were prepared as described previously. 20 Briefly, bone marrow cells were isolated from the tibiae of 4‐ to 44‐week‐old male (n = 15) and 16‐ and 44‐week‐old female G9a fl/fl mice (n = 2) and from C57BL/6 WT mice of similar age. Bone marrow cells from G9a fl/fl and WT mice were cultured for 16 h in alpha‐modified minimal essential medium (αMEM; Life Technologies, Tokyo, Japan) containing 10% fetal bovine serum (FBS), 1% antibiotics (100 mg/mL penicillin and 100 IU/mL streptomycin; Sigma, St. Louis, MO, USA), and 30 ng/mL M‐CSF (Figure 1). Nonadherent cells were collected and used as BMMs. BMMs (1.5 × 105 cells/well) were cultured for 2 or 3 days in the presence of 50 ng/mL M‐CSF in 48‐well culture plates. After removing nonadherent cells by pipetting, adherent BMMs from G9a fl/fl and WT mice were used as osteoclast precursors and infected with an adenovirus expression vector‐harboring Cre recombinase at a multiplicity of infection (MOI = 100) for 24 h (Figure 1), as described previously. 21 The cells were further cultured for 4 days in the presence of 50 ng/mL M‐CSF and 200 ng/mL sRANKL for the osteoclast formation assay (Figure 1).

FIGURE 1.

FIGURE 1

G9a deletion in male BMMs promoted osteoclast differentiation. (A) G9a and osteoclast marker mRNA expressions at day 0 and day 4 after RANKL treatment, as determined by qPCR. The arbitrary value of each transcript was determined by normalized to Gapdh. Significant differences from the respective day 0 group, ***p < .001 (Student's t‐test). (B) Immunofluorescent analysis of G9a protein expression in WT and G9a‐deleted cells on day 4 after RANKL treatment. G9a expression was significantly reduced after Ad‐Cre infection (Ad/Cre) in G9a fl/fl BMMs. Scale bar: 50 μm. (C) TRAP staining on day 4 after RANKL treatment and G9a deletion. Number of TRAP‐positive multinucleated cells (MNCs) increased after G9a deletion (G9a fl/fl BMMs+ Ad/Cre). Scale bar: 50 μm. (D) G9a and osteoclast marker mRNA expressions analyzed by qPCR on days 0, 2, and 4 after RANKL treatment in WT and G9a KO cells. Each column and bar represent the mean and 1 SD. The Acp5 and Ctsk values on day 0 were as follows: Acp5; 0.00741 ± 0.00027 SD (WT), 0.00913 ± SD 0.00002 (KO), and Ctsk; 0.00178 ± SD 0.00006 (WT), 0.0042 ± SD 0.00005 (KO), respectively. Significant differences in the respective BMMs group, *p < .05, **p < .01, ***p < .001. (E) Reintroduction of G9a into G9a KO cells. To create G9a‐deleted BMMs, Ad‐Cre was infected into BMMs isolated from G9a fl/fl mice (Ad/Cre). Then, Ad/G9a was reintroduced into G9a‐deleted BMMs (Ad/Cre and Ad/G9a) to rescue G9a deletion. Infection of GFP‐expressing adenovirus was served as control (Ad/GFP). (a) TRAP staining on day 4 after RANKL treatment. Increase in TRAP‐positive MNC numbers in G9a‐deleted BMMs (Ad/Cre) was abrogated by Ad/G9a infection (Ad/Cre + Ad/Cre). (b) G9a expression shown by qPCR. G9a expression was suppressed by Ad/Cre and restored by Ad/Cre and Ad/G9a infection. Each column and bar represent the mean and 1 SD. Significant differences between any two groups, ***p < .001 (Scheffé method).

2.4. TRAP staining

Cells were fixed with 3.7% neutrally buffered formaldehyde solution for 10 min. The cells were then rinsed with Tris‐buffered saline (TBS) three times, and the fixed cells were permeabilized with a 1:1 solution of acetone and ethanol for 1 min and stained using a TRAP staining kit (Cosmo Bio, Tokyo, Japan). TRAP‐positive multinuclear cells (MNCs) (>3 nuclei/cell) were considered as mature osteoclasts.

2.5. Preparation of adenovirus vectors

The production of adenovirus was performed as described previously 22 according to the manufacturer's protocol (ViraPower Adenoviral Expression System; Thermo Fisher Scientific). An adenovirus vector expressing Cre‐recombinase (AxCANCre; Ad/Cre) was obtained from RIKEN BRC (Tsukuba, Japan). An adenovirus vector expressing enhanced green fluorescent protein (EGFP) (Ad/GFP), which contains CAG‐IRES‐EGFP, was prepared as described. 22 Another adenovirus vector expressing G9a was constructed as follows: a fragment of CAG‐G9a was subcloned into the pENTR vector (Thermo Fisher Scientific). Then, pENTR‐CAG‐G9a‐IRES‐EGFP was subcloned into pAd/PL‐DEST by a site‐specific recombination.

2.6. Bone resorption assay (pit formation assay)

BMMs were further cultured with M‐CSF and sRANKL for 4–6 days in 60 mm RepCell plates (temperature‐responsive polymer‐coated plates; CellSeed Inc., Tokyo, Japan) to allow the development of multinucleated osteoclasts. Osteoclasts were harvested and seeded at approximately 100 osteoclasts per 4‐mm‐square dentin slice in 48‐well plates with 0.5 mL medium supplemented with M‐CSF and sRANKL for 48 h. Dentin slices were stained with TRAP to confirm the presence of osteoclasts. After the removal of cells from the surface, these slices were re‐stained with toluidine blue to detect resorption pits. The stained pit images were binarized, and the resorption area was measured using the ImageJ software (NIH).

2.7. Reverse transcription and real‐time PCR

RNA was isolated using an RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. RNA was reverse transcribed using oligo‐dT primers and reverse transcriptase (Superscript III; Invitrogen) to synthesize cDNA. After cDNA synthesis, qPCR was performed as previously described. 22 Three or more independent cultures were subjected to qPCR and representative results are shown. The primer sequences used for each PCR experiment are listed in Table 1.

TABLE 1.

Primer sequences used in qPCR.

Gene Forward (5′‐>3′) Reverse (5′‐>3′) Accession number
Gapdh GCCAAACGGGTCATCATCTC GTCATGAGCCCTTCCACAAT NM_001411843.1
G9a CTCTACCGGACTGCCAAGAT CTCGGCATCAGAGATCAGC NM_001403869.1
Nfatc1 TCTTCCGAGTTCACATCC ACAGCACCATCTTCTTCC XM_036161029.1
Acp5 GCAGCCAAGGAGGACTAC CCCACTCAGCACATAGCC XM_006509946.3
Ctsk CAGAACGGAGGCATTGAC CGATGGACACAGAGATGG NM_007802.4

2.8. Immunofluorescence staining

BMMs were cultured on coverslips or chamber slides for osteoclast formation as described above. The cells were fixed with 3.7% neutrally buffered formaldehyde solution for 10 min. After rinsing with TBS three times, the fixed cells were permeabilized with a 1:1 solution of acetone and ethanol for 1 min and then blocked with 5% goat serum in TBS. The cells were then incubated with mouse anti‐G9a, rabbit anti‐NFATc1, or mouse anti‐H3K9me2 17 antibodies, with the appropriate Alexa Fluor‐conjugated secondary antibody, and counterstained with DAPI.

NFATc1‐positive nucleus was defined as being its positive immunofluorescent nuclear area occupying more than 50% of its nucleus area. The image thresholding method of Image J (version 1.53k, NIH, USA) was used to count NFATc1‐positive nuclei in fluorescent images and the total number of nuclei in DIC images, respectively. Then, the percentage of NFATc1‐positive nuclei to total nuclei of multinucleated osteoclasts was calculated.

2.9. Chromatin immunoprecipitation (ChIP) assay

The ChIP assay was performed as previously described. 7 BMMs were cultured at 4 × 106 cells per 6 cm tissue culture dish and treated with 1% formaldehyde for 10 min at room temperature. The crosslinked chromatin was sheared to an average length of 200–1000 bp by sonication with a sonicator (Bioruptor UCD‐200; Cosmobio, Tokyo, Japan). The precleaned chromatin was used as the input control. Immunoprecipitation was carried out overnight at 4°C against the following antibodies: NFATc1 (4 μg, sc‐7294X, Santa Cruz, USA), G9a (1 μg, 68851, Cell Signaling Technology, MA, USA), normal mouse IgG (4 μg, 015‐000‐003, Jackson ImmunoResearch, PA, USA), and normal rabbit IgG (1 μg, ab46540, Abcam, Cambridge, UK). Immunoprecipitation was performed using Dynabeads M‐280 sheep anti‐mouse IgG (11202D, Thermo Fisher Scientific) and Dynabeads M‐280 sheep anti‐rabbit IgG (11204D, Thermo Fisher Scientific). After repeated rinses, protein–DNA cross‐linking was reversed, and the precipitated DNA was treated with RNAse at 37°C and with proteinase K at 55°C, and then purified by phenol–chloroform extraction. DNA was analyzed by real‐time PCR (qPCR). Data obtained for each primer pair were expressed as a percentage of the input DNA. Primer sets used for the mouse Ctsk promoters 23 are as follows:

N1 site, 5′‐CCCCCAAAGTCAGTCAGATG‐3′ and 5′‐GGTAAGGATTGCGGAAGTCA‐3′; N3 site, 5′‐CTAGCCTTTCCTCCCCTCTC‐3′ and 5′‐CCCTAGTTGTCTCCATTCTTCC‐3′.

2.10. Microcomputed tomography (μCT) measurement

X‐ray images of the distal end of the femora were obtained using the μCT scanner, ScanXmate‐A080 (Comscantecno, Kanagawa, Japan), with the following parameters: 22 kV and 250 μA. Three‐dimensional images reconstructed from the multi‐slice images, which were taken at the ROI of 1–3 mm (20 μm slices) away from the growth plate, were analyzed using the TRI/3D‐BON software (Ratoc System Engineering, Tokyo, Japan). For cortical bone analysis, images of horizontal cross‐sections 3 mm away from the femoral growth plate were analyzed by ImageJ.

2.11. Bone histomorphometry

Left femora of 10‐week‐old mice were fixed overnight in 4% paraformaldehyde and soaked in a series of sucrose concentrations up to a maximum of 100% prior to mounting in OTC compound in a hexane cooling bath at −70°. Ten μm serial sections were prepared using a cryotome (CM‐300S, Leica, Germany) and stained for tartrate‐resistant acid phosphatase (TRAP) to identify osteoclast cells. The sections were counterstained with methyl green. Images were captured using a microscope (BZ‐X710, Keyence, Osaka, Japan). Osteoclast surface (Oc.S/BS) and osteoclast number (N.Oc/BS) were measured on the trabecular bone in the ROI area (0.75 mm square), 100 μm away from the growth plate cartilage, by tracing cell bodies and contact surfaces on the digitizer in images of tissue sections using ImageJ. Bone measurement parameters were described according to the report of the ASBMR Histomorphometry Nomenclature Committee. 24

2.12. Measurement of serum CTX and osteocalcin

The cross‐linked C‐telopeptide of type I collagen (CTX) in the serum, which reflects degraded collagen derived from resorbed bone, was measured using a RatLap ELISA kit (AC‐06F1, Immunodiagnostic Systems, UK) according to the manufacturer's instructions. Serum Gla‐type osteocalcin levels, an indicator of bone formation, were measured by EIA method according to the manual for the Mouse Gla‐Osteocalcin High Sensitive EIA Kit (MK127, Takara‐Bio, Kyoto, Japan).

2.13. Statistics

The mean values of measures in the two groups of gene expressions, TRAP‐positive MNCs numbers, ratio of NFATC1‐positive nuclei number/total nuclei number of MNOC, and CHIP assay were compared using Student's t‐test with equal variances or Welch test with unequal variances. The mean values in the three or more groups of gene expressions were examined using a one way of variance (ANOVA) and a post hoc test (Scheffé method). Those measures are expressed as the mean and standard deviations (SD).

The statistical differences between the two groups of bone morphometric and pit size data were examined using Mann–Whitney U test. The former and latter measures were expressed as box‐and‐whisker diagrams and mean ± 1 SD, respectively.

3. RESULTS

3.1. G9a deletion in BMMs promoted osteoclast differentiation and G9a overexpression reversed it

First, we examined G9a expression during osteoclastic differentiation. We induced osteoclast differentiation in wild‐type (WT) BMMs in vitro by treatment with M‐CSF and RANKL. G9a mRNA expression slightly diminished after the addition of RANKL, while mRNA for osteoclast‐specific genes, Ctsk, Acp5, and Nfatc1, increased during osteoclastogenesis in vitro (Figure 1A).

To investigate the function of G9a in osteoclastic differentiation, we deleted G9a in BMMs in vitro. BMMs derived from G9a fl/fl mice were cultured with M‐CSF, and an adenovirus expressing Cre (Ad/Cre) was used to infect these cells 24 h before the addition of RANKL (see the experimental outline in Figure S1). We confirmed that G9a expression levels were low in BMMs at 2 days (mRNA) and 4 days (mRNA and protein) after RANKL treatment (Figure 1B,D) in G9a‐deleted cells. TRAP staining revealed that the numbers of TRAP‐positive MNCs increased after G9a deletion (Figure 1C). Furthermore, the mRNA expression of Acp5 and Ctsk, and to a lesser extent, Nfatc1, was increased in G9a KO cells (Figure 1D). Thus, G9a deletion in BMMs generated more mature osteoclasts in vitro. These effects were observed in BMMs derived from female mice (Figure S2) and male mice.

We reintroduced G9a cDNA using a recombinant adenovirus into G9a‐deleted BMMs. Ad/GFP‐expressing cells were used as the control. To create G9a‐deleted BMMs, Ad/Cre was infected into BMMs isolated from G9a fl/fl mice. Then, Ad/G9a was reintroduced into G9a‐deleted BMMs (Ad/Cre and Ad/G9a) to rescue G9a deletion. The reintroduction of G9a suppressed the increase in TRAP‐positive MNCs in G9a‐deleted BMMs (Figure 1E(a)). We confirmed that G9a expression was inhibited by Ad/Cre and restored by Ad/Cre and Ad/G9a infection (Figure 1E(b)).

3.2. Osteoclasts lacking G9a showed reduced H3K9me2 levels

One of the major substrates of G9a is H3K9me2. Next, we investigated the levels of H3K9me2 in G9a‐deleted and non‐deleted BMMs. Immunofluorescent staining of H3K9me2 was uniformly present in macrophage nuclei before induction of osteoclast differentiation (Figure 2A, rightmost panel). On day 2 when osteoclast differentiation was induced by addition of RANKL, H3K9me2 staining intensity was reduced in G9a‐deleted BMMs compared to WT BMMs (Figure 2B, right‐most panels). On day 4 after induction of differentiation, H3K9me2 was localized in the nuclei of WT mononuclear and multinuclear cells, whereas its localization was markedly reduced in G9a‐deleted osteoclasts (Figure 2C, rightmost panels). Quantification of H3K9me2 fluorescence intensity per DAPI‐positive nucleus showed significant reduction in G9a‐deleted cells on days 2 and 4 after differentiation (Figure 2D). Thus, the increased osteoclast differentiation due to G9a deletion was accompanied by a decrease in H3K9me2 level.

FIGURE 2.

FIGURE 2

G9a‐deleted osteoclasts showed decreased H3K9me2 levels. Immunostaining analysis of H3K9me2 in WT and G9a‐deleted male cells. Images of DIC, stained with DAPI and anti‐H3K9me2 on (A) day 0, (B) day 2, and (C) day 4. Scale bar: 20 μm. (D) Values of H3K9me2 fluorescence intensity/DAPI area in multinucleated cells (MNCs) on days 2 and 4 after RANKL treatment are shown. ***p < .001.

3.3. G9a suppressed NFATc1 activity

To address the molecular mechanisms underlying the suppressive roles of G9a in osteoclasts, we investigated the effect of G9a on the function of a key transcription factor of osteoclastogenesis, NFATc1. Because nuclear import is critical for NFATc1 function, we examined NFATc1 localization using immunofluorescent staining. We found that the nuclear localization of NFATc1 protein in multinucleated osteoclasts increased in G9a‐deleted osteoclasts (Figure 3A). We further analyzed the binding of NFATc1 to target regions in the genome. The ChIP experiment confirmed that the binding of NFATc1 to two Ctsk promoter regions, N1 and N3, increased after RANKL‐induced osteoclastic differentiation (data not shown). The binding of NFATc1 to N1 and N3 was further enhanced by G9a deletion (Figure 3B). Enhanced binding to the promoter was more significant in N3 than in N1. These data suggest that G9a suppressed NFATc1 function in the nucleus.

FIGURE 3.

FIGURE 3

G9a regulates NFATc1 activity. (A) Immunostaining analysis of NFATc1 in WT and G9a‐deleted male osteoclasts on day 4. WT: WT BMMs + Ad‐Cre, KO: G9a fl/fl BMMs + Ad‐Cre. Images of anti‐NFATc1, DIC, and anti‐NFATc1 merged with DIC (Merged) are shown. Scale bar. 50 μm. The percentages of NFATc1‐positive nuclei in MNCs containing more than nine nuclei are shown in the right panel. (B) ChIP assay using an anti‐NFATc1 antibody on the Ctsk promoters N1 and N3. Data for DNA obtained after ChIP was shown as percent against input. G9a deletion (KO) enhanced the binding of the NFATc1 protein to the Ctsk promoters N1 and N3 compared to WT osteoclasts differentiated from BMMs on day 4. Each column and bar represent the mean and 1 SD. Significant differences between the WT and KO groups, *p < .05 (Welch test).

3.4. G9a cKO mice showed increased bone resorption

Since in vitro experiments revealed promotion of osteoclast differentiation by G9a deletion, we investigated the effects of G9a deletion on bone resorption in vivo. G9a cKO mice were generated by crossing G9a fl/fl mice with RANK‐Cre (Rank Cre/+ ) mice, in which G9a was deleted from the osteoclast lineage. Micro‐CT analysis of the femur of 10‐week‐old female mice revealed reduced trabecular bone volume in G9a cKO mice as compared to that in WT mice (Figure 4A). G9a cKO mice also displayed reduced bone volume per tissue volume (BV/TV), bone surface per tissue volume (BS/TV), and number of trabecular bones (Tb.N) (Figure 4B). On the other hand, cortical bone area (Ct. Ar) and bone marrow area (Ma. Ar) were unchanged in horizontal cross sections of the femur from control and cKO female mice (Figure S3). No significant differences were evident in the outer and inner circumference length of the cortical bone (Ct.Pm and Ec.Pm, respectively) between the control and cKO mice (Figure S3). Thus, G9a cKO mice demonstrated accelerated bone loss in trabecular bone.

FIGURE 4.

FIGURE 4

Bone loss in G9a cKO female mice. (A) Representative μCT images of the distal end of the femur obtained from 10‐week‐old control (G9a +/+ ;RANK Cre/+ ) and G9a cKO (G9a fl/fl ;RANK Cre/+ ) mice. Scale bar: 1000 μm. (B) Significant differences in trabecular bone morphometric parameters in the femur of control (n = 10) and G9a cKO (n = 8) mice were evident. *p < .05 (Mann–Whitney U test). (C) Significantly increased serum levels of CTx, a bone resorption marker, in G9a cKO than in control mice. Serum levels of osteocalcin, a bone‐forming marker, were not changed. *p < .05 (Mann–Whitney U test). (D) TRAP staining of long bones. Representative sagittal sections of the distal end of the femur were obtained from 10‐week‐old control and G9a cKO mice. Right panels; osteoclast surface and numbers on long bone sections. Not significantly different osteoclast surface on the trabecular bone (Oc.S/BS) and number of osteoclasts per bone (N.Oc/BS) in the femur of G9a cKO and control mice. Scale bar: 100 μm. BS/BV, bone surface per bone volume; BS/TV, bone surface per tissue volume; BV/TV, bone volume per tissue volume; Tb.N, number of trabecular bones; Tb.Sp, space between trabecular bones; Tb.Th, trabecular bone thickness.

Serum levels of cross‐linked C‐telopeptide of type I collagen (CTx), a bone resorption marker, were increased twofold in G9a cKO mice as compared to that in control mice (Figure 4C, upper panel). In contrast, no significant change was observed in serum levels of osteocalcin, a marker of bone formation (Figure 4C, lower panel). We further analyzed the number of osteoclasts using sagittal sections of the distal femur from 10‐week‐old control and G9a cKO female mice. The number of osteoclasts per bone surface (N.Oc/BS) and the surface area of osteoclasts per bone surface (Oc.S/BS) were not significantly different between the two mouse groups (Figure 4D). These results suggest that bone resorption is enhanced in G9a cKO mice, resulting in bone loss. Increased bone resorption in the trabecular bone was observed in female mice but not in male mice (Figure S4A). Additionally, serum levels of cross‐linked C‐telopeptide of type I collagen (CTx) and osteocalcin were also not elevated in male mice (Figure S4B).

3.5. Osteoclasts from G9a cKO mice showed increased bone‐resorbing activity

In accordance with unchanged number of osteoclasts in G9a cKO mice, there was no difference in the number of TRAP+ MNCs differentiated from BMMs in G9a‐cKO and control mice (Figure S5). However, more heavily stained TRAP‐positive cells were observed in G9a cKO mice BMMs than in control mice BMM. Therefore, we examined the bone‐resorbing activity of RANKL‐induced osteoclasts in G9a cKO mice. To investigate the bone‐resorbing activity, we performed a pit formation assay by seeding the same number of TRAP‐positive osteoclasts onto dentin slices (Figure 5A). The pit area per slice and the number of pits increased in osteoclasts derived from G9a cKO than in control mice (Figure 5B,C). These results suggest that G9a cKO BMMs treated with RANKL generated multinucleated TRAP‐positive osteoclasts with increased bone‐resorbing activity.

FIGURE 5.

FIGURE 5

Osteoclasts derived from G9a cKO mice exhibit enhanced bone‐resorbing activity. (A) Bone resorption assay on cultivated osteoclasts derived from G9a cKO and control male mice. Upper panels: bone resorption pits were binarized and are indicated with black spots. Scale bar: 1000 μm. Lower panels: TRAP staining of bone slices before cell removal. Scale bar: 200 μm. (B) Increased pit area per slice in osteoclasts derived from G9a cKO and control mice cultured on dentin slices. *p < .05 (Mann–Whitney U test). (C) G9a cKO‐derived osteoclasts showing a greater area per resorption pit compared to controls as to the distribution of pit size.

4. DISCUSSION

In this study, we investigated the role of G9a in osteoclast differentiation and function. G9a deletion in BMMs yielded more multinucleated osteoclasts in vitro, and this effect was reversed by Ad/G9a overexpression in G9a‐deleted cells. G9a KO osteoclasts showed enhanced gene expression of Acp5 and Ctsk. Previous studies showed that inhibition of G9a using shRNA‐mediated knockdown, or the G9a inhibitor BIX01294, suppresses osteoclastic differentiation. 11 , 12 , 13 The discrepancy of these studies with our results in vitro may be accounted for by differences in the approaches and efficiencies of G9a depletion from G9a fl/fl bone marrow cells, with our results demonstrating more than 95% suppression.

Mechanistically, we showed that G9a KO cells increase the nuclear translocation of NFATc1 and the binding of NFATc1 to the promoter region of NFATc1 target gene, Ctsk. Osteoclast derived from G9a cKO mice showed higher bone‐resorbing activity. In line with the promoted osteoclastogenesis in vitro, G9a cKO mice showed lower bone volume in vivo. Thus, G9a regulates NFATc1 function in the nucleus and suppresses hyperactivated bone resorption.

How does G9a affect NFATc1 function and bone‐resorbing activity? One possible mechanism is through H3K9me2 status by G9a. Immunohistochemical analysis revealed that H3K9me2‐positive signals in multinucleated osteoclasts decreased with G9a deletion. It has been reported that a decrease in H3K9me2‐enriched regions is correlated with compartmental changes from a repressive to an active state, leading to transcriptional activation. 1 The reduction in H3K9me2 by G9a KO could lead to hyperactivated states in the regulatory loci where NFATc1 is functional. Another possibility is that association of G9a, H3K9me2, and B lymphocyte‐induced maturation protein‐1 (Blimp‐1) may affect osteoclastogenesis. G9a‐dependent H3K9me2 and Blimp‐1 maintain a repressive state in regulatory regions in CD8+ T cells. 5 Blimp‐1 and B cell lymphoma 6 (Bcl6) axis targets NFATc1 during osteoclastic differentiation. 25 Recently, a direct interaction between Blimp‐1 and NFATc1 has been shown. 26 Therefore, repressive complexes of G9a, H3K9me2, and Blimp‐1 may suppress NFATc1 function and this possibility needs to be tested in the future.

G9a cKO female mice showed reduced trabecular bone volume and increased serum CTX levels, indicating accelerated bone resorption in trabecular bone, but not cortical bone. However, the number of osteoclasts per bone surface area (N.Oc/BS) and the surface area of osteoclasts per bone surface area (Oc.S/BS) did not change significantly. The discrepancy in the results between RANK‐Cre‐expressed osteoclasts in vivo and adeno‐Cre‐transduced G9a‐deleted osteoclasts in vitro could be due to the deletion efficiency of Cre recombinase. RANK‐Cre‐expressed osteoclasts showed approximately 50% G9a expression of that of WT osteoclasts (not shown), whereas adeno‐Cre‐transduced G9a‐deleted osteoclasts showed over 95% reduction in vitro in G9a expression. Another explanation could be the accelerated cell turnover of osteoclasts in G9a cKO, which could generate more osteoclasts but disappear quickly in vivo, leading to unchanged cell numbers. Osteoclasts in G9a cKO likely gained more active resorbing capacity. In the pit‐forming assay, the total resorption area was significantly increased in the G9a KO than in the WT group when the same number of TRAP‐positive MNCs were seeded on dentine slices. The area per pit, that is, the resorption area per cell, also drastically increased in the osteoclasts derived from the KO group. This suggests that accelerated bone resorption in vivo in the G9a‐deleted mice may be explained by increased function per osteoclast.

Enhanced bone resorption in G9a cKO mice was observed in female mice but not in male mice. This may be because females have a lower bone mass than males, which may account for the differences in susceptibility to osteoporosis in the sexes. A recent study showed that RANKL induces a higher expression of NFATc1 and its target genes in female‐derived osteoclasts than in male‐derived osteoclasts. 27 Therefore, more activated NFATc1 in female osteoclasts may be susceptible to NFATc1 modulators. It has been reported that G9a functions as a coactivator of the endogenous estrogen receptor α (ERα) in breast cancer cells, thus G9a may affect gender‐dependent bone metabolisms through hormonal regulation. 28 In in vitro experiments, we performed G9a deletion using BMM from G9a fl/fl female and male mice and observed increased osteoclastic differentiation in both sexes. Thus, at least in vitro, differences between the sexes were not observed. Therefore, careful and substantial experiments are required to elucidate the mechanisms underlying the differences per sex due to G9a deletion on bone resorption.

A recent study showed that G9a is active in initiating tumors from cancer stem cells in various tumors. 29 The inhibition of G9a activity is a promising approach in preclinical studies on several kinds of tumors. 30 , 31 Nevertheless, our results indicate that G9a inhibition accelerates bone resorption. Thus, these unanticipated effects should be considered before proceeding to clinical trials applying G9a inhibitors.

Our results suggest that BMMs lacking G9a become hyperactivated osteoclasts with more nuclear accumulation of NFATc1. However, we could not show direct evidence of the involvement of NFATc1 in vivo in the phenotype of G9a cKO mice. This can be demonstrated by crossing G9a cKO mice with heterozygous or homozygous osteoclast‐specific NFATc1 knockout mice, which could reverse the phenotype. Future study is needed to explore this point.

5. CONCLUSION

In this study, we demonstrated that G9a deletion during osteoclastogenesis results in hyperactivated osteoclast formation and increased resorption, suggesting that G9a inhibits extravagant (over formation) osteoclast formation and bone resorption. Moreover, our findings suggest that G9a negatively regulates NFATc1 function. These results offer insights into the mechanism underlying bone remodeling, suggesting potential targets for devising therapeutic strategies for bone‐related disorders.

AUTHOR CONTRIBUTIONS

Koichiro Komatsu: performed the research, acquired the data, and manuscript writing; Hisashi Ideno: performed the research and acquired and analyzed the data; Kazuhisa Nakashima, Nobuyuki Udagawa, and Yasuhiro Kobayashi: analyzed and interpreted the data and provision of mechanistic insight; Hiroshi Kimura and Makoto Tachibana: methodology, resources; Teruhito Yamashita and Akira Nifuji: conception and design, experiment performance, provision of mechanistic insight, and manuscript writing.

FUNDING INFORMATION

This study was partially supported by JSPS KAKENHI (grant numbers: 22390344, 15K15679, 18K09515, 20K09897, 21K09823, 23K09407, and 23K09158).

DISCLOSURES

The authors declare that they have no conflicts of interest with the contents of this article.

Supporting information

Figure S1.

FSB2-38-e23779-s001.pptx (1.2MB, pptx)

ACKNOWLEDGMENTS

The authors would like to thank Ms A. Iwadare for her technical assistance with the experiments.

Komatsu K, Ideno H, Nakashima K, et al. The G9a histone methyltransferase represses osteoclastogenesis and bone resorption by regulating NFATc1 function. The FASEB Journal. 2024;38:e23779. doi: 10.1096/fj.202400449RR

Koichiro Komatsu and Hisashi Ideno contributed equally to this work.

Contributor Information

Teruhito Yamashita, Email: teruhito.yamashita@mdu.ac.jp.

Akira Nifuji, Email: nifuji-a@tsurumi-u.ac.jp.

DATA AVAILABILITY STATEMENT

The data that support the findings are available in the Methods section and Supplementary Data of this article.

REFERENCES

  • 1. Fukuda K, Shimura C, Miura H, et al. Regulation of mammalian 3D genome organization and histone H3K9 dimethylation by H3K9 methyltransferases. Commun Biol. 2021;4(1):571. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 2. Jiang Q, Ang JYJ, Lee AY, et al. G9a plays distinct roles in maintaining DNA methylation, retrotransposon silencing, and chromatin looping. Cell Rep. 2020;33(4):108315. [DOI] [PubMed] [Google Scholar]
  • 3. Chen X, Skutt‐Kakaria K, Davison J, et al. G9a/GLP‐dependent histone H3K9me2 patterning during human hematopoietic stem cell lineage commitment. Genes Dev. 2012;26(22):2499‐2511. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 4. Verbaro DJ, Sakurai N, Kim B, Shinkai Y, Egawa T. Cutting edge: the histone methyltransferase G9a is required for silencing of helper T lineage–associated genes in proliferating CD8 T cells. J Immunol. 2018;200(12):3891‐3896. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 5. Shin HM, Kapoor VN, Guan T, Kaech SM, Welsh RM, Berg LJ. Epigenetic modifications induced by blimp‐1 regulate CD8+ T cell memory progression during acute virus infection. Immunity. 2013;39(4):661‐675. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 6. Tu WB, Shiah YJ, Lourenco C, et al. MYC interacts with the G9a histone methyltransferase to drive transcriptional repression and tumorigenesis. Cancer Cell. 2018;34(4):579‐595.e8. [DOI] [PubMed] [Google Scholar]
  • 7. Ideno H, Nakashima K, Komatsu K, et al. G9a is involved in the regulation of cranial bone formation through activation of Runx2 function during development. Bone. 2020;137:115332. [DOI] [PubMed] [Google Scholar]
  • 8. Ishida N, Hayashi K, Hoshijima M, et al. Large scale gene expression analysis of osteoclastogenesis in vitro and elucidation of NFAT2 as a key regulator. J Biol Chem. 2002;277(43):41147‐41156. [DOI] [PubMed] [Google Scholar]
  • 9. Aliprantis AO, Ueki Y, Sulyanto R, et al. NFATc1 in mice represses osteoprotegerin during osteoclastogenesis and dissociates systemic osteopenia from inflammation in cherubism. J Clin Invest. 2008;118(11):3775‐3789. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 10. Udagawa N, Koide M, Nakamura M, et al. Osteoclast differentiation by RANKL and OPG signaling pathways. J Bone Miner Metab. 2021;39(1):19‐26. [DOI] [PubMed] [Google Scholar]
  • 11. Shin Y, Ghate NB, Moon B, Park K, Lu W, An W. DNMT and HDAC inhibitors modulate MMP‐9‐dependent H3 N‐terminal tail proteolysis and osteoclastogenesis. Epigenetics Chromatin. 2019;12(1):25. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 12. Kim K, Shin Y, Kim J, Ulmer TS, An W. H3K27me1 is essential for MMP‐9‐dependent H3N‐terminal tail proteolysis during osteoclastogenesis. Epigenetics Chromatin. 2018;11(1):23. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 13. Tsuda H, Zhao N, Imai K, Ochiai K, Yang P, Suzuki N. BIX01294 suppresses osteoclast diff erentiation on mouse macrophage‐like Raw264.7 cells. Bosn J Basic Med Sci. 2013;13(4):271‐275. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 14. Yasui T, Hirose J, Tsutsumi S, Nakamura K, Aburatani H, Tanaka S. Epigenetic regulation of osteoclast differentiation: possible involvement of Jmjd3 in the histone demethylation of Nfatc1. J Bone Miner Res. 2011;26(11):2665‐2671. [DOI] [PubMed] [Google Scholar]
  • 15. Adamik J, Pulugulla SH, Zhang P, et al. EZH2 supports osteoclast differentiation and bone resorption via epigenetic and cytoplasmic targets. J Bone Miner Res. 2020;35(1):181‐195. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 16. Caviness PC, Gai D, Lazarenko OP, Blackburn ML, Zhan F, Chen JR. Decreased bone resorption in Ezh2 myeloid cell conditional knockout mouse model. FASEB J. 2023;37(7):e23019. [DOI] [PubMed] [Google Scholar]
  • 17. Kimura H, Hayashi‐Takanaka Y, Goto Y, Takizawa N, Nozaki N. The organization of histone H3 modifications as revealed by a panel of specific monoclonal antibodies. Cell Struct Funct. 2008;33(1):61‐73. [DOI] [PubMed] [Google Scholar]
  • 18. Tachibana M, Nozaki M, Takeda N, Shinkai Y. Functional dynamics of H3K9 methylation during meiotic prophase progression. EMBO J. 2007;26(14):3346‐3359. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 19. Maeda K, Kobayashi Y, Udagawa N, et al. Wnt5a‐Ror2 signaling between osteoblast‐lineage cells and osteoclast precursors enhances osteoclastogenesis. Nat Med. 2012;18(3):405‐412. [DOI] [PubMed] [Google Scholar]
  • 20. Mizoguchi T, Muto A, Udagawa N, et al. Identification of cell cycle‐arrested quiescent osteoclast precursors in vivo. J Cell Biol. 2009;184(4):541‐554. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 21. Wada S, Ideno H, Shimada A, et al. H3K9MTase G9a is essential for the differentiation and growth of tenocytes in vitro. Histochem Cell Biol. 2015;144(1):13‐20. [DOI] [PubMed] [Google Scholar]
  • 22. Ideno H, Takanabe R, Shimada A, et al. Protein related to DAN and cerberus (PRDC) inhibits osteoblastic differentiation and its suppression promotes osteogenesis in vitro. Exp Cell Res. 2009;315(3):474‐484. [DOI] [PubMed] [Google Scholar]
  • 23. Balkan W, Martinez AF, Fernandez I, Rodriguez MA, Pang M, Troen BR. Identification of NFAT binding sites that mediate stimulation of cathepsin K promoter activity by RANK ligand. Gene. 2009;446(2):90‐98. [DOI] [PubMed] [Google Scholar]
  • 24. Dempster DW, Compston JE, Drezner MK, et al. Standardized nomenclature, symbols, and units for bone histomorphometry: a 2012 update of the report of the ASBMR histomorphometry nomenclature committee. J Bone Miner Res. 2013;28(1):2‐17. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 25. Miyauchi Y, Ninomiya K, Miyamoto H, et al. The Blimp1‐Bcl6 axis is critical to regulate osteoclast differentiation and bone homeostasis. J Exp Med. 2010;207(4):751‐762. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 26. Koenig A, Vaeth M, Xiao Y, et al. NFATc1/αA and Blimp‐1 support the follicular and effector phenotype of Tregs. Front Immunol. 2022;12:791100. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 27. Mun SH, Jastrzebski S, Kalinowski J, et al. Sexual dimorphism in differentiating osteoclast precursors demonstrates enhanced inflammatory pathway activation in female cells. J Bone Miner Res. 2021;36(6):1104‐1116. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 28. Zhang X, Peng D, Xi Y, et al. G9a‐mediated methylation of ERα links the PHF20/MOF histone acetyltransferase complex to hormonal gene expression. Nat Commun. 2016;7:10810. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 29. Haebe JR, Bergin CJ, Sandouka T, Benoit YD. Emerging role of G9a in cancer stemness and promises as a therapeutic target. Oncogenesis. 2021;10(11):76. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 30. Flesher JL, Fisher DE. G9a: an emerging epigenetic target for melanoma therapy. Epigenomes. 2021;5(4):23. [DOI] [PMC free article] [PubMed] [Google Scholar]
  • 31. Rahman Z, MohdR B, Devabattula G, MohdA K, Godugu C. Targeting H3K9 methyltransferase G9a and its related molecule GLP as a potential therapeutic strategy for cancer. J Biochem Mol Toxicol. 2021;35(3):e22674. [DOI] [PubMed] [Google Scholar]

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Figure S1.

FSB2-38-e23779-s001.pptx (1.2MB, pptx)

Data Availability Statement

The data that support the findings are available in the Methods section and Supplementary Data of this article.


Articles from The FASEB Journal are provided here courtesy of Wiley

RESOURCES