Characterisation of fatty acyl reductases of sunflower (Helianthus annuus L.) seed.

Long and very long chain fatty alcohols are produced from their corresponding acyl-CoAs through the activity of fatty acyl reductases (FARs). Fatty alcohols are important components of the cuticle that protects aerial plant organs, and they are metabolic intermediates in the synthesis of the wax esters in the hull of sunflower (Helianthus annuus) seeds. Genes encoding 4 different FARs (named HaFAR2, HaFAR3, HaFAR4 and HaFAR5) were identified using BLAST, and studies showed that four of the genes were expressed in seed hulls. In this study, the structure and location of sunflower FAR proteins were determined. They were also expressed exogenously in Saccharomyces cerevisiae to evaluate their substrate specificity based on the fatty alcohols synthesized by the transformed yeasts. Three of the four enzymes tested showed activity in yeast. HaFAR3 produced C18, C20 and C22 saturated alcohols, whereas HaFAR4 and HaFAR5 produced C24 and C26 saturated alcohols. The involvement of these genes in the synthesis of sunflower seed wax esters was addressed by considering the results obtained.

Genomics of predictive radiation mutagenesis in oilseed rape: modifying seed oil composition.

Rapeseed is a crop of global importance but there is a need to broaden the genetic diversity available to address breeding objectives. Radiation mutagenesis, supported by genomics, has the potential to supersede genome editing for both gene knockout and copy number increase, but detailed knowledge of the molecular outcomes of radiation treatment is lacking. To address this, we produced a genome re-sequenced panel of 1133 M2 generation rapeseed plants and analysed large-scale deletions, single nucleotide variants and small insertion-deletion variants affecting gene open reading frames. We show that high radiation doses (2000 Gy) are tolerated, gamma radiation and fast neutron radiation have similar impacts and that segments deleted from the genomes of some plants are inherited as additional copies by their siblings, enabling gene dosage decrease. Of relevance for species with larger genomes, we showed that these large-scale impacts can also be detected using transcriptome re-sequencing. To test the utility of the approach for predictive alteration of oil fatty acid composition, we produced lines with both decreased and increased copy numbers of Bna.FAE1 and confirmed the anticipated impacts on erucic acid content. We detected and tested a 21-base deletion expected to abolish function of Bna.FAD2.A5, for which we confirmed the predicted reduction in seed oil polyunsaturated fatty acid content. Our improved understanding of the molecular effects of radiation mutagenesis will underpin genomics-led approaches to more efficient introduction of novel genetic variation into the breeding of this crop and provides an exemplar for the predictive improvement of other crops.

Tackling functional redundancy of Arabidopsis fatty acid elongase complexes.

Very-long-chain fatty acids (VLCFA) are precursors for various lipids playing important physiological and structural roles in plants. Throughout plant tissues, VLCFA are present in multiple lipid classes essential for membrane homeostasis, and also stored in triacylglycerols. VLCFA and their derivatives are also highly abundant in lipid barriers, such as cuticular waxes in aerial epidermal cells and suberin monomers in roots. VLCFA are produced by the fatty acid elongase (FAE), which is an integral endoplasmic reticulum membrane multi-enzymatic complex consisting of four core enzymes. The 3-ketoacyl-CoA synthase (KCS) catalyzes the first reaction of the elongation and determines the chain-length substrate specificity of each elongation cycle, whereas the other three enzymes have broad substrate specificities and are shared by all FAE complexes. Consistent with the co-existence of multiple FAE complexes, performing sequential and/or parallel reactions to produce the broad chain-length-range of VLCFA found in plants, twenty-one KCS genes have been identified in the genome of Arabidopsis thaliana. Using CRISPR-Cas9 technology, we established an expression platform to reconstitute the different Arabidopsis FAE complexes in yeast. The VLCFA produced in these yeast strains were analyzed in detail to characterize the substrate specificity of all KCS candidates. Additionally, Arabidopsis candidate proteins were transiently expressed in Nicotiana benthamiana leaves to explore their activity and localization in planta. This work sheds light on the genetic and biochemical redundancy of fatty acid elongation in plants.

Analysis of Free and Esterified Sterol Content and Composition in Seeds Using GC and ESI-MS/MS.

