Gene expression profiling in porcine mammary gland during lactation and identification of breed- and developmental-stage-specific genes

Characterization and comparative analysis of transcriptional profiles of porcine colostrum and mature milk at different parities

BACKGROUND

Porcine milk is a complex fluid, containing a myriad of immunological, biochemical, and cellular components, made to satisfy the nutritional requirements of the neonate. Whole milk contains many different cell types, including mammary epithelial cells, neutrophils, macrophages, and lymphocytes, as well nanoparticles, such as milk exosomes. To-date, only a limited number of livestock transcriptomic studies have reported sequencing of milk. Moreover, those studies focused only on sequencing somatic cells as a proxy for the mammary gland with the goal of investigating differences in the lactation process. Recent studies have indicated that RNA originating from multiple cell types present in milk can withstand harsh environments, such as the digestive system, and transmit regulatory molecules from maternal to neonate. Transcriptomic profiling of porcine whole milk, which is reflective of the combined cell populations, could help elucidate these mechanisms. To this end, total RNA from colostrum and mature milk samples were sequenced from 65 sows at differing parities. A stringent bioinformatic pipeline was used to identify and characterize 70,841 transcripts.

RESULTS

The 70,841 identified transcripts included 42,733 previously annotated transcripts and 28,108 novel transcripts. Differential gene expression analysis was conducted using a generalized linear model coupled with the Lancaster method for P-value aggregation across transcripts. In total, 1667 differentially expressed genes (DEG) were identified for the milk type main effect, and 33 DEG were identified for the milk type x parity interaction. Several gene ontology (GO) terms related to immune response were significant for the milk type main effect, supporting the well-known fact that immunoglobulins and immune cells are transferred to the neonate via colostrum.

CONCLUSIONS

This is the first study to perform global transcriptome analysis from whole milk samples in sows from different parities. Our results provide important information and insight into synthesis of milk proteins and innate immunity and potential targets for future improvement of swine lactation and piglet development.

Hyperprolactinemia using domperidone in prepubertal gilts: Effects on hormonal status, mammary development and mammary and pituitary gene expression

Objectives of this experiment were to determine if the domperidone protocol previously used for gestating gilts can also lead to hyperprolactinemia in growing gilts, and to assess the effects of such a protocol on hormonal status, mammary development and gene expression in mammary and pituitary tissue of gilts at puberty. The impact on future lactation performance was also determined. At 75 ± 3 kg body weight (BW), gilts were divided between: 1) controls (CTL), receiving daily intramuscular (IM) injections of canola oil (1.1 mL) for 29 d (n = 41), and 2) treated (DOMP), receiving daily IM injections with 0.5 mg/kg BW of the dopamine receptor antagonist domperidone for 29 d (n = 40). In addition to that daily injection, treated gilts also received twice daily IM injections with 0.5 mg/kg BW of domperidone over the first 3 d of treatment. Fifteen gilts per treatment were sacrificed at 210 ± 5 d of age to collect mammary glands (for compositional analysis and gene expression) and the anterior pituitary (for gene expression). Remaining gilts were bred and allowed to farrow. Blood was sampled at the onset of treatment and on days 14 and 30. Gilts that farrowed were also blood sampled on days 3 and 20 of lactation. Blood was assayed for prolactin (PRL), leptin, insulin-like growth factor 1 (IGF-1), urea, free fatty acids and glucose. Concentrations of PRL increased after 14 d and 30 d of treatment (P < 0.01) and were lesser on day 3 of lactation in DOMP than CTL gilts (P < 0.01). At puberty, there were tendencies (P < 0.10) for total parenchymal protein and DNA to be greater in DOMP than CTL gilts. Treatment did not affect mRNA abundance of PRL or the long form of the PRL receptor genes in the pituitary gland at puberty but expression level of the dopamine receptor D2 and PRL genes was much lower in pubertal than late-pregnant gilts (P < 0.001). Furthermore, many genes related with PRL had a much greater expression level in late pregnancy than at puberty. On day 20 of lactation, CTL sows had greater concentrations of urea than DOMP sows (P < 0.01). The growth rate of litters was not affected by treatment nor was milk composition (P > 0.10). Even though PRL concentrations were increased with treatment, the absence of effect on mammary development was either due to timing relative to developmental stage, whereby treatment was initiated when gilts were too young, or was because all PRL receptors may have been saturated thereby preventing biological action of additional PRL.

