Abstract
Periodontitis, one of the most common non-communicable diseases, is characterized by chronic oral inflammation and uncontrolled tooth supporting alveolar bone resorption. Its underlying mechanism to initiate aberrant oral barrier immunity has yet to be delineated. Here, we report a unique fibroblast subpopulation activated to guide oral inflammation (AG fibroblasts) identified in a single-cell RNA sequencing gingival cell atlas constructed from the mouse periodontitis models. AG fibroblasts localized beneath the gingival epithelium and in the cervical periodontal ligament responded to the ligature placement and to the discrete application of Toll-like receptor stimulants to mouse maxillary tissue. The upregulated chemokines and ligands of AG fibroblasts linked to the putative receptors of neutrophils in the early stages of periodontitis. In the established chronic inflammation, neutrophils together with AG fibroblasts appeared to induce type 3 innate lymphoid cells (ILC3s) that were the primary source of interleukin-17 cytokines. The comparative analysis of Rag2-/- and Rag2γc-/-mice suggested that ILC3 contributed to the cervical alveolar bone resorption interfacing the gingival inflammation. We propose that AG fibroblasts function as a previously unrecognized surveillant to initiate gingival inflammation leading to periodontitis through the AG fibroblast-neutrophil-ILC3 axis.
Introduction
The oral mucosa is one of the most active barrier tissues in the human body [1, 2] and oral barrier immunity acts as a crucial surveillance system, which must defend from pathogenic organisms, while ensuring sustained tolerance of commensal microbes and harmless antigens [3]. However, once this gingival inflammation progresses, adjacent alveolar bone supporting the cervical area of the dentition is subjected to a localized, but often severe, osteoclastic bone resorption, resulting in tooth loss and disruption of the maxillofacial structure [4]. Periodontitis is not only the most frequent cause of tooth loss in adults [5], but also, globally, ranks among the most significant contributors to poor health and decreased quality of life, imposing substantial economic and healthcare burdens [6]. It is therefore acutely important to elucidate the molecular mechanisms underlying the initiation of gingival inflammation, in order to facilitate the development of effective preventive and therapeutic measures for this debilitating condition.
The postulated pathological framework has been reconstructed from animal models [7, 8], harvested immune cell analyses, and oral microbial associations [4, 9, 10]. For example, the ligature-induced periodontitis mouse model has shown that fully developed periodontal inflammation is closely associated with the abundant recruitment and excessive activation of neutrophils [11] and pathological induction of interleukin (IL)-17–secreting proinflammatory effector CD4+ T helper (Th)17 cells [12]. Based on our current understanding of osteoimmunology, it was hypothesized that these pathological immune cells activate osteoclasts, leading to induction of periodontal alveolar bone resorption [13, 14]. However, the precise cascade of events leading to this disease outcome has not been fully deciphered.
Investigation of the molecular mechanisms underlying interaction between different inflammatory cell fractions in immunological disorders is an intense area of research. Prior studies have shown that chemokines are involved in the trafficking and activation of inflammatory cells during both homeostasis and disease-associated inflammation [14]. These small proteins, which typically contain disulfide cysteine–cysteine (CC) and cysteine–X–cysteine (CXC) molecular signatures, act as ligands for cellular receptors that modulate inflammatory signaling pathways [15]. Such chemokine–receptor networks direct the migration of inflammatory cells, potentially amplifying the resulting tissue damage. However, increasing evidence suggests that in addition to immune cells, barrier tissue stromal cells are involved in innate immune cell regulation [16–19]. We and others have found that oral fibroblasts secrete cytokines and chemokines in response to microbial stimuli or a proinflammatory environment, suggesting a possible mechanism by which they may regulate immune cell function [8, 20]. Therefore, a comprehensive elucidation of chemokine–receptor signaling networks will likely need to include all types of gingival cells.
In the present study, to better understand the pathways contributing to chronic oral inflammation, we constructed the gingival single-cell transcriptomic atlas of the ligature-induced periodontitis mouse model as well as the newly developed mouse model of maxillary topical application (MTA) of pathological oral microbial biofilm component(s), discretely mimicking the initial events of gingival inflammation and degradation. Here we report a novel subpopulation of fibroblasts activated to guide chronic inflammation (AG fibroblasts). Our findings suggest that chemokines and ligands derived from AG fibroblasts bind to and activate receptors involved in neutrophil and lymphocytes. We further show that in the fully developed chronic gingival inflammation of the ligature-induced periodontitis model, Cd4- innate lymphoid cells (ILCs) predominantly produce proinflammatory IL-17 cytokines and provide evidence that type 3 ILCs (ILC3s) are responsible for cervical alveolar bone resorption in mice—a pathological feature consistent with the human periodontitis disease phenotype. Thus, the results from this study identify the AG fibroblast–neutrophil–ILC3 axis as a previously unrecognized mechanism contributing to pathological inflammation in periodontitis.
Results
Alterations in major cell type proportions during periodontitis development
The alveolar bone and gingival tissue surrounding the maxillary second molar of C57BL/6J mice presented healthy connective tissue supporting the dentition with minimal CD45+ immune cell infiltration and well-developed collagen architecture (Figure 1A – figure supplement 1A). The placement of a ligature around the maxillary left second molar[21] induced a localized and small gingival connective tissue degradation and infiltration of CD45+ immune cells near the cervical area (Figure 1B). The noticeable gingival defect and more prominent though still localized CD45+ immune cell infiltration developed after Day 3 following ligature placement (Figure 1C – figure supplement 1A, 1B). On Day 7, the gingival connective tissue was largely degraded, and chronic inflammation was developed (Figure 1D – figure supplement 1A, 1B). Three-dimensional (3D) reconstruction of micro-computed tomography (microCT) images (figure supplement 1C) further revealed a reduction in alveolar bone height starting from Day 3, which progressively increased on Day 7 (figure supplement 1D, 1E). Overall, the observed pattern of periodontal tissue degradation was consistent with that reported in previous studies [22, 23].
Left-side palatal gingiva tissue was harvested from mice on Day 0 (i.e., healthy gingiva without ligature placement) and on Days 1, 4, and 7 after ligature placement, and gingival cells were dissociated for single-cell RNA sequencing (scRNA-seq) [8, 24]. On Days 0 and 1, the major cell types identified included epithelial cells expressing cadherin 1 (Cdh1) [25] and type XVII collagen (Col17a1) [26], fibroblasts expressing type I collagen (Col1a1) [27] and lumican (Lum) [28], B cells expressing membrane spanning 4 domains A1 (Ms4a1) [29] and cluster of differentiation 79A (Cd79a) [30], T cells expressing epsilon subunit of T cell receptor complex (Cd3e) [31] and cluster of differentiation 5 (Cd5) [32], and myeloid cells expressing lysozyme 2 (Lyz2) [33] and integrin subunit alpha M (Itgam) [34] (Figure 1E and F—figure supplement 2). On Days 4 and 7, an additional endothelial cell fraction expressing selectin P (Selp) [35] and selectin E (Sele) [36] emerged (Figure 1E and F—figure supplement 2), suggesting increased inflammatory neovascularization with the progression of periodontal inflammation. The proportion of B cells was increased on Day 1, and the proportion of myeloid cells increased progressively from Day 4 to 7 (Figure 1F). In addition, the proportion of fibroblasts was increased on Day 7 (Figure 1F).
Fibroblasts activated to guide leukocyte migration in periodontitis development
We previously identified two distinct subpopulations of gingival fibroblasts, differentiated by expression of type XIV collagen (Col14a1) [8], and these were also detected in our current scRNA-seq data from Day 0 to 7 (Figure 2A). Gene ontogeny (GO) enrichment analysis of the Col14a1-expressing fibroblast subpopulation revealed expression of major gene clusters related to immune regulation, including “Regulation of leukocyte migration” (Figure 2B). This immune regulatory phenotype appears to be unique to Col14a1-expressing fibroblasts, which are therefore referred to as ‘fibroblasts activated to guide leukocyte migration’ or AG fibroblasts. The other subpopulation of Col1a1-expressing fibroblasts appeared to keep typical fibroblastic features and is thus referred to as KT fibroblasts (Figure 2C). GO analysis suggested that KT fibroblasts primarily maintained connective tissue remodeling but did not express an immune regulatory phenotype. An additional fibroblast subpopulation expressing smooth muscle aortic actin 2 (Acta2) was detected on Day 7, and these were identified as myofibroblasts [37] (MF) (Figure 2A). We found that the proportions of KT and AG fibroblasts were equal on Day 0 and on Days 1 and 4 after ligature placement (Figure 2E, 2F). However, on Day 7, the proportion of AG fibroblasts decreased, and the MF fraction emerged (Figure 2E, 2F).
AG fibroblasts and immune surveillance in early periodontitis development
To further characterize AG fibroblasts and their “Regulation of leukocyte migration” phenotype, we analyzed expression of CC motif chemokine ligands (CCLs) and CXC motif chemokine ligands (CXCLs) within fibroblast subpopulations in our scRNA-seq dataset [13]. Results show that AG fibroblasts activated expression of Ccl8, Ccl11, Ccl19, Cxcl1, Cxcl10, and Cxcl12 immediately after ligature placement on Days 1 and 4 (Figure 3A). Chemokine expression levels then decreased on Day 7, suggesting that AG fibroblasts might initiate early leukocyte migration into gingival tissue. Given that activation of Toll-like receptors (TLR) is known to increase chemokine expression, we further assessed expression of TLRs in fibroblast subpopulations. We found that AG fibroblasts displayed temporal activation of Tlr2, Tlr3, and Tlr4 expression on Days 1 and 4 (Figure 3B). Similarly, the TLR downstream genes Myd88, Irak1, Map3k7, and RelA were also expressed on Days 1 and 4 (Figure 3B). These data suggest that AG fibroblasts may guide the establishment of an early inflammatory environment within the gingiva and thereby promote periodontitis pathogenesis.
To validate the presence of AG fibroblasts, Day 1 gingiva and periodontal tissue histological sections were subjected to immunohistochemistry with antibodies against COL14A1 and CXCL12 (Figure 3C). We detected COL14A1- and CXCL12-positive AG fibroblasts localized near gingival epithelial cells in the connective tissue papillae and free gingiva, as well as in the cervical zone of the periodontal ligament (PDL) space (Figure 3C). We note that this localization pattern of AG fibroblasts appears to be highly suitable for early immune surveillance during periodontitis pathogenesis. Thus, we have hypothesized that AG fibroblasts initially sense the pathological stress including oral microbial stimuli through TLRs and secrete inflammatory signals through chemokine expression.
AG fibroblasts induced by maxillary topical application (MTA) of unmethylated cytidine phosphate guanosine oligonucleotide (CpG ODN) and of Porphyromonas gingivalis lipopolysaccharide (LPS)
To further characterize the behavior of AG fibroblasts, we used the recently developed maxillary topical application (MTA) model without the placement of a ligature. The MTA model applied oral microbial biofilm directly to the maxillary gingiva and held under a custom-made oral appliance (Figure 4A) for 1 hr. Four days after the topical application, mouse maxillary tissue was harvested for histological analysis and for scRNA-seq (Figure 4B) [8]. We previously reported that oral microbial biofilm and its extracellular polymeric substances such as extracellular DNA, but not planktonic microbes, induced initial periodontal tissue reaction in the MTA model [8]. The present study discretely applied ligands of TLR9 and TLR2/4: unmethylated CpG ODN and P. gingivalis LPS, respectively; using the MTA model. CpG ODN induced the expression of cathepsin K (Ctsk) in the periodontal ligament (PDL) and gingival connective tissue (Figure 4C) validating the previous report [8]. Furthermore, there were signs of localized Ctsk+ osteoclastic activities on the surface of alveolar bone (Figure 4C).
The scRNA-seq gingival cell composition of the MTA of CpG ODN and LPS indicated the early periodontitis pattern (Figure 4E) consistent with Day 1 scRNA-seq of the ligature-induced periodontitis (Figure 1F). The MTA of CpG ODN and LPS identified AG fibroblasts along with KT fibroblasts, whereas the MTA of LPS induced the additional MF (Figure 4E). The further analysis demonstrated the expression of CCL and CXCL chemokines by AG fibroblasts activated by the MTA of CpG ODN and LPS (Figure 4F). CpG ODN and LPS are sensed by Toll-Like Receptor 9 (TLR9) and TLR2/4, respectively. The scRNA-seq of MTA of CpG ODN did not capture the transcriptional activation of Tlr9, while the downstream effector genes such as Myd99, Irak1 and RelA were found upregulated in AG fibroblasts (Figure 4G). The MTA of LPS appeared to upregulate Tlr2 and Tlr4 gene transcription of AG fibroblasts and of MF albeit at lesser degrees (Figure 4G). These results indicated that AG fibroblasts were activated by CpG ODN and LPS through the MTA model, while these microbial elements appeared to exert different reactions.
Myeloid cell composition and activity during periodontitis development
The infiltration of proinflammatory neutrophils into the gingiva has been extensively characterized in this mouse model of periodontitis [11, 38]. Here, we found that macrophages expressing cluster of differentiation 86 (Cd86) [39] and integrin subunit alpha X (Itgax) [40] were predominant on Day 0 in the myeloid cell fraction from healthy gingiva (Figure 5A, 5B—figure supplement 3). In contrast, the proportion of neutrophils expressing CXC motif chemokine receptor (Cxcr2) [41] and G0/G1 switch gene 2 (G0s2) [42] increased after ligature placement and during periodontitis development from Day 1 to 7 (Figure 5A, 5B—figure supplement 3), suggesting continuous neutrophil infiltration into the early and established gingival lesion. Moreover, after ligature placement, gingival neutrophils upregulated the expression of triggering receptor expressed on myeloid cells 1 (Trem1), indicating that these cells are activated and participating in the amplification of inflammatory signals [43, 44] (Figure 5C). Strikingly, the Trem1-expressing activated neutrophils also show upregulation of matrix metalloproteinase 9 (Mmp9; Figure 5C)— a protein associated with extracellular matrix degeneration within gingival tissue [45]. We further detected expression of tumor necrosis factor (Tnf) and transforming growth factor beta 1 (Tgfb1) in both macrophages and neutrophils on Day 0 and after ligature placement (Figure 5D). Collectively, these myeloid cell behaviors are consistent with those reported in prior studies on periodontitis development, thus validating our scRNA-seq data.
Myeloid cells are also known to stimulate other immune cells through the expression of CCL and CXCL chemokines, many of which are associated with periodontitis development [46]. Our scRNA-seq analysis revealed that macrophages expressed Ccl2, Ccl9, Cxcl4, and Cxcl6, and neutrophils expressed Ccl3, Ccl4, Ccl6, Cxcl2, and Cxcl3 throughout periodontitis development (Figure 5E, 5F). These data suggest that chemokines and cytokines produced by macrophages and neutrophils in inflamed tissue may amplify and polarize the immune response toward chronic gingival inflammation.
Role of AG fibroblasts in myeloid cell activation
The interaction between chemokine ligands and their receptors has been extensively studied. Here, to evaluate the interaction between chemokine ligands strongly expressed by AG fibroblasts and chemokine receptors in innate immune cells, we matched AG fibroblast-expressed chemokines to expression of their putative receptors in myeloid cells. Results show that AG fibroblasts may regulate macrophages via the expression Ccl8 and Ccl11, which encode chemokines that can interact with CC chemokine receptors (CCRs) encoded by Ccr2 and Ccr5 in macrophages. Similarly, gene expression signatures suggest that AG fibroblast-mediated neutrophil regulation may occur through CCL8–CCR1, CXCL1–CXCR2, and CXCL12–CXCR4 interactions (Figure 6A). Notably, all chemokine–receptor pairs are expressed throughout periodontitis development, although expression of factors mediating the interaction between AG fibroblasts and macrophages was decreased on Day 7.
Intercellular communication may occur via other ligand–receptor interactions, which induce downstream target gene expression in recipient cells. Here, we performed NicheNet analysis to identify potential interactions between gingival fibroblasts and myeloid cells, revealing a trend toward increasing ligand–target interactions between fibroblasts and myeloid cells from Day 1 to Day 7 (Figure 6B). On Day 1, ligand expression was present primarily in AG fibroblasts, whereas on Day 7, MFs also exhibited prominent ligand expression (Figure 6C). Analysis of target gene expression revealed that both macrophages and neutrophils expressed receptors capable of interacting with AG fibroblast ligands on Day 1. In contrast, on Day 7, receptor expression was almost exclusively present in neutrophils, suggesting these cells are primarily targeted by MFs and AG fibroblasts at later stages of periodontitis development (Figure 6D).
The chemokine-receptor interaction was also suggested between AG fibroblasts and myeloid cells in the MTA of CpG DON (Figure 6E). However, in the MTA of LPS, macrophage lacked the detectable expression of chemokine receptor and the AG fibroblastic chemokine interaction appeared to be limited to neutrophils (Figure 6E).
Expression of osteoclastogenic cytokines
Macrophage-colony stimulating factor (M-CSF), encoded by Csf1, and receptor activator of nuclear factor kappa-Β ligand (RANKL), encoded by Tnfsf11, are known to play critical roles in osteoclast induction [47]. We found that in the ligature-induced periodontitis model, the expression of M-CSF was detected in various cell types including fibroblasts (Figure 7A), whereas the RNAKL expression was more restricted to fibroblasts and T cells (Figure 7B). AG fibroblasts consistently expressed M-CSF from Day 0 to Day 7 and neutrophils started to express M-CSF in Day 4 and Day 7 (Figure 7C). The expression of RNAKL was detected in AG fibroblasts but not in myeloid cells (Figure 7D). T cells and type 3 innate lymphoid cells (ILC3; see below) were also shown to express M-CSF and RNAKL (Figure 7C, 7D).
The MTA of CpG ODN and LPS suggested the AG fibroblasts were one of the predominant cellular sources of M-CSF (Figure 7E). By contrast, the expression of RNAKL was only detected in AG fibroblasts of the MTA of CpG ODN and LPS did not seem to induce RANKL expression in the MTA model (Figure 7F). The osteoclastic activity found in the MTA of CpG ODN (Figure 4D) might be induced by M-CSF and RANKL derived from AG fibroblasts.
Helper T (Th) cells, cytotoxic T (Tc) cells, regulatory T (Treg) cells and innate lymphoid cells (ILCs) in periodontitis development
Our scRNA-seq dataset found that the Day 0 healthy gingiva exclusively contained Th cells expressing Cd4 [48]. However, after ligature placement, Tc cells expressing Cd8 [48] emerged, and on Day 7, Cd4-Zbtb16+ ILCs and Cd4+Foxp3+ Treg cells [49] were also detected (Figure 8A, 8B). The ILCs expressed Nfil3, a basic leucine zipper (bZIP) transcription factor required for ILC development [50] (Figure 8C). However, to our surprise, expression of Rorc, Il17a, and Il17f was detected in ILCs (Figure 8C), indicating that this gingival subpopulation is predominantly composed of type 3 ILCs (ILC3s). ILC3s and Th17 cells share similar regulatory functions. However, the role of ILC3s in the development of periodontitis has not been fully deciphered.
ILC3s are critical for cervical alveolar bone resorption in the mouse periodontitis model
We thus examined the role of ILC3s in periodontitis pathogenesis by measuring ligature-induced gingival defects and alveolar bone resorption in Rag2-/- mice, which lack functional B, Th, and Tc cells, and Rag2–IL-2 receptor common gamma (Il2rg) double-knockout mice (Rag2γc-/-), lacking all lymphocytes including ILCs. After ligature placement, we found that alveolar bone loss was decreased in Rag2-/- mice and nearly eliminated in Rag2γc-/- mice (Figure 8D). MicroCT image analysis indicated a better perseveration of alveolar bone structure in Rag2γc-/- mice, relative to the other groups (Figure 8E—figure supplement 4A, 4B). However, gingival defects developed similarly in wild-type (WT), Rag2-/-, and Rag2γc-/- mice (figure supplement 4C, 4D).
Histologically, osteoclastic resorption lacunae were observed on the alveolar bone surface at the cervical PDL and tooth apex PDL zones in WT mice (Figure 8F, 8G). In addition, we detected a significant decrease in the number of tartrate-resistant acid phosphatase (TRAP)-positive osteoclasts in the cervical PDL zone of Rag2γc-/- mice (Figure 8H) and in the apical PDL zone of both Rag2-/- and Rag2γc-/- mice (Figure 8I). Collectively, these data suggest that ILC3s, not Th17 cells, are responsible for cervical alveolar bone resorption in the mouse periodontitis model, a pathological phenotype consistent with human disease (figure supplement 4E, 4F).
The role of AG fibroblasts and neutrophils in ILC3 development in periodontitis
Based on our present data, we hypothesize that ILC3s within the gingival tissue play a pathological role in cervical alveolar bone resorption in our mouse model. We therefore aimed to identify the cells that promote ILC3 development in mice. Similar to Th17 cells, ILC3 development was shown to be triggered by IL-6 and IL-23a [51, 52]. A survey of our scRNA-seq data identified epithelial cells, fibroblasts, and myeloid cells as the source of Il6 in gingival tissue (Figure 9A). On Days 1 and 4 following ligature placement, AG fibroblasts primarily expressed Il6, whereas neutrophils became the predominant source on Day 7. We further found that gingival epithelial cells comprise at least four different subpopulations, plus an additional group displaying an epithelial–mesenchymal transition (EMT) phenotype on Day 7; Il6 was expressed by several of these epithelial subpopulations, including the EMT subgroup (Figure 9A—figure supplement 5A–D). Il23a was also detected AG fibroblasts and myeloid cells (Figure 9B), with expression present in all epithelial cell subsets at various points in periodontitis development (Figure 9B— figure supplement 5D).
Lastly, we evaluated potential interactions between chemokine ligands expressed by AG fibroblasts and neutrophils and chemokine receptors in innate immune cells. Our data suggest the presence of a chemokine–receptor association between ILC3s and both AG fibroblasts (Figure 9C) and neutrophils (Figure 9D), although interactions with other innate immune cells are also possible. NicheNet analysis further identified potential ligand–target gene associations between lymphocytes, including ILC3s, and both fibroblasts (Figure 9E) and myeloid cells (Figure 9F). Ligand expression was more prominent in AG fibroblasts than in other fibroblast subpopulations (Figure 9G) and elevated in neutrophils relative to macrophages (Figure 9H). Additionally, target gene expression was detected in ILC3s, although it was not specific to these cells (Figure 9I, 9J). Thus, in total, our data suggest a regulatory role for a newly identified AG fibroblast subpopulation in the gingiva, which appears to orchestrate chronic gingival inflammation, at least in the early stages, and to promote alveolar bone resorption via stimulation of neutrophils and ILC3s.
Discussion
The present study revealed that the gene signature of a unique and previously uncharacterized subpopulation of gingival fibroblasts, referred to as AG fibroblasts, could possess the functional capability to serve as an oral immune surveillant and to orchestrate the initiate gingival inflammation. Oral barrier immunity presents a complex interaction between different types of immune cells to protect and maintain the oral environment and structure. Given the frequent exposure to a diverse commensal and pathological microbiome, physical and chemical insults, and dietary and airborne antigens, the oral barrier immune mechanism must resiliently establish the highly tolerant homeostasis [53]. This oral immune homeostasis mechanism is not yet fully understood to date; however, the aberrant oral immune response and gingival chronic inflammation leading to periodontitis have provided an important clue to elucidate the oral barrier immunity.
The previous studies suggested gingival intraepithelial γ8T cells [54], a subset of neutrophils [55], macrophages [56], and dendritic cells [57] should serve the candidate immune surveillant in the oral barrier tissue. These immune cells recognize the pathological signals through the TLR sensing mechanism. TLR2-/-, TLR4-/-, TLR2&4-/- [58] and TLR9-/- [59, 60] mice exhibited the reduced alveolar bone loss in various mouse periodontitis models, suggesting that TLRs indeed played an important role in developing discordant oral barrier immunity. However, these studies did not identify the TLR-carrying surveillant cells centrally involved in initiating the chronic gingival inflammation.
It has been suggested that TLR-expressing gingival fibroblasts may regulate innate immune responses [61, 62]. In the present study, we found that not all gingival fibroblasts acquired an immune-sensing capability, but rather, AG fibroblasts represented a distinct fibroblast subpopulation, capable of responding to microbial and tissue damage signals to serve initiate immune surveillance. To test if TLR ligands stimulate the AG fibroblast activation, we applied a newly developed mouse system employing the discrete MTA model [8]. The present study used TLR9 ligand: CpG ODN; and TLR2/4 ligand: LPS. Topical application of these TLR ligands to the maxillary gingival tissue activated AG fibroblasts and increased expression of CC and CXC chemokines (Figure 4F). It was noted that the CpG ODN application did not upregulate Tlr9 expression, while the LPS application increased Tlr2/4 expression. TLR9 mRNA level has not been fully correlated to the chronic inflammatory diseases; however, the pathological activation of TLR9 caused the discordant downstream inflammation [63]. In the present study, the differential expression of the TLR downstream signaling molecules suggested the ligand specific response by AG fibroblasts. In fact, both stimulants increased the expression of myeloid differentiation marker 88 (Myd88), IL1R-associated kinase (Irak), mitogen-activated protein kinase (Map3k7) and NF-kB subunit RelA (Figure 4G), suggesting that the corresponding TLRs were indeed stimulated [64, 65]. The Tlr expression pattern of Day 1 AG fibroblasts was similar to that of the MTA of CpG ODN, whereas the Day 4 AG fibroblasts resembled the Tlr expression pattern of the MTA of LPS (Figure 3B). It is tempting to speculate that microbial DNA/TLR9 activation of AG fibroblasts may initiate the early pathological process followed by the LPS/TRL2/4 stimulation in the more established periodontal inflammation. Taken together, we hypothesize that AG fibroblasts sensitively detect the microbial signals and play a previously unrecognized role in innate immune regulation during the early periodontitis development.
One of the major observations in the present study was that AG fibroblasts clearly upregulated the chemokine expression only after the pathological stimulation. CCL and CXCL chemokines orchestrate the chronic inflammation through the ligation and activation of their putative chemokine receptors in a specific mechanism. For example, neutrophils express a relatively limited number of chemokine receptors—CXCR2, CXCR4, and CCR1 [66, 67]. Therefore, neutrophil migration and trafficking to the affected gingiva as well as its pathological transformation in chronic periodontitis patients [11, 68, 69] may be guided through the specific chemokine-receptor interaction. Consistent with prior studies, we detected the expression of these chemokine receptors in gingival neutrophils, and they were linked to the chemokine ligands expressed by AG fibroblasts—CXCL1, CXCL12, and CCL8, respectively (Figure 6). In addition, neutrophils can respond to various chemo-attractants that modulate and finetune various cellular behaviors, such as migration direction, adhesion strength, and functional heterogeneity [70]. Here, we observed evidence for substantial intracellular communication between AG fibroblasts and neutrophils through ligand–target gene interactions (Figure 6). Although validation of each putative ligand–receptor interaction is beyond the scope of this study, our data suggest that AG fibroblasts support neutrophil trafficking in both early and established periodontitis lesions.
A primary pathological consequence of periodontitis is uncontrolled alveolar bone resorption, leading to tooth loss, which is thought to be mediated by osteoimmune pathways involving IL-17– expressing Th17 cells [71]. Surprisingly, analysis of our scRNA-seq revealed that the expression of Rorc, Il17a, and Il17f was not detected in CD4+ Th cells, but rather, in lymphocytes with ILC characteristics (Figure 8). ILCs share phenotypic and functional features with CD4+ T cells, although they lack antigen-specific T cell receptors [50]. The gingival ILCs expressed Rorc, Il17a, and Il17f, but not Tbx21 and Gata3, indicating that they were type 3 ILCs (ILC3s). Notably, recent clinical studies have reported the presence of ILCs in gingiva from human periodontitis patients [72, 73] and leptin receptor–deficient mice [74], although their pathological contributions have not been fully elucidated.
In the present study, we further explored the role of ILC3s in periodontitis-associated bone loss using Rag2-/- and Rag2γc-/- mice. Rag2-/- mutation prevents V(D)J recombination required for generating immunoglobulin and T cell receptors, resulting in the production of functionally immature B and T cells, including Th17 cells. However, ILCs do not undergo genomic receptor rearrangements and, thus, are unaffected by Rag2-/- mutation. In contrast, Rag2γc-/- mice have the additional Il2rg-/- mutation, which disables common γ chain cytokines (γc). Therefore, in addition to non-functional B and T cells, these animals also have defective γc-dependent ILCs. The ligature-placed Rag2-/- mice exhibited a reduced number of osteoclasts in the apical PDL area, indicative of reduced bone resorption in this root apical region. By contrast, Rag2γc-/- mice demonstrated a significant loss of bone resorption in both the cervical and apical PDL areas (Figure 8). Human periodontitis induces the localized alveolar bone resorption at the cervical PDL zone interfacing the gingival inflammatory legion. In contrast, apical alveolar bone resorption is observed in clinical cases of root canal infection (apical periodontitis). Therefore, our findings suggest that ILC3s, which are differentially impacted by Rag2γc-/- and Rag2-/- mutations, may be primarily responsible for human periodontitis-like cervical alveolar bone resorption near the site of gingival inflammation.
Differentiation of ILC3s and Th17 cells is mediated by similar environmental cues, with IL-6 and IL-23a playing critical roles in both pathways [75, 76]. Our present scRNA-seq data (Figure 9) suggest that AG fibroblasts and neutrophils are the primary cellular sources of Il6 during the early and later stages, respectively, of periodontitis development. In contrast, Il23a was found to be expressed by a variety of cell types, such as AG fibroblasts, myeloid cells, and multiple subsets of epithelial cells, including those with an EMT phenotype [77]. Further, CC and CXC chemokine– receptor associations between ILC3s and both AG fibroblasts and neutrophils appeared to be non-specific. Therefore, our data suggest that rather than a specific trigger, the collective gingival environment, which includes AG fibroblasts, might contribute to ILC3 differentiation.
In conclusion, based on our present findings, we propose that a previously unrecognized AG fibroblast subpopulation in the gingiva can facilitate immune surveillance and participate in the pathological regulation of innate immune cells, such as proinflammatory neutrophils, within oral barrier tissue (Figure 10). Moreover, we hypothesize that ILC3s in the inflamed gingiva play a critical role in pathological alveolar bone resorption in the mouse model of periodontitis and, potentially, in human disease. Thus, the newly proposed AG fibroblast–neutrophil–ILC3 axis may hold valuable clues for unraveling the pathological mechanisms underlying periodontitis development. Moreover, these findings also provide a basis for investigation of new preventive and therapeutic strategies to contain oral barrier inflammation and potentially sever the link between periodontitis and debilitating non-communicative metabolic and cardiovascular diseases.
Methods
Animal care
All protocols for animal experiments were reviewed and approved by the University of California Los Angeles (UCLA) Animal Research Committee (ARC# 2003-009) and followed the Public Health Service Policy for the Humane Care and Use of Laboratory Animals and the UCLA Animal Care and Use Training Manual guidelines. C57BL/6J WT, Rag2-/-, and Rag2γc-/- mice were purchased from the Jackson Laboratory (Bar Harbor, ME, USA). Animals had free access to regular rodent diet and water and were maintained in standard housing conditions with 12-h light/dark cycles in the Division of Laboratory Animal Medicine at UCLA.
Evaluation of gingival defect and alveolar bone resorption in a ligature-induced mouse model of periodontitis
A silk thread was gently tied around the left maxillary second molar of 8- to 12-week-old female WT, Rag2-/-, and Rag2γc-/- mice under general inhalation anesthesia with isoflurane (Henry Schein, Melville, NY, USA). On Days 1, 3, and 7 after ligature placement, WT mice were randomly chosen and euthanized by 100% CO2 inhalation. WT mice without ligature placement were used as Day 0 pre-periodontitis control. Rag2-/- and Rag2γc-/- mice were euthanized on Day 7. The palatal tissue was digitally photographed, and maxillae were harvested and fixed in 10% buffered formalin (Thermo Fisher Scientific, Waltham, MA, USA). The gingival defect area was measured from digital photographs, using the ImageJ Java-based image-processing program (NIH, Bethesda, MD, USA), and normalized to the circumferential area of the maxillary first molar. Fixed maxillae were subjected to microCT imaging at an energy level of 60 kV and 166 µA, and 3D images were reconstructed from microCT scans (Skyscan 1275: Bruker, Billerica, MA, USA). Alveolar bone loss was assessed at three sites (mesiobuccal cusp, distobuccal cusp, and distal cusp) of the first molar, two sites (mesiobuccal cusp and distobuccal cusp) of the second molar, and one site (buccal cusp) of the third molar by measuring the distance from the cementoenamel junction (CEJ) to the alveolar bone crest (ABC) on the buccal and palatal side of the alveolar bone. Total bone loss was calculated from the six-site total CEJ–ABC distance. The bone volume/total volume (BV/TV) ratio, bone surface, trabecular number (Tb.N), and trabecular thickness (Tb.Th) in the buccal side of alveolar bone of the second molar were determined using the proprietary analysis program (CTan: Bruker, Billerica, MA, USA). Statistical analysis to assess differences among multiple experimental groups was performed using one-way analysis of variance (ANOVA) with Tukey’s multiple-comparison test; p < 0.05 was considered to be statistically significant.
Evaluation of gingival effect in the maxillary topical application (MTA) model
We have developed a method to apply chemical therapeutic agents topically to the mouse maxillary tissue [24, 78]. We used the maxillary topical application (MTA) model to apply oral microbial components [8]. The present study used the MTA model using 1 µg/ml unmethylated CpG oligonucleotide (CpG ODN: InvivoGen, San Diego, CA) or P. gingivalis lipopolysaccharide (LPS: InvivoGen). First, a custom-made oral appliance was fabricated using clear dental resin, covering the maxillary/palatal tissue between the molars. Mice were anesthetized and placed on a supine position. CpG ODN (3 µl) or LPS (3 µl) was pipetted over the maxillary tissue and the oral appliance was placed to hold the solution with a bite block for 1 hr in an anesthetization chamber. Mice were then transferred to the operation table and the oral appliance and bite block were removed. In general, there was no remaining solution in the mouth. Mice were returned to the cage in the vivarium. On Day 4 of the MTA, mice were euthanized by 100% CO2 inhalation and the maxillary tissue was harvested and fixed with 10% buffered formalin for histological evaluation.
Histological analysis
Fixed maxillae were decalcified in 10% EDTA (Sigma–Aldrich) for 3 weeks and then embedded in paraffin. Histological cross-sections were stained with hematoxylin and eosin (HE) and evaluated on a light microscope. Adjacent paraffin sections (4 µm) were subjected to a heat-induced epitope retrieval procedure and then immunohistochemically stained with polyclonal antibodies to CD45 (#PA5-11671, Thermo Fisher Scientific), COL14A1 (#PA5-49916, Thermo Fisher Scientific), CXCL12 (#PA5-30603, Thermo Fisher Scientific), or Ctsk (PA5-14270, Thermo Fisher Scientific) at the suggested dilution, followed by secondary antibody application, diaminobenzidine staining, and methylene blue counterstaining.
Using maxillary cross-sections of WT, Rag2-/-, and Rag2-/-γc-/- mice, osteoclasts were evaluated by TRAP staining using a commercially available kit (Acid Phosphatase TRAP kit, Sigma–Aldrich), according to the manufacturer’s instructions. TRAP-positive cells were counted in the 0.01-mm2 area. Statistical analysis to assess the differences among multiple experimental groups was performed using one-way ANOVA with Tukey’s multiple-comparison test; p < 0.05 was considered to be statistically significant.
Single cell dissociation from mouse maxillary gingiva
On Days 1, 4, and 7 after ligature placement, and on Day 4 of the MTA of CpG ODN or LPS, mice were euthanized by 100% CO2 inhalation. And maxillary gingival tissues (n = 4 per group) were harvested from freshly isolated mouse maxillae.
Collagenase II treatment
The tissues were cut into 1-mm pieces and placed immediately into digestion buffer, containing 1-mg/ml collagenase II (Life Technologies, Thermo Fisher Scientific, Waltham, MA, USA), 10-units/ml DNase I (Sigma–Aldrich, St. Louis, MO, USA), and 1% bovine serum albumin (BSA; Sigma–Aldrich) in Dulbecco’s Modified Eagle Medium (DMEM; Life Technologies). The tissues were incubated in digestion buffer for 20 min at 37 °C on a shaker at 150 rpm and then passed through a 70-µm cell strainer. The collected cells were pelleted at 1,500 rpm for 10 min at 4 °C and resuspended in phosphate-buffered saline (PBS; Life Technologies), supplemented with 0.04% BSA (Cell suspension A).
Trypsin treatment
Immediately following collagenase II treatment, tissues were incubated in 0.25% trypsin (Life Technologies) and 10-units/ml DNase I for 30 min at 37 °C on a 150-rpm shaker. Trypsin was neutralized with fetal bovine serum (Life Technologies), and the tissues were passed through a 70-µm cell strainer and washed with DMEM. The collected cells were then pelleted at 1,500 rpm for 10 min at 4 °C and resuspended in PBS with 0.04% BSA (Cell suspension B). Cell suspension A and Cell suspension B were combined in one tube. An equal number of combined cell suspensions A and B from four animals per group were combined for scRNA-seq analysis (10X Genomics, San Francisco, CA, USA).
Cell clustering and Identification
Cell Ranger was used to align reads, generate feature–barcode matrices, and perform clustering and gene expression analyses on the scRNA-seq data, and the output from this program was analyzed using the R-program Seurat (https://satijalab.org/seurat/). Cells with <2,400 genes detected or >1% mitochondrial gene expression were filtered out as low-quality cells. Individual gene counts for each cell were divided by the total gene counts for that cell and multiplied by a scale factor of 10,000; natural-log transformation was then applied to the counts. The FindVariableFeatures function was used to select 2,000 variable genes with default parameters, and the ScaleData function was used to scale and center the counts in the dataset. Principal component analysis and Uniform Manifold Approximation and Projection dimensional reduction were performed on variably expressed genes. The cluster markers were found using the FindAllMarkers function, and cell types were manually annotated based on the cluster markers. Cell types were assigned based on expression of cell marker genes, and gene expression within different cell types was displayed using dot plots and violin plots.
Functional annotation and pathway enrichment analysis
Annotation and visualization of GO terms were performed by Metascape (http://metascape.org/gp/index.html#/main/step1). The top-100 differentially expressed genes in each population were input and filtered with the term “immune”. Filtered genes were then input, and only “biological process” gene sets were retrieved from the GO database.
Ligand–target matrix prediction
NicheNet (v.1.0.0, https://github.com/saeyslab/nichenetr) was used to predict interactions between cell types. In brief, the integrated Seurat object containing each cell subpopulation was input into the NicheNet Seurat wrapper. Sender cells and receiver cells were determined, and interactions between active ligands expressed by sender cells and target receptors expressed by receiver cells were predicted based on information in signaling and ligand–receptor databases.
Acknowledgements
We thank Dr. Connie Lee of Division of Periodontology at UCLA School of Dentistry for her guidance on the clinical manifestations of human periodontitis. We also thank Dr. Yunfeng Li of the Translational Pathology Core Laboratory, Department of Pathology and Laboratory Medicine, David Geffen School of Medicine at UCLA for her immunohistology work. This study was supported by NIH grants R01DE022550, R44DE025524, and C06RR014529 and by SINTX Technologies.
Competing interests
I.N. is a consultant for FUJI FILM Corp and BioVinc LLC and received a research fund from SINTX Technologies, Inc. I.N. and A.H. received a research fund from Maruho Co. Ltd. The other authors declare no competing interests.
Data availability statement
We specify where the data supporting the results of this study can be found and provide hyperlinks to publicly archived datasets analyzed or generated during the study, where applicable. scRNA-seq data are found in NIH Gene Expression Omnibus (GSE228635).
Supplemental Information
References
- 1.The subgingival microbiome in health and periodontitis and its relationship with community biomass and inflammationISME J 7:1016–25
- 2.Homeostatic Immunity and the MicrobiotaImmunity 46:562–576
- 3.The oral mucosa: A barrier site participating in tissue-specific and systemic immunityOral Dis 24:22–25
- 4.The oral microbiota: dynamic communities and host interactionsNat Rev Microbiol 16:745–759
- 5.Predictors for tooth loss in periodontitis patients: Systematic review and meta-analysisJ Clin Periodontol 46:699–712
- 6.Oral diseases: a global public health challengeLancet 394:249–260
- 7.Application of Ligature-Induced Periodontitis in Mice to Explore the Molecular Mechanism of Periodontal DiseaseInt J Mol Sci 22
- 8.Oral microbial extracellular DNA initiates periodontitis through gingival degradation by fibroblast-derived cathepsin K in miceCommun Biol 5
- 9.Local and systemic mechanisms linking periodontal disease and inflammatory comorbiditiesNat Rev Immunol 21:426–440
- 10.Effect of Calcium Ion Supplementation on Oral Microbial Composition and Biofilm Formation In VitroMicroorganisms 10
- 11.New developments in neutrophil biology and periodontitisPeriodontol 82:78–92
- 12.Calcium Signaling in T Cells and Chronic Inflammatory Disorders of the Oral CavityJ Dent Res 100:693–699
- 13.Inflammatory and immune pathways in the pathogenesis of periodontal diseasePeriodontol 64:57–80
- 14.The chemokine system in innate immunityCold Spring Harb Perspect Biol 7
- 15.A guide to chemokines and their receptorsFEBS J 285:2944–2971
- 16.IL-33 and the intestine: The good, the bad, and the inflammatoryCytokine 100:1–10
- 17.Intestinal fibroblast-derived IL-10 increases survival of mucosal T cells by inhibiting growth factor deprivation- and Fas-mediated apoptosisJ Immunol 175:2000–9
- 18.PD-1 ligand expression by human colonic myofibroblasts/fibroblasts regulates CD4+ T-cell activityGastroenterology 135:1228–1237
- 19.CCL19-producing fibroblastic stromal cells restrain lung carcinoma growth by promoting local antitumor T-cell responsesJ Allergy Clin Immunol 142:1257–1271
- 20.Human oral mucosa cell atlas reveals a stromal-neutrophil axis regulating tissue immunityCell 184:4090–4104
- 21.Optimization of the ligature-induced periodontitis model in miceJ Immunol Methods 394:49–54
- 22.Analysis of Experimental Ligature-Induced Periodontitis Model in MiceMethods Mol Biol 2210:237–250
- 23.An experimental murine model to study periodontitisNat Protoc 13:2247–2267
- 24.Mechanism of bisphosphonate-related osteonecrosis of the jaw (BRONJ) revealed by targeted removal of legacy bisphosphonate from jawbone using competing inert hydroxymethylene diphosphonateElife 11
- 25.Oral Mucosal Epithelial CellsFront Immunol 10
- 26.High Expression of Collagen XVII Compensates for its Depletion Induced by Pemphigoid IgG in the Oral MucosaJ Invest Dermatol 138:1707–1715
- 27.Autocrine regulation of mesenchymal progenitor cell fates orchestrates tooth eruptionProc Natl Acad Sci U S A 116:575–580
- 28.Transcript profiling of periodontal fibroblasts and osteoblastsJ Periodontol 76:1044–55
- 29.Characterization of inflammatory cell infiltrate in human dental pulpitisInt Endod J 43:1013–21
- 30.Investigation of immune cell markers in feline oral squamous cell carcinomaVet Immunol Immunopathol 202:52–62
- 31.Cell Biology of T Cell Receptor Expression and RegulationAnnu Rev Immunol 36:103–125
- 32.Immunoregulatory function of CD3+, CD4-, and CD8-T cells. Gamma delta T cell receptor-positive T cells from nude mice abrogate oral toleranceJ Immunol 143:3415–22
- 33.Mouse lysozyme M gene: isolation, characterization, and expression studiesProc Natl Acad Sci U S A 85:6232–6
- 34.M2-polarized macrophages contribute to neovasculogenesis, leading to relapse of oral cancer following radiationSci Rep 6
- 35.Expression of P-selectin on endothelial cells is upregulated by LPS and TNF-alpha in vivoCell Adhes Commun 2:7–14
- 36.E-selectin mediates Porphyromonas gingivalis adherence to human endothelial cellsInfect Immun 80:2570–6
- 37.Neuronal PAS Domain 2 (Npas2)-Deficient Fibroblasts Accelerate Skin Wound Healing and Dermal Collagen ReconstructionAnat Rec (Hoboken 303:1630–1641
- 38.Primary immunodeficiencies reveal the essential role of tissue neutrophils in periodontitisImmunol Rev 287:226–235
- 39.Macrophage depletion abates Porphyromonas gingivalis-induced alveolar bone resorption in miceJ Immunol 193:2349–62
- 40.Prognostic Relevance of Macrophage Phenotypes in High-grade Oral Tongue Squamous Cell CarcinomasAppl Immunohistochem Mol Morphol 29:359–365
- 41.Loss of Neutrophil Homing to the Periodontal Tissues Modulates the Composition and Disease Potential of the Oral MicrobiotaInfect Immun 89
- 42.Activation of G0S2 is coordinated by recruitment of PML/RARalpha and C/EBPepsilon to its promoter during ATRA-induced APL differentiationJ Leukoc Biol 101:655–664
- 43.Cutting edge: inflammatory responses can be triggered by TREM-1, a novel receptor expressed on neutrophils and monocytesJ Immunol 164:4991–5
- 44.Critical limb ischaemia is characterised by an increased production of whole blood reactive oxygen species and expression of TREM-1 on neutrophilsAtherosclerosis 229:396–403
- 45.Immunolocalization of matrix metalloproteinases-2 and -9 during apical periodontitis developmentArch Oral Biol 54:764–71
- 46.Relationship between chemokines and dendritic cells in human chronic periodontitisJ Periodontol 85:1416–23
- 47.Bone resorption by osteoclastsScience 289:1504–8
- 48.Memory T cell subsets in healthy gingiva and periodontitis tissuesJ Periodontol 89:1121–1130
- 49.An update of knowledge on the regulatory role of Treg cells in apical periodontitisOral Dis 27:1356–1365
- 50.Innate lymphoid cells. Innate lymphoid cells: a new paradigm in immunologyScience 348
- 51.Skin-resident innate lymphoid cells converge on a pathogenic effector stateNature 592:128–132
- 52.Interleukin 6 Increases Production of Cytokines by Colonic Innate Lymphoid Cells in Mice and Patients With Chronic Intestinal InflammationGastroenterology 149:456–67
- 53.Tissue-Specific Immunity at the Oral Mucosal BarrierTrends Immunol 39:276–287
- 54.gammadeltaT cells in oral tissue immune surveillance and pathologyFront Immunol 13
- 55.Distinct Oral Neutrophil Subsets Define Health and Periodontal Disease StatesJ Dent Res 95:931–8
- 56.Innate Phagocyte Polarization in the Oral CavityFront Immunol 12
- 57.Dendritic cells of the oral mucosaMucosal Immunol 7:27–37
- 58.Different engagement of TLR2 and TLR4 in Porphyromonas gingivalis vs. ligature-induced periodontal bone lossBraz Oral Res 31
- 59.Toll-Like Receptor 9-Mediated Inflammation Triggers Alveolar Bone Loss in Experimental Murine PeriodontitisInfect Immun 83:2992–3002
- 60.Interplay of Toll-Like Receptor 9, Myeloid Cells, and Deubiquitinase A20 in Periodontal InflammationInfect Immun 85:e00814–16
- 61.Single-Cell RNA Sequencing Identifies New Inflammation-Promoting Cell Subsets in Asian Patients With Chronic PeriodontitisFront Immunol 12
- 62.Biological Roles of Fibroblasts in Periodontal DiseasesCells 11
- 63.Importance of TLR9-IL23-IL17 axis in inflammatory bowel disease development: Gene expression profiling studyClin Immunol 197:86–95
- 64.Cutting edge: Role of Toll-like receptor 9 in CpG DNA-induced activation of human cellsJ Immunol 167:3555–8
- 65.TLR signaling pathwaysSemin Immunol 16:3–9
- 66.Regulation of human neutrophil chemokine receptor expression and function by activation of Toll-like receptors 2 and 4Immunology 115:90–8
- 67.Solitary pulmonary nodules: evaluation with a CT reference phantomRadiology 170:653–6
- 68.Neutrophil recruitment and function in health and inflammationNat Rev Immunol 13:159–75
- 69.Neutrophil-derived chemokines on the road to immunitySemin Immunol 28:119–28
- 70.Neutrophil chemoattractant receptors in health and disease: double-edged swordsCell Mol Immunol 17:433–450
- 71.The Th17/Treg cell balance: crosstalk among the immune system, bone and microbes in periodontitisJ Periodontal Res 57:246–255
- 72.Enrichment of Innate Lymphoid Cell Populations in Gingival TissueJ Dent Res 97:1399–1405
- 73.Characterization of the human immune cell network at the gingival barrierMucosal Immunol 9:1163–1172
- 74.Single-Cell Transcriptomic Atlas of Gingival Mucosa in Type 2 DiabetesJ Dent Res 101:1654–1664
- 75.A T-bet gradient controls the fate and function of CCR6-RORgammat+ innate lymphoid cellsNature 494:261–5
- 76.Lineage relationship analysis of RORgammat+ innate lymphoid cellsScience 330:665–9
- 77.Epithelial-mesenchymal transition in gingival tissues from chronic periodontitis patients: A case-control studyDent Med Probl 58:311–319
- 78.Fluorescent risedronate analogue 800CW-pRIS improves tooth extraction-associated abnormal wound healing in zoledronate-treated miceCommun Med (Lond 2
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Copyright
© 2023, Kondo et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 1,359
- downloads
- 179
- citations
- 5
Views, downloads and citations are aggregated across all versions of this paper published by eLife.