Sage Journals: Discover world-class research

Abstract

Periodontitis is a common complex inflammatory disease of the oral cavity. It is characterized by inflammation of gingival tissues and alveolar bone loss. Recently, a genome-wide association study and 2 genome-wide association study meta-analyses found 2 associated regions (haplotype blocks) at the inhibitory immune receptor gene SIGLEC5 to increase the risk for periodontitis. The aims of the current study were the identification of the putative causal variants underlying these associations, characterization of their molecular biological effects, and validation of SIGLEC5 as the target gene. We mapped the associated single-nucleotide polymorphisms to DNA elements with predictive features of regulatory functions and screened the associated alleles for transcription factor (TF) binding sites. Antibody electrophoretic mobility shift assays (EMSAs) with allele-specific probes were used to identify TF binding and to quantify allele-specific effects on binding affinities. Luciferase reporter assays were used to quantify the effect directions and allele-specific strength of the associated regulatory elements. We used CRISPR-dCas9 gene activation to validate SIGLEC5 as a target of the association. EMSA in peripheral blood mononuclear cells showed that E-26 transformation–specific TF-related gene (ERG) binds at rs11084095, with almost complete loss of binding at the minor A-allele. Allele-specific reporter genes showed enhancer function of the DNA sequence at rs11084095, which was abrogated in the background of the A-allele. EMSA in B lymphocytes showed that TF MAF bZIP (MAFB) binds at the common G-allele of rs4284742, whereas the minor A-allele reduced TF binding by 69%, corresponding to 9-fold reduction of luciferase reporter gene activity by the A-allele. Using CRISPR-dCas9, we showed that the enhancer at rs4284742 strongly activated SIGLEC5 expression, validating this gene as the target gene of the association. We conclude that rs11084095 and rs4284742 are putatively causal for the genome-wide significant associations with periodontitis at SIGLEC5 that impair ERG and MAFB binding, respectively.

Keywords

oral inflammation ERG transcription factor MAFB transcription factor mouth mucosa enhancer elements genetic

Introduction

Periodontitis is a common complex inflammatory disease of the oral cavity and is characterized by chronic inflammation of the gingiva in conjunction with destruction of connective tissue and, subsequently, alveolar bone and tooth loss. Periodontitis has a range of manifestations that differ among individuals in severity and progression of tissue destruction and age of disease onset. The basis of phenotypic variation of complex diseases largely is the genetic variability among individuals (Timpson et al. 2018). The genetic variants mostly are single-nucleotide polymorphisms (SNPs; 1000 Genomes Project et al. 2015), and the number of functional SNPs and their biological consequences, magnitude of effects, and interactions with one another and with environmental factors shape disease manifestations (Boyle et al. 2017; Wray et al. 2018). To identify genetic susceptibility loci of periodontitis, various genome-wide association studies (GWASs) have been performed. Recently, a GWAS and 2 GWAS meta-analyses independently found SNPs at the inhibitory immune receptor gene sialic acid–binding Ig-like lectin 5 (SIGLEC5) to be associated at a genome-wide significance level (P < 5 × 10^-8) with early-onset generalized stage III, grade C periodontitis (Munz et al. 2017) and more moderate late-onset forms (Munz et al. 2019; Shungin et al. 2019). Additionally, a study that integrated expression quantitative trait locus (eQTL) data with GWAS associations implied SIGLEC5 as periodontitis risk gene (Li et al. 2020). All SNPs of these associations are located in the introns of SIGLEC5. These are regions that do not code for a protein, but variants in introns can change the activity of a gene or cause a protein to be produced in the wrong place or at the wrong time. The closest genes in a distance of the associated SNPs are SIGLEC14 and ZNF175. eQTL effects for rs4284742 and rs11084095 suggested SIGLEC5 as the target gene of the association, with the most significant eQTL effects detected in blood and monocytes, the tissue where SIGLEC5 is almost exclusively expressed (except of placenta), with P = 7.7 × 10^–14 and P = 6.4 × 10^–23, respectively. Reported eQTLs of the associated SNPs were much weaker for other genes, implying SIGLEC5 as the most likely target gene of the association. SIGLEC5 is an inhibitory transmembrane receptor that binds sialic acids and sialic acid–containing glycan ligands, which are expressed on the surfaces of certain pathogens and bacteria. It is discussed that it mediates crosstalk of pathogen-associated molecular patterns (PAMPs) and danger-associated molecular patterns (DAMPs) during infection and wound healing (Pillai et al. 2012).

To elucidate the molecular mechanisms that predispose to increased disease susceptibility, genetic associations need to be leveraged to biological meaning. This poses a challenge because the most significant associated variant, called the GWAS lead SNP or sentinel variant, most often is not identical with the functional variants that caused the association. This is explained by the fact that numerous SNPs are in strong linkage disequilibrium (LD) and coinherited with the GWAS lead SNPs, comprising associated haplotype blocks (Gabriel et al. 2002).

