- 1Department of Ruminant Production, IRTA, Caldes de Montbui, Spain
- 2Department of Animal and Veterinary Sciences, University of Lleida, Lleida, Spain
- 3ICREA, Institut de Recerca i Estudis Avançats, Barcelona, Spain
Colostrum is rich in membranous vesicles of endocytic origin named exosomes, with proteins, lipids, RNA, and/or DNA cargos which can play different roles in physiological processes. Like other colostrum bioactive compounds, exosomes could be also influenced by individual characteristics. The objective of the study was to characterize miRNA cargo of colostrum exosomes from primiparous and multiparous cows in different farms. Twenty-seven colostrum samples of clinically healthy Holstein cows (11 primiparous and 16 multiparous) from 3 different farms were obtained and frozen. After thawing, exosomes were isolated following an ultracentrifugation protocol, and characterized morphologically. Particle size distribution and western immunoblotting were also analyzMaed. After RNA extraction, miRNAs were sequenced and analyzed to assess potential differences in profiles between primiparous and multiparous cows from different farms. Fourteen miRNA were upregulated and 11 miRNAs downregulated in primiparous compared with multiparous cows. Most of the miRNA differences between primiparous and multiparous cows regulate the gene expression of factors involved in mammary gland development and differentiation, and lipogenesis. In addition, miRNAs from one of the farms showed 8 miRNAs downregulated and 12 upregulated compared with the other 2 farms, independently of parity. Differences in miRNA between farms were mainly associated with immune and inflammatory-related genes. In conclusion, miRNA cargos of bovine colostrum exosomes differ in primiparous and multiparous cows, and some on-farm practices might also determine the content and activity of miRNA in colostrum exosomes.
1 Introduction
Colostrum is an important source of nutrients, and it is also a pivotal source of key hormones and growth factors that may act at both gastrointestinal and systemic levels (1). Colostrum quantity and quality depend on management and environmental factors such as season, calving ease, length of dry period, parity, temperature humidity index, and sex of calves (2), but also on nutritional (3) or metabolic health status of the dams (4). Colostrum quality research has mainly focused on immunoglobulins concentration and, to a lesser extent, on other molecules such as cytokines (5), hormones (3), nutrients (6) or immune cells (7). However, colostrum is also rich in membranous vesicles of endocytic origin, named exosomes, secreted by multiple cell types, found in all biofluids, and considered major players in cell–cell communication (8).
Exosomes can contain proteins, lipids, RNA, and DNA, participating in physiological and pathological processes. miRNAs are a group of small (~22 nucleotides) endogenous RNA molecules that regulate gene expression by binding coding RNA and blocking translation or methylation of promoter CpG-islands. There are more miRNAs in colostrum exosomes than in milk, and some of them are related to immune function and cell development miRNAs (9). It has been suggested that exosomes may participate in the development of the rumen (10), and region-specific miRNA have been found throughout the small intestine, whose expression patterns change throughout the first 6 weeks of the calf life (11). Interestingly, expression changes were most diverse during the first week of life, and many of these miRNAs targeted genes involved in the immune response, suggesting that the miRNA could regulate the development of the immune function during early life. In recent years, several studies have characterized bovine colostrum and milk microvesicles. Milk-isolated microvesicles carrying mRNAs and miRNAs have been related with calf gastrointestinal development and immune system (12). Furthermore, there exist differences in bovine colostrum miRNA between cow breeds (13), or between cows with different colostrum IgG concentration (14). We hypothesized that cow parity could influence colostrum miRNA in exosomes. Therefore, the objective of this study was to isolate colostrum exosomes and compare their miRNA cargos from primiparous and multiparous cows of different dairy farms.
2 Materials and methods
2.1 Sample collection
One hundred mL of unpasteurized bovine colostrum (first milking) was obtained after calving from 27 clinically healthy Holstein cows [11 primiparous and 16 multiparous (7 of second lactation, 5 of third lactation and 4 of forth lactation)] in 3 different farms all located in the north-east of Spain from July to November 2020. Samples were collected in sterile 50-mL conical tubes and immediately frozen at −20°C. All farms offered a single diet throughout the dry period, but its ingredient and nutrient composition differed by farm as described in Supplementary Table S1. All farms targeted 60-days for the length of the dry period, but actual duration differed among individuals within farms.
2.2 Exosomes isolation
Colostrum exosomes were isolated using ultracentrifugation based on (15) protocol with some modifications. After thawing the samples at room temperature, they were centrifuged twice at 12,500 × g for 15 min at 4°C and the upper fat layer removed. Then, samples were diluted 1:2 with PBS and centrifuged at 50,000 × g for 1 h at 4°C in a high-speed centrifuge. Then, four fifths of the supernatant were discarded (23 mL), and the remaining supernatant was washed in PBS for a final centrifugation at 50,000 × g for 2 h at 4°C in 38-mL ultracentrifuge tubes in a high-speed centrifuge. Then, the supernatant was discarded, and the smooth pellet formed on the top of the firm casein pellet was recovered and resuspended in PBS for exosome characterization analysis.
2.3 Exosomes characterization
2.3.1 Cryo-transmission electron microscopy (Cryo-TEM)
Ultrastructural morphology of purified exosomes was analyzed by Cryo-Transmission Electron Microscopy (Cryo-TEM). A 3.9 μL drop of the purified exosomes samples were deposited onto holey carbon on a 400-mesh copper grid, previously treated by glow discharge. The grid was mounted on a plunger (Leica EM GP main unit, Leica Microsystems, Vienna, Austria), water excess was removed by blotting with filter paper (Whatman n°1 Filter Paper 90 mm Diameter, Silmid, Birmingham, UK) and exosomes suspension was straightforward vitrified by rapid immersion in liquid ethane (−178°C). Samples were mounted on a Gatan 626 cryo-transfer system and inserted into a Jeol JEM 2011 under cryo conditions. Electron microscope operating energy was 200 kV. Images were recorded using a Gatan Ultrascan US1000 CCD camera.
2.3.2 Nanoparticle tracking analysis (NTA)
Particle size distribution and concentration were further assessed by Nanoparticle Tracking Analysis (NTA) with Nanosight NS300 (Malvern Panalytical, Malvern, UK). A 1:2,000–1:5,000 dilution of the sample in PBS was performed. Data were further processed using NanoSight Software NTA 3.4.
2.3.3 Western immunoblotting
Purified exosomes samples were analyzed by western blot following the protocol previously described (16). Eight μL of purified exosomes were loaded per lane. Primary antibodies used were against CD9 (1:500, Invitrogen, Massachusetts), MFGE8 (1:1,000, Sigma, MO, USA), TSG101 (1:250, Invitrogen, MA, USA), and Alix (1:1,000, Cell Signaling Technology, MA, USA). Secondary anti-mouse Alkaline Phosphatase-conjugated antibody (1:20,000, Sigma, MO, USA) was used for CD9 and TSG101 and anti-rabbit Alkaline Phosphatase-conjugated antibody (1:30,000, Sigma, MO, USA) was used for MFGE8 and Alix. Positive control corresponding to bovine milk exosomes (Lyophilized Exosome Standards (Bovine Milk) ref: ESL-01, Creative Biolabs) was used.
2.3.4 RNA sequencing and profiling
RNA from 200 μL of purified exosomes was extracted through the miRNeasy Serum/Plasma kit (Qiagen). RNA quality and quantity were assessed using an Agilent 2100 Bioanalyzer and RNA ladder. Small RNA libraries were prepared using the Small RNA_NebNext kit (NEB, Ipswich, MA, USA) according to the manufacturer’s instructions. Briefly, 1 μg of total RNA was subjected to DNase I treatment (RNase-Free) to remove genomic DNA contamination. RNA fragments of 10–40 nt in length were generated using the fragmentation buffer. The RNA fragments were then ligated with 3′ and 5′ adapters, followed by reverse transcription and PCR amplification to generate the sequencing library. The final libraries were analyzed using the Agilent 2100 Bioanalyzer. The libraries were sequenced on an Illumina HiSeq2500 with a single-end 50 bp read length according to the manufacturer’s instructions. The sequencing depth for each sample was 10 million reads.
2.4 Data analysis
The quality of the miRNA profiling raw data was checked using FastQC (17). The reads were trimmed for QC and the presence of the Universal Illumina adapter using Skewer v.2.0.0 (18). Trimmed reads were mapped to the Bos taurus reference genome (Ensembl, release 104, file Bos_taurus.ARS-UCD1.2.dna.fa.gz) using Bowtie v.1.2.2 (19) and counts were assigned to annotated miRNA genes (Ensembl, release 104, file Bos_taurus.ARS-UCD1.2.104.chr.gtf.gz) using HTSeq v.2.0 (20). The number of unaligned tags was increased to the refinement made using ShortStack v.3 (21) by removing the multi-mapping reads impossible to reassign. Qualimap v.2.2.1 (22) was used to check the quality of the aligned reads. Differential expression analysis was performed in R (version 4.0.0) (23)/Bioconductor (24) environment, using the DESeq2 package version 1.28.1 (25). Genes that had less than 10 read counts across all samples were filtered out before processing the data. Principal component analysis (PCA) was conducted using the prcomp method from the stats core R package on the regularized log-transformed data and plotted using the ggplot2 package (26). First parity or other parities were the fixed effect, and farms were considered independent among them. Samples 604, 6,598, and 9,280 were removed from the analysis as outliers. The expression of miRNAs genes between type of animals (primiparous vs. multiparous) within farms were selected as differentially expressed if the FDR-adjusted p-value was less than 0.05, and the absolute log2 fold change was more than 1. In addition, bedtools v2.27.1 (27) and samtools 1.8 (28) were used. Similarly, the different expressions of miRNAs among the three farms, compared by pairs, were evaluated following the same procedure.
3 Results
3.1 Exosome characterization
Bovine exosomes were characterized by Cryo-TEM, NTA, and Western immunoblotting. Images of Cryo-TEM showed that colostrum vesicles had a round appearance (Figure 1). No differences were observed in particle size or yield either in primiparous or multiparous cows, or among the different farms (Figure 2). The mean particle size of all samples was 182 ± 8.3 nm with a main peak at 125–135 nm, and they were within the 10 and 90 percentile range of 111 ± 2.4 nm and 276 ± 10.6 nm, respectively, and mean particle yield was of 1.97 ± 0.230×1012 particles/mL. Western blot confirmed that samples were enriched with exosome-specific markers MFGE8 and ALIX (Figure 3). Marker CD9 was not detected in colostrum samples whereas TSG101 was detected in all analyzed samples but not in the positive control used (Supplementary Figure S1).
Figure 2. Examples of nanoparticles size distribution of isolated exomes samples from primiparous cows in Farm I (a), Farm 2 (b), and Farm 3 (c).
Figure 3. Western blot analyses of the presence of two exosome markers (Alix and MFGE8) in samples of either primiparous or multiparous cows of the 3 farms included in the study. Positive control corresponds to commercial purified exosomes from bovine milk (Lyophilized Exosome Standards, Creative Biolabs).
3.2 Characterization of miRNAs
A total of 943 miRNAs were isolated, with 88.4% of them associated with known genes. Twenty miRNAs represented between 80 and 85% of total miRNAs reads, from those, 12 were known (bta-miR-30a, bta-miR-26a-2, bta-miR-99a, bta-let-7f-2, bta-miR-200a, bta-let-7a-3, bta-miR-2285 t, bta-miR-200c, bta-miR-27b, bta-miR-30d, bta-miR-30f, bta-let-7a-2). The comparison between primiparous and multiparous cows revealed 14 miRNAs downregulated (bta-miR-2478, bta-miR-885, bta-miR-503, bta-miR-708, bta-miR-424, bta-miR-9-1, bta-miR-218-1, bta-miR-452, bta-miR-497, bta-miR-195, bta-miR-3431, bta-miR-378c, bta-miR-224, and bta-miR-181c, in order of log2Folds magnitude), and 11 miRNAs upregulated (bta-miR-128-2, bta-miR-143, bta-miR-145, bta-miR-181a-1, bta-miR-181b-2, bta-miR-1388, bta-miR-181a-2, bta-miR-193a, bta-miR-2299, bta-miR-28, and bta-miR-10174, in the order of 2 log-folds) in primiparous in comparison with multiparous cows. The overview of these differences was represented in a heatmap as Supplementary Figure S2.
When assessing potential anecdotical differences in colostrum samples among farms, one farm was distinguished in the PCA plot from the other 2 farms (Figure 4), and it was also revealed in the DESeq2 analysis. Colostrum samples from the outstanding farm had 8 downregulated (bta-miR-99a, bta-miR-181d, bta-miR-126, bta-miR-195, bta-miR-2387, bta-miR-2285e-1, and 2 unknown miRNAs in chromosome 7and 8), and 12 upregulated (bta-miR-221, bta-miR-146a, bta-miR-345, bta-miR-222, bta-miR-142, bta-miR-155, bta-miR-223, bta-miR-2284w, bta-miR-378c, bta-miR-378-2 and 2 unknown miRNAs in chromosome 7 and 19) miRNAs compared with the other 2 farms, which did not differ in the expression of any miRNAs.
Figure 4. Principal Component Analysis plot of miRNA of exomes isolated from colostrum bovine samples from primiparous (parity l) and multiparous (parity 2) cows from three different dairy farms. Parity 2 includes parities from 2 to 4 lactations (animal ID of 2 lactations: 510, 6,129, 6,131, 6,598, 7,404, 10,000, 10,640, animal ID of 3 lactations: 462, 464, 471, 5,683, 7,215, animal ID of4 lactations: 402, 3,055, 5,036, 17,224).
4 Discussion
Herein mean particle size of exosomes (182 ± 8.3 nm) was greater than that described in the literature in colostrum samples (149 ± 8.6 nm) (29) or in milk samples (100 nm) (30). The mean particle yield was of 1.97 ± 0.230×1012 particles/mL, which was also greater than that reported in the literature (1.4×1011 particles/mL) (31), probably mainly due to different exosome isolation protocols. In the current study, several filtration procedures (filtration using cellulose acetate filters (1′2, 0′8 and 0′45 mm) or PVDF filters (0′45 mm), dialysis membranes with a 300 kDa cut-off, and centrifuge filters with a 300 kDa cut-off) were tested without affecting the quality of exosomes (Supplementary Figure S3). By contrast, previous studies reported a benefit of using filter membranes before (29, 31) or after (30) ultracentrifugation, or using a size exclusion chromatography column after centrifugation (31). The lack of filtration steps in the current samples may have increased particle size and number of exosomes in comparison with previous studies. Similar to our study, TSG101 and Alix were found in colostrum samples (8), and MFGE8 was detected in bovine raw milk (32). Colostrum samples herein did not mark for CD9, which it has been previously reported to be present in exosomes from bovine milk (29, 31), and bovine transition milk samples collected during the first and second days in milk of lactating cows (33).
The highly expressed miRNAs found were common among all animals and farms, and they have also been reported as common miRNA in milk or colostrum samples in previous studies (34, 35). Some of the different miRNA contents in primiparous and multiparous cows regulate the expression of genes participating in mammary gland development and differentiation. The bta-miR-424/503 cluster and bta-miR-452 were downregulated in primiparous cows, and they have been described to be up-regulated during mammary gland involution (36, 37). Specifically, bta-miR-424/503 cluster has been reported to be involved in the transcriptional regulation of TGF-ß (36). Furthermore, bta-miR-708 was also downregulated in primiparous cows, and it has been reported to regulate the expression of TNFSF11, EGF, and HOXA5, which are genes involved in the development and differentiation of the mammary gland (38). Lastly, bta-miR-128-2 was upregulated in primiparous cows, and it has been described to be repressed by TGF-ß during mammary epithelial oncogenic transformation (39). Some other miRNAs that differed between primiparous and multiparous cows have been involved in lipogenesis. On one hand, some were downregulated in primiparous cows: bta-miR-497, which has been involved in the regulation the fatty acid synthesis (40), bta-miR-378c, which is related to maternal body conditional index (41), and bta-miR-3431 which has been reported to be upregulated in dairy cows when they were supplemented with linseed oil (42). On the other, the miRNA upregulated in primiparous cows were bta-miR-193a, which has been involved in adipocyte differentiation (43), bta-miR-181a-1, bta-miR-181b-2, bta-miR-143, and bta-miR-145, which have been involved in milk fat synthesis (44–47), and bta-miR-2299 and bta-miR-1388, which they were downregulated and bta-miR-28, which was upregulated, in dairy cows when they were supplemented with linseed or safflower oils (42). Other miRNAs downregulated in primiparous cows were bta-miR-2478, which has been associated with the inhibition of melanin production (48), bta-miR-885 which has been related with retinol biosynthesis (37), and two miRNAs (bta-miR-195, and bta-miR-181c) linked to tumor growth signaling pathways (49, 50).
Although it was not the objective of the study and animals were only checked by clinical diseases signs and clinical mastitis, we observed that miRNA profile of colostrum exosomes of one farm differed from the other 2 in the PCA. The miRNA found to differ upregulate genes related to immune and inflammatory responses like in Staphylococcus aureus infections [bta-miR-378-2 and bta-miR-223; (51, 52), respectively], in the regulation of lymphocyte functions [bta-miR-155; (53)], or as biomarkers in inflammatory diseases [bta-miR-221 and bta-miR-222; (54)]. Those miRNAs that were downregulated in that specific farm were more diverse in their regulation responses: bta-miR-126 was found to regulate vascular integrity (55), bta-miR-181d has been involved in the regulation of inflammatory response in heat-stressed cows (56), or bta-miR-2387 targeted lipid metabolism-related genes in (57) study. Although the study of differences among farms was exploratory, and it was not targeting any specific farm management practice or animal diseases, the fact that 2 farms had similar miRNA profiles (while all farms differed in nutrients and ingredients composition, Supplementary Table S1), may exclude diet composition as an influencing factor on colostrum exosomes miRNAs in the present study. Since most of the upregulated mi-RNAs in the outstanding farm were immune and inflammatory related, and this farm was not vaccinating cows during the dry period while the other 2 farms did, vaccination might be an important factor behind this anectodical difference in miRNA composition of colostrum.
Exosomes in colostrum from primiparous and multiparous cows differed in 25 miRNAs. Most of these miRNAs are involved in the regulation of lipogenic pathways, and development and differentiation of the mammary gland.
Data availability statement
Sequencing data are available in the NCBI BioProject under accession number PRJNA1165141.
Ethics statement
The requirement of ethical approval was waived by Comitè d’ètica i d’experimentació animal de l’IRTA for the studies involving animals because colostrum samples were obtained under the current farm practices. There were no animal invasive practices to obtain the samples. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was not obtained from the owners for the participation of their animals in this study because we considered an oral consent as valid.
Author contributions
MT: Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing. AA: Methodology, Validation, Writing – review & editing. EG-F: Methodology, Supervision, Writing – review & editing. FF: Investigation, Methodology, Resources, Writing – review & editing. AB: Conceptualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the Spanish Ministry of Economy and Competitivity and FEDER (PGC2018-098847-B-100). IRTA thanks the funding support of the CERCA Program and the Consolidated Research Groups of Sustainable Animal Husbandry (ref. 2021 SGR 01552), both from the Generalitat de Catalunya.
Acknowledgments
The authors would like to thank all the farm personnel involved in collecting colostrum samples. The characterization of exosomes with the Nanosight has been performed by Jose Amable Bernabe in the ICTS “NANBIOSIS,” more specifically by the Unit 6. Unit of the CIBER in Bioengineering, Biomaterials & Nanomedicne (CIBER-BBN) at the Barcelona Materials Science Institute. The CryoTEM analysis was done at “Servei de Microscopia” (UAB). The miRNA sequencing and analysis was performed by Berta Fusté and Julia Ponomarenko at Genomics and Bioinformatics units at Center for Genomic Regulation (CRG).
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.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
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/fvets.2024.1463342/full#supplementary-material
References
1. McGrath, BA, and Fox, PF. Composition and properties of bovine colostrum: a review. Dairy Sci Technol. (2016) 96:133–58. doi: 10.1007/s13594-015-0258-x
2. Borchardt, S, Sutter, F, Heuwieser, W, and Venjakob, P. Management-related factors in dry cows and their associations with colostrum quantity and quality on a large commercial dairy farm. J Dairy Sci. (2022) 105:1589–602. doi: 10.3168/jds.2021-20671
3. Mann, S, Leal Yepes, FA, Overton, TR, Lock, AL, Lamb, SV, Wakshlag, JJ, et al. Effect of dry period dietary energy level in dairy cattle on volume, concentration of immunoglobulin G, insulin, and fatty acid composition of colostrum. J Dairy Sci. (2016) 99:1515–26. doi: 10.3168/jds.2015-9926
4. Immler, M, Failing, K, Gärtner, WA, and Donaat, K. Association between the metabolic status of the cow and colostrum quality as determined by brix refractometry. J Dairy Sci. (2021) 104:10131–42. doi: 10.3168/jds.2020-19812
5. Hagiwara, K, Kataoka, S, Yamanaka, H, Kirisawa, R, and Iwai, H. Detection of cytokines in bovine colostrum. Vet Immunol Immunopathol. (2000) 76:183–90. doi: 10.1016/s0165-2427(00)00213-0
6. Valldecabres, A, and Silva-del-Río, N. First-milking colostrum mineral concentrations and yields: comparison to second milking and associations with serum mineral concentrations, parity, and yield in multiparous Jersey cows. J Dairy Sci. (2022) 105:2315–25. doi: 10.3168/jds.2021-21069
7. Langel, SN, Wark, WA, Garst, SN, James, RE, McGilliard, MLL, Petersson-Wollfe, CD, et al. Effect of feeding whole compared with cell-free colostrum on calf immune status: the neonatal period. J Dairy Sci. (2015) 98:3729–40. doi: 10.3168/jds.2014-8422
8. Samuel, M, Chisanga, D, Liem, M, Keerthikumar, S, Anand, S, Ang, CS, et al. Bovine milk-derived exosomes from colostrum are enriched with proteins implicated in immune response and growth. Sci Rep. (2017) 7:5933. doi: 10.1038/s41598-017-06288-8
9. Izumi, H, Kosaka, N, Shimizu, T, Sekine, K, Ochiva, T, and Takase, M. Bovine milk contains microRNA and messenger RNA that are stable under degradative conditions. J Dairy Sci. (2012) 95:4831–41. doi: 10.3168/jds.2012-5489
10. He, L, and Hannon, GJ. MicroRNAs:small RNAs with a big role in gene regulation. Nat Rev Genet. (2004) 5:522–31. doi: 10.1038/nrg1379
11. Liang, G, Malmuthuge, N, McFadden, TB, Bao, H, Griebel, PJ, and Stothard, GLL. Potential regulatory role of MicroRNAs in the development of bovine gastrointestinal tract during early life. PLoS One. (2014) 9:e92592. doi: 10.1371/journal.pone.0092592
12. Hata, T, Murakami, K, Nakatani, H, Yamaamoto, Y, Matsuda, T, and Aoki, N. Isolation of bovine milk-derived microvesicles carrying mRNAs and microRNAs. Biochem Biophys Res Commun. (2010) 396:528–33. doi: 10.1016/j.bbrc.2010.04.135
13. Özdemir, S. Identification and comparison of exosomal microRNAs in the milk and colostrum of two different cow breeds. Gene. (2020) 743:144609. doi: 10.1016/j.gene.2020.144609
14. Ma, T, Li, W, Chen, Y, Cobo, ER, Windeyer, C, Gamsjäger, L, et al. Assessment of microRNA profiles in ssmall extracellular vesicles isolated from bovine colostrum with different immunoglobulin G concentrations. JDS Commun. (2022) 3:328–33. doi: 10.3168/jdsc.2022-0225
15. Ferreira, RF, Blees, T, Shakeri, F, Buness, A, Sylvesterr, M, Savoini, G, et al. Comparative proteome profiling in exosomes derived from porcine colostrum versus mature milk reveals distinct functional proteomes. J Proteome. (2021) 249:104338. doi: 10.1016/j.jprot.2021.104338
16. López-Cano, A, Martíinez-Miquel, M, Guasch, J, Ratera, I, Arís, A, and Garcia-Fruitós, E. Exploring the impact of the recombinant Escherichia coli strain on defensins antimicrobial activity. BL21 versus origami strain. Microb Cell Factories. (2022) 21:77. doi: 10.1186/s12934-022-01803-7
17. Andrews, S. FastQC: a quality control tool for high throughput sequence data. (2010). Available at: http://www.bioinformatics.babraham.ac.uk/projects/fastqc
18. Jiang, H, Lei, R, Ding, SW, and Zhu, S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. (2014) 15:182. doi: 10.1186/1471-2105-15-182
19. Langmead, B, Trapnell, C, Pop, M, and Salzberg, SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. (2009) 10:R25. doi: 10.1186/gb-2009-10-3-r25
20. Putri, G, Anders, S, Pyl, PT, Pimanda, JE, and Zanini, F. Analysing high-throughput sequencing data in Python with HTSeq 2.0. Bioinformatics. (2022) 38:2943–5. doi: 10.1093/bioinformatics/btac166
21. Axtell, MJ. ShortStack: comprehensive annotation and quantification of small RNA genes. RNA. (2013) 19:740–51. doi: 10.1261/rna.035279.112
22. Okonechnikov, K, Conesa, A, and García-Alcalde, F. Qualimap 2: advanced multi-sample quality control for high-throughput sequencing data. Bioinformatics. (2016) 32:292–4. doi: 10.1093/bioinformatics/btv566
23. R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing (2020).
24. Huber, W, Carey, VJ, Gentleman, R, Anders, S, Carlson, M, Carvalho, BS, et al. Orchestrating high-throughput genomic analysis with Bioconductor. Nat Methods. (2015) 12:115–21. doi: 10.1038/nmeth.3252
25. Love, MI, Huber, W, and Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. (2014) 15:550. doi: 10.1186/s13059-014-0550-8
27. Quinlan, AR. BEDTools: the Swiss-Army tool for genome feature analysis. Curr Protoc Bioinformatics 47:11.12.1-34. (2014) 47:11.12.1–34. doi: 10.1002/0471250953.bi1112s47
28. Li, H, Handsaker, B, Wysoker, A, Fennell, T, Ruan, J, Homer, N, et al. 1000 genome project data processing subgroup. The sequence alignment/map format and SAM tools. Bioinformatics. (2009) 25:2078–9. doi: 10.1093/bioinformatics/btp352
29. Wijenayake, S, Eisha, S, Tawhidi, Z, Pitino, MA, Steele, MA, Fleming, AS, et al. Comparison of methods for pre-processing, exosome isolation, and RNA extraction in unpasteurized bovine and human milk. PLoS One. (2021) 16:e0257633. doi: 10.1371/journal.pone.0257633
30. Cai, M, He, H, Jia, X, Chen, S, Wang, J, Shi, Y, et al. Genome-wide microRNA profiling of bovine milk-derived exosomes infected with Staphylococcus aureus. Cell Stress Chaperones. (2018) 23:663–72. doi: 10.1007/s12192-018-0876-3
31. Vaswani, K, Mitchell, MD, Holland, OJ, Koh, YQ, Hill, RJ, Harb, T, et al. Method for the isolation of exosomes from human and bovine milk. J Nutr Metab. (2019) 2019:5764740. doi: 10.1155/2019/5764740
32. Yamauchi, M, Shimizu, K, Rahman, M, Ishikawa, H, Takase, H, Ugawa, S, et al. Efficient method for isolation of exosomes from raw bovine milk. Drug Dev Ind Pharm. (2018) 45:359–64. doi: 10.1080/03639045.2018.1539743
33. Maburutse, BE, Park, MR, Oh, S, and Kim, Y. Evaluation and characterization of Milkderived Microvescicle isolated from bovine colostrum. Korean J Food Sci Ann. (2017) 37:654–62. doi: 10.5851/kosfa.2017.37.5.654
34. Chen, X, Gao, C, Li, H, Huang, L, Sun, Q, Dong, Y, et al. Identification and characterization of microRNAs in raw milk during different periods of lactation, commercial fluid, and powdered milk products. Cell Res. (2010) 20:1128–37. doi: 10.1038/cr.2010.80
35. Yun, B, Kim, Y, Park, DJ, and Oh, S. Comparative analysis of dietary exosome-derived microRNAs from human, bovine and caprine colostrum and mature milk. J Anim Sci Technol. (2021) 63:593–602. doi: 10.5187/jast.2021.e39
36. Llobet-Navas, D, Rodríguez-Barrueci, R, Castro, V, Ugalde, AP, Sumazin, P, Jacob-Sendler, D, et al. The miR-424(322)/503 cluster orchestrates remodeling of the epithelium in the involuting mammary gland. Genes Dev. (2014) 28:765–82. doi: 10.1101/gad.237404.114
37. Do, DN, Li, R, Dudemaine, PL, and Ibeagha-Awemu, EM. MicroRNA roles in signalling during lactation: an insight from differential expression, time course and pathway analyses of deep sequence data. Sci Rep. (2017) 7:44605. doi: 10.1038/srep44605
38. Xuan, R, Chao, T, Wang, A, Zhang, F, Sun, P, Liu, S, et al. Characterization of microRNA profiles in the mammary gland tissue of dairy goats at the late lactation, dry period and late gestation stages. PLoS One. (2020) 15:e0234427. doi: 10.1371/journal.pone.0234427
39. Qian, PX, Banerjee, A, Wu, ZS, Zhang, X, Wang, H, Pandey, V, et al. Loss of SNAIL regulated miR-128-2 on chromosome 3p22.3 targets multiple stem cell factors to promote transformation of mammary epithelial cells. Cancer Res. (2012) 72:6036–50. doi: 10.1158/0008-5472.CAN-12-1507
40. Chen, Z, Xie, Y, Luo, J, Chen, T, Xi, Q, Zhang, Y, et al. Milk exosome-derived miRNAs from water buffalo are implicated in immune response and metabolism process. BMC Vet Res. (2020) 16:123. doi: 10.1186/s12917-020-02339-x
41. Xi, Y, Jiang, X, Li, R, Chen, M, Song, W, and Li, X. The levels of human milk microRNAs and their association with maternal weight characteristics. Eur J Clin Nutr. (2016) 70:445–9. doi: 10.1038/ejcn.2015.168
42. Li, R, Beaudoin, F, Ammah, AA, Bissonnette, N, Benchaar, C, Zhao, X, et al. Deep sequencing shows microRNA involvement in bovine mammary gland adaptation to diets supplemented with linseed oil or safflower oil. BMC Genomics. (2015) 16:884. doi: 10.1186/s12864-015-1965-7
43. Han, F, Zhou, L, Zhao, L, Wang, L, Liu, L, Li, H, et al. Identification of miRNA in sheep intramuscular fat and the role of miR-193a-5p in proliferation and differentiation of 3T3-L1. Front Genet. (2021) 12:633295. doi: 10.3389/fgene.2021.633295
44. Chen, Z, Shi, HP, Sun, S, Xu, HF, Cao, DY, and Luo, J. MicroRNA-181b suppresses TAG via target IRS2 and regulating multiple genes in the hippo pathway. Exp Cell Res. (2016) 348:66–74. doi: 10.1016/j.yexcr.2016.09.004
45. Lian, S, Guo, JR, Nan, XM, Ma, L, Loor, JJ, and Bu, DP. MicroRNA Bta-miR-181a regulates the biosynthesis of bovine milk fat by targeting ACSL1. J Dairy Sci. (2016) 99:3916–24. doi: 10.3168/jds.2015-10484
46. Wang, H, Shi, H, Luo, J, Yi, Y, Yao, D, Zhang, X, et al. MiR-145 regulates lipogenesis in goat mammary cells via targeting INSIG1 and epigenetic regulation of lipid-related genes. J Cell Physiol. (2017) 232:1030–40. doi: 10.1002/jcp.25499
47. Zhang, L, Wu, ZQ, Wang, YJ, Wang, M, and Yang, WC. MiR-143 regulates Milk fat synthesis by targeting Smad3 in bovine mammary epithelial cells. Animals. (2020) 10:1453. doi: 10.3390/ani10091453
48. Bae, IS, and Kim, SH. Milk exosome-derived MicroRNA-2478 suppresses Melanogenesis through the Akt-GSK3 pathway. Cells. (2021) 10:2848. doi: 10.3390/cells10112848
49. Fu, Y, Tang, Y, Wang, J, and Guo, Z. MicroRNA-181c suppresses the biological progression of osteosarcoma via targeting SMAD7 and regulating transforming growth factor-β (TGF-β) signaling pathway. Med Sci Monit. (2019) 25:4801–10. doi: 10.12659/MSM.916939
50. Rao, G, Dwivedi, SKD, Zhang, Y, Dey, A, Shameer, K, Karthik, R, et al. MicroRNA-195 controls MICU1 expression and tumor growth in ovarian cancer. EMBO Rep. (2020) 21:e48483. doi: 10.15252/embr.201948483
51. Ma, S, Tong, C, Ibeagha-Awemu, EM, and Zhao, X. Identification and characterization of differentially expressed exosomal microRNAs in bovine milk infected with Staphylococcus aureus. BMC Genomics. (2019) 20:934. doi: 10.1186/s12864-019-6338-1
52. Han, S, Li, X, Liu, J, Zou, Z, Luo, L, Wu, R, et al. Bta-miR-223 targeting CBLB contributes to resistance to Staphylococcus aureus mastitis through the PI3K/AKT/NF-κB pathway. Front Vet Sci. (2020) 7:529. doi: 10.3389/fvets.2020.00529
53. Lawless, N, Vegh, P, O’Farrelly, C, and Lynn, DJ. The role of microRNAs in bovine infection and immunity. Front Immunol. (2014) 5:611. doi: 10.3389/fimmu.2014.00611
54. Song, J, Ouyang, Y, Che, J, Li, X, Zhao, Y, Yang, K, et al. Potential value of miR-221/222 as diagnostic, prognostic, and therapeutic biomarkers for diseases. Front Immunol. (2017) 8:56. doi: 10.3389/fimmu.2017.00056
55. Fish, JE, Santoro, MM, Morton, SU, Yu, S, Yeh, RF, Wythe, JD, et al. miR-126 regulates angiogenic signalling and vascular integrity. Dev Cell. (2008) 15:272–84. doi: 10.1016/j.devcel.2008.07.008
56. Fan, C, Hu, RR, Fan, H, Yang, Y, Gong, B, Zhang, S, et al. Effects of seasonal ambient heat stress on expression of microRNAs in the mammary gland of Holstein cows. Int J Biometeorol. (2021) 65:235–46. doi: 10.1007/s00484-020-02025-5
Keywords: colostrum, dairy cow, exosome, miRNA, parity
Citation: Terré M, Arís A, Garcia-Fruitós E, Fàbregas F and Bach A (2024) Micro RNA profiles in colostrum exosomes obtained from primiparous or multiparous dairy cows. Front. Vet. Sci. 11:1463342. doi: 10.3389/fvets.2024.1463342
Edited by:
Kai Wang, Chinese Academy of Agricultural Sciences (CAAS), ChinaReviewed by:
Diana Giannuzzi, University of Padua, ItalyLiang Deng, Shenyang Agricultural University, China
Copyright © 2024 Terré, Arís, Garcia-Fruitós, Fàbregas and Bach. 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: Marta Terré, marta.terre@irta.cat