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

Single-Cell Transcription Mapping of Murine and Human Mammary Organoids Responses to Female Hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Single-cell transcription mapping of murine and human mammary organoids responses to female hormones

During female adolescence and pregnancy, rising levels of hormones result in a cyclic source of signals that control the development of mammary tissue. While such alterations are well understood from a whole-gland perspective, the alterations that such hormones bring to organoid cultures derived from mammary glands have yet to be fully mapped. This is of special importance given that organoids are considered suitable systems to understand cross species breast development. Here we utilized single-cell transcriptional profiling to delineate responses of murine and human normal breast organoid systems to female hormones across evolutionary distinct species. Collectively, our study represents a molecular atlas of epithelial dynamics in response to estrogen and pregnancy hormones.

Exosome-Derived microRNAs from Mouthrinse Have the Potential to Be Novel Biomarkers for Sjögren Syndrome

Sjögren syndrome (SS) is diagnosed based on invasive tissue biopsies and blood sampling. Therefore, a novel non-invasive and simple inspection diagnostic marker of SS is required. Here, we identified exosome-derived microRNAs (miRNAs) as biomarkers for SS using non-invasive mouthrinse samples collected from patients with SS and healthy volunteers. We compared miRNAs derived from exosomes in mouthrinse samples from the two groups using microarrays and real-time polymerase chain reaction (PCR) and identified 12 miRNAs as biomarker candidates. The expression ratios of four miRNAs were significantly increased in the SS group compared to the control group. Logistic regression analysis revealed a more significant influence of miR-1290 and let-7b-5p in the SS group than that in the control group. We combined these miRNAs to create a diagnostic prediction formula using logistic regression analysis. The combination of miR-1290 and let-7b-5p distinguished SS from the control samples with an AUC, sensitivity, specificity, positive predictive value, and negative predictive value of 0.856, 91.7%, 83.3%, 84.6%, and 90.9%, respectively. These results indicated that an increased ratio of these miRNAs could serve as a novel and non-invasive diagnostic marker for SS. This is the first report of diagnosis and screening of SS by adopting a non-invasive method using mouthrinse.

The effect of single-nucleotide polymorphism in the promoter region of bovine alpha-lactalbumin (LALBA) gene on LALBA expression in milk cells and milk traits of cows

Polymorphisms of milk protein genes have been proposed as candidate markers for dairy production traits in cattle. In the present study, a polymorphism was detected in the 5'-flanking (promoter) region of the bovine alpha-lactalbumin (LALBA) gene, a T/C transition located at nucleotide -1,001 relative to the transcription start site g.-1001T > C (NC_037332.1:g.31183170T > C), which is recognizable with PstI restriction endonuclease. In silico analyses showed that this mutation created novel retinoid X receptor alpha and vitamin D receptor transcription factor binding sites. Real-time PCR found that cows with different genetic variants of the promoter demonstrated different levels of expression of LALBA mRNA in milk somatic cells (MSCs). The TT genotype cows demonstrated low expression, whereas those with CT demonstrated much higher expression (P < 0.05). ELISA analysis found milk LALBA protein levels also differed between the TT and CT cows (P < 0.05) and that these levels were not correlated with the mRNA abundance in MSC. Association analysis found that the g.-1001T > C polymorphism in the promoter region of the LALBA gene influenced milk production traits in Polish Holstein-Friesian cows. High daily milk yield and dry matter yield, and high lactose yield and concentration were associated with the TT genotype. The TT genotype cows also had a lower number of somatic cells in the milk, considered as an indicator of udder health status. Therefore, the TT genotype could be more desirable from the breeder's perspective.

Characterization of Gene Expression Signatures for the Identification of Cellular Heterogeneity in the Developing Mammary Gland