The aims of the current study were to identify the putative causal variants of the associations with SIGLEC5, to characterize the role of the effect alleles in the disease etiology, and to validate SIGLEC5 as the target gene of the association.

Methods and Materials

Protocols for cell culture, transfection, quantitative reverse transcription polymerase chain reaction, isolation of peripheral blood mononuclear cells (PBMCs), and eQTL analysis are given in the appendix.

Selection of Putative Causal Variants

To identify putative causal variants, we integrated data from ENCODE (ENCODE Project Consortium 2012) with the genomic locations of the associated SNP (Appendix). We determined all SNPs in strong LD (r² > 0.8) to the GWAS lead SNPs (Munz et al. 2017; Shungin et al. 2019) in the northwestern European populations CEU (Utah residents with northern and western European ancestry) and GBR (British in England and Scotland; 1000 Genomes Project et al. 2010) of the International Genome Sample Resource using the online tool LDproxy (Machiela and Chanock 2015). Subsequently, sequence to motif alignments for these SNPs were performed with various libraries of binding matrixes for transcription factor binding sites (TFBSs) from Transfac professional (geneXplain), SNPInspector (Genomatix), and the open access database Jaspar. The transcription factor (TF) binding motif was confirmed via the web interface for position weight matrix (PWM) model generation and evaluation PWMTools (Ambrosini et al. 2018).

Electrophoretic Mobility Shift Assay

To determine allele-specific protein-DNA binding, the Gelshift Chemiluminescent EMSA Kit (Activemotif) was used. The nuclear protein extract of Raji cells was extracted as described in the Appendix. For the supershift binding reaction, Raji cells or PBMC nuclear extract (10 µg) and biotin-labeled double-stranded oligonucleotides (20 fmol) were incubated at room temperature for 20 min with 1× binding buffer and 2 µL of specific antibody. To verify the result of DNA-protein interaction, unlabeled oligonucleotides (4 pmol) were added to the binding reaction. The reactions were loaded onto a 5% native polyacrylamide gel and run in 0.5× TBE buffer at 100 V for 1 to 1.5 h. After electric transfer of the products to a nylon membrane, the membrane was cross-linked at 120 mJ/cm²; the biotin-labeled oligonucleotides were detected by chemiluminescence; and the absolute value area of the shifted bands was quantified by ImageJ. All electrophoretic mobility shift assays (EMSAs) were performed in triplicates. Sequences of oligonucleotide probes, antibodies, and analysis details are given in the Appendix.

Luciferase Reporter Gene Assay

Cloning and transfection of the luciferase reporter gene experiments are described in the appendix. All transfections were performed in 3 independent biological replicates.

CRISPR-dCas9 Activation

CRISPR-dCas9 activation (CRISPRa) provides the possibility to test whether a genomic site serves as a cis-regulatory element for a target gene of interest (Simeonov et al. 2017). It allows specific and efficient quantification of the regulatory potential that a chromatin sequence has on gene expression in the endogenous context, including naturally occurring variants. We used CRISPRa to analyze if the DNA elements at rs4284742 and rs11084095 had regulatory effects on SIGLEC5 expression. Details are given in the Appendix.

Results

Assignment of Putative Causal Associated Variants by Integration of LD and Regulatory Chromatin Features

The haplotype block tagged by the GWAS lead SNPs were relatively narrow and encompassed few intronic cosegregating SNPs (Fig. 1A, B). The closest genes in a distance of the associated SNPs are SIGLEC14 and ZNF175 (Appendix Fig. 1). rs4284742 (Munz et al. 2017) had no other variant in strong LD (r² > 0.8), whereas the adjacent associated haplotype block that was tagged by the GWAS meta-analyses’ lead SNPs rs12461706 (Shungin et al. 2019) and rs11084095 (Munz et al. 2019) comprised 5 SNPs in strong LD. These included the 2 sentinel variants and 3 additional tagging SNPs (Table). Of the 6 SNPs that tagged the 2 haplotype blocks, only rs4284742 mapped to regulatory DNA elements as determined by ChIP-Seq and DNAse I hypersensitivity. rs11084095, rs34984145, and rs11880807 located to a region that showed H3K4Me1 methylation in the B-cell line GM12878, a methylation mark that is enriched at active and primed cell type–specific enhancers. Both haplotype blocks were separated from one another by an insulator element (Fig. 1C).

Figure 1.

