Research Paper Volume 13, Issue 8 pp 11218—11233
Exosomal circRNAs contribute to intestinal development via the VEGF signalling pathway in human term and preterm colostrum
- 1 Department of Pediatrics, Women’s Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital, Nanjing 210004, China
- 2 Department of Pediatrics, Fourth Clinical Medicine College, Nanjing Medical University, Nanjing 210029, China
Received: June 16, 2020 Accepted: February 8, 2021 Published: April 4, 2021
https://doi.org/10.18632/aging.202806How to Cite
Copyright: © 2021 Zhou et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Abstract
Human breast milk (HBM) provides essential nutrients for newborn growth and development, and contains a variety of biologically active ingredients that can affect gastrointestinal tract and immune system development in breastfed infants. HBM also contains mRNAs, microRNAs and lncRNAs, most of which are encapsulated in milk-derived exosomes and exhibit various important infant development related biological functions. While previous studies have shown that exosomal circRNAs are involved in the intestinal epithelial cells’ proliferation and repair. However, the effect of HBM exosomal circRNAs on intestinal development is not clear. In this study, we identified 6756 circRNAs both in preterm colostrum (PC) and term colostrum (TC), of which 66 were upregulated, and 42 were downregulated (|fold change>2|, p < 0.05) in PC. Pathway analysis showed that the VEGF signalling pathway was involved, and network analysis revealed that the differentially expressed circRNAs bound various miRNAs. Further analyses showed that has_circRNA_405708 and has_circRNA_104707 were involved in the VEGF signalling pathway, and that they all bound various mirRNAs. Exosomes found in preterm colostrum (PC) and term colostrum (TC) promoted VEGF protein expression and induced the proliferation and migration of small intestinal epithelial cells (FHCs). Exosomal circRNAs found in human colostrum (HC) binding to related miRNAs may regulate VEGF signalling, and intestinal development.
Introduction
Human breast milk (HBM), contains nutrients needed by infants as well as various non-nutritive bioactive ingredients that can promote growth and development [1]. Preterm milk (PM) is also known to contain more proteins and lipids than term milk (TM) [2, 3], and while there are many studies on the differences in macronutrients between PM and TM, there have only a few reports comparing exosomal content.
The immaturity of the intestines of premature babies means that they are particularly vulnerable to damage. Multiple studies have shown that breastfeeding significantly reduced necrotizing enterocolitis (NEC) [4, 5], and that breast milk (BM) can increase intestinal permeability, downregulate oxidative stress [6], and directly decrease damage intestinal epithelial cells [7].
Milk exosomes have important biological functions [8]. A previous study has shown that HBM exosomes reduce intestinal epithelial cell death [9] and promote intestinal epithelial cells growth [10]. There is emerging evidence that exosomal RNAs are involved in many processes, including cellular growth, antiviral activity, and insulin sensitivity [11].
CircRNAs were initially considered to be the product of abnormal RNA splicing [12], and many circRNAs have been discovered in recent years [13]. circRNAs have been shown to promote intestinal stem cell self-renewal [14], which is related to inflammatory bowel disease and intestinal barrier formation [15]. Recent studies have shown that exosomal circRNAs play critical roles in the development of neonatal tissue and organ, such as the brain [16] and the nervous system [17]. However, the expression and function of HBM exosomal circRNAs in intestinal development still remains unclear.
The aim of this study was to evaluate the differential expression of exosomal circRNAs in PC and TC and analyse key circRNAs that may regulate intestinal development. Differential expression of circRNAs between PC and TC groups was detected by microarray analysis, and 6756 exosomal circRNAs were identified. GO, pathway analyses and analysis of downstream miRNAs were used to identify and predict circRNAs functions. Our study may help to elucidate the role of HC in intestinal development.
Materials and Methods
Ethics statement
HC was obtained from the Woman’s Hospital of Nanjing Medical University. This study was approved by the Institutional Review Board at the Women’s Hospital of Nanjing Medical University, Nanjing Maternity, and Child Health Care Hospital [Permission Number (2013)78].
Human colostrum (HC) collection
Colostrum was collected from lactating mothers (n=18) who had donated their milk at the Woman’s Hospital of Nanjing Medical University. Nine had a full-term birth (≥37 weeks), and nine had preterm births (24-36 weeks). Colostrum was selected as our experimental sample material to reduce any differences that may be seen in human breast milk expressed at different times post birth, and to account for any differences between individual donors, each sample contained colostrum from there different donors. HBM was then obtained from these samples. The hospital provided the mother’s food, minimising the chance of food affecting colostrum composing between samples. Detailed data are shown in Table 1.
Table 1. Data related to the samples.
| Group | Sample | Mother | Age | Week | Weight | Delivery way | Parity |
| 1 | 30 | 40+1 | 3430 | natural birth | G2P2 | ||
| T-1 | 2 | 27 | 39+1 | 3180 | natural birth | G1P1 | |
| 3 | 26 | 38+2 | 3310 | natural birth | G1P1 | ||
| 4 | 30 | 39 | 3600 | Caesarean section | G2P2 | ||
| Term | T-2 | 5 | 37 | 40+1 | 3640 | natural birth | G2P1 |
| 6 | 30 | 40+3 | 3590 | natural birth | G2P2 | ||
| 7 | 25 | 41+3 | 3740 | natural birth | G2P1 | ||
| T-3 | 8 | 28 | 41+1 | 3550 | natural birth | G1P1 | |
| 9 | 26 | 39+3 | 2800 | natural birth | G1P1 | ||
| P-1 | 10 | 24 | 33 | 1870 | Caesarean section | G1P1 | |
| 11 | 25 | 35+2 | 2730 | natural birth | G1P1 | ||
| 12 | 29 | 32 | 2120 | natural birth | G1P1 | ||
| 13 | 31 | 34+3 | 2480 | natural birth | G2P1 | ||
| Preterm | P-2 | 14 | 26 | 35+4 | 2460 | natural birth | G1P1 |
| 15 | 30 | 33+4 | 1870 | natural birth | G2P1 | ||
| 16 | 25 | 33+4 | 2290 | natural birth | G3P1 | ||
| P-3 | 17 | 34 | 32+5 | 1830 | Caesarean section | G3P2 | |
| 18 | 27 | 33+6 | 2790 | natural birth | G4P2 |
Isolation and purification of human colostrum exosomes
Exosomes were isolated immediately after colostrum samples were obtained from both breasts of each donor. The isolation method used in this study has been previously described [18]. In brief, 50 ml of HC was centrifuged twice at 3000 g for 15 min at 4° C to remove cells and fat globules. The supernatant was then transferred to new tubes, filtered through a 0.22 μm filter to remove any remaining debris and centrifuged at 42000 rpm at 4° C for 120 min. Exosomal pellets were stored at -80° C until use. The exosomes were observed with transmission electron microscopy (TEM) (USA, FEI Company, FEI Tecnai T20), and their quantity and size were assessed using the Nano Sight NS (Malvern, UK, 300 system, Nano Sight Technology). The collection time for all samples was within one month.
Exosome labelling for in vitro studies
HC-Exo were labelled with PKH26 according to the manufacturer's protocol. The standard solution was passed through a 0.22μm filter (American microporous, Billerica's) and washed three times to remove excess dye. FHC cells were seeded in a 6-well plate for 6 h, then treated with PKH26 marker solution.
Extraction and reverse transcription of total RNA
Total RNA was extracted from the samples using TRIzol Reagent (Invitrogen, Carlsbad CA, USA), and a pure tissue kit (Tiangen, DP431) was used to isolate the RNA. A 10 μM reverse transcription reaction was performed using 1000 ng of total RNA. The random primer reverse transcription method was used with the ScriptTM RT Master Mix kit, according to manufactures’ instructions, to generate cDNA synthesis for each sample.
CircRNA extraction and microarray analysis
Arraystar Human circRNA Array v2 analysis was performed on the samples (Kangchen, Shanghai, China), and a NanoDrop ND-1000 was used to quantify otal RNA in each sample. The standard l Arraystar-based protocol was used for sample preparation and microarray hybridization. Briefly, RNase R (Epicentre, Inc.) was used to digest total RNA, remove linear RNAs, and enrich circular RNAs. We amplified and transcribed the enriched circRNAs into the fluorescent cRNA (Arraystar Super RNA Labelling Kit; Arraystar) and hybridized labelled cRNAs on the Arraystar Human circRNA Array v2 (8x15K, Arraystar). After washing, the array was scanned using an Agilent G2505C scanner.
The array image obtained by Agilent Feature Extraction software (version 11.0.1.1) was analyzed. Quantile normalization and subsequent data processing were performed using the LIMMA (implemented with the R software package). Volcanic maps were used to identify significant differences in the circRNA expression between the two groups and the hierarchical clustering method shows the distinguishable patterns of circRNA expression between samples. |log2 fold change (FC)| ≥ 2 and a P value < 0.05 by t test were considered to be statistically significant.
qRT-PCR verification of target genes
qRT-PCR was used to evaluate the results of SYBR. The 2-ΔΔCT method was used to analyse the experimental data. All data were averaged from three separate experiments. GAPDH was the internal reference for the target genes. Primer sequences are shown in Table 2.
Table 2. Primer sequences of circRNAs.
| Primer name | Sequence |
| 004239-F | ACCCAACAACCTGGTCCATA |
| 004239-R | TGTTGGGGACCTTGTTCCTA |
| 004239-F | ACCCAACAACCTGGTCCATA |
| 004239-R | TGTTGGGGACCTTGTTCCTA |
| 101018-F | GGCACAGTGAGACAGATGCT |
| 101018-R | CATCATTTACAGTTTTTCCTGGTG |
| 101018-F | GGCACAGTGAGACAGATGCT |
| 101018-R | CATCATTTACAGTTTTTCCTGGTG |
| 103356-F | GCATGATGCCTCATCAACAG |
| 103356-R | GTCAAGTTCCTCCGACAAGC |
| GAPDH-F | GAAGGTGAAGGTCGGAGTC |
| GAPDH-R | GAAGATGGTGATGGGATTTC |
GO and KEGG pathway analysis
CircBase (http://www.circbase.org) was used to retrieve the encoded genes and predict their target genes. DAVID (https://david.ncifcrf.gov) was used for the target gene analysis.
Western blotting analysis
Proteins were extracted for detection using the radioimmunoprecipitation (RIPA) test kit (Sigma-Aldrich, USA) according to the manufacturer's instructions and protein concentration measured using the bicinchoninic acid (BCA) toolkit (Pierce, USA). and SDS-containing polyacrylamide gels (SDS-PAGE) were used to separate equal amounts of protein Samples (30 μg for exosome pellets). These samples were then transferred to a polyvinylidene fluoride (PVDF) membrane (Bio-Rad, USA). Then, 1 x TBST (0.1 M, pH 7.4) containing 5% skim milk was used to block the membrane for 1 h, hen hybridized with CD9 (Proteintech, 1:1000 dilution), CD63 (Proteintech, 1:1000 dilution), CD81 (Proteintech, 1:1000 dilution), PI3K/P-PI3K (Proteintech, 1:1000 dilution), AKT/P-AKT (Proteintech, 1:1000 dilution) and VEGF (Proteintech, 1:1000 dilution) antibodies. The immune complex was incubated with a secondary antibody conjugated to horseradish peroxidase (Applygen, China; 1:5000 dilution) for 1 h at room temperature. Protein bands were examined and the signal intensities of interest were calculated using ImageJ software (National Institutes of Health, USA). The results are presented as the relative fold change ± SD.
Cell proliferation assay
Cell proliferation was assessed using the Cell Counting Kit-8 (CCK-8) assay (Dojindo, Tokyo, Japan) according to the manufacturer's instructions. FHC cells (X-Y Biotechnology) were grown in 96-well plates with a cell density of approximately 1x103 cells per well then co-cultured with TM or PM exosomes (20.0 μg of protein) for 0, 12, 24, or 36 h. The cells were measured at 450 nm (optical density (OD) value) using a microplate reader (BioTek Instruments, Inc., Germany).
Wound closure assay
FHC cells were plated into a 6-well plates and scratched with the tip of a sterile pipette. TM and PM exosomes were added (200 μg). The scratched areas were imaged under a microscope at 0 and 12 h, and ZEN software (Zeiss, Germany) was used to measure the scratches widths.
Statistical analysis
Data were analysed by SPSS, and Student’s t-test or one-way ANOVA was employed for statistical comparisons. P < 0.05 was statistically significant.
Results
Characteristics of the HC exosomes
Exosomes were isolated from human TC and PC using ultracentrifugation and transmission electron microscopy (TEM) was used to investigate them at the nanometre scale (Figure 1A). CD63, CD9 and CD81 are exosome surface markers and all were expressed in PM. CD81 and CD9 were also expressed in the TM exosomes, although CD81 expression was lower in TC exosomes compared to PC exosomes (Figure 1B). Nanosight (NTA) analysis showed that TC and PC exosomes were both 40-50 nm in diameter (Figure 1C) and there was a significantly higher abundance of exosomes in PM compared to TM (Figure 1D). These results suggest that exosomes were have been successfully isolated from the TC and the PC, and that the exosome levels may be higher in the PC than TC.
Figure 1. HC Exosome characterization (A) TEM observation of breast milk exosomes showed a typical cup-shaped structure. (B) WB verification of presence of CD9, CD81 and CD63. (C) NTA analysis of preterm and term milk exosomes diameters. (D) Comparison of preterm and term breast milk exosomes abundance.
Differentially expressed circRNAs in HC exosomes
Exosomal circRNAs from TC and the PC were identified through microarray analysis. A total of 6756 circRNAs were detected and 2845 were downregulated and 3911 circRNAs were upregulated in PC. These differentially expressed circRNAs were distributed across all chromosomes, including sex chromosomes (Figure 2A). Detailed information on the differentially expressed circRNAs is shown in Supplementary Table 1. Detected circRNAs are also shown by scatter plots and volcano maps (Figure 2B, 2C), and a heat map illustrates the differences in circRNAs, 66 of which were upregulated and 42 were downregulated (fold change (FD) > 2, P < 0.05) (Figure 2D).
Figure 2. Differentially expressed circRNAs in HC exosomes. (A) A total of 6756 circRNAs were detected, 3911 were upregulated and 2845 were downregulated. Histogram showing the differentially expressed circRNAs in the human chromosome. (B) The X-axis and Y-axis values are normalized signal values (log2 scaling). The green line represents a broken line. Above and below the upper and lower green lines, the circRNAs exhibited more than 2.0-fold changes between the preterm milk (PM) and term milk (TM) groups. (C) The red dots in the volcano map represent statistically significant differentially expressed circRNAs. (D) Heat map of circRNA results from the PC and TC groups. Red bars indicate higher expression, and green bars indicate lower expression.
Verification of circRNAs by qRT-PCR
To verify the exosomal circRNAs which were significantly expressed in the PC and TC, we randomly selected 6 differentially expressed circRNAs for qRT-PCR verification in the same samples used for the HC microarray analysis; 3 upregulated (hsa_circRNA_104423, hsa_circRNA_101018, hsa_circRNA_044097) and 3 downregulated (hsa_circRNA_004239, hsa_circRNA_084900, hsa_circRNA_103356) circRNAs were selected (Figure 3A, 3B). The results showed that verified circRNAs expression was consistent with the microarray analysis results.
Figure 3. circRNAs expression verification by qPCR. *P-value < 0.05, **P < 0.01 ***P < 0.001. (A) The circRNAs upregulated in PC, (B) The circRNAs downregulated in PC.
GO enrichment and KEGG signal pathway analysis
GO enrichment and signal pathway analysis were performed to assess potential functions. Biological process analysis showed that the circRNAs were mostly involved in organelle localization, post-translational protein modification, glycoprotein metabolic process, immune response-activating signal transduction translation, positive regulation of immune response, carbohydrate metabolic process, apoptotic signalling pathway, and protein complex assembly (Figure 4A). The cellular component analysis demonstrated that the potential functions were significantly related to focal adhesion, cell-substrate junction, Golgi apparatus, cytosol, cytoskeleton, membrane-bound vesicle, vesicle and extracellular vesicular exosomes (Figure 4B). Molecular function indicated that the most significant functions were protein kinase binding, cytoskeletal protein binding, enzyme regulator activity, ribonucleoside binding, nucleoside binding, poly(A) RNA binding, ATP binding purine ribonucleoside triphosphate binding, adenyl nucleotide binding and transferase activity (Figure 4C). Analysis of the KEGG pathway showed that the VEGF signalling pathway was the most significant pathway.
Figure 4. Top 10 GO terms with gene count and enrichment analysis for biological processes. (A) Biological process, (B) Cellular component, (C) Molecular function. (D) The top 10 enriched KEGG pathways of the parental genes of the differentially expressed circRNAs.
Detailed annotation of the interaction between circRNAs and miRNAs
Previous studies found that circRNAs perform biological functions through their downstream miRNAs. To evaluate the target miRNAs, Target Scan and miRanda databases were used to make theoretical predictions based on conservative seed matching sequences. 108 circRNAs (FD>2.0, P<0.05) were identified that could bind to miRNAs (Supplementary Table 2). The relationship between the circRNAs and miRNAs is shown as a network (top 20 differentially expressed circRNAs) (Figure 5A). The KEGG signalling pathway analysis showed that the VEGF signalling pathway was the most significantly enriched. Further analysis identified two circRNAs involved in the VEGF signalling pathway, one of which was upregulated and the other was downregulated. They bound to the miRNAs hsa_miR-4481, hsa_miR-6875-3p, hsa_miR-4745-5p, hsa_miR-6827-38, hsa_miR-4722-5p, hsa_miR-301a-5p, hsa_miR-125a-3p, hsa_miR-134-5p, hsa_miR-222-5p and hsa_miR-501-5p (Figure 5B).
Figure 5. Network of the relationship between circRNAs and miRNAs. (A) Red dots represent upregulated circRNAs in PC, blue dots represent downregulated circRNAs in PC, and yellow dots represent downstream miRNAs that bind to circRNAs; (B) Flow chart of the experiment and two circRNAs involved in the VEGF signalling pathway, one of which was upregulated and the other was downregulated.
HC exosomes function in small intestinal epithelial cells (FHCs) is through the VEGF signalling pathway
The KEGG signalling pathway analysis showed that the VEGF signalling pathway was the most significantly enriched. Further functional experiments have shown that HC-Exo, taken up by FHC cells after PKH26 labelling for 6h, could promote proliferation and migration of intestinal epithelial cells (FHC) (Figure 6A–6C) and that the PM exosomes were significantly more potent than the TM exosomes. Experiments have also shown that PC and TC exosomes promoted the VEGF protein expression (Figure 6C) and phosphorylation of PI3K and AKT in FHC cells, which are the key downstream regulators of VEGF signalling. PI3K phosphorylation was higher with PC than TC exosomes (Figure 6D).
Figure 6. HC exosomes function through the VEGF signalling pathway in small intestinal epithelial cells (FHCs). (A) PKH26 labelled exosomes taken up into the cells, (B) Proliferation of intestinal epithelial cells (FHC); (C) the migration of intestinal epithelial cells (FHC); (D) the VEGF protein expression; (E) P-PI3K/PI3K and P-AKT/AKT protein expression.
VEGFR inhibition blocked the proliferation and migration of FHC cells induced by PC and TC exosomes
AV951 (a VEGFR inhibitor) was added to suppress the VEGR signalling to further verify the HC-Exo involvement in the pathway. It was shown that inhibition of VEGF signalling significantly decreased the proliferation and migration of FHC cells promoted by TC and PC exosomes (Figure 7A, 7B).
Figure 7. VEGFR inhibition blocks PC and TC exosome-induced FHC proliferation and migration. (A) TM and PM exosomes-induced proliferation, and (B) migration of FHC cells after VEGFR inhibitor intervention.
Discussion
Because it contains many active substances, which are particularly important for the growth and development of newborns, HBM quality is especially important for premature infants [1, 19]. PM and TM compositions are significantly different, but they had significant biological function [3]. To an extent, milk's nutritional content is related to the mother's diet [20]. HBM used in this study was from the colostrum of the mothers of newborns During the collection period, our hospital food department provided the mothers' food, which reduced the HBM content variation to some extent. Previous studies have shown that breast milk reduces the risk of NEC [7]. The abundant bioactive factors in HBM are essential for the development of the innate immune system of the gastrointestinal tract [21], and a number of studies have shown that BM exosomes may affect intestinal development. For example, intestinal epithelial cells growth in infants can be promoted by BM exosomes. [10, 22]. Inhibition of the TLR4/NF-κB and p53 pathways attenuated LPS-induced intestinal epithelial cell apoptosis [23] and promoted cell proliferation [24]. In this study, TEM, Nano-sight, and surface markers (CD63, CD9 and CD81) analysis has shown that we successfully extracted HC-Exo. It was observed that CD63 and CD81 levels were lower in TC exosomes than in PC, which may be due to there being fewer exosomes in TC than in PC.
Exosomes contain many active ingredients (including miRNAs, lncRNAs and mRNAs) [25, 26], as well as new types of long noncoding RNAs (circRNAs) [13]. Recent studies have suggested that exosomal miRNAs extracted from milk promote commitment of Treg cells commitment by decreasing the immune response of Il-4/Th2-mediated allergic sensitization and that milk transfer of exogenous miR-155 may induce thymic regulatory T cells to secrete key components required for immune regulation and epigenetic modification [27]. Intestinal immune system maturation is crucial for intestine developing, and these RNAs may strengthen the immune system and, in turn, promote development. The roles of exosomal circRNAs involved in intestinal development remain unclear. In this study, we isolated circRNAs from PC and TC exosomes and analysed them by high-throughput microarray screening.
We performed GO analysis and found that differentially expressed circRNAs in PC and the TC exosomes significantly participated in immune response, activation of signal transduction, translation, and positive regulation of the immune response. The immune response and immune activation play roles in newborns development [28], and by improving the intestinal immune system of newborns, help to promote resistance to bacterial infection [29, 30]. The in vitro synthesis of circRNAs activates the innate immune response [31], which protects against viral infection. In addition, the immune factors NF90/NF110 regulate circRNA biosynthesis and inhibit viral infection by interacting with viral miRNAs [31, 32]. These findings indicate that exosomal circRNAs from HC may be critical in the development of the immune system of newborns and thus promote intestinal development and help defend against neonatal infection.
The VEGF signalling pathway regulates angiogenesis of various organs and tissues and alleviates adaptive cerebral ischaemia/reperfusion injury in rats [33]. Exosomal MMP2 from mature osteoblasts promotes endothelial cell angiogenesis through the VEGF signalling pathway [34]. The intestinal microvascular system is critical in NEC [35]. In neonatal mice, loss of VEGFR2 signalling led to dysplasia of the intestinal microvascular system and necrotizing enterocolitis in neonatal mice. Inhibition of VEGFR2 kinase activity also reduced villous endothelial cells proliferation in newborns mice and increased morbidity and mortality from NEC [36]. In this study, we found that the primary enrichment pathway of exosomal circRNAs of the PC and the TC was the VEGF signalling pathway. We also showed in FHCs, PC and TC exosomes can promote cell proliferation and migration, and that VEGF protein expression was also increased. PC and TC exosomal circRNAs may therefore promote intestinal epithelial cell proliferation and migration through the VEGF signalling pathway.
Recent studies have also shown that circRNAs can bind miRNA and regulate downstream miRNAs, all of which are involved in various pathological mechanisms. For example, studies conducted by Eur Heart J et al. showed that the cardiac circRNA HRCR could bind directly binds to miR-223 and inhibits HRCR expression in human cardiomyocytes, producing a hypertrophic response [37], and ciRS-7 (human circRNA cerebellar degeneration-related protein 1 transcript (CDR1)) binding to related miRNAs can affect the availability of miR-7 to bind to its target mRNAs [38–40]. CircSRY also acts as a sponge of miR-138 and as a platform for the binding of miR-138, regulating the mRNA translation of its targets [41]. Sponging of 9 miRNAs, including miR-124, which regulates cell growth, has been previously reported [42]. In our study, we predicted that exosomal circRNAs derived from the PC and the TC bind various miRNAs, some of which have been shown to regulate the physiology and pathology of intestinal development. For example, hsa_circ_0000726 binds miRNA181 and reduces HK and TNF-α expression, inducing glycometabolism and the inflammatory response, respectively [43]. Hsa_circ_083171 binds to miRNA-455, and miR-455-5p downregulates STRA6, reducing cell proliferation [44] and alleviating apoptosis, oxidative stress and the inflammatory response [45]. MiR-455-5p can also alleviate oxygen-induced myocardial cell injury [46]. This study indicated that the exosomal circRNAs derived from HBM might regulate these miRNAs and promote intestinal development. After further analysis, two exosomal circRNAs derived from PC and TC, hsa_circRNA_405708 and hsa_circRNA_104707, respectively, were found to be involved in the VEGF signalling pathway. It was previous reported that circRNAs, such as circ_001621, can sponges the miR-578 and regulate VEGF expression to promote osteosarcoma cell proliferation and migration [47]. These two exosomal circRNAs also bind to many miRNAs, including hsa-miR-4481, hsa-miR-6875-3p, hsa-miR-4745-5p hsa-miR-6827-3p, hsa-miR-4722-5p, hsa-miR-301a-5p, hsa-miR-125a-3p, hsa-miR-134-5p hsa-miR-222-5p and hsa-miR-501-5p. Our results indicated that the PM and TM derived exosomal circRNAs might bind miRNAs, regulate the VEGF signalling pathway, and then influence intestinal development.
Overall, circRNAs have an important effect on multiple pathways in the intestine and should be the focus of future studies. Further research is also needed to validate the potential function of these exosomal circRNAs and their mechanism in intestinal development.
Conclusions
This study confirmed that there were significant differences in PC and TC derived exosomal circRNAs, and that these circRNAs may play important biological roles through their downstream miRNAs. KEGG pathway analysis identified the significant enrichment of the VEGF signalling pathway, and it was shown that PC and TC exosomes can promote the proliferation and migration of FHC cells and VEGF expression. HC exosomes may promote cell proliferation and migration through the VEGF signalling pathway, thereby promoting intestinal development.
Ethics approval and consent to participate
The study was approved by the Institutional Review Board at Women’s Hospital of Nanjing Medical University, Nanjing Maternity and Child Health Care Hospital [Permission Number (2013)78]. Written consent was obtained from all participants.
Data availability statement
The date that support the findings of this study are available from the corresponding author upon reasonable request.
Abbreviations
BM: Breast Milk; PM: Preterm Milk; TM: Term Milk; PC: Preterm Colostrum; TC: Term Colostrum; HBM: Human Breast Milk; HC: Huma Colostrum; TEM: Transmission electron microscopy; NTA: Nanosight; NEC: Necrotizing enterocolitis; FD: Fold change.
Author Contributions
Shuping Han designed the research; Jing Yin, Wenjuan Chen and Yahui Zhou conducted the research and wrote the paper; Yiwen Liu provided the samples, Xingyun Wang and Zhangbin Yu analyzed data; All authors read and approved the final manuscript.
Acknowledgments
We gratefully acknowledge the women who donated breast milk in this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Funding
The grants of National Natural Science Foundation of China (Grant No. 81971427, 81901512, 81870240), Science and Technology Development Fund of Wuxi Health and Family Planning Commission (Q201742) and Science and Technology Development Fund of Nanjing Medical University (No. NMUB2018127) were supported the study.
Editorial Note
This corresponding author has a verified history of publications using the personal email address for correspondence.
References
- 1. Andreas NJ, Kampmann B, Mehring Le-Doare K. Human breast milk: a review on its composition and bioactivity. Early Hum Dev. 2015; 91:629–35. https://doi.org/10.1016/j.earlhumdev.2015.08.013 [PubMed]
- 2. Bauer J, Gerss J. Longitudinal analysis of macronutrients and minerals in human milk produced by mothers of preterm infants. Clin Nutr. 2011; 30:215–20. https://doi.org/10.1016/j.clnu.2010.08.003 [PubMed]
- 3. Lewallen LP. Breastfeeding is important for cognitive development in term and preterm infants. Evid Based Nurs. 2012; 15:85–86. https://doi.org/10.1136/ebnurs-2012-100619 [PubMed]
- 4. Schanler RJ, Shulman RJ, Lau C. Feeding strategies for premature infants: beneficial outcomes of feeding fortified human milk versus preterm formula. Pediatrics. 1999; 103:1150–57. https://doi.org/10.1542/peds.103.6.1150 [PubMed]
- 5. Altobelli E, Angeletti PM, Verrotti A, Petrocelli R. The impact of human milk on necrotizing enterocolitis: a systematic review and meta-analysis. Nutrients. 2020; 12:1322. https://doi.org/10.3390/nu12051322 [PubMed]
- 6. Shoji H, Shimizu T, Shinohara K, Oguchi S, Shiga S, Yamashiro Y. Suppressive effects of breast milk on oxidative DNA damage in very low birthweight infants. Arch Dis Child Fetal Neonatal Ed. 2004; 89:F136–38. https://doi.org/10.1136/adc.2002.018390 [PubMed]
- 7. Penn AH, Altshuler AE, Small JW, Taylor SF, Dobkins KR, Schmid-Schönbein GW. Digested formula but not digested fresh human milk causes death of intestinal cells in vitro: implications for necrotizing enterocolitis. Pediatr Res. 2012; 72:560–67. https://doi.org/10.1038/pr.2012.125 [PubMed]
- 8. Kahn S, Liao Y, Du X, Xu W, Li J, Lönnerdal B. Exosomal MicroRNAs in milk from mothers delivering preterm infants survive in vitro digestion and are taken up by human intestinal cells. Mol Nutr Food Res. 2018; 62:e1701050. https://doi.org/10.1002/mnfr.201701050 [PubMed]
- 9. Martin C, Patel M, Williams S, Arora H, Brawner K, Sims B. Human breast milk-derived exosomes attenuate cell death in intestinal epithelial cells. Innate Immun. 2018; 24:278–84. https://doi.org/10.1177/1753425918785715 [PubMed]
- 10. Hock A, Miyake H, Li B, Lee C, Ermini L, Koike Y, Chen Y, Määttänen P, Zani A, Pierro A. Breast milk-derived exosomes promote intestinal epithelial cell growth. J Pediatr Surg. 2017; 52:755–59. https://doi.org/10.1016/j.jpedsurg.2017.01.032 [PubMed]
- 11. Zeng B, Chen T, Xie MY, Luo JY, He JJ, Xi QY, Sun JJ, Zhang YL. Exploration of long noncoding RNA in bovine milk exosomes and their stability during digestion in vitro. J Dairy Sci. 2019; 102:6726–37. https://doi.org/10.3168/jds.2019-16257 [PubMed]
- 12. Floris G, Zhang L, Follesa P, Sun T. Regulatory role of circular RNAs and neurological disorders. Mol Neurobiol. 2017; 54:5156–65. https://doi.org/10.1007/s12035-016-0055-4 [PubMed]
- 13. Wang Y, Liu J, Ma J, Sun T, Zhou Q, Wang W, Wang G, Wu P, Wang H, Jiang L, Yuan W, Sun Z, Ming L. Exosomal circRNAs: biogenesis, effect and application in human diseases. Mol Cancer. 2019; 18:116. https://doi.org/10.1186/s12943-019-1041-z [PubMed]
- 14. Zhu P, Zhu X, Wu J, He L, Lu T, Wang Y, Liu B, Ye B, Sun L, Fan D, Wang J, Yang L, Qin X, et al. IL-13 secreted by ILC2s promotes the self-renewal of intestinal stem cells through circular RNA circPan3. Nat Immunol. 2019; 20:183–94. https://doi.org/10.1038/s41590-018-0297-6 [PubMed]
- 15. Rankin CR, Lokhandwala ZA, Huang R, Pekow J, Pothoulakis C, Padua D. Linear and circular CDKN2B-AS1 expression is associated with inflammatory bowel disease and participates in intestinal barrier formation. Life Sci. 2019; 231:116571. https://doi.org/10.1016/j.lfs.2019.116571 [PubMed]
- 16. Marell PS, Blohowiak SE, Evans MD, Georgieff MK, Kling PJ, Tran PV. Cord blood-derived exosomal CNTN2 and BDNF: potential molecular markers for brain health of neonates at risk for iron deficiency. Nutrients. 2019; 11:2478. https://doi.org/10.3390/nu11102478 [PubMed]
- 17. Lu D, Xu AD. Mini review: circular RNAs as potential clinical biomarkers for disorders in the central nervous system. Front Genet. 2016; 7:53. https://doi.org/10.3389/fgene.2016.00053 [PubMed]
- 18. Wang X, Yan X, Zhang L, Cai J, Zhou Y, Liu H, Hu Y, Chen W, Xu S, Liu P, Chen T, Zhang J, Cao Y, et al. Identification and peptidomic profiling of exosomes in preterm human milk: insights into necrotizing enterocolitis prevention. Mol Nutr Food Res. 2019; 63:1801247. https://doi.org/10.1002/mnfr.201801247 [PubMed]
- 19. Ballard O, Morrow AL. Human milk composition: nutrients and bioactive factors. Pediatr Clin North Am. 2013; 60:49–74. https://doi.org/10.1016/j.pcl.2012.10.002 [PubMed]
- 20. Karcz K, Królak-Olejnik B. Vegan or vegetarian diet and breast milk composition - a systematic review. Crit Rev Food Sci Nutr. 2021; 61:1081–98. https://doi.org/10.1080/10408398.2020.1753650 [PubMed]
- 21. Jakaitis BM, Denning PW. Human breast milk and the gastrointestinal innate immune system. Clin Perinatol. 2014; 41:423–35. https://doi.org/10.1016/j.clp.2014.02.011 [PubMed]
- 22. Gao HN, Guo HY, Zhang H, Xie XL, Wen PC, Ren FZ. Yak-milk-derived exosomes promote proliferation of intestinal epithelial cells in an hypoxic environment. J Dairy Sci. 2019; 102:985–96. https://doi.org/10.3168/jds.2018-14946 [PubMed]
- 23. Chen T, Xie MY, Sun JJ, Ye RS, Cheng X, Sun RP, Wei LM, Li M, Lin DL, Jiang QY, Xi QY, Zhang YL. Porcine milk-derived exosomes promote proliferation of intestinal epithelial cells. Sci Rep. 2016; 6:33862. https://doi.org/10.1038/srep33862 [PubMed]
- 24. Li T, Sun X, Chen L. Exosome circ_0044516 promotes prostate cancer cell proliferation and metastasis as a potential biomarker. J Cell Biochem. 2020; 121:2118–26. https://doi.org/10.1002/jcb.28239 [PubMed]
- 25. Tang YT, Huang YY, Zheng L, Qin SH, Xu XP, An TX, Xu Y, Wu YS, Hu XM, Ping BH, Wang Q. Comparison of isolation methods of exosomes and exosomal RNA from cell culture medium and serum. Int J Mol Med. 2017; 40:834–44. https://doi.org/10.3892/ijmm.2017.3080 [PubMed]
- 26. Zhang J, Li S, Li L, Li M, Guo C, Yao J, Mi S. Exosome and exosomal microRNA: trafficking, sorting, and function. Genomics Proteomics Bioinformatics. 2015; 13:17–24. https://doi.org/10.1016/j.gpb.2015.02.001 [PubMed]
- 27. Zheng B, Yin WN, Suzuki T, Zhang XH, Zhang Y, Song LL, Jin LS, Zhan H, Zhang H, Li JS, Wen JK. Exosome-mediated miR-155 transfer from smooth muscle cells to endothelial cells induces endothelial injury and promotes atherosclerosis. Mol Ther. 2017; 25:1279–94. https://doi.org/10.1016/j.ymthe.2017.03.031 [PubMed]
- 28. Stras SF, Werner L, Toothaker JM, Olaloye OO, Oldham AL, McCourt CC, Lee YN, Rechavi E, Shouval DS, Konnikova L. Maturation of the human intestinal immune system occurs early in fetal development. Dev Cell. 2019; 51:357–73.e5. https://doi.org/10.1016/j.devcel.2019.09.008 [PubMed]
- 29. Furman D, Davis MM. New approaches to understanding the immune response to vaccination and infection. Vaccine. 2015; 33:5271–81. https://doi.org/10.1016/j.vaccine.2015.06.117 [PubMed]
- 30. Zimmermann P, Curtis N. Factors that influence the immune response to vaccination. Clin Microbiol Rev. 2019; 32:e00084–18. https://doi.org/10.1128/CMR.00084-18 [PubMed]
- 31. Xiu Y, Jiang G, Zhou S, Diao J, Liu H, Su B, Li C. Identification of potential immune-related circRNA-miRNA-mRNA regulatory network in intestine of Paralichthys olivaceus during Edwardsiella tarda infection. Front Genet. 2019; 10:731. https://doi.org/10.3389/fgene.2019.00731 [PubMed]
- 32. Yang L, Fu J, Zhou Y. Circular RNAs and their emerging roles in immune regulation. Front Immunol. 2018; 9:2977. https://doi.org/10.3389/fimmu.2018.02977 [PubMed]
- 33. Peng L, Yin J, Wang S, Ge M, Han Z, Wang Y, Zhang M, Xie L, Li Y. TGF-β2/Smad3 signaling pathway activation through enhancing VEGF and CD34 ameliorates cerebral ischemia/reperfusion injury after isoflurane post-conditioning in rats. Neurochem Res. 2019; 44:2606–18. https://doi.org/10.1007/s11064-019-02880-8 [PubMed]
- 34. Tang H, He Y, Li L, Mao W, Chen X, Ni H, Dong Y, Lyu F. Exosomal MMP2 derived from mature osteoblasts promotes angiogenesis of endothelial cells via VEGF/Erk1/2 signaling pathway. Exp Cell Res. 2019; 383:111541. https://doi.org/10.1016/j.yexcr.2019.111541 [PubMed]
- 35. Sabnis A, Carrasco R, Liu SX, Yan X, Managlia E, Chou PM, Tan XD, De Plaen IG. Intestinal vascular endothelial growth factor is decreased in necrotizing enterocolitis. Neonatology. 2015; 107:191–98. https://doi.org/10.1159/000368879 [PubMed]
- 36. Yan X, Managlia E, Liu SX, Tan XD, Wang X, Marek C, De Plaen IG. Lack of VEGFR2 signaling causes maldevelopment of the intestinal microvasculature and facilitates necrotizing enterocolitis in neonatal mice. Am J Physiol Gastrointest Liver Physiol. 2016; 310:G716–25. https://doi.org/10.1152/ajpgi.00273.2015 [PubMed]
- 37. Wang K, Long B, Liu F, Wang JX, Liu CY, Zhao B, Zhou LY, Sun T, Wang M, Yu T, Gong Y, Liu J, Dong YH, et al. A circular RNA protects the heart from pathological hypertrophy and heart failure by targeting miR-223. Eur Heart J. 2016; 37:2602–11. https://doi.org/10.1093/eurheartj/ehv713 [PubMed]
- 38. Xu L, Zhang M, Zheng X, Yi P, Lan C, Xu M. The circular RNA ciRS-7 (Cdr1as) acts as a risk factor of hepatic microvascular invasion in hepatocellular carcinoma. J Cancer Res Clin Oncol. 2017; 143:17–27. https://doi.org/10.1007/s00432-016-2256-7 [PubMed]
- 39. Yu L, Gong X, Sun L, Zhou Q, Lu B, Zhu L. The circular RNA Cdr1as act as an oncogene in hepatocellular carcinoma through targeting miR-7 expression. PLoS One. 2016; 11:e0158347. https://doi.org/10.1371/journal.pone.0158347 [PubMed]
- 40. Yang X, Xiong Q, Wu Y, Li S, Ge F. Quantitative proteomics reveals the regulatory networks of circular RNA CDR1as in hepatocellular carcinoma cells. J Proteome Res. 2017; 16:3891–902. https://doi.org/10.1021/acs.jproteome.7b00519 [PubMed]
- 41. Hansen TB, Jensen TI, Clausen BH, Bramsen JB, Finsen B, Damgaard CK, Kjems J. Natural RNA circles function as efficient microRNA sponges. Nature. 2013; 495:384–88. https://doi.org/10.1038/nature11993 [PubMed]
- 42. Li F, Zhang L, Li W, Deng J, Zheng J, An M, Lu J, Zhou Y. Circular RNA ITCH has inhibitory effect on ESCC by suppressing the Wnt/β-catenin pathway. Oncotarget. 2015; 6:6001–13. https://doi.org/10.18632/oncotarget.3469 [PubMed]
- 43. Cui Y, Yin K, Gong Y, Qu Y, Liu H, Lin H. Atrazine induces necroptosis by miR-181-5p targeting inflammation and glycometabolism in carp lymphocytes. Fish Shellfish Immunol. 2019; 94:730–38. https://doi.org/10.1016/j.fsi.2019.09.068 [PubMed]
- 44. Zheng J, He Q, Tang H, Xia H. miR-455-5p overexpression reduces rat lung alveolar type II cell proliferation by downregulating STRA6. Anat Rec (Hoboken). 2019; 302:2062–69. https://doi.org/10.1002/ar.24145 [PubMed]
- 45. Chen P, Miao Y, Yan P, Wang XJ, Jiang C, Lei Y. MiR-455-5p ameliorates HG-induced apoptosis, oxidative stress and inflammatory via targeting SOCS3 in retinal pigment epithelial cells. J Cell Physiol. 2019; 234:21915–24. https://doi.org/10.1002/jcp.28755 [PubMed]
- 46. Huang S, Tao W, Guo Z, Cao J, Huang X. Suppression of long noncoding RNA TTTY15 attenuates hypoxia-induced cardiomyocytes injury by targeting miR-455-5p. Gene. 2019; 701:1–8. https://doi.org/10.1016/j.gene.2019.02.098 [PubMed]
- 47. Ji X, Shan L, Shen P, He M. Circular RNA circ_001621 promotes osteosarcoma cells proliferation and migration by sponging miR-578 and regulating VEGF expression. Cell Death Dis. 2020; 11:18. https://doi.org/10.1038/s41419-019-2204-y [PubMed]