Recent progress of porcine milk components and mammary gland function

As the only nutritional source for newborn piglets, porcine colostrum and milk contain critical nutritional and immunological components including carbohydrates, lipids, and proteins (immunoglobulins). However, porcine milk composition is more complex than these three components. Recently, scientists identified additional and novel components of sow colostrum and milk, including exosomes, oligosaccharides, and bacteria, which possibly act as biological signals and modulate the intestinal environment and immune status in piglets and later in life. Evaluation of these nutritional and non-nutritional components in porcine milk will help better understand the nutritional and biological function of porcine colostrum and milk. Furthermore, some important functions of the porcine mammary gland have been reported in recent published literature. These preliminary studies hypothesized how glucose, amino acids, and fatty acids are transported from maternal blood to the porcine mammary gland for milk synthesis. Therefore, we summarized recent reports on sow milk composition and porcine mammary gland function in this review, with particular emphasis on macronutrient transfer and synthesis mechanisms, which might offer a possible approach for regulation of milk synthesis in the future.

Preparing and Analyzing Expressed Sequence Tags (ESTs) Library for the Mammary Tissue of Local Turkish Kivircik Sheep

Kivircik sheep is an important local Turkish sheep according to its meat quality and milk productivity. The aim of this study was to analyze gene expression profiles of both prenatal and postnatal stages for the Kivircik sheep. Therefore, two different cDNA libraries, which were taken from the same Kivircik sheep mammary gland tissue at prenatal and postnatal stages, were constructed. Total 3072 colonies which were randomly selected from the two libraries were sequenced for developing a sheep ESTs collection. We used Phred/Phrap computer programs for analysis of the raw EST and readable EST sequences were assembled with the CAP3 software. Putative functions of all unique sequences and statistical analysis were determined by Geneious software. Total 422 ESTs have over 80% similarity to known sequences of other organisms in NCBI classified by Panther database for the Gene Ontology (GO) category. By comparing gene expression profiles, we observed some putative genes that may be relative to reproductive performance or play important roles in milk synthesis and secretion. A total of 2414 ESTs have been deposited to the NCBI GenBank database (GW996847–GW999260). EST data in this study have provided a new source of information to functional genome studies of sheep.

Transcriptomic analysis reveals key regulators of mammogenesis and the pregnancy-lactation cycle

An organ unique to mammals, the mammary gland develops 90% of its mass after birth and experiences the pregnancylactation-involution cycle (PL cycle) during reproduction. To understand mammogenesis at the transcriptomic level and using a ribo-minus RNA-seq protocol, we acquired greater than 50 million reads each for the mouse mammary gland during pregnancy (day 12 of pregnancy), lactation (day 14 of lactation), and involution (day 7 of involution). The pregnancy-, lactation- and involution-related sequencing reads were assembled into 17344, 10160, and 13739 protein-coding transcripts and 1803, 828, and 1288 non-coding RNAs (ncRNAs), respectively. Differentially expressed genes (DEGs) were defined in the three samples, which comprised 4843 DEGs (749 up-regulated and 4094 down-regulated) from pregnancy to lactation and 4926 DEGs (4706 up-regulated and 220 down-regulated) from lactation to involution. Besides the obvious and substantive up- and down-regulation of the DEGs, we observe that lysosomal enzymes were highly expressed and that their expression coincided with milk secretion. Further analysis of transcription factors such as Trps1, Gtf2i, Tcf7l2, Nupr1, Vdr, Rb1, and Aebp1, and ncRNAs such as mir-125b, Let7, mir-146a, and mir-15 has enabled us to identify key regulators in mammary gland development and the PL cycle.