Linkage and chromosomal positions of the single-nucleotide polymorphisms (SNPs) that showed genome-wide significant associations with increased periodontitis susceptibility. (A) rs12461706 shows strong linkage disequilibrium (LD) with 3 SNPs (r² > 0.8). The plot was generated with rs12461706 as the index SNP with genotype data from the 1000 Genomes Project for populations CEU (Utah residents with northern and western European ancestry) and GBR (British in England and Scotland) in a 500,000–base pair window (LDproxy Tool). The dashed horizontal line indicates r² = 0.8. The SNPs were aligned to their chromosomal positions (x-axis). The blue vertical line shows the position of the LD SNPs at r² > 0.8. rs3829655 missed the LD threshold with r² = 0.77. This SNP is a benign exonic synonymous SNP and was not included in the transcription factor binding site analysis. (B) rs4284742 is not in strong LD with other SNPs (r² > 0.8). (C) rs4284742 locates at chromatin elements that correlate with regulatory functions of gene expression (open chromatin as determined by DNAse I hypersensitivity and transcription factor binding sites experimentally confirmed by ChIP-Seq), as determined in various cell types (data from ENCODE). Chromatin segmentation in B-lymphocyte cells (GM12878) indicated enhancer elements at rs34984145, rs11880807, and rs11084095. The GWAS lead SNPs rs12461706, rs11084095, and rs4284742 and their LD SNPs (r² > 0.8) are highlighted with blue vertical lines. SNP order from left to right: rs34984145, rs11880807, rs12461706, rs11084095, rs4801882, and rs4284742. Because rs11084095 and rs4801882 have a distance of 24 base pairs, they are indicated as 1 line. This figure is available in color online.

Table.

SNPs in Linkage Disequilibrium (r² > 0.8) with the GWAS Lead SNPs rs4284742, rs12461706, and rs11084095.

	Location (GRCh37/hg19)	Allele			r² with rs12461706
SNP ID	Location (GRCh37/hg19)	Common	Minor	MAF^a	r² with rs12461706	GWAS Lead SNP	Predicted TFBS^b
rs4284742	chr19:52131733	G	A	0.25	0.19	Munz et al. (2017)	MAFB (82%)
rs12461706	chr19:52121235	A	T	0.40	1.00	Shungin et al. (2019)	None
rs11084095	chr19:52127030	G	A	0.40	1.00	Munz et al. (2019)	ERG (93%)
rs11880807	chr19:52120522	A	G	0.30	0.81		None
rs34984145	chr19:52120410	A	T	0.43	0.81		BACH2 (83%)
rs4801882	chr19:52127053	G	A	0.44	0.81		None

Bold indicates GWAS-lead SNPs.

GWAS, genome-wide association study; MAF, minor allele frequency; SNP, single-nucleotide polymorphism; TFBS, transcription factor binding site.

Indicated for CEU (Utah residents with northern and western European ancestry). MAFs for different human ethnicities were similar to northwest Europeans.

Position weight matrix similarity score.

We investigated PWM libraries whether the nucleotide variants of the 6 SNPs changed predicted TFBSs. For the common allele of rs11084095, a TFBS for the ETS transcription factor ERG with a matrix similarity of 93% was predicted (Fig. 2A); for the common allele of rs4284742, a TFBS for the TF MAF bZIP transcription factor B (MAFB) with a matrix similarity of 82% was predicted (Fig. 2B); and for the common allele of rs34984145, a TFBS for BACH2 with a matrix similarity of 83% was predicted (Fig. 2C). PWMTools confirmed the high similarity to the TF binding motifs at these SNPs. At the noneffect alleles of the other SNPs, no TFBS was predicted. We selected the 3 SNPs rs4284742, rs11084095, and rs34984145 for subsequent validation of TF binding in vitro.

Figure 2.

rs11084095 and rs4284742 have allele-specific effects on ERG and MAFB binding. (A) The ERG transcription factor (TF) binding motif has a matrix similarity of 93% with the DNA sequence at the common G-allele of rs11084095. (B) The MAFB TF binding motif shows a matrix similarity of 82% with the DNA sequence at the common G-allele of rs4284742. (C) For the common A-allele of rs34984145, a transcription factor binding site was predicted by SNPInspector for the TF BACH2 with a matrix similarity of 83%, whereas the rare T-allele strongly reduces binding affinity.

TFs ERG and MAFB Show Allele-Specific Binding at rs11084095 and rs4284742

To prove protein binding at the 3 selected SNPs and to give evidence for allele-specific binding of the predicted TFs, we performed an EMSA using ERG, MAFB, and BACH2 antibodies with allele-specific DNA probes. SIGLEC5 is mainly expressed in various lymphocytes of the innate immune system and in B cells (Crocker and Varki 2001). We performed the EMSA for rs11084095 with protein extract from PBMCs because in addition to SIGLEC5, ERG is weakly expressed in this cell type, too. MAFB is expressed in whole blood (including B lymphocytes), and BACH2 is expressed in B lymphocytes. Therefore, we performed the EMSA for rs4284742 and rs34984145 with protein extract from a B-lymphocyte cell line (Raji cells). In each EMSA we observed a band supershift with the TF-specific antibody and allele-specific oligonucleotides (Fig. 3A–F). The background of the rare A-allele of rs11084095 significantly reduced ERG binding with P = 0.005 to 99% as compared with binding at the common G-allele (Fig. 3A, B). The background of the rare A-allele of rs4284742 significantly reduced MAFB binding with P = 0.02 to 69% compared to binding at the common G-allele (Fig. 3C, D). At the common A-allele of rs34984145, BACH2 binding showed a significant reduction with P = 0.007 of TF binding by 59% as compared with the rare T-allele (Fig. 3E, F).