The developing mammary gland depends on several transcription-dependent networks to define cellular identities and differentiation trajectories. Recent technological advancements that allow for single-cell profiling of gene expression have provided an initial picture into the epithelial cellular heterogeneity across the diverse stages of gland maturation. Still, a deeper dive into expanded molecular signatures would improve our understanding of the diversity of mammary epithelial and non-epithelial cellular populations across different tissue developmental stages, mouse strains and mammalian species. Here, we combined differential mammary gland fractionation approaches and transcriptional profiles obtained from FACS-isolated mammary cells to improve our definitions of mammary-resident, cellular identities at the single-cell level. Our approach yielded a series of expression signatures that illustrate the heterogeneity of mammary epithelial cells, specifically those of the luminal fate, and uncovered transcriptional changes to their lineage-defined, cellular states that are induced during gland development. Our analysis also provided molecular signatures that identified non-epithelial mammary cells, including adipocytes, fibroblasts and rare immune cells. Lastly, we extended our study to elucidate expression signatures of human, breast-resident cells, a strategy that allowed for the cross-species comparison of mammary epithelial identities. Collectively, our approach improved the existing signatures of normal mammary epithelial cells, as well as elucidated the diversity of non-epithelial cells in murine and human breast tissue. Our study provides a useful resource for future studies that use single-cell molecular profiling strategies to understand normal and malignant breast development.

Identification of differentially expressed miRNAs in serum extracellular vesicles (EVs) of Kazakh sheep at early pregnancy

MiRNAs-containing extracellular vesicles (EVs) possess the unique function of mediating intercellular communication and participating in many biological processes such as post-transcriptional gene regulation of embryo implantation and placental development. In the present study, Illumina small-RNA sequencing was used to identify differentially expressed (DE) miRNAs in serum EVs of pregnant (P) and non-pregnant (NP) Kazakh sheep at Day 17 from mating. The specifically and differentially expressed miRNAs at early pregnancy in sheep were verified by using RT-PCR. The target genes of DE miRNAs were predicted by bioinformatics software, and the functional and pathway enrichment analysis was performed on Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) terms. A total of 562 miRNAs (210 novel miRNAs) were identified by sequencing, of which 57 miRNAs were differentially expressed, 49 were up-regulated, 8 were down-regulated and 22 novel miRNAs were specifically expressed in the pregnant sheep. Eight highly expressed known miRNA (miR-378-3p, miR-320-3p, miR-22-3p, let-7b, miR-423-3p, miR-221, miR-296-3p, miR-147-3p) in pregnant group were down-regulated in the control group. miRNAs-containing pregnancy-related terms and regulatory pathways regulation were enriched using both GO and KEGG analyses. Moreover, we also envisioned a miRNA-mRNA interaction network to understand the function of miRNAs involved in the early pregnancy serum regulatory network. The results of RT-PCR verification confirmed the reliability of small-RNA sequencing. Among them, miR-22-3p and miR-378-3p were significantly differentially expressed (DE) between pregnant sheep and non-pregnant group (p <  0.01). The site at which oar-miR-22-3p binds MAPK3 was determined with a dual-luciferase system. This is the first integrated analysis of the expression profiles of EV-miRNAs and their targets during early pregnancy in ewes. These data identify key miRNAs that influence the implantation of sheep in the early stage of pregnancy, and provide theoretical basis for further molecular regulatory mechanisms research.

Epigenetics: New Insights into Mammary Gland Biology

The mammary gland undergoes important anatomical and physiological changes from embryogenesis through puberty, pregnancy, lactation and involution. These steps are under the control of a complex network of molecular factors, in which epigenetic mechanisms play a role that is increasingly well described. Recently, studies investigating epigenetic modifications and their impacts on gene expression in the mammary gland have been performed at different physiological stages and in different mammary cell types. This has led to the establishment of a role for epigenetic marks in milk component biosynthesis. This review aims to summarize the available knowledge regarding the involvement of the four main molecular mechanisms in epigenetics: DNA methylation, histone modifications, polycomb protein activity and non-coding RNA functions.

AEBP1 is a Novel Oncogene: Mechanisms of Action and Signaling Pathways

