1. Genetics and Genomics

A novel gene ZNF862 causes hereditary gingival fibromatosis

  1. Juan Wu
  2. Dongna Chen
  3. Hui Huang
  4. Ning Luo
  5. Huishuang Chen
  6. Junjie Zhao
  7. Yanyan Wang
  8. Tian Zhao
  9. Siyuan Huang
  10. Yang Ren
  11. Teng Zhai
  12. Weibin Sun
  13. Houxuan Li  Is a corresponding author
  14. Wei Li  Is a corresponding author
  1. Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, China
  2. BGI Genomics, BGI-Shenzhen, China
  3. Academy for Advanced Interdisciplinary Studies, Peking University, China
https://doi.org/10.7554/eLife.66646
Share this article
Cite this article
  1. Juan Wu
  2. Dongna Chen
  3. Hui Huang
  4. Ning Luo
  5. Huishuang Chen
  6. Junjie Zhao
  7. Yanyan Wang
  8. Tian Zhao
  9. Siyuan Huang
  10. Yang Ren
  11. Teng Zhai
  12. Weibin Sun
  13. Houxuan Li
  14. Wei Li
(2022)
A novel gene ZNF862 causes hereditary gingival fibromatosis
eLife 11:e66646.
https://doi.org/10.7554/eLife.66646
  2. Download BibTeX
  3. Download .RIS

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.