Figure 3.

Transcription factor binding site of BACH2 at rs34984145. (A) ERG antibody electrophoretic mobility shift assay (EMSA) performed with rs11084095 allele-specific oligonucleotide probes and nuclear protein extract from peripheral blood mononuclear cells. (B) Absolute value area of the antibody-specific bands. In the background of the rare A-allele of rs11084095, ERG binding to the allele–specific oligonucleotides was 99% reduced compared to the common G-allele (P = 0.005). (C) MAFB antibody EMSA performed with rs4284742 allele-specific oligonucleotide probes and nuclear protein extract from Raji cells. (D) Absolute value area of the antibody-specific bands. In the background of the rare A-allele of rs4284742, MAFB binding to the allele-specific oligonucleotides was 69% reduced (P = 0.02) as compared with the common G-allele. (E) BACH2 antibody EMSA performed with rs34984145 allele-specific oligonucleotide probes and nuclear protein extract from Raji cells. (F) Absolute value area of the antibody-specific bands. In the background of the common A-allele of rs34984145, BACH2 binding to the allele-specific oligonucleotides was reduced by 59% (P = 0.007) as compared to the rare T-allele. Binding of ERG, MAFB, and BACH2 antibodies to allele-specific probes is shown in lanes 2 and 3. Probes with the effect allele abrogated antibody binding as compared with probes with the common G-allele. Unlabeled DNA was added to verify that the band shift was antibody specific in lanes 1 and 4. Values are presented as mean ± SD. *P < .05. **P < .01.

rs4284742, rs11084095, and rs34984145 Are Located in Transcriptional Activators

Regulatory DNA elements can either activate or repress transcription. To measure the activity of the regulatory elements at rs11084095, rs4284742, and rs34984145 and to discriminate their effect directions and allele-specific effect sizes, we employed luciferase reporter gene assays. ERG is expressed in HeLa cells at similar levels as in PBMCs. Because HeLa cells can be transfected more efficiently than PBMCs, which improves the detection of differences in reporter gene expression, we performed the reporter gene experiments in HeLa cells. Here, luciferase activity showed a significant increase in the background of the common G-allele of rs11084095 by 9.9-fold (P = 0.015) as compared with the empty plasmid. The rare A-allele showed no increase in luciferase activity. The difference between alleles was significant with P = 0.013 (Fig. 4A).

Figure 4.

rs11084095 and rs4284742 have allele-specific effects on luciferase activity. The 65-bp DNA sequence up- and downstream of rs11084095 and the 75-bp DNA sequence spanning rs4284742 showed allele-specific enhancer activity in HeLa cells. (A) In the background of the common rs11084095 G-allele, luciferase activity was 9.9-fold increased (P = 0.015) versus the rare A-allele, which showed no upregulation. The difference between alleles was significant with P = 0.013. (B) The common rs4284742 G-allele increased luciferase activity 13-fold (P = 0.008) as compared with the rare A-allele, which increased the activity 4-fold (P = 0.002). The differences between alleles was significant with P = 0.010. (C) Both alleles of rs34984145 showed significant enhancer activity in HeLa cells. The rare T-allele showed stronger upregulation of 4.2-fold (P = 0.0001) versus the common A-allele, which showed an upregulation of 2.3-fold (P = 0.002). The difference between alleles was not significant. (D) The haplotype block at rs4284742 has a cis-regulatory effect on SIGLEC5 expression. CRISPRa of the genomic region at rs4284742 induced SIGLEC5 expression as compared with the scrambled gRNA. gRNA 1 (chr19:52131631-52131649), gRNA 2 (chr19:52131441-52131459), and gRNA 3 (chr19:52131754-52131773) induced SIGLEC5 expression 380-fold (P = 0.00008), 137-fold (P = 0.03), and 80-fold (P = 0.0009), respectively (positions are given for genome build GRCh37/hg19). (E) CRISPRa of the genomic region at rs11084095 did not significantly induce SIGLEC5 expression (gRNA 1: chr19:52126982-52127000, gRNA 2: chr19:52127105-52127123, gRNA 3: 2126912-52126930). Error bars indicate 95% CI. *P < .05. **P < .01. ***P < .001.

