A novel gene ZNF862 causes hereditary gingival fibromatosis

Abstract
Editor's evaluation
Introduction
Results and discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Hereditary gingival fibromatosis (HGF) is the most common genetic form of gingival fibromatosis which is featured as a localized or generalized overgrowth of gingivae. Currently two genes (SOS1 and REST), as well as four loci (2p22.1, 2p23.3–p22.3, 5q13–q22, and 11p15), have been identified as associated with HGF in a dominant inheritance pattern. Here, we report 13 individuals with autosomal-dominant HGF from a four-generation Chinese family. Whole-exome sequencing followed by further genetic co-segregation analysis was performed for the family members across three generations. A novel heterozygous missense mutation (c.2812G > A) in zinc finger protein 862 gene (ZNF862) was identified, and it is absent among the population as per the Genome Aggregation Database. The functional study supports a biological role of ZNF862 for increasing the profibrotic factors particularly COL1A1 synthesis and hence resulting in HGF. Here, for the first time we identify the physiological role of ZNF862 for the association with the HGF.

Editor's evaluation

This study is of clinical relevance to those interested in the etiology and pathology of hereditary gingival fibromatosis (HGF). The paper discusses two novel findings: identification of a causative role of a missense mutation in the gene encoding the zinc finger protein 862 (ZNF862) that leads to hereditary gingival fibromatosis (HGF), a rare disease characterized by overgrowth of gingivae, in an examined family, and a suggestion of the molecular consequences of that mutation that leads to the disease.

https://doi.org/10.7554/eLife.66646.sa0

Introduction

Gingival fibromatosis (GF) is a rare disorder which is characterized by a benign, nonhemorrhagic, localized, or generalized fibrous enlargement of free and attached gingivae with slow progression. GF could be complicated by epilepsy, hypertrichosis, and mental retardation (Balaji and Balaji, 2017; Snyder, 1965) or it can develop as a part of syndromes like Cowden’s syndrome (Witkop, 1971), Zimmerman-Laband syndrome (Guglielmi et al., 2019), Cross syndrome (Poulopoulos et al., 2011), Rutherford syndrome (Häkkinen and Csiszar, 2007), Ramon syndrome (Suhanya et al., 2010), Jones syndrome (Gita et al., 2014), Costello syndrome (Hennekam, 2003), Ectro-amelia syndrome (Morey and Higgins, 1990), and hyaline fibromatosis syndrome (Hamada et al., 1980); it may also be caused by poor oral hygiene-related gingival inflammation, puberty, pregnancy, or a side effect of the common medications such as calcium channel blockers (nifedipine and verapamil) (Livada and Shiloah, 2012); usually it betrays as an isolated condition as non-syndromic hereditary gingival fibromatosis (HGF) (Jorgenson and Cocker, 1974). HGF is the most common genetic form of GF and usually transmitted as an autosomal-dominant inheritance pattern. Nevertheless, sporadic cases and autosomal-recessive inheritance pedigrees have also been previously described (Majumder et al., 2013).

As reported previously, to date, four loci (2p22.1 [MIM: 135,300], 5q13–q22 [MIM: 605,544], 2p23.3–p22.3 [MIM: 609,955], and 11p15 [MIM: 611,010]) have been identified associated with HGF (Xiao et al., 2001; Ye et al., 2005; Zhu et al., 2007; Xiao et al., 2000); besides a heterozygous frameshift mutation in SOS1 (MIM: 182530) has been identified as causative for the autosomal-dominant HGF in a Brazilian family (Hart et al., 2002); and protein truncating mutations in REST (MIM: 600571) have been reported as the genetic cause of an autosomal-dominant inheritance pattern HGF across three independent families (Bayram et al., 2017). The estimated incidence of HGF is 1 per 175,000 of the population (Ahmed and Ali, 2015), and equal between males and females (Gawron et al., 2016; Odessey et al., 2006). Periodontal surgery including gingivectomy, gingivoplasty, and flap surgery can be applied to patients with HGF due to the aesthetic and functional concern, whereas the recurrence of hyperplasia is relatively high potentially owing to the underlying genetic predisposition (Chaurasia, 2014; Zhou et al., 2011). Therefore, it is important to explore the genetic causes and etiology accounting for HGF with expectations of the precise genetic diagnosis and of the potential gene therapy development, to help individuals who desire to circumvent the sufferings and improve their lives. We report a large Chinese family with 13 affected individuals clinically diagnosed with non-syndromic HGF, and with 12 unaffected family members. Genetic analysis of this large family led to the identification of a heterozygous mutations in a novel gene ZNF862.

Results and discussion

This four-generation pedigree with HGF was shown in Figure 1A. Visual inspection indicated severe or mild generalized GF in all affected members, whereas not in any unaffected member, from this family. No individuals in this family was ever exposed to the medication that may result in GF according to our investigation. No other congenital abnormality has been found for all family members. The clinical data of all individuals in this family are summarized in Table 1, and clinical photos of one control and three patients from each generation are available in Figure 1B. From a clinical perspective, the proband (IV-2) had the mild recurrence of hyperplasia during the second year after the periodontal surgery of gingivectomy, gingivoplasty, and flap surgery; while the patient II-2 had an additional clinical finding (hypertension, with no medication of calcium channel blockers treated).

Figure 1

Download asset Open asset

Pedigree and co-segregation analysis.

