- University College of Medical Sciences, University of Delhi, Delhi, India
Aim: TCF7L2 gene is believed to increase the risk of T2DM by its effects on insulin secretion. However, the exact mechanism of this enhanced risk is not clearly known. While TCF7L2 gene has been shown to affect lipid metabolism, these effects have remained largely unexplored in the context of diabetes risk.
Methods: Postprandial lipid responses to a standardized fat challenge test were performed in 620 Asian Indian subjects (310 with NGT and 310 with T2DM/prediabetes) and compared between the risk and wild genotypes of the rs7903146 TCF7L2 gene. In 30 subjects scheduled to undergo abdominal surgery (10 each with NGT, Prediabetes and T2DM), adipocyte TCF7L2 gene expression was also performed by real time qPCR and confirmed by protein expression in western blot.
Results: T allele of rs7903146 TCF7L2 gene was confirmed as the risk allele for T2DM (OR=1.8(1.2-2.74), p=0.005). TT+CT genotypes of rs7903146 TCF7L2 gene showed significantly higher 4hrTg (p<0.01), TgAUC (p<0.01), peakTg (p<0.01) as well as higher postprandial plasma glucose (p=.006) levels and HOMA-IR (p=0.03) and significantly lower adiponectin levels (p=0.02) as compared to CC genotype. The expression of TCF7L2 gene in VAT was 11-fold higher in prediabetes group as compared to NGT (P<0.01) and 5.7-fold higher in T2DM group as compared to NGT group(P=0.003) and was significantly associated with PPTg and glucose levels.
Conclusion: There is significant PPTg dysmetabolism associated with the risk allele of rs7903146 polymorphism as well as adipocyte expression of TCF7L2 gene. Significant upregulation of TCF7L2 gene expression in VAT that correlates with PPTg and glycaemia is also seen in Asian Indians with glucose intolerance. Modulation of PPTg metabolism by TCF7L2 gene and the resultant PPHTg may be a novel mechanism that contributes to its diabetes risk in them.
Introduction
Transcription factor 7 like 2 (TCF7L2) gene has been shown to have the strongest association with type 2 diabetes mellitus among all diabetogenic genes and this association has been replicated in all countries of the world including India (1–6). The risk of type 2 diabetes mellitus (T2DM) increases by 50% in the presence of the risk allele of TCF7L2 gene at an individual level while the population attributable risk varies between 10-25% depending on the allele frequency (1, 7, 8). TCF7L2 gene variant has also been shown to predict future diabetes risk in subjects with impaired glucose tolerance (9, 10). Despite this strong and consistent association, the exact mechanism of action by which this gene enhances diabetes risk is still unclear. The TCF7L2 gene encodes a high mobility group (HMG) box-containing transcription factor, which is the effector of the canonical WNT signalling pathway that plays a key role in cell differentiation (11–13). TCF7L2 protein also regulates the expression of several genes involved in cell cycle process including cyclin D1 and c-myc, which control the G1 to S phase transition (12, 14). It is generally believed that its effect on insulin secretion is the primary mechanism underlying its association with T2DM. Most studies have focused on the effects of TCF7L2 gene on beta cell proliferation where it was found that altered expression of TCF7L2 gene in beta cells of the pancreas results in impaired insulin secretion (9, 10, 15) Dominant-negative mutant of TCF7L2 abolishes pro glucagon mRNA levels (16, 17). TCF7L2 gene also affects GLP1 production and gluconeogenesis (18). Very few studies have explored its effects on adipogenesis and lipid metabolism as a possible mechanism for the TCF7L2 gene associated risk of development of diabetes (19–22).
Postprandial hypertrglyceridaemia (PPHTg) may play an important role in the pathogenesis of diabetes (23–27). Recently, we have shown that PPHTg can predict the development of insulin resistance and diabetes (28). TCF7L2 gene is known to upregulate the WNT signalling pathway in the adipose tissue and the resultant decrease in adipogenesis can lead to PPHTg (11, 13, 14, 16, 29). In an earlier study, variants of TCF7L2 were found to modulate postprandial triglycerides levels in individuals with polyunsaturated fatty acid (PUFA) intakes of more than 7.36% (21). However, one of the studies conducted at our institute in subjects with normal glucose tolerance (NGT) did not show a significant difference in postprandial hypertriglyceridemia between the risk allele and wild type variants of this gene in first degree relatives of diabetes patients who had NGT and lower BMI (30). Decreased adipogenesis due to altered TCF7L2 gene (11) function may result into impaired Tg trapping and PPHTg. The resultant overflow of triglycerides to ectopic sites such as liver and muscles may lead to lipotoxicity, insulin resistance and diabetes. Therefore, we aimed to study the polymorphisms of TCF7L2 gene and its expression in subcutaneous and visceral adipose tissue to find out its association with postprandial hypertriglyceridaemia and risk of diabetes in subjects with varying degree of glucose tolerance.
Methods
Ethical considerations
The study was approved by Institutional Ethics Committee of University college of medical sciences and Guru Teg Bahadur Hospital. A written informed consent was obtained from all the subjects prior to the initiation of the study which was conducted in accordance with the principles of Declaration of Helsinki (31)
Study population
It was a case control study conducted in 620 Asian Indians subjects (310 NGT, 310 subjects with glucose intolerance) between the age group of 20-60 years. The subjects with normal glucose tolerance, prediabetes and T2DM were identified on the basis of standard oral glucose tolerance tests as per the ADA criteria. Known diabetes subjects were recruited from Diabetic clinic of Department of Endocrinology. The genotype frequency of rs7903146 polymorphic form of TCF7L2 gene in the literature was used to calculate the sample size required (2). The reported genotype frequency and the required number of subjects in each group to detect a significant difference in genotyping pattern using 2/3 chi-square with 5% alpha error and 90% power came out to be a minimum of 295 for each group. The sample size calculated to study the association between postprandial Tg parameters and genotypes of TCF7L2 with 5% alpha error and 90% power came out to be 55 for each group based on existing literature. However, we have considered the higher sample size so that the objectives of this study can be met. Groups were matched for age, sex and BMI. Subjects with any chronic disease such as hyperlipidaemia inherited lipoprotein metabolism disorder, pancreatic disease, obstructive biliary disease or coronary heart disease, liver or kidney disease, hypothyroidism, Cushing’s syndrome, chronic diarrhoea, and malabsorption affecting lipid metabolism were excluded. Flow chart of study design is given in Figue 1.
Fat challenge test and Biochemical estimations
A standardized 75 gm oral glucose tolerance test was performed in all the subjects to characterize them into different categories of glucose intolerance i.e., NGT and subjects with glucose intolerance on the basis of ADA criteria (32). Fasting blood samples was also collected for estimation of insulin levels. One week post Oral glucose tolerance test (OGTT), all subjects were called after overnight fasting for a standardized oral fat challenge test (33). Fasting blood samples were collected for the measurement of lipids, adiponectin and DNA extraction. A standardized fatty meal containing 729 kcal/m2 body surface area of calories, 62.5 gm of fat, 24.75 gm of carbohydrates, 5.3 gm of protein and 240 mg of cholesterol, that was made from whipped cream, banana and sugar was given to the subjects to be consumed over a period of 10 to 15 min. Blood samples were collected for estimation of lipids after the fat challenge 2 hourly upto 4 hours in all subjects and upto 8 hours in 2/3rd of them.
Plasma glucose levels, Glucose Oxidase - Peroxidase (GOD POD method; Randox UK catalogue number: GL8308), serum insulin (IRMA method using 125I isotopes; Backmencoulter, catalouge no.-IM3210) and HbA1C (Biorad D-10 catalouge no: 2200101) were measured. Serum triglyceride, Total cholesterol and HDL were measured using Randox UK kits (Catalouge no; Tg: 120211, Tchol: 120194, HDL: CH2673) and adiponectin levels were estimated using Bio Vendor Elisa kits (Catalouge No. – RD194015200R). All the investigation procedures were initially standardized by repeatedly running the tests on same set of samples under a range of conditions. Homeostatic model assessment of insulin resistance (HOMA-IR) was calculated using fasting insulin and fasting plasma glucose levels.
DNA extraction and genotyping
DNA was extracted from whole blood collected in EDTA tube by spin protocol using Qiagen kit (QIAamp DNA Mini kit, catalouge no. 51306) The rs7903146 polymorphism of TCF7L2 gene was genotyped in 620 subjects using the polymerase chain reaction- restriction fragment length polymorphism method using following primers: forward primer 5’AAGAGAAGATTCCTTTTTAAATGGTG3’Reverse primer 3’CCTCATACGGCAATTAAATTATACA5’ (genei). A final reaction volume of 25 µl for the PCR was made which contained 100 ng of whole DNA, PCR master mix Thermo scientific #k0172 (0.05 U/µL Taq DNA polymerase, 4 mM Mgcl2, 0.4 mM of dNTP) and nuclease free water to make up the volume. The PCR was carried out on Biorad T100 model under following condition 95°C: 3 min 95°C: 30 sec 55°C:40 sec 72°C:45 sec and 72°C: 5 min. HPYCHIV restriction enzyme were used for the digestion of PCR product at 65 °C for 5 min. The digested products were analysed by electrophoresis on 2% agarose gel.
Sequencing
We have further confirmed the Restriction fragment length polymorphism (RFLP) results for each genotype (TT, CC and CT) in a representative sample through sanger sequencing. Three samples of each genotype of rs7903146 polymorphic form of TCF7L2gene PCR product were randomly selected and sequenced. For sequencing PCR products were first purified with a column-based purification kit (Invitrogen™), and sequenced by the Sanger method using an Applied Biosystems 3130 Genetic Analyzer. The recognition sites of the TaaI enzyme within the sequence region were analysed with the Serial Cloner V2.5 program. The homology of DNA sequences to the published sequences was determined by using the BLAST platform on the national centre for biotechnology information (NCBI) website. The figures of the peak obtained for different genotypes are given in Supplementary Figure 1.
Collection of biopsies and Gene expression studies
Thirty subjects (10 subjects each of NGT, Pre diabetes, and T2DM) who were to undergo elective appendicectomy, inguinal hernia repair or laparoscopic cholecystectomy for gall bladder stone removal were recruited for the gene expression study if they fulfilled the inclusion and exclusion criteria and were matched for age sex and BMI. Approximately 1 cm2 biopsies of subcutaneous adipose tissue (SAT) and visceral adipose tissue (VAT) were taken during surgery. The same was washed in physiological saline and immediately transferred in liquid nitrogen and carried to the biochemistry lab. Biopsies were homogenized and fixed in Trizol LS reagent (Invitrogen, Thermo scientific, catalogue no. 15596026) on the same day and were kept at −80°C until analysis. Chloroform-isopropanol method was used for the extraction of RNA. The purity of RNA sample was assessed by measuring the absorption at 230 nm, 260 nm and 280 nm using nano drop. The average yield of RNA in SAT was 384.07 ng/ul and in VAT was 1198.8 ng/ul and the integrity of RNA was checked by running the samples on 3% agarose gel. The ratio of 260:280 between 1.9-2.1 were considered for further analysis. cDNA was synthesized using iscript cDNA synthesis kit Biorad USA (catalogue no. 1708891).
TCF7L2 gene expression in subcutaneous and visceral adipose tissue was quantified by quantitative real time polymerase chain reaction (qPCR) experiment conducted on CFX Connect TM Real-Time PCR Detection System (BioRad, USA). As housekeeping genes, β actin and GAPDH were used. All the samples were run in triplicate along with negative control. The results were calculated using the comparative 2-∆∆CT method (Livak Method) manually (34). PCR amplification condition for TCF7L2 gene was: denaturation at 94°C for 3 min followed by 40 cycles of denaturation for 30 s, annealing at 55°C for 30 s, extension at 72°C, and final extension at 72°C for 9 min using the forward primer 5’AAGAGAAGATTCCTTTTTAAATGGTG3’ Reverse primer 3’CCTCATACGGCAATTAAATTATACA5’. One µg of cDNA was used uniformly along with PCR master mix (SsoFast EvaGreen Supermix, Biorad US catalogue number 1725200) containing 0.4 mM of each dNTPs and 4 mM MgCl2. The length of TCF7L2 amplicon was 267 bp and the mRNA variants were X4-XM011533842. The PCR efficiency was calculated from the slope.
RIPA lysis buffer (G Biosciences cat no. R159) was used to extract total protein from homogenized visceral adipose tissue which was then quantified by Bradford protein assay (Bedford, MA, USA). SDS PAGE (10%) was used to separate proteins which were transferred to PVDF membrane (Biorad cat no. 12-0177). The PVDF was further blocked using 5% non-fat dry milk in TBST. Primary antibody TCF7L2 (1:700) (abbkine, cat no.-Abp52252-1) of rabbit anti goat source and anti-β-actin (abbkine, cat no.-Abp50593-1) (1:1000) and secondary antibody (abbkine, cat no.-Abp 21020-2) were used. The luminescence was recorded under Chemidoc (Thermo scientific ECL Imager, USA) which was further quantified by using Image J software.
Statistical analysis
Differences between the clinical characteristics of the group were calculated by using unpaired t-test (for two groups) and one-way analysis of variance (for more than two groups). Comparisons of distribution of polymorphic forms between the groups were done by Chi-square test. Logistic regression analysis was performed to find out the odds ratio for disease risk associated with genotypes. The fold changes for expression of genes were calculated by using Livak method. Differences in gene expression between the groups were assessed with Mann Whitney test. Spearman correlation was performed to find the association of different variables with each other. Statistical analysis was performed using the software SPSS 20.0
Results
Total number of study subjects were 620 (Subjects with normal glucose tolerance =310; Subjects with glucose Intolerance=310). Mean age of study subjects was 40.67 ± 8.92 years and mean BMI was 27.98 ± 4.89 kg/m2.
Table 1 shows comparison of clinical and biochemical parameters between the study groups. All the groups were matched for age, sex and BMI shown by the absence of any significant difference of these parameters between the groups.
It was observed that the T allele of rs7903146 polymorphic form of TCF7L2 gene was significantly higher in subjects with Glucose Intolerance (44.2%) as compared to NGT groups (35.6%). There was a significant association of the T allele with glucose intolerance (OR=1.66(1.14-2.18) p=0.005) which indicates that T allele is the risk allele (Table 2). When subjects were re-analysed after further dividing those with glucose intolerance into prediabetes and T2DM, it was found that T allele remained significantly higher in the T2DM group when compared to those with NGT (OR=1.8(1.2-2.74) p=0.005) or prediabetes (OR=1.5(1.7-2.25) p=0.03). T allele remained higher in the prediabetes group when compared with those with NGT (OR=1.4(0.99-2.05) p=0.05). Further, T allele was also significantly associated with T2DM (Table 3). The genotypic distribution in our study was as per the Hardy- Weinberg equilibrium.
Table 2 - Genotypic and allelic frequencies and estimates of relative risks for the rs7903146 polymorphic form of TCF7L2 gene in subjects with and without glucose intolerance.
Table 3 - Genotypic and allelic frequencies and estimates of relative risks for the rs7903146 polymorphic form of TCF7L2 gene in subjects with normal glucose tolerance, Prediabetes and Diabetes.
Fat challenge test was carried out until 4 hours in all 620 subjects and lipid samples were collected till 8 hours in 406 subjects of them. Postprandial triglycerides (PPTg) levels and glycaemic parameters were compared between risk genotype and wild genotype of rs7903146 polymorphic form of TCF7L2 gene. It was found that TT+CT genotypes (T allele) of rs7903146 polymorphic form of TCF7L2 gene showed significantly higher Tg 4hr (p<0.01), TgAUC (p<0.01) and peakTg (p<0.01) levels as well as higher postprandial plasma glucose levels (p=0.006) and HOMA-IR(p=0.03), and significantly lower adiponectin levels(p=0.02) as compared to CC genotype (Table 4). Figure 2 shows that there was a significant positive correlation of TgAUC with HOMA-IR (r=0.22, p<0.01) and HbA1c (r=0.21, p<0.01) in risk genotype subjects but not in subjects with wild genotype. However, fasting plasma glucose and postprandial plasma glucose correlated significantly with TgAUC in both the groups.
Table 4 - Comparison of postprandial triglyceride, glycaemic and insulin resistance parameters between Risk and wild type genotypes of rs7903146 polymorphic form of TCF7L2 gene in all the subjects.
Figure 2 Correlation of risk and wild genotypes of rs7903146 polymorphic form of TCF7L2 gene with TgAUC, glycaemic and insulin resistance parameters. PP glucose, postprandial plasma glucose; TgAUC, triglyceride area under curve; HOMA-IR , Homeostatic model assessment of insulin resistance; statical analysis used: Pearson correlation.
Adipocyte expression of TCF7L2 gene in VAT and SAT is shown in Figure 3. The expression of TCF7L2 gene in VAT was 11-fold higher in prediabetes group as compared to NGT (P<0.01) and 5.7-fold higher in T2DM group as compared to NGT group (P<0.01). Protein expression of TCF7L2 gene in VAT by western blot technique confirmed the findings of gene expression viz. 6-fold higher TCF7L2 expression in Prediabetes group and 2.5-fold higher expression in T2DM group as compared to NGT group as shown in Figure 3.
Figure 3 TCF7L2 gene/protein expression in the three study groups. VAT= visceral adipose tissue, SAT= subcutaneous adipose tissue, NGT, normal glucose tolerance; PD, prediabetes; T2DM, type 2 diabetes mellitus; N=30; statical analysis: one way anova; **= Prediabetes vs NGT (p=≤0.01); *= T2DM vs NGT (p=≤0.01); # = Prediabetes vs NGT (p=≤0.01); ## = T2DM vs NGT (p=≤0.01). (A) Fold Change of TCF7L2 gene expression between the three groups. (B) Fold Change of TCF7L2 protein expression between the three groups. (C) Representative western blot image.
There was no significant difference in the expression of TCF7L2 gene between the three groups in SAT.
There was a significant correlation of expression of TCF7L2 gene in visceral adipose tissue with TgAUC (p<0.03), PeakTg (p<0.02) and postprandial plasma glucose levels (p<0.01) (Figure 4). However, there was no correlation of expression of TCF7L2 gene in visceral adipose tissue with fasting plasma glucose (p =0.15), HOMA-IR (p =0.12) or serum adiponectin levels (p =0.3).Also, there was no significant correlation of expression of TCF7L2 gene in subcutaneous adipose tissue with any of the PPTg, glycaemic or insulin resistance parameters.
Figure 4 Correlation of TCF7L2 gene expression in VAT with PPTg and Glycaemic parameters. PPg, postprandial plasma glucose; Tg, Triglycerides; TgAUC, triglyceride area under curve; N=30 ; statical analysis: Pearson correlation.
Discussion
Our study found that individuals with TT/CT genotypes of rs7903146 polymorphism showed significantly higher levels of postprandial triglycerides, fasting and postprandial plasma glucose, and HOMA-IR and significantly lower serum adiponectin levels as compared to CC genotype. Those with TT/CT genotypes were at a significantly higher risk of having glucose intolerance as compared to CC genotypes. Adipocyte expression of TCF7L2 gene in visceral adipose tissue was significantly higher in those with glucose intolerance as compared to NGT. Furthermore, there was a significant positive correlation between TCF7L2 gene expression in visceral adipose tissue and PPTg parameters such as TgAUC, PeakTg as well as postprandial plasma glucose.
The strength of the present study is that it not only seeks to explore the association of the rs7903146 polymorphic form of TCF7L2 gene with altered postprandial triglycerides metabolism but also examines for the first time whether altered adipocyte expression of this gene could also be associated with significant postprandial hypertriglyceridemia and insulin resistance.
The association of rs7903146 polymorphic form of TCF7L2 gene with enhanced risk of diabetes has been consistently demonstrated from different populations (4, 35, 36) in the world including Asian Indians (2, 3, 37). However, the underlying mechanisms are not fully understood. It is believed that these effects are secondary to an impairment in insulin secretion (10, 15). The role of TCF7L2 gene in modulating lipid metabolism is a less well studied area although these are mentioned in literature (19–22, 38, 39)
Our previous studies showed that postprandial hypertriglyceridemia is seen early in the natural history of diabetes and also in first degree relatives of subjects with prediabetes (33) and could lead to the development of prediabetes and diabetes (40). This prompted us to ascertain if postprandial triglycerides dysmetabolism could be an important and novel mechanism by which TCF7L2 associated risk of diabetes is mediated.
The current study demonstrated the association of TT/CT genotypes of the rs7903146 polymorphic form of TCF7L2 gene with PPTg levels in a wide range of subjects from both sexes aged between 20-60 years and included those with all classes of glucose tolerance. This would suggest a consistent and robust association between TCF7L2 gene with PPTg metabolism. The significantly higher PPTg levels in subjects with the T allele who also displayed greater insulin resistance and a higher risk of glucose intolerance clearly indicate that the risk allele of this gene is associated with PPTg dysmetabolism. After thorough research of the literature, only 5-6 studies could be identified which investigated the influence of TCF7L2 gene with respect to postprandial lipemia (19–21, 30, 38, 41). Association of T allele of rs7903146 polymorphic form of TCF7L2 gene and postprandial hypertriglyceridemia was reported by Martinez et al. in young and elder healthy males (19). However, these researchers could not demonstrate altered PPTg metabolism in subjects with metabolic syndrome. In another study, Musso G et al. demonstrated that TT/CT rs7903146 polymorphic form of TCF7L2 gene was involved in modulation of postprandial lipids and displayed higher postprandial lipoprotein levels in patients with NASH (41). Also, healthy subjects with T allele with PUFA intake of more than 7.36% were shown to have higher VLDL concentrations and postprandial triglyceride lipoproteins (21). However, there was no association of postprandial triglyceride levels with risk allele of TCF7L2 gene in healthy first degree relatives of diabetes patients who had NGT and a BMI < 25 Kg/m2 (22). Each of these studies represented specific groups of subjects mostly healthy individuals without diabetes or prediabetes unlike in our study which included subjects with varying degree of glucose tolerance.
Logistic regression analysis revealed that the risk of T2DM was significantly higher in subjects with TT genotype of rs7903146 polymorphic form of TCF7L2 gene as compared to CC genotype which confirms T allele as the risk allele for this polymorphism as reported earlier (2, 3, 37). TT/CT genotypes of the rs7903146 polymorphic form of TCF7L2 gene are also associated with higher HOMA-IR and lower serum adiponectin levels both of which are markers of higher insulin resistance. Similar results were shown in an earlier study (42). Taken together, our findings with respect to rs7903146 polymorphism indicate that postprandial dysmetabolism could be an important contributor to TCF7L2 associated insulin resistance.
Whether TCF7L2 gene associated postprandial lipid dysmetabolism leads to insulin resistance or TCF7L2 associated insulin resistance results is postprandial lipid dysmetabolism cannot be said with certainty from the current study. Association of the risk allele of TCF7L2 gene with insulin resistance (43) and post meal insulin resistance (44) has been reported before. However, when read with previous studies that have suggested that postprandial lipemia is an early phenomenon in diabetes, prediabetes (27), first degree relatives of T2DM (24, 27) as well as our previous study in animals (28), it would appear likely that it is TCF7L2 gene associated PPHTg that results in insulin resistance and not vice versa.
We also found a significant positive correlation between expression of TCF7L2 gene in visceral adipose tissue and PPTg parameters such as TgAUC and PeakTg levels in subjects who had varying levels of glucose intolerance. This is the first demonstration of a significant association between TCF7L2 expression in visceral adipose tissue and PPTg dysmetabolism confirming the findings from our polymorphism studies and points to an important role of TCF7L2 gene in postprandial hypertriglyceridemia. We did not find any study in the literature regarding TCF7L2 gene expression in adipose tissue in the context of postprandial lipemia.
We found a significantly higher expression of TCF7L2 gene in VAT in patients with T2DM while there was no significant change in SAT. This finding in VAT is both novel as well as interesting and may be related to the unique pattern of body fat distribution described in Asian Indians. Results of previous studies on adipose tissue TCF7L2 gene expression in T2DM and impaired glucose metabolism have been conflicting and have reported both upregulation (45) as well as downregulation (46) of TCF7L2 gene in adipose tissue. Few studies reported no significant change in TCF7L2 expression also (47) for expression studies (46, 48–50). The increased adipose tissue expression of TCF7L2 gene in our study in subjects with glucose intolerance may reflect the known phenotypic differences in adipose tissue morphology and distribution in both these populations. It is well known that Asian Indians have significantly higher VAT and significantly lower SAT when compared to white Caucasians (51–53). This, phenotype which has been described as the “thin fat” Indian phenotype, is associated with (54) normal body weight in many cases, and is believed to be a form of partial lipodystrophy that facilitates ectopic lipid accumulation in visceral adipose depots. Central obesity and higher VAT depots are believed to be central to the development of IR and T2DM in south Asians including Asian Indians (51).
Recent genetic and transcriptomic studies using gene set enrichment analysis from India have shown that Adipocyte is a key player in the pathogenesis of insulin resistance and T2DM in Indians (55–57). Asian Indians showed differentially expressed lipodystrophy genes in peripheral subcutaneous tissue of type 2 diabetic patients (55) significantly greater expression of several genes affecting adipogenesis and inflammation, particularly in VAT and visceral adipocyte hypertrophy in diabetic subjects compared to controls whereas in SAT the difference in adipocyte size was not significant (55). Inadequate activation of adipogenesis results in a lower capacity of SAT to store fat and any response to chronic overnutrition in this scenario is by adipocyte hypertrophy in SAT and VAT (58). It has therefore been suggested that genetic predisposition for T2DM in south Asians is associated with restricted adipogenesis and hypertrophic obesity (59). It is well known that TCF7L2 gene activates the WNT pathway through the effects of beta-catenin in the visceral adipose tissue thereby inhibiting adipogenesis (12, 60). It is therefore possible that the TCF7L2-related genetic predisposition to T2DM in Asian Indians is at least partly mediated by its effects on the WNT pathway. The finding of increased TCF7L2 expression in VAT in our study supports this hypothesis. However, it will not be possible from this study to say that it is only the overexpression of TCF7L2 gene in visceral adipose tissue that results in insulin resistance. It could also be possible that insulin resistance in these subjects upregulates this gene. Only longitudinal cohort studies can answer this with clarity.
Recent evidence reveals a more complex regulatory role for TCF7L2 gene and its adipose tissue expression in mediating WNT pathway related metabolic effects that are associated with Insulin resistance and T2DM. Chen et al. have clearly demonstrated that there could be important differences depending on whether TCF7L2 expression is in preadipocytes or in mature adipocyte. In preadipocytes with high levels of beta-catenin, TCF7L2 overexpression is associated with greater beta-catenin binding that facilitates WNT activation and inhibition of adipogenesis whereas in mature adipocytes, with low levels of beta-catenin, it is TCF7L2 gene downregulation that activates WNT pathway leading to similar metabolic effects such as adipocyte hypertrophy, increased fat deposition in the liver and insulin resistance (11). It is very well possible that in Caucasians, who have a much higher subcutaneous fat depot (SAT) and an abundance of mature adipocytes, it is the down regulation of TCF7L2 gene in mature adipocytes that plays a dominant role in causing insulin resistance and glucose intolerance. On the contrary, in Asian Indians, the paucity of mature adipocyte (61) secondary to the partial lipodystrophic phenotype would make any effect of TCF7L2 gene under expression in these adipocytes far less significant. Hence, at least in Asian Indians the overexpression of TCF7L2 gene in preadipocytes is likely to be having a dominant role in the development of adipocyte hypertrophy and all its consequences leading to T2DM particularly in the face of chronic high energy intake. Clearly, more research is required to understand the precise role of TCF7L2 and its adipose tissue expression in adipogenesis and adipocyte function in relation to glucose intolerance and T2DM.
Our study revealed a significant association of TCF7L2 expression in VAT with all PPTg and glycaemic parameters. Activation of the canonical WNT pathway in VAT secondary to upregulation of TCF7L2 gene can also result in impairment of Tg trapping by the visceral adipose tissue and consequent PPHTg. The resultant lipid spill over and ectopic deposition of fat in liver, muscles and pancreas can lead to insulin resistance and type 2 diabetes (49, 62). We propose that dysregulation of TCF7L2 gene can result in postprandial hypertriglyceridemia through its actions on the WNT pathway in VAT and could thereby contribute to the associated diabetes risk. This could be a novel mechanism by which TCF2L2 gene enhances the risk of T2DM
Limitations
There are a few limitations of this study. In the polymorphism study, only rs7903146 SNP of TCF7L2 gene was explored as it was most commonly associated with T2DM. Examining other SNPs could have provided additional information on the association studied. We could analyze only representative samples by Sanger sequencing due to cost limitations; it would have been better to perform sequencing in all the samples. We did not use microarray for gene expression experiments which could have further enhanced the accuracy of the results. However real-time PCR is considered a valid method to perform gene expression studies. Lastly our findings need to be validated through prospective cohort studies.
Conclusion
There is significant PPTg dysmetabolism associated with the risk allele of rs7903146 polymorphism as well as adipocyte expression of TCF7L2 gene. Significant upregulation of TCF7L2 gene expression in VAT that correlates with PPTg and glycaemia is also seen in Asian Indians with glucose intolerance. Modulation of PPTg metabolism by TCF7L2 gene and the resultant PPHTg may be a novel mechanism that contributes to its diabetes risk in them.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: Gene Bank data repository with the accession number 2622729.
Ethics statement
The studies involving human participants were reviewed and approved by Institutional Ethics Committee of University college of medical sciences & Guru Teg Bahadur Hospital. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author contributions
BM: Contributed in conception and design, acquisition ofdata, data analysis, interpretation of data, drafting the article.VM: acquisition of data, writing work. MA: acquisition of data,data analysis. BB: contributed in conception and design, revising article. VA: contributed in conception and design, acquisition of data. SM: Contributed in conception and design, data analysis,interpretation of data, revising article and final approval of theversion to be published.
Funding
This work was supported by Indian Council of Medical Research, New Delhi. (Grant No. 5/4/5-3/Diab-16-NCD-II).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.973718/full#supplementary-material
References
1. Tong Y, Lin Y, Zhang Y, Yang J, Zhang Y, Liu H, et al. Association between TCF7L2 gene polymorphisms and susceptibility to type 2 diabetes mellitus: A large human genome epidemiology (HuGE) review and meta-analysis. BMC Med Genet (2009) 19:1–25. doi: 10.1186/1471-2350-10-15
2. Chauhan G, Spurgeon CJ, Tabassum R, Bhaskar S, Kulkarni SR, Mahajan A, et al. Impact of common variants of PPARG, KCNJ11, TCF7L2, SLC30A8, HHEX, CDKN2A, IGF2BP2, and CDKAL1 on the risk of type 2 diabetes in 5,164 indians. Diabetes. (2010) 59(8):2068–74. doi: 10.2337/db09-1386
3. Chandak GR, Janipalli CS, Bhaskar S, Kulkarni SR, Mohankrishna P, Hattersley AT, et al. Common variants in the TCF7L2 gene are strongly associated with type 2 diabetes mellitus in the Indian population. Diabetologia (2007) 59:2068–74. doi: 10.1007/s00125-006-0502-2
4. Grant SFA, Thorleifsson G, Reynisdottir I, Benediktsson R, Manolescu A, Sainz J, et al. Variant of transcription factor 7-like 2 (TCF7L2) gene confers risk of type 2 diabetes. Nat Genet (2006) 38(3):320–3. doi: 10.1038/ng1732
5. Gamboa-Meléndez MA, Huerta-Chagoya A, Moreno-Macías H, Vázquez-Cárdenas P, Ordóñez-Sánchez ML, Rodríguez-Guillén R, et al. Contribution of common genetic variation to the risk of type 2 diabetes in the Mexican mestizo population. Diabetes (2012) 50:63–7. doi: 10.2337/db11-0550
6. Bodhini D, Radha V, Dhar M, Narayani N, Mohan V. The rs12255372(G/T) and rs7903146(C/T) polymorphisms of the TCF7L2 gene are associated with type 2 diabetes mellitus in Asian indians. Metabolism. (2007) 56(9):1174–8. doi: 10.1016/j.metabol.2007.04.012
7. Zeggini E, McCarthy MI. TCF7L2: The biggest story in diabetes genetics since HLA? Diabetologia (2007) 50(1):1–4. doi: 10.1007/s00125-006-0507-x
8. Cauchi S, Froguel P. TCF7L2 genetic defect and type 2 diabetes. Curr Diabetes Rep (2008) 8(2):1–4. doi: 10.1007/s11892-008-0026-x
9. Florez JC, Jablonski KA, Bayley N, Pollin TI, de Bakker PIW, Shuldiner AR, et al. TCF7L2 polymorphisms and progression to diabetes in the diabetes prevention program. N Engl J Med (2006) 355(3):241–50. doi: 10.1056/NEJMoa062418
10. Lyssenko V, Lupi R, Marchetti P, Del Guerra S, Orho-Melander M, Almgren P, et al. Mechanisms by which common variants in the TCF7L2 gene increase risk of type 2 diabetes. J Clin Invest. (2007) 117(8):2155–63. doi: 10.1172/JCI30706
11. Chen X, Ayala I, Shannon C, Fourcaudot M, Acharya NK, Jenkinson CP, et al. The diabetes gene and WNT pathway effector TCF7L2 regulates adipocyte development and function. Diabetes. (2018) 67(4):554–68. doi: 10.2337/db17-0318
12. Barker N, Morin PJ, Clevers H. The yin-yang of TCF/beta-catenin signaling. Adv Cancer Res (2000) 77:1–24. doi: 10.1016/S0065-230X(08)60783-6
14. Tian X, Liu Z, Niu B, Zhang J, Tan TK, Lee SR, et al. E-cadherin/β-catenin complex and the epithelial barrier. J Biomed Biotechnol (2011) 2011:1–6. doi: 10.1155/2011/567305
15. Zhou Y, Park SY, Su J, Bailey K, Ottosson-Laakso E, Shcherbina L, et al. TCF7L2 is a master regulator of insulin production and processing. Hum Mol Genet (2014) 23(24):6419–31. doi: 10.1101/003202
16. Shao W, Wang D, Chiang YT, Ip W, Zhu L, Xu F, et al. The WNT signaling pathway effector TCF7L2 controls gut and brain proglucagon gene expression and glucose homeostasis. Diabetes (2013) 62(3):789–800. doi: 10.2337/db12-0365
17. Yi F, Brubaker PL, Jin T. TCF-4 mediates cell type-specific regulation of proglucagon gene expression by beta-catenin and glycogen synthase kinase-3beta. J Biol Chem (2005) 280(2):1457–64. doi: 10.1074/jbc.M411487200
18. Boj SF, van Es JH, Huch M, Li VSW, Jose A, Hatzis P, et al. Diabetes risk gene and WNT effector Tcf7l2/TCF4 controls hepatic response to perinatal and adult metabolic demand. Cell. (2012) 151(7):1595–607. doi: 10.1016/j.cell.2012.10.053
19. Perez-Martinez P, Perez-Caballero AI, Garcia-Rios A, Yubero-Serrano EM, Camargo A, Gomez-Luna MJ, et al. Effects of rs7903146 variation in the Tcf7l2 gene in the lipid metabolism of three different populations. PloS One (2012) 7(8):1–8. doi: 10.1371/journal.pone.0043390
20. Engelbrechtsen L, Hansen TH, Mahendran Y, Pyl P, Andersson E, Jonsson A, et al. Homozygous carriers of the TCF7L2 rs7903146 T-allele show altered postprandial response in triglycerides and triglyceride-rich lipoproteins. Sci Rep (2017) 7:1–8. doi: 10.1038/srep43128
21. Warodomwichit D, Arnett DK, Kabagambe EK, Tsai MY, Hixson JE, Straka RJ, et al. Polyunsaturated fatty acids modulate the effect of TCF7L2 gene variants on postprandial lipemia. J Nutr (2009) 139(3):439–46. doi: 10.3945/jn.108.096461
22. Mishra BK, Velmurugan M, Gambhir JK, Madhu SV. Postprandial lipemia and its relation to TCF7L2 gene polymorphisms in normoglycemic first-degree relatives of type 2 diabetes patients. international journal of diabetes in developing countries. International journal of diabetes in developing countries (2019) 39(2):268–72. doi: 10.1007/s13410-018-0678-2
23. Mingrone G, Henriksen FL, Greco AV, Krogh LN, Capristo E, Gastaldelli A, et al. Triglyceride-induced diabetes associated with familial lipoprotein lipase deficiency. Diabetes. (1999) 48(6):1258–63. doi: 10.2337/diabetes.48.6.1258
24. Axelsen M, Smith U, Eriksson JW, Taskinen MR, Jansson PA. Postprandial hypertriglyceridemia and insulin resistance in normoglycemic first-degree relatives of patients with type 2 diabetes. Ann Intern Med (1999) 131(1):27–31. doi: 10.7326/0003-4819-131-1-199907060-00006
25. Mori Y, Nakagiri H, Kondo H, Murase T, Tokimitsu I, Tajima N. Dietary diacylglycerol reduces postprandial hyperlipidemia and ameliorates glucose intolerance in otsuka long-Evans tokushima fatty (OLETF) rats. Nutrition (2005) 21:933–9. doi: 10.1016/j.nut.2005.01.009
26. Kawano K, Hirashima T, Mori S, Natori T. OLETF (Otsuka long-Evans tokushima fatty) rat: a new NIDDM rat strain. Diabetes Res Clin Pract (1994) 24:s317–20. doi: 10.1016/0168-8227(94)90269-0
27. Madhu S, Sinha B, Aslam M, Mehrotra G, Dwivedi S. Postprandial triglyceride responses and endothelial function in prediabetic first-degree relatives of patients with diabetes. J Clin Lipidol (2017) 11(6):1415–20. doi: 10.1016/j.jacl.2017.08.001
28. Aslam M, Aggarwal S, Sharma KK, Galav V. Postprandial hypertriglyceridemia predicts development of insulin resistance glucose intolerance and type 2 diabetes. Plos one (2016) 25(11(1):1–15. doi: 10.1371/journal.pone.0145730
29. Smith U, Kahn BB. Adipose tissue regulates insulin sensitivity: role of adipogenesis, de novo lipogenesis and novel lipids. J Intern Med (2016) 280(5):465–75. doi: 10.1111/joim.12540
30. Mishra B K, Velmurugan M, Gambhir J, Madhu SV. Postprandial lipemia and its relation to TCF7L2 gene polymorphisms in normoglycemic first-degree relatives of type 2 diabetes patients. Int J Diabetes Developing Countries. (2018) 39(2):268–72. doi: 10.1007/s13410-018-0678-2
31. Li SJ, Wu YY, Li W, Wang SJ, Fan YM. Ultrastructural observation in a case of mucinous nevus. JDDG - J German Soc Dermatol (2018) 16(6):778–80. doi: 10.1111/ddg.13528
32. Genuth SM, Palmer JP, Nathan DM. 2.Classification and diagnosis of diabetes. Diabetes Care (2017) 40(Supplement 1):S11 LP–S24. doi: 10.2337/dc17-S005
33. Madhu SV, Kant S, Srivastava S, Kant R, Sharma SB, Bhadoria DP. Postprandial lipaemia in patients with impaired fasting glucose, impaired glucose tolerance and diabetes mellitus. Diabetes Res Clin Pract (2008) 80(3):380–5. doi: 10.1016/j.diabres.2008.01.016
34. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-delta delta C(T)) method. Methods. (2001) 25(4):402–8. doi: 10.1006/meth.2001.1262
35. Cauchi S, Meyre D, Dina C, Choquet H, Samson C, Gallina S, et al. Transcription factor TCF7L2 genetic study in the French population: Expression in human β-cells and adipose tissue and strong association with type 2 diabetes. Diabetes. (2006) 55(10):2903–8. doi: 10.2337/db06-0474
36. Miyake K, Horikawa Y, Hara K, Yasuda K, Osawa H, Furuta H, et al. Association of TCF7L2 polymorphisms with susceptibility to type 2 diabetes in 4,087 Japanese subjects. J Hum Genet (2008) 19:1–25. doi: 10.1007/s10038-007-0231-5
37. Bodhini D, Radha V, Dhar M, Narayani N, Mohan V. The rs12255372(G/T) and rs7903146(C/T) polymorphisms of the TCF7L2 gene are associated with type 2 diabetes mellitus in Asian indians. metabolism: Clinical and experimental. Metabolism (2007) 56(9):1174–8. doi: 10.1016/j.metabol.2007.04.012
38. Perez-Martinez P, Delgado-Lista J, Perez-Jimenez F, Lopez-Miranda J. Update on genetics of postprandial lipemia. Atheroscler Suppl. (2010) 11(1):39–43. doi: 10.1016/j.atherosclerosissup.2010.03.002
39. Sanghera DK, Nath SK, Ortega L, Gambarelli M, Kim-Howard X, Singh JR, et al. TCF7L2 polymorphisms are associated with type 2 diabetes in khatri sikhs from north India: Genetic variation affects lipid levels. Ann Hum Genet (2008) 53(2):174–80. doi: 10.1111/j.1469-1809.2008.00443.x
40. Kumar V, Madhu SV, Singh G, Gambhir JK. Post-prandial hypertriglyceridemia in patients with type 2 diabetes mellitus with and without macrovascular disease. J Assoc Physicians India. (2010) 58:603–7.
41. Musso G, Gambino R, Pacini G, Pagano G, Durazzo M, Cassader M. Transcription factor 7-like 2 polymorphism modulates glucose and lipid homeostasis, adipokine profile, and hepatocyte apoptosis in NASH. Hepatology (2009) 49(2):426–35. doi: 10.1002/hep.22659
42. Zhang L, Zhang M, Wang JJ, Wang CJ, Ren YC, Wang BY, et al. Association of TCF7L2 and GCG gene variants with insulin secretion, insulin resistance, and obesity in new-onset diabetes. BioMed Environ Sci (2016) 29(11):814–7. doi: 10.3967/bes2016.108
43. Damcott CM, Pollin TI, Reinhart LJ, Ott SH, Shen H, Silver KD, et al. Polymorphisms in the transcription factor 7-like 2 (TCF7L2) gene are associated with type 2 diabetes in the Amish: Replication and evidence for a role in both insulin secretion and insulin resistance. Diabetes. (2006) 55(9):2654–9. doi: 10.2337/db06-0338
44. Ferreira MC, da Silva MER, Fukui RT, Arruda-Marques MDC, dos Santos RF. Tcf7l2 correlation in both insulin secretion and postprandial insulin sensitivity. Diabetol Metab Syndrome. (2018) 10(1):1–7. doi: 10.1186/s13098-018-0338-1
45. Ahlzén M, Johansson LE, Cervin C, Tornqvist H, Groop L, Ridderstråle M. Expression of the transcription factor 7-like 2 gene (TCF7L2) in human adipocytes is down regulated by insulin. Biochem Biophys Res Commun (2008) 370(1):49–52. doi: 10.1016/j.bbrc.2008.03.006
46. Huertas-Vazquez A, Plaisier C, Weissglas-Volkov D, Sinsheimer J, Canizales-Quinteros S, Cruz-Bautista I, et al. TCF7L2 is associated with high serum triacylglycerol and differentially expressed in adipose tissue in families with familial combined hyperlipidaemia. Diabetologia. (2008) 51(1):62–9. doi: 10.1007/s00125-007-0850-6
47. Kovacs P, Berndt J, Ruschke K, Kloting N, Schon MR, Korner A, et al. TCF7L2 gene expression in human visceral and subcutaneous adipose tissue is differentially regulated but not associated with type 2 diabetes mellitus. Metabolism. (2008) 57(9):1227–31. doi: 10.1016/j.metabol.2008.04.016
48. Nguyen-Tu MS, Martinez-Sanchez A, Leclerc I, Rutter GA, da Silva Xavier G. Adipocyte-specific deletion of Tcf7l2 induces dysregulated lipid metabolism and impairs glucose tolerance in mice. Diabetologia. (2021) 64(1):129–41. doi: 10.1007/s00125-020-05292-4
49. Karczewska-Kupczewska M, Stefanowicz M, Matulewicz N, Nikołajuk A, Straczkowski M. WNT signaling genes in adipose tissue and skeletal muscle of humans with different degrees of insulin sensitivity. J Clin Endocrinol Metab (2016) 101(8):3079–87. doi: 10.1210/jc.2016-1594
50. Geoghegan G, Simcox J, Seldin MM, Parnell TJ, Stubben C, Just S, et al. Targeted deletion of Tcf7l2 in adipocytes promotes adipocyte hypertrophy and impaired glucose metabolism. Mol Metab (2019) 24:44–63. doi: 10.1016/j.molmet.2019.03.003
51. Raji A, Seely EW, Arky RA, Simonson DC. Body fat distribution and insulin resistance in healthy Asian indians and caucasians (2001). Available at: https://academic.oup.com/jcem/article/86/11/5366/2849495. doi: 10.1210/jcem.86.11.7992
52. Chandalla M, Lin P, Seenivasan T, Livingston EH, Snell PG, Grundy SM, et al. Insulin resistance and body fat distribution in south Asian men compared to Caucasian men. PloS One (2007) 2(8):1–7. doi: 10.1371/journal.pone.0000812
53. Kapoor N, Lotfaliany M, Sathish T, Thankappan KR, Thomas N, Furler J, et al. Prevalence of normal weight obesity and its associated cardio-metabolic risk factors – results from the baseline data of the kerala diabetes prevention program (KDPP). PloS One (2020) 15:1–11. doi: 10.1371/journal.pone.0237974
54. Yajnik CS. A critical evaluation of the fetal origins hypothesis and its implications for developing countries early life origins of insulin resistance and type 2 diabetes in India and other Asian countries 1. J Nutr (2004) 134:205–10. doi: 10.1093/jn/134.1.205
55. Kumar A, Tiwari P, Saxena A, Purwar N, Wahi N, Sharma B, et al. The transcriptomic evidence on the role of abdominal visceral vs. subcutaneous adipose tissue in the pathophysiology of diabetes in Asian indians indicates the involvement of both. Biomolecules. (2020) 10(9):1–23. doi: 10.3390/biom10091230
56. Saxena A, Wahi N, Kumar A, Mathur SK. Functional interactomes of genes showing association with type-2 diabetes and its intermediate phenotypic traits point towards adipo-centric mechanisms in its pathophysiology. Biomolecules (2020) 10(4):1–11. doi: 10.3390/biom10040601
57. Saxena A, Tiwari P, Wahi N, Kumar A, Mathur SK. The common pathophysiologic threads between Asian Indian diabetic’s ‘Thin fat phenotype’ and partial lipodystrophy: the peripheral adipose tissue transcriptomic evidences. Adipocyte. (2020) 9(1):253–63. doi: 10.1080/21623945.2020.1776082
58. Muir LA, Neeley CK, Meyer KA, Baker NA, Brosius AM, Washabaugh AR, et al. Adipose tissue fibrosis, hypertrophy, and hyperplasia: Correlations with diabetes in human obesity. Obesity. (2016) 24(3):597–605. doi: 10.1002/oby.21377
59. Arner P, Arner E, Hammarstedt A, Smith U. Genetic predisposition for type 2 diabetes, but not for overweight/obesity, is associated with a restricted adipogenesis. PloS One (2011) 6(4):1–5. doi: 10.1371/journal.pone.0018284
60. Karczewska-Kupczewska M, Stefanowicz M, Matulewicz N, Nikolajuk A, Straczkowski M. WNT signaling genes in adipose tissue and skeletal muscle of humans with different degrees of insulin sensitivity. J Clin Endocrinol Metab (2016) 101(8):3079–87. doi: 10.1210/jc.2016-1594
61. Bays HE. Adiposopathy: Is “Sick fat” a cardiovascular disease? J Am Coll Cardiol [Internet] (2011) 57(25):2461–73. doi: 10.1016/j.jacc.2011.02.038
Keywords: TCF7L2 gene, postprandial hypertriglyceridemia, type 2 diabetes mellitus, prediabetes, adipose tissue
Citation: Madhu SV, Mishra BK, Mannar V, Aslam M, Banerjee B and Agrawal V (2022) TCF7L2 gene associated postprandial triglyceride dysmetabolism- a novel mechanism for diabetes risk among Asian Indians. Front. Endocrinol. 13:973718. doi: 10.3389/fendo.2022.973718
Received: 20 June 2022; Accepted: 31 August 2022;
Published: 03 October 2022.
Edited by:
Alok Raghav, Ganesh Shankar Vidyarthi Memorial Medical College, IndiaReviewed by:
Ziyi Song, Guangxi University, ChinaRuth Gutierrez-Aguilar, Universidad Nacional Autónoma de México, Mexico
Jamal Ahmad, Aligarh Muslim University, India
Copyright © 2022 Madhu, Mishra, Mannar, Aslam, Banerjee and Agrawal. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sri Venkata Madhu, ZHJzdm1hZGh1QGdtYWlsLmNvbQ==
†ORCID: S V Madhu orcid.org/0000-0003-0018-5984