The DNA sequence with the common G-allele of rs4284742 showed 13-fold upregulation of luciferase activity (P = 0.008) in HeLa cells as compared with the reporter gene without this regulatory element (Fig. 4B). The reporter gene with the minor A-allele showed an upregulation of 4-fold (P = 0.002). The difference between alleles was significant with P = 0.01.

The luciferase gene activity of rs34984145 showed weak but significant upregulation for both alleles (T-allele: 4.2-fold, P = 0.0001; A-allele: 2.3-fold, P = 0.002), but the fold change difference between alleles was <2 (Fig. 4C), indicating allele-specific effects of low biological relevance.

These experiments indicated that the regulatory elements act as enhancers. Additionally, these experiments showed strong allele differences of rs11084095 and rs4284742 on enhancer activity.

CRISPR-dCas9 Activation of rs4284742 Showed Cis-regulation of SIGLEC5

For rs11084095, an eQTL effect on the expression of SIGLEC5 was reported in monocytes with P = 6.4 × 10^–23 (Zeller et al. 2010), and for rs4284742, an eQTL effect on the expression of SIGLEC5 was noted in whole blood with P = 7.7 × 10^–14 (Blood eQTL browser; Westra et al. 2013), suggesting that the disease-associated elements have cis-regulatory effects on the expression of SIGLEC5 (Appendix Tables 1 and 2). To validate the regulatory potential of the associated haplotype blocks tagged by rs11084095 and rs4284742, we performed CRISPRa in HeLa cells because Raji cells showed poor survival after transfection with the CRISPRa plasmids, probably because of DNA toxicity (Kim et al. 2014). CRISPRa of rs4284742 (gRNAs located in 20- to 300-bp distance from the SNP) with the synergistic activation mediator (SAM) system strongly increased the expression of SIGLEC5 as compared with cotransfected unspecific control sgRNAs (FC = 380, Fig. 4D; Appendix Fig. 2). This proved that the DNA element at rs4284742 exerts a strong cis-regulatory effect on SIGLEC5 and implied that SIGLEC5 is the target gene of the disease-associated genetic variant. CRISPRa of rs11084095 with 3 gRNAs (25- to 120-bp distance from the SNP) showed no significant upregulation of SIGLEC5 (Fig. 4E).

Because the alleles of rs34984145 showed no different effects on luciferase activity and no bloodborne eQTLs of rs34984145 on the expression of SIGLEC5 were reported, we considered that rs34984145 is not a causal variant of the association and did not perform CRISPRa for this SNP.

Discussion

Two large GWASs independently found genetic variants at SIGLEC5 to be associated with different forms of periodontitis, indicating broad relevance of this locus for the disease etiology. We identified rs4284742 and rs11084095 as putative causal variants of the associations with early-onset generalized stage III, grade C periodontitis (Munz et al. 2017) and with more moderate forms (Shungin et al. 2019), respectively. Both SNPs were not in LD (r² = 0.2) and located on different haplotype blocks. However, rs12461706 was in complete LD (r² = 1) with rs11084095, the lead SNP of a GWAS-meta-analysis that combined the generalized stage III, grade C periodontitis cases with a GWAS sample of more moderate periodontitis (Teumer et al. 2013), which was also included in the GWAS meta-analysis of Shungin et al. (2019). To identify the putative causal variants of these associations, we analyzed all SNPs that were in strong LD to the GWAS lead SNPs for locating to predicted TFBSs. The 2 complementary TF databases of Transfac professional and SNPInsepctor allowed prediction of binding sites of all human TFs currently known. However, it is possible that we missed a TF with no currently known binding site. Another possible limitation was that we confined the analysis to SNPs that indicated strong linkage according to r² > 0.8. We used this LD measure because the r² coefficient of correlation takes account of allele frequency. Strong LD indicated by D′ but not by r² would include alleles that are inherited with the particular GWAS lead SNP but are not carried by the majority of cases because they are rare. Such alleles would not be suggestive as causative variants because they would not explain the association for most cases. However, this does not exclude the existence of rare susceptibly variants at SIGLEC5, but such variants have no disease-relevant role for the general population. Instead, their effects would become noticeable in individual cases.

At the position of rs11084095, we identified a binding site for the TF ERG with 93% PWM matrix similarity. We showed that the sequence at rs11084095 is a binding site for ERG and has enhancer activity. The background of the rare effect allele significantly reduced ERG binding and enhancer activity. ERG is an essential TF for endothelial homeostasis, a system control state that encompasses acute responses to injury to support repair of damaged endothelium (Heiss et al. 2015). Changes in endothelial homeostasis (i.e., endothelial injury) affect the vascular structure by interacting with extracellular matrix turnover and influence endothelial membrane function and adhesiveness to proteins of the coagulation cascade and platelets, membrane permeability, and integrity (Gulino-Debrac 2013). Corresponding to function, ERG is mainly expressed in endothelial cells but also in leukocytes (Uhlen et al. 2015). The GWAS catalog (Welter et al. 2014) reports that genetic variations of ERG are associated with numerous blood cell traits and blood pressure, as well as bone mineral density and osteoarthritis. Notably, it was recently shown that soluble SIGLEC5 appears efficient in blocking leukocyte rolling over P-selectin (platelet selectin) and E-selectin (endothelial selectin; Pepin et al. 2016). In human endothelial cells, ERG associates with enhancers of von Willebrand factor (Kalna et al. 2019) and von Willebrand factor may act as a ligand for SIGLEC5 (Pegon et al. 2012), suggesting a possible link between endothelial homeostasis and SIGLEC5 function.

