Next Article in Journal
Understanding the Long-Lasting Effects of Fetal Nutrient Restriction versus Exposure to an Obesogenic Diet on Islet-Cell Mass and Function
Next Article in Special Issue
β-Cell Death in Diabetes: Past Discoveries, Present Understanding, and Potential Future Advances
Previous Article in Journal
The Leucine Catabolite and Dietary Supplement β-Hydroxy-β-Methyl Butyrate (HMB) as an Epigenetic Regulator in Muscle Progenitor Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TCF19 Impacts a Network of Inflammatory and DNA Damage Response Genes in the Pancreatic β-Cell

1
Division of Endocrinology, Diabetes, and Metabolism, Department of Medicine, University of Wisconsin-Madison, Madison, WI 53705, USA
2
Department of Biological Sciences, Marquette University, Milwaukee, WI 53233, USA
3
Department of Neuroscience, University of Wisconsin-Madison, Madison, WI 53705, USA
4
William S. Middleton Memorial Veterans Hospital, Madison, WI 53705, USA
*
Author to whom correspondence should be addressed.
Metabolites 2021, 11(8), 513; https://doi.org/10.3390/metabo11080513
Submission received: 18 June 2021 / Revised: 31 July 2021 / Accepted: 2 August 2021 / Published: 4 August 2021
(This article belongs to the Special Issue The Role of β Cells in Diabetes)

Abstract

:
Transcription factor 19 (TCF19) is a gene associated with type 1 diabetes (T1DM) and type 2 diabetes (T2DM) in genome-wide association studies. Prior studies have demonstrated that Tcf19 knockdown impairs β-cell proliferation and increases apoptosis. However, little is known about its role in diabetes pathogenesis or the effects of TCF19 gain-of-function. The aim of this study was to examine the impact of TCF19 overexpression in INS-1 β-cells and human islets on proliferation and gene expression. With TCF19 overexpression, there was an increase in nucleotide incorporation without any change in cell cycle gene expression, alluding to an alternate process of nucleotide incorporation. Analysis of RNA-seq of TCF19 overexpressing cells revealed increased expression of several DNA damage response (DDR) genes, as well as a tightly linked set of genes involved in viral responses, immune system processes, and inflammation. This connectivity between DNA damage and inflammatory gene expression has not been well studied in the β-cell and suggests a novel role for TCF19 in regulating these pathways. Future studies determining how TCF19 may modulate these pathways can provide potential targets for improving β-cell survival.

Graphical Abstract

1. Introduction

The pancreatic β-cells are endocrine cells whose primary role is to synthesize and secrete insulin. Insulin is required to maintain euglycemia. However, the pancreatic β-cell is susceptible to many different stressors including oxidative stress, endoplasmic reticulum (ER) stress, and inflammation [1,2]. These stressors are exacerbated in patients with obesity, insulin resistance, and diabetes [3,4,5]. This can lead to β-cell apoptosis and reduced β-cell mass [5,6]. Pancreatic islets from patients with T2DM have increased ER stress which can lead to β-cell dysfunction and apoptosis [7,8,9]. In addition, increased circulating cytokines and localized islet inflammation are characteristics of T2DM and can contribute to β-cell death [10]. Hyperglycemia, as well as metabolic abnormalities associated with diabetes can lead to oxidative stress, resulting in increased intracellular reactive oxygen species (ROS) that contribute to β-cell dysfunction [11,12]. While many of these sources of β-cell stress have been well studied, there are other factors that can lead to β-cell dysfunction and apoptosis that are less studied. In particular, DNA damage has started gaining attention in recent years as having a role in diabetes pathogenesis. The microenvironment in the islet during diabetes involves oxidative stress and inflammatory insults that can increase DNA damage [13,14,15,16]. Additionally, DNA damage in islets elicited by the β-cell toxin, streptozotocin (STZ), causes an elevation of proinflammatory cytokines [14]. However, this inflammatory response is attenuated after inactivation of the master DNA repair gene, ataxia telangiectasia mutated (ATM) [14]. Horwitz et al. also demonstrated that the β-cell DNA damage response (DDR) was more frequent in islets infiltrated by CD45+ immune cells [14]. This brings to light a fascinating connection between DNA damage, inflammatory, and immune responses in the islet. A better understanding of the intersection between these processes will provide potential regulatory targets to reduce and resolve DNA damage and inflammatory stress on the β-cell that may serve to help maintain adequate β-cell mass and function in diabetes.
In humans, the gene TCF19 (transcription factor 19) is associated with both T1DM and T2DM in genome-wide association studies [17,18,19,20]. TCF19 is expressed in human islets and shows a positive correlation with BMI in nondiabetic subjects [21]. In mice, Tcf19 is widely expressed; however, its expression is highest in the pancreatic islet and increases with obesity when β-cells are known to increase proliferation [21]. Others have similarly identified Tcf19 as a gene upregulated in proliferating β-cells and found that knockdown of TCF19 impairs insulin secretion in a human β-cell line [22,23,24,25,26]. We have previously demonstrated that siRNA-mediated knockdown of Tcf19 in rat insulinoma INS-1 cells reduces β-cell proliferation and survival and impairs cell cycle progression beyond the G1/S checkpoint [21]. Additionally, Tcf19 knockdown increases apoptosis via reduced expression of genes involved in the maintenance of ER homeostasis and increased expression of proapoptotic genes [21].
The TCF19 protein contains a forkhead-associated (FHA) domain, which is a phosphopeptide recognition domain commonly found in many transcription factors that participate in DNA repair and cell cycle regulation [27]. The human TCF19 (hsTCF19) protein, but not the mouse protein, also harbors a plant homeodomain (PHD) finger, allowing it to interact with chromatin. PHD finger proteins are often considered “chromatin readers” that recognize modified histones and can recruit additional transcriptional machinery to these areas [28]. Specifically, the tryptophan residue at position 316 in hsTCF19 has been shown to bind to chromatin via tri-methylated histone H3 and to regulate cell proliferation in liver cells through this interaction [29,30]. Taken together, these characteristics support the role of TCF19 as a transcriptional regulator of β-cell proliferation and survival.
The purpose of this study was to determine the effect of TCF19 overexpression on proliferation and survival in the β-cell. In this study, we overexpressed TCF19 in INS-1 cells and found that TCF19 overexpression does not induce proliferation, cell cycle progression, or impact cell survival. Rather, there was significant upregulation of a tightly interconnected set of genes involved in cell stress, inflammation, and antiviral responses, alluding to a previously unexplored role for TCF19 in the β-cell. Additionally, we find that TCF19 overexpression in human islets leads to significant upregulation of several DDR genes. Using a novel analysis for potential transcriptional co-regulators on these upregulated genes, we identified STAT1, STAT2, and IRF1 as likely drivers of the tight transcriptional gene network. Interestingly, there was no measurable activation of these transcription factors, indicating alternate mechanisms of regulating the inflammatory and DDR gene expression. These findings not only identify an intriguing connection between DNA damage and inflammatory responses in the β-cell but elucidate a novel role for TCF19 in modulating these two pathways.

2. Results

2.1. Human TCF19 Overexpression Increases 3H-Thymidine Incorporation in INS-1 Cells but Does Not Change Cell Cycle Gene Expression

Based on our original studies on Tcf19, we concluded that Tcf19 was necessary for normal β-cell proliferation, as Tcf19 knockdown led to impaired cell cycle progression, reduced 3H-thymidine incorporation, and G1/S cell cycle arrest [21]. We next wanted to determine if increased levels of TCF19 could drive β-cell proliferation, and therefore, we overexpressed hsTCF19 in INS-1 rat insulinoma cells. The human TCF19 protein was chosen for overexpression as it contains the PHD finger domain that is known to mediate interactions with methylated histones (specifically trimethylated histone 3 at lysine 4 (H3K4me3)) [29,30]. The PHD finger domain is not present in the rodent protein. As we have not yet identified a reliable and specific TCF19 antibody, we generated a C-terminal myc-tagged TCF19 to allow for probing on the western blot. TCF19 overexpression was confirmed at both the mRNA and protein level (Figure 1A,B).
As an assay to assess proliferation, we measured 3H-thymidine nucleotide incorporation in cells expressing hsTCF19 vs. empty vector control. INS-1 cells overexpressing hsTCF19 showed a significant two-fold increase in 3H-thymidine nucleotide incorporation suggesting increased cell proliferation (Figure 1C). To confirm that the 3H-thymidine nucleotide incorporation observed correlated with an increase in the expression of cell cycle genes, as would be expected in a dividing cell, we assessed cell cycle gene expression with quantitative real-time PCR (qRT-PCR). Interestingly, there was no significant change in expression of cell cycle genes, including the proliferative marker, Ki67 (Figure 1D). We concluded that overexpression of hsTCF19 in INS-1 cells does not lead to transcriptional activation of cell cycle genes, suggesting an alternate process for nucleotide incorporation that does not result in cell cycle progression. DNA repair may be an alternative pathway that leads to increased 3H-thymidine nucleotide incorporation [31]. DNA damage and repair responses are important in preserving genome integrity, and an accumulation of DNA damage without sufficient repair can result in cell cycle arrest at the G1/S checkpoint [32]. However, qRT-PCR showed no significant change in cell cycle inhibitors Cdkn2c (p18), Cdkn1a (p21), and Cdkn1b (p27) with hsTCF19 overexpression, suggesting that there was no induction of substantial DNA damage leading to cell cycle arrest or activation of checkpoint inhibitors (Figure 1D). We next hypothesized that if hsTCF19 overexpression is affecting DNA repair, it may elicit a change in cell viability. However, after staining cells with trypan blue, we found that the percentage of live cells was not significantly affected by hsTCF19 overexpression (Figure 1E).