Total sterol content and composition in plant tissues can be easily determined by gas chromatography (GC) after saponification of the total lipid extract. However, in oleogenic tissues a significant proportion of the sterol is esterified to fatty acids, with GC methodologies unable to provide information about the proportion and the molecular species composition of intact steryl esters (SEs). Here we describe an electrospray ionization-tandem mass spectrometry (ESI-MS/MS) and Multiple Reaction Monitoring (MRM) method which, in parallel with GC analysis, allows for the accurate determination of both free and esterified sterol content and composition in seeds. After extraction of seed oil with hexane, free sterols are derivatized with undecanoyl chloride, total steryl esters are then purified from triacylglycerol (TAG) by liquid chromatography, infused and ionized as ammonium adducts, with molecular species identified and quantified by fragmentation in the presence of internal standards.

DES2 is a fatty acid Δ11 desaturase capable of synthesizing palmitvaccenic acid in the arbuscular mycorrhizal fungus Rhizophagus irregularis.

Arbuscular mycorrhizal (AM) fungi are oleaginous organisms, and the most abundant fatty acyl moiety usually found in their lipids is palmitvaccenic acid (16:1Δ11cis ). However, it is not known how this uncommon fatty acid species is made. Here, we have cloned two homologues of lepidopteran fatty acyl-coenzyme A Δ11 desaturases from the AM fungus Rhizophagus irregularis. Both enzymes, DES1 and DES2, are expressed in intraradical mycelium and can complement the unsaturated fatty acid-requiring auxotrophic growth phenotype of the Saccharomyces cerevisiae ole1Δ mutant. DES1 expression leads almost exclusively to oleic acid (18:1Δ9cis ) production, whereas DES2 expression results in the production of 16:1Δ11cis and vaccenic acid (18:1Δ11cis ). DES2 therefore encodes a Δ11 desaturase that is likely to be responsible for the synthesis of 16:1Δ11cis in R. irregularis.

Piglet weight gain during the first two weeks of lactation influences the immune system development.

The aim of this study was to investigate the effect of the piglet growth during the first week of life on ileal expression of genes and on development of the immune system. Eight litters adjusted to 12 piglets were used. Within each litter, the piglet that showed the lowest weight gain (LWG; n = 8) and the one that showed the highest weight gain (HWG; n = 8) in their first week of life were enrolled. Peripheral blood mononuclear cells (PBMC) were isolated on days 8 and 16 to characterize cellular population profiles and to assess ex-vivo secretion of interleukin-10 (IL-10), IL-6 and tumor necrosis factor-α (TNF-α). On day 16, piglets were euthanized and ileum samples were collected to extract RNA for microarray analysis and gene expression by qPCR. As expected, growth performance of LWG piglet was impaired compared to HWG piglets (P < 0.05). From day 8 to 16, the percentage of CD21+ B cells significantly increased in blood of heavier HWG piglets while the percentage remained constant in smaller LWG piglets (P weight x day = 0.01). For the CD4+CD8α- Th cells, a marked increase was observed in LWG piglets from 8 to 16 days of age (P = 0.002) whereas no significant change occurred in HWG piglets. Percentages of CD14+ monocytes and other MHC-II+ cells were respectively higher and lower on day 8 compared to day 16 for both groups of piglets (P < 0.01). On day 8, LPS-activated PBMC from LWG piglets produced less IL-6 compared to HWG piglets (P < 0.05). Microarray analysis of gene expression in piglets' ileum tissue indicated that several genes involed in defense response and response to oxidative stress were modulated differently in LWG compared to HWG. Gene analysis by Q-PCR confirmed microarray results and revealed that IL-10, SOD1, NOS2, NOD2, TLR4, TLR9, CD40 and CD74 expressions were significantly decreased (P < 0.05) in LWG in comparison to HWG piglets, while MYD88 and NFkBiA showed a tendency to decrease (0.05 ≤ P < 0.07). These results suggest that birth weight and milk intake affect the growth performances and the development of immunity by modulating the expression of genes associated with immunity and oxidative stress in piglets' intestinal tissue, and by affecting the leukocyte populations involved in innate and cell-mediated immunity in nursing piglets. Therefore, impaired development of immune system in LWG piglets might have an impact on their resistance to infections later in life.

New insights in the composition of wax and sterol esters in common and mutant sunflower oils revealed by ESI-MS/MS.

