Skip to main content
  • Research article
  • Open access
  • Published:

Starch intake, amylase gene copy number variation, plasma proteins, and risk of cardiovascular disease and mortality

Abstract

Background

Salivary amylase, encoded by the AMY1 gene, initiate the digestion of starch. Whether starch intake or AMY1 copy number is related to disease risk is currently rather unknown. The aim was to investigate the association between starch intake and AMY1 copy number and risk of cardiovascular disease (CVD) and mortality and whether there is an interaction. In addition, we aim to identify CVD-related plasma proteins associated with starch intake and AMY1 copy number.

Methods

This prospective cohort study used data from 21,268 participants from the Malmö Diet and Cancer Study. Dietary data were collected through a modified diet history method and incident CVD and mortality were ascertained through registers. AMY1 gene copy number was determined by droplet digital polymerase chain reaction, a risk score of 10 genetic variants in AMY1 was measured, and a total of 88 selected CVD-related proteins were measured. Cox proportional hazards regression was used to analyze the associations of starch intake and AMY1 copy number with disease risk. Linear regression was used to identify plasma proteins associated with starch intake and AMY1 copy number.

Results

Over a median of 23 years’ follow-up, 4443 individuals developed CVD event and 8125 died. After adjusting for potential confounders, a U-shape association between starch intake and risk of CVD (P-nonlinearity = 0.001) and all-cause mortality (P-nonlinearity = 0.03) was observed. No significant association was found between AMY1 copy number and risk of CVD and mortality, and there were no interactions between starch intake and AMY1 copy number (P interaction > 0.23). Among the 88 plasma proteins, adrenomedullin, interleukin-1 receptor antagonist protein, fatty acid-binding protein, leptin, and C-C motif chemokine 20 were associated with starch intake after adjusting for multiple testing.

Conclusions

In this large prospective study among Swedish adults, a U-shaped association between starch intake and risk of CVD and all-cause mortality was found. Several plasma proteins were identified which might provide information on potential pathways for such association. AMY1 copy number was not associated with CVD risk or any of the plasma proteins, and there was no interaction between starch intake and AMY1 copy number on disease risk.

Peer Review reports

Background

Cardiovascular disease (CVD), a major cause of disability and premature death throughout the world, can be prevented to a large extent by a healthy diet [1]. Evidence on the role of carbohydrate intake on risk of CVD and mortality have been inconsistent [2,3,4,5,6], which could due to opposite effects by different types of carbohydrates. For example, sugar intake has been positively associated, while fiber intake has been inversely associated, with the various disease outcomes [7, 8]. Compared with sugar and fiber intake, the association between starch intake and risk of CVD and mortality have been investigated to a limited extent and the results are inconsistent [9,10,11]. The inconsistency could be due to genetic variation across populations.

Starch digestion is initiated in the oral cavity by the action of salivary amylase, which is encoded by the salivary amylase gene (AMY1) [12]. Although foods remain in the mouth for a relatively short time, salivary amylase bound to starch can also be detected in the duodenum [13]. The copy number variation (CNV) of the AMY1 gene, ranging from 2 to 17 diploid copies, is correlated to salivary amylase concentrations and was found to be higher in populations that historically have consumed large amounts of starch [12, 14]. A modestly higher postprandial glycemia after a starchy meal has been shown among individuals with high AMY1 copy number compared to those with low AMY1 copy number [15, 16], although some studies found no difference [17, 18]. In addition, the association between AMY1 copy number and metabolic traits is inconsistent [17, 19,20,21], although a negative association between AMY1 copy number and several cardiometabolic risk factors, including obesity and inflammation which have been shown in some studies, mainly will be small sample sizes [22,23,24]. Thus, there is a need to further explore the association between AMY1 copy number and cardiometabolic risk markers, especially to simultaneously examine a broad range of risk markers, including inflammatory markers. In addition, to the best of our knowledge, there are no studies that have investigated the association between AMY1 copy number and risk of CVD and mortality. Although we did not find significant associations between AMY1 copy number and body mass index (BMI) or glucose homeostasis traits in our previous study, we observed an interaction between AMY1 copy number and starch intake on BMI and insulin resistance [20, 21]. For example, the inverse association of high-starch diet and fasting insulin levels was stronger with higher AMY1 copy numbers [20].

Therefore, the aim of this study was to investigate the association between starch intake, AMY1 copy number, and risk of CVD and mortality and to examine the interaction between starch intake and AMY1 copy number. In addition, we aim to identify whether any of 88 CVD-related plasma proteins were associated with starch intake or AMY1 copy number.

Methods

Study population

The Malmö Diet and Cancer Study (MDCS) is a population-based prospective cohort study initiated to examine associations between diet and health outcomes [25]. It includes 30,446 individuals from Malmö, in the south of Sweden, aged 44–74 years at enrollment between 1991 and 1996. At baseline, all the MDCS participants underwent physical examinations, had blood taken, and answered questionnaires regarding lifestyle factors. Between October 1991 and February 1994, 6103 individuals was randomly selected to be studied for the epidemiology of carotid artery disease, and this sample is referred to as the Malmö Diet and Cancer-Cardiovascular Cohort (MDC-CC) [26]. All participants provided written informed consent. The study was approved by the Ethical Review Board at the Faculty of Medicine at Lund University, Sweden (LU 51/90), and was carried out in accordance with the Helsinki Declaration.

