Comparison of goat and cow milk-derived extracellular vesicle miRNomes.

miRNAs present in milk are mainly found in extracellular vesicles (EVs), which are nanosized membrane vesicles released by most of the cell types to ensure intercellular communication. The majority of the studies performed so far on these vesicles have been conducted on human and cow's milk and focused on their miRNA content. The objectives of this study were to profile the miRNA content of purified EVs from five healthy goats and to compare their miRNome to those obtained from five healthy cows, at an early stage of lactation. EV populations were morphologically characterized using Transmission Electron Microscopy and Nanoparticle Tracking Analysis. The presence of EV protein markers checked by Western blotting and the absence of contamination of preparations by milk proteins. The size distribution and concentration of bovine and goat milk-derived EVs were similar. RNA-sequencing were performed, and all sequences were mapped to the cow genome identifying a total of 295 miRNAs. This study reports for the first-time a goat miRNome from milk EVs and its validation using cow miRNomes.

Top-Down proteomics based on LC-MS combined with cDNA sequencing to characterize multiple proteoforms of Amiata donkey milk proteins.

An in-depth molecular characterization of the main milk proteins, caseins (CNs) and whey proteins, from Amiata donkey combining top-down proteomic analysis (LC-MS) and cDNA sequencing revealed multiple proteoforms arising from complex splicing patterns, including cryptic splice site usage and exon skipping events. Post-translational modifications, in particular phosphorylation, increased the variety and complexity of proteoforms. αs2-CN perfectly exemplifies such a complexity. With 2 functional genes, CSN1S2 I and CSN1S2 II, made of 20 and 16 exons respectively, nearly 30 different molecules of this CN were detected in the milk of one Amiata donkey. A cryptic splice site usage, leading to a singular shift of the open reading frame and generating two αs2-CN I isoforms with different C-terminal sequences, was brought to light. Twenty different αs1-CN molecules with different phosphorylation levels ranging between 4 and 9P were identified in a single milk sample, most of them resulting from exon skipping events and cryptic splice site usage. Novel genetic polymorphisms were detected for CNs (ß- and αs-CN) as well as for whey proteins (lysozyme C and ß-LG I). The probable new ß-LG I variant, with a significantly higher mass than known variants, appears to display an N-terminal extension possibly related to the signal peptide sequence. This represents the most comprehensive report to date detailing the complexity of donkey milk protein micro-heterogeneity, a prerequisite for discovering new elements to objectify the original properties of donkey's milk.

Comparative transcriptome analysis of goat (Capra hircus) adipose tissue reveals physiological regulation of body reserve recovery after the peak of lactation.

Adipose tissue is the energy storage organ providing energy to other tissues, including mammary gland, that supports the achievement of successive lactation cycles. Our objective was to investigate the ability of goats to restore body fat reserves by comparing lipogenic enzyme activities and by transcriptomic RNA-Seq data at two different physiological stages, mid- and post-lactation. Key lipogenic enzyme activities were higher in goat omental adipose tissue during mid-lactation (74 days in milk) than during the post-lactation period (300 days postpartum). RNA-Sequencing analysis revealed 19,271 expressed genes in the omental adipose tissue. The comparison between adipose transcriptome analysis from mid- and post-lactation goats highlighted 252 differentially expressed genes (padj < 0.05) between these two physiological stages. The differential expression of 11 genes was confirmed by RT-qPCR. Functional genomic analysis revealed that 31% were involved in metabolic processes among which 38% in lipid metabolism. Most of the genes involved in lipid synthesis and those in lipid transport and storage were upregulated in adipose tissue of mid- compared to post-lactation goats. In addition, adipose tissue plasticity was emphasized by genes involved in cellular signaling and tissue integrity. Network analyses also highlighted three key regulators of lipid metabolism (LEP, APOE and HNF4A) and a key target gene (VCAM1). The greatest lipogenic enzyme activities with the upregulation of genes involved in lipid metabolism highlighted a higher recovery of lipid reserves after the lactation peak than 4 months post-lactation. This study contributes to a better understanding of the molecular mechanisms controlling the body lipid reserves management during the successive lactations.