At the position of SNP rs4284742, we identified a predicted TFBS of the TF MAFB and showed that the effect allele, which is the common allele, provided strong MAFB binding affinity as compared with the noneffect allele. Additionally, we showed that CRISPR-dCas9 activation of the genomic sequences at rs4284742 activated SIGLEC5 expression. These experiments indicated that MAFB regulation is linked with activation of SIGLEC5 expression and with increased risk for early-onset periodontitis. MAFB belongs to the subfamily of the large Maf transcription factors (Santos-Gallego 2016). It is expressed by monocytes and is required for differentiation to macrophages (Kelly et al. 2000; Gemelli et al. 2006). Furthermore, it negatively regulates osteoclast generation via inhibition of the transcription factor NFATc1 and the osteoclast-associated receptor (OSCAR; Kim et al. 2007), implying a functional context of the genetic association with osteoclast differentiation. Additionally, MAFB was shown to promote sprouting angiogenesis (Jeong et al. 2017). During healing of tissue injuries, angiogenic capillary sprouts invade the wound clot and organize into a microvascular network throughout the granulation tissue (reviewed by Tonnesen et al. 2000).A challenge for healing of aseptic tissue injuries is discrimination from infections. This is achieved by the innate immune system through danger- or pathogens-associated molecular patterns, DAMPs and PAMPs. Sialoside-based pattern recognition by SIGLEC5 receptors was shown to selectively suppress the immune response to DAMPs, suggesting a mechanism by which aseptic tissue injury and infection are distinguished (Chen et al. 2009). This mechanism is normally well controlled to avoid autoimmune destruction. However, many viruses and pathogenic bacteria express sialidase as virulence factors (Crennell et al. 1993; reviewed by Drzeniek 1972; Crennell et al. 1993), and microbial-expressed sialidases might have potential to abrogate the SIGLEC-mediated inhibitory effects of DAMPs on the innate immune system during healing of injured tissues. As a result, DAMPs and PAMPs would become indistinguishable, which could provide an explanation for massive inflammation and tissue destruction (Liu et al. 2009), as typically seen in severe periodontitis phenotypes (stage III) with a rapid rate of progression (grade C). This may explain why, in specific situations, SIGLEC5 expression may increase the risk for alveolar bone loss.

Notably, the maximal PWM score of MAFB was 82%. We interpret the incomplete similarity score in a way that the specificity of protein-DNA binding depends not solely on the DNA sequence but also on the 3-dimensional structure of DNA and TF protein macromolecules (Rohs et al. 2010). This result in variation of functional binding motifs and, accordingly, the predicative accuracy of a PWM should be interpreted with caution (Weirauch et al. 2013). However, the rare allele further reduced the matrix similarity to 60%, which corresponded with reduced MAFB binding in the EMSA and reduced reporter gene activity. Using CRISPRa, we could show that the enhancer at rs4284742 regulates SIGLEC5 expression. However, we lacked experimental evidence for such effects at rs11084095. This does not signify that SIGLEC5 is not the target gene of this association. Notably, the strongest eQTL of rs11084095 affected SIGLEC5 expression (e.g., P = 6.4 × 10^–23 in monocytes; Zeller et al. 2010) but no other genes. We consider technical problems that impeded the function of the CRISPRa system at this haplotype block—for example, poor targeting or binding of the gRNAs at the selected PAM sequences.

For the common allele of rs34984145, a TFBS for BACH2 was predicted with a PWM matrix similarity of 83%. The antibody-specific EMSA gave evidence for BACH2 binding at the sequence of this SNP, but the luciferase reporter gene showed no differences of enhancer activity between alleles. This is why we considered that rs34984145 is not a causal variant of the association.

We conclude that rs4284742 and rs11084095 are functional variants, that the risk alleles reduce enhancer activity and binding affinity of the TFs MAFB and ERG, and that SIGLEC5 is the target gene of these associations. We suggest a functional context of the associations with impaired regulation of endothelial homeostasis and healing of aseptic tissue injuries.

Author Contributions

