Detection of fecal flora
Research object:In this experiment, we enrolled patients who were hospitalized in the Department of Neurology and Cardiology of the Second Affiliated Hospital of Shandong First Medical University in April 2021, including 7 patients with hypertension complicated with ischemic stroke and only 7 patients with hypertension.
The admission criteria of hypertension complicated with ischemic stroke group are: (1) Age over 18 years old; (2) Clinical diagnosis of acute cerebral infarction, and the onset time is less than 10 days; (3) Volunteer to participate in this research. Exclusion criteria: (1) Connective tissue disease; (2) diabetes; (3) suffering from acute and chronic digestive diseases; (4) Patients who are delirious and unable to collect stool samples correctly.
The admission criteria of the hypertension control group: are (1) Age over 18 years old; (2) Volunteer to participate in this research. Exclusion criteria: (1) Previous history of cerebral infarction; (2) Connective tissue disease; (3) diabetes; (4) suffering from acute and chronic digestive diseases; (5) Patients who are delirious and unable to collect stool samples correctly.
During the hospitalization, the feces were collected, and the blood lipid, blood glucose, blood routine, and other test indexes were recorded. Collect the clinical information of the experimental group and the control group, including gender, age, current medical history, past history, etc. This study was approved by the ethics committee of our hospital, and all participating researchers were informed.
Collection of feces:Use a sterile cotton swab to dip 50mg of the middle stool of each of the patients with hypertension complicated with stroke and those with hypertension only, immerse it in the sampling tube of cereal flora storage buffer, rotate and stir it for 30 seconds, then tighten the cover of the sampling tube, shake it left and right for 10 seconds to complete the mixing, record the name and collection time of the sample, store it at room temperature, and send it for inspection within one week.
DNA extraction:Use the DNA extraction kit GHFDE100 (Zhejiang Hangzhou Equipment Preparation 20190952) to extract the metagenomic DNA from the sample according to the manufacturer's instructions.
PCR amplification and 16S rRNA detection method:The V4 region of bacterial 16S rRNA gene was amplified by PCR using forward primer 515f (515F(5‘-GTGCCAGCMGCCGCGGTAA-3'-3') and reverse primer 806r (5' ggactachvgggtwtctaat-3'). Sample-specific paired-end 6-bp barcodes are integrated into the TrueSeq adapter for multiple sequencing. The polymerase chain reaction consists of 25μl Phusion High-Fidelity PCR Master Mix, 3μ l (10μ m) forward and reverse primers, 10μl template DNA,3μl dimethyl sulfoxide, and 6μl DDH_2O. Heat cycle includes initial denaturation at 98℃ for the 30s, then 25 cycles, including denaturation at 98℃ for 15s, annealing at 58℃ for 15s, renaturation at 72℃ for 15s, and extension at 72℃ for 1min. PCR amplification products were purified with Agencourt AMPure XP magnetic beads and quantified with PicoGreen dsDNA Assay Kit. After the individual quantitative step, the amplified fragments were mixed in equal amounts, and the double-ended 2×150bp sequencing was carried out using the IllLumina NovaSeq6000 platform of Guhe Information Technology Co., Ltd.
Statistics and analysis of OTU data:The original sequencing readings that exactly match the barcode are assigned to the corresponding samples and identified as valid sequences. Low-quality sequences were screened by the following criteria: sequences with a length of < 150 bp, sequences with an average Phred score of < 20, sequences with fuzzy bases, and sequences with single nucleotide repeats of > 8 bp. Use Vsearch v2.4.4 (-fastq _ mergepairs-fastq _ minovlen5) to assemble end-to-end reading, and then use Vsearch v2.15.0 to select the operational taxon (OTU), including de-replication, clustering and chimera detection. A microbial quantitative insight pipeline is used to process OTUs and sequencing data. According to the representative sequence set of SILVA138 database, OTU taxonomy was carried out.
Further, the OTU table is generated to record the abundance of each OTU in each sample and the classification of these OTUs. All samples were discarded, which contained less than 0.001% OTU of the total sequence. In order to minimize the difference of sequencing depth among samples, an average and rounded thin OTU table is generated by averaging 100 uniformly resampled OTU subsets at the minimum sequencing depth of 90% for further analysis.
Bioinformatics and Statistical Analysis:Sequence analysis is mainly carried out by using QIIME2 and R software packages (version 3.2.0). Use OTU table in QIIME2 to calculate OTU-level α diversity index, such as Chao1 richness estimate, ACE measure (coverage estimate based on abundance), PD_whole_tree, Shannon diversity index and Simpson index. Generate the ordered abundance curve of OTU level to compare the abundance and uniformity of OTUs among samples. UniFrac distance measurement was used for β diversity analysis to investigate the structural changes of microbial communities in the samples, and it was visualized by principal coordinate analysis (PCoA) and non-metric multidimensional scale (NMDS). Students' T-test and Monte Carlo permutation test with 1000 permutations were used to determine the difference of Unifrac distance between groups, and visualized by box diagram. Principal component analysis (PCA) is also based on the composition profile of genus level. R package "vegan" was used to evaluate the significance of microbial community structure differences among groups by PERMANOVA (displacement multivariate analysis of variance). Venn diagram is generated by using R package Venn diagram according to the occurrence of OTU in samples/groups (regardless of its relative abundance) to show the shared and unique OTU between samples or groups.
The abundance of taxa at the level of phylum, class, order, family, genus and species in a sample or group was statistically compared by Kruskal.test in R stats software package. LEfSe (Linear Discriminant Analysis Effect Size) is executed with default parameters to detect the difference between different groups and enrich the classification groups. Random forest analysis is used to distinguish samples from different groups using R package "Random Forest", which has 1,000 trees and all default settings. Use 10x cross-validation to estimate generalization error. The expected "baseline" error of is also included, which is obtained by simply predicting the classifiers of the most common category labels. Co-occurrence analysis was carried out by calculating Spearman rank correlation among dominant taxonomic groups.
Based on high-quality sequences, the microbial function was predicted by PICRUST 2 (phylogenetic investigation of the community by reconstructing unobserved state, https://github.com/picrust/picrust2/). The bacterial function database and gut-brain module (GBMs) were also analyzed by using human gut metabolism module. For each sample, the function is analyzed using OMixer-RPM version 1.0 (https://github.com/raeslab/omixer-rpm) with KO redundancy predicted from PICRUSt 2.
The blood serum of patients hospitalized in the Second Affiliated Hospital of Shandong First Medical University in September 2022 was collected, including 2 patients with hypertension complicated with ischemic stroke and only 3 patients with hypertension. Add 2 L serum sample into 96-well plate, and further dilute 0.1M nh4hc 3 (Sigma, Missouri, USA), 6 M urea (Sigma, Missouri, USA), and 2 M thiourea (Sigma, Missouri, USA) with lysis buffer containing the following components. The disulfide bond was reduced at 33℃ for 45 minutes at 600 rpm, and 10 mM tris (2- carboxyethyl) phosphine hydrochloride (TCEP) was added and alkylated with iodoacetamide. All crude protein extracts were precipitated with acetone at 20℃ for 2 h, and then re-dissolved with 200μl Teab buffer. Then it was digested with trypsin at 37℃ overnight, and then TFA was added to extinguish it. Sep-Pak column of Waters Co was used to desalt these peptides.
Orbitrap Eclipse mass spectrometer and Easy-NLC 1200 system were used for data acquisition, and the chromatographic column was ACRADEST PepMapTM 100 C18 column (75μ m× 2cm, nanviper, 3μ m, 100, Thermoscience). The flow rate is 300nL/min. Mobile phase A is 5% acetonitrile +0.1% formic acid (v/v, the best LC/MS grade), mobile phase B is 80% acetonitrile +0.1% formic acid (v/v, the best LC/MS grade), mobile phase B is 3% ~ 5% B, 0.05 min; 5%~15%B,23.55 m in; 15%~28%B,21 m in; 28%~38%B,7.3 m in; 38%~100%B,0.05 m in; 100%B 7.25 m in。 The MS conditions are as follows: ion funnel Rf is 40,Ms1 Orbitrap resolution is 60000 (m/z 200), and MS1 maximum sampling time is 20ms. The charge state of the polypeptide is 2-7. Data-dependent acquisition mode was used to trigger precursor separation and sequencing. The MSOrbit RAP resolution is set to 15000(m/z200), and the high energy collision induced dissociation (HCD) with normalized collision energy is 30%. The MS2 isolation window is set to 1.6 Da.
Thermoproteome Discoverer (PD, 2.4.1.15) was used to carry out feature detection, database search and protein/polypeptide quantification on all original documents. The identity of MS polypeptide sequence and protein was determined by matching fragment patterns in the UniProt mouse database (downloaded on March 23, 2021, containing 17 040 sequences). Aminomethylation of cysteine residues is set as fixed modification. There are only two missing cleavages per peptide. The fragment quality tolerance is 0.02Da. Label-free quantification (LFQ) was calculated in each parameter group. Select unique polypeptide and razor polypeptide for protein quantification. Mass spectrometry protein omics data has been stored in ProteomeXchange consortium (http://proteomesscentral. proteomexc ange.org) through Pride partner repository with data set identifier PXD030026. Bioinformatics analysis was carried out using David, Cytoscape and OMICStudio tools. Only protein with p value < 0.05 was included in GO concentration analysis and network formation.
2.3 Analysis of hematological indicators
research object:This experiment enrolled the inpatients in the Department of Applied Neurology and Cardiology of the Second Affiliated Hospital of Shandong First Medical University from 2021 to 2022, including 80 patients with hypertension complicated with ischemic stroke and 48 patients with hypertension only.
The admission criteria of hypertensive stroke group are: (1) Age over 18 years old; (2) The clinical diagnosis is acute cerebral infarction, and the onset time is less than 10 days. Exclusion criteria: (1) Connective tissue disease; (2) diabetes; (3) suffering from acute and chronic digestive diseases; (4) Patients who are delirious and unable to collect stool samples correctly. The admission criteria of hypertension control group: (1) Age over 18 years old. Exclusion criteria: (1) Previous history of cerebral infarction; (2) Connective tissue disease; (3) Diabetes. The blood lipid, blood sugar, blood routine and other test indexes were recorded. Collect the clinical information of the experimental group and the control group, including gender, age, current medical history, past history, etc. This study was approved by the ethics committee of our hospital, and all participating researchers were informed.
detection methods:After fasting for 12 hours, 5ml venous blood was collected at 7-9am the next morning, serum was centrifuged, and seven blood lipids were analyzed by immunoturbidimetry. BECKMAN COULTER AU5800 automatic biochemical analyzer is adopted, and the operation procedures are strictly followed. The daily maintenance and quality control of the instrument are well done.
Reference range immunoturbidimetry :
Total cholesterol (CHOL): 2. 33 ~ 5.69 mmol/L.
TG (triglyceride): 0. 58 ~ 1.88 mmol/L.
High-density lipoprotein (HDL): 1. 03 ~ 1. 85mmol /L
Low-density lipoprotein (LDL): 2. 07-3. 33 mmol/L.
Apolipoprotein A (ApoA): 1 ~ 1. 6g /L
Apolipoprotein B (ApoB): 0. 6 ~ 1. 1g /L
Free fatty acid (NEFA): 0. 1 ~ 0. 6mmol /L
Statistical analysis:SPSS25.0 software is used, and the data is expressed by the mean standard deviation. If the data is normal, two independent samples T-test is adopted, but it is not normal, and a nonparametric rank sum test is adopted.
ELISA detection
This experiment enrolled the inpatients in the Department of Applied Neurology and Cardiology of the Second Affiliated Hospital of Shandong First Medical University in 2022, including 6 hypertensive patients with ischemic stroke and 4 hypertensive patients.
The admission criteria of the hypertensive stroke group are: (1) Age over 18 years old; (2) The clinical diagnosis is acute cerebral infarction, and the onset time is less than 10 days. Exclusion criteria: (1) Connective tissue disease; (2) diabetes; (3) suffering from acute and chronic digestive diseases; (4) Patients who are delirious and unable to collect stool samples correctly.The admission criteria of the hypertension control group: (1) Age over 18 years old. Exclusion criteria: (1) Previous history of cerebral infarction; (2) Connective tissue disease; (3) Diabetes. After fasting all night (at least 10 hours), blood samples were collected from the subjects' anterior elbow vein from 7:30 am to 9:00 am. The collected blood was coagulated at 4 C for 30 min, and then centrifuged at 3000 rpm at 4 C for 20 min. Immediately store the analysis sample at-20 C to avoid repeated freeze-thaw cycles. ELISA kit (Jiangsu Enzyme-free Industry Co., Ltd., Jiangsu, China, catalogue number E1766H1) was used to determine the concentration of APOC2 and APOC3, and the operation was carried out according to the manufacturer's plan.