2.2. RNA-Seq Analysis Reveals a Role for TCF19 in Regulating Viral, Inflammatory, and DNA Damage Genes

To obtain a more global perspective on what genes TCF19 could be regulating, we performed RNA-seq analysis on INS-1 cells overexpressing hsTCF19. Notably, this revealed only a relatively small number of differentially expressed genes. Of the 160 genes differentially expressed between the groups (false discovery rate (FDR) < 5%), 136 genes were upregulated and 24 were downregulated (Table S1), suggesting that TCF19 likely acts as a positive regulator of transcription. 92 gene IDs were identified as being upregulated >2-fold from the original list of 136 upregulated genes, and of these, 85 were uniquely mapped and included in the PANTHER Fisher’s Exact overrepresentation test [33,34,35]. This analysis revealed 199 significantly overrepresented Biological Process Gene Ontology (GO) terms (FDR < 5%). These GO terms were sorted hierarchically and the most specific subclasses of GO terms are listed in Figure 2A. This analysis shows that genes upregulated with hsTCF19 overexpression are highly enriched for biological processes relating to double-strand break repair, apoptosis in response to ER stress, antigen presentation, interferon signaling, immune system processes, and viral responses (Figure 2A).
To determine the relationship between the significantly upregulated genes, we performed Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) analysis on the 92 genes that were upregulated >2-fold. STRING uses an algorithm to reflect co-occurrences of genes in literature to predict associations for a particular group of genes [36]. STRING analysis revealed highly significant connectivity between almost all input genes (enrichment p-value < 1 × 10−16), suggesting that TCF19 may be regulating one cluster of interconnected genes (Figure 2B). We hypothesized that this cluster of genes may have roles in viral and interferon responses, as well as the DDR.
Among the upregulated genes (Table S1), several are known to be involved in DDR and repair pathways (Parp9, Parp10, Parp12, Parp14) [37]. In particular, Parp9 and another gene from the dataset, Dtx3l, have been shown to work as a complex to promote DNA repair [38]. Other significantly upregulated genes include those from the oligoadenylate synthase (Oas) family (Oas1i, Oasl2, Oas2, Oas1a, Oas1g, Oas1f), which are stimulated by type 1 interferons in response to viral infections [39]. However, they can also be activated by DNA damage, where they may have roles in Poly(ADP-ribose) (PAR) synthesis and interacting with PARP1 during DNA repair [39,40]. Mx1, Ddx60, and Usp18 are genes with known antiviral roles and were also significantly upregulated [41,42,43]. These observations suggest that TCF19 may play a previously unreported role in the DDR and viral and inflammatory response pathways.
To assess the extent to which the findings in this overexpression model could be translated to human islets, we overexpressed hsTCF19 in human islets (Figure 3A) and assessed several of the differentially expressed genes from the RNA-seq dataset (Figure 3B). Notably, transcript levels for DDR genes PARP9 and DTX3L were significantly upregulated in human islets overexpressing hsTCF19 compared to the empty vector control islets (Figure 3B). Antiviral genes MX1 and DDX60 were also significantly upregulated (Figure 3B).

2.3. Mining Algorithm for GenetIc Controllers (MAGIC) Analysis for Common Transcriptional Regulators

To look for common transcriptional regulators associated with the promoters of the 92 genes that were upregulated >2-fold in INS-1 cells overexpressing hsTCF19, we performed MAGIC analysis [44]. MAGIC compares input gene lists to ChIP seq tracks archived in Encyclopedia of DNA Elements (ENCODE) and predicts which nuclear proteins are enriched at the promoters/regulatory regions of the input gene list [44]. These analyses revealed Signal Transducer and Activator of Transcription (STAT)1 and STAT2 as positive drivers of the gene set (Figure 4). The associations were striking with p-values of 7.81 × 10−19 and 3.23 × 10−20, respectively. Specifically, STAT1 and STAT2 are known to interact with the promoter of 17 genes out of the 92 from the gene set. There was strong enrichment within the gene set compared to overall promoter interactions for these STAT proteins across the genome, suggesting that TCF19 leads to upregulation of genes that can also be regulated by the STAT proteins. However, these associations in ENCODE were not determined in β-cells or islets and were often based on experiments involving interferon stimulation. Additionally, Interferon Response Factor (IRF)1, a transcription factor important in both innate and adaptive immunity, also showed striking enrichment for promoter interactions with the upregulated gene set (p = 2.77 × 10−16).
As these transcription factors could be potential regulators of the upregulated genes in this dataset, we assessed the activation level of these transcription factors. Full activation of STAT1 involves both phosphorylation at tyrosine (Y701) and serine (S727) residues [45]. Y701 phosphorylation is required for nuclear accumulation of STAT1, while full transcriptional activity of STAT1 requires phosphorylation at S727 [45]. STAT1 is phosphorylated at the S727 residue in response to interferons [45,46]. We did not observe a significant difference in S727 levels (Figure 5A). Interestingly, we observed a significant decrease in Y701 phosphorylation (Figure 5B). There was no change in IRF1 levels (Figure 5C). STAT2 was not detectable in the INS-1 cells. Taken together, this suggests that TCF19 does not directly modulate the levels of STAT1, STAT2, or IRF1 in β-cells, nor does it impact the phosphorylation of S727 on STAT1 that is required for full activity. We actually see a decrease in Y701 STAT1 phosphorylation, suggesting that less STAT1 is capable of moving to the nucleus for transcriptional activity in the presence of TCF19.
Although not identified as a potential co-regulator in MAGIC analysis, Nuclear Factor Kappa-B (NF-κB) has a well-characterized role in mediating inflammation and is also activated by the cyclic GMP-AMP synthase (cGAS)-Stimulator of Interferon Genes (STING) pathway, which is a component of the innate immune system that functions to detect cytosolic DNA and leads to the production of type 1 interferons [47]. Additionally, positive regulation of the NF-κB pathway signaling was a biological process that was significantly overrepresented in PANTHER GO analysis (Figure 2A). Therefore, we predicted that NF-κB may be a possible regulator of the upregulated gene set. However, we found no increase in the phosphorylation of NF-κB with TCF19 overexpression (Figure 5D).

3. Discussion