R. Mueller, A.S. Schaefer, contributed to conception, design, data acquisition, analysis, and interpretation, drafted and critically revised the manuscript; A. Chopra, contributed to conception, design, data analysis, and interpretation, critically revised the manuscript; H. Dommisch, contributed to data interpretation, critically revised the manuscript. All authors gave final approval and agree to be accountable for all aspects of the work.

Supplemental Material

sj-docx-1-jdr-10.1177_00220345211049984 – Supplemental material for Periodontitis Risk Variants at SIGLEC5 Impair ERG and MAFB Binding

Supplemental material, sj-docx-1-jdr-10.1177_00220345211049984 for Periodontitis Risk Variants at SIGLEC5 Impair ERG and MAFB Binding by R. Mueller, A. Chopra, H. Dommisch and A.S. Schaefer in Journal of Dental Research

Footnotes

A supplemental appendix to this article is available online.

Declaration of Conflicting Interests

The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was funded by the German Research Foundation (SCHA1582 5-1).

ORCID iDs

R. Mueller

A. Chopra

H. Dommisch

References

1000 Genomes Project Consortium; Abecasis

Altshuler

Auton

Brooks

Durbin

Gibbs

Hurles

McVean

. 2010. A map of human genome variation from population-scale sequencing. Nature. 467(7319):1061–1073.

1000 Genomes Project Consortium; Auton

Brooks

Durbin

Garrison

Kang

Korbel

Marchini

McCarthy

McVean

, et al. 2015. A global reference for human genetic variation. Nature. 526(7571):68–74.

Ambrosini

Groux

Bucher

2018. PWMScan: a fast tool for scanning entire genomes with a position-specific weight matrix. Bioinformatics. 34(14):2483–2484.

Boyle

Pritchard

. 2017. An expanded view of complex traits: from polygenic to omnigenic. Cell. 169(7):1177–1186.

Chen

Tang

Zheng

Liu

2009. CD24 and Siglec-10 selectively repress tissue damage-induced immune responses. Science. 323(5922):1722–1725.

Crennell

Garman

Laver

Vimr

Taylor

. 1993. Crystal structure of a bacterial sialidase (from Salmonella typhimurium LT2) shows the same fold as an influenza virus neuraminidase. Proc Natl Acad Sci U S A. 90(21):9852–9856.

Crocker

Varki

2001. Siglecs, sialic acids and innate immunity. Trends Immunol. 22(6):337–342.

Drzeniek

. 1972. Viral and bacterial neuraminidases. Curr Top Microbiol Immunol. 59:35–74.

ENCODE Project Consortium. 2012. An integrated encyclopedia of DNA elements in the human genome. Nature. 489(7414):57–74.

10.

Gabriel

Schaffner

Nguyen

Moore

Roy

Blumenstiel

Higgins

DeFelice

Lochner

Faggart

, et al. 2002. The structure of haplotype blocks in the human genome. Science. 296(5576):2225–2229.

11.

Gemelli

Montanari

Tenedini

Zanocco Marani

Vignudelli

Siena

Zini

Salati

Tagliafico

Manfredini

, et al. 2006. Virally mediated MafB transduction induces the monocyte commitment of human CD34+ hematopoietic stem/progenitor cells. Cell Death Differ. 13(10):1686–1696.

12.

Gulino-Debrac

. 2013. Mechanotransduction at the basis of endothelial barrier function. Tissue Barriers. 1(2):e24180.

13.

Heiss

Rodriguez-Mateos

Kelm

2015. Central role of eNOS in the maintenance of endothelial homeostasis. Antioxid Redox Signal. 22(14):1230–1242.

14.

Jeong

Hernandez-Rodriguez

Kim

Enriquez-Gasca

Yoon

Adams

Scholer

Vaquerizas

Adams

. 2017. Transcriptional regulation of endothelial cell behavior during sprouting angiogenesis. Nat Commun. 8(1):726.

15.

Kalna

Yang

Peghaire

Frudd

Hannah

Shah

Osuna Almagro

Boyle

Gottgens

Ferrer

, et al. 2019. The transcription factor ERG regulates super-enhancers associated with an endothelial-specific gene expression program. Cir Res. 124(9):1337–1349.

16.

Kelly

Englmeier

Lafon

Sieweke

Graf

2000. MafB is an inducer of monocytic differentiation. EMBO J. 19(9):1987–1997.

17.

Kim

Lee

Jin

Kook

Kim

Lee

Kim

2007. MafB negatively regulates RANKL-mediated osteoclast differentiation. Blood. 109(8):3253–3259.

18.

Kim

Cho

Kim

. 2014. Highly efficient RNA-guided genome editing in human cells via delivery of purified Cas9 ribonucleoproteins. Genome Res. 24(6):1012–1019.

19.

Zheng

Meng

Chen

2020. Integration of genome-wide association study and expression quantitative trait loci data identifies AIM2 as a risk gene of periodontitis. J Clin Periodontol. 47(5):583–593.

20.

Liu

Chen

Zheng