(A) The black arrow indicates the proband. The affected individuals are indicated with black-filled boxes in this family. Whole-exome sequencing (WES) was performed for the proband and nine other members (numbers indicated in red color). (B) The photographs of gingival overgrowth are revealed by intraoral examination in affected members compared with in the unaffected member in this family. (C) In this family, all affected individuals harbor the heterozygous variant (c.2812G > A) whereas the unaffected individuals are wild-type. (D) Schematic structure with putative domains of ZNF862 protein and the localization of the novel variant (red arrows). Green rectangles indicate KRAB (krueppel-associated box) domains; yellow rectangles indicate ZnF_TTF (zinc finger in transposases and transcription factors) domains; blue rectangle indicates Dimer_Tnp_hAT (hAT family C-terminal dimerization region) domain.

Table 1

Clinical Characteristics of the individuals in this family.

NA = Not available.

Members	Identified variant	Age at last exam	Gender	Gingival fibromatosis	Age of onset	Exposure to medication	Gingivectomy/recurrence	Other clinical finding
Ⅱ–1	Wild type	67 years	Male	−	–	No	-/-	No
Ⅱ–2	p.A938T (ZNF862)	70 years	Female	+	NA	No	-/-	Hypertension
II-3	p.A938T (ZNF862)	66 years	Male	+	NA	No	-/-	No
II-4	Wild type	66 years	Female	−	–	No	-/-	No
II-5	Wild type	64 years	Male	−	–	No	-/-	No
II-6	Wild type	64 years	Female	−	–	No	-/-	No
Ⅱ–7	Wild type	62 years	Male	−	–	No	-/-	No
Ⅱ–8	p.A938T (ZNF862)	60 years	Female	+	NA	No	-/-	No
Ⅱ–9	Wild type	54 years	Female	−	–	No	-/-	No
Ⅲ–1	Wild type	49 years	Female	−	–	No	-/-	No
Ⅲ–2	p.A938T (ZNF862)	47 years	Male	+	NA	No	-/-	No
Ⅲ–3	p.A938T (ZNF862)	45 years	Female	+ (mild)	NA	No	-/-	No
Ⅲ–4	p.A938T (ZNF862)	41 years	Male	+	NA	No	-/-	No
Ⅲ–5	Wild type	36 years	Female	−	–	No	-/-	No
Ⅲ–6	p.A938T (ZNF862)	36 years	Male	+	NA	No	-/-	No
III-7	p.A938T (ZNF862)	34 years	Female	+	NA	No	-/-	No
III-8	Wild type	32 years	Female	−	–	No	-/-	No
Ⅲ–9	Wild type	30 years	Male	−	–	No	-/-	No
Ⅳ–1	p.A938T (ZNF862)	26 years	Male	+	6–7 years	No	+/-	No
Ⅳ–2	p.A938T (ZNF862)	22 years	Female	+	6–7 years	No	+/+ (mild)	No
Ⅳ–3	p.A938T (ZNF862)	13 years	Female	+	2–3 years	No	-/-	No
Ⅳ–4	Wild type	11 years	Female	−	–	No	-/-	No
Ⅳ–5	p.A938T (ZNF862)	7 years	Male	+	2–3 years	No	-/-	No

Whole-exome sequencing (WES) was performed for the proband and nine other members (numbers indicated in red color in Figure 1A) in the family according to previously described protocols (Zhang et al., 2015). After analyzing all single nucleotide variants (SNVs) and indels for each gene including the promoter region, exons, splicing sites, introns, and untranslated region (UTR), only three variants ATP7B (c.3403G > A), CDADC1 (c.83–13G > T), ZNF862 (c.2812G > A) were identified to be co-segregated with the phenotype among the 10 family members who underwent WES. To further screen and validate the potential disease-causing variants co-segregation with the phenotype in this pedigree, conventional PCR was performed for sanger sequencing evaluation. Eventually, only ZNF862 (c.2812G > A), no any other variant, was validated to be co-segregated with the phenotype in all 23 alive members in this family, as is shown in Figure 1C.

ZNF862 gene harbors eight exons, encodes a 1169-amino acid protein, and is located on chromosome 7q36.1. It does not escape our notice that ZNF862 does not map to any one of the loci that were reported previously as be associated with HGF. ZNF862 is a predicted intracellular protein which function, to the best of our knowledge, is not yet clearly identified, as a zinc finger protein it may be involved in transcriptional regulation. ZNF862 is expressed ubiquitously across tissues (Fagerberg et al., 2014), it may play various roles under different physiological condition. Schwartz et al., 2019, suggested a plausible role of ZNF862 in matrix-producing metaplastic carcinoma. Peng et al., 2018, reported ZNF862 associated with IgE-mediated type-I hypersensitivity in children. The WES-identified heterozygous variant (p.A938T) in our study located close to and upstream of the Dimer_Tnp_hAT (hAT family C-terminal dimerization region) domain of ZNF862 protein (Figure 1).

The pathophysiological mechanisms underlying HGF remain largely elusive. Nevertheless, overproduction of extracellular matrix, particularly the major component collagen type I (COL1A1), may account for the gingival fibroblast overgrowth phenotype; meanwhile increased synthesis of TIMP-1 seems to be associated with COL1A1 excessive accumulation in HGF gingival fibroblasts (Gawron et al., 2018; Roman-Malo et al., 2019). Besides, Martelli-Junior et al. reported that overexpression of TGF-β1 and IL-6 in gingival fibroblasts play pivotal roles in increasing the synthesis of the COL1A1 along with other specific growth factors (Martelli-Junior et al., 2003). Systematically, transcriptomic analysis of HGF patients and controls performed by Han et al., 2019, demonstrated that regulatory network connection of TGF-β/SMAD signaling pathway and craniofacial development processes contributes to the molecular mechanism of clinical-pathological manifestations.

