- 1College of Light Industry and Food Engineering, Guangxi University, Nanning, China
- 2Medical College, Guangxi University, Nanning, China
- 3Department of Pharmacology and Nutritional Sciences, University of Kentucky, Lexington, KY, United States
Introduction: Recent evidence supports a role for the gut microbe-metabolites in longevity. However, the phenomenon of hypertension is more common in the longevity area and whether hypertension is associated with longevity remains unclear. Here, we hypothesize that the levels of gut microbiota, SCFAs, and urine metabolites were different between hypertension elderly and hypertension longevity.
Methods: We recruited 46 elderly volunteers from Donglan County, Guangxi, and 32 were selected and included in the experiment. The subjects with hypertension were divided into two groups according to age, Hypertension Elderly (HTE, aged 70.5 ± 8.59, n = 19) and Hypertension Longevity (HTL, aged 100 ± 5.72, n = 13). The gut microbiota, SCFAs, and urine metabolites were determined by three-generation 16S rRNA full-length sequencing, GC-MS, and 1H-NMR, respectively.
Results: Compared with the HTL group, the HTE group had higher levels of hypertension-related genera Klebsiella and Streptococcus, while having lower levels of the SCFA-producing genera Bacteroides, Faecalibacterium, and Alistipes. Based on LEFse analysis, Klebsiella pneumoniae, Lactobacillus gasseri, Streptococcus salivarius, Ruminococcus, Actinomyces, Rikenellaceae, f_Saccharimonadaceae, Clostridium perfringens, and Bacteroids, Faecalibacterium prausnitzii, Parabacteroides, Alistipes were biomarkers that showed significant differences between the groups. In addition, the microbial pathways associated with K. pneumoniae and E. coli may promote hypertension, while A. muciniphila may play a role in reversing the development of hypertension in long-lived elderly. Metabolomics revealed that HTL contained a lower concentration of fecal acetate and propionate than HTE, while it contained a higher concentration of serum acetate and urine acetate. Furthermore, their immune cells exhibited no significant changes in SCFAs receptors.
Conclusion: Although long-lived elderly have extremely high systolic blood pressure, their unique gut microbiota composition and efficient acetate absorption in the colon may offset the damages caused by hypertension and maintain healthy homeostasis.
Introduction
Hypertension is a major cause of cardiovascular disease and has become a public health problem worldwide. The pathogenesis of hypertension involves multiple factors, including genetics, environment, hormones, and inflammation (1). There is increasing evidence that the gut microbiota plays an essential role in the development and pathogenesis of hypertension (2). The gastrointestinal tract is the largest immune interface with the environment in the human body. An association between hypertension and gut microbiota has been found in human and animal models (3). In hypertensive patients, the α-diversity of gut microbiota was decreased and highly correlated with gram-negative flora, including increased genera Klebsiella, Streptococcus, Parabacteria, Desulfovibrio, and Prevotella (4), while some short-chain fatty acid-producing bacteria, such as Faecalibacterium, Roseburia spp will be reduced (5). It was also found that high blood pressure is correlated to gut dysbiosis in spontaneously hypertensive rats (6). Germ-free mice transplanted with fecal bacteria from hypertensive humans developed gut microbiota similar to their donors and increased systolic and diastolic blood pressure after 8 weeks (2). Together, these studies suggest that hypertensive patients have disturbed gut microbiota and an increase in some of the characteristic genera associated with hypertension.
Changes in gut microbiota can lead to changes in microbial-derived metabolites. Among them, short-chain fatty acids (SCFAs), as the product of bacterial fermentation of dietary fiber, are closely related to health (7). SCFAs are produced and absorbed in the colon and excreted in the feces, so fecal SCFAs reflect levels that cannot be absorbed in the colon (8). Previous research indicated that the efficient absorption of acetate, propionate, and butyrate in the colon decreased the blood pressure in hypertensive experimental animals (9, 10), which was further confirmed by the higher levels of fecal SCFAs in hypertensive patients (11). Therefore, fecal SCFAs can be used as a surrogate index of the absorption of SCFAs in the body.
Numerous studies have used the longevity population to study the relationship between gut microbiota and longevity (12–14). The gut microbiota of centenarians is rich in health-related bacteria such as Akkermansia, Christensenellaceae, and Bifidobacterium (15), which produce high levels of specific secondary bile acids and effectively inhibit gram-positive pathogenic bacteria (16). It has a higher ability of glycolysis and fermentation to generate short-chain fatty acids (17). However, long-lived elderly and their offspring have widespread hypertension (18). Most of them live in remote mountain villages, such as Bama in Guangxi (19), Sardinia (17), rural Korea (20), and Moscow (21), etc., are rarely affected by the modern urbanization environment. Bama, Donglan, and Fengshan are the three most famous longevity populations in Guangxi, and which prevalence of the three counties was 35.6, 31.4, and 28.5 centenarians/100,000 citizens in 2019, respectively (13). In our previous study, 66.7% of the population aged >60 in this longevity area had SBP≥140 mmHg (22). This co-existence of longevity and hypertension is well worth investigating. We hypothesize that these long-lived populations have unique gut microbiota and microbial-derived metabolites that offsets hypertension damage and thus maintain healthy homeostasis.
Based on the above, this study conducts a cross-sectional study of gut microbiota, urine metabolites, SCFAs, and their receptors in longevity areas in Donglan County, Guangxi, through three generations of 16S full-length sequencing, metabolomics, and qPCR, respectively. Analyzing the particularity of gut microbiota composition of long-lived people with hypertension, and combing the corresponding metabolites to explain the relationship between longevity and hypertension, can provide a reference for the prevention and treatment of hypertension.
Materials and methods
Participant recruitment and study groups
This research was approved by the Ethics Committee of Guangxi University (approval number: gxdxyxll05). For the experiment design, 46 subjects from Donglan County, Hechi City, Guangxi, were recruited, and all subjects signed written informed consent. The general conditions of the research subjects were recorded, including their age, gender, height, weight, and blood pressure of the subjects. The age information was verified by comparing the ID card or household registration. Clinical and medical history were collected based on the self-report and physical examination reports (Supplementary Table 1). Blood pressure (BP) was measured using a calibrated arm electronic sphygmomanometer (Omron). When measuring blood pressure, subjects were asked to sit still for 5 min with their feet flat on the floor without crossing them, and their upper arms were relaxed and placed on a table at the same level as the heart. Two measurements were taken 1 min apart, and if systolic BP was >15 mm Hg apart or diastolic BP > 10 mm Hg apart, a third measurement was performed and averaged. Inclusion criteria: (i) age ≥ 60 years, and either sex, (ii) body mass index (BMI) 17–25 kg/m2, (iii) systolic blood pressure (SBP) ≥ 140 mmHg, and not using BP-lowering medication, (iv) from autochthonous families, none was coming from other regions originally. Exclusion criteria: (i) history of chronic medical conditions (diabetes, gastroenteritis, chronic kidney disease, morbid obesity, chronic pancreatitis, or other malabsorption disorder), (ii) use of probiotics, prebiotics, or antimicrobial medication (antibiotic or antifungal treatments) 1 year before sampling. Four participants were excluded due to high BMI, four were excluded due to the SBP lower than 140 mmHg, and five were excluded due to chronic medical conditions. One participant dropped out due to unwillingness to continue (Supplementary Figure 1). Following criteria, a total of 32 participants remained: 19 participants aged 60–90 were divided into Hypertension-Elderly (HTE, average 70.5) group, and 13 participants aged 92+ were divided into Hypertension-Longevity (HTL, average 100). The characteristics of the participants are shown in Table 1.
Dietary assessment
Participants’ dietary and nutritional status was assessed using a semi-quantitative food frequency questionnaire (FFQ), which was recorded for consecutive 7-days at home (23, 24). All participants were asked to adopt the mode of individual dining and were urged to maintain their usual eating patterns at home during the dietary survey. Trained investigators used electronic food scales, measuring cups, and spoons to weigh and record each food and drink the participants consumed. For unmeasurable food, the standard portion size was used. Qualified dietitians evaluated and verified the dietary data at a dietary survey location for quality control purposes. Dietary data was converted to nutritional consumption levels according to the 2022 Dietary Guidelines for Chinese Residents. The average daily nutritional intake was determined by multiplying the food consumed (in grams) or the portion size consumed by the nutrient content per 100 g of the food listed in the guideline (25).
16S rRNA third-generation full-length sequencing of gut microbiota
Our study followed guidelines for gut microbiota studies in hypertension (26) and the STORMS (Strengthening the Organization and Reporting of Microbiome Studies) reporting (27) (Supplementary Table 2). Fecal samples of the participants were collected by themselves or their family members with an icebox at home in sterilized empty tubes (for SCFA determination) or tubes containing RNAlater (Thermo Scientific) for microbial DNA extraction. Morning urine samples were collected in sterilized empty tubes with an icebox. Tubes were immediately stored in a small car refrigerator (Alpicool, A30) at −20°C for <24 h. Then, the samples were transported to the laboratory by professionals and stored at −80°C after dispensing into 2 mL sterilized EP tubes until analysis. DNA was extracted using the Stool Genomic DNA Extraction Kit (Solarbio) according to the manufacturer’s instructions. After extraction, PCR amplification of the full-length region of bacterial 16S rRNA was performed using universal primers 27F (5′-AGRGTTTGATYNTGGCTCAG-3′), 1492R (5′-TASGGHTACCTTGTTASGACTT-3′), and purified with magicpure size selection DNA beads (TANGEN Biotech Corporation Ltd., Beijing, China).
Subsequently, the quality inspection was performed on the formed sequencing library, and processing, including barcode recognition, was performed on the high-quality circular consensus sequencing (CCS) sequence obtained. The generated optimization CCS was clustered at 97% similarity (USEARCH, version 10.0), and its species classification was obtained based on the operational taxonomic unit (OTU) sequence composition. The platforms of 16S: Silva database and RDP Classifier were used to analyze species annotation and taxonomy and gut microbiota diversity. Alpha diversity analysis was performed to examine the species richness and diversity within a single sample, and each sample’s ACE, Chao, Shannon, and Simpson indices were calculated. The differences in the community composition and structure of different samples were compared through Beta diversity analysis. Finally, Metastats analysis was carried out to compare the significant differences between the groups at the genus level, and linear discriminant analysis effect size (LEfSe) analysis was used to screen biomarkers that were statistically different between the groups (criteria: LDA score > 3.5). The raw sequencing data generated in this study have been deposited to the NCBI SRA database (BioProject: PRJNA888351).
Fecal short-chain fatty acid measurement
Serum SCFAs were measured in 400 μL, and fecal SCFAs were measured from 0.2 g of fecal sample, all in triplicates. For the fecal sample extraction, 0.2 g of −80°C frozen feces were mixed with 2 mL ultra-pure water in a centrifuge tube. Samples were shaken for 2 min and centrifuged at 12,000 rpm for 20 min (4°C). After centrifugation, 1 mL supernatant was mixed with 0.25 mL of 25% metaphosphoric acid, kept in ice bath at −20°C for 2 h. After thawed, the samples were centrifuged at 12,000 rpm for 10 min (4°C). The supernatant (1 mL) was mixed with 1 mL ethyl acetate, vortexed and centrifuged for 12,000 rpm for 5 min (4°C) to obtain the clear organic layer, which was aspirated and filtered through an organic green membrane into the cold GC glass vials for analysis. For the serum sample extraction, 400 μL serum was mixed with 50 μL of 25% metaphosphoric acid and kept in ice bath at −20°C for 2 h. Follow steps were similar to fecal extraction and took care to add proportion. The analysis of acetate, propionate, isobutyrate, butyrate, isovalerate, and valerate was performed by gas chromatography-mass spectrometry (Agilent 7890A Series) using a capillary DB-Wax column (30 m, 0.25 μm, 0.25 μm; SGE, Cromlab SL, Barcelona, Spain), coupled with a flame ionization detector (GC-FID). The column temperature was programmed at 80°C, reached 110°C at a rate of 10°C/min, maintained for 1 min, then increased to 160°C at 4°C/min, and maintained 1 min (total run time 17.5 min). Helium was the carrier gas (0.5 mL/min). The injection was carried out with a split injector (1:100) at 220°C, the detector temperature was 250°C, and 1 μL of the solution was injected into the GC-FID system. The SCFAs was identified according to standard compounds’ retention time of standard compounds (acetate, propionate, isobutyrate, butyrate, isovalerate, and valerate; Sigma-Aldrich).
Analysis of urine metabolites based on 1H-nuclear magnetic resonance
Aliquots of 500 μL of urine were transferred into a 2.0 mL centrifuge tube and added with 500 μL of 100 mM phosphate buffer containing 0.05% TSP-labeled heavy water (10% heavy water, K2HPO4/NaH2PO4, PH = 7.4). The mixtures were vortexed and centrifuged (15 min, 12,000 × rpm, 4°C) to obtain the urine supernatant.
1H-NMR analysis
The prepared urine supernatant (550 uL) was placed on the Bruker AVANCE 500 MHZ nuclear magnetic resonance spectrometer of the Prodigy liquid nitrogen cryogenic probe for NMR spectrum detection. The pre-saturation method suppressed the water peak, and the NOESY pulse sequence was used. The specific setting parameters are sampling times NS = 64, temperature 25°C, spectrum width SWH = 10,000 Hz, measurement frequency SF = 500.13 MHz, relaxation delay D1 = 2 s, number of sampling points TD = 65,536, sampling time AQ = 3.277 s, mixing time 0.1 s, O1P = 4.7ppm, FID resolution 0.245.
1H-NMR spectrum processing and data analysis
The Fourier transform of the fecal 1H-NMR spectrum was performed using MestReNova NMR data processing software (Mestrelab Research, Spain), and the phases and baseline were adjusted. The chemical shift of TSP (δ 0.0) was taken as the standard to calibrate the chemical shift. The water peak within the range of δ 4.70–5.00 was excised, and the NMR spectra with δ 0.60–9.00 were segmented with equal width in the unit of 0.01 for segment integration, and the data were normalized. The above data were stored in Excel tables for statistical analysis.
Blood expression of SCFA receptors and transporters
The three main SCFA-sensing receptors GPR41 (FFAR3), GPR43 (FFAR2), and GPR109A (HCAR2) have a role in cardiovascular dysfunction (28). As these receptors are highly expressed in immune cells (29), we quantified the mRNA expression of the 3 receptors in circulating immune cells in 30 participants. Since no blood samples were collected from the participants, we selected 30 samples(15 in HTE and 15 in HTL) from long-lived populations with hypertension characteristics collected earlier (23). The basic information of these samples is presented in Supplementary Table 3. Whole blood was treated with Red Blood Cell Lysis Buffer (ThermoFisher Scientific), and RNA was extracted using the Total RNA Extraction Kit (Solarbio). RNA was quantified in a Nanodrop Spectrophotometer, and the first-strand complementary synthesis reaction (cDNA) was made using the BeyoRT™ III First Strand cDNA Synthesis Kit (Beyotime). SYBR Green I was used in a Roche LIGHTCYCLER 96 Real-time quantitative PCR (qPCR) system (Roche Diagnostics Co., Ltd., Basel, Switzerland), with GAPDH as housekeeping genes (Supplementary Table 4). All expression experiments were run in duplicates, and significance was assessed by the 2–ΔΔCT method.
Statistics and analysis
Microbiome data were analyzed as explained above. Graphs were generated using GraphPad Prism version 9.0. Adobe Illustrator version 2020 was used to figure combinations and graphical abstract drawings. All the analyses were performed using SPSS version 26.0 (SPSS Inc., United States) and a 5% significance level. The data were presented as mean ± SD values. The Mann–Whitney U test was used to determine the statistical difference between the two groups gut microbiota, SCFAs, and urine metabolites. A spearman correlation test was conducted to explore correlations between the differential genus and significant microbial pathways and nutrient intake. Pearson correlation analysis was performed to examine correlations between blood pressure and levels of SCFAs, and urine metabolites, after adjustment for age and BMI. Further analyses were conducted using simple linear regression models for BP, acetate, butyrate, and propionate.
Results
Nutrient intake in subjects
The nutrient intake levels in the two groups are presented in Table 2. The two groups had significantly different nutrient intake levels. Notably, the level of energy, fat, saturated fatty acid (SFA), monounsaturated fatty acids (MUFA), polyunsaturated fatty acids (PUFA), carbohydrate, thiamine, riboflavin, vitamin E, folic acid, phosphorus, potassium, calcium, zinc, and manganese were significantly lower in subjects in the HTL compared with the group HTE (p < 0.05). Furthermore, the intake of protein, dietary fiber, cholesterol, vitamin C, vitamin K, nicotinic acid, phosphorus, sodium, magnesium, iron, and selenium were at the same level between the groups. In addition, there were no significant differences in the number of food items consumed by participants, probably due to the similar dietary habits of the participants in the longevity areas. The spearman correlation was used to calculate the relationship between genus and nutrient intake (Supplementary Figure 2). The genus of Faecalibacterium was positively correlated with sodium, whereas the genus of Subdoligranulum was positively correlated with energy, fat, SFA, PUFA, vitamin A, thiamine, riboflavin, vitamin E, phosphorus, magnesium, calcium, iron, and zinc. Clostridium was negatively correlated with vitamin A, iron, and zinc. Christensenellaceae_R7_group had a negative correlation with vitamin C.
Analysis of gut microbiota characteristics between group
Through high-throughput sequencing and sequence screening, 400,724 CCS sequences were obtained through the Barcode identification of 32 samples after sequencing. Each sample generated at least 10,254 CCS sequences, with an average of 12,523 CCS sequences. Operational taxonomic units (OTUs) were divided according to 97% sequence similarity, and a total of 360 OTUs were obtained for subsequent analysis. Through the Shannon curve and the rank abund curve, it was found that the number of OTUs in all samples increased with the increase in sequencing amount and eventually plateaued, indicating that the sequencing volume is sufficient to achieve the ideal sequencing depth (Supplementary Figure 3). The shared and unique OTUs within and between groups were analyzed and plotted as Venn diagrams. HTL had six shared OTU features, HTE had three shared OTU features, and 306 shared OTU features between the two groups (Supplementary Figure 4). Alpha diversity was assessed by five parameters and showed no significant difference between HTE and HTL (Supplementary Table 5).
A total of 11 phylum were dominated by Firmicutes, Proteobacteria, Bacteroidetes, and Verrucomicrobiota, which accounted for 97.32% (HTE) and 99.01% (HTL), respectively. The ratio of Firmicutes to Bacteroidetes(F/B) of HTL (2.14) was significantly lower than HTE (3.71). Figure 1 shows species with relative abundances greater than 1% at the family and genus levels. Both groups were dominated by five families, Enterobacteriaceae, Ruminococcaceae, Akkermansiaceae, Bacteroidaceae, and Lachnospiraceae, with cumulative abundances of 72.61% (HTE) and 75.90% (HTL), respectively. Compared with HTL, the HTE group had much higher levels of the hypertension-related bacteria genera Klebsiella and Streptococcus (+655% and +242%), while it had lower levels of SCFA-producing bacteria Bacteroides (–60%), Roseburia (–80%), Faecalibacterium (–50%), and Alistipes (–79%), respectively. Table 3 further presents the significant difference between the groups in taxa at the genus and species level. The above results showed that the HTE experienced gut dysbiosis, while HTL were abundant in SCFA-producing bacteria. These results are consistent with the commonly reported hypertensive gut genera in addition to identifying two new health-associated bacterial genera (Supplementary Table 6).
Figure 1. Bar chart of the relative abundance (%) gut microbiota profiles. The taxa with average relative abundance above 1% are displayed at the family (A), and genus (B) levels.
The LEfSe method was used to analyze and screen out the biomarker with an LDA score > 3.5. A total of 23 differential species were found between the two groups, Klebsiella pneumoniae, Lactobacillus gasseri, Streptococcus salivarius, Ruminococcus, Actinomyces, Rikenellaceae, f_Saccharimonadaceae, Clostridium perfringens in the HTE and Bacteroids, Faecalibacterium prausnitzii, Parabacteroides, Alistipes in the HTL group, respectively (Figure 2).
Figure 2. The gut microbiome taxa linear discriminant analysis effect sizes. Predicted gut microbiome taxa that are different between hypertension elderly and hypertension longevity, with a linear discriminant analysis (LDA) score of least 3.5.
Hypertension-related microbial pathways
PICRUSTs2 was performed to predict the functional prediction of KEGG microbial pathways in gut microbiota, and a total of 248 Class3 microbial pathways were obtained. The biosynthesis of secondary metabolites, biosynthesis of antibiotics, microbial metabolism in diverse environments, and biosynthesis of amino acids were significantly enriched in both HTE and HTL groups. There were 114 Class3 microbial pathways with relative abundance greater than 0.1% and 16 microbial pathways were significantly different between HTE and HTL groups, including ABC transporters, alanine, aspartate and glutamate metabolism, arginine biosynthesis, benzoate degradation, degradation of aromatic compounds, GABAergic synapse, glycosaminoglycan metabolism, inositol phosphate metabolism, lysosome, other glycan degradation, phosphotransferase system (PTS), propanoate degradation, quorum sensing, selenocompound metabolism, sphingolipid metabolism, and tyrosine metabolism (Supplementary Table 7). The correlations between the top 20 species and 16 significant microbial pathways were further analyzed by Spearman correlation (Figure 3). Interestingly, the association of the hypertensive marker species Klebsiella pneumoniae with microbial pathways was consistent with Escherichia coli, but the correlation was reversed for Akkermansia muciniphila.
Figure 3. Heat map of Spearman correlation between significantly microbial pathways and top 20 abundance species. The X-axis indicates the significantly different microbial pathways, and the Y-axis indicates the top 20 levels of species. Red represents positive correlation, and blue represents negative correlation. The stronger the color, the closer the number is to the Spearman correlation value of 1 or –1, *indicates significant correlation, *p < 0.05, **p < 0.01.
Fecal short-chain fatty acid levels
Six short-chain fatty acids (SCFAs) in fecal samples, namely acetate, propionate, isobutyrate, butyrate, isovalerate, and valerate, were qualitatively detected by GC-MS. Acetate (p = 0.008) and propionate (p = 0.033) were significantly lower in HTL feces than in HTE (Figures 4A, C). Correlation analysis showed that there was no correlation between SBP and SCFAs, but DBP was negatively correlated with acetate (r = −0.368, p = 0.045) and propionate (r = −0.387, p = 0.035) (Figures 4B, D). There were significant correlations between SCFAs, and acetate was significantly positively correlated with propionate (r = 0.708, p = 0.000), butyrate (r = 0.545, p = 0.002), and valerate (r = 0.453, p = 0.012) (correlation results after adjusted for age and BMI).
Figure 4. Fecal levels of gut microbial metabolites short-chain fatty acids (SCFAs) in association with hypertension. (A) Acetate, (B) the correlation analysis between acetate and diastolic blood pressure (DBP), (C) propionate, (D) the correlation analysis between propionate and diastolic blood pressure (DBP). Sample size: n = 19 hypertension elderly and 13 hypertension longevity. *Indicates p < 0.05, **p < 0.01.
Statistical analysis of urine metabolites
Through urine 1H-NMR metabolomic analysis, 36 metabolites were identified (Supplementary Figure 5). Compared with HTE, urine acetate (p < 0.001), pyruvate (p = 0.022), alanine (p = 0.014), and acetoacetate (p = 0.010) were significantly higher in HTL, while creatinine was significantly lower (p = 0.041) (Figure 5). Supplementary Table 8 shows the urine metabolites significantly associated with SBP and DBP. In particular, N-acetyl glycoprotein was negatively correlated with both SBP and DBP (correlation results were adjusted by age and BMI).
Figure 5. Urine metabolites with significant differences between groups. Sample size: n = 19 hypertension elderly and 13 hypertension longevity. *Indicates p < 0.05. **P < 0.01.
SCFAs and receptors
To further explore the relationship between SCFAs and hypertension, we then directly measured the levels of serum SCFAs and plasma SCFAs receptors. Since no blood samples were collected from the participants, 30 samples (15 in HTE and 15 in HTL) from long-lived populations with hypertension characteristics collected earlier were selected (23). Hypertension long-lived elderly had higher serum acetate levels than hypertension elderly (p = 0.011) (Figure 6A). There were no significant differences between the groups’ levels of serum propionate, butyrate (Figures 6B, C), isobutyrate, valerate, and isovalerate. We then quantified the mRNA levels of SCFA-sensing receptors, GPR41, GPR43, and GPR109A, in white blood cells from these samples. There were no changes in the expression of GPR41, GPR43, and GPR019A (Figures 6D–F).
Figure 6. Plasma levels of short-chain fatty acids (SCFAs) and relative expression of short-chain fatty acids (SCFAs)-sensing receptors. (A) Acetate, (B) propionate, (C) butyrate, the mRNA expression in white blood cells between groups, while there are no change in GPR41 (D), GPR43 (E), and GPR109A (F). Sample size: n = 15 hypertension elderly and 15 hypertension longevity. *Indicates p < 0.05.
Discussion
This study is the first cross-sectional study to explore the phenomenon of hypertension in long-lived populations. Microbiome and metabolomic techniques were conducted to characterize the gut microbiota, urine metabolites, and SCFAs of longevity populations with hypertension. The HTE group experienced gut dysbiosis, similar to those reported in hypertensive patients (2, 28), such as higher levels of Klebsiella and Streptococcus, but not in the HTL group. The HTL group had higher Bacteroides, Faecalibacterium, and Alistipes, which may protect the long-lived elderly against hypertension-induced damage. In addition, the microbial pathways associated with K. pneumoniae and E. coli may promote hypertension, while A. muciniphila may play a role in reversing the development of hypertension in long-lived elderly. Fecal SCFAs levels showed that fecal acetate and propionate in the HTL group were significantly lower, while urine acetate was significantly higher than in the HTE group. Furthermore, HTL exhibited higher serum acetate, but their immune cells expressed no significant changes in SCFAs receptors. In conclusion, our study revealed that hypertension elderly experienced gut dysbiosis, while hypertension long-lived elderly had a unique microbiome and efficient acetate absorption capability in the colon, which can offset the damage of hypertension and maintain healthy homeostasis.
Compared to HTL, the HTE group had significantly higher levels of F/B, and more hypertension-related genera increased in the HTE group, such as Klebsiella and Streptococcus. F/B has been used to assess the stability of gut microbiota. A higher value indicates a more disordered intestinal flora and a greater the possibility of hypertension. These conclusions have been verified in hypertensive patients (30) and spontaneously hypertensive rats (6). From the characteristics of gut microbiota, compared with the control population, hypertensive patients always have a higher abundance of Klebsiella, Clostridium, Streptococcus, Porphyromonas, and Actinomyces, and a lower abundance of Bacteroides, Faecalibacterium, Roseburia than the control population (2). In this study, the disturbance of gut microbiota in the HTE group was similar to the previous report, but the HTL group was quite different. The HTL group had higher Bacteroides, Faecalibacterium, and Alistipes than the HTE group. Bacteroides and Alistipes are important acetate and propionate producing bacteria (31), and Faecalibacterium are important butyrate-producing bacteria (32), which are biomarkers to distinguish hypertensive patients in the average population (11). These genera are essential for maintaining health, suggesting that the HTL group has a health-relevant gut microbiome profile.
The search for biomarkers with significant differences between groups (Figure 3) has implications for exploring the phenomenon of hypertension in long-lived people and can explore the relationship between gut microbiota and hypertension at the species level. Bacteroides, F. prausnitzii, Parabacteroides, and Alistipes were enriched in HTL. F.prausnitzii can produce butyrate, which affects blood pressure through vasodilation or plasminogen activator inhibitor-1 (33). Alistipes are a producer of acetate and propionate and play a protective role in cardiovascular diseases (34). Alistipes finegoldii belonging to the genus Alistipes is a species that protects the integrity of the intestinal barrier against colitis (35), which contributes to the absorption of energy substances in the colon. Compared to the HEL group, the HTE had higher levels of Klebsiella pneumoniae, and Streptococcus salivarius, a further extension of the genus-level species (Klebsiella and Streptococcus) reported in previous studies (4). The severity of hypertension was associated with an increased abundance of these two species (36). The enrichment of K. pneumoniae directly contributes to blood pressure elevation and hypertension pathogenesis, and it causes intestinal damage, fecal metabolic changes, and renal shifts may be integrated mediators (23).
The correlation analysis showed that Klebsiella pneumoniae and Escherichia coli might be associated with the microbial pathways of hypertension, while Akkermansia muciniphila was the opposite. As we know, E. coli is an opportunistic pathogen, the most abundant intestinal flora in the longevity population in this study, which was reported to be associated with the development of hypertension during follow-up (37). Moreover, both E. coli and K. pneumoniae have been reported to be significantly enriched in carotid atherosclerosis and hypertension patients, which generated lipopolysaccharide (LPS) (38) and Trimethylamine (TMA) (39). LPS can be transferred to the bloodstream when the intestinal barrier function is impaired, causing asymptomatic hypoendotoxemia (40). TMA is converted by liver enzymes to TMAO, which causes cardiovascular damage (41) and exacerbates angiotensin II-induced hypertension (42). Based on these findings, we inferred that E. coli and K. pneumoniae have similar microbial pathways and may be involved in the development of hypertension. On the contrary, the abundance of A. muciniphila is negatively associated with obesity, type 2 diabetes, and hypertension (43, 44). Therefore, as a new generation of probiotics, the microbial pathway of A.muciniphila is significantly different from that of K. pneumoniae and E. coli. It may play a role in reversing the development of hypertension, and its specific mechanism needs further exploration.
Besides gut microbiota, the changes in SCFAs were also correlated to blood pressure (45). The concentrations of fecal acetate and propionate in the HTL group were significantly lower, while urine acetate was significantly higher than HTE. Furthermore, HTL exhibited higher serum acetate, but their immune cells expressed no significant changes in SCFAs receptors. A common explanation is that the lower fecal SCFAs in the long-lived elderly may be due to the reduction of SCFA-producing bacteria and the decreased metabolic activity of the gut microbiota (46). However, the relative abundance of SCFA-producing bacteria such as Bacteroides, Faecalibacterium, and Alistipes was higher in HTL than in HTE, and there was no significant difference in dietary fiber intake between the two groups, so it would not be lower from the perspective of SCFAs production. Everything gets interesting when explained in terms of SCFA uptake levels. SCFAs are efficiently absorbed in the colon, and less than 5% are excreted in the feces. Fecal SCFAs reflect SCFAs that cannot be absorbed in the host colon (3) and are surrogate indicators of colonic SCFAs absorption. Fecal SCFAs concentrations are more representative of SCFAs uptake than their production (8). Fecal SCFAs are significantly increased, and serum SCFAs are significantly decreased in hypertensive patients, indicating that the adequate absorption level of SCFAs in hypertensive patients is lower (11). Higher levels of fecal SCFAs are associated with gut dysbiosis, obesity, and hypertension, presumably due to the lower absorption rate of SCFAs in the gut (47). In contrast, participants with lower blood pressure had lower levels of fecal SCFAs and higher levels of SCFA-producing bacteria, which may promote intestinal absorption of SCFAs (5). Above all, these SCFA-producing bacteria have two functions: they produce SCFAs and promote the absorption of SCFAs in the human intestine, with their activity contributing to the maintenance of low levels of fecal SCFAs and higher levels of blood SCFAs and even urine of the host.
More importantly, we noticed that the acetate level in the HTL group was lower in the fecal and higher in the blood and urine compared to the HTE group. It is important to note that most SCFAs are consumed in the gastrointestinal tract and liver. The remaining SCFAs (mainly acetate) can be transported to other organs via the bloodstream to regulate the host’s metabolism, physiology and energy balance (48). Acetate and propionate both enter the portal venous circulation and peripheral blood, but they play different roles. The propionate is mainly used for glycoisogenesis in the liver, which is then metabolized to CO2 through the tricarboxylic acid cycle and other metabolic pathways. In contrast, acetate enters the systemic circulation and reaches peripheral tissues, which activates signaling via binding to 3 G-protein–coupled receptors, with GPR43 being the most predominant of these receptors (9). These receptors are highly expressed in immune cells, including T, B, and innate lymphoid cells, and activate anti-inflammatory downstream pathways to decrease blood pressure (49). The high absorption of acetate in the colon and more urine excretion of the HTL group supports the hypothesis that the HTL group have unique gut microbiota and efficient acetate absorption capacity, which can offset hypertension damage and thus maintain healthy homeostasis.
Although this research used a relatively small sample size and did not include a healthy control group, our study took advantage of the only longevity and hypertension cohort published to date with well-characterized BP monitoring in both men and women, all untreated for BP-lowering medication. This study is also the only cohort with detailed information regarding diet, urine metabolites, plasma, fecal SCFAs, and their receptors, which allowed us to explore the interplay between the gut microbiome, diet, urine metabolites, SCFAs, and their receptors regard to longevity. Future studies should expand the sample size and use ambulatory BP monitoring to minimize intra-individual variation in samples. Meanwhile, due to the limitation of 16S rRNA sequencing, future studies will use metagenome to uncover the species that play a crucial role in the association between longevity and hypertension.
Conclusion
In conclusion, we found that people in longevity areas of Guangxi, China, generally with hypertension. Hypertension elderly had much higher levels of hypertension-related bacteria Klebsiella and Streptococcus, while hypertension long-lived elderly had higher levels of short-chain fatty acid-producing bacteria (Bacteroides, Faecalibacterium, Alistipes). Moreover, K. pneumoniae and E. coli may play a role in promoting hypertension, while A. muciniphila may play a role in reversing the development of hypertension. SCFAs further suggest that acetate can be effectively absorbed in the colon of long-lived elderly with hypertension to offset the negative effects of hypertension. An important mechanism may be driven by SCFA-producing bacteria and their main metabolite, acetate, associated with hypertension longevity. The enrichment of SCFA-producing bacteria, such as Bacteroides spp, Faecalibacterium spp, and Alistipes spp may promote efficient acetate absorption in the colon and represent a new target for BP therapy in the future.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA888351.
Ethics statement
The studies involving human participants were reviewed and approved by the Medical Ethics Committee, Guangxi University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
QL, ZQ, and QZ conceived and designed the experiments. QZ, NM, YL, HZ, RL, KH, HL, XY, and JM contributed sample collection, investigation, and data analysis. QZ wrote the original draft. QL, DH, ZZ, YL, and ZQ contributed to review and edit. QL, ZQ, and YL contributed to final approval of the publication. QL, ZQ, YL, and QZ are accountable for all the aspects of the work. All authors have read and agreed to the published version of the manuscript.
Funding
This study was funded by the National Natural Science Foundation of China (31871802 and 81971505), the Guangxi Key Research and Development Program (AB18221065), and the National Key R&D Program of China (2018YFA0108304).
Acknowledgments
We are very grateful to all the participants and the technical support from BMKCloud (www.biocloud.net).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2022.1070223/full#supplementary-material
References
1. Giri A, Hellwege JN, Keaton JM, Park J, Qiu C, Warren HR, et al. Trans-ethnic association study of blood pressure determinants in over 750,000 individuals. Nat Genet. (2019) 51:51–62. doi: 10.1038/s41588-018-0303-9
2. Li J, Zhao F, Wang Y, Chen J, Tao J, Tian G, et al. Gut microbiota dysbiosis contributes to the development of hypertension. Microbiome. (2017) 5:14. doi: 10.1186/s40168-016-0222-x
3. Avery EG, Bartolomaeus H, Maifeld A, Marko L, Wiig H, Wilck N, et al. The gut microbiome in hypertension. Circ Res. (2021) 128:934–50.
4. Yan Q, Gu Y, Li X, Yang W, Jia L, Chen C, et al. Alterations of the gut microbiome in hypertension. Front Cell Infect Microbiol. (2017) 7:381. doi: 10.3389/fcimb.2017.00381
5. Verhaar BJH, Collard D, Prodan A, Levels JHM, Zwinderman AH, Bäckhed F, et al. Associations between gut microbiota, faecal short-chain fatty acids, and blood pressure across ethnic groups: the HELIUS study. Eur Heart J. (2020) 41:4259. doi: 10.1093/eurheartj/ehaa704
6. Yang T, Santisteban MM, Rodriguez V, Li E, Ahmari N, Carvajal JM, et al. Gut dysbiosis is linked to hypertension. Hypertension. (2015) 65:1331–40.
7. Huart J, Leenders J, Taminiau B, Descy J, Saint-Remy A, Daube G, et al. Gut microbiota and fecal levels of short-chain fatty acids differ upon 24-Hour blood pressure levels in men. Hypertension. (2019) 74:1005–13. doi: 10.1161/HYPERTENSIONAHA.118.12588
8. Vogt J, Wolever T. Fecal acetate is inversely related to acetate absorption from the human rectum and distal colon. J Nutr. (2003) 133:3145–8. doi: 10.1093/jn/133.10.3145
9. Kaye DM, Shihata WA, Jama HA, Tsyganov K, Ziemann M, Kiriazis H, et al. Deficiency of prebiotic fiber and insufficient signaling through gut metabolite-sensing receptors leads to Cardiovascular disease. Circulation. (2020) 141:1393–403. doi: 10.1161/CIRCULATIONAHA.119.043081
10. Kim S, Goel R, Kumar A, Qi Y, Lobaton G, Hosaka K, et al. Imbalance of gut microbiome and intestinal epithelial barrier dysfunction in patients with high blood pressure. Clin Sci (Lond). (2018) 132:701–18.
11. Calderón-Pérez L, Gosalbes MJ, Yuste S, Valls RM, Pedret A, Llauradó E, et al. Gut metagenomic and short chain fatty acids signature in hypertension: a cross-sectional study. Sci Rep. (2020) 10:6436. doi: 10.1038/s41598-020-63475-w
12. Wilmanski T, Diener C, Rappaport N, Patwardhan S, Wiedrick J, Lapidus J, et al. Gut microbiome pattern reflects healthy ageing and predicts survival in humans. Nat Metab. (2021) 3:274.
13. Ren M, Li H, Fu Z, Li Q. Succession analysis of gut microbiota structure of participants from long-lived families in hechi. Guangxi, China. Microorganisms. (2021) 9:2524. doi: 10.3390/microorganisms9122524
14. Odamaki T, Kato K, Sugahara H, Hashikura N, Takahashi S, Xiao JZ, et al. Age-related changes in gut microbiota composition from newborn to centenarian: a cross-sectional study. BMC Microbiol. (2016) 16:90. doi: 10.1186/s12866-016-0708-5
15. Biagi E, Franceschi C, Rampelli S, Severgnini M, Ostan R, Turroni S, et al. Gut microbiota and extreme longevity. Curr Biol. (2016) 26:1480–5.
16. Sato Y, Atarashi K, Plichta DR, Arai Y, Sasajima S, Kearney SM, et al. Novel bile acid biosynthetic pathways are enriched in the microbiome of centenarians. Nature. (2021) 599:458–64. doi: 10.1038/s41586-021-03832-5
17. Wu L, Zeng T, Zinellu A, Rubino S, Kelvin DJ, Carru CA. Cross-Sectional study of compositional and functional profiles of gut microbiota in sardinian centenarians. Msystems. (2019) 4:e319–25. doi: 10.1128/mSystems.00325-19
18. Hai P, Yao D, Zhao R, Dong C, Saymuah S, Pan Y, et al. BMI, blood pressure, and plasma lipids among centenarians and their offspring. Evid-Based Compl Alt. (2022) 2022:3836247. doi: 10.1155/2022/3836247
19. Wang F, Yu T, Huang G, Cai D, Liang X, Su H, et al. Gut microbiota community and its assembly associated with age and diet in chinese centenarians. J Microbiol Biotechn. (2015) 25:1195. doi: 10.4014/jmb.1410.10014
20. Kim B, Choi CW, Shin H, Jin S, Bae J, Han M, et al. Comparison of the gut microbiota of centenarians in longevity villages of South Korea with those of other age groups. J Microbiol Biotechn. (2019) 29:429–40. doi: 10.4014/jmb.1811.11023
21. Eruslanova K, Runichina N, Kotovskay J, Tkacheva O. Abstract P119: prevalence and coincidence of arterial hypertension and orthostatic hypotension in pre and centenarians: preliminary results of moscow centenarian study. Hypertension. (2019) 74(Suppl. 1):A119.
22. Li H, Ren M, He Q, Gao J, Li Q. Revealing the longevity code of humans with up to extreme longevity in guangxi based on physical examination indicators and personalized biomarkers of aging. Biomed Res Int. (2022) 2022:1–10. doi: 10.1155/2022/2810379
23. Li H, Ren M, Li Q. 1H NMR-Based metabolomics reveals the intrinsic interaction of age, plasma signature metabolites, and nutrient intake in the longevity population in Guangxi, China. Nutrients. (2022) 14:2539. doi: 10.3390/nu14122539
24. Ke L, Toshiro T, Fengyan S, Ping Y, Xiaoling D, Kazuo T. Relative validity of a semi-quantitative food frequency questionnaire versus 3 day weighed diet records in middle-aged inhabitants in Chaoshan area, China. Asian Pac J Cancer Prev. (2005) 6:376–81.
25. Wang Y, Feng L, Zeng G, Zhu H, Sun J, Gao P, et al. Effects of cuisine-based Chinese heart-healthy diet in lowering blood pressure among adults in china: multicenter, single-blind, randomized, parallel controlled feeding trial. Circulation. (2022) 146:303–15. doi: 10.1161/CIRCULATIONAHA.122.059045
26. Marques FZ, Jama HA, Tsyganov K, Gill PA, Rhys-Jones D, Muralitharan RR, et al. Guidelines for transparency on gut microbiome studies in essential and experimental hypertension. Hypertension. (2019) 74:1279–93. doi: 10.1161/HYPERTENSIONAHA.119.13079
27. Mirzayi C, Renson A, Furlanello C, Sansone S, Zohra F, Elsafoury S, et al. Reporting guidelines for human microbiome research: the STORMS checklist. Nat Med. (2021) 27:1885–92. doi: 10.1038/s41591-021-01552-x
28. Marques FZ, Mackay CR, Kaye DM. Beyond gut feelings: how the gut microbiota regulates blood pressure. Nat Rev Cardiol. (2018) 15:20–32. doi: 10.1038/nrcardio.2017.120
29. Nakai M, Ribeiro RV, Stevens BR, Gill P, Muralitharan RR, Yiallourou S, et al. Essential hypertension is associated with changes in gut microbial metabolic pathways: a multisite analysis of ambulatory blood pressure. Hypertension. (2021) 78:804–15. doi: 10.1161/HYPERTENSIONAHA.121.17288
30. Antza C, Stabouli S, Kotsis V. Gut microbiota in kidney disease and hypertension. Pharmacol Res. (2018) 130:198–203.
31. Kim S, Rigatto K, Gazzana MB, Knorst MM, Richards EM, Pepine CJ, et al. Altered gut microbiome profile in patients with pulmonary arterial hypertension. Hypertension. (2020) 75:1063–71. doi: 10.1161/HYPERTENSIONAHA.119.14294
32. Machiels K, Joossens M, Sabino J, De Preter V, Arijs I, Eeckhaut V, et al. A decrease of the butyrate-producing species Roseburia hominis and Faecalibacterium prausnitzii defines dysbiosis in patients with ulcerative colitis. Gut. (2014) 63:1275–83. doi: 10.1136/gutjnl-2013-304833
33. Gomez-Arango LF, Barrett HL, McIntyre HD, Callaway LK, Morrison M, Dekker NM. Increased systolic and diastolic blood pressure is associated with altered gut microbiota composition and butyrate production in early pregnancy. Hypertension. (2016) 68:974–81. doi: 10.1161/HYPERTENSIONAHA.116.07910
34. Parker BJ, Wearsch PA, Veloo ACM, Rodriguez-Palacios A. The genus alistipes: gut bacteria with emerging implications to inflammation, cancer, and mental health. Front Immunol. (2020) 11:906. doi: 10.3389/fimmu.2020.00906
35. Dziarski R, Park SY, Kashyap DR, Dowd SE, Gupta D. Pglyrp-Regulated gut microflora Prevotella falsenii, Parabacteroides distasonis and Bacteroides eggerthii Enhance and Alistipes finegoldii attenuates colitis in mice. PLoS One. (2016) 11:e0146162. doi: 10.1371/journal.pone.0146162
36. Jin L, Shi X, Yang J, Zhao Y, Xue L, Xu L, et al. Gut microbes in cardiovascular diseases and their potential therapeutic applications. Protein Cell. (2021) 12:346–59.
37. Meiland R, Geerlings SE, Stolk RP, Hoepelman AI, Peeters PH, Coenjaerts FE, et al. Escherichia coli bacteriuria in female adults is associated with the development of hypertension. Int J Infect Dis. (2010) 14:e304–7. doi: 10.1016/j.ijid.2009.05.015
38. Chen J, Qin Q, Yan S, Yang Y, Yan H, Li T, et al. Gut microbiome alterations in patients with carotid atherosclerosis. Front Cardiovasc Med. (2021) 8:739093. doi: 10.3389/fcvm.2021.739093
39. Zheng S, Piao C, Liu Y, Liu X, Liu T, Zhang X, et al. Glycan biosynthesis ability of gut microbiota increased in primary hypertension patients taking antihypertension medications and potentially promoted by macrophage-adenosine monophosphate-activated protein kinase. Front Microbiol. (2021) 12:719599. doi: 10.3389/fmicb.2021.719599
40. Violi F, Cammisotto V, Bartimoccia S, Pignatelli P, Carnevale R, Nocella C. Gut-derived low-grade endotoxaemia, atherothrombosis and cardiovascular disease. Nat Rev Cardiol. (2022) 20:1–14. doi: 10.1038/s41569-022-00737-2
41. Tang WHW, Li DY, Hazen SL. Dietary metabolism, the gut microbiome, and heart failure. Nat Rev Cardiol. (2019) 16:137–54. doi: 10.1038/s41569-018-0108-7
42. Jiang S, Shui Y, Cui Y, Tang C, Wang X, Qiu X, et al. Gut microbiota dependent trimethylamine N-oxide aggravates angiotensin II–induced hypertension. Redox Biol. (2021) 46:102115. doi: 10.1016/j.redox.2021.102115
43. Xu Y, Wang N, Tan HY, Li S, Zhang C, Feng Y. Function of Akkermansia muciniphila in obesity: interactions with lipid metabolism, immune response and gut systems. Front Microbiol. (2020) 11:219. doi: 10.3389/fmicb.2020.00219
44. Depommier C, Everard A, Druart C, Plovier H, Van Hul M, Vieira-Silva S, et al. Supplementation with Akkermansia muciniphila in overweight and obese human volunteers: a proof-of-concept exploratory study. Nat Med. (2019) 25:1096–103. doi: 10.1038/s41591-019-0495-2
45. Yang F, Chen H, Gao Y, An N, Li X, Pan X, et al. Gut microbiota-derived short-chain fatty acids and hypertension: mechanism and treatment. Biomed Pharmacother. (2020) 130:110503.
46. Badal VD, Vaccariello, Murray ER, Yu KE, Knight R, Jeste DV, et al. The gut microbiome, aging, and longevity: a systematic review. Nutrients. (2020) 12: 3759.
47. de la Cuesta-Zuluaga J, Mueller N, Álvarez-Quintero R, Velásquez-Mejía E, Sierra J, Corrales-Agudelo V, et al. Higher fecal short-chain fatty acid levels are associated with gut microbiome dysbiosis, obesity, hypertension and cardiometabolic disease risk factors. Nutrients. (2019) 11:51. doi: 10.3390/nu11010051
48. Kong D, Schipper L, van Dijk G. Distinct effects of short chain fatty acids on host energy balance and fuel homeostasis with focus on route of administration and host species. Front Neurosci. (2021) 15:755845. doi: 10.3389/fnins.2021.755845
Keywords: gut microbiota, hypertension, longevity, short-chain fatty acids, urine metabolites
Citation: Zhang Q, Meng N, Liu Y, Zhao H, Zhao Z, Hao D, Li R, Han K, Li H, Ma J, Yu X, Qi Z and Li Q (2023) Protection effect of gut microbiota composition and acetate absorption against hypertension-induced damages on the longevity population in Guangxi, China. Front. Nutr. 9:1070223. doi: 10.3389/fnut.2022.1070223
Received: 14 October 2022; Accepted: 20 December 2022;
Published: 16 January 2023.
Edited by:
Daniela Caetano Gonçalves, Federal University of São Paulo, BrazilReviewed by:
Geovana Leite, University of São Paulo, BrazilTianhao Liu, Affiliated Hospital of Jiangnan University, China
Claudia Alves, Universidade Federal de São Paulo, Brazil
Copyright © 2023 Zhang, Meng, Liu, Zhao, Zhao, Hao, Li, Han, Li, Ma, Yu, Qi and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Quanyang Li, bGlxdWFueWFuZ0BneHUuZWR1LmNu; Zhongquan Qi, eXh5eXpAZ3h1LmVkdS5jbg==