2009. CD24-Siglec G/10 discriminates danger- from pathogen-associated molecular patterns. Trends Immunol. 30(12):557–561.

21.

Machiela

Chanock

. 2015. LDlink: a web-based application for exploring population-specific haplotype structure and linking correlated alleles of possible functional variants. Bioinformatics. 31(21):3555–3557.

22.

Munz

Richter

Loos

Jepsen

Divaris

Offenbacher

Teumer

Holtfreter

Kocher

Bruckmann

, et al. 2019. Meta-analysis of genome-wide association studies of aggressive and chronic periodontitis identifies two novel risk loci. Eur J Hum Genet. 27(1):102–113.

23.

Munz

Willenborg

Richter

Jockel-Schneider

Graetz

Staufenbiel

Wellmann

Berger

Krone

Hoffmann

, et al. 2017. A genome-wide association study identifies nucleotide variants at SIGLEC5 and DEFA1A3 as risk loci for periodontitis. Hum Mol Genet. 26(13):2577–2588.

24.

Pegon

Kurdi

Casari

Odouard

Denis

Christophe

Lenting

. 2012. Factor VIII and von Willebrand factor are ligands for the carbohydrate-receptor Siglec-5. Haematologica. 97(12):1855–1863.

25.

Pepin

Mezouar

Pegon

Muczynski

Adam

Bianchini

Bazaa

Proulle

Rupin

Paysant

, et al. 2016. Soluble Siglec-5 associates to PSGL-1 and displays anti-inflammatory activity. Sci Rep. 6:37953.

26.

Pillai

Netravali

Cariappa

Mattoo

2012. Siglecs and immune regulation. Annu Rev Immunol. 30:357–392.

27.

Rohs

Jin

West

Joshi

Honig

Mann

. 2010. Origins of specificity in protein-DNA recognition. Annu Rev Biochem. 79:233–269.

28.

Santos-Gallego CG. 2016. MafB and the role of macrophage apoptosis in atherosclerosis: a time to kill, a time to heal. Atherosclerosis. 252:194–196.

29.

Shungin

Haworth

Divaris

Agler

Kamatani

Keun Lee

Grinde

Hindy

Alaraudanjoki

Pesonen

, et al. 2019. Genome-wide analysis of dental caries and periodontitis combining clinical and self-reported data. Nat Commun. 10(1):2773.

30.

Simeonov

Gowen

Boontanrart

Roth

Gagnon

Mumbach

Satpathy

Lee

Bray

Chan

, et al. 2017. Discovery of stimulation-responsive immune enhancers with crispr activation. Nature. 549(7670):111–115.

31.

Teumer

Holtfreter

Volker

Petersmann

Nauck

Biffar

Volzke

Kroemer

Meisel

Homuth

, et al. 2013. Genome-wide association study of chronic periodontitis in a general German population. J Clin Periodontol. 40(11):977–985.

32.

Timpson

Greenwood

CMT

Soranzo

Lawson

Richards

. 2018. Genetic architecture: the shape of the genetic contribution to human traits and disease. Nat Rev Genet. 19(2):110–124.

33.

Tonnesen

Feng

Clark

. 2000. Angiogenesis in wound healing.J Investig Dermatol Symp Proc. 5(1):40–46.

34.

Uhlen

Fagerberg

Hallstrom

Lindskog

Oksvold

Mardinoglu

Sivertsson

Kampf

Sjostedt

Asplund

, et al. 2015. Proteomics. Tissue-based map of the human proteome. Science. 347(6220):1260419.

35.

Weirauch

Cote

Norel

Annala

Zhao

Riley

Saez-Rodriguez

Cokelaer

Vedenko

Talukder

, et al. 2013. Evaluation of methods for modeling transcription factor sequence specificity. Nat Biotechnol. 31(2):126–134.

36.

Welter

MacArthur

Morales

Burdett

Hall

Junkins

Klemm

Flicek

Manolio

Hindorff

, et al. 2014. The NHGRI GWAS Catalog, a curated resource of SNP-trait associations. Nucleic Acids Res. 42(database issue):D1001–D1006.

37.

Westra

Peters

Esko

Yaghootkar

Schurmann

Kettunen

Christiansen

Fairfax

Schramm

Powell

, et al. 2013. Systematic identification of trans eQTLs as putative drivers of known disease associations. Nat Genet. 45(10):1238–1243.

38.

Wray

Wijmenga

Sullivan

Yang

Visscher

. 2018. Common disease is more complex than implied by the core gene omnigenic model. Cell. 173(7):1573–1580.

39.

Zeller

Wild

Szymczak

Rotival

Schillert

Castagne

Maouche

Germain

Lackner

Rossmann

, et al. 2010. Genetics and beyond—the transcriptome of human monocytes and disease susceptibility. PLoS One. 5(5):e10693.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

1.12 MB