Adipocyte enhancer-binding protein 1 (AEBP1) is a transcriptional repressor involved in the regulation of critical biological processes including adipogenesis, mammary gland development, inflammation, macrophage cholesterol homeostasis, and atherogenesis. Several years ago, we first reported the ability of AEBP1 to exert a positive control over the canonical NF-κB pathway. Indeed, AEBP1 positively regulates NF-κB activity via its direct interaction with IκBα, a key NF-κB inhibitor. AEBP1 overexpression results in uncontrollable activation of NF-κB, which may have severe pathogenic outcomes. Recently, the regulatory relationship between AEBP1 and NF-κB pathway has been of great interest to many researchers primarily due to the implication of NF-κB signaling in critical cellular processes such as inflammation and cancer. Since constitutive activation of NF-κB is widely implicated in carcinogenesis, AEBP1 overexpression is associated with tumor development and progression. Recent studies sought to explore the effects of the overexpression of AEBP1, as a potential oncogene, in different types of cancer. In this review, we analyze the effects of AEBP1 overexpression in a variety of malignancies (e.g., breast cancer, glioblastoma, bladder cancer, gastric cancer, colorectal cancer, ovarian cancer, and skin cancer), with a specific focus on the AEBP1-mediated control over the canonical NF-κB pathway. We also underscore the ability of AEBP1 to regulate crucial cancer-related events like cell proliferation and apoptosis in light of other key pathways (e.g., PI3K-Akt, sonic hedgehog (Shh), p53, parthanatos (PARP-1), and PTEN). Identifying AEBP1 as a potential biomarker for cancer prognosis may lead to a novel therapeutic target for the prevention and/or treatment of various types of cancer.

Transcriptional profiling of swine mammary gland during the transition from colostrogenesis to lactogenesis using RNA sequencing

Background

Colostrum and milk are essential sources of antibodies and nutrients for the neonate, playing a key role in their survival and growth. Slight abnormalities in the timing of colostrogenesis/lactogenesis potentially threaten piglet survival. To further delineate the genes and transcription regulators implicated in the control of the transition from colostrogenesis to lactogenesis, we applied RNA-seq analysis of swine mammary gland tissue from late-gestation to farrowing. Three 2nd parity sows were used for mammary tissue biopsies on days 14, 10, 6 and 2 before (−) parturition and on day 1 after (+) parturition. A total of 15 mRNA libraries were sequenced on a HiSeq2500 (Illumina Inc.). The Dynamic Impact Approach and the Ingenuity Pathway Analysis were used for pathway analysis and gene network analysis, respectively.

Results

A large number of differentially expressed genes were detected very close to parturition (−2d) and at farrowing (+ 1d). The results reflect the extraordinary metabolic changes in the swine mammary gland once it enters into the crucial phases of lactogenesis and underscore a strong transcriptional component in the control of colostrogenesis. There was marked upregulation of genes involved in synthesis of colostrum and main milk components (i.e. proteins, fat, lactose and antimicrobial factors) with a pivotal role of CSN1S2, LALBA, WAP, SAA2, and BTN1A1. The sustained activation of transcription regulators such as SREBP1 and XBP1 suggested they help coordinate these adaptations.

Conclusions

Overall, the precise timing for the transition from colostrogenesis to lactogenesis in swine mammary gland remains uncharacterized. However, our transcriptomic data support the hypothesis that the transition occurs before parturition. This is likely attributable to upregulation of a wide array of genes including those involved in ‘Protein and Carbohydrate Metabolism’, ‘Immune System’, ‘Lipid Metabolism’, ‘PPAR signaling pathway’ and ‘Prolactin signaling pathway’ along with the activation of transcription regulators controlling lipid synthesis and endoplasmic reticulum biogenesis and stress response.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4719-5) contains supplementary material, which is available to authorized users.

Vitamin D Receptor Signaling and Cancer

The vitamin D receptor (VDR) binds the secosteroid hormone 1,25(OH)₂D₃ with high affinity and regulates gene programs that control a serum calcium levels, as well as cell proliferation and differentiation. A significant focus has been to exploit the VDR in cancer settings. Although preclinical studies have been strongly encouraging, to date clinical trials have delivered equivocal findings that have paused the clinical translation of these compounds. However, it is entirely possible that mining of genomic data will help to refine precisely what are the key anticancer actions of vitamin D compounds and where these can be used most effectively.

MiR-139 suppresses β-casein synthesis and proliferation in bovine mammary epithelial cells by targeting the GHR and IGF1R signaling pathways

