Animal experiment
Six macaque monkeys were obtained from Kunming Institute of Zoology, Chinese Academy of Sciences [SCXK (Yunnan) K2017-0003]. The macaque monkeys were aged 22-26 years old and were randomly divided into two groups: the senescence model group (n=3) and the senescence treatment group (n=3). P4 generation UCSSCs were diluted with 0.9% sterile sodium chloride solution at a concentration of 2 × 106 cells/ml. After the macaque monkeys in the treatment group were immobilized, UCSSCs were injected into the femoral vein at 1 × 107/kg per monkey once every other day 3 times. The macaque monkeys in the elderly model group were given the same amount of 0.9% normal saline. Experimental animal certificate number: SYXK (military) 2017-0051. Ethical approval for the use of animal models was provided by the Experimental Animal Committee of the 920th Hospital of the Joint Logistics Support Force.
Thymus index measurement
After three months of treatment, macaque monkeys were anaesthetized with more than 3% pentobarbital sodium and euthanized. Then, the following steps were followed: the chest cavity was opened, the heart was located, then the thymus was identified along the great blood vessels, removed, and weighed (g) with an electronic balance; finally, the thymus organ index of rhesus monkeys after treatment was determined.
Thymus tissue HE staining
The thymus tissue was fixed in 4% paraformaldehyde solution for 24 hours and then dehydrated and embedded in paraffin. The thickness of the sections was approximately 4μm. After being treated with xylene-anhydrous ethanol-75% ethanol, the sections were sealed with neutral glue. Sections were observed and photographed with a Nikon microscope (NikonDS-U3, Nikon, Japan).
Determination of thymosin α and thymopoietin levels
After three months of treatment, the serum was isolated from the peripheral blood of macaque monkeys, and the expression levels of thymosin α and thymopoietin in serum were measured by ELISA kits (Millipore).
Detection of P16 expression by immunofluorescence staining
The thymus tissue was fixed in 4% paraformaldehyde solution for 24 hours and then dehydrated and embedded in paraffin. The thickness of the section was approximately 4μm. An anti-rabbit anti-monkey antibody P16 was added for staining. FITC goat anti-rabbit was incubated at room temperature for 50 minutes, and the nucleus was restained with DAPI for 15 minutes. Finally, the sections were sealed with anti-fluorescence quenching sealant and photographed under a microscope.
Immunohistochemical staining of the thymic epithelial cell-specific marker molecule CK14
Thymus tissue section with a thickness of 4 μm were dewaxed in water for antigen repair and blocking of endogenous peroxidase. One anti-mouse anti-monkey CK14 antibody was used, and the second antibody was labelled with HRP. After incubation with the secondary antibody, the nuclei were stained with DAPI. Sections were observed and photographed with a Nikon microscope (Nikon DS-U3, Nikon, Japan). The positive area signal was analysed by ImageJ software (NIH, Bethesda, MD).
TECs
TECs were purchased from Hangzhou Meisen Co. and highly expressed the epithelial cell-specific marker EPCAM. They were cultured in a medium (DMEM/F12) containing 10% calf serum (Invitrogen).
Establishment and identification of the TEC ageing model
Fourth-generation TECs were collected, the cell concentration was adjusted to 5 × 105 L-1, and 100 μL of cell suspension was inoculated in a 6-well culture plate. When the cell fusion rate was approximately 80%, 200 µmol/LH2O2 (China Nanguo Pharmaceutical Co., Ltd.) was added and cocultured for 72 hours.
Molecular analysis of ageing markers
β-Galactosidase (SA- β-Gal) staining TECs were collected, washed with PBS once, operated according to the kit instructions (Solarbio, G1580), and observed under an ordinary optical microscope (ZEN2011), and the number of blue-positive cells in three visual fields was counted.
Expression level analysis of P21 and P16 TECs were collected for Western blotting analysis and lysed in RIPA buffer (Solarbio, China, Beijing). The total proteins of TEC cells were separated by Tris-Gly and 10% gel electrophoresis and then transferred to a nitrocellulose membrane. The cell membrane was sealed with 5% skim milk at room temperature for 1.5 hours and then incubated with an antibody overnight. The next day, the membrane was washed with PBST and incubated for 30 min with horseradish peroxidase-labelled secondary antibody. Protein bands were detected after treatment with SuperSignal West Femto (Thermo Fisher Science, Waltham, MA, USA). The primary anti-rabbit anti-human P16 (1RV 1000, ab32503), P21 (1RV 1000, ab32351), SRC (1RV 300000 ab179473 ab109381), and ITGB3 (1RV 1000 ab179473) antibodies were obtained from Abcam. The primary anti-rabbit anti-human β-actin antibody (1RV 25000~66009-1-Ig) and the secondary anti-HRP-Goat IgG antibody (1RV 1000, SA00001-2) were purchased from Proteintech.
Observation of the cell substructure by fluoroscopic electron microscopy
TECs were collected, digested with 0.25% trypsin and made into a cell suspension. After washing with PBS, the cell precipitate was centrifuged and fixed for 4 hours. After immersion in acetone, ultrathin sections of 80 nm were stained with uranium and lead, observed under a transmission electron microscope (Hitachi, HT7700), and photographed for image analysis.
Establishment and identification of a coculture model of UCSSCs and senescent TECs
Fourth-generation TECs were collected, the cell concentration was adjusted to 5 × 105 L-1, and 100 μL of cell suspension was inoculated in a 6-well culture plate. When the cell fusion rate was approximately 80%, 200 μmol/L H2O2 was added and treated for 72 h, and the supernatant was discarded. Fourth-generation UCSSCs were collected, the cell concentration was adjusted to 5 × 105 L-1, a 100μL cell suspension was inoculated into the Transwell chamber, the Transwell chamber was embedded on the upper layer of the 6-well culture plate, and DMEM/F12 containing 10% foetal bovine serum was used for coculture of cells for 48 hours.
Analysis of ageing-related indexes of TECs after coculture: After coculture, the senescence marker molecules of TECs were analysed, and the cell substructure was observed by fluoroscopic electron microscopy.
Continuous observation of the morphology of cocultured TECs by a Nanolive 4D microscope
P4 generation TECs at a cell concentration of 5000/ml were added to a Nanolive 4D microscope petri dish, and then 1.5 ml DMEM/F12 (containing 10% foetal bovine serum) medium was added and mixed slowly; the cells were then cultured for 24 h. The culture medium was discarded, and the culture medium was discarded after the TECs were treated with 200 μmol/L H2O2 for 72 h. Next, the cell concentration was adjusted to 5000 cells/ml of P4 generation UCSSCs, the Nanolive 4D microscope was started and connected to the computer, the focal length was adjusted to take pictures under the appropriate field of vision for 48 hours, and pictures were taken every 4 min.
DNA methylation and transcriptome sequencing
The TECs of UCSSCs were 1.3x107 in coculture and 1.2x107 in the model group. The cell precipitates were frozen in liquid nitrogen and stored at -80 °C for 24 hours. The whole-genome methylation and transcriptome data were obtained by whole-genome methylation and transcriptome sequencing at Guangzhou Jidiao Biotechnology Co., Ltd.
Analysis of DNA methylation sequencing data
The obtained clean data were mapped to the species reference genome using BSMAP software (version 2.90). Methylated cytosine level was then detected using a custom Perl script and assessed using the correction algorithm described by Lister R, et al. The level of methylation was calculated based on the percentage of methylation of each sequence (CG, CHG, and CHH) in the entire genome, each chromosome, and different regions of the genome. To explain the changes in DNA methylation in TECs, methylation maps of the upstream 1 kb region and gene body were drawn according to the average methylation level of each window. Changes in GC content were screened according to the number of differentially methylated regions (DMRs) (≥ 5) and the absolute value of methylation ratio difference (≥ 0.1). The senescence genes were counted according to GeneOntologyConsortium.
RNA-seq analysis
The transcriptome database was analysed by DESeq2 software, and the mRNAs satisfying the conditions of log2FC > 1.2 and P value < 0.05 were regarded as differentially expressed genes. The differentially expressed genes in ageing TECs and ageing TECs after coculture were counted and analysed by gene function annotation.
Analysis of the relationship between gene methylation and transcriptional level
To calculate the correlation between gene methylation and transcription, the genes were divided into four groups according to the transcription level in RNA-seq data (FPKM ≤ 1 for no transcription, 1 < FPKM ≤ 10 for low transcription, 10 < FPKM ≤ 100 for medium transcription and FPKM > 100 for high transcription). The Pearson correlation coefficients between upstream/gene body methylation levels and corresponding gene expression levels in TECs of the control, model and coculture groups were calculated.
Establishment of the transcription factor-methylation regulatory network
To explore the relationship between specific gene methylation and transcription factor binding, the TFs of ageing-related genes were predicted by CistromeDB, and the TF-DNA regulatory network was constructed by Cytoscape.
Enrichment analysis
Through the Gene Ontology database (GO item in http://www.geneontology.org/)), the number of genes in each item was calculated, and the GO items significantly enriched relative to the genomic background were defined by hypergeometric test. KEGG is the main public pathway-related database (http://www.kegg.jp/kegg/)).
Functional verification
Expression analysis of RAP1 signal pathway molecules RAP1, SRC and ITGB3 For Western blotting analysis, the primary anti-rabbit anti-human RAP1 (1RV 1000, ab14404), SRC (1RV 30000 ab179473 ab109381) and ITGB3 (1RV 1000 ab179473) were all obtained from Abcam, as described above.
Analysis of the relative expression of the ageing regulatory genes NGF, KRT17 and FOXJ1 Total RNA was extracted according to the production manual (Carlsbad Invitrogen Company, California) by RT-QPCR using TRIzol reagent. cDNA was amplified by reverse transcriptase (GoldenstarTMRT6, China) and random primers. According to the SYBR Green Master Mix specification (TsingKe Biotech Co., China), a CFX96TM Real-Time System (Applied 2720) was used for qPCR. The glyceraldehyde 3-phosphate dehydrogenase (GAPDH) gene was used as an internal control. Each experiment was repeated 3 times, and the gene expression level was normalized by the 2-T square method. All qPCR primers are listed in( Table 1).
Immunofluorescence detection of NGF and KI67 expression TECs were collected and fixed with paraformaldehyde, and the fixation solution was removed. Anti-mouse anti-human antibodies against Ki67 and NGF (1R1000 ~ (th) ServiceBio) were added for staining. The nuclei at 15 min were restrained with CY3 goat anti-rabbit DAPI at room temperature for 50 minutes and then sealed with anti-fluorescence quenching sealant. Microscopic examination: The sections were observed under a fluorescence microscope, and images were collected (Nikon DS-U3, Nikon, Japan; DAPI UV excitation wavelength 380 nm, emission wavelength 420 nm, blue light; CY3 excitation wavelength 510-560 nm, emission wavelength 590 nm, red light). The number of positive cells was analysed by ImageJ software (NIH, Bethesda, MD).
Statistical analysis
SPSS 21.0 statistical software was used for analysis, and the data are expressed as the mean ± standard deviation. A t test was used to compare the two groups, and one-way analysis of variance (one-way ANOVA) was used to analyse the significance of differences among the three groups.