Wax esters (WEs) and steryl esters (SEs) are minor components of sunflower oils formed by the esterification of long chain fatty alcohols and sterols to fatty acids. These compounds have similar carbon numbers and polarities making them difficult to separate using conventional chromatographic methods. In this study, electrospray ionisation-tandem mass spectrometry (ESI-MS/MS) allowed the rapid and accurate profiling of WEs and SEs acyl moieties in total ester fractions of common and mutant sunflower oils with different fatty acid profiles. The acyl composition of both WEs and SEs partially reflected that of the oil and the high oleic background displayed the lowest level of crystallisable waxes. ESI-MS/MS complemented by GC-MS analyses revealed that SEs contain 17-30% of previously unreported methylsterol moieties. We demonstrated that these compounds are overlooked by official sterol analytical methods which may have consequences for quality control and authentication of vegetable oils prior to commercialisation.

Tailoring the composition of novel wax esters in the seeds of transgenic Camelina sativa through systematic metabolic engineering.

The functional characterization of wax biosynthetic enzymes in transgenic plants has opened the possibility of producing tailored wax esters (WEs) in the seeds of a suitable host crop. In this study, in addition to systematically evaluating a panel of WE biosynthetic activities, we have also modulated the acyl-CoA substrate pool, through the co-expression of acyl-ACP thioesterases, to direct the accumulation of medium-chain fatty acids. Using this combinatorial approach, we determined the additive contribution of both the varied acyl-CoA pool and biosynthetic enzyme substrate specificity to the accumulation of non-native WEs in the seeds of transgenic Camelina plants. A total of fourteen constructs were prepared containing selected FAR and WS genes in combination with an acyl-ACP thioesterase. All enzyme combinations led to the successful production of wax esters, of differing compositions. The impact of acyl-CoA thioesterase expression on wax ester accumulation varied depending on the substrate specificity of the WS. Hence, co-expression of acyl-ACP thioesterases with Marinobacter hydrocarbonoclasticus WS and Marinobacter aquaeolei FAR resulted in the production of WEs with reduced chain lengths, whereas the co-expression of the same acyl-ACP thioesterases in combination with Mus musculus WS and M. aquaeolei FAR had little impact on the overall final wax composition. This was despite substantial remodelling of the acyl-CoA pool, suggesting that these substrates were not efficiently incorporated into WEs. These results indicate that modification of the substrate pool requires careful selection of the WS and FAR activities for the successful high accumulation of these novel wax ester species in Camelina seeds.

A Palmitic Acid Elongase Affects Eicosapentaenoic Acid and Plastidial Monogalactosyldiacylglycerol Levels in Nannochloropsis.

Nannochloropsis species are oleaginous eukaryotes containing a plastid limited by four membranes, deriving from a secondary endosymbiosis. In Nannochloropsis, thylakoid lipids, including monogalactosyldiacylglycerol (MGDG), are enriched in eicosapentaenoic acid (EPA). The need for EPA in MGDG is not understood. Fatty acids are de novo synthesized in the stroma, then converted into very-long-chain polyunsaturated fatty acids (FAs) at the endoplasmic reticulum (ER). The production of MGDG relies therefore on an EPA supply from the ER to the plastid, following an unknown process. We identified seven elongases and five desaturases possibly involved in EPA production in Nannochloropsis gaditana Among the six heterokont-specific saturated FA elongases possibly acting upstream in this pathway, we characterized the highly expressed isoform Δ0-ELO1 Heterologous expression in yeast (Saccharomyces cerevisiae) showed that NgΔ0-ELO1 could elongate palmitic acid. Nannochloropsis Δ0-elo1 mutants exhibited a reduced EPA level and a specific decrease in MGDG In NgΔ0-elo1 lines, the impairment of photosynthesis is consistent with a role of EPA-rich MGDG in nonphotochemical quenching control, possibly providing an appropriate MGDG platform for the xanthophyll cycle. Concomitantly with MGDG decrease, the level of triacylglycerol (TAG) containing medium chain FAs increased. In Nannochloropsis, part of EPA used for MGDG production is therefore biosynthesized by a channeled process initiated at the elongation step of palmitic acid by Δ0-ELO1, thus acting as a committing enzyme for galactolipid production. Based on the MGDG/TAG balance controlled by Δ0-ELO1, this study also provides novel prospects for the engineering of oleaginous microalgae for biotechnological applications.

Dual Fatty Acid Elongase Complex Interactions in Arabidopsis.

