Next Article in Journal / Special Issue
Traceability of Plant Diet Contents in Raw Cow Milk Samples
Previous Article in Journal
Bovine Colostrum Supplementation During Running Training Increases Intestinal Permeability
Previous Article in Special Issue
Development and Optimization of a Fluorescent Differential Display PCR System for Analyzing the Stress Response in Lactobacillus sakei Strains
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nutrients
Volume 1
Issue 2
10.3390/nu1020235
nutrients-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Tool for Sheep Product Quality: Custom Microarrays from Public Databases

by
Silvia Bongiorni
1,*,
Giovanni Chillemi
2,
Gianluca Prosperini
2,
Susana Bueno
2,
Alessio Valentini
1 and
Lorraine Pariset
1
1
Dipartimento di Produzioni Animali, Università della Tuscia, 01100, Viterbo, Italy
2
Consorzio per le Applicazioni di Supercalcolo Per Università e Ricerca, 00185, Roma, Italy
*
Author to whom correspondence should be addressed.
Nutrients 2009, 1(2), 235-250; https://doi.org/10.3390/nu1020235
Submission received: 20 October 2009 / Accepted: 1 December 2009 / Published: 4 December 2009
(This article belongs to the Special Issue Foodomics 2009)
Browse Figures
Versions Notes

Abstract

:
Milk and dairy products are an essential food and an economic resource in many countries. Milk component synthesis and secretion by the mammary gland involve expression of a large number of genes whose nutritional regulation remains poorly defined. The purpose of this study was to gain an understanding of the genomic influence on milk quality and synthesis by comparing two sheep breeds with different milking attitude (Sarda and Gentile di Puglia) using sheep-specific microarray technology. From sheep ESTs deposited at NCBI, we have generated the first annotated microarray developed for sheep with a coverage of most of the genome.

1. Introduction