BACKGROUND

MicroRNAs have important roles in many biological processes. However, the role of miR-139 in healthy mammary gland remains unclear. The objective of this study was to investigate the effects of miR-139 on lactation in dairy cows.

RESULTS

Here, we found that miR-139 was down-regulated in mid-lactation dairy cow mammary tissues compared with mid-pregnancy tissues. Then, we prioritized two of potential target genes of miR-139 in cow, growth hormone receptor (GHR) and type I insulin-like growth factor receptor (IGF1R) for further functional studies based on their roles in lactation processes. Dual luciferase reporter assays validated direct binding of miR-139 to the 3'- untranslated region (UTR) of GHR and IGF1R. Moreover, over-expression or silencing of miR-139 affected mRNA levels of GHR and IGF1R in cultured bovine mammary epithelial cells (BMECs). Furthermore, over-expression of miR-139 decreased protein levels of β-casein, proliferation in mammary epithelial cell, and the protein levels of IGF1R and key members of the GHR or IGF1R pathways as well, whereas silencing miR-139 produced the opposite result. Among these signal molecules, signal transducer and activator of transcription-5 (STAT5), protein kinase B (also known as AKT1), mammalian target of rapamycin (mTOR), and p70-S6 Kinase (p70S6K) are involed in β-casein synthesis, and Cyclin D1 is involved in cell proliferation. In addition, silencing GHR decreased protein levels of β-casein, IGF1R, and key members of the IGF1R pathway, whereas co-silencing miR-139 and GHR rescued the expression of GHR and reversed GHR silencing effects.

CONCLUSIONS

Our results demonstrate that GHR and IGF1R are target genes of miR-139 in dairy cow. MiR-139 suppresses β-casein synthesis and proliferation in BMECs by targeting the GHR and IGF1R signaling pathways.

Bioinformatic approaches to interrogating vitamin D receptor signaling

Bioinformatics applies unbiased approaches to develop statistically-robust insight into health and disease. At the global, or "20,000 foot" view bioinformatic analyses of vitamin D receptor (NR1I1/VDR) signaling can measure where the VDR gene or protein exerts a genome-wide significant impact on biology; VDR is significantly implicated in bone biology and immune systems, but not in cancer. With a more VDR-centric, or "2000 foot" view, bioinformatic approaches can interrogate events downstream of VDR activity. Integrative approaches can combine VDR ChIP-Seq in cell systems where significant volumes of publically available data are available. For example, VDR ChIP-Seq studies can be combined with genome-wide association studies to reveal significant associations to immune phenotypes. Similarly, VDR ChIP-Seq can be combined with data from Cancer Genome Atlas (TCGA) to infer the impact of VDR target genes in cancer progression. Therefore, bioinformatic approaches can reveal what aspects of VDR downstream networks are significantly related to disease or phenotype.

Molecular conservation of marsupial and eutherian placentation and lactation

Eutherians are often mistakenly termed 'placental mammals', but marsupials also have a placenta to mediate early embryonic development. Lactation is necessary for both infant and fetal development in eutherians and marsupials, although marsupials have a far more complex milk repertoire that facilitates morphogenesis of developmentally immature young. In this study, we demonstrate that the anatomically simple tammar placenta expresses a dynamic molecular program that is reminiscent of eutherian placentation, including both fetal and maternal signals. Further, we provide evidence that genes facilitating fetal development and nutrient transport display convergent co-option by placental and mammary gland cell types to optimize offspring success.

Biosynthesis of milk fat, protein, and lactose: roles of transcriptional and posttranscriptional regulation