Different miRNA contents between mammary epithelial cells and milk fat globules: a random or a targeted process?

MicroRNAs (miRNAs) are small noncoding RNAs present in milk-derived extracellular vesicles and milk fat globules (MFG). Nucleic acid content between the lactating mammary tissue (MT) and MFG are quite similar but discrepancies exist in their miRNA content. Our objective was to identify the origin of these discrepancies, and to evaluate the existence of a possible mechanism sorting miRNAs that will or will not be exported from the mammary epithelial cells (MECs) in bovine MFG. miR-125b-5p, miR-126-3p, miR-141-3p, and miR-204-5p, chosen on the basis of their abundance in the MT, were quantified using RT-qPCR in lactating cow MT, MFG, and laser capture-microdissected MECs. Two miRNAs (miR-125b-5p and miR-141-3p) were detected in the MT as well as in MFG and MECs. miR-204-5p was detected only in the MT, suggesting that it is very likely expressed in a cell type other than MECs. miR-126-3p was detected both in the MT and in MECs but not in MFG, suggesting a targeting mechanism for miRNAs in MECs. These results highlights differences in miRNA content between MECs and MFG may be due to a possibly not random mechanism for loading MFG with miRNA cargos that could involve a variable distribution in MECs or a sorting mechanism.

Extracellular vesicles produced by human and animal Staphylococcus aureus strains share a highly conserved core proteome.

Staphylococcus aureus is an important opportunistic pathogen of humans and animals. It produces extracellular vesicles (EVs) that are involved in cellular communication and enable inter-kingdom crosstalk, the delivery of virulence factors and modulation of the host immune response. The protein content of EVs determines their biological functions. Clarifying which proteins are selected, and how, is of crucial value to understanding the role of EVs in pathogenesis and the development of molecular delivery systems. Here, we postulated that S. aureus EVs share a common proteome containing components involved in cargo sorting. The EV proteomes of five S. aureus strains originating from human, bovine, and ovine hosts were characterised. The clustering of EV proteomes reflected the diversity of the producing strains. A total of 253 proteins were identified, 119 of which composed a core EV proteome with functions in bacterial survival, pathogenesis, and putatively in EV biology. We also identified features in the sequences of EV proteins and the corresponding genes that could account for their packaging into EVs. Our findings corroborate the hypothesis of a selective sorting of proteins into EVs and offer new perspectives concerning the roles of EVs in S. aureus pathogenesis in specific host niches.

An improved LC-MS method to profile molecular diversity and quantify the six main bovine milk proteins, including genetic and splicing variants as well as post-translationally modified isoforms.

Here we describe a method based on Liquid Chromatography coupled with Mass Spectrometry (LC-MS) that provides an accurate determination of the six main bovine milk proteins, including allelic and splicing variants, as well as isoforms resulting from post-translational modifications, with an unprecedented level of resolution. Proteins are identified from observed molecular masses in comparison with theoretical masses of intact proteins indexed in an "in-house" database that includes nearly 3000 entries. Quantification was performed either from UV (214 nm) or mass signals. Thus, up to one hundred molecules, derived from the six major milk proteins, can be identified and quantified from an individual milk sample. This powerful and reliable method, initially developed as an anchoring method to estimate the composition of the six main bovine milk proteins from MIR spectra, is transferable to several mammalian species, including small ruminants, camels, equines, rabbits, etc., for which specific mass databases are available.

Alternative splicing events expand molecular diversity of camel CSN1S2 increasing its ability to generate potentially bioactive peptides.