Among the 30,446 participants, we excluded those with missing information on dietary intake (n = 2213) and other relevant covariates (n = 563). We also excluded participants with a history of CVD (n = 829), diabetes (n = 1123), and cancer (n = 1574) at baseline, leaving 24,144 individuals. To reduce the bias of population stratification for the genetic analyses, we further excluded participants who were not born in Sweden or with unknown country of birth (n = 2876), leaving 21,268 individuals. A total of 20,264 individuals had available data on specific genetic variants within the AMY1 gene, and 3932 had information on AMY1 copy number.

Dietary assessment and calculation of starch intake

Dietary intake was assessed using a modified diet history method consisting of (1) a 7-day food diary covering all cooked meals, cold beverages, and dietary supplements; (2) a 168-item food-frequency questionnaire (FFQ) for assessment of consumption frequencies and portion sizes of regularly eaten foods that were not covered by the 7-day food diary; and (3) an interview to ask for cooking methods and usual portion sizes for foods recorded in the food diary and to check for overlap between intakes reported by the food diary and the FFQ. Daily food intake was calculated by combing the intake data from the FFQ and the food diary. The energy and nutrient intake values were obtained from the MDCS Food and Nutrient Database that originated from PC KOST2-93, based on the Swedish National Food Agency. Starch intake was calculated by subtracting the intake of total sugars (i.e., mono- and disaccharides) from the total carbohydrate intake and was expressed as the percentage of total non-alcoholic energy intake (E%). The participants were divided into sex-specific quartiles based on their starch intake. Data on the validity [27, 28] and reproducibility [29] of the method have been published. The energy-adjusted Pearson correlation coefficients for carbohydrate intake were 0.66 and 0.70 for men and women, respectively.

Genotyping

AMY1 gene copy number was determined by droplet digital polymerase chain reaction (ddPCR) using the QX200 AutoDG ddPCR System (Bio-Rad Laboratories, Hercules, CA) following the manufacturer’s instructions, as previously described in detail by Rukh et al. [21]. Samples were analyzed along with negative control (RNase free water) on each plate. A high copy number control (NA18972; HapMap DNA: 17 copies of AMY1) and 2 samples from female donors with known AMY1 copy numbers (4 and 8 copies of AMY1) were included as controls in each run. Samples with high copy numbers were diluted and reanalyzed to exclude the possibility of the presence of too much DNA or a technical error. After amplification, the reactions were stored at 4 °C until the plates were read on a Bio-Rad QX200 Droplet Reader, and copy numbers were calculated through the use of the QuantaSoft software version 1.7.4 (Bio-Rad Laboratories). Copy number of AMY1 was divided into four groups: 1–4, 5–6, 7–9, and 10 and above copies.

We created AMY1-genetic risk scores (GRS) using 10 single nucleotide polymorphisms (SNPs) associated with copy number variation in AMY1 gene loci based on a previous publication of European data [30]. The relationship of these SNPs to AMY1 copy number was consistent across 2 independent cohorts of European-ancestry individuals sampled in the USA and Europe [30]. Individual SNP was coded as 0, 1, and 2 according to the number of increasing alleles. We calculated a weighted GRS using the following equation: AMY1-GRS = (β1 × SNP1 + β2 × SNP2 +…+ β3 × SNP3) × (n/sum of the β coefficients), where β is the β coefficient of each SNP for change in AMY1 copy number in the MDCS, and n was 10. The spearman correlation coefficient between AMY1 copy number and AMY1-GRS is 0.43 in our cohort. Details of SNPs and β coefficients of each SNP are presented in Additional file 1: Table S1.

Plasma protein quantification

