- Research article
- Open Access
Areca catechu-(Betel-nut)-induced whole transcriptome changes in a human monocyte cell line that may have relevance to diabetes and obesity; a pilot study
BMC Endocrine Disorders volume 21, Article number: 165 (2021)
Betel-nut consumption is the fourth most common addictive habit globally and there is good evidence linking the habit to obesity, type 2 diabetes (T2D) and the metabolic syndrome. The aim of our pilot study was to identify gene expression relevant to obesity, T2D and the metabolic syndrome using a genome-wide transcriptomic approach in a human monocyte cell line incubated with arecoline and its nitrosated products.
The THP1 monocyte cell line was incubated separately with arecoline and 3-methylnitrosaminopropionaldehyde (MNPA) in triplicate for 24 h and pooled cDNA indexed paired-end libraries were sequenced (Illumina NextSeq 500). After incubation with arecoline and MNPA, 15 and 39 genes respectively had significant changes in their expression (q < 0.05, log fold change 1.5). Eighteen of those genes have reported associations with T2D and obesity in humans; of these genes there was most marked evidence for CLEC10A, MAPK8IP1, NEGR1, NQ01 and INHBE genes.
Our preliminary studies have identified a large number of genes relevant to obesity, T2D and metabolic syndrome whose expression was changed significantly in human TPH1 cells following incubation with betel-nut derived arecoline or with MNPA. These findings require validation by further cell-based work and investigation amongst betel-chewing communities.
Obesity and T2D are reaching epidemic proportions worldwide, but particularly so in South Asian communities living in the Indian-subcontinent or who have migrated to other countries . In the UK there is a three to four fold increase in T2D prevalence in South Asians compared to the general population; furthermore, the disease manifests at a 10–15 years younger age and is strongly associated with the metabolic syndrome and cardiovascular disease. Apart from lifestyle, potentially reversible environmental factors driving this disease are largely unknown.
Betel quid consumption is the fourth most common addictive habit, used by 10 % of the global population and very common in South Asians. The link between cancer risks (oropharyngeal, oesophageal and hepatocellular) and the ‘betel-chewing’ habit is well established [2,3,4]. Evidence also demonstrates associations between betel consumption and obesity, the metabolic syndrome and T2D. In a meta-analysis of 17 Asian studies, betel quid chewing was found to be significantly associated with obesity (relative risk (RR) = 1.47), metabolic syndrome (RR = 1.51), diabetes (RR = 1.47), hypertension (RR = 1.45) and cardiovascular disease (cardiovascular disease: RR = 1.2) . Furthermore, in the Keelung Community Integrated Screening program from Taiwan, paternal betel-chewing was associated, dose-wise, with increases in early onset metabolic syndrome in their never-chewing offspring, while betel-chewing in adults increased their risks, dose-wise, of early onset T2D, cardiovascular disease and metabolic syndrome [6,7,8]. These data in humans support earlier data reported in CD1 mice, where it was found that a proportion of betel-fed adult mice developed hyperglycaemia and obesity and, most remarkably, that amongst their never-betel-fed offspring hyperglycaemia was detected in up to the 4th generation and that the vertical transmission of hyperglycaemia was associated with paternal, but not maternal, hyperglycaemia .
Betel quids (also known as paan) usually contain betel (areca catechu) nut, slaked lime and sometimes tobacco wrapped in piper betle leaf  though in Taiwan tobacco is not used. Nitrosation of the major arecal alkaloid, arecoline, generates MNPA and 3-methylnitrosaminopropionitrile (MNPN); both these compounds have been accepted as carcinogens . Many nitroso-compounds have been reported as being diabetogenic. For instance, all reported survivors of poisoning by the nitroso-rat poison Vacor (Pyrinuron (Pyriminil)) in the US have developed ketosis-prone insulin-dependent diabetes while ~ 20 % of survivors of Vacor poisoning in Korea developed ketosis-prone type 1 diabetes (T1D) but the other ~ 80 % developed non-insulin dependent diabetes; the difference being attributed to the lower content of the active ingredient in the Vacor packs sold in Korea than in those sold in America . Furthermore, studies in mouse strains normally free of diabetes have regularly shown induction of ‘non-immune T1D’ following single large doses of nitroso compounds, classically streptozotocin (STZ) . In contrast, repeated small doses of STZ induce a non-insulin dependent form of diabetes with enlarged pancreatic islets and beta cell abnormalities whose appearance matches those seen in pancreatic islets in humans developing T2D. With direct relevance to the current study, these same abnormalities were also seen in the pancreatic islets of the 8.3 % of young adult CD1 mice that developed hyperglycaemia with obesity following being fed standard RM1 rodent chow containing ground Areca catechu nuts for 2–6 days; interestingly, those changes also appeared in the islets of the 21.7 % of the never-betel-fed F1 offspring of the originally betel-fed mice that themselves developed obesity and hyperglycaemia .
These studies provides evidence that betel-chewing might be one of the aetiological factors for the increases in T2D and associated metabolic disorders in South Asians [12, 14]. Mechanisms that might link betel chewing with these disorders include inflammation, increases in hepatic synthesis of lipids and glucose, in adipogenesis, in adipose tissue glucose uptake, reductions in lipolysis and glycolysis, or neurological, hepatic or intestinal effects on appetite and adverse effects on vitamin D metabolism [15, 16].
In the present ‘proof of principle’ study we sought to investigate possible gene effects that might explain the links between T2D, obesity and related disorders and exposure to arecoline and its nitrosated products in a human monocyte derived cell line (THP1) using a whole transcriptome sequencing approach. THP1 was chosen due the central role that low-grade inflammation plays in the underlying causes and progression of obesity (for example through adipose tissue macrophages), metabolic syndrome and related disorders including both cardiovascular disease and T2D [17, 18]. It was also hoped to determine whether the strength of the evidence for any of the genes significantly affected might warrant further cell-based studies and investigation amongst betel-chewing communities.
Determination of concentration of arecoline, MNPA and MNPN to be used in cell based assays
Cell viability using serial dilutions of arecoline extracts in culture media was assessed by phase contrast light microscopy. Representative images are seen in Additional Fig. 1. Control THP-1 cells are clear, round and uniform in shape, size and clarity. When THP1 cells are incubated for 48 h or more with arecoline concentrations of either 10 or 50 μg/ml the cells look healthy with similar regular appearance by phase contrast microscopy to controls. In contrast, cells incubated with 200 or 1000 μg/ml, appear unhealthy being shrunken and of irregular shape with loss of clarity and release of much disorganised cellular debris. No cell death or damage was observed in incubates containing MNPA or MNPN up to a concentration of 200 μg/ml, though at very high concentrations, media containing more than 1 % of the solvent (methanol and ethyl acetate) was toxic to the cells as assessed by phase light microscopy. A MTT assay for cell viability and proliferation was also performed (Additional Fig. 2), that demonstrates the THP-1 cells are healthy after incubation with a concentration of 50 μg/ml of arecoline but that toxicity was induced at 200 μg/ml and above, supporting the phase contrast light microscopy data.
Selection of samples for whole transcriptome analysis
Whilst an inflammatory cellular response was seen after incubation of the THP1 cells with either arecoline, MNPA or Phorbol 12-myristate 13-acetate (PMA), minimal responses were found at 6, 24 or 48 h with MNPN and therefore, we did not proceed further with investigation of the last of these compounds (Additional Fig. 3).
PMA induces a macrophage-like phenotype differentiation in THP-1 cells, functionally mimicking many aspects of primary human macrophages. PMA, as a specific activator of Protein Kinase C, activates nuclear factor-kappa B and is regularly used as an internal control to test for cytokine and biological responsiveness of THP1 cells. In the absence of any response, one would have to conclude that data from the planned experiments on that particular cell line could not be used. We therefore decided to proceed with the incubations at 24 h. All transcriptome expression experiments (including whole transcriptome) were run in triplicate from independent rounds of cell culture analysis.
Whole transcriptome analysis
Incubation with arecoline
Two hundred seventy-five gene hits were identified with a q < 0.05 (Additional Table 1) reducing to 15 with a log-fold change in either direction of 1.5 (Table 1). The relevance of these genes to diabetes, obesity and/or metabolic syndrome was assessed independently by 2 of us (BJB and GAH) with reference to GeneCards, The Type 2 Diabetes Knowledge Portal and PubMed as described in the methods. Four genes were identified of potential interest:
Insulin Like Growth Factor Binding Protein 3 (IGFB3 non-logged fold change 0.08),
C-Type Lectin Domain Containing 10 A (CLEC10A fold change 0.14),
Junction Plakoglobin (JUP fold change 0.21).
Mitogen-Activated Protein Kinase 8 Interacting Protein 1 (MAPK8IP1 fold change 0.21.
Pathway (Metascape) analysis of the 275 genes (Fig. 1) listed in the Additional Table 1, revealed 5 significant pathways after statistical correction by the false discovery rate (FDR): namely, myeloid cell activation involved in immune response, cellular response to thyroid hormone stimulus, responses to toxic substances and hematopoietic cell lineage.
Incubation with MNPA
359 gene hits were identified after incubation with MNPA with a q < 0.05 (Additional Table 2) reducing to 39 with a log-fold change of +/-1.5 (Table 2). The relevance of these genes to diabetes, obesity and/or metabolic syndrome was assessed independently by 2 of us (BJB and GAH) with reference to GeneCards, The Type 2 Diabetes Knowledge Portal and PubMed as described in the methods. Fourteen genes were identified of potential interest:
Gliomedin (GLDN fold change 0.11)
Glutamate Receptor Interacting Protein 1 (GRIP1 fold change 0.15).
Neuronal Growth Regulator 1 (NEGR1 fold change 0.14).
Potassium Voltage-Gated Channel Subfamily Q Member 5 (KCNQ5 fold change 0.21).
Cytotoxic And Regulatory T-Cell Molecule (CRTAM fold change 4.9).
NAD(P)H Quinone Dehydrogenase 1 (NQO1 fold change 4.9).
Semaphorin 6B (SEMA6B fold change 4.9).
Inhibin Subunit Beta E (INHBE fold change 5.4).
Clusterin (CLU fold change 6.7).
Spectrin Alpha, Erythrocytic 1 (SPTA1 fold change 9.3).
Heme Oxygenase 1 (HMOX1 fold change 8.4).
Transmembrane Protein 140 (TMEM140 fold change 10.0).
Sarcoglycan Gamma (SGCG fold change 24.8).
Triggering Receptor Expressed On Myeloid Cells Like 4 (TREML4 fold change 26.7).
Pathway (Metascape) analysis of the 359 genes (Fig. 2) listed in the Additional Table 2, revealed 5 significant pathways after statistical correction by the FDR: namely, regulation of cell adhesion, response to inorganic substances, apoptotic signalling pathway, response to toxic substances and regulation of the innate immune response.
Gene expression changes in common with arecoline and MNPA
Using the strict criteria of q < 0.05 and log fold change of 1.5 there is only one overlapping gene, H3F3AP4, which is the H3 Histone Pseudogene 6 (chromosome location 2q31.1) of unknown function, although it might regulate gene function as a short interfering RNA or microRNA decoy with probable relevance to betel-induced carcinogenesis. Using the less stringent criteria of ‘q < 0.05 only’ then there are 69 of 275 (25 %) genes identified with arecoline that are also affected by its nitrosated product MNPA (Additional Fig. 4; data available from authors on request).
The aim of our study was to identify genes known to be associated with diabetes, obesity and metabolic syndrome whose expression was significantly altered by exposure to the arecal compounds arecoline and its nitrosated metabolite, MNPA. Whole transcriptome analysis by RNA-Seq of the human monocyte line THP1 did reveal a significant number of genes that are either downregulated or upregulated in response to incubation with arecoline or with MNPA. It was also hoped to determine whether the strength of the evidence for any of the genes found to be significantly affected might warrant further cell-based work and investigation amongst betel-chewing communities.
Consistent with the established effects of betel-quid chewing in humans, a number of cellular pathways and genes have been identified as being significantly affected by the arecal compounds used in our approach; these genes are known to relate to immune responses, to cell differentiation and lineage, to cellular responses to toxic and inorganic substances and to the development of obesity and T2D in humans. Other genes known to be associated with carcinogenesis were also found to have altered expression, as expected; some genes affected are associated with immune function and a smaller number of affected genes are related to neural development and could prove to be associated with susceptibility to addiction.
Genes of specific relevance to obesity, diabetes and metabolic syndrome, based on published studies and genetic evidence, that showed good evidence suggesting that arecoline and/or MNPA may alter their expression, included (Table 3): -.
C-Type Lectin Domain Containing 10 A (CLEC10A), a calcium dependent endocytic receptor also known as the macrophage galactose-type lectin (MGL or CD301). It has been demonstrated to have a role in regulating adaptive and innate immune responses and is expressed in adipose tissue macrophages where it is associated with phenotypic switching of ATM subclasses in mice that then demonstrate either a lean or an obese phenotype . Furthermore, evidence in humans demonstrates that missense and protein truncating mutations of the CLEC10A gene are strongly associated with the development of T2D (http://www.type2diabetesgenetics.org/gene/geneInfo/CLEC10A). In earlier rodent experiments MGL1 was found to be a novel regulator of monocyte trafficking in adipose tissue in response to dietary induced obesity [19, 20].
Mitogen-Activated Protein Kinase 8 Interacting Protein 1 (MAPK8IP1) gene encodes a regulator of pancreatic beta-cell function; it is also expressed in a large number of tissues including many associated with immune function and is a known trans-activator of the glucose transporter GLUT2. MAPK8IP1 has a strong association with T2D with a missense mutation found in one family and, in vitro, that specific mutation has been found to be a key down-regulator of beta cell function .
Neuronal Growth Regulator 1 (NEGR1) is involved in cell adhesion. Multiple genome-wide association studies demonstrate strong genome-wide association (GWAS) signals for this gene with BMI, waist circumference and T2D  and certain mutations of this gene lead to Niemannn-Pick disease, a rare inherited metabolic disorder. (http://www.type2diabetesgenetics.org/gene/geneInfo/NEGR1). Furthermore, NEGR1 knockout mice develop increased adiposity, including increased hepatocyte fat deposition, together with increases in glycaemia and in fasting serum insulin levels .
NAD(P)H Quinone Dehydrogenase 1 (NQO1) gene is a member of the NAD(P)H dehydrogenase (quinone) family, encodes a cytoplasmic 2-electron reductase (Entrez Gene) and is part of the antioxidant defence system. There is strong genetic evidence to support an association between NQO1 variants by GWAS and increased risks of T2D and of an increased BMI (http://www.type2diabetesgenetics.org/gene/geneInfo/NQO1). NQO1 is highly expressed in human adipose tissue and its expression is reduced during diet-induced weight loss; furthermore its expression correlates directly with adiposity, glycaemia and markers of liver dysfunction . Together, these findings indicate a role for NQO1 in the aetiology of obesity and T2D.
Inhibin Subunit Beta E (INHBE) gene is a member of the Transforming Growth Factor (TGF) beta superfamily. The transcribed peptide Activin E is ubiquitously expressed in a large number of normal tissues, many being known to be especially active during cell proliferation, apoptosis, immune response and hormone secretion. The highest expression is found in the liver where it acts as a hepatokine with effects on energy homeostasis in both brown and white adipose tissue . The candidacy of the INHBE gene is further supported by strong GWAS signals associating it with cardiometabolic traits including raised serum triglycerides and overt coronary heart disease (http://www.type2diabetesgenetics.org/gene/geneInfo/INHBE).
Other genes with evidence for relevance to T2D, Metabolic syndrome and obesity include the following (also included in Table 3):
Glutamate receptor interacting protein 1 (GRIP1) mediates the trafficking and membrane organisation of a number of trans-membrane proteins in various cells including neurons and macrophages. Obese mice with a conditional knockout of GRIP1 in macrophages develop massive macrophage infiltration and inflammation in many metabolically active tissues leading to many of the features that are associated with the metabolic syndrome such as hepatic steatosis, hyperglycaemia and increased insulin resistance . Clusterin (CLU) is a molecular chaperone. Secretory clusterin is also known as ApoJ which has recently been identified as a novel hepatokine and deletion of hepatic ApoJ leads to insulin resistance and glucose tolerance . Furthermore, in humans, serum ApoJ levels correlate directly with increases in insulin resistance but are decreased by rosiglitazone treatment . Insulin Like Growth Factor Binding Protein 3 (IGFBP3) is the most abundant of 6 IGF-binding proteins. Important interactions have been observed between IGFBP3, vitamin D metabolism and obesity . Furthermore, in people with T2D IGFB3 levels may inversely contribute to accelerated cerebrovascular disease . Potassium Voltage-Gated Channel Subfamily Q Member 5 (KCNQ5) is a component of potassium channels and a strong GWAS association is seen between KCNQ5 and body mass index (http://www.type2diabetesgenetics.org/gene/geneInfo/KCNQ5). Cytotoxic And Regulatory T-Cell Molecule (CRTAM) is a Protein Coding gene that has a role in the innate immune system and has also been proposed as a potential determinant of insulin secretion. A strong GWAS association is seen between CRTAM and both body mass index and systolic blood pressure (http://www.type2diabetesgenetics.org/gene/geneInfo/CRTAM). Spectrin (SPTA1) is a component of the erythrocyte plasma membrane. A strong association is seen between SPTA1 and separately with HbA1c (GWAS) and with T2D after adjustment for BMI (mainly missense mutations) (http://www.type2diabetesgenetics.org/gene/geneInfo/SPTA1).
A weakness of our approach is that the THP1 cell line is derived from a child with monocytic leukemia, however, in its favour is that these cells have been widely used to study monocyte and macrophage function . Furthermore, their stable genetic background (compared to primary cells) with reproducible treatment responses makes this transformed cell line ideal for use in preliminary pilot studies aiming to justify further cellular and in vivo studies . However, after incubation of THP1 cells with arecoline or MNPA we cannot, therefore, be sure if we are dealing with monocytes or macrophages as is described in the literature after PMA stimulation. Various compounds have been isolated and identified from Areca catechu nuts including alkaloids, tannins, flavones, triterpenes, steroids, and fatty acids . We chose to focus on arecoline and two of the many nitrosated products of the arecal alkaloids (MNPA and MNPN), both because they are known to be the most carcinogenic of them and because low-dose nitrosamines cause T2D experimentally and in humans [12, 14]. Tissue and cell specific levels of arecoline derived from betel nut and its nitrosated products including MNPA are difficult to gauge. The average arecoline concentration in plasma in those chewing betel nut is 7.0 ng/ml (+/- SD of 10.7); in contrast the concentrations of MNPA are not reported . The range of arecoline concentrations affecting different tissues adversely is clearly very wide as reviewed by Liu and colleagues  with a minimal active concentration/dose influencing macrophage expression in vitro of 10− 5 mol/l. Similarly, there is a wide range of concentrations (20–111 μg/ml) leading to minimal toxic effects on various cell types; 50 μg/ml had a minimal toxic effect on mice bone marrow and human peripheral blood lymphocytes generating reactive oxygen species . Hence the administered concentrations in our experiments are consistent with the literature, but nonetheless, we cannot be certain that the concentrations we used are definitely physiologically relevant, which is a limitation of our approach.
Unfortunately, for technical reasons, we did not get results with MNPN. Areca catechu chewing quids often contain various other additives such as slaked lime, spices, sweeteners, and are wrapped in leaves of the Piper betle vine; furthermore, in many countries other than Taiwan they often contain chewing tobacco . We cannot therefore exclude the possibility that many major effects of chewing betel quids in humans may be due to ingestion of betel quid components other than those from the Areca catechu nut. However, obesity and hyperglycaemia were induced in CD1 mice fed Areca catechu nut without any other betel quid component  and this data contributed to our focus on the findings for genes associated with those particular disorders in humans.
In the biostatistical analysis we have been conservative in applying a q value of < 0.05 and a log fold change of x1.5 in order to maximise the chances of identifying relevant transcripts from our exploratory genome-wide experiments. Further experiments will be required to assess the effect of either upregulation or downregulation of these genes in suitable models for detecting changes in function of the particular genes identified in this way.
In summary, the gene expression profiles described, induced by arecoline and MNPA derived from Areca catechu-(Betel-nut) would act through a number of mechanisms including regulation of adaptive and immune responses, transactivation of GLUT2, anti-oxidant defences and secretion of hepatocytokines. We would, therefore, hypothesise that the cumulative effect of betel-nut ingestion in humans, through the described changes in gene expression, would contributes to a predisposition to obesity, T2D and thereby, hepatic steatosis and thereby to an increased risk of cardiovascular disease. This would be supported by the population and rodent studies as noted in the introduction.
This pilot study has identified a large number of genes whose expression was changed significantly in human TPH1 cells following incubation with arecoline and MNPA and that are known to be associated with increased risks of obesity and T2D in humans. These findings suggest that further investigation of these genes is warranted both in vitro and in humans from betel-chewing communities, and that our findings should be placed in the public domain so that they are available to the scientific community.
The aim of our pilot study was to identify gene expression relevant to obesity, T2D and the metabolic syndrome using a genome-wide transcriptomic approach in a human monocyte cell line incubated with arecoline and its nitrosated products.
The THP1 (human acute monocytic leukemia derived; ATCC number TIB-202 purchased from Thermofisher) cell line  was regularly maintained in RPMI 1640 medium containing GlutaMAX™, supplemented with 10 % FCS (Gibco™ Newborn Calf Serum [heat inactivated], of New Zealand origin; Thermo Fisher Scientific), 5 % AA (Gibco® MEM), Non-Essential Amino Acids, 100 U/mL penicillin and 100 µg/mL streptomycin (Thermo Fisher Scientific). Cells were grown at 37 °C in a humidified atmosphere of 5 % CO2 in air, and sub-passaged with fresh complete RPMI medium every three days. The cell line was regularly checked to be mycoplasma free using the VenorGeM Mycoplasma detection kit (Cambio Ltd, UK) according to manufacturer’s instructions.
Stock solutions of Arecoline (100 mg/ml of methanol), MNPA (2.5 mg/ml of methanol), MNPN (20 mg/ml of 100 % ethyl acetate) and PMA((Phorbol-12-Myrsitate-13-Acetate (2 mg/ml of DMSO)) were made, aliquoted and stored at -20oC. To assess the toxicity of these preparations, THP-1 cells was incubated with 10-fold serial dilutions of arecoline, MNPN as well as MNPA in culture media. The cells were then observed under phase contrast light microscopy for up to 6 days. We then selected the highest non-toxic concentration listed below for our experiments. 1 × 106 THP1 cells were treated with either 50 μg/ml Arecoline, 25 μg/ml MNPA, 200 μg/ml MNPN or 200 ng/ml PMA as a positive control in 6-well plates and cells were harvested after 6 h, 24 h and 48 h of treatment. Methanol and ethyl acetate were used as negative controls. Three independent experiments were performed for each exposure. Cell viability and proliferation was assessed by phase contrast light microscopy and using a MTT assay. All chemicals were purchased from Sigma-Aldrich.
RNA extraction and RT-qPCR for gene expression of human TNFa, IL-6 and IL-8 analysis
Prior to whole genome sequencing we sought to confirm the biological viability of the THP1 cells studied following their incubation with arecoline, MNPA or MNPN by measuring their cytokine responses. Total RNA was extracted from treated cells using QIAGEN RNeasy Kit according to manufacturer’s instructions. cDNA was synthesized using 1 µg of the extracted RNA with an Oligo (dT) primer using a SuperScript® IV First-Strand Synthesis System (Thermo Fisher Scientific) as follows: primer annealing at 65 °C for 5 min; RNA reverse transcriptase at 50 °C for 1 h 10 min and at 70 °C for 15 min. The cDNA was used as a template to determine the expression level of human TNFa, IL-6, IL-8 and 18 S  treated (arecoline, MNPA or MNPN) /untreated THP1 cells. The RT-PCR was performed on StepOne™ Real-Time PCR System thermal cycler (Applied Biosystems™). Each PCR reaction consisted of 2 µl of cDNA, 2X SYBR® Green JumpStart™ Taq ReadyMix™ (Sigma-Aldrich), 0.2 µM of each forward and reverse primers (Additional Table 3). qPCR reaction conditions were: cDNA denaturation at 95 °C for 5 min, cDNA amplification at 95 °C for 15 s, primer annealing at 62 °C for 1 min for 45 cycles, then melt curves were obtained at 95 °C for 15 s, 60 °C for 1 min and a final step at 95 °C for 15 s. All target genes were normalised to 18 S RNA using the standard ΔΔCt method. Results were analysed with Thermo Fisher StepOne software v2.3. Each experiment was performed in triplicate and fold change expression level was reported relative to 18 S level.
RNA-sequencing and bioinformatics
Fragmented cDNA Sequencing libraries were generated from 100ng of Total RNA using NEBNext Ultra II with polyA isolation module (Illumina, San Diego, California, USA) according to manufacturer’s protocol. cDNA quantity and quality were evaluated using the Qubit dsDNA HS assay kit (Thermo Fisher Scientific, Erembodegem-Aalst, Belgium). Size distribution of our library was determined using an Agilent 2100 Bioanalyser. Pooled indexed paired-end libraries were sequenced on the Illumina NextSeq 500 (Illumina, USA) using the manufacturer’s instructions. Sequencing was performed at Queen Mary University of London Genome Centre core facility in the Blizard Institute London.
Sequenced reads were mapped using Kallisto  with default settings. Mean insert sizes and standard deviations were provided as input. Analysis of differential gene expression was performed using sleuth, applying a generalized linear model and utilising bootstraps on reads to estimate inferential variance. Genome-wide corrected p-values were calculated using the Bonferroni multiple testing adjustment procedure.
Functional annotation as well as pathway enrichment analyses were performed using DAVID, Reactome and Metascape (https://metascape.org/gp/index.html#/main/step1).
Candidacy of genes identified were assessed by look-ups in:
GeneCards: The human genome database (https://www.genecards.org/) to check for alias’s, gene summary (Entrez, Genecards and UniProtKB/Swiss-Prot) and mRNA expression in normal human tissues (GTEx, Illumina, BioGPS).
Type 2 diabetes knowledge portal (http://www.type2diabetesgenetics.org/). Genes considered were only those with ‘strong’ evidence for signal defined by either at least one variant within the coding sequence ± 100 kb that is associated with at least one phenotype with a p value < 5e-8 identified by a genome-wide association scan (GWAS), or at least one known variant with a missense or protein truncating mutation in the encoded protein that is associated with at least one phenotype with a p value < 5e-6.
PubMed (GAH and BJB independently) by searching biology of the identified gene and biological relevance to obesity, T2D or the metabolic syndrome before collation of results.
Availability of data and materials
All genes identified after incubation with significant expression changes after Bonferroni correction, irrespective of candidacy, are listed in Additional Tables 2 and 3. The full dataset generated and/or analysed during the current study are available in the GEO repository with an accession number of GSE179143, website: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE179143.
Type 2 diabetes
Type 1 diabetes
Hills AP, Arena R, Khunti K, Yajnik CS, Jayawardena R, Henry CJ, et al. Epidemiology and determinants of type 2 diabetes in south Asia. Lancet Diabetes Endocrinol. 2018;6(12):966–78.
Ko YC, Huang YL, Lee CH, Chen MJ, Lin LM, Tsai CC. Betel quid chewing, cigarette smoking and alcohol consumption related to oral cancer in Taiwan. J Oral Pathol Med. 1995;24(10):450–3.
Tsai JF, Chuang LY, Jeng JE, Ho MS, Hsieh MY, Lin ZY, et al. Betel quid chewing as a risk factor for hepatocellular carcinoma: a case-control study. Br J Cancer. 2001;84(5):709–13.
Wu MT, Lee YC, Chen CJ, Yang PW, Lee CJ, Wu DC, et al. Risk of betel chewing for oesophageal cancer in Taiwan. Br J Cancer. 2001;85(5):658–60.
Yamada T, Hara K, Kadowaki T. Chewing betel quid and the risk of metabolic disease, cardiovascular disease, and all-cause mortality: a meta-analysis. PLoS One. 2013;8(8):e70679.
Tung TH, Chiu YH, Chen LS, Wu HM, Boucher BJ, Chen TH, et al. A population-based study of the association between areca nut chewing and type 2 diabetes mellitus in men (Keelung Community-based Integrated Screening programme No. 2). Diabetologia. 2004;47(10):1776–81.
Yen AM, Chiu YH, Chen LS, Wu HM, Huang CC, Boucher BJ, et al. A population-based study of the association between betel-quid chewing and the metabolic syndrome in men. Am J Clin Nutr. 2006;83(5):1153–60.
Yen AM, Boucher BJ, Chiu SY, Fann JC, Chen SL, Huang KC, et al. Longer Duration and Earlier Age of Onset of Paternal Betel Chewing and Smoking Increase Metabolic Syndrome Risk in Human Offspring, Independently, in a Community-Based Screening Program in Taiwan. Circulation. 2016;134(5):392–404.
Boucher BJ, Ewen SW, Stowers JM. Betel nut (Areca catechu) consumption and the induction of glucose intolerance in adult CD1 mice and in their F1 and F2 offspring. Diabetologia. 1994;37(1):49–55.
Boucher BJ, Mannan N. Metabolic effects of the consumption of Areca catechu. Addict Biol. 2002;7(1):103–10.
IARC Working Group on the Evaluation of Carcinogenic Risks to Humans. Betel-quid and areca-nut chewing and some areca-nut derived nitrosamines. IARC Monogr Eval Carcinog Risks Hum. 2004;85:160–2.
Karam JH, Lewitt PA, Young CW, Nowlain RE, Frankel BJ, Fujiya H, et al. Insulinopenic diabetes after rodenticide (Vacor) ingestion: a unique model of acquired diabetes in man. Diabetes. 1980;29(12):971–8.
Zhang M, Lv XY, Li J, Xu ZG, Chen L. The characterization of high-fat diet and multiple low-dose streptozotocin induced type 2 diabetes rat model. Exp Diabetes Res. 2008;2008:704045.
Portha B, Giroix MH, Cros JC, Picon L. Diabetogenic effect of N-nitrosomethylurea and N-nitrosomethylurethane in the adult rat. Ann Nutr Aliment. 1980;34(5–6):1143–51.
Garg A, Chaturvedi P, Gupta PC. A review of the systemic adverse effects of areca nut or betel nut. Indian J Med Paediatr Oncol. 2014;35(1):3–9.
Ogunkolade WB, Boucher BJ, Bustin SA, Burrin JM, Noonan K, Mannan N, et al. Vitamin D metabolism in peripheral blood mononuclear cells is influenced by chewing “betel nut” (Areca catechu) and vitamin D status. J Clin Endocrinol Metab. 2006;91(7):2612–7.
Kasikara C, Doran AC, Cai B, Tabas I. The role of non-resolving inflammation in atherosclerosis. J Clin Invest. 2018;128(7):2713–23.
Saltiel AR, Olefsky JM. Inflammatory mechanisms linking obesity and metabolic disease. J Clin Invest. 2017;127(1):1–4.
Lumeng CN, DelProposto JB, Westcott DJ, Saltiel AR. Phenotypic switching of adipose tissue macrophages with obesity is generated by spatiotemporal differences in macrophage subtypes. Diabetes. 2008;57(12):3239–46.
Westcott DJ, Delproposto JB, Geletka LM, Wang T, Singer K, Saltiel AR, et al. MGL1 promotes adipose tissue inflammation and insulin resistance by regulating 7/4hi monocytes in obesity. J Exp Med. 2009;206(13):3143–56.
Waeber G, Delplanque J, Bonny C, Mooser V, Steinmann M, Widmann C, et al. The gene MAPK8IP1, encoding islet-brain-1, is a candidate for type 2 diabetes. Nat Genet. 2000;24(3):291–5.
Willer CJ, Speliotes EK, Loos RJ, Li S, Lindgren CM, Heid IM, et al. Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet. 2009;41(1):25–34.
Joo Y, Kim H, Lee S, Lee S. Neuronal growth regulator 1-deficient mice show increased adiposity and decreased muscle mass. Int J Obes (Lond). 2019;43(9):1769–82.
Palming J, Sjoholm K, Jernas M, Lystig TC, Gummesson A, Romeo S, et al. The expression of NAD(P)H:quinone oxidoreductase 1 is high in human adipose tissue, reduced by weight loss, and correlates with adiposity, insulin sensitivity, and markers of liver dysfunction. J Clin Endocrinol Metab. 2007;92(6):2346–52.
Hashimoto O, Funaba M, Sekiyama K, Doi S, Shindo D, Satoh R, et al. Activin E Controls Energy Homeostasis in Both Brown and White Adipose Tissues as a Hepatokine. Cell Rep. 2018;25(5):1193–203.
Coppo M, Chinenov Y, Sacta MA, Rogatsky I. The transcriptional coregulator GRIP1 controls macrophage polarization and metabolic homeostasis. Nat Commun. 2016;7:12254.
Seo JA, Kang MC, Yang WM, Hwang WM, Kim SS, Hong SH, et al. Apolipoprotein J is a hepatokine regulating muscle glucose metabolism and insulin sensitivity. Nat Commun. 2020;11(1):2024.
Seo JA, Kang MC, Ciaraldi TP, Kim SS, Park KS, Choe C, et al. Circulating ApoJ is closely associated with insulin resistance in human subjects. Metabolism. 2018;78:155–66.
Al-Daghri NM, Manousopoulou A, Alokail MS, Yakout S, Alenad A, Garay-Baquero DJ, et al. Sex-specific correlation of IGFBP-2 and IGFBP-3 with vitamin D status in adults with obesity: a cross-sectional serum proteomics study. Nutr Diabetes. 2018;8(1):54.
Hjortebjerg R, Laugesen E, Hoyem P, Oxvig C, Stausbol-Gron B, Knudsen ST, et al. The IGF system in patients with type 2 diabetes: associations with markers of cardiovascular target organ damage. Eur J Endocrinol. 2017;176(5):521–31.
Curran AM, Scott-Boyer MP, Kaput J, Ryan MF, Drummond E, Gibney ER, et al. A proteomic signature that reflects pancreatic beta-cell function. PLoS One. 2018;13(8):e0202727.
Chanput W, Mes JJ, Wichers HJ. THP-1 cell line: an in vitro cell model for immune modulation approach. Int Immunopharmacol. 2014;23(1):37–45.
Bosshart H, Heinzelmann M. THP-1 cells as a model for human monocytes. Ann Transl Med. 2016;4(21):438.
Wu IC, Chen PH, Wang CJ, Wu DC, Tsai SM, Chao MR, et al. Quantification of blood betel quid alkaloids and urinary 8-hydroxydeoxyguanosine in humans and their association with betel chewing habits. J Anal Toxicol. 2010;34(6):325–31.
Liu YJ, Peng W, Hu MB, Xu M, Wu CJ. The pharmacology, toxicology and potential applications of arecoline: a review. Pharm Biol. 2016;54(11):2753–60.
Kumpawat K, Deb S, Ray S, Chatterjee A. Genotoxic effect of raw betel-nut extract in relation to endogenous glutathione levels and its mechanism of action in mammalian cells. Mutat Res. 2003;538(1–2):1–12.
IARC Working Group on the Evaluation of Carcinogenic Risks to Humans. Betel-quid and areca-nut chewing and some areca-nut-derived nirosamines. IARC Monogr Eval Carcinog Risks Hum. 2004;85:51–2.
Tsuchiya S, Yamabe M, Yamaguchi Y, Kobayashi Y, Konno T, Tada K. Establishment and characterization of a human acute monocytic leukemia cell line (THP-1). Int J Cancer. 1980;26(2):171–6.
Nhu QM, Shirey K, Teijaro JR, Farber DL, Netzel-Arnett S, Antalis TM, et al. Novel signaling interactions between proteinase-activated receptor 2 and Toll-like receptors in vitro and in vivo. Mucosal Immunol. 2010;3(1):29–39.
Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol. 2016;34(5):525–7.
The authors thank Nisha Patel (Borchukova laboratory, Blizard Institute, Queen Mary University of London) for technical support.
The project was supported by a grant from the Barts and the London Charity (grant number MRC0183). The Charity approved the funding but did not influence the planning of the studies, or the interpretation of the results; they did receive an interim and final report of the study which they approved.
Ethics approval and consent to participate
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Phase light microscopy of THP-1 cells. Individual panels are of control cells or various concentration of arecoline studied at 48 h.
MTT-Cell proliferation assay. Bar chart with 95% confidence limits plotting % cell viability by control and after incubation of arecoline/well at varying concentrations.
Betel nut compounds treatment in THP1 cells. Bar chart with 95% confidence limits plotting a time course for TNFa, IL6 and IL8 mRNA (normalised relative to 18S RNA) in response to incubation with arecoline, MNPA, MNPN and PMA.
Shared genes expressed (q < 0.05) after incubation with either arecoline or MNPA. Illustrated as a Venn diagram.
All genes identified after incubation with arecoline with significant increased expression after Bonferroni correction. Table 2. Genes identified after incubation with MNPA with significant increased expression after Bonferroni correction. Table S3. Primer sequences used in qPCR analysis https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2851245/
About this article
Cite this article
Cardosa, S.R., Ogunkolade, B.W., Lowe, R. et al. Areca catechu-(Betel-nut)-induced whole transcriptome changes in a human monocyte cell line that may have relevance to diabetes and obesity; a pilot study. BMC Endocr Disord 21, 165 (2021). https://doi.org/10.1186/s12902-021-00827-1
- Type 2 diabetes