The demand for high-quality milk is increasing worldwide. The efficiency of milk synthesis can be improved by taking advantage of the accumulated knowledge of the transcriptional and posttranscriptional regulation of genes coding for proteins involved in the synthesis of fat, protein, and lactose in the mammary gland. Research in this area is relatively new, but data accumulated in the last 10 years provide a relatively clear picture. Milk fat synthesis appears to be regulated, at least in bovines, by an interactive network between SREBP1, PPARγ, and LXRα, with a potential role for other transcription factors, such as Spot14, ChREBP, and Sp1. Milk protein synthesis is highly regulated by insulin, amino acids, and amino acid transporters via transcriptional and posttranscriptional routes, with the insulin-mTOR pathway playing a central role. The transcriptional regulation of lactose synthesis is still poorly understood, but it is clear that glucose transporters play an important role. They can also cooperatively interact with amino acid transporters and the mTOR pathway. Recent data indicate the possibility of nutrigenomic interventions to increase milk fat synthesis by feeding long-chain fatty acids and milk protein synthesis by feeding amino acids. We propose a transcriptional network model to account for all available findings. This model encompasses a complex network of proteins that control milk synthesis with a cross talk between milk fat, protein, and lactose regulation, with mTOR functioning as a central hub.

Circulating miRNA signatures of early pregnancy in cattle

Background

Low fertility remains a leading cause of poor productivity in dairy cattle. In this context, there is significant interest in developing novel tools for accurate early diagnosis of pregnancy. MicroRNAs (miRNAs) are short RNA molecules which are critically involved in regulating gene expression during both health and disease. MiRNAs have been shown to regulate ovarian function, uterine receptivity, embryonic development and placental function. Circulating miRNAs can provide useful biomarkers of tissue function and disease; importantly, differential miRNA profiles have been linked to pregnancy and preeclampsia in humans. This study sought to establish the potential of circulating miRNAs as biomarkers of early pregnancy in cattle.

Results

We applied Illumina small-RNA sequencing to profile miRNAs in plasma samples collected from eight non-pregnant heifers on Days 0, 8 and 16 of the oestrous cycle and 11 heifers on Days 16 and 24 of pregnancy. We sequenced a total of 46 samples and generated 9.2 million miRNA reads per sample. There were no differences in miRNA read abundance between any of the pregnant and non-pregnant time-points (FDR > 0.1). As a complementary approach, we analysed sample pools (3–4 samples/pool) corresponding to Days 0, 8 and 16 of the oestrous cycle and Day 24 of pregnancy (n = 3 pools/group) using Qiagen PCR arrays. A total of 16 miRNAs were differentially expressed (FDR < 0.1) in plasma between pregnant and non-pregnant animals. RT-qPCR validation using the same plasma samples confirmed that miR-26a was differentially upregulated on Day 16 pregnant relative to non-pregnant heifers (1.7-fold; P = 0.043), whereas miR-1249 tended to be upregulated in Day 16 pregnant heifers (1.6-fold; P = 0.081). Further validation in an independent group of heifers confirmed an increase in plasma miR-26a levels during early pregnancy, which was significant only on Day 24 (2.0-fold; P = 0.027).

Conclusions

Through genome-wide analyses we have successfully profiled plasma miRNA populations associated with early pregnancy in cattle. We have identified miR-26a as a potential circulating biomarker of early pregnancy.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2529-1) contains supplementary material, which is available to authorized users.

An epigenetic memory of pregnancy in the mouse mammary gland

Pregnancy is the major modulator of mammary gland activity. It induces a tremendous expansion of the mammary epithelium and the generation of alveolar structures for milk production. Anecdotal evidence from multiparous humans indicates that the mammary gland may react less strongly to the first pregnancy than it does to subsequent pregnancies. Here, we verify that the mouse mammary gland responds more robustly to a second pregnancy, indicating that the gland retains a long-term memory of pregnancy. A comparison of genome-wide profiles of DNA methylation in isolated mammary cell types reveals substantial and long-lasting alterations. Since these alterations are maintained in the absence of the signal that induced them, we term them epigenetic. The majority of alterations in DNA methylation affect sites occupied by the Stat5a transcription factor and mark specific genes that are upregulated during pregnancy. We postulate that the epigenetic memory of a first pregnancy primes the activation of gene expression networks that promote mammary gland function in subsequent reproductive cycles. More broadly, our data indicate that physiological experience can broadly alter epigenetic states, functionally modifying the capacity of the affected cells to respond to later stimulatory events.

Vitamin D receptor and RXR in the post-genomic era