Inflammation is a pathophysiological state associated with both T1DM and T2DM. In T1DM, immune cells are critical mediators of islet inflammation through their secretion of cytokines such as interleukin 1 beta (IL-1beta) and tumor necrosis factor-alpha (TNF-alpha) [48]. Additionally, substantial evidence suggests that triggering events such as a viral infection may initiate the β-cell damaging process [49]. In T2DM, obesity induces chronic, low grade inflammation which activates inflammatory pathways and the release of pro-inflammatory cytokines and adipokines [50,51]. Inflammation not only exacerbates insulin resistance and promotes β-cell death but can also contribute to DNA damage [15,52].
In this study, we delineate a role for TCF19, in the inflammatory and DNA damage pathways. We find that TCF19 overexpression significantly increases expression of inflammatory and DDR genes, suggesting a novel role for TCF19 in regulating these two pathways. We find that the significantly upregulated genes from TCF19 overexpression are tightly associated, and we describe potential transcription factor co-regulators of these genes. This brings to light an interesting crosstalk between the inflammatory and DNA damage pathways in the β-cell and suggests how alterations in TCF19 expression or function may contribute to diabetes pathogenesis in both T1DM and T2DM. Although one of many genes associated with diabetes in GWAS, TCF19 is unusual in having associations with both types of diabetes. This shared association suggests that TCF19 may regulate a mechanism involving shared pathophysiology in both T1DM and T2DM, such as β-cell damage or inflammatory responses.
Knockdown of Tcf19 has been shown to result in cell cycle arrest [21]. While we show here that overexpression of hsTCF19 does not result in significant changes in cell cycle genes, hsTCF19 overexpression does result in increased expression of DDR genes. The DDR is made up of DNA damage sensing proteins, transducers, and effectors [53]. Once an aberrant DNA structure is recognized, downstream phosphorylation cascades within the DDR network are initiated with many of the downstream effector proteins having roles in promoting cell cycle arrest [53]. This allows time for the cell to repair the damaged DNA. Other effector proteins upregulate DNA damage repair genes or promote senescence or apoptosis in the face of unrepairable DNA damage [32]. With TCF19 overexpression, we find an increase in genes involved in the DDR but no decrease in cell viability, suggesting that these cells are not undergoing apoptosis. Additionally, the lack of significant change in cell cycle genes including cell cycle inhibitors suggests there is no DNA damage-induced cell-cycle arrest. Notably, these experiments were all performed in the absence of any inducers of DNA damage or interferons, yet we observed upregulation of classic interferon-response genes. Therefore, enhanced TCF19 expression alone is sufficient to independently activate these pathways. Since overexpression of TCF19 led to upregulation of many DNA damage repair genes, this suggests that within the DDR network TCF19 most likely plays a role as a transcriptional regulator that may promote DNA damage repair.
Tcf19 knockdown leads to cell cycle arrest at the G1/S transition [21]. This is consistent with the cell cycle arrest that occurs upon DNA damage to cells in the G1 phase to prevent entry into the S phase [32]. Sustained DNA damage can eventually result in cellular apoptosis [54]. We previously showed that 3–7 days of Tcf19 knockdown led to an increase in cells undergoing apoptosis and a decrease in cell viability. Combining these prior results with current data, we propose that cells lacking TCF19 are inefficient at repairing DNA damage, ultimately leading to cell cycle arrest or cell death due to accumulated DNA damage. We hypothesize that with TCF19 overexpression, DNA damage repair is upregulated.
Interestingly, many of the genes upregulated by TCF19 overexpression are also involved in interferon and immune responses. Additionally, GO analysis revealed an overrepresentation of genes involved in viral, DNA damage, and stress response processes. This signature of viral, inflammatory, and DNA damage responses brings to light an interesting and emerging field connecting DNA damage and the interferon response. Treatment of cells with etoposide, an agent that induces double stranded DNA breaks, leads to the induction of interferon-stimulated genes regulated by NF-κB [55]. The cGAS-STING pathway is a component of the innate immune system that functions to detect cytosolic DNA and, upon activation of STING, results in the production of type 1 interferons [56]. However, after etoposide treatment, there is noncanonical activation of the STING pathway by the DNA repair proteins, ATM and PARP1 [56]. Additionally, the DNA sensor, cGAS, has been shown to be shuttled to the nucleus under conditions of DNA damage [57]. Given these connections, we also looked for an increase in phospho-STING after TCF19 overexpression but did not see any significant changes (data not shown). Further exploration of possible connections between the cGAS-STING pathway and DNA damage and inflammatory responses in the β-cell remain intriguing new directions for future study.
While these studies show that DNA damage can lead to inflammatory gene expression, inflammation can also induce DNA damage. Chronic inflammation can lead to the production of ROS, which are capable of DNA damage through the formation of free radicals and DNA lesions [58]. In further support of the coordinate regulation between these two pathways, viruses can activate the DDR network and also inhibit several DDR proteins [59]. As viral infection is an important initiating factor in T1DM, this could serve as a potential link between the two pathways where the immune system’s viral response may trigger DNA damage and progression to T1DM. It is likely that the DDR and inflammatory pathways are part of a positive feedback loop [60]. We see a dual response gene signature of viral/interferon and DNA damage processes with TCF19 overexpression, suggesting that TCF19 may regulate both these processes. However, we acknowledge the possibility that TCF19 may regulate just one of these processes, and in turn, be indirectly affecting the other.
STRING analysis further supports the tight association between the DNA damage and inflammatory genes in this dataset. MAGIC analysis revealed STAT1, STAT2, and IRF1 as common regulators of this gene set. These transcription factors have well-characterized roles in response to interleukins and interferons, specifically type 1 interferons [61]. However, there have also been studies showing a role for these transcription factors in the DDR and repair pathway [59]. A few of these transcription factors have been found to be responsible for the induction of interferon alpha and gamma genes in response to DNA damage or have roles in regulating DNA damage repair proteins [59,62].
While we did not observe direct increases in phosphorylation or protein levels of these transcription factors, phosphorylation events can be transient and tightly regulated. It is possible that the time point of harvest (48 h post-transfection) may have been too late to capture the phosphorylation event. We did observe a significant decrease of STAT1 Y701 phosphorylation which is required for STAT1 dimerization, nuclear translocation, and DNA binding [45]. This alludes to the possibility that TCF19 may actually be attenuating STAT1 activity. However, there has been debate that dimerization of STAT1 may not be necessary for initiation of interferon-dependent signaling, and that positive and negative transcriptional control may also be modulated by unphosphorylated STAT1 [63]. Additionally, studies have shown that S727 of STAT1 can be phosphorylated independently of Y701, and that Y701 is necessary but not sufficient for interferon-induced S727 phosphorylation [45]. The time scale by which each phosphorylation event reaches its maximal activity is different, with Y701 reaching its maximal level earlier than S727, after cytokine stimulation [45]. This demonstrates the complexity of STAT1 regulation and bodes for further investigation as to how TCF19 overexpression may be modulating its different phosphorylation status.
While we chose to look at phosphorylation events for activation of the transcription factors identified through MAGIC analysis, other types of post-translational modifications, such as those that may work to alter chromatin structure or recruit histone modifiers cannot be ruled out. Notably, TCF19 has been shown to interact with H3K4me3 through its PHD finger to repress gluconeogenic gene expression and to modulate proliferation in HepG2 cells [29,30]. Therefore, it is likely that TCF19 is not directly activating these transcription factors through phosphorylation events, but instead may bind to H3K4me3 at a transcriptionally active promoter and thereby impact transcriptional activation. Additionally, the TCF19 protein harbors an FHA domain, which may allow binding to phosphor-epitopes on proteins [64,65]. FHA domains are often found in proteins that are critical in the cell cycle and regulated through phosphorylation events but are also found in proteins that are involved in the DDR [64]. The FHA domain of TCF19 contains a serine residue at position 78 (Ser78) that has been shown to be phosphorylated after DNA damage [66]. Ser78 in TCF19 is located within a Ser-Gln motif, which is recognized by kinases involved in the DDR such as ATM and ATR [66]. Therefore, it is possible that TCF19 is a downstream target of ATM or can alter gene expression by acting as a co-regulator to other kinases.
We hypothesize that TCF19 affects DDR gene expression through interactions with modified histones via the PHD finger and/or acts as a co-activator to DNA damage proteins by recruiting other DNA damage transcription factors to areas of active chromatin. While we did not directly measure an interaction of any of the transcription factors from the MAGIC analysis with relevant promoters in response to TCF19 overexpression, our data suggest that TCF19 either modulates their ability to activate transcription or may in fact simply be regulating the expression of these genes independently of these transcription factors. The exact mechanism of how TCF19 modulates these inflammatory and DDR genes to promote diabetes susceptibility requires further investigation.
Overall, our work highlights the complexity of regulation of gene expression involved in DNA damage and inflammatory response genes and alludes to the interesting crosstalk between these processes in the context of TCF19. With respect to diabetes susceptibility, individuals with genetic variants of TCF19 may be unable to properly regulate β-cell responses to DNA damage and inflammatory insults, therefore predisposing them to increased β-cell apoptosis. Future experiments will explore the nature by which TCF19 modulates DNA damage repair and inflammatory genes under conditions of stress. Furthermore, this will provide for potential therapeutic targets to prevent or attenuate DNA damage and inflammation to preserve functioning β-cells in at-risk individuals.

4. Materials and Methods

4.1. Human Islets and INS-1 Cell Culture

INS-1E rat insulinoma cells were cultured in RPMI 1640 supplemented with 1% antibiotic-antimycotic (Gibco, 15240–062), 1% L-glutamine, 1% sodium pyruvate, and 10% fetal bovine serum. 2-Mercaptoethanol was added to a final concentration of 50 μM to supplemented media before each use. Human islets were obtained from nondiabetic organ donors through the Integrated Islet Distribution Program. An exemption was granted for human islet work by the Institutional Review Board at the University of Wisconsin. Human islets were cultured in uncoated petri dishes with RPMI 1640 containing 8mM glucose, 10% heat-inactivated fetal bovine serum, and 1% penicillin/streptomycin. INS-1 cells and islets were cultured at 37 °C and 5% CO2 in a humidified atmosphere.