Very long chain fatty acids (VLCFAs) are involved in plant development and particularly in several cellular processes such as membrane trafficking, cell division and cell differentiation. However, the precise role of VLCFAs in these different cellular processes is still poorly understood in plants. In order to identify new factors associated with the biosynthesis or function of VLCFAs, a yeast multicopy suppressor screen was carried out in a yeast mutant strain defective for fatty acid elongation. Loss of function of the elongase 3 hydroxyacyl-CoA dehydratase PHS1 in yeast and PASTICCINO2 in plants prevents growth and induces cytokinesis defects. PROTEIN TYROSIN PHOSPHATASE-LIKE (PTPLA) previously characterized as an inactive dehydratase was able to restore yeast phs1 growth and VLCFAs elongation but not the plant pas2-1 defects. PTPLA interacted with elongase subunits in the Endoplasmic Reticulum (ER) and its absence induced the accumulation of 3-hydroxyacyl-CoA as expected from a dehydratase involved in fatty acid (FA) elongation. However, loss of PTPLA function increased VLCFA levels, an effect that was dependent on the presence of PAS2 indicating that PTPLA activity repressed FA elongation. The two dehydratases have specific expression profiles in the root with PAS2, mostly restricted to the endodermis, while PTPLA was confined in the vascular tissue and pericycle cells. Comparative ectopic expression of PTPLA and PAS2 in their respective domains confirmed the existence of two independent elongase complexes based on PAS2 or PTPLA dehydratase that are functionally interacting.

ACYL-ACYL CARRIER PROTEIN DESATURASE2 and 3 Are Responsible for Making Omega-7 Fatty Acids in the Arabidopsis Aleurone.

Omega-7 monounsaturated fatty acids (ω-7s) are specifically enriched in the aleurone of Arabidopsis (Arabidopsis thaliana) seeds. We found significant natural variation in seed ω-7 content and used a Multiparent Advanced Generation Inter-Cross population to fine-map a major quantitative trait loci to a region containing ACYL-ACYL CARRIER PROTEIN DESATURASE1 (AAD1) and AAD3 We found that AAD3 expression is localized to the aleurone where mutants show an approximately 50% reduction in ω-7 content. By contrast, AAD1 is localized to the embryo where mutants show a small reduction in ω-9 content. Enzymatic analysis has previously shown that AAD family members possess both stearoyl- and palmitoyl-ACP Δ(9) desaturase activity, including the predominant isoform SUPPRESSOR OF SALICYLIC ACID INSENSITIVE2. However, aad3 ssi2 aleurone contained the same amount of ω-7s as aad3 Within the AAD family, AAD3 shares the highest degree of sequence similarity with AAD2 and AAD4. Mutant analysis showed that AAD2 also contributes to ω-7 production in the aleurone, and aad3 aad2 exhibits an approximately 85% reduction in ω-7s Mutant analysis also showed that FATTY ACID ELONGASE1 is required for the production of very long chain ω-7s in the aleurone. Together, these data provide genetic evidence that the ω-7 pathway proceeds via Δ(9) desaturation of palmitoyl-ACP followed by elongation of the product. Interestingly, significant variation was also identified in the ω-7 content of Brassica napus aleurone, with the highest level detected being approximately 47% of total fatty acids.

The Zinc-Finger Protein SOP1 Is Required for a Subset of the Nuclear Exosome Functions in Arabidopsis.

Correct gene expression requires tight RNA quality control both at transcriptional and post-transcriptional levels. Using a splicing-defective allele of PASTICCINO2 (PAS2), a gene essential for plant development, we isolated suppressor mutations modifying pas2-1 mRNA profiles and restoring wild-type growth. Three suppressor of pas2 (sop) mutations modified the degradation of mis-spliced pas2-1 mRNA species, allowing the synthesis of a functional protein. Cloning of the suppressor mutations identified the core subunit of the exosome SOP2/RRP4, the exosome nucleoplasmic cofactor SOP3/HEN2 and a novel zinc-finger protein SOP1 that colocalizes with HEN2 in nucleoplasmic foci. The three SOP proteins counteract post-transcriptional (trans)gene silencing (PTGS), which suggests that they all act in RNA quality control. In addition, sop1 mutants accumulate some, but not all of the misprocessed mRNAs and other types of RNAs that are observed in exosome mutants. Taken together, our data show that SOP1 is a new component of nuclear RNA surveillance that is required for the degradation of a specific subset of nuclear exosome targets.

C1 Metabolism Inhibition and Nitrogen Deprivation Trigger Triacylglycerol Accumulation in Arabidopsis thaliana Cell Cultures and Highlight a Role of NPC in Phosphatidylcholine-to-Triacylglycerol Pathway.