In a previous study on camel milk from Kazakhstan, we reported the occurrence of two unknown proteins (UP1 and UP2) with different levels of phosphorylation. Here we show that UP1 and UP2 are isoforms of camel αs2-CN (αs2-CNsv1 and αs2-CNsv2, respectively) arising from alternative splicing events. First described as a 178 amino-acids long protein carrying eight phosphate groups, the major camel αs2-CN isoform (called here αs2-CN) has a molecular mass of 21,906 Da. αs2-CNsv1, a rather frequent (35%) isoform displaying a higher molecular mass (+1,033 Da), is present at four phosphorylation levels (8P to 11P). Using cDNA-sequencing, αs2-CNsv1 was shown to be a variant arising from the splicing-in of an in-frame 27-nucleotide sequence encoding the nonapeptide ENSKKTVDM, for which the presence at the genome level was confirmed. αs2-CNsv2, which appeared to be present at 8P to 12P, was shown to include an additional decapeptide (VKAYQIIPNL) revealed by LC-MS/MS, encoded by a 3'-extension of exon 16. Since milk proteins represent a reservoir of biologically active peptides, the molecular diversity generated by differential splicing might increase its content. To evaluate this possibility, we searched for bioactive peptides encrypted in the different camel αs2-CN isoforms, using an in silico approach. Several peptides, putatively released from the C-terminal part of camel αs2-CN isoforms after in silico digestion by proteases from the digestive tract, were predicted to display anti-bacterial and antihypertensive activities.

The main WAP isoform usually found in camel milk arises from the usage of an improbable intron cryptic splice site in the precursor to mRNA in which a GC-AG intron occurs.

BACKGROUND: Whey acidic protein (WAP) is a major protein identified in the milk of several mammalian species with cysteine-rich domains known as four-disulfide cores (4-DSC). The organization of the eutherian WAP genes is highly conserved through evolution. It has been proposed that WAP could play an important role in regulating the proliferation of mammary epithelial cells. A bacteriostatic activity was also reported. Conversely to the other mammalian species expressing WAP in their milk, camel WAP contains 4 additional amino acid residues at the beginning of the second 4-DSC domain, introducing a phosphorylation site. The aim of this study was to elucidate the origin of this specificity, which possibly impacts its physiological functions. RESULTS: Using LC-ESI-MS, we identified in Camelus bactrianus from Kazakhstan a phosphorylated whey protein, exhibiting a molecular mass (12,596 Da), 32 Da higher than the original WAP (12,564 Da) and co-eluting with WAP. cDNA sequencing revealed a transition G/A, which modifies an amino acid residue of the mature protein (V12 M), accounting for the mass difference observed between WAP genetic variants. We also report the existence of two splicing variants of camel WAP precursors to mRNA, arising from an alternative usage of the canonical splice site recognized as such in the other mammalian species. However, the major camel WAP isoform results from the usage of an unlikely intron cryptic splice site, extending camel exon 3 upstream by 12-nucleotides encoding the 4 additional amino acid residues (VSSP) in which a potentially phosphorylable Serine residue occurs. Combining protein and cDNA sequences with genome data available (NCBI database), we report another feature of the camel WAP gene which displays a very rare GC-AG type intron. This result was confirmed by sequencing a genomic DNA fragment encompassing exon 3 to exon 4, suggesting for the GC donor site a compensatory effect in terms of consensus at the acceptor exon position. CONCLUSIONS: Combining proteomic and molecular biology approaches we report: the characterization of a new genetic variant of camel WAP, the usage of an unlikely intron cryptic splice site, and the occurrence of an extremely rare GC-AG type of intron.

A New High Hydrostatic Pressure Process to Assure the Microbial Safety of Human Milk While Preserving the Biological Activity of Its Main Components.

Background: The main process used to pasteurize human milk is the low-temperature, long-time Holder method. More recently, the high-temperature, short-time method has been investigated. Both processes lead to the appropriate inactivation of vegetative bacterial forms but are ineffective against bacterial spores. Research Aims/Questions: We aimed to accomplish two main objectives: inactivation of all pathogens, including spores; and preservation of the activity of milk components. Design/Methods: Recently, a novel high-hydrostatic pressure process has been developed by HPBioTECH. Using the same raw human milk samples, we compared the effects of this method with those of the Holder method on vegetative and spore forms of pathogens and on bioactive components (lipase activity, immunoproteins). Results: Two main microbial strains were selected: Staphylococcus aureus (as a reference for vegetative forms) and Bacillus cereus (as a reference for spores). Use of the high-hydrostatic pressure process led to microbial decontamination of 6 log for both S. aureus and B. cereus. Additionally, the bioactivity of the main components of human milk was preserved, with activities of lipase, α-lactalbumin, casein, lysozyme, lactoferrin, and sIgA of ~80, 96-99, 98-100, 95-100, 93-97, and 63-64%, respectively. Conclusions: Use of this novel high-hydrostatic pressure process to generate microbiologically safe human milk may provide important benefits for preterm infants, including improved assimilation of human milk (leading increased weight gain) and improved resistance to infections. Because 10% of all human milk collected is contaminated by B. cereus, use of this method will also prevent waste.