Aiming to explore the underlying mechanism related to the ZNF862 mutation, RNA sequencing was performed according to previously described protocols (Liang et al., 2019). The primary fibroblasts were isolated from fibromatic gingival specimens obtained from two patients who underwent gingivectomy (IV-1 and IV-2) using standard explant culture as described previously (de Andrade et al., 2001), while control cells were obtained from four independent age- and gender-matched controls who underwent crown lengthening surgery for the restorative purpose. The genes expression profiling and changes in patients over controls were summarized in Supplementary file 1, the samples average counts per million mapped reads (CPM) for each gene is shown, the expression fold-change (FC) and their statistical significance false discovery rate (FDR) are also provided. The expression profiling of attested profibrotic genes (Gao et al., 2019), which may be involved in HGF, was interrogated (Figure 2A). A part of such genes, including COL1A1, TGFB1/2, and SMAD1, was up-regulated in HGF fibroblasts compared with controls; whereas another part was down-regulated, suggesting that the special HGF in our study was attributed to a corresponding part of profibrotic factors, including TGF-β/SMAD1 signaling pathway and COL1A1. Besides RNA sequencing, the in vitro cultured gingival fibroblasts from control 1 underwent adenovirus-mediated delivery of short hairpin RNA (shRNA) targeting ZNF862. The expression of four most significant genes ACTA2, FOS, SMAD1, AGT, as well as COL1A1, under shRNA treatment were investigated (Figure 2B), the knockdown of ZNF862 in control gingival fibroblasts leads to the profibrotic genes up-regulation, reminiscent of the expression in the patients. Here, we speculate that the ZNF862 mutation in patients attenuate the protein function, similar with the ZNF862 shRNA, and promote special profibrotic genes expression to result in the HGF trait. Meanwhile, it did not escape our notice that the causative gene SOS1 was mildly down-regulated, which may be associated with the phenotype (Figure 2 and Figure 2—figure supplement 3). Nevertheless, the proliferation in control gingival fibroblasts was not significantly stimulated by the knockdown of ZNF862 during 5 days’ culture period (Figure 2—figure supplement 2), which may be partially resulting from the too short time course to observe the sophisticated physiological effect.

Figure 2 with 3 supplements see all

Download asset Open asset

RNA sequencing and knockdown analysis.

(A) Heatmap of the 70 attested profibrotic genes expression profiling in RNA sequencing. Rows represent genes and columns represent samples. The heatmap is color-coded based on an average normalization; red represents high expression value and green represents low expression value. (B) Real-time PCR analysis of *ACTA2*, *FOS*, *SMAD1*, *AGT*, *COL1A1*, and *ZNF862* mRNA abundance in untreated (Un), adenovirus delivered scrambler (Sc), adenovirus delivered *ZNF862* short hairpin RNA (shRNA) 1 (S1), and adenovirus delivered *ZNF862* shRNA 2 (S2) groups, respectively. Results are given as means ± SEM, dots indicate the relative values of gene expression levels. *p < 0.05, **p < 0.01 in comparison with Sc group. *N.S.* means not significant. (C) Heatmap of the 100 most up-regulated differentially expressed (DE) genes and 100 most down-regulated DE genes expression profiling in RNA sequencing. The color key is the same with (A). (D) Functional annotation chart of DE genes, the 20 most significantly enriched functions of DE genes were plotted. Functional pathways enrichment tests were based on KEGG database. Enrichment denotes the proportion of DE genes among the indicated pathway.

To screen all the genes differentially expressed (DE) between patients and controls, we quantified genes whose expression changed over twofold (|log₂ FC| ≥ 1.0), and we employed a cut-off of FDR < 0.05 due to the statistical power limitation of the small number of entire samples (N = 6). Finally, the transcriptomic comparison of HGF patients and controls yielded 597 DE transcripts, 355 up-regulated while 242 down-regulated (Figure S1). Of which, the 100 most up-regulated DE genes and 100 most down-regulated DE genes (Supplementary file 2) were shown in Figure 2C, furthermore the top 10 genes were illustrated accordingly. These DE genes probably correlated with the HGF trait. The 20 most significantly enriched functional pathways of all DE genes were shown (Figure 2D). Thereinto, lipoic acid metabolism is the most enriched cluster; notably the MAPK signaling pathway is associated with IL-6 pathway, besides several infection-related clusters in this chart is related with IL-6-induced autophagy pathway.

Taken together we can speculate that ZNF862 may execute as transcriptional repressors since ZNF862 has the zinc-finger DNA-binding domain, meanwhile this missense mutation may attenuate biological function of ZNF862 and hereby enhance the TGF-β/SMAD1 signaling pathway that increases profibrotic factors, particularly COL1A1 accumulation in gingival tissue and results in HGF. Nevertheless, the underlying mechanism of this novel mutation leading to HGF remains enigmatic, with no more specific experiment designed for exploring the biological function of ZNF862. Further studies are expected to expand the mutational spectrum of ZNF862 and the associations with phenotypes, and to investigate the physiological role of ZNF862 in the process of HGF.