Triacylglycerol (TAG) accumulation often occurs in growth limiting conditions such as nutrient deprivations. We analyzed and compared the lipid contents of Arabidopsis cells grown under two conditions that inhibited growth as a way to study interactions between membrane and storage lipids. In order to inhibit C1 metabolism, the first condition utilized methotrexate (MTX), a drug that inhibits methyl transfer reactions and potentially reduces Pi-choline synthesis, the polar head of phosphatidylcholine (PC). MTX-treated cells displayed a 10- to 15-fold increase in TAG compared to that found in control cells. This corresponded to a net increase of lipids as the total amount of membrane glycerolipids was minimally affected. Under this condition, PC homeostasis appeared tightly regulated and not strictly dependent on the rate of Pi-choline synthesis. The second condition we investigated involved nitrogen deprivation. Here, we observed a 40-fold increase of TAG. In these cells, the overall lipid content remained unchanged, but membrane lipids decreased by a factor of two suggesting a reduction of the membrane network and a rerouting of membrane lipids to storage lipids. Under all conditions, fatty acid (FA) analyses showed that the FA composition of TAG was comparable to that in PC, but different from that in acyl-CoA, suggesting that TAG accumulation involved PC-derived DAG moieties. In agreement, analyses by qPCR of genes coding for TAG synthesis showed a strong increase of non-specific phospholipase C (NPC) expressions, and experiments using labeled (fluorescent) PC indicated higher rates of PC-to-TAG conversion under both situations. These results highlight a role for NPC in plant cell oil production.

Deep sequencing shows microRNA involvement in bovine mammary gland adaptation to diets supplemented with linseed oil or safflower oil.

BACKGROUND: Bovine milk fat composition is responsive to dietary manipulation providing an avenue to modify the content of fatty acids and especially some specific unsaturated fatty acid (USFA) isomers of benefit to human health. MicroRNAs (miRNAs) regulate gene expression but their specific roles in bovine mammary gland lipogenesis are unclear. The objective of this study was to determine the expression pattern of miRNAs following mammary gland adaptation to dietary supplementation with 5 % linseed or safflower oil using next generation RNA-sequencing. METHODS: Twenty-four Canadian Holstein dairy cows (twelve per treatment) in mid lactation were fed a control diet (total mixed ration of corn:grass silages) for 28 days followed by a treatment period (control diet supplemented with 5 % linseed or safflower oil) of 28 days. Milk samples were collected weekly for fat and individual fatty acid determination. RNA from mammary gland biopsies harvested on day-14 (control period) and on days +7 and +28 (treatment period) from six randomly selected cows per treatment was subjected to small RNA sequencing. RESULTS: Milk fat percentage decreased significantly (P < 0.001) during treatment with the two diets as compared to the control period. The individual saturated fatty acids C4:0, C6:0, C8:0, C14:0 and C16:0 decreased significantly (P < 0.05) while five USFAs (C14:1, C18:1n11t, C20:3n3, C20:5n3 and CLA:t10c12) increased remarkably (P < 0.05) in response to both treatments. Analysis of 361 million sequence reads generated 321 known bovine miRNAs and 176 novel miRNAs. The expression of fourteen and twenty-two miRNAs was affected (P < 0.05) by linseed and safflower oil treatments, respectively. Seven miRNAs including six up-regulated (bta-miR-199c, miR-199a-3p, miR-98, miR-378, miR-148b and miR-21-5p) and one down-regulated (bta-miR-200a) were found to be regulated (P < 0.05) by both treatments, and thus considered core differentially expressed (DE) miRNAs. The gene targets of core DE miRNAs have functions related to gene expression and general cellular metabolism (P < 0.05) and are enriched in four pathways of lipid metabolism (3-phosphoinositide biosynthesis, 3-phosphoinositide degradation, D-myo-inisitol-5-phosphate metabolism and the superpathway of inositol phosphate compounds). CONCLUSION: Our results suggest that DE miRNAs in this study might be important regulators of bovine mammary lipogenesis and metabolism. The novel miRNAs identified in this study will further enrich the bovine miRNome repertoire and contribute to understanding mammary gland biology.

ECERIFERUM2-LIKE proteins have unique biochemical and physiological functions in very-long-chain fatty acid elongation.