4.2. Creation of TCF19 Overexpression Vector

The human TCF19 clone HsCD00002769 was purchased in the pDNR-Dual vector backbone (DNASU Plasmid repository). The pcDNA4-TO-myc/his B backbone vector (Invitrogen) was chosen for overexpression. This vector utilizes a CMV promoter, which ensures robust expression of the inserted gene of interest. Following the inserted TCF19 sequence is both a C-terminal c-myc tag as well as six histidine residues to allow for identification of the overexpressed protein in the absence of reliable TCF19 antibodies. The hsTCF19-pcDNA4 vector was created with In-Fusion HD cloning (Clontech) following kit instructions. Colonies were screened with PCR for insert size and then sequenced to confirm TCF19 insertions and sequence integrity.

4.3. Transfection with hsTCF19-His/Myc-pcDNA4

INS-1 cells and islets were transfected with either hsTCF19 or pcDNA4 control, using Lipofectamine 2000 (Invitrogen, Waltham, MA, USA). INS-1 cells were trypsinized and resuspended in transfection medium (RPMI 1640 supplemented with 1% L-glutamine, 1% sodium pyruvate, and 10% fetal bovine serum). Cells in transfection medium were then added to an hsTCF19 or control plasmid-Lipofectamine mixture at 2–5 µg DNA/5 × 106 cells and plated. Transfection medium was removed 12–18 h post-transfection and replaced with complete growth medium. These conditions were the same for all INS-1 overexpression studies, including RNA-Seq sample preparation.
Human islets were washed in 1x PBS and resuspended in Accutase (Sigma, St. Louis, MO, USA) dissociation solution for 3 min at 37 °C, with tube inversions every 30 s. Islets were then resuspended in 2 mL transfection medium and plated into dishes. hsTCF19 or control plasmid-Lipofectamine 2000 mixture was added at 2 µg DNA/1000 islets. Transfection medium was removed 12–18 h post-transfection and replaced with complete growth media.

4.4. Western Blotting

INS-1 cells were harvested 48 h after transfection and washed in ice-cold PBS. Cells were lysed in protein lysis buffer (0.05 M HEPES, 1% NP-40, 2 mM activated sodium orthovanadate, 0.1 M sodium fluoride, 0.01 M sodium pyrophosphate, 4 mM PMSF, 1 mM leupeptin, 2 μM okadaic acid and Sigma Protease inhibitor cocktail). Cells were incubated in the lysis buffer on ice for 15 min with vortexing every 5 min. The protein concentrations were determined using Bradford protein assay. The protein samples were run on 4–10% SDS-PAGE gradient gel and transferred to polyvinylidene difluoride (PVDF) membrane. Membranes were blocked in 5% milk in Tris-buffered saline with 0.1% Tween 20 (TBST) for 1 h at room temperature and were incubated overnight in primary antibody, washed 3X in TBST and incubated 1 h in secondary antibody. Blots were developed with Pierce ECL Western Blotting Substrate (Thermo Fisher, Waltham, MA, USA) or Supersignal West Femto Maximum Sensitivity Substrate (Thermo Fisher), imaged with a GE ImageQuant charge-coupled device camera, and then quantified by densitometry with Image J 1.44o (http://imagej.nih.gov/ij, accessed 28 February 2020). Primary antibodies and dilutions were as follows: Myc antibody (9E10:sc-40, Santa Cruz Biotechnology, 1:1000), Beta actin (8H10D10, Cell Signaling, 1:1000), STAT1 (#9172, Cell Signaling, 1:1000), phospho- STAT1 Ser727 (#9177, Cell Signaling, 1:1000), phospho-STAT1 Tyr701 (D4A7 #7649, Cell Signaling, 1:1000), phospho- NF-κB p65 Ser536 (93H1 #3033, Cell Signaling, 1:1000), NF-κB p65 (D14E12 XP #8242, Cell Signaling, 1:1000), IRF1 (D5E4, Cell Signaling, 1:1000), GAPDH (14C10, Cell Signaling, 1:1000), Beta-tubulin (#2146, Cell Signaling, 1:1000) in 5% BSA-TBST or 5% non-fat milk. Secondary antibodies and dilutions were as follows: Anti-rabbit IgG, HRP-linked antibody (#7074, Cell Signaling, 1:2000), anti-mouse IgG, HRP-linked antibody (#7076, Cell Signaling, 1:2000).

4.5. Quantitative Real-Time PCR

RNA was isolated from INS-1 and human islets 48 h post-transfection using RNeasy cleanup kit (Qiagen, Hilden, Germany) according to manufacturer’s protocol. Concentration and purity of RNA were determined using a NanoDrop ND-2000c Spectrophotometer, and 100–250 ng of RNA was reverse transcribed to make cDNA with Applied Biosystem High-Capacity cDNA synthesis kit. Quantitative real-time PCR (qRT-PCR) reactions were carried out using Power SYBR green PCR Master Mix (Applied Biosystems) and the StepOnePlus Real Time PCR System (Applied Biosystems). Reverse transcriptase free samples were used as negative controls. All samples were run in triplicates with Cycle threshold (Ct) values normalized to β-actin to yield ∆Ct. Fold changes were then calculated between experimental and control samples: fold change 2(∆Cexperimental—∆Ctcontrolt). For gene expression in INS-1 cells, results were analyzed by non-paired t-test of the ∆Ct t values, while human islets were analyzed by paired t-test. Significance was determined by p < 0.05. Primer sequences used are in Table S2.

4.6. Viability

Transiently transfected INS-1 cells were harvested at 48 h post-transfection by using a cell scraper to dislodge all cells, and 10 μL of cells were collected from each well. Cell viability was determined using trypan blue (Corning, Corning, NY, USA) staining using the TC-10 Automated Cell Counter (BioRad, Hercules, CA, USA). Comparisons were made by paired t-test, including all technical and biological replicates; statistical significance was determined by p < 0.05.

4.7. Proliferation/3H-Thymidine Incorporation

To measure cell proliferation, transiently transfected INS-1 cells were incubated with 3H-thymidine (Perkin Elmer, Waltham, MA, USA, NET0270001MC) at a final concentration of 1 µL 3H-thymidine/mL of supplemented RPMI media for 4 h. Cells were then trypsinized and washed three times with ice-cold PBS. DNA and protein were precipitated by the addition of ice-cold 10% trichloroacetic acid (TCA) and incubated for 30 min on ice. The precipitate was then pelleted at 18,000× g for 10 min at 4 °C. Pelleted precipitate was solubilized in 0.3 N NaOH and vortexed for 15 min. Radioactivity was measured using a liquid scintillation counter, and a fraction of the solubilized product was kept to measure total protein by the Bradford assay. Sample counts per minute were individually normalized to protein, and an average for each transfection was determined. Results were analyzed by unpaired t-test, and statistical significance was determined by p < 0.05.

4.8. RNA Sequencing

INS-1 cells were transfected with either hsTCF19-pcDNA4 or pcDNA4 control vector as stated in the methods above. Cells were cultured 48 h post-transfection before being collected for RNA using the RNeasy Kit (Qiagen). Total RNA was verified for concentration and purity using a NanoDrop ND-2000c Spectrophotometer and Agilent 2100 BioAnalyzer. Samples that met the Illumina TruSeq Stranded Total RNA (Human/Mouse/Rat) (Illumina Inc., San Diego, CA, USA) sample input guidelines were prepared according to the kit’s protocol. Cytoplasmic ribosomal RNA reduction of each sample was accomplished by using complementary DNA probe sequences attached to paramagnetic beads. Subsequently, each mRNA sample was fragmented using divalent cations under elevated temperature, and purified with Agencourt RNA Clean Beads (Beckman Coulter, Pasadena, CA, USA). First strand cDNA synthesis was performed using SuperScript II Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) and random primers. Second strand cDNAs were synthesized using DNA Polymerase I and RNAse H for removal of mRNA. Double-stranded cDNA was purified using Agencourt AMPure XP beads (Beckman Coulter, Pasadena, CA, USA). cDNAs were end-repaired by T4 DNA polymerase and Klenow DNA Polymerase and phosphorylated by T4 polynucleotide kinase. The blunt ended cDNA was purified using Agencourt AMPure XP beads. The cDNA products were incubated with Klenow DNA Polymerase to add an ‘A’ base (Adenine) to the 3′ end of the blunt phosphorylated DNA fragments and then purified using Agencourt AMPure XP beads. DNA fragments are ligated to Illumina adapters, which have a single ‘T’ base (Thymine) overhang at their 3′end. The adapter-ligated products are purified using Agencourt AMPure XP beads. Adapter ligated DNA was amplified in a Linker Mediated PCR reaction (LM-PCR) for 12 cycles using Phusion™ DNA Polymerase and Illumina’s PE genomic DNA primer set followed by purification using Agencourt AMPure XP beads. Quality and quantity of finished libraries were assessed using an Agilent DNA1000 series chip assay (Agilent Technologies, Santa Clara, CA, USA) and Invitrogen Qubit HS cDNA Kit (Invitrogen, Carlsbad, CA, USA), respectively. Libraries were standardized to 2 nM. Cluster generation was performed using the Illumina cBot. Paired-end, 100bp sequencing was performed, using standard SBS chemistry on an Illumina HiSeq2500 sequencer. Images were analyzed using the standard Illumina Pipeline, version 1.8.2. RNA Library preparation and RNA Sequencing was performed by the University of Wisconsin-Madison Biotechnology Center.
Sequencing reads were adapted, and quality trimmed using the Skewer trimming program [67]. Quality reads were subsequently aligned to the annotated reference genome (Rnor_6.0) using the STAR aligner [68]. Quantification of expression for each gene was calculated by RSEM [69]. The expected read counts from RSEM were filtered for low/empty values and used for differential gene expression analysis using EdgeR [70] using a FDR cut off of <5%.