Host factors determine the evolution of infection with Staphylococcus aureus to gangrenous mastitis in goats.

Staphylococcus aureus is the major cause of very severe mastitis of dairy goats. The initial objective of our study was to fine-tune an experimental model of infection of the goat mammary gland with two strains of S. aureus and two lines of goats (low and high somatic cell score lines). Following the challenge, the 10 infected goats divided in two clear-cut severity groups, independently of the S. aureus strain and the goat line. Five goats developed very severe mastitis (of which four were gangrenous) characterized by uncontrolled infection (UI group), whereas the other five kept the infection under control (CI group). The outcome of the infection was determined by 18 h post-infection (hpi), as heralded by the bacterial milk concentration at 18 hpi: more than 107/mL in the UI group, about 106/mL in the CI group. Leukocyte recruitment and composition did not differ between the groups, but the phagocytic killing at 18 hpi efficiency did. Contributing factors involved milk concentrations of α-toxin and LukMF' leukotoxin, but not early expression of the genes encoding the pentraxin PTX3, the cytokines IL-1α and IL-1ß, and the chemokines IL-8 and CCL5. Concentrations of TNF-α, IFN-γ, IL-17A, and IL-22 rose sharply in the milk of UI goats when infection was out of control. The results indicate that defenses mobilized by the mammary gland at an early stage of infection were essential to prevent staphylococci from reaching critical concentrations. Staphylococcal exotoxin production appeared to be a consequent event inducing the evolution to gangrenous mastitis.

Combining different proteomic approaches to resolve complexity of the milk protein fraction of dromedary, Bactrian camels and hybrids, from different regions of Kazakhstan.

Nutritional suitability of milk is not only related to gross composition, but is also strongly affected by the microheterogeniety of the protein fraction. Hence, to go further into the evaluation of the potential suitability of non-bovine milks in human/infant nutrition it is necessary to have a detailed characterization of their protein components. Combining proven proteomic approaches (SDS-PAGE, LC-MS/MS and LC-ESI-MS) and cDNA sequencing, we provide here in depth characterization of the milk protein fraction of dromedary and Bactrian camels, and their hybrids, from different regions of Kazakhstan. A total 391 functional groups of proteins were identified from 8 camel milk samples. A detailed characterization of 50 protein molecules, relating to genetic variants and isoforms arising from post-translational modifications and alternative splicing events, belonging to nine protein families (κ-, αs1-, αs2-, ß-; and γ-CN, WAP, α-LAC, PGRP, CSA/LPO) was achieved by LC-ESI-MS. The presence of two unknown proteins UP1 (22,939 Da) and UP2 (23,046 Da) was also reported as well as the existence of a ß-CN short isoform (946 Da lighter than the full-length ß-CN), arising very likely in both genetic variants (A and B) from proteolysis by plasmin. In addition, we report, for the first time to our knowledge, the occurrence of a αs2-CN phosphorylation isoform with 12P groups within two recognition motifs, suggesting thereby the existence of two kinase systems involved in the phosphorylation of caseins in the mammary gland. Finally, we demonstrate that genetic variants, which hitherto seemed to be species- specific (e.g. ß-CN A for Bactrian and ß-CN B for dromedary), are in fact present both in Camel dromedarius and C. bactrianus.

Within-breed and multi-breed GWAS on imputed whole-genome sequence variants reveal candidate mutations affecting milk protein composition in dairy cattle.

BACKGROUND: Genome-wide association studies (GWAS) were performed at the sequence level to identify candidate mutations that affect the expression of six major milk proteins in Montbéliarde (MON), Normande (NOR), and Holstein (HOL) dairy cattle. Whey protein (α-lactalbumin and ß-lactoglobulin) and casein (αs1, αs2, ß, and κ) contents were estimated by mid-infrared (MIR) spectrometry, with medium to high accuracy (0.59 ≤ R2 ≤ 0.92), for 848,068 test-day milk samples from 156,660 cows in the first three lactations. Milk composition was evaluated as average test-day measurements adjusted for environmental effects. Next, we genotyped a subset of 8080 cows (2967 MON, 2737 NOR, and 2306 HOL) with the BovineSNP50 Beadchip. For each breed, genotypes were first imputed to high-density (HD) using HD single nucleotide polymorphisms (SNPs) genotypes of 522 MON, 546 NOR, and 776 HOL bulls. The resulting HD SNP genotypes were subsequently imputed to the sequence level using 27 million high-quality sequence variants selected from Run4 of the 1000 Bull Genomes consortium (1147 bulls). Within-breed, multi-breed, and conditional GWAS were performed. RESULTS: Thirty-four distinct genomic regions were identified. Three regions on chromosomes 6, 11, and 20 had very significant effects on milk composition and were shared across the three breeds. Other significant effects, which partially overlapped across breeds, were found on almost all the autosomes. Multi-breed analyses provided a larger number of significant genomic regions with smaller confidence intervals than within-breed analyses. Combinations of within-breed, multi-breed, and conditional analyses led to the identification of putative causative variants in several candidate genes that presented significant protein-protein interactions enrichment, including those with previously described effects on milk composition (SLC37A1, MGST1, ABCG2, CSN1S1, CSN2, CSN1S2, CSN3, PAEP, DGAT1, AGPAT6) and those with effects reported for the first time here (ALPL, ANKH, PICALM). CONCLUSIONS: GWAS applied to fine-scale phenotypes, multiple breeds, and whole-genome sequences seems to be effective to identify candidate gene variants. However, although we identified functional links between some candidate genes and milk phenotypes, the causality between candidate variants and milk protein composition remains to be demonstrated. Nevertheless, the identification of potential causative mutations that underlie milk protein composition may have immediate applications for improvements in cheese-making.

Characterisation of adipocyte-derived extracellular vesicle subtypes identifies distinct protein and lipid signatures for large and small extracellular vesicles.

Extracellular vesicles (EVs) are biological vectors that can modulate the metabolism of target cells by conveying signalling proteins and genomic material. The level of EVs in plasma is significantly increased in cardiometabolic diseases associated with obesity, suggesting their possible participation in the development of metabolic dysfunction. With regard to the poor definition of adipocyte-derived EVs, the purpose of this study was to characterise both qualitatively and quantitatively EVs subpopulations secreted by fat cells. Adipocyte-derived EVs were isolated by differential centrifugation of conditioned media collected from 3T3-L1 adipocytes cultured for 24 h in serum-free conditions. Based on morphological and biochemical properties, as well as quantification of secreted EVs, we distinguished two subpopulations of adipocyte-derived EVs, namely small extracellular vesicles (sEVs) and large extracellular vesicles (lEVs). Proteomic analyses revealed that lEVs and sEVs exhibit specific protein signatures, allowing us not only to define novel markers of each population, but also to predict their biological functions. Despite similar phospholipid patterns, the comparative lipidomic analysis performed on these EV subclasses revealed a specific cholesterol enrichment of the sEV population, whereas lEVs were characterised by high amounts of externalised phosphatidylserine. Enhanced secretion of lEVs and sEVs is achievable following exposure to different biological stimuli related to the chronic low-grade inflammation state associated with obesity. Finally, we demonstrate the ability of primary murine adipocytes to secrete sEVs and lEVs, which display physical and biological characteristics similar to those described for 3T3-L1. Our study provides additional information and elements to define EV subtypes based on the characterisation of adipocyte-derived EV populations. It also underscores the need to distinguish EV subpopulations, through a combination of multiple approaches and markers, since their specific composition may cause distinct metabolic responses in recipient cells and tissues.

Thermal manipulation of the chicken embryo triggers differential gene expression in response to a later heat challenge.

BACKGROUND: Meat type chickens have limited capacities to cope with high environmental temperatures, this sometimes leading to mortality on farms and subsequent economic losses. A strategy to alleviate this problem is to enhance adaptive capacities to face heat exposure using thermal manipulation (TM) during embryogenesis. This strategy was shown to improve thermotolerance during their life span. The aim of this study was to determine the effects of TM (39.5 °C, 12 h/24 vs 37.8 °C from d7 to d16 of embryogenesis) and of a subsequent heat challenge (32 °C for 5 h) applied on d34 on gene expression in the Pectoralis major muscle (PM). A chicken gene expression microarray (8 × 60 K) was used to compare muscle gene expression profiles of Control (C characterized by relatively high body temperatures, Tb) and TM chickens (characterized by a relatively low Tb) reared at 21 °C and at 32 °C (CHC and TMHC, respectively) in a dye-swap design with four comparisons and 8 broilers per treatment. Real-time quantitative PCR (RT-qPCR) was subsequently performed to validate differential expression in each comparison. Gene ontology, clustering and network building strategies were then used to identify pathways affected by TM and heat challenge. RESULTS: Among the genes differentially expressed (DE) in the PM (1.5 % of total probes), 28 were found to be differentially expressed between C and TM, 128 between CHC and C, and 759 between TMHC and TM. No DE gene was found between TMHC and CHC broilers. The majority of DE genes analyzed by RT-qPCR were validated. In the TM/C comparison, DE genes were involved in energy metabolism and mitochondrial function, cell proliferation, vascularization and muscle growth; when comparing heat-exposed chickens to their own controls, TM broilers developed more specific pathways than C, especially involving genes related to metabolism, stress response, vascularization, anti-apoptotic and epigenetic processes. CONCLUSIONS: This study improved the understanding of the long-term effects of TM on PM muscle. TM broilers displaying low Tb may have lower metabolic intensity in the muscle, resulting in decreased metabolic heat production, whereas modifications in vascularization may enhance heat loss. These specific changes could in part explain the better adaptation of TM broilers to heat.

Protein profiling of isolated uterine AA amyloidosis causing fetal death in goats.

Pathologic amyloid accumulates in the CNS or in peripheral organs, yet the mechanism underlying the targeting of systemic amyloid deposits is unclear. Serum amyloid A (SAA) 1 and 2 are produced predominantly by the liver and form amyloid most commonly in the spleen, liver, and kidney. In contrast, SAA3 is produced primarily extrahepatically and has no causal link to amyloid formation. Here, we identified 8 amyloidosis cases with amyloid composed of SAA3 expanding the uterine wall of goats with near-term fetuses. Uterine amyloid accumulated in the endometrium, only at the site of placental attachment, compromising maternal-fetal gas and nutrient exchange and leading to fetal ischemia and death. No other organ contained amyloid. SAA3 mRNA levels in the uterine endometrium were as high as SAA2 in the liver, yet mass spectrometry of the insoluble uterine peptides identified SAA3 as the predominant protein, and not SAA1 or SAA2. These findings suggest that high local SAA3 production led to deposition at this unusual site. Although amyloid A (AA) amyloid deposits typically consist of an N-terminal fragment of SAA1 or SAA2, here, abundant C-terminal peptides indicated that the uterine amyloid was largely composed of full-length SAA3. The exclusive deposition of SAA3 amyloid in the uterus, together with elevated uterine SAA3 transcripts, suggests that the uterine amyloid deposits were due to locally produced SAA3. This is the first report of SAA3 as a cause of amyloidosis and of AA amyloid deposited exclusively in the uterus.

Early bovine embryos regulate oviduct epithelial cell gene expression during in vitro co-culture.

In mammals, the oviduct may participate to the regulation of early embryo development. In vitro co-culture of early bovine embryos with bovine oviduct epithelial cells (BOEC) has been largely used to mimic the maternal environment. However, the mechanisms of BOEC action have not been clearly elucidated yet. The aim of this study was to determine the response of BOEC cultures to the presence of developing bovine embryos. A 21,581-element bovine oligonucleotide array was used compare the gene expression profiles of confluent BOEC cultured for 8 days with or without embryos. This study revealed 34 differentially expressed genes (DEG). Of these 34 genes, IFI6, ISG15, MX1, IFI27, IFI44, RSAD2, IFITM1, EPSTI1, USP18, IFIT5, and STAT1 expression increased to the greatest extent due to the presence of embryos with a major impact on antiviral and immune response. Among the mRNAs at least 25 are already described as induced by interferons. In addition, transcript levels of new candidate genes involved in the regulation of transcription, modulation of the maternal immune system and endometrial remodeling were found to be increased. We selected 7 genes and confirmed their differential expression by quantitative RT-PCR. The immunofluorescence imaging of cellular localization of STAT1 protein in BOEC showed a nuclear translocation in the presence of embryos, suggesting the activation of interferon signaling pathway. This first systematic study of BOEC transcriptome changes in response to the presence of embryos in cattle provides some evidences that these cells are able to adapt their transcriptomic profile in response to embryo signaling.

Comparison of five different RNA sources to examine the lactating bovine mammary gland transcriptome using RNA-Sequencing.

The objective of this study was to examine five different sources of RNA, namely mammary gland tissue (MGT), milk somatic cells (SC), laser microdissected mammary epithelial cells (LCMEC), milk fat globules (MFG) and antibody-captured milk mammary epithelial cells (mMEC) to analyze the bovine mammary gland transcriptome using RNA-Sequencing. Our results provide a comparison between different sampling methods (invasive and non-invasive) to define the transcriptome of mammary gland tissue and milk cells. This information will be of value to investigators in choosing the most appropriate sampling method for different research applications to study specific physiological states during lactation. One of the simplest procedures to study the transcriptome associated with milk appears to be the isolation of total RNA directly from SC or MFG released into milk during lactation. Our results indicate that the SC and MFG transcriptome are representative of MGT and LCMEC and can be used as effective and alternative samples to study mammary gland expression without the need to perform a tissue biopsy.

Characterization of goat milk lactoferrin N-glycans and comparison with the N-glycomes of human and bovine milk.

Numerous milk components, such as lactoferrin, are recognized as health-promoting compounds. A growing body of evidence suggests that glycans could mediate lactoferrin's bioactivity. Goat milk lactoferrin is a candidate for infant formula supplementation because of its high homology with its human counterpart. The aim of this study was to characterize the glycosylation pattern of goat milk lactoferrin. After the protein was isolated from milk by affinity chromatography, N-glycans were enzymatically released and a complete characterization of glycan composition was carried out by advanced MS. The glycosylation of goat milk lactoferrin was compared with that of human and bovine milk glycoproteins. Nano-LC-Chip-Q-TOF MS data identified 65 structures, including high mannose, hybrid, and complex N-glycans. Among the N-glycan compositions, 37% were sialylated and 34% were fucosylated. The results demonstrated the existence of similar glycans in human and goat milk but also identified novel glycans in goat milk that were not present in human milk. These data suggest that goat milk could be a source of bioactive compounds, including lactoferrin that could be used as functional ingredients for food products beneficial to human nutrition.

Combining proteomic tools to characterize the protein fraction of llama (Lama glama) milk.

Llamas belong to the Camelidae family along with camels. While dromedary camel milk has been broadly characterized, data on llama milk proteins are scarce. The objective of this study was thus to investigate the protein composition of llama milk. Skimmed llama milk proteins were first characterized by a 2D separation technique coupling RP-HPLC in the first dimension with SDS-PAGE in the second dimension (RP-HPLC/SDS-PAGE). Llama milk proteins, namely caseins (αs1 -, αs2 -, ß-, and κ-caseins), α-lactalbumin, lactoferrin, and serum albumin, were identified using PMF. Llama milk proteins were also characterized by online LC-ESI-MS analysis. This approach allowed attributing precise molecular masses for most of the previously MS-identified llama milk proteins. Interestingly, α-lactalbumin exhibits distinct chromatographic behaviors between llama and dromedary camel milk. De novo sequencing of the llama α-lactalbumin protein by LC coupled with MS/MS (LC-MS/MS) showed the occurrence of two amino acid substitutions (R62L/I and K89L/I) that partly explained the higher hydrophobicity of llama α-lactalbumin compared with its dromedary counterpart. Taken together, these results provide for the first time a thorough description of the protein fraction of Lama glama milk.