The extension of very-long-chain fatty acids (VLCFAs) for the synthesis of specialized apoplastic lipids requires unique biochemical machinery. Condensing enzymes catalyze the first reaction in fatty acid elongation and determine the chain length of fatty acids accepted and produced by the fatty acid elongation complex. Although necessary for the elongation of all VLCFAs, known condensing enzymes cannot efficiently synthesize VLCFAs longer than 28 carbons, despite the prevalence of C28 to C34 acyl lipids in cuticular wax and the pollen coat. The eceriferum2 (cer2) mutant of Arabidopsis (Arabidopsis thaliana) was previously shown to have a specific deficiency in cuticular waxes longer than 28 carbons, and heterologous expression of CER2 in yeast (Saccharomyces cerevisiae) demonstrated that it can modify the acyl chain length produced by a condensing enzyme from 28 to 30 carbon atoms. Here, we report the physiological functions and biochemical specificities of the CER2 homologs CER2-LIKE1 and CER2-LIKE2 by mutant analysis and heterologous expression in yeast. We demonstrate that all three CER2-LIKEs function with the same small subset of condensing enzymes, and that they have different effects on the substrate specificity of the same condensing enzyme. Finally, we show that the changes in acyl chain length caused by each CER2-LIKE protein are of substantial importance for cuticle formation and pollen coat function.

Understanding and manipulating plant lipid composition: Metabolic engineering leads the way.

The manipulation of plant seed oil composition so as to deliver enhanced fatty acid compositions suitable for feed or fuel has long been a goal of metabolic engineers. Recent advances in our understanding of the flux of acyl-changes through different key metabolic pools such as phosphatidylcholine and diacylglycerol have allowed for more targeted interventions. When combined in iterative fashion with further lipidomic analyses, significant breakthroughs in our capacity to generate plants with novel oils have been achieved. Collectively these studies, working at the interface between metabolic engineering and synthetic biology, demonstrate the positive fundamental and applied outcomes derived from such research.

Transcriptome microRNA profiling of bovine mammary epithelial cells challenged with Escherichia coli or Staphylococcus aureus bacteria reveals pathogen directed microRNA expression profiles.

BACKGROUND: MicroRNAs (miRNAs) can post-transcriptionally regulate gene expression and have been shown to be critical regulators to the fine-tuning of epithelial immune responses. However, the role of miRNAs in bovine responses to E. coli and S. aureus, two mastitis causing pathogens, is not well understood. RESULTS: The global expression of miRNAs in bovine mammary epithelial cells (MAC-T cells) challenged with and without heat-inactivated Staphylococcus aureus (S. aureus) or Escherichia coli (E. coli) bacteria at 0, 6, 12, 24, and 48 hr was profiled using RNA-Seq. A total of 231 known bovine miRNAs were identified with more than 10 counts per million in at least one of 13 libraries and 5 miRNAs including bta-miR-21-5p, miR-27b, miR-22-3p, miR-184 and let-7f represented more than 50% of the abundance. One hundred and thirteen novel miRNAs were also identified and more than one third of them belong to the bta-miR-2284 family. Seventeen miRNAs were significantly (P < 0.05) differentially regulated by the presence of pathogens. E. coli initiated an earlier regulation of miRNAs (6 miRNAs differentially regulated within the first 6 hrs post challenge as compared to 1 miRNA for S. aureus) while S. aureus presented a delayed response. Five differentially expressed miRNAs (bta-miR-184, miR-24-3p, miR-148, miR-486 and let-7a-5p) were unique to E. coli while four (bta-miR-2339, miR-499, miR-23a and miR-99b) were unique to S. aureus. In addition, our study revealed a temporal differential regulation of five miRNAs (bta-miR-193a-3p, miR-423-5p, miR-30b-5p, miR-29c and miR-un116) in unchallenged cells. Target gene predictions of pathogen differentially expressed miRNAs indicate a significant enrichment in gene ontology functional categories in development/cellular processes, biological regulation as well as cell growth and death. Furthermore, target genes were significantly enriched in several KEGG pathways including immune system, signal transduction, cellular process, nervous system, development and human diseases. CONCLUSION: Using next-generation sequencing, our study identified a pathogen directed differential regulation of miRNAs in MAC-T cells with roles in immunity and development. Our study provides a further confirmation of the involvement of mammary epithelia cells in contributing to the immune response to infecting pathogens and suggests the potential of miRNAs to serve as biomarkers for diagnosis and development of control measures.