Milk and dairy products are an essential food and economic resource in many countries. Milk provides the primary source of nourishment for mammals’ offspring before their adult diet and contains the principal nutrients plus a huge number of micronutrient molecules, some of them with still unknown properties [1,2]. Therefore, the quality of milk and its control is becoming increasingly important. Milk component synthesis and secretion by the mammary gland varies dramatically across species and involves the expression of a large number of genes whose nutritional regulation remains poorly defined [3]. Nutritional genomics is an integrated science which studies gene expression to identify genetic and nutritional effects of a diet (the nutrient influence) on a single individual; while nutrigenetics seeks to understand the individual genetic differences which affect response to diet [4]. Knowledge of mammary uptake of nutrients, biosynthesis pathways, and the relation between diet and milk composition have been achieved in many studies [5,6,7,8]. Although much is known about the biochemistry of milk synthesis, the regulatory and cellular signaling systems of mammary gland are not well understood [9]. In dairy animals, mammary gland undergoes huge functional and metabolic adaptation to prepare lactogenesis. In all mammals, lactogenesis is characterized by two stages [10,11,12,13]. During the first stage (stage 01), which starts few weeks before parturition, the mammary gland differentiates for secreting colostrum and milk proteins. After parturition (stage 02), the metabolic activity increases the levels of milk production. Milk yield significantly rises during the first few weeks of lactation. During this period a well-studied set of genes, involved in milk synthesis, also increases its expression [12,14,15,16]. After the lactation peak, milk synthesis and qualified gene expression gradually decrease [14,15]. The end of milking activates the involution of the mammary gland which is characterized by epithelial cell death and by the mammary adipose tissue remodeling [17,18]. In dairy animals the nonlacting period, commonly referred to as the dry period between two lactations, is very important for milk production. A dry period of 40–60 days is necessary for optimal milk production during the next lactation [19].
In Italy, sheep is the second species in economic importance as a milk supply. Milk yield and composition, as well as lactation length, can fluctuate between breeds and within breeds. In normal sheep milk, fat ranges from 6% to 9%, protein from 4% to 7%, total solids from 17% to 21% and lactose from 4% to 6% [20]. Also other milk components implicated in human health vary considerably among breeds. Recently, Signorelli and collaborators [21] analyzed milk quality parameters and milk fatty acid profiles of three Italian breeds, Altamurana, Gentile di Puglia and Sarda, finding significant differences between breeds. The lowest content of saturated fatty acids (SFAs) was estimated in Gentile di Puglia breed, while mono-unsaturated FAs (MUFAs) were lowest in the Altamurana. No differences between breeds were found for conjugated linoleic acid (CLA) and poly-unsaturated FAs (PUFAs). Cheese quality is expected to be influenced by the differences between breeds in milk fatty acid contents [21].
The comparative analysis of some sheep breeds with different attitude to milk production could demonstrate the association between genetic variants and milk quality [21]. Candidate genes responsible for milk composition were intensively analyzed by Moioli and collaborators [22] to identify the molecular mechanisms underlying sheep milk quality. Among milk protein genes, the major effects were assessed for the αs1-casein, k-casein, β-lactoglobulin. Other important genes are those implicated in fatty acid metabolism, such as ACACA, SCD, LPL and DGAT1 [22]. However, in order to improve the overall picture, many more genes need to be deeply investigated.
Microarray technology is a powerful tool that helps to explore an organism transcriptome by measuring, in a particular cell or tissue, the expression levels of thousands of genes simultaneously. In livestock species, the microarray technology was discussed and reviewed as potential nutrigenomics tools, in the context of its economic benefits and improvement of food quality and safety in dairy and meat industries [23,24,25]. However, microarrays have been designed for very few livestock species. Moreover, the few devices so far developed, feature a largely incomplete coverage of the genome [26].
The objective of this study was to evaluate temporal changes in mammary gene network expression profiles by comparing two sheep breeds with different milking attitude. Gentile di Puglia (or Merino di Puglia, Pugliese Migliorata, Merino d'Italia, Merino Gentile) is a fine wooled breed from southern Italy. Development of this breed began in the 15th century, but the primary improvement began from the 18th century onwards. The breed was developed by crossing Spanish Merino with the local breeds. Today the selection objective of Gentile is focused onto meat production. Sarda is an Italian breed with high usefullness in milk production. It is widespread, mostly in Sardinia and in Central Italy, and representing 40% of the Italian ovine population.
In this study, we used a sheep-specific microarray chip technology covering most of the species’ transcriptome, representing the first annotated microarray developed for sheep with a covering of 50% of the genome [26]. The chip was generated from sheep ESTs deposited at NCBI and carries 21,743 non-redundant features in quadruplicate, 73.4% of which are fully annotated and corresponding to 10,190 genes. We analyzed the mammary transcriptome using biopsies from individuals of Gentile di Puglia and Sarda at two lactation stages to assess the differences between breeds, with the aim to identify genes controlling milk composition and their metabolic pathways.

2. Results and Discussion

We succefully hybridized eight microarray slides (four slides per lactation stage). Since every spot was replicated four times, for each lactation stage we performed 16 gene replicates (see Table 1).
In wet lab experiments 213 genes resulted differentially expressed between the two breeds at stages 01 (Table 2, Appendix) and 36 genes at stage 02 (Table 3, Appendix), with |FC| > 1.3, and p-value ≤ 0.05. The patterns of differentially expressed genes in Ovis aries were fully reproducible (see experimental section). At stage 01, 70 genes are upregulated in Gentile while 143 are upregulated in Sarda (Table 2, Appendix). At stage 02, only 8 genes are upregulated in Gentile while 28 are upregulated in Sarda, which is by far the most productive dairy breed (Table 3, Appendix).
We performed an analysis to show the most represented KEGG pathways among the differentially expressed genes (Figure 1, Figure 2), in order to identify molecular differences in milk synthesis between breeds and to identify genes controlling milk production and correlated metabolic pathways.
Nutrients 01 00235 g001 1024
Figure 1. KEGG pathway analysis of the differentially expressed genes at stage 01.
Figure 1. KEGG pathway analysis of the differentially expressed genes at stage 01.
Nutrients 01 00235 g001
At stage 01, among the 143 upregulated genes in Sarda, we could recognize caseins αS2, β and K. In addition to milk protein genes, we identified genes involved in processes linked to both lactation and mammary involution, as oxidative metabolism, apoptosis, cell cycle control, oncogenes, ubiquitination pathway and cell communication (focal adhesion adherens junction), endocrine system (insulin signaling pathway, adipocytokine signaling pathway) [27]. The KEGG pathways (like amino acid and carbohydrate metabolism, glycan biosynthesis, cell communication, cell growth and death and the immune system) were significantly (p < 0.05) enriched. The molecular events underlying mammary development during pregnancy, lactation and involution are incompletely understood. The processes of lactation include the development of mammary tissue, as well as the synthesis and secretion of milk. After parturition, the proliferation and differentiation of mammary secretory cells lead to an increase in milk secretion, whereas after lactation peak, milk production declines largely because of apoptotic mammary cell death, which exceeds cell proliferation. The development of mammary gland is spatially regulated by the communication of the mammary epithelium with the extracellular matrix (ECM) through a family of adhesion receptors called integrins. Integrins, in response to both hormones and growth factors, support cells in proliferation, accurate morphological organisation, as well as in milk secretion. Cell adhesion to the ECM plays a key role in alveolar survival, morphogenesis and function [28]. In this context, we could observe a significant difference, between the two breeds, in expression of genes involved in extracellular matrix formation and cell adhesion (TJP1 upregulated in Gentile, CDH5 and TNXB upregulated in Sarda). Remarkably, during stage 01, the expression of the oncogene VAV3 is higher in Gentile, while one of the initiators of apoptosis CFLAR (CASP8 and FADD-like apoptosis regulator) increases in Sarda. Apoptosis, in fact, occurs during involution of mammary gland in cattle [29,30], and an overexpression of many apoptosis-related genes during lactation was recently reported [31]. At stage 01 (early lactation), we found a differential expression of genes, like USP9X, involved in the ubiquitination pathway in Sarda. The protein ubiquitination pathway is the most significantly enriched pathway during both lactation and involution [32]. Another category of genes found differentially expressed between the two breeds encompasses genes involved in oxidoreductase activity, like cytochrome C oxidase, NADH dehydrogenase and ferritin. The activity of cytochrome C oxidase was found to increase from late pregnancy to the first days of lactation [33,34]. The overall expansion of oxidative metabolism is a response to the increased energy demands of the lactation period. At stage 01 we observed an upregulation of cytochrome C oxidase, NADH dehydrogenase and ferritin in Gentile di Puglia. This may reflect the different lactation persistence, which is lower in Gentile di Puglia (60–150 days) as compared to Sarda (210 days).
Nutrients 01 00235 g002 1024
Figure 2. KEGG pathway analysis of the differentially expressed genes at stage 02.
Figure 2. KEGG pathway analysis of the differentially expressed genes at stage 02.
Nutrients 01 00235 g002
At stage 02, some interesting genes as those encoding casein K, and proteins involved in oxidoreductase activity (like TGOLN2 and FTH1) and in ECM-interaction (like COL1A2), resulted overexpressed in Sarda. Finally, we can observe in Sarda an overexpression of genes implicated in lipolysis, like lipase (DAGLB) and phospholipase (PLD3). Several studies on different kind of cheeses have demonstrated that the fatty acid (FA) profile of raw milk influences cheese characteristics [22]. Lipolysis is particularly important in sheep cheeses due to the high fat content and lipase activity [35]. In this perspective, we like to stress that in the World the main output of sheep husbandry is cheese making.
Nutrients 01 00235 g003 1024
Figure 3. The two breeds analyzed: Gentile di Puglia (left) and Sarda (right).
Figure 3. The two breeds analyzed: Gentile di Puglia (left) and Sarda (right).
Nutrients 01 00235 g003

3. Experimental Section

3.1. Animals and Sampling

Whole mammary gland tissue samples were collected from four lactating individuals of two sheep (Ovis aries) breeds, Gentile di Puglia and Sarda (Figure 3). Lactating mammary tissue were taken at two lactation stages (first record, stage 01: 6 days after lambing; second record, stage 02: 44 days after lambing) in both breeds. Tissues from mammary gland were immersed in RNAlater (Sigma) and stored at −20 °C.

3.2. RNA Extraction

Tissues were subjected to RNA extraction with ice-cold TRIzol (Invitrogen) and using RNeasy Midi Kit columns (Qiagen). RNA integrity was assessed by electrophoretic analysis of 28S and 18S rRNA subunits. The purity of RNA and preliminary concentration were assessed with a spectrophotometer (GeneQuantpro). A260/A280 ratio was >1.9.

3.3. RNA Amplification and Labeling

RNA was quantified using a DTX fluorimeter (Beckman Coulter) using the Quant-iT kit (Invitrogen). Aliquots of 1 μg were amplified and Cy3/Cy5 labeled using the Kreatech Diagnostics kit.

3.4. Microarray Study Design and Hybridization

We designed an oligonucleotide chip from sheep ESTs deposited at NCBI. The oligonucleotide microarray platform is electrochemically synthesized and contains 21,743 non-redundant features in quadruplicate, 73.4% of which are fully annotated corresponding to 10,190 genes. A genome assembly for Ovis aries is not yet available, but considering the number of genes in the bovine genome (22,000), we estimate to have a 50% coverage of the sheep genome. We achieved very good technical outcomes, as reproducible patterns of differentially expressed genes (in each slide, replicates show a coefficient of variation <0.25 for differentially expressed genes with P < 0.01) [26]. Oligos were generated in situ on the chip using the Combimatrix (Seattle, WA, USA) equipment. Platform and microarray data have been deposited in the NCBI GEO database (Platform Accession GPL9461; Series Accession GSE18619).
The labeled aRNA was fragmented into 35-200-base fragments and then hybridized onto the slide according to Combimatrix’s instructions. Hybridization was performed overnight at 50 °C. After hybridization, arrays were washed and scanned with a ScanArray Lite (Perkin Elmer) laser scanner. Microarray Imager 5.9.3 software was used to extract feature data from microarray fluorescence images.
The microarray study was designed as described in Table 1. Separate microarrays were used for individual samples. At each stage, RNA of from one lactating individual of each breed was labeled in dual color using Cy3 and Cy5 fluorochromes. RNA aliquots from the same stage of distinct breeds, labeled with different fluorocromes, were hybridized on the same microarray slide. We performed two technical microarray replicates per stage (four slides). The entire microarray experiment was repeated starting from a new RNA extraction of the same tissue sample (see Table 1), for a total of eight slides.
Table
Table 1. Microarray experimental design.
Table 1. Microarray experimental design.
Nutrients 01 00235 i001

3.5. Microarray Data Analysis

Single microarray step

Saturated (foreground median intensity, FMI, over the limits) and bad spots were flagged using the software Microarray Imager 5.9.3 (CombiMatrix). For each channel (Cy5, red and Cy3, green) the mean of empty spots FMI (E) was calculated. Non-empty spots with FMI to E ratio below 1.5 at least in one channel were filtered, together with flagged spots. We calculated the M and the A value (M = log2(R/G), A = (log2(R*G))/2, R = red FMI, G = green FMI) for each spot in order to obtain a measure of the differential expression for the two conditions analysed (M) and the log-mean of the spot intensity (A).

Paired dye-swap microarrays step and normalization

Systematic bias in the data was removed by applying the dye-swap normalization, that makes use of the reverse labelling in the two microarray replicates in order to remove the intrinsic difference of the two fluorochromes output [36]. In particular we paired each microarray with its dye-swap by calculating a new M value, MD = (M1 − M2)/2, where M1 is the M value for an experiment and M2 is the M value calculated for the correspondent experiment with inverse fluorochromes, and a new A value, AD = A1 + A2)/2, where A1 is the A value for an experiment and and A2 is the A value calculated for the correspondent experiment with inverse fluorochromes. In order to remove intensity based bias we also applied the lowess normalization [37], obtaining a new M value (ML) from MD.