4.9. GO Term Enrichment

Upregulated genes with fold change >2 were input into PANTHER 16.0 and compared against Rattus norvegicus reference genome list using the Overrepresentation Test to test for enrichment using the GO biological process complete Ontology database DOI:10.5281/zenodo.4735677 (Released 2021-05-01, http://www.pantherdb.org, accessed 8 June 2021) [33,34,35]. Fisher’s exact test and FDR correction was used to determine statistical significance. We had 92 gene IDs on the initial input list, and 85 were uniquely mapped to their corresponding PANTHER ID.

4.10. Protein-Protein Interaction Network Construction

Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database was used to construct the protein-protein interaction network on genes that were upregulated >2 fold (https://string-db.org/, access date 8 June 2021) [36]. Interaction score of >0.4 was used as the cutoff criterion.

4.11. Mining Algorithm for GenetIc Controllers (MAGIC) Analysis

MAGIC analysis uses Encyclopedia of DNA Elements (ENCODE) ChIPseq data to look for statistical enrichment of transcription factors (TFs) that are predicted to bind to regions in a gene set. It determines if genes in a list are associated with higher ChIP values than expected by chance for a given transcription factor or cofactor based on ENCODE data. Detailed methods are found in Roopra et al. [44] All genes that were induced more than 2-fold with an associated FDR < 5% were used as input and tested against the 5Kb_Gene.mtx matrix.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/metabo11080513/s1, Table S1: Differentially expressed genes with FDR < 5% from RNA-seq on hsTCF19 overexpression in INS-1 cells. Table S2: Primer sequences used for quantitative real-time PCR experiment.

Author Contributions

Conceptualization, G.H.Y., D.A.F., S.L. and D.B.D.; Data curation, G.H.Y., D.A.F., S.L. and J.T.B.; Formal analysis, G.H.Y., D.A.F., S.L., J.T.B. and A.R.; Funding acquisition, D.B.D.; Methodology, G.H.Y., D.A.F., S.L., J.T.B. and A.R.; Writing—original draft, G.H.Y., D.A.F. and S.L.; Writing—review & editing, G.H.Y., J.T.B., A.R. and D.B.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Department of Veterans Affairs Merit Awards I01BX001880 and I01BX004715 to DBD. DAF was supported by NIA T32AG000213 and the University of Wisconsin SciMed program. GHY was supported by NRSA Training Core TL1 TR002375 of the NCATS grant UL11TR002373. JTB received support from NIDDK T32 DK007665. Human pancreatic islets were provided by the NIDDK-funded Integrated Islet Distribution Program (IIDP) 2UC4DK098085.

Institutional Review Board Statement

Not applicable. An exemption was granted for human islet work by the Institutional Review Board at the University of Wisconsin.

Informed Consent Statement

Not applicable.

Data Availability Statement

The Python script and information on how to access required Matrix files are available at https://github.com/asroopra/MAGIC.

Acknowledgments

We thank Matthew Flowers for critical review of this manuscript. The authors utilized the University of Wisconsin—Madison Biotechnology Center’s Gene Expression Center Core Facility (Research Resource Identifier—RRID:SCR_017757) for RNA library preparation, the DNA Sequencing Facility (RRID:SCR_017759) for sequencing, and the Bioinformatics Resource Center (RRID:SCR_017799) for RNASeq quality control and analysis. This work was performed with facilities and resources from the William S. Middleton Memorial Veterans Hospital. This work does not represent the views of the Department of Veterans Affairs or the United States government.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hinokio, Y.; Suzuki, S.; Hirai, M.; Chiba, M.; Hirai, A.; Toyota, T. Oxidative DNA Damage in Diabetes Mellitus: Its Association with Diabetic Complications. Diabetologia 1999, 42, 995–998. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhong, Y.; Li, J.; Chen, Y.; Wang, J.J.; Ratan, R.; Zhang, S.X. Activation of Endoplasmic Reticulum Stress by Hyperglycemia Is Essential for Müller Cell–Derived Inflammatory Cytokine Production in Diabetes. Diabetes 2012, 6, 492–504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Djuric, Z.; Lu, M.H.; Lewis, S.M.; Luongo, D.A.; Chen, X.W.; Heilbrun, L.K.; Reading, B.A.; Duffy, P.H.; Hart, R.W. Oxidative DNA Damage Levels in Rats Fed Low-Fat, High-Fat, or Calorie-Restricted Diets. Toxicol. Appl. Pharm. 1992, 115, 156–160. [Google Scholar] [CrossRef]
  4. Robertson, R.P.; Harmon, J.; Tran, P.O.; Tanaka, Y.; Takahashi, H. Glucose Toxicity in β-Cells: Type 2 Diabetes, Good Radicals Gone Bad, and the Glutathione Connection. Diabetes 2003, 52, 581–587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Cnop, M.; Foufelle, F.; Velloso, L.A. Endoplasmic Reticulum Stress, Obesity and Diabetes. Trends Mol. Med. 2012, 18, 59–68. [Google Scholar] [CrossRef]
  6. Cnop, M.; Welsh, N.; Jonas, J.-C.; Jörns, A.; Lenzen, S.; Eizirik, D.L. Mechanisms of Pancreatic β-Cell Death in Type 1 and Type 2 Diabetes Many Differences, Few Similarities. Diabetes 2005, 54, S97–S107. [Google Scholar] [CrossRef] [Green Version]
  7. Eizirik, D.L.; Cardozo, A.K.; Cnop, M. The Role for Endoplasmic Reticulum Stress in Diabetes Mellitus. Endocr. Rev. 2008, 29, 42–61. [Google Scholar] [CrossRef] [Green Version]
  8. Donath, M.Y.; Dalmas, É.; Sauter, N.S.; Böni-Schnetzler, M. Inflammation in Obesity and Diabetes: Islet Dysfunction and Therapeutic Opportunity. Cell Metab. 2013, 17, 860–872. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Harding, H.P.; Ron, D. Endoplasmic Reticulum Stress and the Development of Diabetes: A Review. Diabetes 2002, 51, S455–S461. [Google Scholar] [CrossRef] [Green Version]
  10. Butcher, M.J.; Hallinger, D.; Garcia, E.; Machida, Y.; Chakrabarti, S.; Nadler, J.; Galkina, E.V.; Imai, Y. Association of Proinflammatory Cytokines and Islet Resident Leucocytes with Islet Dysfunction in Type 2 Diabetes. Diabetologia 2014, 57, 491–501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Goodarzi, M.T.; Navidi, A.A.; Rezaei, M.; Babahmadi-Rezaei, H. Oxidative Damage to DNA and Lipids: Correlation with Protein Glycation in Patients with Type 1 Diabetes. J. Clin. Lab. Anal. 2010, 24, 72–76. [Google Scholar] [CrossRef] [PubMed]
  12. Kaneto, H.; Kawamori, D.; Matsuoka, T.; Kajimoto, Y.; Yamasaki, Y. Oxidative Stress and Pancreatic β-Cell Dysfunction. Am. J. Ther. 2005, 12, 529–533. [Google Scholar] [CrossRef]
  13. Blasiak, J.; Arabski, M.; Krupa, R.; Wozniak, K.; Zadrozny, M.; Kasznicki, J.; Zurawska, M.; Drzewoski, J. DNA Damage and Repair in Type 2 Diabetes Mellitus. Mutat. Res. Fundam. Mol. Mech. Mutagen. 2004, 554, 297–304. [Google Scholar] [CrossRef] [PubMed]
  14. Horwitz, E.; Krogvold, L.; Zhitomirsky, S.; Swisa, A.; Fischman, M.; Lax, T.; Dahan, T.; Hurvitz, N.; Weinberg-Corem, N.; Klochendler, A.; et al. Beta Cell DNA Damage Response Promotes Islet Inflammation in Type 1 Diabetes. Diabetes 2018, 67, 2305–2318. [Google Scholar] [CrossRef]
  15. Oleson, B.J.; Broniowska, K.A.; Schreiber, K.H.; Tarakanova, V.L.; Corbett, J.A. Nitric Oxide Induces Ataxia Telangiectasia Mutated (ATM) Protein-Dependent ΓH2AX Protein Formation in Pancreatic β Cells. J. Biol. Chem. 2014, 289, 11454–11464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Oleson, B.J.; Broniowska, K.A.; Naatz, A.; Hogg, N.; Tarakanova, V.L.; Corbett, J.A. Nitric Oxide Suppresses β-Cell Apoptosis by Inhibiting the DNA Damage Response. Mol. Cell Biol. 2016, 36, 2067–2077. [Google Scholar] [CrossRef] [Green Version]
  17. Cheung, Y.H.; Watkinson, J.; Anastassiou, D. Conditional Meta-Analysis Stratifying on Detailed HLA Genotypes Identifies a Novel Type 1 Diabetes Locus around TCF19 in the MHC. Hum. Genet. 2011, 129, 161–176. [Google Scholar] [CrossRef] [Green Version]
  18. Harder, M.N.; Appel, E.V.R.; Grarup, N.; Gjesing, A.P.; Ahluwalia, T.S.; Jørgensen, T.; Christensen, C.; Brandslund, I.; Linneberg, A.; Sørensen, T.I.A.; et al. The Type 2 Diabetes Risk Allele of TMEM154-Rs6813195 Associates with Decreased Beta Cell Function in a Study of 6,486 Danes. PLoS ONE 2015, 10, e0120890. [Google Scholar] [CrossRef] [Green Version]
  19. Nair, A.K.; Baier, L.J. Complex Genetics of Type 2 Diabetes and Effect Size: What Have We Learned from Isolated Populations? Rev. Diabet. Stud. 2015, 12, 299–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Mahajan, A.; Go, M.J.; Zhang, W.; Below, J.E.; Gaulton, K.J.; Ferreira, T.; Horikoshi, M.; Johnson, A.D.; Ng, M.C.; Prokopenko, I.; et al. Genome-Wide Trans-Ancestry Meta-Analysis Provides Insight into the Genetic Architecture of Type 2 Diabetes Susceptibility. Nat. Genet. 2014, 46, 234–244. [Google Scholar] [CrossRef]
  21. Krautkramer, K.A.; Linnemann, A.K.; Fontaine, D.A.; Whillock, A.L.; Harris, T.W.; Schleis, G.J.; Truchan, N.A.; Marty-Santos, L.; Lavine, J.A.; Cleaver, O.; et al. Tcf19 Is a Novel Islet Factor Necessary for Proliferation and Survival in the INS-1 β-Cell Line. Am. J. Physiol. Endoc. 2013, 305, E600–E610. [Google Scholar] [CrossRef]
  22. Yang, Y.; Gurung, B.; Wu, T.; Wang, H.; Stoffers, D.A.; Hua, X. Reversal of Preexisting Hyperglycemia in Diabetic Mice by Acute Deletion of the Men1 Gene. Proc. Natl. Acad. Sci. USA 2010, 107, 20358–20363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Keller, M.P.; Paul, P.K.; Rabaglia, M.E.; Stapleton, D.S.; Schueler, K.L.; Broman, A.T.; Ye, S.I.; Leng, N.; Brandon, C.J.; Neto, E.C.; et al. The Transcription Factor Nfatc2 Regulates β-Cell Proliferation and Genes Associated with Type 2 Diabetes in Mouse and Human Islets. Plos Genet. 2016, 12, e1006466. [Google Scholar] [CrossRef]
  24. Klochendler, A.; Caspi, I.; Corem, N.; Moran, M.; Friedlich, O.; Elgavish, S.; Nevo, Y.; Helman, A.; Glaser, B.; Eden, A.; et al. The Genetic Program of Pancreatic β-Cell Replication In Vivo. Diabetes 2016, 65, 2081–2093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Barbosa-Sampaio, H.C.; Liu, B.; Drynda, R.; Ledesma, A.M.R.D.; King, A.J.; Bowe, J.E.; Malicet, C.; Iovanna, J.L.; Jones, P.M.; Persaud, S.J.; et al. Nupr1 Deletion Protects against Glucose Intolerance by Increasing Beta Cell Mass. Diabetologia 2013, 56, 2477–2486. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ndiaye, F.K.; Ortalli, A.; Canouil, M.; Huyvaert, M.; Salazar-Cardozo, C.; Lecoeur, C.; Verbanck, M.; Pawlowski, V.; Boutry, R.; Durand, E.; et al. Expression and Functional Assessment of Candidate Type 2 Diabetes Susceptibility Genes Identify Four New Genes Contributing to Human Insulin Secretion. Mol. Metab. 2017, 6, 459–470. [Google Scholar] [CrossRef] [PubMed]
  27. Mohammad, D.H.; Yaffe, M.B. 14-3-3 Proteins, FHA Domains and BRCT Domains in the DNA Damage Response. DNA Repair 2009, 8, 1009–1017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Morrison, E.A.; Musselman, C.A. Chromatin Signaling and Diseases; Academic Press Inc.: Cambridge, MA, USA, 2016; pp. 127–147. [Google Scholar]
  29. Sen, S.; Sanyal, S.; Srivastava, D.K.; Dasgupta, D.; Roy, S.; Das, C. Transcription Factor 19 Interacts with Histone 3 Lysine 4 Trimethylation and Controls Gluconeogenesis via the Nucleosome-Remodeling-Deacetylase Complex. J. Biol. Chem. 2017, 292, 20362–20378. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Mondal, P.; Sen, S.; Klein, B.J.; Tiwary, N.; Gadad, S.S.; Kutateladze, T.G.; Roy, S.; Das, C. TCF19 Promotes Cell Proliferation through Binding to the Histone H3K4me3 Mark. Biochemistry 2019, 59, 389–399. [Google Scholar] [CrossRef] [PubMed]
  31. Limsirichaikul, S.; Niimi, A.; Fawcett, H.; Lehmann, A.; Yamashita, S.; Ogi, T. A Rapid Non-Radioactive Technique for Measurement of Repair Synthesis in Primary Human Fibroblasts by Incorporation of Ethynyl Deoxyuridine (EdU). Nucleic Acids Res. 2009, 37, e31. [Google Scholar] [CrossRef]
  32. Dasika, G.K.; Lin, S.-C.J.; Zhao, S.; Sung, P.; Tomkinson, A.; Lee, E.Y.-H.P. DNA Damage-Induced Cell Cycle Checkpoints and DNA Strand Break Repair in Development and Tumorigenesis. Oncogene 1999, 18, 7883–7899. [Google Scholar] [CrossRef] [Green Version]
  33. Mi, H.; Ebert, D.; Muruganujan, A.; Mills, C.; Albou, L.-P.; Mushayamaha, T.; Thomas, P.D. PANTHER Version 16: A Revised Family Classification, Tree-Based Classification Tool, Enhancer Regions and Extensive API. Nucleic Acids Res. 2020, 49, D394–D403. [Google Scholar] [CrossRef] [PubMed]
  34. Thomas, P.D.; Campbell, M.J.; Kejariwal, A.; Mi, H.; Karlak, B.; Daverman, R.; Diemer, K.; Muruganujan, A.; Narechania, A. PANTHER: A Library of Protein Families and Subfamilies Indexed by Function. Genome Res. 2003, 13, 2129–2141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Mi, H.; Thomas, P. Protein Networks and Pathway Analysis. Methods Mol. Biol. 2009, 563, 123–140. [Google Scholar] [CrossRef]
  36. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING V11: Protein–Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res. 2018, 47, D607–D613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Morales, J.; Li, L.; Fattah, F.J.; Dong, Y.; Bey, E.A.; Patel, M.; Gao, J.; Boothman, D.A. Review of Poly (ADP-Ribose) Polymerase (PARP) Mechanisms of Action and Rationale for Targeting in Cancer and Other Diseases. Crit. Rev. Eukar Gene 2014, 24, 15–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Yang, C.-S.; Jividen, K.; Spencer, A.; Dworak, N.; Ni, L.; Oostdyk, L.T.; Chatterjee, M.; Kuśmider, B.; Reon, B.; Parlak, M.; et al. Ubiquitin Modification by the E3 Ligase/ADP-Ribosyltransferase Dtx3L/Parp9. Mol. Cell 2017, 66, 503–516 e5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Hornung, V.; Hartmann, R.; Ablasser, A.; Hopfner, K.-P. OAS Proteins and CGAS: Unifying Concepts in Sensing and Responding to Cytosolic Nucleic Acids. Nat. Rev. Immunol. 2014, 14, 521–528. [Google Scholar] [CrossRef] [PubMed]
  40. Kondratova, A.A.; Cheon, H.; Dong, B.; Holvey-Bates, E.G.; Hasipek, M.; Taran, I.; Gaughan, C.; Jha, B.K.; Silverman, R.H.; Stark, G.R. Suppressing PARylation by 2′,5′-oligoadenylate Synthetase 1 Inhibits DNA Damage-induced Cell Death. Embo. J. 2020, 39, e101573. [Google Scholar] [CrossRef]
  41. Verhelst, J.; Parthoens, E.; Schepens, B.; Fiers, W.; Saelens, X. Interferon-Inducible Protein Mx1 Inhibits Influenza Virus by Interfering with Functional Viral Ribonucleoprotein Complex Assembly. J. Virol. 2012, 86, 13445–13455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Miyashita, M.; Oshiumi, H.; Matsumoto, M.; Seya, T. DDX60, a DEXD/H Box Helicase, Is a Novel Antiviral Factor Promoting RIG-I-like Receptor-Mediated Signaling. Mol. Cell Biol. 2011, 31, 3802–3819. [Google Scholar] [CrossRef] [Green Version]
  43. Honke, N.; Shaabani, N.; Zhang, D.-E.; Hardt, C.; Lang, K.S. Multiple Functions of USP18. Cell Death Dis. 2016, 7, e2444. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Roopra, A. MAGIC: A Tool for Predicting Transcription Factors and Cofactors Driving Gene Sets Using ENCODE Data. PLoS Comput. Biol. 2020, 16, e1007800. [Google Scholar] [CrossRef]
  45. Sadzak, I.; Schiff, M.; Gattermeier, I.; Glinitzer, R.; Sauer, I.; Saalmüller, A.; Yang, E.; Schaljo, B.; Kovarik, P. Recruitment of Stat1 to Chromatin Is Required for Interferon-Induced Serine Phosphorylation of Stat1 Transactivation Domain. Proc. Natl. Acad. Sci. USA 2008, 105, 8944–8949. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Varinou, L.; Ramsauer, K.; Karaghiosoff, M.; Kolbe, T.; Pfeffer, K.; Müller, M.; Decker, T. Phosphorylation of the Stat1 Transactivation Domain Is Required for Full-Fledged IFN-γ-Dependent Innate Immunity. Immunity 2003, 19, 793–802. [Google Scholar] [CrossRef] [Green Version]
  47. Cai, X.; Chiu, Y.-H.; Chen, Z.J. The CGAS-CGAMP-STING Pathway of Cytosolic DNA Sensing and Signaling. Mol. Cell 2014, 54, 289–296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Willcox, A.; Richardson, S.J.; Bone, A.J.; Foulis, A.K.; Morgan, N.G. Analysis of Islet Inflammation in Human Type 1 Diabetes. Clin. Exp. Immunol. 2009, 155, 173–181. [Google Scholar] [CrossRef]
  49. Åkerblom, H.K.; Vaarala, O.; Hyöty, H.; Ilonen, J.; Knip, M. Environmental Factors in the Etiology of Type 1 Diabetes. Am. J. Med. Genet. 2002, 115, 18–29. [Google Scholar] [CrossRef] [PubMed]
  50. Solinas, G.; Becattini, B. JNK at the Crossroad of Obesity, Insulin Resistance, and Cell Stress Response. Mol. Metab. 2017, 6, 174–184. [Google Scholar] [CrossRef]
  51. Nishimura, S.; Manabe, I.; Nagasaki, M.; Eto, K.; Yamashita, H.; Ohsugi, M.; Otsu, M.; Hara, K.; Ueki, K.; Sugiura, S.; et al. CD8+ Effector T Cells Contribute to Macrophage Recruitment and Adipose Tissue Inflammation in Obesity. Nat. Med. 2009, 15, 914–920. [Google Scholar] [CrossRef] [PubMed]
  52. Tatsch, E.; Carvalho, J.A.M.D.; Hausen, B.S.; Bollick, Y.S.; Torbitz, V.D.; Duarte, T.; Scolari, R.; Duarte, M.M.M.F.; Londero, S.W.K.; Vaucher, R.A.; et al. Oxidative DNA Damage Is Associated with Inflammatory Response, Insulin Resistance and Microvascular Complications in Type 2 Diabetes. Mutat. Res. Fundam. Mol. Mech. Mutagen. 2015, 782, 17–22. [Google Scholar] [CrossRef]
  53. Nakad, R.; Schumacher, B. DNA Damage Response and Immune Defense: Links and Mechanisms. Front. Genet. 2016, 7, 147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Tavana, O.; Puebla-Osorio, N.; Sang, M.; Zhu, C. Absence of P53-Dependent Apoptosis Combined With Nonhomologous End-Joining Deficiency Leads to a Severe Diabetic Phenotype in Mice. Diabetes 2010, 59, 135–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Brzostek-Racine, S.; Gordon, C.; Scoy, S.V.; Reich, N.C. The DNA Damage Response Induces IFN. J. Immunol. 2011, 187, 5336–5345. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Unterholzner, L.; Dunphy, G. CGAS-Independent STING Activation in Response to DNA Damage. Mol. Cell Oncol. 2019, 6, 1–3. [Google Scholar] [CrossRef] [Green Version]
  57. Bai, J.; Liu, F. The CGAS-CGAMP-STING Pathway: A Molecular Link between Immunity and Metabolism. Diabetes 2019, 68, 1099–1108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Wållberg, M.; Cooke, A. Immune Mechanisms in Type 1 Diabetes. Trends Immunol. 2013, 34, 583–591. [Google Scholar] [CrossRef]
  59. Mboko, W.P.; Mounce, B.C.; Wood, B.M.; Kulinski, J.M.; Corbett, J.A.; Tarakanova, V.L. Coordinate Regulation of DNA Damage and Type I Interferon Responses Imposes an Antiviral State That Attenuates Mouse Gammaherpesvirus Type 68 Replication in Primary Macrophages. J. Virol. 2012, 86, 6899–6912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Kay, J.; Thadhani, E.; Samson, L.; Engelward, B. Inflammation-Induced DNA Damage, Mutations and Cancer. DNA Repair 2019, 83, 102673. [Google Scholar] [CrossRef]
  61. Michalska, A.; Blaszczyk, K.; Wesoly, J.; Bluyssen, H.A.R. A Positive Feedback Amplifier Circuit That Regulates Interferon (IFN)-Stimulated Gene Expression and Controls Type I and Type II IFN Responses. Front. Immunol. 2018, 9, 1135. [Google Scholar] [CrossRef] [Green Version]
  62. Frontini, M.; Vijayakumar, M.; Garvin, A.; Clarke, N. A ChIP–Chip Approach Reveals a Novel Role for Transcription Factor IRF1 in the DNA Damage Response. Nucleic Acids Res. 2009, 37, 1073–1085. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Zhong, M.; Henriksen, M.A.; Takeuchi, K.; Schaefer, O.; Liu, B.; Hoeve, J.t.; Ren, Z.; Mao, X.; Chen, X.; Shuai, K.; et al. Implications of an Antiparallel Dimeric Structure of Nonphosphorylated STAT1 for the Activation–Inactivation Cycle. Proc. Natl. Acad. Sci. USA. 2005, 102, 3966–3971. [Google Scholar] [CrossRef] [Green Version]
  64. Durocher, D.; Jackson, S.P. The FHA Domain. Febs. Lett. 2002, 513, 58–66. [Google Scholar] [CrossRef] [Green Version]
  65. Li, J.; Lee, G.I.; Doren, S.R.V.; Walker, J.C. The FHA Domain Mediates Phosphoprotein Interactions. J. Cell Sci. 2000, 23, 4143–4149. [Google Scholar] [CrossRef]
  66. Matsuoka, S.; Ballif, B.A.; Smogorzewska, A.; McDonald, E.R.; Hurov, K.E.; Luo, J.; Bakalarski, C.E.; Zhao, Z.; Solimini, N.; Lerenthal, Y.; et al. ATM and ATR Substrate Analysis Reveals Extensive Protein Networks Responsive to DNA Damage. Science 2007, 316, 1160–1166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Li, H. Aligning Sequence Reads, Clone Sequences and Assembly Contigs with BWA-MEM. arXiv 2013, arXiv:1303.3997v2. [Google Scholar]
  68. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  69. Li, B.; Dewey, C.N. RSEM: Accurate Transcript Quantification from RNA-Seq Data with or without a Reference Genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 2009, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) Overexpression of human transcription factor 19 (TCF19) in INS-1 cells was confirmed by qRT-PCR. (B) Western blot against anti-myc tag confirms overexpression of human TCF19 (hsTCF19) (C) Overexpression of hsTCF19 in INS-1 cells leads to increased 3H-thymidine incorporation (n = 5). (D) hsTCF19 overexpression does not lead to any significant changes in cell cycle gene expression (n = 5) (E) Overexpression of hsTCF19 in INS-1 cells does not affect cell viability (n = 5). Data are means ± SEM * p < 0.05.
Figure 1. (A) Overexpression of human transcription factor 19 (TCF19) in INS-1 cells was confirmed by qRT-PCR. (B) Western blot against anti-myc tag confirms overexpression of human TCF19 (hsTCF19) (C) Overexpression of hsTCF19 in INS-1 cells leads to increased 3H-thymidine incorporation (n = 5). (D) hsTCF19 overexpression does not lead to any significant changes in cell cycle gene expression (n = 5) (E) Overexpression of hsTCF19 in INS-1 cells does not affect cell viability (n = 5). Data are means ± SEM * p < 0.05.
Metabolites 11 00513 g001
Figure 2. RNA-seq analysis of INS-1 cells overexpressing hsTCF19 identifies upregulated genes that form a tight node of interconnected genes with roles in viral, immune, and interferon response pathways (A) PANTHER overrepresentation test for GO biological processes on >2- fold upregulated genes show enrichment for genes involved in viral responses (red), regulation of immune response processes (blue), regulation of cytokines (purple), and antigen presentation (green). Additionally, other processes such as regulation of double strand break repair, apoptotic signaling pathway in response to ER stress, and positive regulation of NF-κB signaling were also overrepresented (black) (FDR < 5%). (B) STRING analysis on differentially expressed genes with >2-fold upregulation shows tightly interconnected network of genes.
Figure 2. RNA-seq analysis of INS-1 cells overexpressing hsTCF19 identifies upregulated genes that form a tight node of interconnected genes with roles in viral, immune, and interferon response pathways (A) PANTHER overrepresentation test for GO biological processes on >2- fold upregulated genes show enrichment for genes involved in viral responses (red), regulation of immune response processes (blue), regulation of cytokines (purple), and antigen presentation (green). Additionally, other processes such as regulation of double strand break repair, apoptotic signaling pathway in response to ER stress, and positive regulation of NF-κB signaling were also overrepresented (black) (FDR < 5%). (B) STRING analysis on differentially expressed genes with >2-fold upregulation shows tightly interconnected network of genes.
Metabolites 11 00513 g002aMetabolites 11 00513 g002b
Figure 3. Human TCF19 overexpression in human islets upregulates key DNA damage repair genes as well as viral response genes. (A) hsTCF19 overexpression in human islets was confirmed with qRT-PCR. Control represents human islets transfected with an empty vector and hsTCF19 represents human islets transfected with TCF19 plasmid. (B) hsTCF19 overexpression in human islets leads to upregulation of genes (qRT-PCR) that were also upregulated in INS-1 cells (n = 5). Gene names are shown on the horizontal axis. Data are means ± SEM * p < 0.05.
Figure 3. Human TCF19 overexpression in human islets upregulates key DNA damage repair genes as well as viral response genes. (A) hsTCF19 overexpression in human islets was confirmed with qRT-PCR. Control represents human islets transfected with an empty vector and hsTCF19 represents human islets transfected with TCF19 plasmid. (B) hsTCF19 overexpression in human islets leads to upregulation of genes (qRT-PCR) that were also upregulated in INS-1 cells (n = 5). Gene names are shown on the horizontal axis. Data are means ± SEM * p < 0.05.
Metabolites 11 00513 g003
Figure 4. MAGIC analysis on the list of upregulated genes after TCF19 expression in INS-1 cells identifies significant enrichment for genes with known ChIP signals for STAT1, STAT2, and IRF1 in their promoters (FDR < 10%). A higher score reflects increased likelihood that the factor is enriched for binding gene regions in the list of upregulated genes. Genes included in analysis were those that were upregulated more than 2-fold with an associated FDR < 5%.
Figure 4. MAGIC analysis on the list of upregulated genes after TCF19 expression in INS-1 cells identifies significant enrichment for genes with known ChIP signals for STAT1, STAT2, and IRF1 in their promoters (FDR < 10%). A higher score reflects increased likelihood that the factor is enriched for binding gene regions in the list of upregulated genes. Genes included in analysis were those that were upregulated more than 2-fold with an associated FDR < 5%.
Metabolites 11 00513 g004
Figure 5. Overexpression of TCF19 in INS-1 cells does not lead to increased activation of transcription factor targets but leads to decreased activation of STAT1 Y701 phosphorylation. Representative images of two western blot replicates along with analysis of all replicates are shown. Densitometry with Image J 1.44o was used to quantify the bands on the western blots, which were then normalized to the housekeeper protein band. (A) Serine 727 phospho-STAT1/STAT1 protein expression does not show a statistically significant difference between control and hsTCF19 overexpressing cells (n = 5) (phospho-STAT1~91 kDa, STAT1~84, 91 kDa). (B) Tyrosine 701 phospho-STAT1/STAT1 protein expression shows statistically significant decrease in hsTCF19 overexpressing cells compared to control (n = 3) (phospho-STAT1~84, 91 kDa). (C) IRF1 protein levels are not significantly different. Representative western blots in the figure have IRF1 levels normalized to GAPDH. Three other replicates are normalized to beta tubulin (n = 5) (IRF1~48 kDa). (D) There is no difference in phospho NF-κB/NF-κB levels with hsTCF19 overexpression (n = 3) (phospho NF-κB~65 kDa, NF-κB~65 kDa). All data are means ± SEM * p < 0.05.
Figure 5. Overexpression of TCF19 in INS-1 cells does not lead to increased activation of transcription factor targets but leads to decreased activation of STAT1 Y701 phosphorylation. Representative images of two western blot replicates along with analysis of all replicates are shown. Densitometry with Image J 1.44o was used to quantify the bands on the western blots, which were then normalized to the housekeeper protein band. (A) Serine 727 phospho-STAT1/STAT1 protein expression does not show a statistically significant difference between control and hsTCF19 overexpressing cells (n = 5) (phospho-STAT1~91 kDa, STAT1~84, 91 kDa). (B) Tyrosine 701 phospho-STAT1/STAT1 protein expression shows statistically significant decrease in hsTCF19 overexpressing cells compared to control (n = 3) (phospho-STAT1~84, 91 kDa). (C) IRF1 protein levels are not significantly different. Representative western blots in the figure have IRF1 levels normalized to GAPDH. Three other replicates are normalized to beta tubulin (n = 5) (IRF1~48 kDa). (D) There is no difference in phospho NF-κB/NF-κB levels with hsTCF19 overexpression (n = 3) (phospho NF-κB~65 kDa, NF-κB~65 kDa). All data are means ± SEM * p < 0.05.
Metabolites 11 00513 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, G.H.; Fontaine, D.A.; Lodh, S.; Blumer, J.T.; Roopra, A.; Davis, D.B. TCF19 Impacts a Network of Inflammatory and DNA Damage Response Genes in the Pancreatic β-Cell. Metabolites 2021, 11, 513. https://doi.org/10.3390/metabo11080513

AMA Style

Yang GH, Fontaine DA, Lodh S, Blumer JT, Roopra A, Davis DB. TCF19 Impacts a Network of Inflammatory and DNA Damage Response Genes in the Pancreatic β-Cell. Metabolites. 2021; 11(8):513. https://doi.org/10.3390/metabo11080513

Chicago/Turabian Style

Yang, Grace H., Danielle A. Fontaine, Sukanya Lodh, Joseph T. Blumer, Avtar Roopra, and Dawn Belt Davis. 2021. "TCF19 Impacts a Network of Inflammatory and DNA Damage Response Genes in the Pancreatic β-Cell" Metabolites 11, no. 8: 513. https://doi.org/10.3390/metabo11080513

APA Style

Yang, G. H., Fontaine, D. A., Lodh, S., Blumer, J. T., Roopra, A., & Davis, D. B. (2021). TCF19 Impacts a Network of Inflammatory and DNA Damage Response Genes in the Pancreatic β-Cell. Metabolites, 11(8), 513. https://doi.org/10.3390/metabo11080513

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop