Breast cancer is one of the most common cancers affecting women worldwide, with an age-standardised incidence rate of 47·8 cases per 100 000 women(Reference Sung, Ferlay and Siegel1). Based on the estimate from the WHO, 416 371 women were newly diagnosed with breast cancer in China in 2020(2). Several lifestyle factors, including a healthy diet, impact breast cancer risk(Reference Coughlin3,Reference Kruk4) . Over the past few decades, interest has increased in the associations between fatty acids and carcinogenesis, particularly n-3 PUFA. Accumulated animal and in vitro experimental evidence has indicated that n-3 PUFA inhibit the pathogenesis and progression of breast cancer(Reference Witte and Hardman5,Reference Fabian, Kimler and Hursting6) .
Erythrocyte membrane n-3 PUFA as a biomarker of n-3 PUFA status(Reference Ferreri, Sansone and Ferreri7,Reference Harris and Thomas8) reflects dietary intake over several months and represents an integrative measure of the interaction concerning dietary, metabolic and genetic factors(Reference Arab9). Despite growing experimental evidence supporting the anti-tumourigenic effects of n-3 PUFA, the epidemiological evidence on the potential protective role of n-3 PUFA in breast cancer risk reduction has been inconclusive. Thus far, only two case–control studies(Reference Kuriki, Hirose and Wakai10,Reference Shannon, King and Moshofsky11) and three nested case–control studies(Reference Pala, Krogh and Muti12–Reference Hirko, Chai and Spiegelman14) have evaluated erythrocyte membrane n-3 PUFA in relation to the risk of breast cancer and the results remain inconsistent. Inverse associations of erythrocyte membrane n-3 PUFA with the odds of breast cancer were observed in a case–control study from Japan(Reference Kuriki, Hirose and Wakai10), a case–control study from China(Reference Shannon, King and Moshofsky11) and a nested case–control study from Italy(Reference Pala, Krogh and Muti12). Other nested case–control studies from the USA and Europe found no significant relationships(Reference Wirfält, Vessby and Mattisson13,Reference Hirko, Chai and Spiegelman14) . Considering the vast differences in the dietary habits and intake pattern between the Chinese and Western populations, and the changes over time, the relationship between n-3 PUFA and risk of breast cancer in Chinese women still demands investigation.
We aimed to investigate the associations of individual and total erythrocyte membrane n-3 PUFA with the odds of being a breast cancer patient among Chinese women. We hypothesised that higher levels of erythrocyte membrane total n-3 PUFA and individual n-3 PUFA were inversely associated with odds of breast cancer.
Methods
Study population
The detailed study methods for this hospital-based case–control study have been published elsewhere(Reference Feng, Zhan and Zuo15). Briefly, eligible breast cancer cases were recruited from Sun Yat-sen University First Affiliated Hospital and Guangdong Women and Children Hospital from September 2011 to December 2019. The inclusion criteria of the breast cancer cases were as follows: (1) aged between 25 and 70 years, (2) histologically confirmed breast cancer and diagnosed within 3 months before the interview and (3) Guangdong natives or residing in Guangdong for at least 5 years. Participants who could not understand or speak Mandarin or Cantonese were excluded. Overall, 1675 of 1882 eligible cases were recruited, with a participation rate of 89·00 %. Among them, 869 provided blood samples. Finally, 853 cases were included in the analysis after excluding sixteen participants who had insufficient blood samples for analysis.
Control subjects were recruited from other departments of the same hospitals during the same period. The same inclusion criteria were applied for the control subjects, and they must have no history of cancer. Finally, we identified 1932 eligible control subjects, of whom 1749 were successfully interviewed (with a participation rate of 90·53 %) and 1169 provided blood samples. Among them, 892 control subjects were frequency matched to the cases by 5-year age intervals and completed the laboratory analysis.
The study was approved by the Ethical Boards of the School of Public Health, Sun Yat-sen University (approval number: 2011–18) and was conducted according to the Declaration of Helsinki. Written informed consent was obtained from all participants.
Data collection
A face-to-face interview was conducted by trained interviewers to collect information on potential confounding factors, including demographic characteristics, prior history of disease, menstrual and reproductive history, anthropometric measurements and lifestyle factors. A FFQ with eighty-one food items(Reference Zhang and Ho16) was used to collect information on dietary intake. Total energy intake was assessed according to the China Food Composition Table(Reference Yang17). The BMI was calculated as the current body weight (kg) divided by the height squared (m2). The metabolic equivalent (MET) hours per day was used to estimate physical activity, and the detailed methods of calculating MET have been described previously(Reference Ainsworth, Haskell and Herrmann18).
Information on the oestrogen receptor (ER) status was available for 803 cases (94·14 %); among these, 581 cases were ER+ and 222 cases were ER−. Information on the progesterone receptor (PR) status was available for 800 cases (93·79 %); among these, 494 cases were PR+ and 306 cases were PR−.
Laboratory analysis of erythrocyte membrane n-3 PUFA
Venous blood samples were obtained following an overnight fasting period (8–12 h) on the second day of the participants’ admission and collected in EDTA tubes. The breast cancer cases and control subjects did not receive medication or surgery before blood collection. Erythrocytes were washed and separated within 2 h of collection and were stored at −80°C for subsequent analysis. The fatty acid compositions in the erythrocytes were measured by GC(Reference Ding, Li and Xiao19). The extraction of fatty acids was conducted with 100 μl erythrocyte sample and 750 μl chloroform/methanol (2:1, v/v)(Reference Folch, Lee and Sloane-Stanley20) and transmethylation with 750 μl 14 % boron-trifluoride ether/methanol (1:3, v/v) solution(Reference Lepage and Roy21) for 60 min at 90°C. The fatty acid methyl esters were analysed using an Agilent 7890A GC system (Agilent) equipped with a DB-23 capillary column (60 m × 0·25 mm internal diameter × 0·15 μm film; Agilent) and a flame ionisation detector. Considering nitrogen as carrier gas, the GC method included a split ratio of 5:1 with the injection temperature at 250°C. The oven temperature started at 50°C for 1 min and was programmed from 50 to 185°C at a rate of 25°C/min followed by a 5-min hold period, and the temperature was continuously increased to 230°C at a rate of 2·4°C/min followed by a 20-min hold period. Comparing the retention time of the samples with commercially available standards to identify individual fatty acids, the amount of each fatty acid was expressed as a percentage of the total erythrocyte membrane fatty acids (relative, %). The following n-3 PUFA was included in the analysis: α-linolenic acid (ALA), EPA, docosapentaenoic acid (DPA), DHA and total n-3 PUFA. Intra-assay and inter-assay CV were determined using replicate assays of samples on the same day (n 12–16) and different days (n 100), respectively. The intra-assay CV were 5·37 %, 6·42 %, 6·05 % and 4·56 % for ALA, EPA, DPA and DHA, respectively. The inter-assay CV were < 20 % for all measured fatty acids.
Statistical analysis
Differences in characteristics between the case and control groups were examined using the Wilcoxon rank-sum test for continuous variables with skewed distribution and χ 2 test for categorical variables. Logistic regression analysis was used to examine the association between erythrocyte membrane n-3 PUFA and odds of breast cancer, and then the OR and 95 % CI were estimated. We classified the participants’ erythrocyte membrane n-3 PUFA into quartiles (Q) based on the proportion of n-3 PUFA in control subjects and assumed that the median values of each quartile were continuous variables to perform the P-value for trend test.
Multivariable regression models were applied to adjust for confounding factors, and covariates were selected based on preliminary analyses and the literature(Reference McPherson, Steel and Dixon22,Reference Turner23) , such as age, educational level, BMI, passive smoking, alcohol drinking, physical activity, age at menarche, history of breast cancer in first-degree relatives, history of benign breast disease and total energy intake. The level of erythrocyte membrane total n-6 PUFA was further adjusted as a covariate. In addition, we performed restricted cubic spline(Reference Desquilbet and Mariotti24) with knots placed at the 5th, 35th, 65th and 95th percentiles to reveal the potential non-linear relationship between erythrocyte membrane n-3 PUFA and odds of breast cancer.
We conducted stratified analysis by menopausal status (pre- v. postmenopausal women). Stratified analysis by BMI (< 24 kg/m2 v. ≥ 24 kg/m2) was also performed to examine the consistency of overall associations. Additionally, we considered the potential heterogeneity of tumour subtypes, defined by hormone receptors (ER + v. ER−; PR + v. PR−). The Wald test was used to test for heterogeneity, comparing the coefficients across cancer subtypes.
All the analyses were performed using SPSS version 22.0 (SPSS; IBM) and STATA version 15.0 (StataCorp). A two-sided P value < 0·05 was considered statistically significant.
Results
Participant characteristics
A detailed comparison of characteristics between the case and control subjects is shown in Table 1. Compared with the control subjects, the cases had higher BMI values, lower energy intake, lower educational level and participated in lower levels of physical activity. Additionally, the cases were more likely to drink alcohol, to be passive smokers, to be younger at menarche, to have a history of breast cancer in first-degree relatives and benign breast disease than the control subjects. The levels of erythrocyte membrane ALA, EPA, DPA, DHA and total n-3 PUFA were significantly lower in cases than in control subjects.
ALA, α-linolenic acid; DPA, docosapentaenoic acid.
* P value for between-group differences between cases and controls is based on Wilcoxon rank-sum test for continuous variables and χ 2 test for categorical variables.
† Among women with children.
Associations between erythrocyte membrane n-3 PUFA and odds of breast cancer
As shown in Table 2, after adjustment for potential confounding factors, individual and total erythrocyte membrane n-3 PUFA were inversely associated with the odds of breast cancer (all P trend < 0·001). Compared with the lowest quartile (Q1) of erythrocyte membrane ALA, the highest quartile (Q4) had a 43 % reduction in the odds of breast cancer. Likewise, compared with Q1 of erythrocyte membrane EPA, DPA, DHA and total n-3 PUFA, Q4 had, respectively, a 40 %, 51 %, 47 % and 61 % reduction in the odds of breast cancer. These associations were unaltered after further adjusting for erythrocyte membrane total n-6 PUFA.
ALA, α-linolenic acid; DPA, docosapentaenoic acid.
* Adjusted for age, educational level, BMI, passive smoking, alcohol drinking, physical activity, age at menarche, history of breast cancer in first-degree relatives, history of benign breast disease and total energy intake.
† Additionally adjusted for erythrocyte membrane total n-6 PUFA.
Dose–response relationships between erythrocyte membrane n-3 PUFA and odds of breast cancer
Figure 1 demonstrates the non-linear association between erythrocyte membrane n-3 PUFA and odds of breast cancer using restricted cubic spline. Non-linear associations of ALA, DPA and total n-3 PUFA with odds of breast cancer were observed (P non-linear < 0·05). Overall, odds of having breast cancer decreased with increased erythrocyte membrane EPA and DHA (P non-linear > 0·05). Specifically, the odds of breast cancer decreased rapidly as erythrocyte membrane ALA increased to 0·62 % and increased slowly above 0·62 %. However, the increased odds of breast cancer were not significant and may be driven by a wide 95 % CI when erythrocyte membrane ALA above 0·91 %. We also noted a similar pattern between odds of breast cancer and erythrocyte membrane DPA, with the cut-off point of DPA of 1·82 %.
Associations of erythrocyte membrane n-3 PUFA with odds of breast cancer stratified by menopausal status
Stratification by menopausal status showed that the associations of erythrocyte membrane total and individual n-3 PUFA with odds of breast cancer were similar to the overall association, except for erythrocyte membrane ALA. The decreased odds of breast cancer in association with increasing erythrocyte membrane ALA appeared in postmenopausal women (ORQ4 v. Q1 = 0·30; 95 % CI 0·18, 0·50; P trend < 0·001) but not in premenopausal women (ORQ4 v. Q1 = 0·80; 95 % CI 0·56, 1·13; P trend = 0·094; P interaction = 0·015; see Table 3).
ALA, α-linolenic acid; DPA, docosapentaenoic acid.
* Adjusted for age, educational level, BMI, alcohol drinking, passive smoking, physical activity, age at menarche, history of breast cancer in first-degree relatives, history of benign breast disease, total energy intake and erythrocyte membrane total n-6 PUFA.
Associations of erythrocyte membrane n-3 PUFA with odds of breast cancer stratified by BMI
Stratified analysis according to BMI showed that the levels of erythrocyte membrane ALA, EPA, DPA, DHA and total n-3 PUFA were consistently inversely associated with odds of breast cancer in BMI-stratified models (all P interaction > 0·05; see Fig. 2).
Associations of erythrocyte membrane n-3 PUFA with odds of breast cancer stratified by oestrogen receptor and progesterone receptor status
Subgroup analysis by hormone receptor status (ER and PR) showed consistent associations with the overall association. Notably, women with higher proportions of erythrocyte membrane DHA had a lower likelihood of being ER+ breast cancer patient than ER− breast cancer (ER+: ORQ4 v. Q1 = 0·46; 95 % CI 0·33, 0·63; ER−: ORQ4 v. Q1 = 0·68; 95 % CI 0·44, 1·06; see Table 4), but the P value for the test of heterogeneity was not statistically significant (P heterogeneity = 0·167).
ALA, α-linolenic acid; DPA, docosapentaenoic acid; ER, oestrogen receptor; PR, progesterone receptor.
* Adjusted for age, educational level, BMI, passive smoking, alcohol drinking, physical activity, age at menarche, history of breast cancer in first-degree relatives, history of benign breast disease, total energy intake and erythrocyte membrane total n-6 PUFA.
† P heterogeneity for ER + v. ER−.
‡ P heterogeneity for PR + v. PR−.
Discussion
This study aimed to evaluate the relationships between erythrocyte membrane n-3 PUFA and the odds of being a breast cancer patient. Higher erythrocyte membrane individual and total n-3 PUFA were inversely associated with the odds of breast cancer. Additionally, our findings suggested significant non-linear relationships between erythrocyte membrane ALA, DPA, total n-3 PUFA and odds of breast cancer.
Most previous studies used plasma(Reference Chajès, Assi and Biessy25–Reference Matta, Deubler and Chajes29) or serum(Reference Lope, Guerrero-Zotano and Casas30–Reference Saadatian-Elahi, Toniolo and Ferrari32) n-3 PUFA as biomarkers to evaluate the association of n-3 PUFA with breast cancer risk. Although one study from plasma(Reference Newell, Ghosh and Goruk28) supports our findings, most studies using plasma(Reference Chajès, Assi and Biessy25–Reference Pouchieu, Chajès and Laporte27,Reference Matta, Deubler and Chajes29) or serum(Reference Lope, Guerrero-Zotano and Casas30–Reference Saadatian-Elahi, Toniolo and Ferrari32) n-3 PUFA as biomarkers reported null findings on the relationship between n-3 PUFA and breast cancer risk. Compared with the plasma or serum fatty acid, the erythrocyte membrane fatty acid profile shows a longer time for fatty acid intake, approximately 4 months, and is not susceptible to the individual postprandial status(Reference Harris and Thomas8,Reference Harris33) . Our observation of the inverse associations between erythrocyte membrane total and individual n-3 PUFA and breast cancer risk agreed with some previous study findings(Reference Kuriki, Hirose and Wakai10–Reference Pala, Krogh and Muti12). A case–control study of Japanese women showed that erythrocyte membrane EPA, DPA, DHA and total n-3 PUFA were inversely associated with odds of breast cancer(Reference Kuriki, Hirose and Wakai10). Another case–control study from Shanghai, China noted that erythrocyte membrane EPA and total n-3 PUFA were inversely associated with odds of breast cancer(Reference Shannon, King and Moshofsky11). By contrast, erythrocyte membrane total n-3 PUFA, EPA, DPA and DHA were not associated with odds of having breast cancer in two nested case–control studies within the Malm Diet Cancer cohort (Sweden)(Reference Wirfält, Vessby and Mattisson13) and within the Nurses’ Health Study II(Reference Hirko, Chai and Spiegelman14). One study in Italy only observed a marginal inverse relationship between erythrocyte membrane DHA and odds of breast cancer(Reference Pala, Krogh and Muti12). To our best knowledge, no study has reported that the erythrocyte membrane ALA is significantly associated with odds of breast cancer(Reference Kuriki, Hirose and Wakai10–Reference Hirko, Chai and Spiegelman14). Our observation of an inverse association between the erythrocyte membrane ALA and odds of having breast cancer must be confirmed in further studies.
The conflicting results from different studies may be due to the variation in the contents of erythrocyte membrane n-3 PUFA (total and subclasses). Of the four studies mentioned above(Reference Shannon, King and Moshofsky11–Reference Hirko, Chai and Spiegelman14), the levels of erythrocyte membrane total n-3 PUFA were higher in the Shanghai(Reference Shannon, King and Moshofsky11), Italian(Reference Pala, Krogh and Muti12) and Swedish(Reference Wirfält, Vessby and Mattisson13) studies than in ours but were lower in the US study(Reference Hirko, Chai and Spiegelman14). Regarding the individual levels of erythrocyte membrane n-3 PUFA, the EPA level was lower in three studies(Reference Shannon, King and Moshofsky11,Reference Pala, Krogh and Muti12,Reference Hirko, Chai and Spiegelman14) than in ours except for the Swedish study(Reference Wirfält, Vessby and Mattisson13). The content of erythrocyte membrane DHA was higher in three studies(Reference Shannon, King and Moshofsky11–Reference Wirfält, Vessby and Mattisson13) than in ours and except the US study(Reference Hirko, Chai and Spiegelman14). Concerning the erythrocyte membrane DPA level, all the studies had higher levels than ours(Reference Shannon, King and Moshofsky11–Reference Hirko, Chai and Spiegelman14). However, the Japanese study(Reference Kuriki, Hirose and Wakai10) revealed a similar median proportion of erythrocyte membrane EPA, DPA, DHA and total n-3 PUFA to ours, likely explaining the consistent results. Additionally, the relatively higher level of erythrocyte membrane ALA in our study compared with that in other studies may help explain our observed reduction of breast cancer with erythrocyte membrane ALA.
Our study demonstrated non-linear associations of the erythrocyte membrane ALA, DPA and total n-3 PUFA and odds of breast cancer. The findings revealed that the erythrocyte membrane ALA level of 0·48–0·91 % may provide the maximum benefit for odds of breast cancer. We also showed a similar pattern of erythrocyte membrane DPA at the level from 0·86 % to 4·85 %. Our observation of non-linear associations of erythrocyte membrane ALA, DPA and DHA with odds of breast cancer might also help to explain the inconsistent results across studies. Further investigation is warranted regarding the optimal level of erythrocyte membrane n-3 PUFA in different populations and the validation of cut-off points.
Possible mechanisms explaining the apparent benefit of n-3 PUFA on breast cancer risk are as follows. n-3 PUFA is precursor to a range of biologically active metabolites, such as resolvins(Reference Serhan, Hong and Gronert34) and n-3 PUFA-derived eicosanoids(Reference Calder35) that act as anti-inflammatory agents to alter cell signalling(Reference Al-Jawadi, Moussa and Ramalingam36). Animal studies have documented that n-3 PUFA modulates the tumour microenvironment, resulting in a decrease of neutrophils and macrophages infiltration and an increase of CD3+ lymphocytes infiltration, tumour cell apoptosis and IL-10 expression(Reference Khadge, Thiele and Sharp37,Reference Garay, Comba and Vara Messler38) . Additionally, n-3 PUFA reduce the production of pro-inflammatory PGE2 metabolites by inhibiting arachidonic acid synthesis, reducing breast cancer risk(Reference Liu and Ma39). Overall, laboratory evidence demonstrates that the role of n-3 PUFA in inhibiting breast cancer occurs at all stages of cancer(Reference Witte and Hardman5).
Breast cancer is a heterogeneous hormone-related disease with varying aetiologies depending on the menopausal status(Reference Harbeck, Penault-Llorca and Cortes40,Reference Parsa and Parsa41) . Our analysis showed an inverse association of erythrocyte membrane ALA with breast cancer in postmenopausal women but not in premenopausal women. Consistent with our results, a Canadian cohort showed a decreased risk for breast cancer in postmenopausal women with higher concentrations of plasma ALA and DHA(Reference Newell, Ghosh and Goruk28). This Canadian study and another meta-analysis(Reference Zheng, Hu and Zhao42) were concordant in observing a trend of a reduced risk of breast cancer associated with total n-3 PUFA in postmenopausal women, driven by a reduced risk associated with ALA. However, two prior studies conducted in postmenopausal women had inconsistent results when examining the relationships between erythrocyte membrane n-3 PUFA and odds of breast cancer(Reference Pala, Krogh and Muti12,Reference Wirfält, Vessby and Mattisson13) . One nested case–control study in Italy illustrated a borderline negative association between the erythrocyte membrane DHA and odds of breast cancer(Reference Pala, Krogh and Muti12). Another nested case–control study in Sweden found no associations(Reference Wirfält, Vessby and Mattisson13). Further studies are needed to explore the associations between erythrocyte membrane n-3 PUFA and odds of breast cancer by menopausal status.
Likewise, breast cancer as a heterogeneous disease is classified into different subtypes according to the expression of ER and PR(Reference Perou, Sørlie and Eisen43). Genes and pathways differentially express on ER− and ER+ cells in response to EPA and arachidonic acid by microarrays(Reference Alquobaili, Miller and Muhie44). However, no significant heterogeneity existed in the ER/PR status between erythrocyte membrane n-3 PUFA and odds of breast cancer, a finding that was consistent with other study findings concerning the relationships of plasma(Reference Chajès, Assi and Biessy25), serum(Reference Lope, Guerrero-Zotano and Casas30) and erythrocyte membrane(Reference Hirko, Chai and Spiegelman14) n-3 PUFA with breast cancer risk. Given the small number of ER− and PR− patients in the available studies, suggestive results must be confirmed with larger case numbers.
We did not find that the obesity status significantly modified the inverse associations of erythrocyte membrane individual and total n-3 PUFA with odds of breast cancer. However, the inverse associations between erythrocyte membrane n-3 PUFA and odds of breast cancer, except for erythrocyte membrane DHA, were stronger in overweight/obese women than in women with BMI < 24 kg/m2. Similarly, a nested case–control study(Reference Hirko, Chai and Spiegelman14) observed significant inverse associations of erythrocyte membrane total n-3 PUFA, ALA, EPA and DPA in overweight/obese women, while no associations were observed among women with BMI < 25 kg/m2. Obesity is a critical risk factor for breast cancer(Reference Suzuki, Kojima and Tokudome45). Inflammation has been proposed as a potential link between obesity and breast cancer(Reference Roberts, Dive and Renehan46), and n-3 PUFA have anti-inflammatory effects. Thus, more investigations are warranted to carefully assess how the BMI may modulate the association of total n-3 PUFA and n-3 PUFA subclasses with breast cancer.
Our study has several strengths. The level of erythrocyte membrane n-3 PUFA is considered a good biomarker to reflect a relatively long-term intake of dietary n-3 PUFA(Reference Sun, Ma and Campos47). Additionally, a relatively large number of study subjects allowed us to explore the relationship between erythrocyte membrane n-3 PUFA and odds of breast cancer stratified by BMI or menopausal status.
This study has some limitations that must be noted. First, the cases and control subjects were recruited from large tertiary hospitals; thus, selection bias may occur. To minimise this bias, we recruited cases from two major hospitals consecutively as well as control subjects with several conditions that were not related to either dietary causes or breast cancer. Additionally, the level of erythrocyte membrane total n-3 PUFA in the control subjects in our study was similar to that in a previous study in a community-dwelling Chinese population (participants without incident metabolic syndrome: 7·09 (sd 1·93) %)(Reference Chen, Wu and Zhu48). Second, the level of n-3 PUFA in the blood may change over time during blood banking/storage. However, one study showed that no statistically significant increase was observed in the level of n-3 PUFA over time in the Cardiovascular Health Study across 13 consecutive years of measures(Reference Lai, de Oliveira Otto and Lemaitre49). Furthermore, the half-life of erythrocytes is 4 months. Therefore, a single pre-operative blood collection can still represent the long-term status of n-3 PUFA.
In summary, our findings showed that total and individual erythrocyte membrane n-3 PUFA were inversely associated with odds of having breast cancer. Other factors, such as menopause and hormone receptor status, may warrant further investigation when examining the association between n-3 PUFA and breast cancer. However, due to the lack of measuring other biomarkers of dietary intake, increased n-3 PUFA levels could be related more to healthier dietary habits in general. Therefore, increased n-3 PUFA may suggest an overall healthier eating which might result in a reduced risk of breast cancer. It is needed to measure more biomarkers of food intake in future studies and to examine their relationship with breast cancer.
Acknowledgements
We gratefully acknowledge the contribution of the study participants; without them, the study would not have been finished.
This work was supported by the National Natural Science Foundation of China (No: 81973020, No. 81102188). The funders had no role in the design of the study, analysis of the data or writing of this manuscript.
Conceptualisation, C. X. Z. and Z. L. Z.; data curation, C. X. Z. and Z. L. Z.; formal analysis, Z. L. Z. and D. D. S.; funding acquisition, C. X. Z.; investigation, Z. L. Z., D. D. S., L. X. and Q. X. W.; methodology, C. X. Z. and Z. L. Z.; project administration, C. X. Z.; resources, X. X. Z. and C. X. Z.; supervision, C. X. Z. and S. C. H.; visualisation, C. X. Z. and Z. L. Z.; writing – original draft, Z. L. Z.; writing – review and editing, C. X. Z. and S. C. H.
The authors declare that they have no competing interests.