Abstract
Introduction: Despite the prevalence of depression and anxiety worldwide, their aetiologies remain unclear, and they can be difficult to diagnose and treat. Changes in salt-taste perception have been found in both conditions. Single-nucleotide polymorphisms (SNPs) in the salt-taste-related gene, TRPV1, have been associated with alterations to salt-taste perception, preference, and sodium consumption. Diet quality is a known modifier of depression and anxiety and recently, sodium intake has been studied in mental health. However, the relationships between salt-taste genetics, depression, anxiety, and these dietary factors are yet to be elucidated. Methods: Data from the well-characterized cross-sectional Retirement Health and Lifestyle Study (n = 536, ≥65 y) were used to explore the relationships between the salt-taste SNP TRPV1-rs8065080, levels of depression and anxiety (Hospital Anxiety and Depression Scale, HADS), estimated sodium intake, and diet quality in this secondary analysis. Standard least-squares regression and nominal logistic regression modelling were used to compare continuous and categorical variables, respectively, with analyses stratified by sex. Results: Presence of the TRPV1-rs8065080 variant allele (C) was found to increase the likelihood of having depression (HADS) in the total population and in males. The associations remained significant after adjusting for sodium intake, three diet quality indices, and demographic variables, suggesting that TRPV1-rs8065080 genotype is driving the association with depression. Discussion/Conclusion: Future studies should explore extra-oral functions of the SNP and salt-taste receptors in the brain and the roles of neurotransmitters common to both depression and salt taste to improve the management of this increasingly prevalent and difficult-to-treat condition.
Introduction
Depression and anxiety are common mental health conditions, impacting more than 280 million people worldwide [1]. Both conditions contribute to significant reductions in disability-adjusted life years [1] and cost the global economy an estimated USD 1 trillion each year [2]. Concerns over worldwide prevalence in older adults (depression 7%; anxiety 3.8%) are heightened by the impact of rapidly ageing populations [3]. Furthermore, depression severity increases, and prognosis worsens with age [4, 5]. The aetiologies of depression and anxiety are unclear, and diagnoses can be difficult [6, 7]. Treatment effectiveness varies for both conditions [8‒15] and medication side effects lead to non-adherence [14‒17]. Therefore, additional information that assists with detection, prevention, and treatment is needed.
There is evidence indicating a healthy diet that includes a variety of nutrient-rich foods can have positive effects on mental health [18‒20]. Nutrient adequacy and diet quality are inversely associated with depression [21‒27] and anxiety [25, 28, 29]. Poor quality diets are typically characterized by higher intakes of discretionary and processed foods, which are often high in sodium [30‒32]. Data suggesting a potential link between sodium intake and mental health have begun to emerge with a large cross-sectional study (n = ∼10,000) finding people with depression were more likely to add salt to their food [33]. In contrast, women with low salt intakes were more likely to be depressed [33]. In this study by Goldstein et al. [33], sodium intake was higher in men, which was hypothesized to explain their lower rates of depression. However, other studies have shown that high sodium intake is associated with increased inflammation and oxidative stress [34‒37], which are thought to play a role in the development of depression and anxiety [38, 39].
Variations in taste can influence dietary preferences and intakes of foods [40‒44] important to mental health [18, 45, 46]. Changes in the way we perceive taste have been found in mental health conditions [47‒49] such as increased salt-taste thresholds in chronic anxiety [48] and depression [49]. Taste sensitivity, dietary preferences, and food intakes are linked to variance in genes coding for taste receptors [50‒54]. The TRPV1 (transient receptor potential cation subfamily V member 1) gene codes for a cation channel involved in the detection and modulation of salt taste [51, 55, 56]. Single-nucleotide polymorphisms (SNPs) in this gene have been associated with changes to salt-taste perception [57‒61]. The most frequently implicated SNP is TRPV1-rs8065080 for which T allele carriage is associated with higher taste sensitivity to salt [57‒59], reduced liking of salty foods [59], and higher consumption of sodium [57]. Other TRPV1 SNPs have been studied with varying outcomes. For example, the rs4790522 and rs150908 SNPs have been associated with alterations in salt taste sensitivity [60, 61] while an earlier study found no associations [58]. Whether the relationships between salt-taste genetics, sodium intake, and diet quality are involved in depression and anxiety remains to be explored.
Taste receptors are also expressed extra-orally and have additional non-gustatory functions [62]. The presence of taste receptors in brain regions [63, 64] that demonstrate functional and structural abnormalities in anxiety [65, 66] and depression [67, 68] further highlights possible involvement of these receptors in the pathology of these conditions. Anxiolytic [65, 69] and antidepressant behaviours have been observed in TRPV1-knockout mice [69]. Changes in serotonin (5-HT) and gamma-amino-butyric-acid (GABA) receptors of the brain in the absence of TRPV1 receptors have also been observed [69]. Other neurotransmitters are modulated by TRPV1 channel activation such as the release of glutamate from synapses [70, 71] and the reduction of GABA, dopamine, and other catecholamines [69, 71]. Neurotransmitters released after salt stimulation of taste receptors, such as 5-HT, GABA, and norepinephrine (NE) [72, 73], are altered in anxiety and depression [74, 75] and have also been found to modulate taste thresholds [76]. Together, these studies indicate possible involvement of TRPV1 receptors in anxiety and depression; however, the role of SNPs in genes coding for salt-taste receptors is yet to be investigated.
Exploring associations between taste genetics, depression, and anxiety may identify new markers to support risk assessment and diagnosis and develop new dietary approaches useful to their management and prevention. Therefore, this cross-sectional study explored the associations between the TRPV1-rs8065080 salt-taste genotype, sodium intake, diet quality, and indicative levels of depression and anxiety in an elderly cohort. This study used the Hospital Anxiety and Depression Scale (HADS) to obtain data on these mental health states and a food frequency questionnaire (FFQ) to generate three diet quality indices and estimated sodium intake data.
Materials and Methods
Subjects
Data for this secondary cross-sectional analysis were obtained from the Retirement Health and Lifestyle Study (RHLS). Participants were eligible to take part in the RHLS if they were aged 65 years or older, were residents of the Wyong and Gosford local government areas of the NSW Central Coast of Australia, and were either living independently in private dwellings or retirement villages [77]. Participants were excluded if they were unable to provide written informed consent. They were not excluded based on the presence of pre-existing health conditions, including depression or anxiety [77]. This sub-study (n = 536) included those who had provided blood samples, were successfully genotyped for TRPV1-rs865080, and completed a valid FFQ. Ethics approval for the RHLS research was granted by the Human Research Ethics Committee of the University of Newcastle (Reference No. H-2008-0431) [78, 79].
Demographics and Anthropometrics
Prior studies have found age, sex, and diet quality are associated with taste genetics and/or mental health [53, 80‒84]. Smoking is known to modify genetic expression [85] and has been linked to the presence of bitter TAS2R38 haplotypes [86]. Depression and anxiety have been correlated to income, education [87], and body mass index (BMI) [88]. Therefore, these variables were included as possible confounders to the associations between TRPV1-rs8065080 and HADS. Data on age, sex, income, education, and smoking history were obtained through interviewer-administered questionnaires. The measurement of height and weight followed the International Society for the Advancement of Kinanthropometry (ISAK) guidelines [89]. Digital scales (Wedderburn© UWPM150 Platform) were used to measure weight, which was recorded to the nearest 0.01 kg [77, 89]. The stretch stature method was used to measure height, which was recorded to the nearest 0.01 cm [77, 89]. BMI (BMI = weight [kg]/height [m2]) calculations were derived from the height and weight measures.
Genotyping
Fasted whole blood samples were collected from participants into EDTA-lined tubes, then stored at −20°C. The QIAGEN QIAmp DNA mini kit was used to isolate DNA from peripheral blood cells following the manufacturer’s instructions. TaqMan™ assay code C___11679656_10 (Applied Biosystems™, ThermoFisher Scientific, CA, USA) was used for allelic discrimination of TRPV1-rs8065080. Further detail on the genotyping protocol can be found in a prior publication [90].
Assessment of Depression and Anxiety
Indicative depression and anxiety scores were obtained for participants through the self-administered HADS assessment tool. The HADS has been widely used in hospital (non-psychiatric) and general population research [91] and is recommended as a reliable measure of anxiety and depression [91, 92]. The HADS is composed of 14 questions equally divided to produce scores measuring each subscale of depression (HADS-D) or anxiety (HADS-A) [92]. The subscales are each scored between 0 and 21 with a score of ≥8 within a subscale indicating the presence of either depression or anxiety [92].
Dietary Assessment
Three diet quality indices and daily sodium intake estimates were obtained via data from an FFQ containing 225 items [93]. Each of the indices focuses on different aspects of diet quality. The Dietary Guideline Index (DGI) provides an assessment of key nutrient intakes and the proportions obtained from healthy food types, food diversity, and unhealthy foods [94]. The Australian Recommended Food Score (ARFS) assesses dietary variety as the indicator of diet quality [95], and the Australian Healthy Eating Index (Aust-HEI) includes measures of dietary variety, healthy food choices, and consumption of fruits, vegetables, fats, and discretionary foods [96]. Inverse associations have been found in previous studies between depression and DGI [27], ARFS [22, 24], and Aust-HEI [26] diet quality scores and between anxiety and ARFS [22, 97] and Aust-HEI [26] diet quality scores. Further information about each of the indices and the scoring methods can be found in a prior publication [90].
Statistical Analyses
The data analyses were undertaken using JMP (Pro v.14.2.0; SAS Institute Inc., Cary, NC, USA 27513). The cohort characteristics are described for continuous variable distributions (means, 95% confidence intervals, and standard deviations) and categorical variable distributions (number and percentage of cohort). The TRPV1-rs8065080 polymorphism was analysed by presence or absence of the variant allele (C) and genotype occurrence. Both genotype and HADS data were reported as the number and percentage of the study cohort. Due to sex differences in salt taste [52, 59], in TRPV1 findings [58], and in anxiety and depression [1, 98], results were further stratified by sex.
Continuous variables were compared using standard least-squares regression analyses (least-squares mean, 95% confidence intervals, p values) and categorical variables through nominal logistic regression analyses (odds ratios, χ2, p values). p values are presented to one significant number and a threshold of <0.05 was considered statistically significant. Where appropriate, results were adjusted for potential confounding factors such as age, sex, education, income, smoking status, and BMI. In these adjustment models, values missing at random conditional on variables were omitted from the analyses. As there are relationships between diet quality, sodium intake, and depression [21‒27, 33] and anxiety [26, 28, 29], further adjustments for these dietary measures were made.
Results
Participant Characteristics
The average age of the 536 participants was 77.4 (SD ± 6.8) years, the mean BMI was 28.5 kg/m2 (SD ± 4.8, Table 1), and there were no differences by sex (55% female, Table 2) for either variable (Table 1). The mean diet-quality scores were 97.0/150 points (DGI), 29.2/74 points (ARFS), and 30.5/60 points (Aust-HEI) (Table 1). Females had higher mean scores than males across all indices (DGI: p = 0.0004, ARFS: p = 0.002, and Aust-HEI: p < 0.0001, Table 1). Estimated mean sodium intake in the total population was 2,053.0 (SD ± 842.8) mg/day (Table 1) and was higher in males than females (2,226.9 mg/day vs. 1,910.8 mg/day, respectively, p < 0.0001, Table 1).
Distributions of continuous variables
Variable . | Total . | Females . | Males . | p value . | ||||||
---|---|---|---|---|---|---|---|---|---|---|
mean (SD) . | min . | max . | mean (SD) . | min . | max . | mean (SD) . | min . | max . | ||
Age, years | 77.4 (±6.8) | 65.0 | 94.0 | 77.5 (±6.9) | 65 | 94 | 77.3 (±6.7) | 65 | 93 | 0.8 |
BMI, kg/m2 | 28.5 (±4.8) | 17.1 | 46.3 | 28.5 (±5.1) | 17.6 | 46.3 | 28.5 (±4.5) | 17.1 | 45.4 | 0.9 |
DGI | 97.0 (±15.8) | 30.9 | 132.6 | 99.2 (±16.1) | 30.9 | 132.6 | 94.3 (±15.0) | 51.8 | 130.4 | 0.0004 |
ARFS | 29.2 (±8.1) | 6.0 | 50.0 | 30.2 (±8.2) | 6 | 50 | 28.0 (±7.9) | 10 | 49 | 0.002 |
Aust-HEI | 30.5 (±9.5) | 4.9 | 50.8 | 32.1 (±9.0) | 6.4 | 50.8 | 28.5 (±9.7) | 4.9 | 46.5 | <0.0001 |
Sodium, mg/day | 2,053.0 (±842.8) | 505.7 | 8,250.2 | 1,910.8 (±840.1) | 505.7 | 8,250.2 | 2,226.9 (±814.6) | 535.2 | 5,481.8.2 | <0.0001 |
Variable . | Total . | Females . | Males . | p value . | ||||||
---|---|---|---|---|---|---|---|---|---|---|
mean (SD) . | min . | max . | mean (SD) . | min . | max . | mean (SD) . | min . | max . | ||
Age, years | 77.4 (±6.8) | 65.0 | 94.0 | 77.5 (±6.9) | 65 | 94 | 77.3 (±6.7) | 65 | 93 | 0.8 |
BMI, kg/m2 | 28.5 (±4.8) | 17.1 | 46.3 | 28.5 (±5.1) | 17.6 | 46.3 | 28.5 (±4.5) | 17.1 | 45.4 | 0.9 |
DGI | 97.0 (±15.8) | 30.9 | 132.6 | 99.2 (±16.1) | 30.9 | 132.6 | 94.3 (±15.0) | 51.8 | 130.4 | 0.0004 |
ARFS | 29.2 (±8.1) | 6.0 | 50.0 | 30.2 (±8.2) | 6 | 50 | 28.0 (±7.9) | 10 | 49 | 0.002 |
Aust-HEI | 30.5 (±9.5) | 4.9 | 50.8 | 32.1 (±9.0) | 6.4 | 50.8 | 28.5 (±9.7) | 4.9 | 46.5 | <0.0001 |
Sodium, mg/day | 2,053.0 (±842.8) | 505.7 | 8,250.2 | 1,910.8 (±840.1) | 505.7 | 8,250.2 | 2,226.9 (±814.6) | 535.2 | 5,481.8.2 | <0.0001 |
SD, standard deviation; BMI, body mass index; DGI, Dietary Guideline Index (150 points); ARFS, Australian Recommended Food Score (74 points); Aust-HEI, Australian Health Eating Index (60 points).
Distributions of categorical variables
Variable . | Total, n (%) . | Females, n (%) . | Males, n (%) . | p value . |
---|---|---|---|---|
Sex | ||||
Males | 241 (45.0) | |||
Females | 295 (55.0) | |||
Income, AUD/year | <0.0001 | |||
<10,000 | 2 (0.4) | 2 (0.7) | 0 (0.0) | |
10,000–20,000 | 163 (30.5) | 128 (43.6) | 35 (14.5) | |
20,000–40,000 | 263 (49.2) | 118 (40.1) | 145 (60.2) | |
40,000–60,000 | 59 (11.0) | 22 (7.5) | 37 (15.4) | |
60,000–80,000 | 22 (4.1) | 8 (2.7) | 14 (5.8) | |
80,000–100,000 | 10 (1.9) | 5 (1.7) | 5 (2.1) | |
≥100,000 | 4 (0.8) | 2 (0.7) | 2 (1.9) | |
Education | <0.0001 | |||
No formal schooling | 0 (0.0) | 0 (0.0) | 0 (0.0) | |
School | 177 (33.1) | 119 (40.5) | 58 (24.1) | |
Trade/certificate/apprenticeship | 102 (19.1) | 25 (8.5) | 77 (32.0) | |
Other certificate (TAFE or business college) | 106 (19.8) | 76 (25.9) | 30 (12.5) | |
Dip. or adv. dip. (TAFE or business college) | 87 (16.3) | 46 (15.7) | 41 (17.0) | |
Bachelor degree (incl. honours) | 31 (5.8) | 14 (4.8) | 17 (7.1) | |
Graduate dip. or graduate certificate | 16 (3.0) | 9 (3.1) | 7 (2.9) | |
Postgraduate doctoral or master’s degree | 16 (3.0) | 5 (1.7) | 11 (4.6) | |
Smoking | <0.0001 | |||
Daily | 14 (2.6) | 7 (2.4) | 7 (2.9) | |
Weekly but not daily | 1 (0.2) | 1 (0.03) | 0 (0.0) | |
Less often than weekly | 0 (0.0) | 0 (0.0) | 0 (0.0) | |
Do not smoke now but used to tried | 227 (42.4) | 82 (27.8) | 145 (60.2) | |
Smoking but never smoked regularly | 26 (4.9) | 16 (5.4) | 10 (4.1) | |
Never smoked | 268 (50.1) | 189 (64.1) | 79 (32.8) |
Variable . | Total, n (%) . | Females, n (%) . | Males, n (%) . | p value . |
---|---|---|---|---|
Sex | ||||
Males | 241 (45.0) | |||
Females | 295 (55.0) | |||
Income, AUD/year | <0.0001 | |||
<10,000 | 2 (0.4) | 2 (0.7) | 0 (0.0) | |
10,000–20,000 | 163 (30.5) | 128 (43.6) | 35 (14.5) | |
20,000–40,000 | 263 (49.2) | 118 (40.1) | 145 (60.2) | |
40,000–60,000 | 59 (11.0) | 22 (7.5) | 37 (15.4) | |
60,000–80,000 | 22 (4.1) | 8 (2.7) | 14 (5.8) | |
80,000–100,000 | 10 (1.9) | 5 (1.7) | 5 (2.1) | |
≥100,000 | 4 (0.8) | 2 (0.7) | 2 (1.9) | |
Education | <0.0001 | |||
No formal schooling | 0 (0.0) | 0 (0.0) | 0 (0.0) | |
School | 177 (33.1) | 119 (40.5) | 58 (24.1) | |
Trade/certificate/apprenticeship | 102 (19.1) | 25 (8.5) | 77 (32.0) | |
Other certificate (TAFE or business college) | 106 (19.8) | 76 (25.9) | 30 (12.5) | |
Dip. or adv. dip. (TAFE or business college) | 87 (16.3) | 46 (15.7) | 41 (17.0) | |
Bachelor degree (incl. honours) | 31 (5.8) | 14 (4.8) | 17 (7.1) | |
Graduate dip. or graduate certificate | 16 (3.0) | 9 (3.1) | 7 (2.9) | |
Postgraduate doctoral or master’s degree | 16 (3.0) | 5 (1.7) | 11 (4.6) | |
Smoking | <0.0001 | |||
Daily | 14 (2.6) | 7 (2.4) | 7 (2.9) | |
Weekly but not daily | 1 (0.2) | 1 (0.03) | 0 (0.0) | |
Less often than weekly | 0 (0.0) | 0 (0.0) | 0 (0.0) | |
Do not smoke now but used to tried | 227 (42.4) | 82 (27.8) | 145 (60.2) | |
Smoking but never smoked regularly | 26 (4.9) | 16 (5.4) | 10 (4.1) | |
Never smoked | 268 (50.1) | 189 (64.1) | 79 (32.8) |
AUD, Australian dollars; dip., diploma; adv. dip., advanced diploma; TAFE, Technical and Further Education.
The household income of the largest proportion of participants was between AUD 20,000–AUD 40,000/year (49.2%) with males reporting higher incomes than females (p <0.0001, Table 2). Males were more likely to be educated at trade level or higher (75.9% vs. 59.6%, p < 0.0001) and to have a history of smoking (67.2% vs. 36.9%, p < 0.0001) (Table 2). Where groups had insufficient numbers for statistical analyses as adjustment variables, they were consolidated. Income was collapsed to three categories (“<AUD 20,000”; “AUD 20,000–40,000”; “>AUD 40,000”), education to three categories (“≤school”; “trade/TAFE/certificate”; “≥diploma/bachelor”), and smoking to two categories (“history of smoking”; “never smoked”).
Genotype Distributions
The TRPV1-rs8065080 variant allele (C) occurred at a frequency of 0.36 and the ancestral allele (T) at 0.64. 41.0% of participants were homozygous for TT, 46.7% were heterozygous for CT, and 12.3% were homozygous for CC (Table 3). The genotype frequencies in the sample (CC 12.3%, CT 46.7%, TT 41.0%) did not deviate from expected Hardy-Weinberg equilibrium frequencies (CC 16.8%, CT 46.7%, TT 34.7%) (p = 0.16). The TRPV1-C allele (CC or CT genotypes) was carried by 58.9% of participants, with no difference between female and male subjects (Table 3). Due to low numbers in the CC homozygous group, analyses were undertaken by presence versus absence of the variant allele for more robust statistical power.
TRPV1-rs8065080 variant (C) allele and genotype distributions
. | Total, n (%) . | Females, n (%) . | Males, n (%) . | P value . |
---|---|---|---|---|
TRPV1-C allele present | 316 (58.9) | 167 (56.6) | 149 (61.8) | 0.2 |
TRPV1-C allele absent | 220 (41.1) | 128 (43.4) | 92 (38.2) | |
CC | 66 (12.3) | 33 (11.2) | 33 (137 | 0.4 |
CT | 250 (46.7) | 134 (45.4) | 116 (48.1) | |
TT | 220 (41.0) | 128 (43.3) | 92 (38.2) |
. | Total, n (%) . | Females, n (%) . | Males, n (%) . | P value . |
---|---|---|---|---|
TRPV1-C allele present | 316 (58.9) | 167 (56.6) | 149 (61.8) | 0.2 |
TRPV1-C allele absent | 220 (41.1) | 128 (43.4) | 92 (38.2) | |
CC | 66 (12.3) | 33 (11.2) | 33 (137 | 0.4 |
CT | 250 (46.7) | 134 (45.4) | 116 (48.1) | |
TT | 220 (41.0) | 128 (43.3) | 92 (38.2) |
HADS Distributions
HADS-D scores indicative of depression occurred in 8% of total participants, 8.1% of females, and 8.3% of males (Table 4). HADS-A scores indicative of anxiety occurred in 13.4% of total participants, 15.9% of females, and 10.4% of males (Table 4). There were no statistically significant differences between the sexes for either HADS-D or HADS-A scores.
Threshold frequencies for HADS-D and HADS-A
. | Total . | Female . | Male . | p value . | |||
---|---|---|---|---|---|---|---|
≥8, n (%) . | <8, n (%) . | ≥8, n (%) . | <8, n (%) . | ≥8, n (%) . | <8, n (%) . | ||
HADS-D | 44 (8.0) | 492 (91.8) | 24 (8.1) | 271 (91.9) | 20 (8.3) | 221 (91.7) | 0.9 |
HADS-A | 72 (13.4) | 464 (86.6) | 47 (15.9) | 248 (84.1) | 25 (10.4) | 216 (89.6) | 0.06 |
. | Total . | Female . | Male . | p value . | |||
---|---|---|---|---|---|---|---|
≥8, n (%) . | <8, n (%) . | ≥8, n (%) . | <8, n (%) . | ≥8, n (%) . | <8, n (%) . | ||
HADS-D | 44 (8.0) | 492 (91.8) | 24 (8.1) | 271 (91.9) | 20 (8.3) | 221 (91.7) | 0.9 |
HADS-A | 72 (13.4) | 464 (86.6) | 47 (15.9) | 248 (84.1) | 25 (10.4) | 216 (89.6) | 0.06 |
HADS-D threshold ≥8, depression; HADS-A threshold ≥8, anxiety.
TRPV1-rs8065080, HADS, and Dietary Measures
There were no statistically significant associations between TRPV1-C allele presence and any of the dietary measures in the total population, in females, or males (Table 5). Participants with HADS-D scores indicative of depression had lower diet quality scores on both the DGI (89.7.3 vs. 97.7, p = 0.001) and ARFS (26.3 vs. 29.5, p = 0.02) than those with normal scores, but there was no association with the Aust-HEI index (Table 6). These significant associations were also present in females (DGI 90.0 vs. 100.0, p = 0.003; ARFS 26.8 vs. 30.5, p = 0.04) but not in males (Table 6). In the total population and by sex, there were no associations between HADS-A scores and the diet quality indices (Table 7) or between HADS-D or HADS-A scores and sodium intake (Tables 6, 7).
Association between TRPV1-rs8065080 variant (C) allele presence and dietary measures
. | Total . | p value . | Female . | p value . | Male . | p value . | |||
---|---|---|---|---|---|---|---|---|---|
C allele present . | C allele absent . | C allele present . | C allele absent . | C allele present . | C allele absent . | ||||
LSM (95% CI) . | LSM (95% CI) . | LSM (95% CI) . | |||||||
Sodium, mg/day | 2,073.9 (1,980.7–2,167.1) | 2,022.9 (1,911.2–2,134.6) | 0.5 | 1,958.8 (1,830.9–2,086.7) | 1,848.2 (1,702.2–1,994.3) | 0.3 | 2,202.9 (2,071.2–2,334.5) | 2,265.9 (2,098.4–2,334.5) | 0.6 |
DGI | 96.3 (94.5–98.0) | 98.0 (95.9–100.1) | 0.2 | 98.8 (96.4–101.3) | 99.6 (96.8–102.4) | 0.7 | 93.4 (90.9–95.8) | 95.8 (92.7–98.9) | 0.2 |
ARFS | 29.5 (28.6–30.4) | 28.7 (27.6–29.8) | 0.3 | 30.8 (29.6–32.0) | 29.3 (27.9–30.8) | 0.1 | 28.1 (26.9–29.4) | 27.9 (26.2–29.5) | 0.8 |
Aust-HEI | 30.2 (29.2–31.3) | 30.8 (29.5–32.0) | 0.2 | 32.3 (30.9–33.7) | 31.8 (30.2–33.3) | 0.6 | 27.9 (26.3–29.4) | 29.4 (27.5–31.4) | 0.2 |
. | Total . | p value . | Female . | p value . | Male . | p value . | |||
---|---|---|---|---|---|---|---|---|---|
C allele present . | C allele absent . | C allele present . | C allele absent . | C allele present . | C allele absent . | ||||
LSM (95% CI) . | LSM (95% CI) . | LSM (95% CI) . | |||||||
Sodium, mg/day | 2,073.9 (1,980.7–2,167.1) | 2,022.9 (1,911.2–2,134.6) | 0.5 | 1,958.8 (1,830.9–2,086.7) | 1,848.2 (1,702.2–1,994.3) | 0.3 | 2,202.9 (2,071.2–2,334.5) | 2,265.9 (2,098.4–2,334.5) | 0.6 |
DGI | 96.3 (94.5–98.0) | 98.0 (95.9–100.1) | 0.2 | 98.8 (96.4–101.3) | 99.6 (96.8–102.4) | 0.7 | 93.4 (90.9–95.8) | 95.8 (92.7–98.9) | 0.2 |
ARFS | 29.5 (28.6–30.4) | 28.7 (27.6–29.8) | 0.3 | 30.8 (29.6–32.0) | 29.3 (27.9–30.8) | 0.1 | 28.1 (26.9–29.4) | 27.9 (26.2–29.5) | 0.8 |
Aust-HEI | 30.2 (29.2–31.3) | 30.8 (29.5–32.0) | 0.2 | 32.3 (30.9–33.7) | 31.8 (30.2–33.3) | 0.6 | 27.9 (26.3–29.4) | 29.4 (27.5–31.4) | 0.2 |
LSM, least-squares mean; CI, confidence interval; DGI, Dietary Guideline Index; ARFS, Australian Recommended Food Score; Aust-HEI, Australian Health Eating Index.
Association between depression (HADS-D) and dietary measures
. | Total . | p value . | Female . | p value . | Male . | p value . | |||
---|---|---|---|---|---|---|---|---|---|
≥8 . | <8 . | ≥8 . | <8 . | ≥8 . | <8 . | ||||
LSM (95% CI) . | LSM (95% CI) . | LSM (95% CI) . | |||||||
Sodium, mg/day | 2,108.0 (1,858.2–2,357.7) | 2,048.0 (1,973.3–2,122.7) | 0.7 | 1,972.6 (1,634.6–2,310.6) | 1,905.4 (1,804.8–2,006.0) | 0.7 | 2,270.4 (1,910.9–2,630.0) | 2,223.0 (2,114.8–2,331.1) | 0.8 |
DGI | 89.7 (85.0–94.3) | 97.7 (96.3–99.0) | 0.001 | 90.0 (83.6–96.3) | 100.0 (98.1–101.9) | 0.003 | 89.3 (82.7–95.9) | 94.8 (92.8–96.8) | 0.1 |
ARFS | 26.3 (24.0–28.7) | 29.5 (28.7–30.2) | 0.02 | 26.8 (23.5–30.1) | 30.5 (29.5–31.4) | 0.04 | 25.8 (22.3–29.3) | 28.2 (27.2–29.3) | 0.2 |
Aust-HEI | 28.5 (25.7–21.3) | 30.6 (29.8–31.5) | 0.2 | 29.4 (25.8–33.0) | 32.2 (31.2–33.4) | 0.1 | 27.6 (23.3–31.9) | 28.6 (27.3–29.9) | 0.7 |
. | Total . | p value . | Female . | p value . | Male . | p value . | |||
---|---|---|---|---|---|---|---|---|---|
≥8 . | <8 . | ≥8 . | <8 . | ≥8 . | <8 . | ||||
LSM (95% CI) . | LSM (95% CI) . | LSM (95% CI) . | |||||||
Sodium, mg/day | 2,108.0 (1,858.2–2,357.7) | 2,048.0 (1,973.3–2,122.7) | 0.7 | 1,972.6 (1,634.6–2,310.6) | 1,905.4 (1,804.8–2,006.0) | 0.7 | 2,270.4 (1,910.9–2,630.0) | 2,223.0 (2,114.8–2,331.1) | 0.8 |
DGI | 89.7 (85.0–94.3) | 97.7 (96.3–99.0) | 0.001 | 90.0 (83.6–96.3) | 100.0 (98.1–101.9) | 0.003 | 89.3 (82.7–95.9) | 94.8 (92.8–96.8) | 0.1 |
ARFS | 26.3 (24.0–28.7) | 29.5 (28.7–30.2) | 0.02 | 26.8 (23.5–30.1) | 30.5 (29.5–31.4) | 0.04 | 25.8 (22.3–29.3) | 28.2 (27.2–29.3) | 0.2 |
Aust-HEI | 28.5 (25.7–21.3) | 30.6 (29.8–31.5) | 0.2 | 29.4 (25.8–33.0) | 32.2 (31.2–33.4) | 0.1 | 27.6 (23.3–31.9) | 28.6 (27.3–29.9) | 0.7 |
LSM, least-squares mean; CI, confidence interval; HADS-D threshold score ≥8, depression; DGI, Dietary Guideline Index; ARFS, Australian Recommended Food Score; Aust-HEI, Australian Health Eating Index.
Association between anxiety (HADS-A) and dietary measures
. | Total . | p value . | Female . | p value . | Male . | p value . | |||
---|---|---|---|---|---|---|---|---|---|
≥8 . | <8 . | ≥8 . | <8 . | ≥8 . | <8 . | ||||
LSM (95% CI) . | LSM (95% CI) . | LSM (95% CI) . | |||||||
Sodium, mg/day | 2,049.3 (1,854.0–2,244.6) | 2,053.5 (1,976.6–2,130.5) | 1.0 | 1,990.7 (1,749.3–2,232.0) | 1,895.7 (1,790.6–2,000.8) | 0.4 | 2,159.5 (1,838.0–2,480.9) | 2,234.7 (2,125.4–2,344.1) | 0.7 |
DGI | 94.0 (90.4–97.7) | 97.5 (96.0–98.9) | 0.09 | 96.4 (91.8–97.7) | 99.7 (97.7–101.7) | 0.2 | 89.6 (83.7–95.5) | 94.9 (92.9–96.9) | 0.1 |
ARFS | 28.9 (27.0–30.7) | 29.3 (28.5–30.0) | 0.7 | 30.0 (27.6–32.3) | 30.2 (29.2–31.2) | 0.9 | 26.7 (23.6–29.8) | 28.2 (27.1–29.2) | 0.4 |
Aust–HEI | 30.1 (27.9–32.3) | 30.5 (29.6–31.4) | 0.7 | 31.7 (29.1–34.3) | 32.1 (31.0–33.3) | 0.8 | 27.1 (23.3–30.9) | 28.6 (27.3–29.9) | 0.5 |
. | Total . | p value . | Female . | p value . | Male . | p value . | |||
---|---|---|---|---|---|---|---|---|---|
≥8 . | <8 . | ≥8 . | <8 . | ≥8 . | <8 . | ||||
LSM (95% CI) . | LSM (95% CI) . | LSM (95% CI) . | |||||||
Sodium, mg/day | 2,049.3 (1,854.0–2,244.6) | 2,053.5 (1,976.6–2,130.5) | 1.0 | 1,990.7 (1,749.3–2,232.0) | 1,895.7 (1,790.6–2,000.8) | 0.4 | 2,159.5 (1,838.0–2,480.9) | 2,234.7 (2,125.4–2,344.1) | 0.7 |
DGI | 94.0 (90.4–97.7) | 97.5 (96.0–98.9) | 0.09 | 96.4 (91.8–97.7) | 99.7 (97.7–101.7) | 0.2 | 89.6 (83.7–95.5) | 94.9 (92.9–96.9) | 0.1 |
ARFS | 28.9 (27.0–30.7) | 29.3 (28.5–30.0) | 0.7 | 30.0 (27.6–32.3) | 30.2 (29.2–31.2) | 0.9 | 26.7 (23.6–29.8) | 28.2 (27.1–29.2) | 0.4 |
Aust–HEI | 30.1 (27.9–32.3) | 30.5 (29.6–31.4) | 0.7 | 31.7 (29.1–34.3) | 32.1 (31.0–33.3) | 0.8 | 27.1 (23.3–30.9) | 28.6 (27.3–29.9) | 0.5 |
LSM, least-squares mean; CI, confidence interval; HADS-A threshold ≥8, anxiety; DGI, Dietary Guideline Index; ARFS, Australian Recommended Food Score; Aust-HEI, Australian Health Eating Index.
TRPV1-rs8065080 and Depression (HADS-D)
Associations between TRPV1-rs8065080 and Depression (HADS-D)
Participants carrying the TRPV1-C allele were more likely to have HADS-D scores indicative of depression in the unadjusted model (p = 0.003) and in the models adjusting for age and sex (p = 0.003) and age, sex, income, education, smoking, and BMI (p = 0.01) than those not carrying the TRPV1-C allele (Table 8). With lower certainty, males carrying the TRPV1-C allele were more likely to have HADS-D scores indicative of depression in the unadjusted model (p = 0.003) and in the models adjusting for age (p = 0.003) and age, income, education, smoking, and BMI (p = 0.002) (Table 8). No associations between presence of the TRPV1-C allele and depression (HADS-D) scores were found in females (Table 8).
Association between TRPV1-rs8065080 variant (C) allele presence and depression (HADS-D)
HADS-D ≥8 . | Unadjusted . | Model 1 . | Model 2 . | |||
---|---|---|---|---|---|---|
χ2 (p value) . | OR (95%CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | |
TRPV1-C allele present | ||||||
Total cohort | 9.1 (0.003) | 2.9 (1.4–6.2) | 8.7 (0.003) | 2.9 (1.4–6.1) | 6.2 (0.01) | 2.8 (1.2–6.5) |
Females | 2.2 (0.1) | 2.0 (0.8–4.9) | 1.9 (0.2) | 1.9 (0.7–4.6) | 0.4 (0.5) | 1.4 (0.5–4.0) |
Males | 8.8 (0.003) | 6.2 (1.4–27.3) | 9.0 (0.003) | 6.3 (1.4–28.1) | 9.8 (0.002) | 11.1 (1.4–86.7) |
HADS-D ≥8 . | Unadjusted . | Model 1 . | Model 2 . | |||
---|---|---|---|---|---|---|
χ2 (p value) . | OR (95%CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | |
TRPV1-C allele present | ||||||
Total cohort | 9.1 (0.003) | 2.9 (1.4–6.2) | 8.7 (0.003) | 2.9 (1.4–6.1) | 6.2 (0.01) | 2.8 (1.2–6.5) |
Females | 2.2 (0.1) | 2.0 (0.8–4.9) | 1.9 (0.2) | 1.9 (0.7–4.6) | 0.4 (0.5) | 1.4 (0.5–4.0) |
Males | 8.8 (0.003) | 6.2 (1.4–27.3) | 9.0 (0.003) | 6.3 (1.4–28.1) | 9.8 (0.002) | 11.1 (1.4–86.7) |
HADS-D, Hospital Anxiety and Depression Scale-Depression; OR, odds ratio; CI, confidence interval; model 1, adjusted for age and sex (total cohort) or adjusted for age (by sex); model 2, adjusted for age, sex, income, education, smoking, and BMI (total cohort, n = 475) or adjusted for age, income, education, smoking, and BMI (by sex; females, n = 278; males, n = 197).
Associations between TRPV1-rs8065080 and Depression (HADS-D), Adjusting for Dietary Measures
The association between presence of the TRPV1-C allele and HADS-D scores indicative of depression remained statistically significant after adjusting for the dietary measures in the total population (DGI, p = 0.004; ARFS, p = 0.001; Aust-HEI, p = 0.003; sodium, p = 0.003) and in males (DGI, p = 0.004; ARFS, p = 0.003; Aust-HEI, p = 0.003; sodium, p = 0.003) but not in females (Table 9).
Association between TRPV1-rs8065080 variant (C) allele presence and depression (HADS-D), adjusting for dietary measures
HADS-D ≥8 . | Model 1 . | Model 2 . | Model 3 . | Model 4 . | ||||
---|---|---|---|---|---|---|---|---|
χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | |
TRPV1-C allele present | ||||||||
Total | 8.3 (0.004) | 2.8 (1.3–6.0) | 10.2 (0.001) | 3.1 (1.5–6.7) | 8.9 (0.003) | 2.9 (1.4–6.2) | 9.1 (0.003) | 2.9 (1.4–6.2) |
Females | 2.0 (0.2) | 1.9 (0.8–4.8) | 3.0 (0.08) | 2.2 (0.9–5.6) | 2.5 (0.1) | 2.0 (0.8–5.1) | 2.2 (0.1) | 1.9 (0.8–4.8) |
Males | 8.3 (0.004) | 6.0 (1.4–26.6) | 9.0 (0.003) | 6.4 (1.4–28.4) | 8.6 (0.003) | 6.1 (1.4–27.1) | 8.8 (0.003) | 6.2 (1.4–27.6) |
HADS-D ≥8 . | Model 1 . | Model 2 . | Model 3 . | Model 4 . | ||||
---|---|---|---|---|---|---|---|---|
χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | |
TRPV1-C allele present | ||||||||
Total | 8.3 (0.004) | 2.8 (1.3–6.0) | 10.2 (0.001) | 3.1 (1.5–6.7) | 8.9 (0.003) | 2.9 (1.4–6.2) | 9.1 (0.003) | 2.9 (1.4–6.2) |
Females | 2.0 (0.2) | 1.9 (0.8–4.8) | 3.0 (0.08) | 2.2 (0.9–5.6) | 2.5 (0.1) | 2.0 (0.8–5.1) | 2.2 (0.1) | 1.9 (0.8–4.8) |
Males | 8.3 (0.004) | 6.0 (1.4–26.6) | 9.0 (0.003) | 6.4 (1.4–28.4) | 8.6 (0.003) | 6.1 (1.4–27.1) | 8.8 (0.003) | 6.2 (1.4–27.6) |
HADS-D, Hospital Anxiety and Depression Scale-Depression; OR, odds ratio; CI, confidence interval; model 1: DGI, Dietary Guideline Index; model 2: ARFS, Australian Recommended Food Score; model 3: Aust-HEI, Australian Health Eating Index; model 4: sodium intake (mg/day).
TRPV1-rs8065080 and Anxiety (HADS-A)
Associations between TRPV1-rs8065080 and Anxiety (HADS-A)
There were no statistically significant associations between presence of the TRPV1-C allele and HADS-A scores indicative of anxiety in the total cohort or by sex, in unadjusted and adjusted models (Table 10). Therefore, no further analyses were undertaken.
Association between TRPV1-rs8065080 variant (C) allele presence and anxiety (HADS-A)
HADS-A ≥8 . | Unadjusted . | Model 1 . | Model 2 . | |||
---|---|---|---|---|---|---|
χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | |
TRPV1-C allele present | ||||||
Total cohort | 2.1 (0.2) | 1.5 (0.9–2.5) | 2.3 (0.1) | 1.5 (0.9–2.5) | 1.4 (0.2) | 1.4 (0.8–2.5) |
Females | 2.0 (0.2) | 1.6 (0.8–3.1) | 1.8 (0.2) | 1.6 (0.8–3.0) | 1.7 (0.2) | 1.6 (0.8–3.2) |
Males | 0.5 (0.5) | 1.4 (0.6–3.3) | 0.5 (0.5) | 1.4 (0.6–3.3) | 0.3 (0.6) | 1.3 (0.5–3.5) |
HADS-A ≥8 . | Unadjusted . | Model 1 . | Model 2 . | |||
---|---|---|---|---|---|---|
χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | χ2 (p value) . | OR (95% CI) . | |
TRPV1-C allele present | ||||||
Total cohort | 2.1 (0.2) | 1.5 (0.9–2.5) | 2.3 (0.1) | 1.5 (0.9–2.5) | 1.4 (0.2) | 1.4 (0.8–2.5) |
Females | 2.0 (0.2) | 1.6 (0.8–3.1) | 1.8 (0.2) | 1.6 (0.8–3.0) | 1.7 (0.2) | 1.6 (0.8–3.2) |
Males | 0.5 (0.5) | 1.4 (0.6–3.3) | 0.5 (0.5) | 1.4 (0.6–3.3) | 0.3 (0.6) | 1.3 (0.5–3.5) |
HADS-A, Hospital Anxiety and Depression Scale-Anxiety; OR, odds ratio; CI, confidence interval; model 1: adjusted for age and sex (total cohort) or adjusted for age (by sex); model 2: adjusted for age, sex, income, education, smoking, and BMI (total cohort; n = 475) or adjusted for age, income, education, smoking, and BMI (by sex; females, n = 278; males, n = 197).
Discussion
The key finding in this cross-sectional analysis is that presence of the TRPV1-rs8065080 variant allele (C) was found to increase the likelihood of having a HADS-D score indicative of depression in an elderly cohort. The association remained significant after adjusting for the demographic variables of age, sex, income, education, BMI, and smoking and after adjusting for the three diet quality indices and sodium intake. This suggests that the association may be related to the function of the salt-taste receptor, rather than the receptor being linked to depression through changes in diet quality or sodium intake.
The association between TRPV1-C allele presence and HADS-D scores indicative of depression (with and without adjustments) was sex specific, with significant findings in males but not females. This is not due to differences in sex-specific distributions in TRPV1-C allele carriage or HADS-D scores as these did not vary between groups. However, females had higher diet quality scores on all three indices than males, and females with lower diet quality scores (DGI and ARFS) were more likely to have indicative depression while there were no associations between diet quality and depression in males. This indicates there may be a threshold effect where the impacts of genotypic variance manifest in the context of poor diet quality in males, or a sex hormone-specific interaction may be occurring.
As a missense mutation [99], TRPV1-rs8065080 may be modifying protein structure and function [57, 58, 100‒103]. Neuropathic pain irregularities such as higher tolerance to cold pain [101, 102], lower risk of wheezing in childhood asthma [100], and protective effects against hypertension (TRPV1-C allele) [103] are associated with the SNP. TRPV1-rs8065080 has been shown to decrease TRPV1 channel activity by approximately 20–30% after heat and capsaicin stimulation [100]. Decreased salt sensitivity has been found in TRPV1-C allele carriers [57] and they have been found to rate salt solutions as less intense than non-carriers [58]. Therefore, as TRPV1-rs8065080 likely causes a reduction in taste function [57, 58], and the current study suggests there is a greater likelihood of having HADS-D scores indicative of depression in the presence of the SNP, the results indicate that people with genes that reduce salt-taste function are more likely to experience depression.
The findings in this study indicate that the hypothesized functional changes in the TRPV1-rs8065080 receptor may have extra-oral implications aside from changing taste. Sodium intake was not a significant determinant of HADS-D threshold scores, TRPV1-C allele presence did not influence sodium intake, and sodium intake did not impact on the association between presence of the TRPV1-C allele and HADS-D scores. Therefore, the analyses suggest that TRPV1-rs8065080 genotype is driving the association with depression independent of sodium intake and that the functional changes in the SNP may be occurring extra-orally. TRPV1 mRNA is expressed in the human brain including in the insula, amygdala, orbitofrontal cortex, thalamus, and hippocampus [63, 64]. These same brain regions demonstrate structural and functional abnormalities in human depression [67, 68]. The missense mutation may be creating functional changes in TRPV1 and the presence of these receptors in regions of the brain affected by depression could be an avenue for further exploration.
There is an additional biological pathway to consider as a possible mechanism involved in the findings of this study and for further research. Salt triggers the release of 5-HT and NE from type III taste cells [104, 105]. These neurotransmitters are found to be low in depression [106] and to downregulate depression [74, 75]. Salt also stimulates the release of anandamide which activates TRPV1 [107, 108] and restores 5-HT and NE levels in the brain [109, 110]. Absence of TRPV1 in knockout mice shows changes in 5-HT receptors and decreased depressive behaviours [69]. Therefore, the association between TRPV1-C allele carriage and higher likelihood of scores indicative of depression may be a result of SNP-related functional changes to levels of neurotransmitters linked to depression.
As the nature of this study is observational and mental health conditions and genetic expression are multifactorial, the findings are limited to association. Additional limitations include those around the collection of subjective data such as the HADS test and FFQ. For example, rates of depression and anxiety in the study sample differed to those of similar age in the general population at the time of data collection (2010–2012). HADS-A threshold scores occurred in 13.4% and HADS-D threshold scores in 8.0% of the study population. In comparison, the Australian Health Survey data from 2011 to 2012 reported that approximately 3.3% of people over 65 experienced anxiety and 10.6% experienced depression [111]. Therefore, significant associations to depression may be under-stated.
The study’s strengths include the large cohort size and availability of broad demographic and dietary data that enabled adjustments to be made for sex, age, income, education, smoking, BMI, sodium intake, and three diet quality indices. As age was not a confounding factor in the significant associations, the results can be more widely considered in other age groups. The sex distribution is reflective of the Australian statistics for people aged 65 years and older [112] and HADS is a reliable and valid measure of levels of depression and anxiety in research settings [91]. While cross-sectional studies have limitations, they are a necessary first step towards experimental studies to test the results under controlled conditions. Future studies should investigate how salt-taste receptors, TRPV1 channels, and SNPs in the brain are involved in the biology of depression, which may generate new approaches to prevention and treatment. Studies that elucidate the role neurotransmitters have in salt taste and mental health and that determine if salt-taste SNPs alter neurotransmitter levels may also contribute to improved management of depression.
Conclusions
Overall, the novel findings of this exploration of associations between genetic variance of a salt-taste receptor (TRPV1) and levels of anxiety and depression in an elderly population suggest that the TRPV1-rs8065080 SNP is more likely to be present when depression is indicated. This association appears to be independent of sodium intake and diet quality in the total sample. However, poor diet quality in males may provide the context of a threshold effect for the impacts of the gene. This knowledge may aid in the development of improved risk assessment or diagnosis of depression and the development of preventative strategies and new therapies. Given the significant and increasing disease burden of depression, further research is needed to explore these associations to improve outcomes and the quality of life of sufferers.
Acknowledgments
Authors acknowledge the role of Paul Roach in the design of the original study and the roles of Charlotte Martin, Zoe Yates, Katrina King, and Suzanne Niblett in sample collection and data management.
Statement of Ethics
Ethics approval for the Retirement Health and Lifestyle Study was granted by the Human Research Ethics Committee of the University of Newcastle (Reference No. H-2008-0431). Written informed consent was provided by all subjects included in this study.
Conflict of Interest Statement
There are no conflicts of interest to declare.
Funding Sources
This research was conducted as part of the Retirement Health and Lifestyle Study, with initial and ongoing funding provided by the Australian Research Council (G0188386), Central Coast Local Health District Public Health Unit (G0190658/G1700259), UnitingCare Ageing NSW/ACT (G0189230), Urbis Pty Ltd. (G0189232), Valhalla Village Pty Ltd. (G1000936), and Hunter Valley Research Foundation.
Author Contributions
Conceptualization, formal analysis, methodology, and writing – original draft: C.F. and E.L.B.; data curation: C.F., M.L., and E.L.B.; funding acquisition, project administration, and resources: M.V., M.L., and E.L.B.; and investigation and writing – review and editing: C.F., C.S., M.V., M.L., T.B., and E.L.B. All authors have read and agreed to the published version of the manuscript.
Data Availability Statement
All data generated or analysed during this study are included in this article. Further enquiries can be directed to the corresponding author.