Following the elucidation of the human genome and components of the epigenome, it is timely to revisit what is known of vitamin D receptor (VDR) function. Early transcriptomic studies using microarray approaches focused on the protein coding mRNA that were regulated by the VDR, usually following treatment with ligand. These studies quickly established the approximate size and surprising diversity of the VDR transcriptome, revealing it to be highly heterogenous and cell type and time dependent. Investigators also considered VDR regulation of non-protein coding RNA and again, cell and time dependency was observed. Attempts to integrate mRNA and miRNA regulation patterns are beginning to reveal patterns of co-regulation and interaction that allow for greater control of mRNA expression, and the capacity to govern more complex cellular events. Alternative splicing in the trasncriptome has emerged as a critical process in transcriptional control and there is evidence of the VDR interacting with components of the splicesome. ChIP-Seq approaches have proved to be pivotal to reveal the diversity of the VDR binding choices across cell types and following treatment, and have revealed that the majority of these are non-canonical in nature. The underlying causes driving the diversity of VDR binding choices remain enigmatic. Finally, genetic variation has emerged as important to impact the transcription factor affinity towards genomic binding sites, and recently the impact of this on VDR function has begun to be considered.

Vitamin D and the RNA transcriptome: more than mRNA regulation

The GRCh37.p13 primary assembly of the human genome contains 20805 protein coding mRNA, and 37147 non-protein coding genes and pseudogenes that as a result of RNA processing and editing generate 196501 gene transcripts. Given the size and diversity of the human transcriptome, it is timely to revisit what is known of VDR function in the regulation and targeting of transcription. Early transcriptomic studies using microarray approaches focused on the protein coding mRNA that were regulated by the VDR, usually following treatment with ligand. These studies quickly established the approximate size, and surprising diversity of the VDR transcriptome, revealing it to be highly heterogenous and cell type and time dependent. With the discovery of microRNA, investigators also considered VDR regulation of these non-protein coding RNA. Again, cell and time dependency has emerged. Attempts to integrate mRNA and miRNA regulation patterns are beginning to reveal patterns of co-regulation and interaction that allow for greater control of mRNA expression, and the capacity to govern more complex cellular events. As the awareness of the diversity of non-coding RNA increases, it is increasingly likely it will be revealed that VDR actions are mediated through these molecules also. Key knowledge gaps remain over the VDR transcriptome. The causes for the cell and type dependent transcriptional heterogenetiy remain enigmatic. ChIP-Seq approaches have confirmed that VDR binding choices differ very significantly by cell type, but as yet the underlying causes distilling VDR binding choices are unclear. Similarly, it is clear that many of the VDR binding sites are non-canonical in nature but again the mechanisms underlying these interactions are unclear. Finally, although alternative splicing is clearly a very significant process in cellular transcriptional control, the lack of RNA-Seq data centered on VDR function are currently limiting the global assessment of the VDR transcriptome. VDR focused research that complements publically available data (e.g., ENCODE Birney et al., 2007; Birney, 2012), TCGA (Strausberg et al., 2002), GTEx (Consortium, 2013) will enable these questions to be addressed through large-scale data integration efforts.

Ribogenomics: the science and knowledge of RNA

Ribonucleic acid (RNA) deserves not only a dedicated field of biological research — a discipline or branch of knowledge — but also explicit definitions of its roles in cellular processes and molecular mechanisms. Ribogenomics is to study the biology of cellular RNAs, including their origin, biogenesis, structure and function. On the informational track, messenger RNAs (mRNAs) are the major component of ribogenomes, which encode proteins and serve as one of the four major components of the translation machinery and whose expression is regulated at multiple levels by other operational RNAs. On the operational track, there are several diverse types of RNAs — their length distribution is perhaps the most simplistic stratification — involving in major cellular activities, such as chromosomal structure and organization, DNA replication and repair, transcriptional/post-transcriptional regulation, RNA processing and routing, translation and cellular energy/metabolism regulation. An all-out effort exceeding the magnitude of the Human Genome Project is of essence to construct just mammalian transcriptomes in multiple contexts including embryonic development, circadian and seasonal rhythms, defined life-span stages, pathological conditions and anatomy-driven tissue/organ/cell types.