In this study, we identify a missense variant (c.2812G > A) in a novel gene ZNF862 that causes autosomal-dominant HGF in a four-generation Chinese family. The functional study supports a biological role of ZNF862 for increasing the profibrotic factors, particularly COL1A1 synthesis in gingiva and hence resulting in HGF. This variant has been absent among the large population according to the report of Exome Aggregation Consortium (ExAC) database, 1000 Genomes, and Genome Aggregation Database. The log of the odds (LOD) score is calculated as 3.6 based on the dominant segregations. According to the combining criteria set guidelines by the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen) (Oza et al., 2018; Richards et al., 2015), the missense variant is rated as ‘Pathogenic’ with the following evidence rules: PP1_Strong (LOD > 3), PS3 (in vitro functional study support), and PM2 (absent from controls). We can speculate that this finding may shed lights on the precise genetic diagnosis and on the potential therapy development, helping individuals continue their lives better.

Materials and methods

Samples collection and WES

Request a detailed protocol

The participants from this family were recruited from Nanjing Stomatological Hospital. The genomic DNA was extracted from the peripheral blood. WES was performed for probands and nine other members (numbers indicated in red color in Figure 1A) in the family. Average WES sequencing depth of target region for each sample was higher than 100, and 20× coverage for more than 95% targeted bases was achieved for each sample. Sequencing reads were aligned to the human genome reference sequence (hg19) through the in-house programs. Genome Analysis Toolkit (GATK) was used to call SNVs and indels; and Annotation of Genetic Variants (ANNOVAR) was used for the annotation. A custom Perl script (Source code 1) was used for retrieving the variants corresponding to the inheritance model of this family and that do not appear in unaffected individuals. All the WES data underwent the same quality control filtering and pruning procedures to maximize parity. The following primers are used to amplify products for Sanger sequencing: ZNF862 (NM_001099220.3) forward: 5'-TCATGCGCGGCTTCCACTTTGT; reverse: 5'-AGCACTCGAAATACCTGGCCAG; ATP7B (NM_000053.4) forward: 5'-TCCAGTCGGTAACCTGTTCAC; reverse: 5'-GCACAGCAGAGGCAATCAC; CDADC1 (NM_030911.4) forward: 5'-CGACTAAGGCTGAAGCGTCT; reverse: 5'-AATTGGGGGAAATTATGTGG.

Transcriptome sequencing

Request a detailed protocol

The primary fibroblasts were isolated from fibromatic gingival specimens obtained from the patients who underwent gingivectomy (IV-1 as patient 1 and IV-2 as patient 2) using standard explant culture, while control cells were obtained from independent age- and gender-matched controls (two males at 25 and 26 years of age as controls 1 and 2, respectively, and two females at 21 and 23 years of age as controls 3 and 4, respectively) who underwent crown lengthening surgery for the restorative purpose.

Total RNA samples were extracted from the cultured fibroblast cells using Trizol, RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA). A total amount of 1 μg RNA per sample was used as input material for the RNA sample preparations. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in First Strand Synthesis Reaction Buffer. First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase (RNase H-). Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3’ ends of DNA fragments, adaptor with hairpin loop structure was ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 370–420 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, MA). Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers, and Index Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. The library preparations were sequenced on an Illumina Novaseq platform and 150 bp paired-end reads were generated.