For plasma proteins, the analysis was restricted to the MDC-CC participants. The available sample was 3680 and 3254 participants for identifying the proteins associated with starch intake and AMY1 copy number, respectively. Serum samples were stored frozen until proteomic profiling was performed using the Olink platform (Olink Proteomics, Uppsala, Sweden) (https://www.olink.com/content/uploads/2015/12/0696-v1.3-Proseek-Multiplex-CVD-I-Validation-Data_final.pdf). A total of 92 plasma proteins were measured by the SciLifeLab with the Olink Proseek Multiplex proximity extension assays Cardiovascular. Values were expressed as normalized protein expression values as arbitrary units on a log2 scale. For statistical analysis, we excluded four proteins that were available in less than 75% of the individuals in the present study sample. Hence, 88 proteins were available for the analysis.

Outcome

All participants were followed up from baseline examination to the diagnosis of CVD, emigration from Sweden, death, or December 31, 2018, whichever came first. Total CVD and mortality were obtained through the Swedish Hospital Discharge Register [31], the Swedish Cause of Death Register, and the Stroke Register of Malmö [32]. Total CVD was defined as a hospital admission or death using International Classification of Diseases (ICD) codes, including coronary heart disease (ICD-9 codes 410A-410X and ICD-10 code I21), death attributable to ischemic heart disease (ICD-9 codes: 410-414; ICD-10 codes: I20-I25), and ischemic stroke (ICD-9 code 434 and ICD-10 code I63). CVD mortality was defined as ICD-9 codes 390-459 and ICD-10 I00-I99.

Assessment of covariates

Data on smoking habits (current, former, or never), educational level (the highest qualification attained), and heredity for diseases (including cancer, myocardial infarction, stroke, and diabetes) and other lifestyle factors were obtained using a structured questionnaire. Four separate heredity scores were constructed and participants were categorized as having low heredity (e.g., with no first-degree relatives with the disease), medium heredity (with one first-degree relative with the disease), and high heredity (with at least two first-degree relatives with the disease) for disease [33]. Alcohol consumption was divided into six categories. Zero consumers reported no consumption of alcohol in the 7-day food diary or during the previous year in the questionnaire. The other individuals were divided into sex-specific quintiles based on the reported alcohol intake from the 7-day food diary. Leisure-time physical activity was categorized into five predefined groups based on metabolic equivalent task (MET) hours per week: < 7.5, 7.5–15, 15–25, 25–50, and > 50 MET hour/week. BMI (kg/m2) was calculated as weight in kilograms divided by the square of the height in meters. Systolic and diastolic blood pressures were measured after 10 min of rest in the supine position. Hypertension was defined as systolic/diastolic blood pressures ≥ 140/90 mmHg and/or current use of antihypertensive medications. The variable for dietary assessment method was created based on the interview time of 60 min before September 1994 and 45 min from 1 September 1994. The season of diet data collection was defined as spring, summer, autumn, and winter. A diet quality index based on the Swedish dietary guidelines was calculated using the following factors: fiber (> 2.4 g/MJ), fruit and vegetables (> 400 g/day), fish (> 300 g/week), added sugar (< 10% energy), and red and processed meat (< 500 g/week). Each point was given for each favorable diet factor, with the total score ranging from 0 to 5 [34].

Statistical analysis

To evaluate differences in baseline characteristics for participants with different quartiles of starch intake, 𝜒2 test was conducted for categorical variables, and analysis of variance was conducted for continuous variables. Cox proportional hazards models were conducted with quartiles of sex-specific starch intake as the exposures of interest, incident CVD and mortality as the outcome, and time-to-event as the time metric. We tested for linear trends by using the ordinal score on quartiles of starch intake. Three models were applied in our analyses. Model 1 was adjusted for age, sex, dietary assessment method, season, and total energy intake. Model 2 was additionally adjusted for leisure-time physical activity, smoking status, alcohol consumption, education, heredity score, and diet quality index. Model 3 was further adjusted for hypertension and BMI. Interaction by cereal intake, bread intake, potato intake, and sex were performed by adding a multiplicative factor between these starch food source (quartiles) and starch (continues) in model 3.

Age- and sex-adjusted Cox proportional hazards models were used to examine the associations of AMY1 copy number and AMY1-GRS with risk of CVD and mortality. Interaction by age were performed by adding a multiplicative factor between age and AMY1 copy number in age- and sex-adjusted model. Furthermore, multiplicative interactions between AMY1 copy number or AMY1-GRS and starch intake were evaluated using the Wald test on cross-product terms into multivariable models (the above model 3). The association between starch intake and risk of CVD and mortality were evaluated in strata of AMY1 copy number (1–4 copies, 5–6 copies 7–9 copies, 10 and above copies) and AMY1-GRS (quintiles). We verified the assumption of linearity between starch intake, AMY1 copy number, AMY1-GRS, and risk of CVD and mortality using restricted cubic spline functions with the SAS macro written by Desquilbet and Mariotti [35].

Associations between starch intake and plasma proteins were studied with three multivariable linear regression models: model 1 adjusted for age, sex, season, and total energy intake, model 2 with further adjustment for education, smoking status, alcohol consumption, and leisure-time physical activity, and model 3 with additional adjustment for BMI. Associations between AMY1 copy number and plasma proteins were studied with multivariable linear regression, adjusted for age and sex. Tests of statistical significance were corrected for multiple comparisons using the Bonferroni method for 88 tests (P < 0.00057).

Cox proportional hazards models were used to study associations between proteins that were associated with starch intake or AMY1 copy number (after Bonferroni correction) and risks of CVD and mortality. Models 1 were adjusted for age and sex, and model 2 were further adjusted for physical activity, smoking status, alcohol consumption, educational level, and heredity scores. To test for the potential influence of obesity in the association between plasma protein and risk of CVD and mortality, we additionally adjusted this model for BMI (model 3). All tests were two-sided, and we considered P < 0.05 to be statistically significant. SAS version 9.4 (SAS Institute Inc., Cary, NC, USA) was used for the analyses.

Results

A total of 21,268 participants (8187 (38.5%) men and 13,081 (61.5%) women) were included in the study. At baseline, the mean intake of starch in the total population was 24E% or 134 g/day.

Table 1 shows the baseline characteristics of participants according to quartiles of sex-specific starch intake. Compared with participants in the lowest quartile, participants in the highest quartile of starch intake tended to be younger, more often ex-smokers and non-smokers, and less likely to have hypertension, but were more likely to have family history of myocardial infarction and diabetes. They had lower intake of energy and fat and higher intake of carbohydrates, protein, potato, and bread and higher diet quality score.

Table 1 Baseline characteristics of the participants according to sex specific quartiles of starch intake (n = 21,268)

During median follow-up of 23 years, 4443 (20.9%) of the participants had a cardiovascular event and 8125 (38.2%) died. Table 2 shows the association between starch intake and the risk of CVD and mortality. The restricted cubic spline analyses indicated U-shape associations between starch intake and CVD risk (P-nonlinearity = 0.0012) and all-cause mortality (P-nonlinearity = 0.03) (Fig. 1). Compared with those in the lowest quartile, those in the second to the third quartile had a 9–11% lower risk of CVD, while there was no statistically reduced risk in the highest quartile (Table 2). Participants in the highest vs. lowest quartile had a hazard ratio of 0.91 (95% confidence interval 0.85 to 0.97) for all-cause mortality. Similar non-linear associations (P-nonlinearity = 0.08–0.11), although not statistically significant were observed between starch intake and CHD, ischemic stroke, and CVD mortality (Fig. 1). There was no heterogeneity of effect for cereal intake, bread intake, potato intake, and sex on association of starch intake and CVD (all P-interaction > 0.07) and mortality risk (all P-interaction > 0.07).

Table 2 Associations between intake of starch and CVD and mortality risk (n = 21,268)
Fig. 1
figure 1

Restricted cubic spline plots to assess association between starch intake and CVD and mortality. The HRs and 95% CIs above were adjusted for age, sex, smoking status, drinking status, education, season, method, physical activity, heredity scores (including cancer, infarct, stroke, and diabetes), total energy intake, the modified diet index, hypertension, and body mass index. CVD, cardiovascular disease; CHD, coronary heart disease

We did not observe statistically significant associations between AMY1 copy number and any of the outcomes (Table 3 and Additional file 1: Fig. S2). There was no heterogeneity of effect for age on association of AMY1 copy number and CVD (P-interaction = 0.07) and mortality risk (P-interaction = 0.14). In addition, no significant associations were found between AMY1-GRS and any of the outcomes (Additional file 1: Table S2 and Fig. S3). There was no evidence of interaction between starch intake and AMY1 copy number or AMY1-GRS on CVD and mortality risk (Fig. 2 and Additional file 1: Fig. S4). Neither AMY1 copy number (P = 0.58) nor AMY1-GRS (P = 0.85) were associated with starch intake.

Table 3 Associations between intake of AMY1 copy number and CVD and mortality risk (n = 3932)
Fig. 2
figure 2

Association between intake of starch and risk of CVD and mortality by AMY1 copy numbera. aThe multivariable analysis was adjusted for age, sex, smoking status, drinking status, education, season, method, physical activity, heredity scores (including cancer, infarct, stroke, and diabetes), total energy intake, the modified diet index, hypertension, and body mass index. bHazard ratio for per increase of 10% in proportion of starch energy in total energy. CVD, cardiovascular disease; CHD, coronary heart disease; CI, confidence interval; HR, hazards ratio

After adjustment for potential confounders and correction for multiple testing, concentrations of the following five among the 88 plasma proteins were found to be significantly and inversely associated with starch intake: adrenomedullin, interleukin-1 receptor antagonist protein (IL1ra), fatty acid-binding protein (FABP4), leptin, and C-C motif chemokine 20 (CCL20) (Fig. 3). The associations were attenuated and not statistically significant after further adjustment for BMI (Additional file 1: Table S3). All five starch-related proteins were positively associated with risk of CVD, CHD, and all-cause mortality in model 2, adjusting for lifestyle factors (P < 0.05). Higher levels of leptin, IL1ra, and adrenomedullin were associated with increased ischemic stroke risk. Higher levels of leptin, FABP4, and IL1ra were associated with increased CVD mortality. Overall, the associations were not significant after further adjustment for BMI (Additional file 1: Table S4). None of the plasma proteins were associated with the AMY1 copy number with adjustment for age and sex (Additional file 1: Table S5).

Fig. 3
figure 3

Volcano plot of association between starch intake and 88 plasma proteins in full sample analysis (n = 3680). Linear regressions were adjusted for age, sex, season, total energy intake, educational level, smoking status, alcohol consumption, and physical activity

Discussion

In this large prospective cohort of Swedish middle aged and elderly individuals, we found non-linear associations between starch intake and risk of CVD and all-cause mortality. We identified five plasma proteins linked to lower starch intake that were also associated with increased CVD and mortality risk, providing clues to underlying biological mechanisms. AMY1 copy number was not associated with plasma proteins or CVD and mortality risk, and we did not observe any interaction between AMY1 copy number and starch intake on risk of CVD and mortality.

The observed U-shaped associations between starch intake and risk of CVD and all-cause mortality observed in our study was similar to previous results in the prospective UK Biobank cohort, which reported lowest risk of CVD and death among participants consuming 20–30E% from starch [10]. In contrast, a prospective cohort study conducted in Japan reported an inverse association between starch intake and all-cause mortality in men [11]. These studies, including ours, have adjusted for total energy intake in the analyses and thereby examining the isocaloric exchange of starch with other macronutrients. This discrepancy might be due to the different levels and sources of starch in Japan and Europe. For example, the mean intake of starch is approximately 39E% in the Japanese study and 24E% both in the UK Biobank study and our study. Furthermore, the main source of starch in Japan is rice, while in Europe, it is potato and bread. In a meta-analysis, high rice consumption was related to a modest reduction in risk of mortality in men [36]. The mean intake of starch in our population was 134 g/day, which is approximately equal to the starch content in 300 g of bread or 1 kg of potatoes [37]. Previous studies have shown a protective association of total bread with mortality [38]. However, a recent meta-analysis did not find a significant association between potato consumption and risk of all-cause mortality [39]. In addition, although not statistically significant, a U-shaped association was also found between starch intake and risk of CHD and stroke in our study. Similarly, previous prospective studies have not observed the association between starch intake and risk of CHD among US, Chinese, and European populations [9, 40,41,42] and risk of stroke among Europeans [43].

We identified adrenomedullin, IL1ra, FABP4, leptin, and CCL20 to be inversely associated with starch intake, and all these proteins were positively associated with increased risk of CVD and mortality. Plasma FABP4, leptin, IL1ra, and CCL20 play a regulatory role in inflammation [44, 45], which is an underlying mechanism of CVD and death. Furthermore, adrenomedullin, leptin, and FABP4 are secreted by the adipose tissue in direct relation to the amount of body fat [46,47,48], and IL1RA and CCL20 has been shown to be highly elevated in obese animals and humans [49,50,51]. Animal experiments and human trials showed that resistant starch (i.e., starch that is not digested and absorbed in the small intestine) had an anti-inflammatory effect [52]. In addition, several studies have shown that higher intakes of slowly digestible starch and resistant starch intake are associated with increased satiety, reduced hunger, and/or reduced body weight [53, 54]. Although it is not easy to explain the association with starch, high intakes may inhibited the expression and function of plasma proteins by increasing insulin resistance and plasma-free fatty acids [55], as well as delaying adipocyte differentiation [56].

In the present study, AMY1 copy number was not associated with any of the outcomes, and as far as we know, there are no previous study that has investigated this. AMY1 copy number was only measured in 3932 of the participants and AMY1-GRS was measured in 20,264 of the participants. Thus, the size of our cohort may have been underpowered to detect small differences and future study needs to confirm such association in other populations. Because several previous studies showed no or weak associations with CVD risk factors (e.g., inflammation, obesity, and insulin sensitivity) [22, 24, 57], we might not expect any major association of AMY1 copy number and CVD risk.

We did not identify any of the included plasma proteins, including several inflammatory markers, to be associated with AMY1 copy number after correction for multiple testing. To our knowledge, the current study is the first to examine a broad range of CVD-related plasma proteins, including inflammatory markers (e.g., interleukin (IL), tumor necrosis factor, and CCL) associated with AMY1 copy number. A previous study of 57 non-diabetic adults found that individuals with low AMY1 CNV had higher concentrations of serum cytokines, including IL-6, IL-1β, monocyte chemoattractant protein-1, and tumor necrosis factor-alpha [22]. Another study of 76 children found that AMY1 CNV were negatively associated with the C-reactive protein, resistin, and CCL2 [58].

The strengths of the study include its prospective study design, long follow-up time, detailed dietary assessment method, and more precise measurement of AMY1 copy number using ddPCR. Nevertheless, several limitations must be acknowledged. Firstly, dietary intakes were collected only at baseline and the dietary habits may have changed over time, which could have potentially underestimated the diet-disease associations. Second, our study only includes individuals of European descent. Thus, caution should be taken when generalizing the findings to other populations. In particular, the consumption of starch could be substantially different in other regions, and there might be diversity in the distribution of genes across different populations. Third, the lack of replication in external cohorts of the plasma proteins in relation to starch intake is a potential concern. Fourth, we were unable to investigate associations with different structures of starch because these data were not available, which could be an effect modifier. Therefore, future studies should focus on the influence of different types of starch on CVD and mortality. Finally, as in other observational studies, even though we have adjusted for a wide variety of covariates that relate to CVD and mortality, residual confounding by unidentified confounders is still possible.

Conclusions

In conclusion, we observed U-shaped associations of starch intake with risk of CVD and mortality, and the plasma proteomic analysis which identified several plasma proteins associated with starch intake and CVD risk could provide clues to underlying biological mechanisms. AMY1 copy number was not associated with risk of CVD, mortality, or CVD-related plasma proteins, and there was no interaction between starch intake and AMY1 copy number.

Availability of data and materials

The data that support the findings of this study are available from the Malmö Population-Based Cohorts Joint Database but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of the Malmö Population-Based Cohorts Joint Database.

Abbreviations

BMI:

Body mass index

CCL20:

C-C motif chemokine 20

CNV:

Copy number variation

CVD:

Cardiovascular disease

E%:

Percentage of total non-alcoholic energy intake

FABP4:

Fatty acid-binding protein

FFQ:

Food frequency questionnaire

GRS:

Genetic risk score

ICD:

International Classification of Diseases

IL1ra:

Interleukin-1 receptor antagonist protein

IL:

Interleukin

MDCS:

The Malmö Diet and Cancer Study

MDC-CC:

The Malmö Diet and Cancer-Cardiovascular Cohort

MET:

Metabolic equivalent task

SNPs:

Single nucleotide polymorphisms

References

  1. Yu E, Malik VS, Hu FB. Cardiovascular disease prevention by diet modification: JACC Health Promotion Series. J Am Coll Cardiol. 2018;72(8):914–26.

    Article  Google Scholar 

  2. Gribbin S, Enticott J, Hodge AM, Moran L, Thong E, Joham A, et al. Association of carbohydrate and saturated fat intake with cardiovascular disease and mortality in Australian women. Heart. 2022;108(12):932–9.

    Article  CAS  Google Scholar 

  3. Trichopoulou A, Psaltopoulou T, Orfanos P, Hsieh CC, Trichopoulos D. Low-carbohydrate-high-protein diet and long-term survival in a general population cohort. Eur J Clin Nutr. 2007;61(5):575–81.

    Article  CAS  Google Scholar 

  4. Lagiou P, Sandin S, Lof M, Trichopoulos D, Adami HO, Weiderpass E. Low carbohydrate-high protein diet and incidence of cardiovascular diseases in Swedish women: prospective cohort study. BMJ. 2012;344:e4026.

    Article  Google Scholar 

  5. Sjogren P, Becker W, Warensjo E, Olsson E, Byberg L, Gustafsson IB, et al. Mediterranean and carbohydrate-restricted diets and mortality among elderly men: a cohort study in Sweden. Am J Clin Nutr. 2010;92(4):967–74.

    Article  Google Scholar 

  6. Dehghan M, Mente A, Zhang X, Swaminathan S, Li W, Mohan V, et al. Associations of fats and carbohydrate intake with cardiovascular disease and mortality in 18 countries from five continents (PURE): a prospective cohort study. Lancet. 2017;390(10107):2050–62.

    Article  CAS  Google Scholar 

  7. Anderson JJ, Gray SR, Welsh P, Mackay DF, Celis-Morales CA, Lyall DM, et al. The associations of sugar-sweetened, artificially sweetened and naturally sweet juices with all-cause mortality in 198,285 UK Biobank participants: a prospective cohort study. BMC Med. 2020;18(1):97.

    Article  CAS  Google Scholar 

  8. Reynolds A, Mann J, Cummings J, Winter N, Mete E, Te Morenga L. Carbohydrate quality and human health: a series of systematic reviews and meta-analyses. Lancet. 2019;393(10170):434–45.

    Article  CAS  Google Scholar 

  9. Liu S, Willett WC, Stampfer MJ, Hu FB, Franz M, Sampson L, et al. A prospective study of dietary glycemic load, carbohydrate intake, and risk of coronary heart disease in US women. Am J Clin Nutr. 2000;71(6):1455–61.

    Article  CAS  Google Scholar 

  10. Ho FK, Gray SR, Welsh P, Petermann-Rocha F, Foster H, Waddell H, et al. Associations of fat and carbohydrate intake with cardiovascular disease and mortality: prospective cohort study of UK Biobank participants. BMJ. 2020;368:m688.

    Article  Google Scholar 

  11. Nagata C, Wada K, Yamakawa M, Konishi K, Goto Y, Koda S, et al. Intake of starch and sugars and total and cause-specific mortality in a Japanese community: the Takayama Study. Br J Nutr. 2019;122(7):820–8.

    Article  CAS  Google Scholar 

  12. Perry GH, Dominy NJ, Claw KG, Lee AS, Fiegler H, Redon R, et al. Diet and the evolution of human amylase gene copy number variation. Nat Genet. 2007;39(10):1256–60.

    Article  CAS  Google Scholar 

  13. Fried M, Abramson S, Meyer JH. Passage of salivary amylase through the stomach in humans. Dig Dis Sci. 1987;32(10):1097–103.

    Article  CAS  Google Scholar 

  14. Mandel AL, C Peyrot des Gachons, Plank KL, Alarcon S, Breslin PA: Individual differences in AMY1 gene copy number, salivary alpha-amylase levels, and the perception of oral starch. PLoS One 2010, 5(10):e13352.

  15. Atkinson FS, Hancock D, Petocz P, Brand-Miller JC. The physiologic and phenotypic significance of variation in human amylase gene copy number. Am J Clin Nutr. 2018;108(4):737–48.

    Article  Google Scholar 

  16. Carpenter D, Dhar S, Mitchell LM, Fu B, Tyson J, Shwan NA, et al. Obesity, starch digestion and amylase: association between copy number variants at human salivary (AMY1) and pancreatic (AMY2) amylase genes. Hum Mol Genet. 2015;24(12):3472–80.

    Article  CAS  Google Scholar 

  17. Mandel AL, Breslin PA. High endogenous salivary amylase activity is associated with improved glycemic homeostasis following starch ingestion in adults. J Nutr. 2012;142(5):853–8.

    Article  CAS  Google Scholar 

  18. Gigliola Alberti JP, Cataldo LR, Vega J, Aguilera CM, Alvarez-Mercado AI, Hodgson MI, et al. Glycemic response after starch consumption in relation to salivary amylase activity and copy-number variation of AMY1 gene. J Food Nutr Res. 2015;3(8, 2015):558–63.

    Google Scholar 

  19. Nakajima K, Nemoto T, Muneyuki T, Kakei M, Fuchigami H, Munakata H. Low serum amylase in association with metabolic syndrome and diabetes: A community-based study. Cardiovasc Diabetol. 2011;10:34.

    Article  CAS  Google Scholar 

  20. Hamid AK, Andersson-Assarsson J, Ericson U, Sonestedt E. Interaction effect between copy number variation in salivary amylase locus (AMY1) and starch intake on glucose homeostasis in the Malmo Diet and Cancer cohort. Front Nutr. 2020;7:598850.

    Article  Google Scholar 

  21. Rukh G, Ericson U, Andersson-Assarsson J, Orho-Melander M, Sonestedt E. Dietary starch intake modifies the relation between copy number variation in the salivary amylase gene and BMI. Am J Clin Nutr. 2017;106(1):256–62.

    Article  CAS  Google Scholar 

  22. Marquina C, Mousa A, Belski R, Banaharis H, Naderpoor N, de Courten B. Increased inflammation and cardiometabolic risk in individuals with low AMY1 copy numbers. J Clin Med. 2019;8(3):382.

    Article  CAS  Google Scholar 

  23. Mayneris-Perxachs J, Mousa A, Naderpoor N, Fernandez-Real JM, de Courten B. Low AMY1 copy number is cross-sectionally associated to an inflammation-related lipidomics signature in overweight and obese individuals. Mol Nutr Food Res. 2020;64(11):e1901151.

    Article  Google Scholar 

  24. Grammatikopoulou MG, Gkiouras K, Markaki AG, Gkouskou KK, Aivaliotis M, Stylianou K, et al. AMY1 diploid copy number among end-stage renal disease patients. Hormones (Athens). 2020;19(3):369–76.

    Article  Google Scholar 

  25. Manjer J, Carlsson S, Elmstahl S, Gullberg B, Janzon L, Lindstrom M, et al. The Malmo Diet and Cancer Study: representativity, cancer incidence and mortality in participants and non-participants. Eur J Cancer Prev. 2001;10(6):489–99.

    Article  CAS  Google Scholar 

  26. Muhammad IF, Borne Y, Bao X, Melander O, Orho-Melander M, Nilsson PM, et al. Circulating HER2/ErbB2 levels are associated with increased incidence of diabetes: a population-based cohort study. Diabetes Care. 2019;42(8):1582–8.

    Article  CAS  Google Scholar 

  27. Elmstahl S, Riboli E, Lindgarde F, Gullberg B, Saracci R. The Malmo Food Study: the relative validity of a modified diet history method and an extensive food frequency questionnaire for measuring food intake. Eur J Clin Nutr. 1996;50(3):143–51.

    CAS  Google Scholar 

  28. Riboli E, Elmstahl S, Saracci R, Gullberg B, Lindgarde F. The Malmo Food Study: validity of two dietary assessment methods for measuring nutrient intake. Int J Epidemiol. 1997;26(Suppl 1):S161–73.

    Article  Google Scholar 

  29. Elmstahl S, Gullberg B, Riboli E, Saracci R, Lindgarde F. The Malmo Food Study: the reproducibility of a novel diet history method and an extensive food frequency questionnaire. Eur J Clin Nutr. 1996;50(3):134–42.

    CAS  Google Scholar 

  30. Usher CL, Handsaker RE, Esko T, Tuke MA, Weedon MN, Hastie AR, et al. Structural forms of the human amylase locus and their relationships to SNPs, haplotypes and obesity. Nat Genet. 2015;47(8):921–5.

    Article  CAS  Google Scholar 

  31. Ludvigsson JF, Andersson E, Ekbom A, Feychting M, Kim JL, Reuterwall C, et al. External review and validation of the Swedish national inpatient register. BMC Public Health. 2011;11:450.

    Article  Google Scholar 

  32. Jerntorp P, Berglund G. Stroke registry in Malmo, Sweden. Stroke. 1992;23(3):357–61.

    Article  CAS  Google Scholar 

  33. Johansson A, Drake I, Engstrom G, Acosta S. Modifiable and non-modifiable risk factors for atherothrombotic ischemic stroke among subjects in the Malmo Diet and Cancer study. Nutrients. 2021;13(6):1952.

    Article  CAS  Google Scholar 

  34. Gonzalez-Padilla E, Tao Z, Sanchez-Villegas A, Alvarez-Perez J, Borne Y, Sonestedt E. Association between adherence to Swedish dietary guidelines and Mediterranean diet and risk of stroke in a Swedish population. Nutrients. 2022;14(6):1253.

    Article  Google Scholar 

  35. Desquilbet L, Mariotti F. Dose-response analyses using restricted cubic spline functions in public health research. Stat Med. 2010;29(9):1037–57.

    Google Scholar 

  36. Saneei P, Larijani B, Esmaillzadeh A. Rice consumption, incidence of chronic diseases and risk of mortality: meta-analysis of cohort studies. Public Health Nutr. 2017;20(2):233–44.

    Article  Google Scholar 

  37. Moraeus L, Lemming EW, Hursti UK, Arnemo M, Sipinen JP, Lindroos AK. Riksmaten Adolescents 2016-17: a national dietary survey in Sweden - design, methods, and participation. Food Nutr Res. 2018:62.

  38. Aune D, Keum N, Giovannucci E, Fadnes LT, Boffetta P, Greenwood DC, et al. Whole grain consumption and risk of cardiovascular disease, cancer, and all cause and cause specific mortality: systematic review and dose-response meta-analysis of prospective studies. BMJ. 2016;353:i2716.

    Article  Google Scholar 

  39. Darooghegi Mofrad M, Milajerdi A, Sheikhi A, Azadbakht L. Potato consumption and risk of all cause, cancer and cardiovascular mortality: a systematic review and dose-response meta-analysis of prospective cohort studies. Crit Rev Food Sci Nutr. 2020;60(7):1063–76.

    Article  Google Scholar 

  40. AlEssa HB, Cohen R, Malik VS, Adebamowo SN, Rimm EB, Manson JE, et al. Carbohydrate quality and quantity and risk of coronary heart disease among US women and men. Am J Clin Nutr. 2018;107(2):257–67.

    Article  Google Scholar 

  41. Rebello SA, Koh H, Chen C, Naidoo N, Odegaard AO, Koh WP, et al. Amount, type, and sources of carbohydrates in relation to ischemic heart disease mortality in a Chinese population: a prospective cohort study. Am J Clin Nutr. 2014;100(1):53–64.

    Article  CAS  Google Scholar 

  42. Sieri S, Krogh V, Berrino F, Evangelista A, Agnoli C, Brighenti F, et al. Dietary glycemic load and index and risk of coronary heart disease in a large italian cohort: the EPICOR study. Arch Intern Med. 2010;170(7):640–7.

    Article  CAS  Google Scholar 

  43. Sieri S, Brighenti F, Agnoli C, Grioni S, Masala G, Bendinelli B, et al. Dietary glycemic load and glycemic index and risk of cerebrovascular disease in the EPICOR cohort. PLoS One. 2013;8(5):e62625.

    Article  CAS  Google Scholar 

  44. Fantuzzi G, Faggioni R. Leptin in the regulation of immunity, inflammation, and hematopoiesis. J Leukoc Biol. 2000;68(4):437–46.

    Article  CAS  Google Scholar 

  45. Ge XN, Bastan I, Dileepan M, Greenberg Y, Ha SG, Steen KA, et al. FABP4 regulates eosinophil recruitment and activation in allergic airway inflammation. Am J Phys Lung Cell Mol Phys. 2018;315(2):L227–40.

    CAS  Google Scholar 

  46. Correia ML, Haynes WG. Leptin, obesity and cardiovascular disease. Curr Opin Nephrol Hypertens. 2004;13(2):215–23.

    Article  Google Scholar 

  47. Farooqi IS, Jebb SA, Langmack G, Lawrence E, Cheetham CH, Prentice AM, et al. Effects of recombinant leptin therapy in a child with congenital leptin deficiency. N Engl J Med. 1999;341(12):879–84.

    Article  CAS  Google Scholar 

  48. Maeda K, Cao H, Kono K, Gorgun CZ, Furuhashi M, Uysal KT, et al. Adipocyte/macrophage fatty acid binding proteins control integrated metabolic responses in obesity and diabetes. Cell Metab. 2005;1(2):107–19.

    Article  CAS  Google Scholar 

  49. Meier CA, Bobbioni E, Gabay C, Assimacopoulos-Jeannet F, Golay A, Dayer JM. IL-1 receptor antagonist serum levels are increased in human obesity: a possible link to the resistance to leptin? J Clin Endocrinol Metab. 2002;87(3):1184–8.

    Article  CAS  Google Scholar 

  50. Franck N, Maris M, Nalbandian S, Talukdar S, Schenk S, Hofmann HP, et al. Knock-down of IL-1Ra in obese mice decreases liver inflammation and improves insulin sensitivity. PLoS One. 2014;9(9):e107487.

    Article  Google Scholar 

  51. Zhao L, Xia J, Wang X, Xu F. Transcriptional regulation of CCL20 expression. Microbes Infect. 2014;16(10):864–70.

    Article  CAS  Google Scholar 

  52. Lu J, Ma B, Qiu X, Sun Z, Xiong K. Effects of resistant starch supplementation on oxidative stress and inflammation biomarkers: a systematic review and meta-analysis of randomized controlled trials. Asia Pac J Clin Nutr. 2021;30(4):614–23.

    CAS  Google Scholar 

  53. Willis HJ, Eldridge AL, Beiseigel J, Thomas W, Slavin JL. Greater satiety response with resistant starch and corn bran in human subjects. Nutr Res. 2009;29(2):100–5.

    Article  CAS  Google Scholar 

  54. Aller EE, Abete I, Astrup A, Martinez JA, van Baak MA. Starches, sugars and obesity. Nutrients. 2011;3(3):341–69.

    Article  CAS  Google Scholar 

  55. Kabir M, Guerre-Millo M, Laromiguiere M, Slama G, Rizkalla SW. Negative regulation of leptin by chronic high-glycemic index starch diet. Metabolism. 2000;49(6):764–9.

    Article  CAS  Google Scholar 

  56. Loor J. Transcriptomics of muscle and fat development in growing cattle. Energy and protein metabolism and nutrition. Wageningen: Wageningen Academic Publishers; 2010. p. 59–68.

    Google Scholar 

  57. Mayneris-Perxachs J, Mousa A, Naderpoor N, Fernandez-Real JM, de Courten B. Low AMY1 copy number is cross-sectionally associated to an inflammation-related lipidomics signature in overweight and obese individuals. Mol Nutr Food Res. 2021;65(18):e2170044.

    Article  Google Scholar 

  58. Selvaraju V, Venkatapoorna CMK, Babu JR, Geetha T. Salivary amylase gene copy number is associated with the obesity and inflammatory markers in children. Diabetes Metab Syndr Obes. 2020;13:1695–701.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors acknowledge the staff and participants of the Malmö Diet and Cancer study.

Funding

Open access funding provided by Lund University. This study was funded by the Swedish Research Council (2020-01412), Heart and Lung Foundation (20190555), Crafoord Foundation (20201032), Håkansson Foundation, and Hjelt Foundation. Huiping Li holds a China Scholarship Council Studentship with the Lund University.

Author information

Authors and Affiliations

Authors

Contributions

HL, YB, and ES designed the study. HL, YB, and ES contributed to the analysis and interpretation of the data. HL wrote the manuscript. HL, YB, YW, and ES reviewed and edited the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Huiping Li or Emily Sonestedt.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Ethical Review Board at the Faculty of Medicine at Lund University, Sweden (LU 51/90), and was carried out in accordance with the Helsinki Declaration. All participants provided written informed consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table S1.

A total of 10 single nucleotide polymorphisms associated with AMY1 copy number. Table S2. Association between intake of AMY1-GRS and CVD and mortality risk (n = 20,264). Table S3. Association between starch intake and plasma proteins (n = 3680). Table S4. Association between plasma proteins and CVD and mortality risk (n = 3680). Table S5. Association between AMY1 copy number and plasma proteins (n = 3254). Figure S1. Flowchart of participant selection from the Malmö Diet and Cancer Study. Figure S2. Restricted cubic spline plots to assess association between AMY1 copy number and CVD and mortality. Figure S3. Restricted cubic spline plots to assess association between AMY1-GRS and CVD and mortality. Figure S4. Association between intake of starch and risk of CVD and mortality by quintiles of AMY1-GRS (n = 20,264).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, H., Borné, Y., Wang, Y. et al. Starch intake, amylase gene copy number variation, plasma proteins, and risk of cardiovascular disease and mortality. BMC Med 21, 27 (2023). https://doi.org/10.1186/s12916-022-02706-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-022-02706-5

Keywords