Identification of Putative Ortholog Gene Blocks Involved in Gestant and Lactating Mammary Gland Development: A Rodent Cross-Species Microarray Transcriptomics Approach
The mammary gland (MG) undergoes functional and metabolic changes during the transition from pregnancy to lactation, possibly by regulation of conserved genes. The objective was to elucidate orthologous genes, chromosome clusters and putative conserved transcriptional modules during MG development. We analyzed expression of 22,000 transcripts using murine microarrays and RNA samples of MG from virgin, pregnant, and lactating rats by cross-species hybridization. We identified 521 transcripts differentially expressed; upregulated in early (78%) and midpregnancy (89%) and early lactation (64%), but downregulated in mid-lactation (61%). Putative orthologous genes were identified. We mapped the altered genes to orthologous chromosomal locations in human and mouse. Eighteen sets of conserved genes associated with key cellular functions were revealed and conserved transcription factor binding site search entailed possible coregulation among all eight block sets of genes. This study demonstrates that the use of heterologous array hybridization for screening of orthologous gene expression from rat revealed sets of conserved genes arranged in chromosomal order implicated in signaling pathways and functional ontology. Results demonstrate the utilization power of comparative genomics and prove the feasibility of using rodent microarrays to identification of putative coexpressed orthologous genes involved in the control of human mammary gland development. 1. Introduction Mammals are the only animals that secrete a complex fluid from an elaborated skin gland to provide both innate protection and nourishment for their newborn. There are more than 4,000 species of mammals with striking similarities in the structure and function of their mammary glands as well as in their unique milk components such as the caseins, α-lactalbumin, lactoferrin, lactose, and milk fat. Nevertheless, variations are exhibited in the arrangement and numbers of mammary gland, milk composition, and suckling strategies. Mammary gland development begins at puberty and is maintained throughout pregnancy until lactation. During these last stages, development compromises numerous overlapping programs such as branching morphogenesis, inductive stromal-epithelial interactions, programmed cell death, extracellular matrix remodeling, and hormone action [1]. Current knowledge of the molecular regulation of mammary development and lactation has largely been derived from the dissection of signaling networks in cell culture systems and phenotypic characterization of genetically altered mice as well as
References
[1]
T. W. Hale and P. Hartmann, Textbook on Human Lactation, Hale Publishing, Amarillo, Tex, USA, 2007.
[2]
D. G. Lemay, M. C. Neville, M. C. Rudolph, K. S. Pollard, and J. B. German, “Gene regulatory networks in lactation: identification of global principles using bioinformatics,” BMC Systems Biology, vol. 1, article 56, 2007.
[3]
X. J. Zhou and G. Gibson, “Cross-species comparison of genome-wide expression patterns,” Genome Biology, vol. 5, no. 7, article 232, 2004.
[4]
R. A. Gibbs, G. M. Weinstock, M. L. Metzker, et al., “Genome sequence of the Brown Norway rat yields insights into mammalian evolution,” Nature, vol. 428, no. 6982, pp. 493–521, 2004.
[5]
D. L. Wheeler, T. Barrett, D. A. Benson et al., “Database resources of the National Center for Biotechnology Information,” Nucleic Acids Research, vol. 35, no. 1, pp. D5–D12, 2007.
[6]
H. Koltai and C. Weingarten-Baror, “Specificity of DNA microarray hybridization: characterization, effectors and approaches for data correction,” Nucleic Acids Research, vol. 36, no. 7, pp. 2395–2405, 2008.
[7]
M. Rodríguez-Cruz, R. Sánchez, M. Bernabe-Garcia, J. Maldonado, M. del Prado, and M. López-Alarcón, “Effect of dietary levels of corn oil on maternal arachidonic acid synthesis and fatty acid composition in lactating rats,” Nutrition, vol. 25, no. 2, pp. 209–215, 2009.
[8]
P. Chomczynski and N. Sacchi, “Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction,” Analytical Biochemistry, vol. 162, no. 1, pp. 156–159, 1987.
[9]
J. H. Kim, “Chapter 8: biological knowledge assembly and interpretation,” PLOS Computational Biology, vol. 8, no. 12, Article ID e1002858, 2012.
[10]
M. C. Rudolph, J. L. McManaman, L. Hunter, T. Phang, and M. C. Neville, “Functional development of the mammary gland: Uuse of expression rofiling and trajectory clustering to reveal changes in gene expression during pregnancy, lactation, and involution,” Journal of Mammary Gland Biology and Neoplasia, vol. 8, no. 3, pp. 287–307, 2003.
[11]
T. F. DeLuca, I.-H. Wu, J. Pu et al., “Roundup: a multi-genome repository of orthologs and evolutionary distances,” Bioinformatics, vol. 22, no. 16, pp. 2044–2046, 2006.
[12]
P. Flicek, B. L. Aken, K. Beal et al., “Ensembl 2008,” Nucleic Acids Research, vol. 36, no. 1, pp. D707–D714, 2008.
[13]
S. M. Anderson, M. C. Rudolph, J. L. McManaman, and M. C. Neville, “Key stages in mammary gland development. Secretory activation in the mammary gland: it's not just about milk protein synthesis!,” Breast Cancer Research, vol. 9, no. 1, article 204, 2007.
[14]
D. J. Flint, M. Boutinaud, C. B. A. Whitelaw, G. J. Allan, and A. F. Kolb, “Prolactin inhibits cell loss and decreases matrix metalloproteinase expression in the involuting mouse mammary gland but fails to prevent cell loss in the mammary glands of mice expressing IGFBP-5 as a mammary transgene,” Journal of Molecular Endocrinology, vol. 36, no. 3, pp. 435–448, 2006.
[15]
Y. Ning, B. Hoang, A. G. P. Schuller et al., “Delayed mammary gland involution in mice with mutation of the insulin-like growth factor binding protein 5 gene,” Endocrinology, vol. 148, no. 5, pp. 2138–2147, 2007.
[16]
A. Morani, M. Warner, and J.-?. Gustafsson, “Biological functions and clinical implications of oestrogen receptors alfa and beta in epithelial tissues,” Journal of Internal Medicine, vol. 264, no. 2, pp. 128–142, 2008.
[17]
S. B. Aicha, J. Lessard, M. Pelletier, A. Fournier, E. Calvo, and C. Labrie, “Transcriptional profiling of genes that are regulated by the endoplasmic reticulum-bound transcription factor AIbZIP/CREB3L4 in prostate cells,” Physiological Genomics, vol. 31, no. 2, pp. 295–305, 2007.
[18]
J. ?lander, J. Lengqvist, P. J. Holm et al., “Microsomal glutathione transferase 1 exhibits one-third-of-the-sites-reactivity towards glutathione,” Archives of Biochemistry and Biophysics, vol. 487, no. 1, pp. 42–48, 2009.
[19]
H. Nakshatri and S. Badve, “FOXA1 in breast cancer,” Expert Reviews in Molecular Medicine, vol. 11, p. e8, 2009.
[20]
D. R. Magnan, D. V. Spacek, N. Ye, Y.-C. Lu, and T. R. King, “The male sterility and histoincompatibility (mshi) mutation in mice is a natural variant of microtubule-associated protein 7 (Mtap7),” Molecular Genetics and Metabolism, vol. 97, no. 2, pp. 155–162, 2009.
[21]
C. T. Rankin, T. Bunton, A. M. Lawler, and S.-J. Lee, “Regulation of left-right patterning in mice by growth/differentiation factor-1,” Nature Genetics, vol. 24, no. 3, pp. 262–265, 2000.
[22]
J. C. Martinez, M. M. Müller, H. Turley et al., “Nuclear and membrane expression of the angiogenesis regulator delta-like ligand 4 (DLL4) in normal and malignant human tissues,” Histopathology, vol. 54, no. 5, pp. 598–606, 2009.
[23]
C. Zhu, Y.-R. Qin, D. Xie et al., “Characterization of tumor suppressive function of P300/CBP-associated factor at frequently deleted region 3p24 in esophageal squamous cell carcinoma,” Oncogene, vol. 28, no. 31, pp. 2821–2828, 2009.
[24]
A. C. Zemke, J. C. Snyder, B. L. Brockway et al., “Molecular staging of epithelial maturation using secretory cell-specific genes as markers,” American Journal of Respiratory Cell and Molecular Biology, vol. 40, no. 3, pp. 340–348, 2009.
[25]
C. Nicolae, Y.-P. Ko, N. Miosge et al., “Abnormal collagen fibrils in cartilage of matrilin-1/matrilin-3-deficient mice,” Journal of Biological Chemistry, vol. 282, no. 30, pp. 22163–22175, 2007.
[26]
D. L. Hadsell, W. Olea, N. Lawrence et al., “Decreased lactation capacity and altered milk composition in insulin receptor substrate null mice is associated with decreased maternal body mass and reduced insulin-dependent phosphorylation of mammary Akt,” Journal of Endocrinology, vol. 194, no. 2, pp. 327–336, 2007.
[27]
C. J. Watson, “Immune cell regulators in mouse mammary development and involution,” Journal of Animal Science, vol. 87, no. 13, pp. 35–42, 2009.
[28]
I. Dalkilic, J. Schienda, T. G. Thompson, and L. M. Kunkel, “Loss of filaminC (FLNc) results in severe defects in myogenesis and myotube structure,” Molecular and Cellular Biology, vol. 26, no. 17, pp. 6522–6534, 2006.
[29]
Z. Luan, Y. Liu, T. J. Stuhlmiller, J. Marquez, and M. I. García-Castro, “SUMOylation of Pax7 is essential for neural crest and muscle development,” Cellular and Molecular Life Sciences, vol. 70, no. 10, pp. 1793–1806, 2013.
[30]
R. C. Wek and D. R. Cavener, “Translational control and the unfolded protein response,” Antioxidants and Redox Signaling, vol. 9, no. 12, pp. 2357–2371, 2007.
[31]
M. C. Thomas and C. M. Chiang, “The general transcription machinery and general cofactors,” Critical Reviews in Biochemistry and Molecular Biology, vol. 41, no. 3, pp. 105–178, 2006.
[32]
I. A. Belyantseva, B. J. Perrin, K. J. Sonnemann et al., “γ-Actin is required for cytoskeletal maintenance but not development,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 24, pp. 9703–9708, 2009.
[33]
N. Saito, Y. Suhara, D. Abe et al., “Synthesis of 2α-propoxy-1α,25-dihydroxyvitamin D3 and comparison of its metabolism by human CYP24A1 and rat CYP24A1,” Bioorganic and Medicinal Chemistry, vol. 17, no. 13, pp. 4296–4301, 2009.
[34]
C. J. Watson and K. Neoh, “The Stat family of transcription factors have diverse roles in mammary gland development,” Seminars in Cell and Developmental Biology, vol. 19, no. 4, pp. 401–406, 2008.
[35]
Y. Li, Y. Chang, L. Zhang et al., “High glucose upregulates pantothenate kinase 4 (PanK4) and thus affects M2-type pyruvate kinase (Pkm2),” Molecular and Cellular Biochemistry, vol. 277, no. 1-2, pp. 117–125, 2005.
[36]
D. Horst, S. Ustanina, C. Sergi et al., “Comparative expression analysis of Pax3 and Pax7 during mouse myogenesis,” International Journal of Developmental Biology, vol. 50, no. 1, pp. 47–54, 2006.
[37]
H. Fang, W. Tong, R. Perkins et al., “Bioinformatics approaches for cross-species liver cancer analysis based on microarray gene expression profiling,” BMC Bioinformatics, vol. 6, no. 2, article S6, 2005.
[38]
M. D. Aupperlee, A. A. Drolet, S. Durairaj, W. Wang, R. C. Schwartz, and S. Z. Haslam, “Strain-specific differences in the mechanisms of progesterone regulation of murine mammary gland development,” Endocrinology, vol. 150, no. 3, pp. 1485–1494, 2009.
[39]
C. J. Watson and W. T. Khaled, “Mammary development in the embryo and adult: a journey of morphogenesis and commitment,” Development, vol. 135, no. 6, pp. 995–1003, 2008.
[40]
M. S. Hestand, M. van Galen, M. P. Villerius, G.-J. B. van Ommen, J. T. den Dunnen, and P. A. C. Hoen, “CORE_TF: a user-friendly interface to identify evolutionary conserved transcription factor binding sites in sets of co-regulated genes,” BMC Bioinformatics, vol. 9, article 495, 2008.
[41]
C. F. Singer, G. Hudelist, W. Lamm, R. Mueller, K. Czerwenka, and E. Kubista, “Expression of tyrosine kinases in human malignancies as potential targets for kinase-specific inhibitors,” Endocrine-Related Cancer, vol. 11, no. 4, pp. 861–869, 2004.
[42]
S. C. P. Renn, N. Aubin-Horth, and H. A. Hofmann, “Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray,” BMC Genomics, vol. 5, article 42, 2004.
[43]
C. Bar-Or, M. Bar-Eyal, T. Z. Gal, Y. Kapulnik, H. Czosnek, and H. Koltai, “Derivation of species-specific hybridization-like knowledge out of cross-species hybridization results,” BMC Genomics, vol. 7, article 110, 2006.
[44]
J. Adjaye, R. Herwig, D. Herrmann et al., “Cross-species hybridisation of human and bovine orthologous genes on high density cDNA microarrays,” BMC Genomics, vol. 5, article 83, 2004.
[45]
M. D. Kane, T. A. Jatkoe, C. R. Stumpf, J. Lu, J. D. Thomas, and S. J. Madore, “Assessment of the sensitivity and specificity of oligonucleotide (50mer) microarrays,” Nucleic Acids Research, vol. 28, no. 22, pp. 4552–4557, 2000.
[46]
K. A. Frazer, L. Elnitski, D. M. Church, I. Dubchak, and R. C. Hardison, “Cross-species sequence comparisons: a review of methods and available resources,” Genome Research, vol. 13, no. 1, pp. 1–12, 2003.