Abstract
Areca (Areca catechu L.) is one of the most important cash crops in China and is considered the fourth most widely used addictive substance. In addition, areca is widely used in traditional and herbal medicines. The major characteristics of the fruit are affected by its genetic background and growth environment. The growing environment in different regions will impact the quality of agricultural products and the processing quality. The quality of areca is not only the basis of its commercialization development and processing quality, but also is an important basis for the scientific planting of areca. Therefore, determining the quality of areca will provide evidence for scientific planting and more optimal applications. We evaluated the quality of areca by comparing the differences in physicochemical characteristics using principal component analysis (PCA) and hierarchical cluster analysis. A total of 165 arecas, in the same growth period, were collected from 11 main producing regions in Hainan Province. Our results illustrate that the physicochemical characteristics of areca in different regions were significantly different. The PCA was conducted using 10 quality indexes, and three principal components were extracted to reflect 80% of the original variables. The first principal component mainly reflected the fruit shape quality, the second principal component mainly reflected the hardness quality, and the third principal component mainly reflected the functional component quality. The relationship between each producing region and the principal component could be obtained intuitively from the principal component score plots. The arecas in Wanning and Wenchang were larger and their cellulose content was greater than in other areas, indicating that they were more suitable for processing. In contrast, the arecas in Baoting, Wuzhishan, Danzhou, Tunchang, and Dongfang had a greater arecoline content than the other areas, making them more suitable for use as medicinal materials. Hierarchical cluster analysis classified the 11 producing regions into five categories based on the measured parameters, which was consistent with the results of the PCA score plots. These results could provide information to improve the use of areca in China.
Areca (Areca catechu L.) is widely distributed throughout South and Southeast Asia, including China, India, Pakistan, Indonesia, Malaysia, the Philippines, and New Guinea (Heatubun et al. 2012; Verma 2011). In China, 95% of areca is planted in Hainan, of which 95% to 97% is processed into chewing products (Li et al. 2019; Sarode et al. 2013). During the past 5 years, Chinese areca production has increased steadily, with an annual output of 276,200 tons, and the annual output value of the areca industry has reached 10 billion US dollars. In the consumer market, the consumption of areca products has become the fourth largest worldwide, after cigarettes, alcohol, and coffee (Li and Zhang 2011). The major characteristics of the fruit are affected by its genetic background and growth environment. The growing environment in different regions will impact the quality of agricultural products and the processing quality (Liu et al. 2022; Wang et al. 2022). Therefore, the quality of areca is not only the basis of its commercialization development and processing quality, but also is an important basis for the scientific planting of areca. For example, the phenotypic traits of the areca fruit, such as weight, transverse diameter, and longitudinal diameter, can determine whether it is easy to eat. The cellulose, lignin, and pectin contents determine the hardness and feeling of chewing; hard fibers will cause oral mucosal damage. Alkaloids and polyphenols, which have refreshing effects, are the main active substances and medicinal ingredients in areca (Amudhan et al. 2012; Peng et al. 2015). These qualities affect sensory taste and consumer choice (Dai et al. 2021; Huang et al. 2021; Peng et al. 2021).
Limited studies have been conducted on the quality of areca. In this study, 165 arecas were collected during the same growth period from 11 main producing regions in Hainan Province. The quality was evaluated using PCA and hierarchical cluster analysis (HCA) and 10 indices, including weight, transverse diameter, longitudinal diameter, fiber layer thickness, cellulose, hemicellulose, pectin, and alkaloids. Consumers are paying increasing attention to the quality and safety of areca products. Therefore, this study provides information about areca quality to support its use, and improves the analysis of the impact of different regions on quality.
Materials and Methods
Plant material.
Arecas with a growth period of 150 d were selected from 11 producing regions in Hainan Province. Arecas that were green or dark green in coloration, free of any blemishes or deformities, and with a growth period of 150 d were selected. They included Wenchang (WC), Qionghai (QH), Wanning (WN), Baoting (BT), Sanya (SY), Dongfang (DF), Ledong (LD), Qiongzhong (QZ), Tunchang (TC), Danzhou (DZ), and Wuzhishan (WZS) (Table 1). Farmers planted the selected arecas. Fifteen arecas with appropriate maturity, no pests or diseases, and no mechanical damage were selected randomly from each region. Immediately after collection, the samples were transported to the laboratory for cold storage.
Sampling location.
Selection and treatment of samples.
Fifteen arecas were selected randomly from each region, washed with clean water, and dried. Each index was measured in triplicate. All chemical tests used 1 g of an areca-treated sample—meaning, after drying, crushing, and passing through a 0.25-mm mesh sieve.
Determination of weight, size, and fiber layer thickness.
A single areca was weighed using an analytical balance, and its width and length were measured using a Vernier caliper. The areca was cut down the middle, its core was removed, and the thickness of the fiber layer was measured using a Vernier caliper.
Determination of cellulose, hemicellulose, and lignin contents.
The content of each lignocellulose component was determined using a cellulose tester. Determination of the content of neutral detergent fiber and acid detergent fiber in the sample, the acid-washing lignin content of the sample, was determined using a 72% concentrated sulfuric acid hydrolysis method, and the ash content was determined after dry ashing at 550 °C for 2 h in a muffle furnace. Last, the cellulose, hemicellulose, and lignin contents of the samples were calculated.
Determination of pectin content.
Pectin content was measured using spectrophotometry according to the Chinese Agricultural Industry Standard (Ministry of Agriculture of China 2008). We used a spectrophotometer to measure the absorbance at a wavelength of 525 nm. We then calculated the pectin content using the standard curve.
Determination of arecoline.
The arecoline content was measured using high-performance liquid chromatography (HPLC; HITACHI Chromaster, Japan). The chromatographic parameters were as follows: XB-SCX strong cation-exchange column, 4.6 mm × 250 mm × 5 μm; mobile phase, acetonitrile, 0.2% phosphoric acid (35:65, pH adjusted to 3.8 with ammonia); detection wavelength, 215 nm; column temperature, 25 °C; injection volume, 20 μL; and flow rate, 1.0 mL⋅min–1.
Determination of tannin content.
The tannin content was measured using HPLC (HITACHI Chromaster). The chromatographic parameters were as follows: chromatographic column, C18 column of 4.6 mm × 150 mm × 5 μm; mobile phase, methanol, 0.1% phosphoric acid (25:75); detection wavelength, 275 nm; column temperature, 30 °C; injection volume 10 μL; and flow rate, 1.0 mL⋅min–1.
Statistical analysis.
Excel and the Statistical Package for the Social Sciences (SPSS) were used for statistical analysis. Single factor analysis of variance, least significant difference (LSD), and Waller–Duncan methods were used for the analysis of variance and multiple comparisons (α = 0.05); correlation analysis, PCA, and HCA were performed on the data of the measured samples. The original 2019b software (OriginLab, Northampton, MA, USA) was used for drawing, and all experiments were repeated in triplicate.
Results
Physical characteristics.
The weight, transverse diameter, and longitudinal diameter are the most important physical characteristics that determine the economic value of areca (Ren et al. 2020). As shown in Table 2, the weights, transverse diameters, and longitudinal diameters in the different regions were significantly different (P < 0.05). The weight, transverse diameter, and longitudinal diameter were the greatest in WN (28.40 g, 32.06 mm, and 55.40 mm, respectively), which were significantly greater than those in the other regions. The thickness of the fiber layer determines the chewing quality of the areca. A thick fiber layer increases the hardness of the areca, which damages the oral mucosa during chewing. Wanning had the smallest fiber layer thickness, with an average thickness of 5.36 mm, which was significantly less than that in the other regions (P < 0.05).
Physical characteristics of areca in different regions.
Chemical characteristics.
The chewing sensation is an important indicator of areca processing. Lignin determines the softness and hardness of chewing, whereas cellulose, hemicellulose, and pectin determine the elasticity and toughness of chewing (Bauer et al. 1973). Lignin is the main factor affecting the hardness of areca fibers (Fahlén and Salmén 2005; Jiang et al. 2018). During plant growth, lignin acts as a matrix to fill the remaining space in the cell wall, increasing the mechanical strength and enhancing the hardness of the plant (Li et al. 2019). Therefore, areca with a low lignin content should be selected for processing. Because of the hard texture of areca fibers, long-term chewing areca will lead to wear and tear and abrasion of the teeth, longitudinal fracture of teeth roots, and, in serious cases, it will damage the oral tissue, resulting in oral submucosal fibrous lesions (Arakeri and Brennan 2013; Brandt et al. 2013; Santos et al. 2020). Pectin plays important roles in cell wall hydration, adhesion, growth, extensibility, and elasticity (Cao et al. 2019; Ferreira-Lazarte et al. 2018; Li et al. 2021; Liao et al. 2021). Cellulose, hemicellulose, and pectin give areca products elasticity and toughness; therefore, fruit with high contents of cellulose, hemicellulose, and pectin should be used primarily when processing areca.
As shown in Table 3, the lignin, cellulose, and hemicellulose contents of areca in the different regions were significantly different (P < 0.05), but there was no significant difference in pectin content. In terms of lignin content, there was a decreasing trend from north to south, with the greatest lignin content in WC (115.80 mg⋅g–1) and the least in SY (102.75 mg⋅g–1). In terms of cellulose, there is a growth pattern from north to south, with the greatest cellulose content in SY (214.34 mg⋅g–1) and the least in QH (191.28 mg⋅g–1). Wanning had the greatest hemicellulose content (252.58 mg⋅g–1), which was significantly greater than those in the other regions.
Chemical characteristics of areca in different regions.
Alkaloids and polyphenols, which have antioxidative, antiviral, antitumor, antibacterial, and antidepressant effects, are the active components of areca (Papke et al. 2015; Yamada et al. 2015), and alkaloids are the main source of its stimulating effects (Bhat et al. 2017; Mehrtash et al. 2017). Therefore, areca, which has high alkaloid and polyphenol contents, is suitable for traditional and herbal medicines. As shown in Table 3, the areca alkaloid contents in the different regions were significantly different (P < 0.05); however, there was no significant difference in polyphenol content. The content was the greatest in the eastern region (5.42 mg⋅g–1), which was significantly greater than that in the other regions.
Principal component analysis.
Principal component analysis is a statistical analysis method that simplifies multiple indicators into a few comprehensive indicators and uses a few variables to reflect the information of the original variables as much as possible (Gao 2005). It can be seen from Table 4 that the eigenvalue of the first four principal components were greater than one. The cumulative variance contribution rates of the first three principal components were 80.00%, which integrated most of the area information, and three principal components were extracted.
Eigenvalues and cumulative variance contribution rates of autonomous components.
Table 5 lists the information contents of the different principal components. The first principal component (PC1) mainly synthesized the information on weight, longitudinal diameter, and fiber layer thickness, and the three indicators all presented a positive distribution on the first principal component. Therefore, PC1 mainly explained the fruit shape of the areca, which can be referred to as the fruit shape index, and the second principal component (PC2) mainly synthesized the information of cellulose and lignin. Lignin presented a negative distribution in PC2; cellulose presented a positive distribution. The second principal component mainly explained the hardness of areca, which can be referred to as the hardness index, and the third principal component (PC3) mainly synthesized information on alkaloids and polyphenols, in which alkaloids are distributed in a negative direction. Hence, PC3 was named the functional component index.
Component load score coefficient matrix.
In our study, PC1 and PC2 contained 40.91% and 25.97% of the original information, respectively. Researchers have used PCA score charts to reflect the relationships among indicators (Baraldi et al. 2007; Karlsen et al. 1999; Patras et al. 2011). The relationship between each region and PC1 vs. PC2 is shown in Fig. 1. Wanning and WC were located in the first quadrant, indicating arecas were larger and the cellulose content was greater. Qiongzhong, BT, and WZS, in the second quadrant, illustrated that arecas in these three regions were smaller and the cellulose content was greater. Six regions (LD, DF, SY, DZ, TC, and QH) were in the negative range of PC2, indicating the arecas in these six regions had a greater lignin content and were unsuitable for processing. Qiongzhong, BT, and WZS were in the negative range of PC1 and the positive range of PC2, explaining why the cellulose content of areca was greater and the fruit were smaller, and thus were more suitable for fresh consumption.
Principal component 3 comprised 13.1% of the original information. Figure 2 shows the relationship between each region and PC1 and PC3. Baoting, WZS, DZ, TC, and DF were in the negative range of PC3, indicating that the arecoline content in these four regions was greater and the arecas were more suitable for use as medicinal materials.
Hierarchical cluster analysis.
Based on physicochemical quality, the HCA of areca in 11 regions was determined. According to the HCA dendrogram (Fig. 3), samples from the 11 regions were divided into three clusters (distance between classes = 15).
The first cluster included WC and WN, which had larger fruit, higher cellulose content, and lower lignin content, and were suitable for processing. The second cluster contained BT, QZ, WZS, SY, DF, and LD. This group mainly included areca with smaller fruit, which are suitable for fresh consumption. Baoting, QZ, and WZS were smaller and had greater cellulose and arecoline contents, making them suitable for fresh consumption and medicine. The third cluster, including QH, TC, and DZ, had a higher lignin content and lower cellulose content, which is not recommended for fresh consumption and processing. In summary, the results of HCA and the principal component comprehensive evaluation were relatively consistent, indicating that both methods could be used to analyze the quality indicators of areca and to determine whether the arecas was suitable for fresh consumption, processing, or medicine.
Discussion
Principal component analysis and HCA are used increasingly by researchers in the quality analysis of multiple indicators of varieties (Geőcze et al. 2013; Keenan et al. 2012; Schnackenberg et al. 2010). In this study, 10 quality indices of areca from 11 producing regions were determined and investigated using statistics, PCA, and HCA. This can provide a foundation for areca use. For example, we can select regions in the third and fourth intervals of the PC1 and PC3 score charts for producing medicinal materials; fresh-consumption areca can be selected using the second and third intervals of the PC1 and PC2 score charts. Interestingly, the results showed that climatic conditions had an important effect on the quality of the arecas. Hainan’s climate is generally divided into three regions: the northern region, the central mountainous region, and the southern region (Chen et al. 2022). The PCA and HCA of the central mountainous region (BT, QZ, and WZS), southern region (SY, LD, and DF), and northern region (WC, WN, QH, TC, and DZ) reflected regional characteristics. In future studies, it will be necessary to consider whether the climatic characteristics of different regions lead to different accumulations of areca metabolites.
Conclusions
In this study, 10 quality indices were determined and analyzed for 165 arecas collected from 11 main areca-producing regions in Hainan. The results showed significant differences in eight quality indices of areca from different regions without pectin or polyphenols. Arecas in WN were the most suitable for processing because they had the largest fruit, the greatest hemicellulose content, and the thinnest fiber layer. The areca alkaloid content of DF was the greatest, making it suitable for medicinal use.
Principal component analysis was performed on the 10 quality indices, and three principal components were extracted to reflect 80% of the original variables. The first principal component mainly reflected the fruit shape quality, PC2 mainly reflected the hardness quality, and PC3 mainly reflected the functional component quality. The relationships among each production region and the principal component can be obtained intuitively from the principal component score plot. The arecas in WN and WC were larger, and the cellulose content was greater, making them more suitable for processing material. Baoting, WZS, DZ, TC, and DF had a greater arecoline content, which is more suitable for use as medicinal materials.
Based on the quality characteristics of arecas, the 11 producing regions were divided into three categories using HCA, and appropriate utilization fields were determined preliminarily. The results of the HCA and PCA were relatively consistent. The correlation between areca quality and the producing region was relatively significant. Principal component analysis and HCA can be used to analyze the quality index of areca and its appropriate use.
References Cited
Amudhan MS, Begum VH & Hebbar KB. 2012 A review on phytochemical and pharmacological potential of Areca catechu L. seed Int J Pharm Sci Res.3 11 4151 4157
Arakeri G & Brennan P. 2013 Dietary copper: A novel predisposing factor for oral submucous fibrosis? Med Hypotheses.80 3 241 243 https://doi.org/10.1016/j.mehy.2012.11.038
Baraldi C, Bodecchi LM, Cocchi M, Durante C, Ferrari G, Foca G, Grandi M, Marchetti A, Tassi L & Ulrici A. 2007 Chemical composition and characterisation of seeds from two varieties (pure and hybrid) of Aesculus hippocastanum Food Chem.104 1 229 236 https://doi.org/10.1016/j.foodchem.2006.11.032
Bauer WD, Talmadge KW, Keegstra K & Albersheim P. 1973 The structure of plant cell walls: II. The hemicellulose of the walls of suspension-cultured sycamore cells Plant Physiol.51 1 174 187 https://doi.org/10.1104/pp.51.1.174
Bhat SK, Ashwin D, Mythri S & Bhat S. 2017 Arecanut (Areca catechu L.) decreases Alzheimer’s disease symptoms: Compilation of research works J Med Plant Studies.5 5 4 9
Brandt A, Gräsvik J, Hallett JP & Welton T. 2013 Deconstruction of lignocellulosic biomass with ionic liquids Green Chem.15 3 550 583 https://doi.org/10.1039/c2gc36364j
Cao M-R, Yuan H-W, Daniyal M, Yu H-H, Xie Q-L, Liu Y-K, Li B, Jian Y-Q, Peng C-Y, Tan D-B, Peng Y-M, Choudhary M-L, Rahman A & Wang W. 2019 Two new alkaloids isolated from traditional Chinese medicine Binglang the fruit of Areca catechu Fitoterapia.138 104276 https://doi.org/10.1016/j
Chen Y-X, Ren L-X, Lou Y-S, Tang L-L, Yang J-Z & Su L. 2022 Climate suitability analysis of green orange cultivation in Hainan island under future climate change scenarios Zhonghua Nongye Qixiang.43 10 786 797 https://doi.org/10.3969/j.issn.1000-6362.2022.10.002
Dai W-T, Ji J-B, Kang X-N, Ma J-S, Dai J-H & Wang S-P. 2021 Research progress of areca drying technology Food Ind.42 3 258 262
Dammström S, Salmén L & Gatenholm P. 2008 On the interactions between cellulose and xylan, a biomimetic simulation of the hardwood cell wall BioResources.4 1 3 14 https://doi.org/10.1515/HF.2009.036
Fahlén J & Salmén L. 2005 Pore and matrix distribution in the fiber wall revealed by atomic force microscopy and image analysis Biomacromolecules.6 1 433 438 https://doi.org/10.1021/bm040068x
Ferreira-Lazarte A, Kachrimanidou V, Villamiel M, Rastall RA & Moreno F-J. 2018 In vitro fermentation properties of pectins and enzymatic-modified pectins obtained from different renewable bioresources Carbohydr Polym.199 482 491 https://doi.org/10.1016/j.carbpol.2018.07.041
Gao H-X. 2005 Applied multivariate statistical analysis Beijing University Press Beijing, China
Geőcze KC, Barbosa LCA, Fidêncio PH, Silvério F-O, Lima C-F, Barbosa M-C-A & Ismail F-M-D. 2013 Essential oils from pequi fruits from the Brazilian Cerrado ecosystem Food Res Int.54 1 1 8 https://doi.org/10.1016/j.foodres.2013.06.005
Heatubun CD, Dransfield J, Flynn T, Tjitrosoedirdjo SS, Mogea JP & Baker WJ. 2012 A monograph of the betel nut palms (Areca: Arecaceae) of East Malesia Bot J Linn Soc.168 2 147 173 https://doi.org/10.1111/j.1095-8339.2011.01199.x
Huang Q, Liu K-K & Chen K. 2021 Effect of processing technology on the content of arecoline in areca product Food Ind.42 9 57 60
Jiang Z, Zhao P & Hu C. 2018 Controlling the cleavage of the inter- and intra-molecular linkages in lignocellulosic biomass for further biorefining: A review Bioresour Technol.256 466 477 https://doi.org/10.1016/j.biortech.2018.02.061
Karlsen AM, Aaby K, Sivertsen H, Baardseth P & Ellekjær M-R. 1999 Instrumental and sensory analysis of fresh Norwegian and imported apples Food Qual Prefer.10 4–5 305 314 https://doi.org/10.1016/S0950-3293(99)00030-0
Keenan DF, Valverde J, Gormley R, Butler F & Brunton NP. 2012 Selecting apple cultivars for use in ready-to-eat desserts based on multivariate analyses of physico-chemical properties Lebensm Wiss Technol.48 2 308 315 https://doi.org/10.1016/j.lwt.2012.04.005
Li Y-C, Cheng A-J, Lee L-Y, Huang Y-C & Chang J-T-C. 2019 Multifaceted mechanisms of areca nuts in oral carcinogenesis: The molecular pathology from precancerous condition to malignant transformation J Cancer.10 17 4054 https://doi.org/10.7150/jca.29765
Li X, Wei Z, Wang X, Duan F-X, Xiong L-D, Li J-W, Tian J, Jia L-R & Gao H. 2021 Premna microphylla Turcz leaf pectin exhibited antioxidant and anti-inflammatory activities in LPS-stimulated RAW 264.7 macrophages Food Chem.349 129164 https://doi.org/10.1016/j.foodchem.2021.129164
Li Z & Zhang H-D. 2011 Present situation and prospect of areca processing Farm Products Processing.8 3 5
Liao T, Zhou L, Liu J-P, Zhou L-Q, Liu C-M & Liu W. 2021 Effects of pectin on the activity and thermal stability of polyphenoloxidase J Chin Inst Food Sci Technol.21 5 67 74 https://doi.org/10.16429/j.1009-7848.2021.05.008
Liu Y-X, Liao Y-X, Guo M-R, Zhang W-D, Sang Y-Y, Wang H, Cheng S-B & Chen G-G. 2022 Comparative elucidation of bioactive and volatile components in dry mature jujube fruit (Ziziphus jujuba Mill.) subjected to different drying methods Food Chem X.14 100311 https://doi.org/10.1016/j.fochx.2022.100311
Mehrtash H, Duncan K, Parascandola M, David A, Gritz ER, Gupta PC, Mehrotra R, Nordin ASA, Pearlman PC, Warnakulasuriya S, Wen CP, Zain RB & Trimble EL. 2017 Defining a global research and policy agenda for betel quid and areca nut Lancet Oncol.18 12 e767 e775 https://doi.org/10.1016/S1470-2045(17)30460-6
Ministry of Agriculture of China 2008 Determination of pectin content in fruits and derived product-spectrophotometry method, NY/T2016-2011 Ministry of Agriculture of China Beijing, China
Papke R-L, Horenstein NA & Stokes C. 2015 Nicotinic activity of arecoline, the psychoactive element of “betel nuts”, suggests a basis for habitual use and anti-inflammatory activity PLoS One.10 10 e0140907 https://doi.org/10.1371/journal.pone.0140907
Patras A, Brunton NP, Downey G, Rawson A, Warriner K & Gernigon G. 2011 Application of principal component and hierarchical cluster analysis to classify fruits and vegetables commonly consumed in Ireland based on in vitro antioxidant activity J Food Compos Anal.24 2 250 256 https://doi.org/10.1016/j.jfca.2010.09.012
Peng Y-M, Gong N-C & Li P-H. 2021 Changes and multivariate statistical analysis of components in edible areca during processing China Food Ind.21 96 101
Peng W, Liu Y-J, Wu N, Sun T, He X-Y, Gao Y-X & Wu C-J. 2015 Areca catechu L. (Arecaceae): A review of its traditional uses, botany, phytochemistry, pharmacology and toxicology J Ethnopharmacol.164 340 356 https://doi.org/10.1016/j.jep.2015.02.010
Ren J-L, Liu H-Y & Wang X-H. 2020 Various approaches of lignocellulose utilization and strategies for hemicellulose-prepared lignocellulose biorefinery Shengwu Jiagong Guocheng.18 1 1 12 https://doi.org/10.3969/j.issn.1672-3678.2020.01.001
Santos EE, Amaro RC, Bustamante CCC, Guerra MHA, Soares LC & Froes RES. 2020 Extraction of pectin from agroindustrial residue with an ecofriendly solvent: Use of FTIR and chemometrics to differentiate pectins according to degree of methyl esterification Food Hydrocoll.107 105921 https://doi.org/10.1016/j.foodhyd.2020.105921
Sarode SC, Mahuli A, Sarode GS & Mahuli S. 2013 Why only areca nut chewing cannot cause oral submucous fibrosis? Med Hypotheses.81 1 47 49 https://doi.org/10.1016/j.mehy.2013.02.025
Schnackenberg BJ, Saini UT, Robinson BL, Ali SF & Patterson TA. 2010 An acute dose of gamma-hydroxybutyric acid alters gene expression in multiple mouse brain regions Neuroscience.170 2 523 541 https://doi.org/10.1016/j.neuroscience.2010.06.049
Verma S. 2011 Areca nut (betel nut) chewing: A popular Indian cultural practice and its mucosal implications Int J Dermatol.50 2 229 232 https://doi.org/10.1111/j.1365-4632.2010.04710.x
Wang S-H, Luo L, Zhang R, Guo K, Zhang X-B, Kong L-J & Zhao Z-G. 2022 Biochemical composition and quality of wild-caught adult mitten crabs from three river basins J Food Compos Anal.110 104574 https://doi.org/10.1016/j.jfca.2022.104574
Yamada T, Hara K & Kadowaki T. 2015 Betel chewing and risk of metabolic disease, cardiovascular disease, and all-cause mortality: A meta-analysis Int J Epidemiol.44 184 https://doi.org/10.1093/ije/dyv096.278