RNA sequencing of the human milk fat layer transcriptome reveals distinct gene expression profiles at three stages of lactation
Danielle G Lemay, Olivia A Ballard, Maria A Hughes, Ardythe L Morrow, Nelson D Horseman, Laurie A Nommsen-Rivers, Danielle G Lemay, Olivia A Ballard, Maria A Hughes, Ardythe L Morrow, Nelson D Horseman, Laurie A Nommsen-Rivers
Abstract
Aware of the important benefits of human milk, most U.S. women initiate breastfeeding but difficulties with milk supply lead some to quit earlier than intended. Yet, the contribution of maternal physiology to lactation difficulties remains poorly understood. Human milk fat globules, by enveloping cell contents during their secretion into milk, are a rich source of mammary cell RNA. Here, we pair this non-invasive mRNA source with RNA-sequencing to probe the milk fat layer transcriptome during three stages of lactation: colostral, transitional, and mature milk production. The resulting transcriptomes paint an exquisite portrait of human lactation. The resulting transcriptional profiles cluster not by postpartum day, but by milk Na:K ratio, indicating that women sampled during similar postpartum time frames could be at markedly different stages of gene expression. Each stage of lactation is characterized by a dynamic range (10(5)-fold) in transcript abundances not previously observed with microarray technology. We discovered that transcripts for isoferritins and cathepsins are strikingly abundant during colostrum production, highlighting the potential importance of these proteins for neonatal health. Two transcripts, encoding β-casein (CSN2) and α-lactalbumin (LALBA), make up 45% of the total pool of mRNA in mature lactation. Genes significantly expressed across all stages of lactation are associated with making, modifying, transporting, and packaging milk proteins. Stage-specific transcripts are associated with immune defense during the colostral stage, up-regulation of the machinery needed for milk protein synthesis during the transitional stage, and the production of lipids during mature lactation. We observed strong modulation of key genes involved in lactose synthesis and insulin signaling. In particular, protein tyrosine phosphatase, receptor type, F (PTPRF) may serve as a biomarker linking insulin resistance with insufficient milk supply. This study provides the methodology and reference data set to enable future targeted research on the physiological contributors of sub-optimal lactation in humans.
Conflict of interest statement
Competing Interests: The authors have declared that no competing interests exist.
Figures
References
- Ip S, Chung M, Raman G, Chew P, Magula N, et al. (2007) Breastfeeding and maternal and infant health outcomes in developed countries. Agency for Healthcare Research and Quality.
- Hauck FR, Thompson JM, Tanabe KO, Moon RY, Vennemann MM (2011) Breastfeeding and reduced risk of sudden infant death syndrome: a meta-analysis. Pediatrics 128: 103–110.
- Gunderson EP, Jacobs DR Jr, Chiang V, Lewis CE, Feng J, et al. (2010) Duration of lactation and incidence of the metabolic syndrome in women of reproductive age according to gestational diabetes mellitus status: a 20-Year prospective study in CARDIA (Coronary Artery Risk Development in Young Adults). Diabetes 59: 495–504.
- US Department of Health and Human Services (2011) The Surgeon General’s Call to Action to Support Breastfeeding. Washington, DC: Office of the Surgeon General.
- Institute of Medicine Committee on Obesity Prevention Policies for Young Children (2011) Early Childhood Obesity Prevention Policies: Goals, Recommendations, and Potential Actions. Washington, D.C.: Institute of Medicine.
- Gartner LM, Morton J, Lawrence RA, Naylor AJ, O’Hare D, et al. (2005) Breastfeeding and the use of human milk. Pediatrics 115: 496–506.
- Li R, Fein SB, Chen J, Grummer-Strawn LM (2008) Why mothers stop breastfeeding: mothers’ self-reported reasons for stopping during the first year. Pediatrics 122 Suppl 2S69–76.
- Declercq E, Labbok MH, Sakala C, O’Hara M (2009) Hospital practices and women’s likelihood of fulfilling their intention to exclusively breastfeed. Am J Public Health 99: 929–935.
- Patton S, Huston GE (1988) Incidence and characteristics of cell pieces on human milk fat globules. Biochim Biophys Acta 965: 146–153.
- Maningat PD, Sen P, Rijnkels M, Sunehag AL, Hadsell DL, et al. (2009) Gene expression in the human mammary epithelium during lactation: the milk fat globule transcriptome. Physiol Genomics 37: 12–22.
- Wang Z, Gerstein M, Snyder M (2009) RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 10: 57–63.
- Nommsen-Rivers LA, Chantry CJ, Peerson JM, Cohen RJ, Dewey KG (2010) Delayed onset of lactogenesis among first-time mothers is related to maternal obesity and factors associated with ineffective breastfeeding. Am J Clin Nutr 92: 574–584.
- Nommsen-Rivers LA, Dolan LM, Huang B (2011) Timing of Stage II Lactogenesis Is Predicted by Antenatal Metabolic Health in a Cohort of Primiparas. Breastfeed Med.
- Neville MC, Allen JC, Archer PC, Casey CE, Seacat J, et al. (1991) Studies in human lactation: milk volume and nutrient composition during weaning and lactogenesis. Am J Clin Nutr 54: 81–92.
- Filteau SM, Lietz G, Mulokozi G, Bilotta S, Henry CJ, et al. (1999) Milk cytokines and subclinical breast inflammation in Tanzanian women: effects of dietary red palm oil or sunflower oil supplementation. Immunology 97: 595–600.
- Patton S, Huston GE (1986) A method for isolation of milk fat globules. Lipids 21: 170–174.
- Roncada P, Stipetic LH, Bonizzi L, Burchmore RJ, Kennedy MW (2013) Proteomics as a tool to explore human milk in health and disease. J Proteomics.
- Jensen RG (1995) Handbook of Milk Composition; Jensen RG, editor. San Diego, CA: Academic Press, Inc., p. 919.
- Jackson JG, Janszen DB, Lonnerdal B, Lien EL, Pramuk KP, et al. (2004) A multinational study of alpha-lactalbumin concentrations in human milk. J Nutr Biochem 15: 517–521.
- Manso MA, Miguel M, Lopez-Fandino R (2007) Application of capillary zone electrophoresis to the characterisation of the human milk protein profile and its evolution throughout lactation. J Chromatogr A 1146: 110–117.
- Arosio P, Ponzone A, Ferrero R, Renoldi I, Levi S (1986) Characteristics of ferritins in human milk secretions: similarities to serum and tissue isoferritins. Clin Chim Acta 161: 201–208.
- Barboza M, Pinzon J, Wickramasinghe S, Froehlich JW, Moeller I, et al. (2012) Glycosylation of human milk lactoferrin exhibits dynamic changes during early lactation enhancing its role in pathogenic bacteria-host interactions. Mol Cell Proteomics 11: M111 015248.
- Christensen B, Schack L, Klaning E, Sorensen ES (2010) Osteopontin is cleaved at multiple sites close to its integrin-binding motifs in milk and is a novel substrate for plasmin and cathepsin D. J Biol Chem. 285: 7929–7937.
- Dallas DC, Guerrero A, Khaldi N, Castillo PA, Martin WF, et al. (2013) Extensive in vivo Human Milk Peptidomics Reveals Specific Proteolysis Yielding Protective Antimicrobial Peptides. J Proteome Res 12: 2295–2304.
- Edmond KM, Zandoh C, Quigley MA, Amenga-Etego S, Owusu-Agyei S, et al. (2006) Delayed breastfeeding initiation increases risk of neonatal mortality. Pediatrics 117: e380–386.
- Garrido D, Nwosu C, Ruiz-Moyano S, Aldredge D, German JB, et al. (2012) Endo-beta-N-acetylglucosaminidases from infant gut-associated bifidobacteria release complex N-glycans from human milk glycoproteins. Mol Cell Proteomics 11: 775–785.
- Sanchez-Juanes F, Alonso JM, Zancada L, Hueso P (2009) Glycosphingolipids from bovine milk and milk fat globule membranes: a comparative study. Adhesion to enterotoxigenic Escherichia coli strains. Biol Chem 390: 31–40.
- Lemay DG, Neville MC, Rudolph MC, Pollard KS, German JB (2007) Gene regulatory networks in lactation: identification of global principles using bioinformatics. BMC Syst Biol 1: 56.
- Bionaz M, Loor JJ (2011) Gene networks driving bovine mammary protein synthesis during the lactation cycle. Bioinform Biol Insights 5: 83–98.
- Bernier-Dodier P, Delbecchi L, Wagner GF, Talbot BG, Lacasse P (2010) Effect of milking frequency on lactation persistency and mammary gland remodeling in mid-lactation cows. J Dairy Sci 93: 555–564.
- Menzies KK, Lefevre C, Macmillan KL, Nicholas KR (2009) Insulin regulates milk protein synthesis at multiple levels in the bovine mammary gland. Funct Integr Genomics 9: 197–217.
- Berlato C, Doppler W (2009) Selective response to insulin versus insulin-like growth factor-I and -II and up-regulation of insulin receptor splice variant B in the differentiated mouse mammary epithelium. Endocrinology 150: 2924–2933.
- Belfiore A, Frasca F, Pandini G, Sciacca L, Vigneri R (2009) Insulin receptor isoforms and insulin receptor/insulin-like growth factor receptor hybrids in physiology and disease. Endocr Rev 30: 586–623.
- Koren S, Fantus IG (2007) Inhibition of the protein tyrosine phosphatase PTP1B: potential therapy for obesity, insulin resistance and type-2 diabetes mellitus. Best Pract Res Clin Endocrinol Metab 21: 621–640.
- Hill PD, Humenick SS (1994) The occurrence of breast engorgement. J Hum Lact 10: 79–86.
- Garcia-San Frutos M, Fernandez-Agullo T, Carrascosa JM, Horrillo D, Barrus MT, et al. (2012) Involvement of protein tyrosine phosphatases and inflammation in hypothalamic insulin resistance associated with ageing: effect of caloric restriction. Mech Ageing Dev 133: 489–497.
- Chagnon MJ, Uetani N, Tremblay ML (2004) Functional significance of the LAR receptor protein tyrosine phosphatase family in development and diseases. Biochem Cell Biol 82: 664–675.
- Macheda ML, Williams ED, Best JD, Wlodek ME, Rogers S (2003) Expression and localisation of GLUT1 and GLUT12 glucose transporters in the pregnant and lactating rat mammary gland. Cell Tissue Res 311: 91–97.
- Shennan DB, Peaker M (2000) Transport of milk constituents by the mammary gland. Physiol Rev 80: 925–951.
- Shields SL, Rezamand P, Sevier DL, Seo KS, Price W, et al. (2011) Effects of increased milking frequency for the first 21 days post partum on selected measures of mammary gland health, milk yield and milk composition. J Dairy Res 78: 301–307.
- Wall EH, McFadden TB (2010) The effects of milk removal or four-times-daily milking on mammary expression of genes involved in the insulin-like growth factor-I axis. J Dairy Sci 93: 4062–4070.
- Littlejohn MD, Walker CG, Ward HE, Lehnert KB, Snell RG, et al. (2010) Effects of reduced frequency of milk removal on gene expression in the bovine mammary gland. Physiol Genomics 41: 21–32.
- Geraghty SR, List BA, Morrow GB (2010) Guidelines for establishing a donor human milk depot. J Hum Lact 26: 49–52.
- Chapman DJ, Perez-Escamilla R (1999) Identification of risk factors for delayed onset of lactation. J Am Diet Assoc 99: 450–454; quiz 455–456.
- Dewey KG, Nommsen-Rivers LA, Heinig MJ, Cohen RJ (2003) Risk factors for suboptimal infant breastfeeding behavior, delayed onset of lactation, and excess neonatal weight loss. Pediatrics 112: 607–619.
- Matsuda M, DeFronzo RA (1999) Insulin sensitivity indices obtained from oral glucose tolerance testing: comparison with the euglycemic insulin clamp. Diabetes Care 22: 1462–1470.
- Retnakaran R, Qi Y, Goran MI, Hamilton JK (2009) Evaluation of proposed oral disposition index measures in relation to the actual disposition index. Diabet Med 26: 1198–1203.
- Trapnell C, Pachter L, Salzberg SL (2009) TopHat: discovering splice junctions with RNA-Seq. Bioinformatics 25: 1105–1111.
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, et al. (2010) Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol 28: 511–515.
- Flicek P, Amode MR, Barrell D, Beal K, Brent S, et al. (2012) Ensembl 2012. Nucleic Acids Res 40: D84–90.
- Anders S, Huber W (2010) Differential expression analysis for sequence count data. Genome Biol 11: R106.
- Kumar L, M EF (2007) Mfuzz: a software package for soft clustering of microarray data. Bioinformation 2: 5–7.
- Huang da W, Sherman BT, Tan Q, Kir J, Liu D, et al. (2007) DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res 35: W169–175.
- Huang da W, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 44–57.
- Bionaz M, Loor JJ (2007) Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle. Physiol Genomics 29: 312–319.
Source: PubMed