Abstract
Background:
Type 2 diabetes mellitus (T2DM) has been shown to increase the risks of cognitive decline and dementia. Paired box gene 4 (PAX4), a transcription factor for beta cell development and function, has recently been implicated in pathways intersecting Alzheimer’s disease and T2DM.
Objective:
In this report, we evaluated the association of the ethnic-specific PAX4 R192H variant, a T2DM risk factor for East Asians which contributes to earlier diabetes onset, and cognitive function of Chinese T2DM patients.
Methods:
590 Chinese patients aged 45–86 from the SMART2D study were genotyped for PAX4 R192H variation using Illumina OmniExpress-24 Array. The Repeatable Battery for the Assessment of Neuropsychological Status (RBANS) which had been validated in the Singapore population was administered to assess five cognitive domains: immediate memory, visuospatial/constructional, language, attention, and delayed memory. Multiple linear regression was used to assess the association of the R192H risk allele and cognitive domains.
Results:
Patients with two PAX4 R192H risk alleles showed significantly lower attention index score (β= –8.46, 95% CI [–13.71, –3.21], p = 0.002) than patients with wild-type alleles after adjusting for age, gender, diabetes onset age, HbA1c, body-mass index, renal function, lipid profiles, systolic blood pressure, metformin usage, smoking history, education level, Geriatric Depression Scale score, and presence of APOE ɛ4 allele.
Conclusion:
Ethnic-specific R192H variation in PAX4 is associated with attention-specific cognitive impairment in Chinese with T2DM. Pending further validation studies, determining PAX4 R192H genotype may be helpful for early risk assessment of early-onset T2DM and cognitive impairment to improve diabetes care.
INTRODUCTION
Diabetes has been identified as a risk factor for cognitive dysfunction and dementia [1]. Cognitive decrements can present during pre-diabetic stage and progress throughout the disease course, at a rate that exceeds that of normal cognitive aging [2]. Three stages of cognitive deficits in type 2 diabetes mellitus (T2DM) can be broadly defined in order of severity: diabetes-associated cognitive decrements, mild cognitive impairment (MCI), and dementia [3]. A recent study identified Paired box gene 4 (PAX4) as the novel missing link between T2DM and Alzheimer’s disease (AD) [4].
PAX4 is a transcription factor that is essential for beta cell development and differentiation. In experimental studies, Pax4-null mice did not survive beyond 3 days of birth and were found to possess more alpha cells than beta and delta cells [5]. Homozygous PAX4 knockout (KO) rabbits generated using CRISPR/Cas9 system exhibited growth retardation, persistent hyperglycemia, reduced insulin-producing beta cells and increased glucagon-producing alpha cells. In addition, they also developed diabetic nephropathy, hepatopathy, myopathy, and cardiomyopathy [6]. PAX4 is also implicated in regulation of beta cell survival and proliferation in the mature pancreatic islets [7, 8]. While PAX4 has been widely linked to beta cell physiology and implicated in diabetes, there is sparse evidence to date to suggest its involvement in diabetes-associated cognitive impairment.
Clinically, PAX4 R192H (rs2233580) genetic variation was identified in several large-scale genome-wide association analysis to be associated with the risk of T2D at non-trivial odds-ratios ranging from 1.48 to 1.79 [9–12]. We also observed this PAX4 R192H genetic variant in our prior studies of early-onset T2DM patients with suspected monogenic diabetes [13] and were able to show that this variant contributes to earlier onset of T2DM subsequently in Chinese patients from two other large T2DM cohorts [14]. Functional and structural studies suggest that this variant is pathogenic and the defects were consistent with PAX4’s roles in beta cell regulation [10, 16].
Neurologically, there has been increasing evidence of PAX4’s involvement in neurological pathways. Volodin et al. demonstrated that PAX4 is critical for activating expression of genes required for myofibril breakdown after denervation [17]. PAX4 was also subsequently identified as one of the transcription factors affecting most of the downstream differentially expressed genes in the Parkinson’s disease regulatory network [18]. More recently, PAX4 was also found to regulate cytoskeletal changes through common effector proteins in both T2DM and AD pathways [4]. Given the link between PAX4 with T2DM and its potential role for mediating neurological functions, we sought to evaluate the association of the ethnic-specific T2DM risk factor, PAX4 R192H genetic variation, with cognitive impairment in Chinese patients with T2DM.
MATERIALS AND METHODS
Study population
Singapore Study of Macro-angiopathy and Micro-vascular Reactivity in Type 2 Diabetes (SMART2D) Phase 2 was a cross-sectional study on adults clinically diagnosed with T2DM between 21–93 years old who attended the Diabetes Centre in the Alexandra Health Private Limited (AHPL) regional hospital or a primary-care polyclinic in Northern Singapore from September 2014 to October 2018. Subjects were excluded if they had one of the following: active inflammation (e.g., systemic lupus erythromatosis or overt untreated cancer, use of non-steroidal anti-inflammatory medications with the same day of assessment, fasting glucose < 4.5 mM after phlebotomy, use of oral steroids equivalent to prednisolone > 7.5 mg/ day, contraindications to bio-impedance analysis (e.g., cardiac pacemaker), existing neuropsychiatric disease which could impair neurocognitive ability, and inability to give written informed consent. These study participants also did not have any cognitive diseases documented in their medical records. Out of a total of 1,547 patients recruited into this cohort during this period, 1,185 patients aged 45 years old and above underwent neurocognitive assessments, and a sub-set of 590 Chinese patients with no missing data or outliers (N = 2 in RBANS total score) in the study variables of interest were analyzed (Fig. 1). Written informed consent was obtained from all study participants. The study was approved by the National Healthcare Group Domain Specific Review Board in Singapore (DSRB 2014/00667). All experiments were performed in accordance with relevant guidelines and regulations.

Study sub-cohort assembly. Participants’ data was included in the analysis if they fulfilled the criteria of Chinese ethnicity with no missing data and no outliers in study variables of interest. SD, standard deviation.
Demographics and clinical data collection
Trained research nurses collected data from patients on demographics, educational level, medical conditions, and medications. Blood pressure (BP) was measured with a standard sphygmomanometer for the participants in the sitting position after resting for at least 10 min (2 readings were obtained). Blood and urine samples were sent to KTPH Department of Laboratory Medicine with accreditation from the College of the American Pathologists (CAP) for clinical measurements. Clinical parameters measured included: Haemoglobin A1c (HbA1c) with the Tina-quant Haemoglobin A1c Gen.3 (Roche Cobas® c501, Switzerland); serum low density lipoprotein cholesterol (LDL-C) with enzymatic colorimetric test (Roche Cobas® c501, Switzerland); serum creatinine with enzymatic colorimetric test (Roche Cobas® c501, Switzerland); and urinary albumin with immunoturbidimetric assay (Roche Cobas® c501, Switzerland). The Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation formula was used to estimate glomerular filtration rate (eGFR) [19]. A 15-item Geriatric Depression Scale (GDS) was used to assess for depressive symptoms [20]. Individuals were considered as having depressive symptoms if they scored ≥5 in GDS [21]. Apolipoprotein E (APOE) genotype (presence of APOE ɛ4 allele) and PAX4 R192H genotype were determined using Illumina OmniExpress-24 array (Illumina, California, USA). PAX4 R192H genotype comprises of: CC genotype (0 risk allele), CT genotype (1 risk allele), and TT genotype (2 risk alleles).
Cognitive function assessments
Repeatable Battery for the Assessment of Neuropsychological Status (RBANS) was used to assess cognitive function in 12 subsets (list learning, story memory, figure copy, line orientation, picture naming, semantic fluency, digit span, coding, list recall, list recognition, story recall and figure recall). Five domains of cognition were derived from these subsets and they were 1) attention, 2) language, 3) visuospatial or constructional abilities, 4) immediate memory, and 5) delayed memory. An overall score (total score) was derived from the five sub-domains [22]. RBANS has been widely used for cognitive assessments worldwide and has also been extensively validated and recognized in our Singapore population [23–25]. The index scores were standardized to a mean of 100 and standard deviation of 15. A cut-off of 1-2 SD is typically used to indicate mild cognitive impairment. In addition, the locally translated and validated version of the Mini-Mental State Examination (MMSE) was also used [26, 27] to assess the global cognitive function of our study participants. MMSE measures cognitive functioning in the domains of memory, attention, language, praxis, and visuospatial ability [28]. The summed scores range from 0–30 and a cut-off of < 24 is typically used to indicate dementia. For mild cognitive impairment, cut-offs ranging from 24 to 29 have been used [27]. Both RBANS and MMSE instruments were administered by trained research nurses.
Statistical analysis
Data was analyzed using SPSS 27.0 for Windows (SPSS Inc., Chicago, IL, USA). Continuous variables were presented as median and interquartile range (IQR). The Kruskal-Wallis test was used to compare continuous variables. Categorical variables were reported as frequencies or percentages. The Fisher-Freeman-Halton test was used to compare categorical variables. A two-tailed value of p < 0.05 was considered as statistically significant. PAX4 R192H genotype (CC (n = 418), CT (n = 162), and TT (n = 10)) fulfilled Hardy-Weinberg equilibrium with Chi-squared value of 1.62 and p value of 0.20. Multiple linear regression was used to assess the association of PAX4 R192H with RBANS sub-domain and total scores. Covariates known to affect cognitive function and scoring were adjusted for and they consist of age [29], gender [30], diabetes onset age (< 40, 41–64, ≥65) [31], body mass index (BMI) [32], HbA1c [31], systolic blood pressure (SBP) [33], high-density lipoprotein cholesterol (HDL-C) [34], triglycerides [35], eGFR < 60 and ≥60 ml/min/1.73 m2 [36], urinary albumin/creatinine ratio (uACR) < 30, 30–300 > 300 mg/g [37], metformin usage [21], smoking history (non-, ex-, current-) [38], education level (≥7 years) [39], GDS score (≥5) [40], and presence of APOE ɛ4 allele [41].
RESULTS
Clinical characteristics of study participants analyzed in this project (overall as well as stratified according to their PAX4 R192H genotype) are given in Table 1. This sub-cohort of 590 study participants comprised of 54.2% male participants. The median age for this cohort was 63.0 (IQR 57.0–69.0) years old and age of diabetes onset was 48.0 (IQR 40.0–55.0) years old with median duration of 14.5 (IQR 8.0–20.0) years. Patients with 2 risk alleles (TT genotype) were younger in index age (median age: 58.5 years old), diabetes onset age (median age: 42.5 years) and had longer duration of diabetes (median years: 18.0 years) as compared to those with CC or CT genotypes. This is consistent with our reported observation that presence of the minor risk (T) allele contributes to earlier onset of T2DM [14].
Clinical characteristics of study participants
Continuous variables were presented as median and interquartile range (IQR). The Kruskal-Wallis test was used to compare continuous variables. Categorical variables were reported as percentages. The Fisher-Freeman-Halton exact test was used to compare categorical variables. A two-tailed value of p < 0.05 was considered as statistically significant (bold). BMI, body mass index; HbA1c, glycated hemoglobin; HDL-C, high density lipoprotein cholesterol; LDL-C, low density lipoprotein cholesterol; uACR, urinary albumin-creatinine ratio; eGFR, estimated glomerular filtration rate; APOE, Apolipoprotein E.
To determine the association of PAX4 R192H genotype with cognitive impairment, multiple linear regression was performed with PAX4 R192H genotype as independent variable and RBANS total score and sub-domain scores as dependent variable while adjusting for age, gender, diabetes onset age, BMI, HbA1c, SBP, HDL-C, triglycerides, eGFR, uACR, metformin usage, smoking history, education level, GDS score, and presence of APOE ɛ4 allele. We observed that the presence of 2 risk alleles (TT genotype) was negatively associated with RBANS total score (β= –3.539, 95% CI [–6.932, –0.146], p = 0.041) as well as the attention index domain score (β= –8.457, 95% CI [–13.705, –3.209], p = 0.002) (Table 2). Comparable associations of the TT genotype with individual component test scores of the attention index domain, namely digit span (β= –9.564, 95% CI [–18.447, –0.682] p = 0.035) and coding (β= –7.350, 95% CI [–11.683, –3.018], p = 0.001) were also observed. MMSE total score (N = 579), however, did not show any significant association with CT genotype (β= 0.195, 95% CI [–0.145, 0.536], p = 0.261) or TT genotype (β= 0.071, 95% CI [–1.092, 1.234], p = 0.905) after full adjustment.
Association between PAX4 genotype and RBANS total and sub-domain index scores in fully adjusted multiple linear regression model
To further evaluate whether the association of PAX4 R192H genotype with reduced RBANS total score and attention index score is independent of diabetes onset age, we compared the results from 3 different models of multiple linear regression with different combinations of PAX4 R192H genotype and diabetes onset age (Table 3). As shown, adding diabetes onset age (Model 1) to the model containing PAX4 R192H genotype (Model 2) did not result in significant changes to the association in terms of coefficients and p-values. Diabetes onset age only (Model 3) did not result in significant association with the scores although the trend is consistently negative. Interaction terms of onset age and PAX4 R192H genotype added into the regression models did not give statistically significant results, suggesting they are independent of each other. In addition, stratified analysis was performed to determine the association of PAX4 R192H genotype with the cognitive scores in each of the onset age groups with adjustment for all relevant covariates. As shown in Table 4, significant associations were observed in the middle-aged (≤40) and elderly-aged (≥65) onset groups, with the greatest magnitude observed in the elderly-aged onset group with the lowest median diabetes duration of 5 years old. Collectively, these observations suggest that the association of PAX4 R192H genotype with total score and attention index score is independent of diabetes onset age.
Relationship between PAX4 R192H TT genotype and diabetes onset age on RBANS total score and attention index score
Multiple linear regression was used to assess the association of PAX4 R192H genotype and/or diabetes onset age with RBANS total score and attention index score, adjusting for Model 1: full set of covariates including both PAX4 R192H genotype and diabetes onset age, Model 2: full set of covariates including PAX4 R192H genotype but excluding diabetes onset age, Model 3: full set of covariates including diabetes onset age but excluding PAX4 R192H genotype. *Coefficient shown is for diabetes onset age≤40 years old against reference of ≥65 years old.
Association of PAX4 R192H TT genotype with RBANS total score and attention index score stratified by diabetes onset age
Multiple linear regression was used to assess the association of PAX4 R192H genotype and/or diabetes onset age with RBANS total score and attention index score, adjusting for full set of covariates except diabetes onset age.
DISCUSSION
In this paper, we report a novel observation that PAX4 R192H genotype is associated with reduced cognitive function in a Singaporean cohort of Chinese patients with T2DM. For patients with TT genotype, the RBANS total score and attention index score were reduced by 3.5 and 8.5 points (0.2 SD and 0.6 SD) respectively. This suggests that the PAX4 R192H genetic variation is associated with cognitive decrements specifically in areas related to executive functions for T2DM patients. As there is no existing evidence to explain the selective association of PAX4 R192H genetic variation with attention, it would be useful to corroborate this observation using brain imaging to further elucidate the specific pathophysiological pathways underlying such cognitive decrements. No significant associations between PAX4 R192H genetic variation and MMSE were observed. This is not unexpected as MMSE has been reported to be less sensitive towards mild cognitive impairment [27, 42].
As PAX4 R192H genetic variation is also known to associate with earlier onset of T2DM in Chinese, it is of interest to determine if this association with cognitive impairment was indirectly caused by the earlier T2DM onset resulting in longer disease duration and increased likelihood of sub-optimal glycemic control resulting in secondary complications such as cognitive impairment [43]. Our analyses support the independent association of PAX4 R192H genotype with cognitive function. Interestingly, our stratified analyses revealed that the association was most prominent in the elderly-onset group with the shortest diabetes duration. This calls to attention the possibility of multiple roles of PAX4 in both T2DM and also cognitive function. While PAX4 has been implicated in beta cell function, there has been no prior reports to suggest a direct role in cognitive function. Nonetheless, our observations contribute to growing evidence pointing to PAX4’s potential involvement in the neurocognitive domain [17, 18]. Given that PAX4 regulates expression of Grb2 and NOX4 which are the effector proteins shared by T2DM and AD pathways, further functional studies to decipher the relationships between R192H variation and these 2 key players would be important to establish the mechanistic link between PAX4 R192H genotype with T2DM and cognitive function. In addition, it would be of interest to explore the interaction of PAX4 R192H T allele with APOE ɛ4 allele in a larger cohort sufficiently powered to determine their contribution to the risk of cognitive impairment. Our stratified analyses have demonstrated that the association of this variant with cognitive impairment is independent of its association with earlier onset of T2DM, suggesting its potential dual roles in T2DM and cognitive impairment. Therefore, it would be interesting to determine if indeed the association of this PAX4 variant with cognitive impairment is also present in the healthy population in future studies.
Our study has several limitations. First, since this is a cross-sectional study, we are unable to establish the causal contribution of PAX4 R192H genotype to cognitive function. Secondly, we were unable to corroborate our findings with brain imaging data as advanced magnetic resonance imaging and positron emission tomography techniques [44–46] required to visualize mild cognitive impairments were beyond the accessibility of this study. Thirdly, concurrent clinical diagnosis for mild cognitive impairment was not available for our cohort. Finally, our findings are preliminary and would benefit from validation in an alternative cohort to further establish the association.
Conclusion
In summary, we have shown that PAX4 R192H TT genotype is associated with increased cognitive impairment specifically in the attention domain. This novel but preliminary finding suggests the potential use of PAX4 R192H genotype for screening and risk assessment of Chinese T2DM patients for early intervention and management of diabetes-associated cognitive dysfunction.
Footnotes
ACKNOWLEDGMENTS
We would like to thank the patients for their participation and doctors for their referrals without which this study would not have been possible.
Su Chi Lim is supported by the Singapore Ministry of Health’s National Medical Research Council under its Clinician Scientist Award (MOH-000714-01). The SMART2D cohort is supported by the Singapore Ministry of Health’s National Medical Research Council under its CS-IRG (MOH-000066).