MembersIdentified variantAge at last examGenderGingival fibromatosisAge of onsetExposure to medicationGingivectomy/recurrenceOther clinical finding
Ⅱ–1Wild type67 yearsMale−–No-/-No
Ⅱ–2p.A938T (ZNF862)70 yearsFemale+NANo-/-Hypertension
II-3p.A938T (ZNF862)66 yearsMale+NANo-/-No
II-4Wild type66 yearsFemale−–No-/-No
II-5Wild type64 yearsMale−–No-/-No
II-6Wild type64 yearsFemale−–No-/-No
Ⅱ–7Wild type62 yearsMale−–No-/-No
Ⅱ–8p.A938T (ZNF862)60 yearsFemale+NANo-/-No
Ⅱ–9Wild type54 yearsFemale−–No-/-No
Ⅲ–1Wild type49 yearsFemale−–No-/-No
Ⅲ–2p.A938T (ZNF862)47 yearsMale+NANo-/-No
Ⅲ–3p.A938T (ZNF862)45 yearsFemale+ (mild)NANo-/-No
Ⅲ–4p.A938T (ZNF862)41 yearsMale+NANo-/-No
Ⅲ–5Wild type36 yearsFemale−–No-/-No
Ⅲ–6p.A938T (ZNF862)36 yearsMale+NANo-/-No
III-7p.A938T (ZNF862)34 yearsFemale+NANo-/-No
III-8Wild type32 yearsFemale−–No-/-No
Ⅲ–9Wild type30 yearsMale−–No-/-No
Ⅳ–1p.A938T (ZNF862)26 yearsMale+6–7 yearsNo+/-No
Ⅳ–2p.A938T (ZNF862)22 yearsFemale+6–7 yearsNo+/+ (mild)No
Ⅳ–3p.A938T (ZNF862)13 yearsFemale+2–3 yearsNo-/-No
Ⅳ–4Wild type11 yearsFemale−–No-/-No
Ⅳ–5p.A938T (ZNF862)7 yearsMale+2–3 yearsNo-/-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 (|log2 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

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

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.
    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
    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.
    1. Witkop CJ
    (1971)
    Heterogeneity in gingival fibromatosis
    Birth Defects Original Article Series 7:210–221.
    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.

Article and author information

Author details

  1. 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
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4494-6570

  2. 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
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7957-8379

  3. 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

  4. Ning Luo

    Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China
    Contribution
    Data curation
    Competing interests
    No competing interests declared

  5. Huishuang Chen

    BGI Genomics, BGI-Shenzhen, Shenzhen, China
    Contribution
    Data curation, Methodology, Validation, Writing – review and editing
    Competing interests
    is employee of BGI Genomics

  6. 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

  7. Yanyan Wang

    BGI Genomics, BGI-Shenzhen, Shenzhen, China
    Contribution
    Formal analysis, Methodology, Writing – review and editing
    Competing interests
    is employee of BGI Genomics

  8. 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

  9. 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

  10. Yang Ren

    Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China
    Contribution
    Investigation
    Competing interests
    No competing interests declared

  11. Teng Zhai

    BGI Genomics, BGI-Shenzhen, Shenzhen, China
    Contribution
    Resources
    Competing interests
    is employee of BGI Genomics

  12. Weibin Sun

    Department of Periodontology, Nanjing Stomatological Hospital, Medical School of Nanjing University, Nanjing, China
    Contribution
    Investigation
    Competing interests
    No competing interests declared

  13. 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

  14. 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
    ORCID icon "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.

Version history

  1. Preprint posted: May 1, 2020 (view preprint)
  2. Received: January 19, 2021
  3. Accepted: February 9, 2022
  4. Accepted Manuscript published: February 10, 2022 (version 1)
  5. Accepted Manuscript updated: February 15, 2022 (version 2)
  6. Version of Record published: February 18, 2022 (version 3)

Copyright

© 2022, Wu 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

  • 712
    views
  • 121
    downloads
  • 6
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Juan Wu
  2. Dongna Chen
  3. Hui Huang
  4. Ning Luo
  5. Huishuang Chen
  6. Junjie Zhao
  7. Yanyan Wang
  8. Tian Zhao
  9. Siyuan Huang
  10. Yang Ren
  11. Teng Zhai
  12. Weibin Sun
  13. Houxuan Li
  14. Wei Li
(2022)
A novel gene ZNF862 causes hereditary gingival fibromatosis
eLife 11:e66646.
https://doi.org/10.7554/eLife.66646

Share this article

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

Categories and tags

Research organism

Further reading

    1. Computational and Systems Biology
    2. Genetics and Genomics

    Imputation of 3D genome structure by genetic–epigenetic interaction modeling in mice

    Lauren Kuffler, Daniel A Skelly ... Gregory W Carter
    Research Article

    Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic–epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic–epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

    1. Genetics and Genomics
    2. Neuroscience

    Genetic Chimerism: Marmosets contain multitudes

    Kenneth Chiou, Noah Snyder-Mackler
    Insight

    Single-cell RNA sequencing reveals the extent to which marmosets carry genetically distinct cells from their siblings.

    1. Genetics and Genomics

    Elimination of subtelomeric repeat sequences exerts little effect on telomere essential functions in Saccharomyces cerevisiae

    Can Hu, Xue-Ting Zhu ... Jin-Qiu Zhou
    Research Article

    Telomeres, which are chromosomal end structures, play a crucial role in maintaining genome stability and integrity in eukaryotes. In the baker’s yeast Saccharomyces cerevisiae, the X- and Y’-elements are subtelomeric repetitive sequences found in all 32 and 17 telomeres, respectively. While the Y’-elements serve as a backup for telomere functions in cells lacking telomerase, the function of the X-elements remains unclear. This study utilized the S. cerevisiae strain SY12, which has three chromosomes and six telomeres, to investigate the role of X-elements (as well as Y’-elements) in telomere maintenance. Deletion of Y’-elements (SY12YΔ), X-elements (SY12XYΔ+Y), or both X- and Y’-elements (SY12XYΔ) did not impact the length of the terminal TG1-3 tracks or telomere silencing. However, inactivation of telomerase in SY12YΔ, SY12XYΔ+Y, and SY12XYΔ cells resulted in cellular senescence and the generation of survivors. These survivors either maintained their telomeres through homologous recombination-dependent TG1-3 track elongation or underwent microhomology-mediated intra-chromosomal end-to-end joining. Our findings indicate the non-essential role of subtelomeric X- and Y’-elements in telomere regulation in both telomerase-proficient and telomerase-null cells and suggest that these elements may represent remnants of S. cerevisiae genome evolution. Furthermore, strains with fewer or no subtelomeric elements exhibit more concise telomere structures and offer potential models for future studies in telomere biology.