Raw data were processed through an in-house bioinformatics workflow. Paired-end clean reads were aligned to the reference genome ensembl release-97 (http://asia.ensembl.org/Homo_sapiens/Info/Index) using Hisat2 v2.0.5, featureCounts v1.5.0-p3 was used to count the reads numbers mapped to each gene. And then the CPM and fragments per kilobase of transcript per million mapped reads (FPKM) of each gene was calculated based on the length of the gene and reads count mapped to this gene. At least 45 million clean mapped reads per sample were yielded. The clean bases after filtration for each sample are higher than 6.7 Giga-base, with the effective rate higher than 99.7%, the Q20 bases is higher than 94% for each sample. The GC content for each sample is ranging from 51% to 55%. DE analysis of two groups (patients and controls) was performed using the edgeR R package (3.22.5). The resulting p-values were adjusted using the Benjamini and Hochberg’s approach for controlling the FDR. Corrected p-value of 0.05 and absolute FC of 2 were set as the threshold for significantly DE. We used clusterProfiler R package to test the statistical enrichment of DE genes in KEGG pathways, corrected p-value less than 0.05 was considered significantly enriched. KEGG is a database resource for understanding high-level functions and utilities of the biological system (http://www.genome.jp/kegg/).

ShRNA interference

Request a detailed protocol

Besides RNA sequencing, the in vitro cultured fibroblasts (from control 1, control 2, and control 3, respectively) underwent adenovirus-mediated delivery of shRNA targeting ZNF862, the following sequences were used for the ZNF862 shRNA interference experiments, ZNF862 shRNA 1: GCTCTGTTCTGCTCTGCTTGC; ZNF862 shRNA 2: GGATTTACATCCGCTACTTCA; scrambler: GCACCCAGTCCGCCCTGAGCAAA, at 72 hr after the adenovirus-mediated U6 promoter-driven shRNA delivered into the cells (MOI = 70) real-time PCR was performed using 18s rRNA as an internal control. Cells between the third and sixth passages were used for all abovementioned experiments. The following primers are used to amplify products for the real-time PCR: FOS (NM_005252) forward: 5'-GCCGGGGATAGCCTCTCTTACT; reverse: 5'-CCAGGTCCGTGCAGAAGTC; SMAD1 (NM_001003688) forward: 5'-AGAGACTTCTTGGGTGGAAACA; reverse: 5'-ATGGTGACACAGTTACTCGGT; AGT (NM_000029) forward: 5'-GCTGACAGGCTACAGGCAATC; reverse: 5'-TGTGAACACGCCCACCACC; CCL2 (NM_002982) forward: 5'-CAGCCAGATGCAATCAATGCC; reverse: 5'-TGGAATCCTGAACCCACTTCT; COL1A1 (NM_000088.3) forward: 5'-GAGTCTACATGTCTAGGGTCT; reverse: 5'-CACGTCATCGCACAACACCT; ZNF862 (NM_001099220.3) forward: 5'-CTTACTCCAGGAGGAATGGGT; reverse: 5'-CTCCCATGTAGCCCATCTGT; ACTA2 (NM_001613.3) forward: 5'-CCAGACATCAGGGGGTGAT; reverse: 5'-TGGTGCCAGATCTTTTCCAT; FN1 (NM_001365524.2) forward: 5'-CCATAAAGGGCAACCAAGAG; reverse: 5'-AAACCAATTCTTGGAGCAGG; SOS1 (NM_005633.4) forward: 5'- GAAACCCTTTATCTCTCCCAGT; reverse: 5'- CTTGTCAGCACACATTGCCACT; 18srRNA (NR_145822.1) forward: 5'-CGAACGTCTGCCCTATCAACT; reverse: 5'-CAGACTTGCCCTCCAATGGATCCTCGTT.

Proliferation assay

Request a detailed protocol

Cell Counting Kit 8 (CCK8, Dojindo Molecular Technologies) proliferation assay was performed with the above-mentioned in vitro cultured fibroblasts. The cells in the logarithmic growth phase were plated overnight in a 96-well plate at an initial density of 2000 per well. Starting with adenovirus infection (MOI = 70), 10 μl CCK8 per well was added to the culture medium on the days 1, 2, 3, 4, and 5. The fluorescent optical density value was measured using a microplate reader. The in vitro cultured fibroblasts underwent adenovirus-mediated delivery of ZNF862 siRNA (NM_001099220.3).

Statistical analysis

Request a detailed protocol

All data from the semi-quantitative real-time PCR analyses are representative of at least three independent experiments. Data are presented as the mean ± SEM of at least three independent experiments and differences (in comparison with Sc group) are considered statistically significant at p < 0.05 using Student’s t-test. *p < 0.05, **p < 0.01. N.S. means not significant. All analyses were performed using R v4.0.2.

Web resources

Request a detailed protocol

1000 Genomes, http://www.internationalgenome.org/
ClinVar, https://www.ncbi.nlm.nih.gov/clinvar/
ExAC, http://exac.broadinstitute.org
Genome Aggregation Database, https://gnomad.broadinstitute.org/
GenBank, http://www.ncbi.nlm.nih.gov/genbank/
OMIM, http://www.omim.org/
UniProt, http://www.uniprot.org/

Data availability

The sequencing data supporting this study have been deposited in the China Genebank Nucleotide Sequence Archive (https://db.cngb.org/cnsa, accession number CNP0000995).

The following data sets were generated

1. Wu J
2. Chen D
3. Huang H
4. Luo N
5. Chen H
6. Zhao J
7. Wang Y
8. Zhao T
9. Huang S
10. Ren Y
11. Zhai T
12. Sun W
13. Li W
(2021) China National GeneBank
Study on the causative gene of hereditary gingival fibromatosis.

https://doi.org/10.26036/CNP0000995

References

1. Ahmed S
2. Ali Z
(2015)
Rare case of idiopathic gingival fibromatosis affecting primary dentition

Journal of Ayub Medical College, Abbottabad 27:933–935.
- PubMed
- Google Scholar
1. Balaji P
2. Balaji SM
(2017) Gingival fibromatosis with hypertrichosis syndrome: Case series of rare syndrome
Indian Journal of Dental Research 28:457–460.

https://doi.org/10.4103/ijdr.IJDR_367_17
- PubMed
- Google Scholar
(2017) REST Final-Exon-Truncating Mutations Cause Hereditary Gingival Fibromatosis
American Journal of Human Genetics 101:149–156.

https://doi.org/10.1016/j.ajhg.2017.06.006
- PubMed
- Google Scholar
1. Chaurasia A
(2014) Hereditary gingival fibromatosis
National Journal of Maxillofacial Surgery 5:42–46.

https://doi.org/10.4103/0975-5950.140171
- PubMed
- Google Scholar
1. de Andrade CR
2. Cotrin P
3. Graner E
4. Almeida OP
5. Sauk JJ
6. Coletta RD
(2001) Transforming growth factor-beta1 autocrine stimulation regulates fibroblast proliferation in hereditary gingival fibromatosis
Journal of Periodontology 72:1726–1733.

https://doi.org/10.1902/jop.2001.72.12.1726
- PubMed
- Google Scholar
1. Fagerberg L
2. Hallström BM
3. Oksvold P
4. Kampf C
5. Djureinovic D
6. Odeberg J
7. Habuka M
8. Tahmasebpoor S
9. Danielsson A
10. Edlund K
11. Asplund A
12. Sjöstedt E
13. Lundberg E
14. Szigyarto CA-K
15. Skogs M
16. Takanen JO
17. Berling H
18. Tegel H
19. Mulder J
20. Nilsson P
21. Schwenk JM
22. Lindskog C
23. Danielsson F
24. Mardinoglu A
25. Sivertsson A
26. von Feilitzen K
27. Forsberg M
28. Zwahlen M
29. Olsson I
30. Navani S
31. Huss M
32. Nielsen J
33. Ponten F
34. Uhlén M
(2014) Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics
Molecular & Cellular Proteomics 13:397–406.

https://doi.org/10.1074/mcp.M113.035600
- PubMed
- Google Scholar
1. Gao Q
2. Yang K
3. Chen D
4. Song Y
5. Qiao W
6. Sun X
7. Meng L
8. Bian Z
(2019) Antifibrotic Potential of MiR-335-3p in Hereditary Gingival Fibromatosis
Journal of Dental Research 98:1140–1149.

https://doi.org/10.1177/0022034519863300
- PubMed
- Google Scholar
(2016) Gingival fibromatosis: clinical, molecular and therapeutic issues
Orphanet Journal of Rare Diseases 11:9.

https://doi.org/10.1186/s13023-016-0395-1
- PubMed
- Google Scholar
(2018) TIMP-1 association with collagen type I overproduction in hereditary gingival fibromatosis
Oral Diseases 24:1581–1590.

https://doi.org/10.1111/odi.12938
- PubMed
- Google Scholar
(2014) Idiopathic gingival fibromatosis associated with progressive hearing loss: A nonfamilial variant of Jones syndrome
Contemporary Clinical Dentistry 5:260–263.

https://doi.org/10.4103/0976-237X.132387
- PubMed
- Google Scholar
(2019) Zimmermann-Laband-1 Syndrome: Clinical, Histological, and Proteomic Findings of a 3-Year-Old Patient with Hereditary Gingival Fibromatosis
Biomedicines 7:E48.

https://doi.org/10.3390/biomedicines7030048
- PubMed
- Google Scholar
1. Häkkinen L
2. Csiszar A
(2007) Hereditary gingival fibromatosis: characteristics and novel putative pathogenic mechanisms
Journal of Dental Research 86:25–34.

https://doi.org/10.1177/154405910708600104
- PubMed
- Google Scholar
1. Hamada K
2. Ohdo S
3. Hayakawa K
4. Kikuchi I
5. Koono M
(1980) Puretić syndrome--gingival fibromatosis with hyaline fibromas
Jinrui Idengaku Zasshi. The Japanese Journal of Human Genetics 25:249–255.

https://doi.org/10.1007/BF01997703
- PubMed
- Google Scholar
1. Han SK
2. Kong J
3. Kim S
4. Lee JH
5. Han DH
(2019) Exomic and transcriptomic alterations of hereditary gingival fibromatosis
Oral Diseases 25:1374–1383.

https://doi.org/10.1111/odi.13093
- PubMed
- Google Scholar
1. Hart TC
2. Zhang Y
3. Gorry MC
4. Hart PS
5. Cooper M
6. Marazita ML
7. Marks JM
8. Cortelli JR
9. Pallos D
(2002) A mutation in the SOS1 gene causes hereditary gingival fibromatosis type 1
American Journal of Human Genetics 70:943–954.

https://doi.org/10.1086/339689
- PubMed
- Google Scholar
1. Hennekam RCM
(2003) Costello syndrome: an overview
American Journal of Medical Genetics. Part C, Seminars in Medical Genetics 117C:42–48.

https://doi.org/10.1002/ajmg.c.10019
- PubMed
- Google Scholar
1. Jorgenson RJ
2. Cocker ME
(1974) Variation in the Inheritance and Expression of Gingival Fibromatosis
Journal of Periodontology 45:472–477.

https://doi.org/10.1902/jop.1974.45.7.472
- PubMed
- Google Scholar
1. Liang WC
2. Wong CW
3. Liang PP
4. Shi M
5. Cao Y
6. Rao ST
7. Tsui SKW
8. Waye MMY
9. Zhang Q
10. Fu WM
11. Zhang JF
(2019) Translation of the circular RNA circβ-catenin promotes liver cancer cell growth through activation of the Wnt pathway
Genome Biology 20:84.

https://doi.org/10.1186/s13059-019-1685-4
- PubMed
- Google Scholar
1. Livada R
2. Shiloah J
(2012)
Gummy smile: could it be genetic? Hereditary gingival fibromatosis

The Journal of the Tennessee Dental Association 92:23–26.
- PubMed
- Google Scholar
1. Majumder P
2. Nair V
3. Mukherjee M
4. Ghosh S
5. Dey SK
(2013) The autosomal recessive inheritance of hereditary gingival fibromatosis
Case Reports in Dentistry 2013:432864.

https://doi.org/10.1155/2013/432864
- PubMed
- Google Scholar
(2003) Effect of transforming growth factor-beta1, interleukin-6, and interferon-gamma on the expression of type I collagen, heat shock protein 47, matrix metalloproteinase (MMP)-1 and MMP-2 by fibroblasts from normal gingiva and hereditary gingival fibromatosis
Journal of Periodontology 74:296–306.

https://doi.org/10.1902/jop.2003.74.3.296
- PubMed
- Google Scholar
1. Morey MA
2. Higgins RR
(1990) Ectro-amelia syndrome associated with an interstitial deletion of 7q
American Journal of Medical Genetics 35:95–99.

https://doi.org/10.1002/ajmg.1320350118
- PubMed
- Google Scholar
(2006) Hereditary gingival fibromatosis: aggressive 2-stage surgical resection in lieu of traditional therapy
Annals of Plastic Surgery 57:557–560.

https://doi.org/10.1097/01.sap.0000229059.20539.0e
- PubMed
- Google Scholar
1. Oza AM
2. DiStefano MT
3. Hemphill SE
4. Cushman BJ
5. Grant AR
6. Siegert RK
7. Shen J
8. Chapin A
9. Boczek NJ
10. Schimmenti LA
11. Murry JB
12. Hasadsri L
13. Nara K
14. Kenna M
15. Booth KT
16. Azaiez H
17. Griffith A
18. Avraham KB
19. Kremer H
20. Rehm HL
21. Amr SS
22. Abou Tayoun AN
23. ClinGen Hearing Loss Clinical Domain Working Group
(2018) Expert specification of the ACMG/AMP variant interpretation guidelines for genetic hearing loss
Human Mutation 39:1593–1613.

https://doi.org/10.1002/humu.23630
- PubMed
- Google Scholar
1. Peng C
2. Cardenas A
3. Rifas-Shiman SL
4. Hivert M-F
5. Gold DR
6. Platts-Mills TA
7. Lin X
8. Oken E
9. Baccarelli AA
10. Litonjua AA
11. DeMeo DL
(2018) Epigenome-wide association study of total serum immunoglobulin E in children: a life course approach
Clinical Epigenetics 10:55.

https://doi.org/10.1186/s13148-018-0488-x
- PubMed
- Google Scholar
(2011) Current concepts on gingival fibromatosis-related syndromes
Journal of Investigative and Clinical Dentistry 2:156–161.

https://doi.org/10.1111/j.2041-1626.2011.00054.x
- PubMed
- Google Scholar
(2015) Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology
Genetics in Medicine 17:405–424.

https://doi.org/10.1038/gim.2015.30
- PubMed
- Google Scholar
(2019) Fibroblasts Collagen Production and Histological Alterations in Hereditary Gingival Fibromatosis
Diseases (Basel, Switzerland) 7:E39.

https://doi.org/10.3390/diseases7020039
- PubMed
- Google Scholar
1. Schwartz CJ
2. Dolgalev I
3. Yoon E
4. Osman I
5. Heguy A
6. Vega-Saenz de Miera EC
7. Nimeh D
8. Jour G
9. Darvishian F
(2019) Microglandular adenosis is an advanced precursor breast lesion with evidence of molecular progression to matrix-producing metaplastic carcinoma
Human Pathology 85:65–71.

https://doi.org/10.1016/j.humpath.2018.10.021
- PubMed
- Google Scholar
1. Snyder CH
(1965) Syndrome of gingival hyperplasia, hirsutism, and convulsions; dilantin intoxication without dilantin
The Journal of Pediatrics 67:499–502.

https://doi.org/10.1016/s0022-3476(65)80413-9
- PubMed
- Google Scholar
(2010) Cherubism combined with epilepsy, mental retardation and gingival fibromatosis (Ramon syndrome): a case report
Head and Neck Pathology 4:126–131.

https://doi.org/10.1007/s12105-009-0155-9
- PubMed
- Google Scholar
1. Witkop CJ
(1971)
Heterogeneity in gingival fibromatosis

Birth Defects Original Article Series 7:210–221.
- PubMed
- Google Scholar
1. Xiao S
2. Wang X
3. Qu B
4. Yang M
5. Liu G
6. Bu L
7. Wang Y
8. Zhu L
9. Lei H
10. Hu L
11. Zhang X
12. Liu J
13. Zhao G
14. Kong X
(2000) Refinement of the locus for autosomal dominant hereditary gingival fibromatosis (GINGF) to a 3.8-cM region on 2p21
Genomics 68:247–252.

https://doi.org/10.1006/geno.2000.6285
- PubMed
- Google Scholar
1. Xiao S
2. Bu L
3. Zhu L
4. Zheng G
5. Yang M
6. Qian M
7. Hu L
8. Liu J
9. Zhao G
10. Kong X
(2001) A new locus for hereditary gingival fibromatosis (GINGF2) maps to 5q13-q22
Genomics 74:180–185.

https://doi.org/10.1006/geno.2001.6542
- PubMed
- Google Scholar
1. Ye X
2. Shi L
3. Cheng Y
4. Peng Q
5. Huang S
6. Liu J
7. Huang M
8. Peng B
9. Bian Z
(2005) A novel locus for autosomal dominant hereditary gingival fibromatosis, GINGF3, maps to chromosome 2p22.3-p23.3
Clinical Genetics 68:239–244.

https://doi.org/10.1111/j.1399-0004.2005.00488.x
- PubMed
- Google Scholar
1. Zhang L
2. Zhou Y
3. Cheng C
4. Cui H
5. Cheng L
6. Kong P
7. Wang J
8. Li Y
9. Chen W
10. Song B
11. Wang F
12. Jia Z
13. Li L
14. Li Y
15. Yang B
16. Liu J
17. Shi R
18. Bi Y
19. Zhang Y
20. Wang J
21. Zhao Z
22. Hu X
23. Yang J
24. Li H
25. Gao Z
26. Chen G
27. Huang X
28. Yang X
29. Wan S
30. Chen C
31. Li B
32. Tan Y
33. Chen L
34. He M
35. Xie S
36. Li X
37. Zhuang X
38. Wang M
39. Xia Z
40. Luo L
41. Ma J
42. Dong B
43. Zhao J
44. Song Y
45. Ou Y
46. Li E
47. Xu L
48. Wang J
49. Xi Y
50. Li G
51. Xu E
52. Liang J
53. Yang X
54. Guo J
55. Chen X
56. Zhang Y
57. Li Q
58. Liu L
59. Li Y
60. Zhang X
61. Yang H
62. Lin D
63. Cheng X
64. Guo Y
65. Wang J
66. Zhan Q
67. Cui Y
(2015) Genomic analyses reveal mutational signatures and frequently altered genes in esophageal squamous cell carcinoma
American Journal of Human Genetics 96:597–611.

https://doi.org/10.1016/j.ajhg.2015.02.017
- PubMed
- Google Scholar
1. Zhou M
2. Xu L
3. Meng HX
(2011)
Diagnosis and treatment of a hereditary gingival fibromatosis case

The Chinese Journal of Dental Research 14:155–158.
- PubMed
- Google Scholar
1. Zhu Y
2. Zhang W
3. Huo Z
4. Zhang Y
5. Xia Y
6. Li B
7. Kong X
8. Hu L
(2007) A novel locus for maternally inherited human gingival fibromatosis at chromosome 11p15
Human Genetics 121:113–123.

https://doi.org/10.1007/s00439-006-0283-1
- PubMed
- Google Scholar

Article and author information

Author details

Juan Wu

Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China

Contribution
Investigation, Resources, Validation, Writing – review and editing

Contributed equally with
Dongna Chen and Hui Huang

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4494-6570
Dongna Chen

BGI Genomics, BGI-Shenzhen, Shenzhen, China

Contribution
Formal analysis, Resources, Validation

Contributed equally with
Juan Wu and Hui Huang

Competing interests
is employee of BGI Genomics

"This ORCID iD identifies the author of this article:" 0000-0002-7957-8379
Hui Huang

BGI Genomics, BGI-Shenzhen, Shenzhen, China

Contribution
Resources, Validation

Contributed equally with
Juan Wu and Dongna Chen

Competing interests
is employee of BGI Genomics
Ning Luo

Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China

Contribution
Data curation

Competing interests
No competing interests declared
Huishuang Chen

BGI Genomics, BGI-Shenzhen, Shenzhen, China

Contribution
Data curation, Methodology, Validation, Writing – review and editing

Competing interests
is employee of BGI Genomics
Junjie Zhao

Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China

Contribution
Formal analysis, Methodology

Competing interests
No competing interests declared
Yanyan Wang

BGI Genomics, BGI-Shenzhen, Shenzhen, China

Contribution
Formal analysis, Methodology, Writing – review and editing

Competing interests
is employee of BGI Genomics
Tian Zhao

Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China

Contribution
Data curation, Investigation

Competing interests
No competing interests declared
Siyuan Huang
1. BGI Genomics, BGI-Shenzhen, Shenzhen, China
2. Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
Contribution
Data curation, Formal analysis, Methodology, Software

Competing interests
No competing interests declared
Yang Ren

Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China

Contribution
Investigation

Competing interests
No competing interests declared
Teng Zhai

BGI Genomics, BGI-Shenzhen, Shenzhen, China

Contribution
Resources

Competing interests
is employee of BGI Genomics
Weibin Sun

Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China

Contribution
Investigation

Competing interests
No competing interests declared
Houxuan Li

Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China

Contribution
Conceptualization, Funding acquisition, Project administration, Resources, Writing – review and editing

For correspondence
lihouxuan3435_0@163.com

Competing interests
No competing interests declared
Wei Li

BGI Genomics, BGI-Shenzhen, Shenzhen, China

Contribution
Conceptualization, Data curation, Formal analysis, Project administration, Supervision, Visualization, Writing – original draft, Writing – review and editing

For correspondence
chorway@pku.edu.cn

Competing interests
is employee of BGI Genomics

"This ORCID iD identifies the author of this article:" 0000-0003-4475-531X

Funding

National Natural Science Foundation of China (51772144)

Houxuan Li

Medical Science and Technology Development Foundation, Nanjing Department of Health (YKK20152)

Juan Wu

Nanjing Clinical Research Center for Oral Diseases, Nanjing Department of Health (2019060009)

Juan Wu

Jiangsu Province Natural Science Foundation of China (BK20200149)

Juan Wu

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank all the participants in this study, and the support by the funding from National Natural Science Foundation of China (51772144); the Medical Science and Technology Development Foundation, Nanjing Department of Health (YKK18121); Nanjing Clinical Research Center for Oral Diseases, Nanjing Department of Health (No. 2019060009); Jiangsu Province Natural Science Foundation of China (BK20200149).

Ethics

Human subjects: The usage and handling of human samples in this study was approved by the Institutional Review Board on Bioethics and Biosafety of BGI (IRB No. 19059) and the written informed consent obtained from each participant. Clinical investigation was performed in accordance with the Declaration of Helsinki.

Copyright

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.