Ethics statement
All animal experiments were conducted according to the recommendations in the Guide for the Care and Use of Laboratory Animals published by the US National Institutes of Health and approved by the Ethics Committee of Hebei Medical University (No.IACUC-Hebmu-2021035).
Experimental Animals
Male C57BL/6J mice of 6-7-weeks age were purchased from Beijing HFK Bio-Technology Co., Ltd. and were housed in a specific pathogen free (SPF) animal facility at Hebei Medical University Laboratory Animal Center. All animals were housed at constant temperature (20 ± 2oC) and humidity (45–55%), with a 12 hours light/12 dark cycle and free diet and water ad libitum.
Establishment of mouse model of DM and hUC-MSC administration: The groups of experimental animals described here are also the subject of a separate manuscript describing the effects of hUC-MSCs on diabetic nephropathy (He et al., currently under review). A STZ model of DM with low toxicity was established as previously described [53] by intraperitoneal (i.p.) injection of STZ (Cayman Chemicals, catalog no. 18883664) 80 mg/kg body weight in citrate buffer pH 4.5 daily for 5 consecutive days. For each experiment, a control group of non-diabetic (Non-DM) mice was generated by 5 consecutive i.p. injections of citrate buffer alone. Blood glucose concentrations were measured once a week in all mice using a glucometer (Roche, Accu-Chek) and blood glucose test strips (Roche, Excellence). Mice with fasting blood glucose (FBG) levels ≥ 16.7 mM were considered to have DM. Two in vivo experiments were performed of 10 weeks and 18 weeks duration following induction of DM. For each, mice were divided into three groups: Non-DM group, DM group (treated with i.v. injection of sterile saline 2 weeks before termination of the experiment) and DM + MSC group (treated with i.v. injection of hUC-MSCs in sterile saline 2 weeks before termination of the experiment). The experimental designs are illustrated in Fig. 1A. For i.v. treatments, 5 x 105 freshly cultured hUC-MSCs suspended in 0.2 ml of sterile saline or 0.2 ml of sterile saline alone were injected via the tail vein into mice with confirmed DM. Mice were monitored weekly for fasting blood glucose and body weight. An Animal General Distress Scoring Sheet was completed on a daily basis with humane interventions undertaken as appropriate (see Table S1 at the end of the manuscript).
Culture, Identification and Characterization of hUC-MSCs
Cryopreserved hUC-MSCs were purchased from Qilu Cell Therapy Engineering Technology Co., Ltd (Shandong, China). Upon receipt, the hUC-MSCs were thawed, transferred to complete medium (DMEM, low glucose medium with 15% FBS,100 units/ml penicillin, and 100 mg/ml streptomycin) and cultured in a humidified 37°C, 5% CO2 incubator. The medium was changed every 3 days. When the hUC-MSCs reached 80–90% confluence, the medium was discarded, the cell culture flask was washed three times with PBS, and the cells were lifted with MSC digestion solution (Jing Meng, Beijing, China).
To confirm the expected surface marker phenotype, freshly-lifted hUC-MSCs were stained with fluorochrome-labelled mouse anti-human antibodies against CD73, CD44, CD29, CD105, CD90, CD45, and HLA-DR and with appropriate mouse isotype control antibodies using the BD Biosciences Human MSC Analysis Kit (catalog no. 562245, BD Biosciences, Franklin Lakes, NJ, USA) according to the manufacturer’s protocol. Flow cytometry was performed on a BD FACS Calibur (BD Biosciences) and the resulting data files were analyzed using FlowJo Software (Treestar, Ashland, OR).
To confirm tri-lineage differentiation capacity, hUC-MSCs were cultured in 6-well plates in adipogenic differentiation (Cyagen Biosciences, Guangzhou, China, catalog no. HUXUC-90031) or osteogenic differentiation media (Cyagen Biosciences, Guangzhou, China, catalog no. HUXUC-90021). For chondrogenic differentiation, hUC-MSCs were seeded at passage 4 in 15 ml sterile centrifuge tubes in chondrogenic differentiation medium (Cyagen Biosciences, Guangzhou, China, catalog no. HUXUC-90041). The induction of differentiation was confirmed by standard protocols from the manufacturer for adipogenesis (Oil red O staining), osteogenesis (Alizarin red staining) and chondrogenesis (Alcian blue staining) compared to hUC-MSCs cultured under baseline conditions at the same passage number served as negative controls.
Echocardiography: Two-dimensional, M-mode echocardiogram and tissue Doppler were performed using a Vevo 2100 ultrasound system under anesthesia with isoflurane at one day before euthanasia. In order to comprehensively evaluate the systolic and diastolic functions of the heart, the parasternal short-axis scan and the apical four-chamber scan were selected for assessment. Using the parasternal short-axis scan, the following parameters were recorded: heart rate (HR), left ventricular internal diameter at end-diastole (LVIDd), left ventricular internal diameter at end-systole (LVIDs), left ventricular anterior wall thickness at end-diastole (LVAWd), left ventricular anterior wall thickness at end-systole (LVAWs), left ventricular posterior wall thickness at end-diastole (LVPWd), left ventricular posterior wall thickness at end-systole (LVPWs), left ventricular weight (LVW), stroke volume (SV), cardiac output (CO), left ventricular ejection fraction (LVEF) and fractional shortening (FS). The first six indicators measured by ultrasound images can reflect the shape of the heart at the end of systole or diastole. The remaining indicators were calculated based on the first six indicators. LVEF and FS are the most typical indicators to evaluate cardiac systolic function [54]. By the apical four-chamber scan, we measured the ratio of peak early diastolic filling velocity to late atrial filling velocity (E/A) and isovolumic relaxation time (IVRT) which mainly reflect the diastolic function of the heart [55].
Tissue harvest, histopathology staining and image analysis
Following humane euthanasia, mouse hearts were isolated, weighed, fixed in paraformaldehyde and then embedded in paraffin. For each paraffin-embedded heart, whole-organ 4-µm thick sections were prepared using a microtome. For initial assessment of tissue quality and structure, a single heart section from one animal per group was stained with hematoxylin and eosin (H&E) and inspected by light microscopy. Next, for each experiment, whole heart sections from 3 different levels of the organ were stained with Masson’s trichrome (MT) for 6 mice per group. The protocols used for H&E and MT staining are provided in Supplemental Methods. For image analysis of MT-stained sections, 6 fields from each section were selected for the interstitial area, and 3 fields were selected for the perivascular area. Images were collected at 200×magnification and the MT-stained sections were analyzed in blinded fashion using a light microscope (Olympus BX53, Japan). Quantitative image analysis was performed using Image-Pro Plus 6.0 software (Image-pro Plus, Media Cybernetics, Inc., USA). For quantification of interstitial fibrosis, collagen volume fraction (CVF) was calculated for each field as the ratio of the blue-stained area to the total interstitial area. For quantification of perivascular fibrosis, the ratio of the blue-stained perivascular collagen area to the total luminal area (PVCA/LA) was calculated for each vessel within the field [56]. In total 90 fields per heart were analyzed for interstitial fibrosis and 45 fields per heart were analyzed for perivascular fibrosis. The final CVF and PVCA/LA values, expressed as %, were derived for each heart from the average values of all analyzed fields.
Selection of candidate differentially expressed miRNAs and mRNAs from publically available sources: In order to identify miRNAs with altered expression in diabetic mouse heart, we selected Non-DM and DM group sequencing data from a miRNA high-throughput sequencing (series GSE210036 [15]) and generated a heatmap of the top 50 differentially-expressed miRNAs. From this, we selected six fibrosis-associated miRNAs: let-7f, miR-26a, miR-29a, miR-29b, miR-29c and miR-133a. Based on the recent review by Z-Q Jin [57] we also selected three additional candidate miRNAs: miR-34a, miR-155 and miR-326. The relative expression of these nine miRNAs, and subsequently, the relative expression of candidate mRNAs in heart tissue from the experimental groups of the current study were then determined by qRT-PCR (methods described below).
To identify candidate mRNAs with altered expression in the heart in diabetes, we also selected publically available Non-DM and DM group data from a recently performed mRNA high-throughput sequencing analysis (series GSE161052 [17]) and generated a volcano plot with annotation of the top 18 differentially expressed mRNAs. We noted that these mRNAs included transcripts encoding the alpha 1 chains of collagen I and collagen III and selected these mRNAs for qRT-PCR analysis in hearts from the current study due to their relevance to fibrosis and prior identification of COL1A1 as a target for miRNA-133a. In addition, 6 other fibrosis-related mRNAs potentially regulated by miRNA-133a - α-SMA, TGF-β, Smad2, Smad3, Smad4 and FGF1 were also selected for quantitative analysis in heart tissue from our experimental groups.
RNA extractions and quantitative reverse transcription and polymerase chain reaction (qRT-PCR) analyses
After euthanasia, one-half of the fresh heart tissue was stored in a -80 ° C freezer for further use (n = 6 per group). For qRT-PCR analyses, total miRNA or mRNA were extracted from heart tissue using miRNA extraction and isolation kit (DP501, TIANGEN, Beijing) and mRNA extraction and isolation kit (DP424, TIANGEN, Beijing) respectively according to the manufacturers’ instructions. The miRNA and mRNA preparations were transcribed to cDNA by miRNA First Strand cDNA Synthesis (Tailing Reaction, B532451, Sango Biotech, Shanghai) and cDNA first strand synthesis kit (ZS-M14003, ZHONGSHI TONTAU, Tianjin), respectively according to the manufacturers’ instructions. Quantitative PCR was performed with a SYBR Green PCR master mix (ZS-M13010, ZHONGSHI TONTAU, Tianjin). Because miRNA was generally composed of more than 20 bases, different from general mRNA, we adopted the poly(A) tailing method for reverse transcription of miRNA. The tailed reverse transcription assay utilizes poly (A) polymerase for mature miRNA plus poly (A) tail, followed by 5 '- end universal tagged oligo dt as the reverse transcription primer to obtain the human elongated miRNA cDNA first strand, and finally fluorescence quantitative PCR detection with a reverse primer complementary to the universal tag sequence by dye method or probe method. Universal PCR primer R (B661601-0002, Sango Biotech, Shanghai) were the same for all miRNAs and internal control U6. The primer F sequences of miRNA are listed in Table S2 and primer sequences of mRNA are listed in Supplemental Table S3. A melt curve was included to ensure primer specificity. Experiments were performed in triplicate, and results were normalized to U6 or β-actin expression (2−ΔΔCT method).
Statistical analysis
For all experiments, group data were expressed as mean ± SD and analyzed using GraphPad Prism 6. The differences between the groups were compared using ordinary one-way analysis of variance (ANOVA) or multiple comparisons of 2-way ANOVA. P < 0.05 was considered statistically significant. For bioinformatics analyses of publically available miRNA and mRNA datasets, R programming language (version 4.1.3) and R studio (version 3.3.0) were used.