Significance analysis

For each lactation stage, we performed one sample t-test to establish, for each gene, if the mean of the ML values (uML) was significantly different from 0, and corrected the p-value for multiple comparisons with the Benjamini and Hochberg False Discovery Rate [38]. Finally, only genes with a satisfactory effect (absolute value of the fold change |FC*| > 1.3, FC = 2^ uML, if FC<0, FC* = −1/FC, if FC > 0, FC* = FC) and a significant p-values were considered. The statistical significance of the enrichment for the KEGG pathways of interest was computed using the hypergeometric test [39].

4. Conclusions

Sheep farming is very important for cheese production and the fatty acid and protein composition of raw milk is crucial in the cheese making process. The fatty acid profile of raw milk has been demonstrated to affect cheese characteristics and differentiate new types of cheese [21,40]. The genetic differences between breeds on milk quality are likely to affect also cheese quality and could be a marker to carry out genetic improvement plans of local and endangered sheep breeds. However, the number of studies on gene expression analysis between breeds aimed at understanding how genetic variations affect milk quality is quite limited in sheep as compared to cow. The main difficulty has been, to date, the absence of devices like microarrays. In this paper we have proven that a homologous chip, generated from sheep ESTs, is a valuable tool which can be employed in gene expression analysis. Furthermore, this approach can be easily extended to other species of which genetic sequences are present in public databases [26].

Acknowledgements

This work has been supported by the GENZOOT project, funded by the Italian Ministry of Agricultural and Forestry Policies, (MiPAF). We wish to thank Dr. Bianca Moioli, GENZOOT project coordinator; the CRA in Monterotondo (Roma) for samples supplying and stimulating discussion; Paolo Ciorba and Gabriella Porcai for technical assistance.

Appendix

Table
Table 2. Differentially expressed genes, stage 01.
Table 2. Differentially expressed genes, stage 01.
P value FC EST RefSeq Gene UPDOWN
0.0089536-4.677EE874449.1----GentileSarda
0.0399319-2.705CD287057.1----GentileSarda
0.00629989-2.578EE788795.1----GentileSarda
0.0148422-2.489EE873971.1NM_001083800IGLL1GentileSarda
0.0107907-2.339DY479414.1----GentileSarda
0.00171376-2.275EE874479.1NM_001025317RPS8GentileSarda
0.00263317-2.249EE814758.1NM_001015592PFN1GentileSarda
0.00360217-2.164EE781281.1----GentileSarda
0.0110063-2.07EE866465.1NM_001033614RPS18GentileSarda
0.0164299-2.067EE873426.1NM_001105455RPL39GentileSarda
0.0354098-2.015DY497824.1--FTH1GentileSarda
0.0306295-1.918EE788254.1----GentileSarda
0.0184745-1.753EE856480.1--FTH1GentileSarda
0.0185062-1.697EE824030.1XM_615127ASPHD2GentileSarda
0.00474836-1.682EE874471.1NM_001034438RPS20GentileSarda
0.0434131-1.632EE865486.1NM_001045929MPV17GentileSarda
0.0250763-1.588EE831966.1----GentileSarda
0.019839-1.582EE866322.1XM_615898VAV3GentileSarda
0.000864936-1.57EE848574.1----GentileSarda
0.0424348-1.568EE774810.1XM_001373571LOC100021448GentileSarda
0.0043326-1.546EE866312.1XM_583460RPUSD2GentileSarda
0.0440448-1.522EE862726.1XM_001787912PLEKHA2GentileSarda
0.0263277-1.52DY485409.1----GentileSarda
0.04731-1.514EE862521.1NM_001076831COL3A1GentileSarda
0.0135288-1.509EE750744.1XM_865008ERBB2IPGentileSarda
0.00245474-1.492EE831293.1----GentileSarda
0.0390691-1.458EE831369.1NM_001046226FEM1AGentileSarda
0.0449235-1.448CO202828.1----GentileSarda
0.0107052-1.445DY484904.1----GentileSarda
0.0275188-1.437EE871094.1----GentileSarda
0.0279427-1.436EE864971.1NM_174590PTPN13GentileSarda
0.0267974-1.432EE792695.1NM_001101894FBXO11GentileSarda
0.000833469-1.419EE816833.1NM_001075714LOC515452GentileSarda
0.000176532-1.418EE874214.1NM_003032ST6GAL1GentileSarda
0.0471172-1.404EE777222.1NM_001001441TNNT3GentileSarda
0.0401322-1.401EE770300.1NM_001034333SAT1GentileSarda
0.00551675-1.4DY513982.1NM_001076011ITPKCGentileSarda
0.0103614-1.398EE802605.1NM_001078161LOC777786GentileSarda
0.0244517-1.389EE849733.1----GentileSarda
0.00575702-1.386EE830948.1XM_001790594LOC100139162GentileSarda
0.0251925-1.377EE873625.1NM_001079783UBE2E3GentileSarda
0.00381358-1.374EE816897.1----GentileSarda
0.00653387-1.373EE837747.1NM_001075367ALDH2GentileSarda
0.0100565-1.373EE846386.1NM_001101154ALAS1GentileSarda
0.0496638-1.37DY500795.1NM_001035376MEOX1GentileSarda
0.0090203-1.367EE825590.1XM_001109283LOC712430GentileSarda
0.015555-1.361DY512463.1NM_001101929HK3GentileSarda
0.0104062-1.358EE807632.1XM_612129ATP2A2GentileSarda
0.00419809-1.357EE798194.1XM_001250724PDCLGentileSarda
0.0129818-1.347EE755345.1XM_865894MORC3GentileSarda
0.0138409-1.344EE795474.1NM_001105041RUSC1GentileSarda
0.0268807-1.342EE864960.1NM_001045971SPINT2GentileSarda
0.0404747-1.339EE812525.1NM_001102149MGC159954GentileSarda
0.0464112-1.339EE869934.1NM_001100314PIK3R4GentileSarda
0.027626-1.338EE792100.1----GentileSarda
0.00996729-1.337EE864412.1NM_001113764TYK2GentileSarda
0.0121423-1.331EE832631.1NM_001034743RAB5CGentileSarda
0.0318396-1.329EE831575.1NM_001046043ANGPTL4GentileSarda
0.035704-1.325EE767566.1XM_001089198beta- galactosidaseGentileSarda
0.000597784-1.32EE857107.1----GentileSarda
0.0323903-1.319EE872862.1NM_001075700RAP2CGentileSarda
0.0421994-1.319EE841857.1NM_001042682REREGentileSarda
0.0208483-1.309EE810668.1NM_001099167TMEM149GentileSarda
0.0279321-1.309EE798509.1----GentileSarda
0.0326586-1.307EE842831.1----GentileSarda
0.00133131-1.304DY521037.1XM_001788744LOC513508GentileSarda
0.0189923-1.304EE867731.1----GentileSarda
0.0237616-1.302EE816414.1XM_614378SPRYD3GentileSarda
0.043391-1.302EE858582.1NM_001077971GRSF1GentileSarda
0.00992968-1.301EE862241.1NM_175610TJP1GentileSarda
0.04873211.301EE830752.1XM_001789168LOC100138505SardaGentile
0.01538921.302EE770765.1XM_001787762HEXBSardaGentile
0.02178681.303EE851241.1NM_138782FCHO2SardaGentile
0.01676421.303EE801658.1NM_001015555AUP1SardaGentile
0.04304051.304EE831852.1NM_001102287FANCGSardaGentile
0.03375511.305EE819075.1NM_001075848RASGRP3SardaGentile
0.01529041.307EE759878.1----SardaGentile
0.007926881.307EE869797.1XM_001789365LOC783484SardaGentile
0.03053441.308EE839816.1NM_001038561RPUSD3SardaGentile
0.007278231.309DY520684.1XM_001720318LOC100129623SardaGentile
0.01030581.31EE845343.1XM_001787123CALCASardaGentile
0.0002128561.311EE867028.1NM_001101171ABI2SardaGentile
0.01981071.312EE782033.1XM_001256327LOC789629SardaGentile
0.04052721.313EE783717.1----SardaGentile
0.002399571.313CF117405.1XM_599530LOC521270SardaGentile
0.01335751.315EE849952.1----SardaGentile
0.01745281.315EE770796.1----SardaGentile
0.01963971.315EE871964.1----SardaGentile
0.02111851.315EE873622.1NM_174718PNNSardaGentile
0.01749241.317EE866027.1NM_053043RBM33SardaGentile
0.002246321.319EE833196.1----SardaGentile
0.00772931.319EE746291.1XR_042970LOC532848SardaGentile
0.04531951.321EE830242.1----SardaGentile
0.001522611.321EE848826.1XM_001788615ANAPC11SardaGentile
0.03713481.321EE861720.1----SardaGentile
0.01247631.322EE840376.1XR_027670LOC539015SardaGentile
0.01106951.323EE836611.1NM_001101080ADAMTS1SardaGentile
0.008695261.324EE854385.1----SardaGentile
0.04708651.325EE824979.1NM_001286CLCN6SardaGentile
0.0005916041.326EE834601.1----SardaGentile
0.01648171.326EE872298.1XM_875686CLDND1SardaGentile
0.01294291.327EE865545.1NM_001102062WDR75SardaGentile
0.009360871.33EE856712.1XM_001788055LOC100138621SardaGentile
0.04221281.332EE836140.1XM_001789542LOC100139498SardaGentile
0.000991391.334EE865157.1NM_018509LRRC59SardaGentile
0.04199621.335EE773093.1----SardaGentile
0.00413491.335EE822639.1NM_001034368ABHD4SardaGentile
0.01911051.337EE842849.1XR_042867ZNF134SardaGentile
0.01513361.338EE825761.1XM_001256069LOC789273SardaGentile
0.02875271.338EE827215.1NM_001002892ST3GAL2SardaGentile
0.01068971.339EE861357.1XM_001255930LOC789066SardaGentile
0.005474981.339EE796870.1XM_614279LGI2SardaGentile
0.03251181.341EE784162.1----SardaGentile
0.01302261.342EE795222.1--IFNAR1ESardaGentile
0.03107871.342EE793726.1----SardaGentile
0.02255651.343EE780347.1----SardaGentile
0.00411031.343EE826336.1NM_001076049EFEMP2SardaGentile
0.008472581.345EE847787.1NM_001075670SMAP2SardaGentile
0.01448381.348EE848020.1XM_593447CCDC61SardaGentile
0.02519651.348DY491137.1----SardaGentile
0.014721.35EE747969.1----SardaGentile
0.02305451.355EE872727.1NM_001078102DNAJC24SardaGentile
0.03403451.357EE797130.1NM_020772NUFIP2SardaGentile
0.04000431.357EE782700.1----SardaGentile
0.01994721.358EE852157.1NM_172127CAMK2DSardaGentile
0.02841781.36EE837214.1NM_001098003TMCO3SardaGentile
0.01522161.36EE844290.1NM_001102498NKAPLSardaGentile
0.0006640051.36EE816910.1----SardaGentile
0.01296561.365EE813986.1XM_584123LOC538993SardaGentile
0.002395961.368EE833638.1----SardaGentile
0.004730051.371EE776285.1XM_600379AHDC1SardaGentile
0.03133931.372EE837130.1NM_001025345MCM7SardaGentile
0.002578391.373DY522523.1NM_001102074QSOX1SardaGentile
0.03088611.373EE868991.1NM_001034633SLC3A1SardaGentile
0.04076811.375EE746824.1XM_865771RUNX1SardaGentile
0.03932771.379EE794485.1----SardaGentile
0.03078121.38EE815257.1NM_001046011CD37SardaGentile
0.01845871.384EE816699.1NM_001083429PANK4SardaGentile
0.0006063511.384EE830551.1NM_001037607ARFRP1SardaGentile
0.04245281.385EE753287.1----SardaGentile
0.02131581.385EE826333.1XM_590179DENND2ASardaGentile
0.02645851.386EE851940.1NM_012319SLC39A6SardaGentile
0.006471711.387EE790636.1NM_001045969ALDH7A1SardaGentile
0.02268061.388CN824197.1NM_053064GNG2SardaGentile
0.03876531.388EE813255.1NM_030935TSC22D4SardaGentile
0.001749151.393EE791965.1NM_001046497MAT1ASardaGentile
0.02648751.397DY490978.1--NDUFS2SardaGentile
0.002957381.398EE829627.1NM_001077854MS4A1SardaGentile
0.02240811.399EE854207.1NM_001077104KRTAP3-1SardaGentile
0.04702491.402EE858765.1NM_001083703RNF185SardaGentile
0.003718681.402EE832497.1NM_001078041PLD3SardaGentile
0.02165811.402EE866310.1NM_001099072UNC45ASardaGentile
0.0103591.403EE808144.1XM_001250150TAF2SardaGentile
0.03178681.407EE833224.1--CDH5SardaGentile
0.005063131.408EE794195.1----SardaGentile
0.03237861.408DY478310.1NM_000088COL1A1SardaGentile
0.006253661.41EE825745.1NM_001102546PDE4BSardaGentile
0.03966241.412EE764887.1----SardaGentile
0.02406261.413EE818092.1NM_001102238METTL3SardaGentile
0.0004523411.416EE843072.1NM_178140PDZD2SardaGentile
0.007725431.418EE856741.1NM_001075176RPA1SardaGentile
0.00559551.421EE823634.1----SardaGentile
0.01402541.422EE864563.1NM_001034339IL11RASardaGentile
4.02e-051.424DY480261.1----SardaGentile
0.04449031.428EE816651.1NM_001083414WWP2SardaGentile
0.0008771211.432EE870673.1----SardaGentile
0.01542731.442EE766357.1XM_585246KIF11SardaGentile
0.007507411.446EE825953.1XM_865238EXOC2SardaGentile
0.002222571.446EE788062.1NM_019863F8SardaGentile
0.007433561.446EE747821.1----SardaGentile
0.03727421.447EE871308.1----SardaGentile
0.0154381.448EE851254.1----SardaGentile
0.00811821.451DY520937.1NM_001102035CERCAMSardaGentile
0.003718391.458EE873738.1NM_001035283TALDO1SardaGentile
0.006129511.458EE765691.1--GALNSardaGentile
0.02112051.46EE780060.1NM_001012281CFLAR SardaGentile
0.01492151.461EE821548.1NM_001046390TESSardaGentile
0.03376441.464BG874259.1----SardaGentile
0.03475331.466DY514734.1XM_589271FBLN2SardaGentile
0.006527371.469EE783524.1----SardaGentile
0.01216991.469EE828474.1XM_870386SH3BP5SardaGentile
0.001710571.477DY504539.1----SardaGentile
0.004928991.478EE856049.1XM_865072PIGTSardaGentile
0.001546591.479EE871601.1----SardaGentile
0.01381371.48EE831423.1XM_001788882LOC789539SardaGentile
0.009525631.481EE829269.1NM_001014865DRG2SardaGentile
0.008085431.489EE818202.1XM_602855GALNT7SardaGentile
0.02124841.491DY485302.1NM_001024824RIPK5SardaGentile
0.02048471.498EE871638.1XM_870378CAPN7SardaGentile
0.01850141.511DY480855.1XM_001787789LOC534471SardaGentile
0.01611771.514EE748326.1NM_001046194CDCA2SardaGentile
0.008259151.514DY496078.1XM_001789157LOC508459SardaGentile
0.01938491.517EE841331.1NM_001039591USP9XSardaGentile
0.005440021.54EE767854.1NM_001076330PRSS16SardaGentile
0.006943911.541EE820835.1NM_001045866BRD2SardaGentile
0.0481831.544EE760222.1XM_001252101LOC784704SardaGentile
0.005852421.551DY491388.1XM_596546SH3PXD2BSardaGentile
0.02422481.554EE749164.1NM_174676RASA3SardaGentile
0.01797861.561EE750374.1----SardaGentile
0.04528291.564EE857270.1XM_596854PLEKSardaGentile
0.0362561.566EE798015.1----SardaGentile
0.000573571.579EE856496.1NM_001046252LITAFSardaGentile
0.01232831.599EE825369.1NM_001081602STK38SardaGentile
0.008330881.603EE812445.1XM_001254977LRRCC1SardaGentile
0.04288151.612EE826810.1NM_001035012RIPK1SardaGentile
0.001790871.642DY500867.1XM_001254158ZNF318SardaGentile
0.006312941.682EE849102.1----SardaGentile
7.62e-051.877EE799074.1NM_174703TNXBSardaGentile
0.01443251.931EE791703.1----SardaGentile
0.01245892.018EE874444.1--CSN2SardaGentile
0.04684522.039EE816347.1NM_001075485LYSMD2SardaGentile
0.02284272.078EE874341.1--LOC443383SardaGentile
0.01760612.197EE874443.1--CSN3SardaGentile
Table
Table 3. Differentially expressed genes, stage 02.
Table 3. Differentially expressed genes, stage 02.
pvalueFC EST RefSeq Gene UPDOWN
0.0262055-1.467EE855884.1NM_001038689C5H12orf45 GentileSarda
0.0440456-1.391EE804883.1NM_172216CAMKK2 GentileSarda
0.0295858-1.37EE752417.1NM_174438PROS1 GentileSarda
0.00912544-1.344EE805588.1NM_001105615CENPP GentileSarda
0.0356366-1.341CF117455.1XM_001254445FBXO48 GentileSarda
0.0153359-1.332EE831205.1NM_001046346WDR1 GentileSarda
0.0167085-1.304EE872615.1NM_001099859EIF4G2 GentileSarda
0.00236899-1.302EE803274.1NM_001102100LENG8 GentileSarda
0.02203971.303EE794531.1NM_001083487DAGLB SardaGentile
0.01175021.305CD288635.1NM_001083793SMAGP SardaGentile
0.009153091.308EE811805.1---- SardaGentile
0.04809771.314DY516851.1XR_028016LOC539596 SardaGentile
0.005974851.316EE869114.1---- SardaGentile
0.01590311.316EE749850.1XM_001249987EEF1A1 SardaGentile
0.04832091.32EE810406.1NM_001045938AP1B1 SardaGentile
0.01644871.322EE840745.1NM_194442LBR SardaGentile
0.01940881.323EE752798.1NM_001083462SEC22A SardaGentile
0.02054221.324DY514991.1XM_001251051LOC782414 SardaGentile
0.008254891.325EE872139.1XR_027898LOC784355 SardaGentile
0.0008613531.327EE858050.1XM_580387PRKD3 SardaGentile
0.03864971.33EE788388.1---- SardaGentile
0.0004856621.335EE844249.1NM_001075142IL4R SardaGentile
0.008427721.336EE816344.1---- SardaGentile
0.02906681.338EE874443.1--CSN3 SardaGentile
0.0328051.342EE834071.1XM_583748TBX21 SardaGentile
0.003633291.359EE860008.1XM_613386ZNF532 SardaGentile
0.0004517451.361DY520339.1NM_174520COL1A2 SardaGentile
0.0282631.362DY480878.1---- SardaGentile
0.0462961.392EE820416.1---- SardaGentile
0.005298471.405DY520664.1NM_001075474MRPS16 SardaGentile
0.01957541.417CF117857.1NM_001098860RAMP2 SardaGentile
0.000721331.456CD287057.1---- SardaGentile
0.001984831.465DY479000.1---- SardaGentile
0.0222871.47EE801908.1NM_001046214ACBD4 SardaGentile
0.01055841.484EE815399.1XM_589331TGOLN2 SardaGentile
0.006576551.518DY497824.1--FTH1 SardaGentile

References and Notes

  1. Allen, J.C.; Keller, R.P.; Archer, P.; Neville, M.C. Studies in human lactation: Milk composition and daily secretion rates of macronutrients in the first year of lactation. Am. J. Clin. Nutr. 1991, 54, 69–80. [Google Scholar]
  2. Neville, M.C. Anatomy and physiology of lactation. Pediatr. Clin. N. Am. 2001, 48, 13–34. [Google Scholar]
  3. Capuco, A.V.; Akers, R.M. The origin and evolution of lactation. J. Biol. 2009, 8 art. no. 37. [Google Scholar]
  4. Chadwick, R. Nutrigenomics, individualism and public health. Proc. Nutr. Soc. 2004, 63, 161–166. [Google Scholar] [Green Version]
  5. Van Knegsel, A.T.M.; Van Den Brand, H.; Dijkstra, J.; Van Straalen, W.M.; Heetkamp, M.J.W.; Tamminga, S.; Kemp, B. Dietary energy source in dairy cows in early lactation: Energy partitioning and milk composition. J. Dairy Sci. 2007, 90, 1467–1476. [Google Scholar]
  6. Van Knegsel, A.T.M.; Van Den Brand, H.; Graat, E.A.M.; Dijkstra, J.; Jorritsma, R.; Decuypere, E.; Tamminga, S.; Kemp, B. Dietary energy source in dairy cows in early lactation: Metabolites and metabolic hormones. J. Dairy Sci. 2007, 90, 1477–1485. [Google Scholar]
  7. Palmquist, D.L.; Beaulieu, A.D.; Barbano, D.M. Feed and animal factors influencing milk fat composition. J. Dairy Sci. 1993, 76, 1753–1771. [Google Scholar]
  8. Kim, S.W.; Wu, G. Regulatory role for amino acids in mammary gland growth and milk synthesis. Amino Acids 2009, 37, 89–95. [Google Scholar]
  9. Andrechek, E.R.; Mori, S.; Rempel, R.E.; Chang, J.T.; Nevins, J.R. Patterns of cell signaling pathway activation that characterize mammary development. Development 2008, 135, 2403–2413. [Google Scholar]
  10. Hartmann, P.E. Changes in the composition and yield of the mammary secretion of cows during the initiation of lactation. J. Endocrinol. 1973, 59, 231–247. [Google Scholar]
  11. Hartmann, P.E.; Trevethan, P.; Shelton, J.N. Progesterone and oestrogen and the initiation of lactation in ewes. J. Endocrinol. 1973, 59, 249–259. [Google Scholar]
  12. Neville, M.C.; McFadden, T.B.; Forsyth, I. Hormonal regulation of mammary differentiation and milk secretion. J. Mammary Gland Biol. 2002, 7, 49–66. [Google Scholar]
  13. McManaman, J.L.; Neville, M.C. Mammary physiology and milk secretion. Ad. Drug Deliver Rev. 2003, 55, 629–641. [Google Scholar]
  14. Capuco, A.V.; Wood, D.L.; Baldwin, R.; Mcleod, K.; Paape, M.J. Mammary cell number, proliferation, and apoptosis during a bovine lactation: Relation to milk production and effect of bST. J. Dairy Sci. 2001, 84, 2177–2187. [Google Scholar] [CrossRef] [PubMed]
  15. Sorensen, M.T.; Nørgaard, J.V.; Theil, P.K.; Vestergaard, M.; Sejrsen, K. Cell turnover and activity in mammary tissue during lactation and the dry period in dairy cows. J. Dairy Sci. 2006, 89, 4632–4639. [Google Scholar]
  16. Finucane, K.A.; McFadden, T.B.; Bond, J.P.; Kennelly, J.J.; Zhao, F.Q. Onset of lactation in the bovine mammary gland: Gene expression profiling indicates a strong inhibition of gene expression in cell proliferation. Funct. Integr. Genomic. 2008, 8, 251–264. [Google Scholar]
  17. Baxter, F.O.; Neoh, K.; Tevendale, M.C. The beginning of the end: Death signaling in early involution. J. Mammary Gland Biol. 2007, 12, 3–13. [Google Scholar]
  18. Stein, T.; Salomonis, N.; Gusterson, B.A. Mammary gland involution as a multi-step process. J. Mammary Gland Biol. 2007, 12, 25–35. [Google Scholar] [CrossRef]
  19. Capuco, A.V.; Akers, R.M. Mammary Involution in Dairy Animals. J. Mammary Gland Biol 1999, 4, 37–144. [Google Scholar]
  20. Haenlein, G.F.W. Past, present, and future perspectives of small ruminant dairy research. J. Dairy Sci. 2001, 84, 2097–2115. [Google Scholar] [CrossRef] [PubMed]
  21. Signorelli, F.; Contarini, G.; Annicchiarico, G.; Napolitano, F.; Orrù, L.; Catillo, G.; Haenlein, G.F.W.; Moioli, B. Breed differences in sheep milk fatty acid profiles: Opportunities for sustainable use of animal genetic resources. Small Ruminant Res. 2008, 78, 24–31. [Google Scholar]
  22. Moioli, B.; D'Andrea, M.; Pilla, F. Candidate genes affecting sheep and goat milk quality. Small Ruminant Res. 2007, 68, 179–192. [Google Scholar]
  23. Kaput, J.; Rodriguez, R.L. Nutritional genomics: The next frontier in the postgenomic era. Phys. Genomic. 2004, 16, 166–177. [Google Scholar]
  24. Kato, H.; Saito, K.; Kimura, T. A perspective on DNA microarray technology in food and nutritional science. Curr. Opin. Clin. Nutr. 2005, 8, 516–522. [Google Scholar]
  25. Zduńczyk, Z.; Pareek, Ch.S. Application of nutrigenomics tools in animal feeding and nutritional research. J. Anim. Feed Sci. 2009, 18, 13–16. [Google Scholar]
  26. Pariset, L.; Chillemi, G.; Bongiorni, S.; Spica, V.R.; Valentini, A. Microarrays and high-throughput transcriptomic analysis in species with incomplete availability of genomic sequences. New Biotech. 2009, 25, 272–279. [Google Scholar]
  27. Bongiorni, S.; Chillemi, G.; Prosperini, G.; Bueno, S.; Signorelli, F.; Moioli, B.; Pariset, L. Transcriptomic analysis of two sheep breeds during lactation, using a new custom microarray platform. Ital. J. Anim. Sci. 2009, 8, 33–35. [Google Scholar]
  28. Katz, E.; Streuli, C.H. The extracellular matrix as an adhesion checkpoint for mammary epithelial function. Inter. J. Biochem. Cell B. 2007, 39, 715–726. [Google Scholar]
  29. Wilde, C.J.; Quarrie, L.H.; Tonner, E.; Flint, D.J.; Peaker, M. Mammary apoptosis. Livest. Prod. Sci. 1997, 50, 29–37. [Google Scholar]
  30. Wilde, C.J.; Addey, C.V.P.; Li, P.; Fernig, D.G. Programmed cell death in bovine mammary tissue during lactation and involution. Exp. Physiol. 1997, 82, 943–953. [Google Scholar]
  31. Suchyta, S.P.; Sipkovsky, S.; Halgren, R.G.; Kruska, R.; Elftman, M.; Weber-Nielsen, M.; Vandehaar, M.J.B; Coussens, P.M. Bovine mammary gene expression profiling using a cDNA microarray enhanced for mammary-specific transcripts. Phys. Genomic 2004, 16, 8–18. [Google Scholar]
  32. Lemay, D.G.; Neville, M.C.; Rudolph, M.C.; Pollard, K.S.; German, J.B. Gene regulatory networks in lactation: Identification of global principles using bioinformatics. BMC Syst. Biol. 2007, 1 art. no. 56. [Google Scholar]
  33. Rosano, T.G.; Jones, D.H. Developmental changes in mitochondria during the transition into lactation in the mouse mammary gland. I. Behavior on isopycnic gradient centrifugation. J. Cell Biol. 1976, 69, 573–580. [Google Scholar]
  34. Rosano, T.G.; Lee, S.K.; Jones, D.H. Developmental changes in mitochondria during the transition into lactation in the mouse mammary gland. II. Membrane marker enzymes and membrane ultrastructure. J. Cell Biol. 1976, 69, 581–588. [Google Scholar]
  35. Fontecha, J.; Mayo, I.; Toledano, G.; Juárez, M. Use of changes in triacylglycerols during ripening of cheeses with high lipolysis levels for detection of milk fat authenticity. Int. Dairy J. 2006, 16, 1498–1504. [Google Scholar]
  36. Fang, Y.; Brass, A.; Hoyle, D.C.; Hayes, A.; Bashein, A.; Oliver, S.G.; Waddingtong, D.; Rattray, M. A model-based analysis of microarray experimental error and normalisation. Nucleic Acids Res. 2003, 31, e96. [Google Scholar]
  37. Yang, Y.H.; Dudoit, S.; Luu, P.; Speed, T.P. Normalization for cDNAmicroarray data. In Microarrays: optical technologies and informatics; Bittner, M., Chen, Y., Dorsel, A., Dougherty, E.R., Eds.; SPIE: San Jose, CA, USA, 2001; pp. 141–152. [Google Scholar]
  38. Benjamini, Y.; Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. 1995, 57, 289–300. [Google Scholar]
  39. Khatri, P.; Drǎghici, S. Ontological analysis of gene expression data: Current tools, limitations, and open problems. Bioinformatics 2005, 21, 3587–3595. [Google Scholar] [CrossRef] [PubMed]
  40. Ha, J.K.; Lindsay, R.C. Release of volatile branched-chain and other fatty acids from ruminant milk fats by various lipases. J. Dairy Sci. 1993, 76, 677–690. [Google Scholar]

Share and Cite

MDPI and ACS Style

Bongiorni, S.; Chillemi, G.; Prosperini, G.; Bueno, S.; Valentini, A.; Pariset, L. A Tool for Sheep Product Quality: Custom Microarrays from Public Databases. Nutrients 2009, 1, 235-250. https://doi.org/10.3390/nu1020235

AMA Style

Bongiorni S, Chillemi G, Prosperini G, Bueno S, Valentini A, Pariset L. A Tool for Sheep Product Quality: Custom Microarrays from Public Databases. Nutrients. 2009; 1(2):235-250. https://doi.org/10.3390/nu1020235

Chicago/Turabian Style

Bongiorni, Silvia, Giovanni Chillemi, Gianluca Prosperini, Susana Bueno, Alessio Valentini, and Lorraine Pariset. 2009. "A Tool for Sheep Product Quality: Custom Microarrays from Public Databases" Nutrients 1, no. 2: 235-250. https://doi.org/10.3390/nu1020235

APA Style

Bongiorni, S., Chillemi, G., Prosperini, G., Bueno, S., Valentini, A., & Pariset, L. (2009). A Tool for Sheep Product Quality: Custom Microarrays from Public Databases. Nutrients, 1(2), 235-250. https://doi.org/10.3390/nu1020235

Article Metrics

Back to TopTop