ErbB3 drives mammary epithelial survival and differentiation during pregnancy and lactation

Lactogenesis factors in the Asian population

The increasing incidence of disrupted lactogenesis in Asian populations underscores the importance of timely identification and efficient intervention. This study acknowledges the influence of ethnicity on genetic variations and aims to investigate the genetic mechanisms that contribute to lactogenesis in individuals of Asian descent. This study examines the possibilities of genetic screening as a means of applying preventive measures, with a particular focus on epigenetic techniques. Additionally, the analysis looks into the underlying mechanisms involved in milk production. This review discusses the intricate mechanisms underlying breast milk production and the potential influence of genetic variables on lactogenesis. Specifically, it explores the association between lactogenesis issues and genetic conditions such as depression, obesity, gestational hypertension, and gestational diabetes. These genetic factors could potentially be found by genetic screening as contributors to disruption in lactogenesis. This study aims to promote future investigation in the subject of genetic testing and its potential association with lactogenesis issues, despite the current scarcity of direct research on this topic. The statement posits that the progress made in genetic testing has the potential to provide novel insights into the timely identification and treatment of disrupted lactogenesis.

Molecular complexity of mammary glands development: a review of lactogenic differentiation in epithelial cells

The mammary gland is a dynamic organ with various physiological processes like cellular proliferation, differentiation, and apoptosis during the pregnancy-lactation-involution cycle. It is essential to understand the molecular changes during the lactogenic differentiation of mammary epithelial cells (MECs, the milk-synthesizing cells). The MECs are organized as luminal milk-secreting cells and basal myoepithelial cells (responsible for milk ejection by contraction) that form the alveoli. The branching morphogenesis and lactogenic differentiation of the MECs prepare the gland for lactation. This process is governed by many molecular mediators including hormones, growth factors, cytokines, miRNAs, regulatory proteins, etc. Interestingly, various signalling pathways guide lactation and understanding these molecular transitions from pregnancy to lactation will help researchers design further research. Manipulation of genes responsible for milk synthesis and secretion will promote augmentation of milk yield in dairy animals. Identifying protein signatures of lactation will help develop strategies for persistent lactation and shortening the dry period in farm animals. The present review article discusses in details the physiological and molecular changes occurring during lactogenic differentiation of MECs and the associated hormones, regulatory proteins, miRNAs, and signalling pathways. An in-depth knowledge of the molecular events will aid in developing engineered cellular models for studies related to mammary gland diseases of humans and animals.

ERBB Receptors and Their Ligands in the Developing Mammary Glands of Different Species: Fifteen Characters in Search of an Author

The ERBB tyrosine kinase receptors and their ligands belong to a complex family that has diverse biological effects and expression profiles in the developing mammary glands, where its members play an essential role in translating hormone signals into local effects. While our understanding of these processes stems mostly from mouse models, there is the potential for differences in how this family functions in the mammary glands of other species, particularly in light of their unique histomorphological features. Herein we review the postnatal distribution and function of ERBB receptors and their ligands in the mammary glands of rodents and humans, as well as for livestock and companion animals. Our analysis highlights the diverse biology for this family and its members across species, the regulation of their expression, and how their roles and functions might be modulated by varying stromal composition and hormone interactions. Given that ERBB receptors and their ligands have the potential to influence processes ranging from normal mammary development to diseased states such as cancer and/or mastitis, both in human and veterinary medicine, a more complete understanding of their biological functions should help to direct future research and the identification of new therapeutic targets.

Per- and Polyfluoroalkyl Substances and Breastfeeding as a Vulnerable Function: A Systematic Review of Epidemiological Studies

Milk formation in the breast during breastfeeding is a complex hormonally regulated process, potentially sensitive to the effects of endocrine-disrupting chemical exposures. The environmental chemicals, per- and polyfluoroalkyl substances (PFAS) are known endocrine disruptors. PFAS exposure have been associated with insufficient mammary gland development in mice and reduced breastfeeding duration in humans. The aim of this review was to gather the epidemiological evidence on the association between PFAS exposure and breastfeeding duration. Using PubMed and Embase, we performed a systematic literature search (on 23 January 2023) to identify epidemiological studies examining the association between maternal PFAS exposure and breastfeeding duration. Animal studies, reviews, and non-English studies were excluded. The risk of bias was assessed using the risk of bias in non-randomized studies of exposures tool. Estimates describing the association between PFAS exposure and the duration of breastfeeding were identified, and the data were synthesized separately for each type of PFAS and for the duration of exclusive and total breastfeeding. Six studies with between 336 and 2374 participants each were identified. PFAS exposure was assessed in serum samples (five studies) or based on residential address (one study). Five out of six studies found shorter total duration of breastfeeding with higher PFAS exposure. The most consistent associations were seen for perfluorooctane sulfonate (PFOS), perfluorooctanoic acid (PFOA), and perfluorononanoic acid (PFNA). The finding of a potential causal association between PFAS exposure and breastfeeding duration is in agreement with findings from experimental studies.

Neuronal-epithelial cross-talk drives acinar specification via NRG1-ERBB3-mTORC2 signaling

Acinar cells are the principal secretory units of multiple exocrine organs. A single-cell, layered, lumenized acinus forms from a large cohort of epithelial progenitors that must initiate and coordinate three cellular programs of acinar specification, namely, lineage progression, secretion, and polarization. Despite this well-known outcome, the mechanism(s) that regulate these complex programs are unknown. Here, we demonstrate that neuronal-epithelial cross-talk drives acinar specification through neuregulin (NRG1)-ERBB3-mTORC2 signaling. Using single-cell and global RNA sequencing of developing murine salivary glands, we identified NRG1-ERBB3 to precisely overlap with acinar specification during gland development. Genetic deletion of Erbb3 prevented cell lineage progression and the establishment of lumenized, secretory acini. Conversely, NRG1 treatment of isolated epithelia was sufficient to recapitulate the development of secretory acini. Mechanistically, we found that NRG1-ERBB3 regulates each developmental program through an mTORC2 signaling pathway. Thus, we reveal that a neuronal-epithelial (NRG1/ERBB3/mTORC2) mechanism orchestrates the creation of functional acini.

Defining mammary basal cell transcriptional states using single-cell RNA-sequencing

Breast cancer is a heterogenous disease that can be classified into multiple subtypes including the most aggressive basal-like and triple-negative subtypes. Understanding the heterogeneity within the normal mammary basal epithelial cells holds the key to inform us about basal-like cancer cell differentiation dynamics as well as potential cells of origin. Although it is known that the mammary basal compartment contains small pools of stem cells that fuel normal tissue morphogenesis and regeneration, a comprehensive yet focused analysis of the transcriptional makeup of the basal cells is lacking. We used single-cell RNA-sequencing and multiplexed RNA in-situ hybridization to characterize mammary basal cell heterogeneity. We used bioinformatic and computational pipelines to characterize the molecular features as well as predict differentiation dynamics and cell-cell communications of the newly identified basal cell states. We used genetic cell labeling to map the in vivo fates of cells in one of these states. We identified four major distinct transcriptional states within the mammary basal cells that exhibit gene expression signatures suggestive of different functional activity and metabolic preference. Our in vivo labeling and ex vivo organoid culture data suggest that one of these states, marked by Egr2 expression, represents a dynamic transcriptional state that all basal cells transit through during pubertal mammary morphogenesis. Our study provides a systematic approach to understanding the molecular heterogeneity of mammary basal cells and identifies previously unknown dynamics of basal cell transcriptional states.

Mammary epithelial cell transcriptome reveals potential roles of lncRNAs in regulating milk synthesis pathways in Jersey and Kashmiri cattle

Background

Long noncoding RNAs (lncRNAs) are now proven as essential regulatory elements, playing diverse roles in many biological processes including mammary gland development. However, little is known about their roles in the bovine lactation process.

Results

To identify and characterize the roles of lncRNAs in bovine lactation, high throughput RNA sequencing data from Jersey (high milk yield producer), and Kashmiri cattle (low milk yield producer) were utilized. Transcriptome data from three Kashmiri and three Jersey cattle throughout their lactation stages were utilized for differential expression analysis. At each stage (early, mid and late) three samples were taken from each breed. A total of 45 differentially expressed lncRNAs were identified between the three stages of lactation. The differentially expressed lncRNAs were found co-expressed with genes involved in the milk synthesis processes such as GPAM, LPL, and ABCG2 indicating their potential regulatory effects on milk quality genes. KEGG pathways analysis of potential cis and trans target genes of differentially expressed lncRNAs indicated that 27 and 48 pathways were significantly enriched between the three stages of lactation in Kashmiri and Jersey respectively, including mTOR signaling, PI3K-Akt signaling, and RAP1 signaling pathways. These pathways are known to play key roles in lactation biology and mammary gland development.

Conclusions

Expression profiles of lncRNAs across different lactation stages in Jersey and Kashmiri cattle provide a valuable resource for the study of the regulatory mechanisms involved in the lactation process as well as facilitate understanding of the role of lncRNAs in bovine lactation biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08406-x.

Screening and Conjoint Analysis of Key lncRNAs for Milk Fat Metabolism in Dairy Cows

Long noncoding RNAs (lncRNAs) play an important regulatory role in various biological processes as a key regulatory factor. However, the complete expression profile of lncRNAs in dairy cows and its function in milk fat synthesis are unknown. In this study, RNA sequencing (RNA-seq) was used to research the whole genome expression of lncRNAs and mRNA transcripts in high and low milk fat percentage (MFP) bovine mammary epithelial cells (BMECs), and joint analysis was carried out. We identified a total of 47 differentially expressed genes (DEGs) and 38 differentially expressed lncRNAs (DELs, Padj <0.05), enrichment analysis screened out 11 candidate DEGs that may regulate milk fat metabolism. Downregulated differential gene ENPP2 (The expression level in BMECs of high milk fat dairy cows was lower than that of low milk fat cows) and upregulated differential gene BCAT1 are more likely to participate in the milk fat metabolism, and its function needs further experiments verification. The enrichment analysis of target genes predicted by DELs identified 7 cis (co-localization) and 10 trans (co-expression) candidate target genes related to milk lipid metabolism, corresponding to a total of 18 DELs. Among them, the targeting relationship between long intervening/intergenic noncoding RNA (lincRNA) TCONS_00082721 and FABP4 is worthy of attention. One hundred and fifty-six competing endogenous RNAs (ceRNAs) interaction regulation networks related to milk fat metabolism were constructed based on the expression information of DELs, differential microRNAs (miRNAs), and lipid metabolism-related target genes. The regulatory network centered on miR-145 will be the focus of subsequent experimental research. The ceRNAs regulatory network related to TCONS_00082721 and TCONS_00172817 are more likely to be involved in milk fat synthesis. These results will provide new ways to understand the complex biology of dairy cow milk fat synthesis and provide valuable information for breed improvement of Chinese Holstein cow.

RNA-seq and nuclear proteomics provide insights into the lactation regulation mechanism of goat transfected IGF-I and GH recombinant vectors

There exists little available information on the mechanisms of lactation regulation until now. In order to explore the underlying mechanism, we injected IGF-I and GH recombinant vectors into the mammary gland, then RNA-seq analysis and nuclear proteomics were used for rapid high-throughput screening of DEGs and DEPs in the two groups linked to lactation regulation. KEGG analysis of 206 DEGs showed that the same 4 of top 10 enrichment pathways (ECM receptor interaction, protein digestion and absorption, focal adhesion and phagosome) involved in 4 co-expressed genes (IDO, BTG1, ITGB6 and keratin 83), the two groups enriched different metabolic pathways yet. Nuclear proteomics analysis showed 75 and 36 DEPs in the IGF-I and GH group respectively; Sixteen common proteins were identified between the IGF-I group and GH group, four of which (ALB, TPT1, CXXC-5 and ACTR2) significantly decreased and three of which (PRP1, PAG-9 and Hsp70) significantly increased. Similarly, DEPs in the two groups were enriched in same one of top 10 enrichment pathways (PI3K-Akt signaling pathway). Protein-protein interaction networks highlighted the contribution of glycosphingolipid biosynthesis, porphyrin and chlorophyll metabolism and the Jak-STAT signaling pathway to lactation regulation of GH and IGFI. GH and IGF-I improve milk yield, which may be linked to important nodal proteins (ALB and ACTB). Our research advances the understanding of the mammary gland transcriptome and nuclear proteomics during GH and IGF-I overexpression. Individual genes, proteins and pathways in this study point towards potential targets for lactation regulation.

Phytochemicals as PI3K/ Akt/ mTOR Inhibitors and Their Role in Breast Cancer Treatment

BACKGROUND

Breast cancer is the predominant form of cancer in women; various cellular pathways are involved in the initiation and progression of breast cancer. Among the various types of breast cancer that differ in their growth factor receptor status, PI3K/Akt signaling is a common pathway where all these converge. Thus, the PI3K signaling is of great interest as a target for breast cancer prevention; however, it is less explored.

OBJECTIVE

The present review is aimed to provide a concise outline of the role of PI3K/Akt/mTOR pathway in breast carcinogenesis and its progression events, including metastasis, drug resistance and stemness. The review emphasizes the role of natural and synthetic inhibitors of PI3K/Akt/m- TOR pathway in breast cancer prevention.

METHODS

The data were obtained from PubMed/Medline databases, Scopus and Google patent literature.

RESULTS

PI3K/Akt/mTOR signaling plays an important role in human breast carcinogenesis; it acts on the initiation and progression events associated with it. Numerous molecules have been isolated and identified as promising drug candidates by targeting the signaling pathway. Results from clinical studies confirm their application in the treatment of human breast cancer alone and in combination with classical chemotherapeutics as well as monoclonal antibodies.

CONCLUSION

PI3K/mTOR signaling blockers have evolved as promising anticancer agents by interfering breast cancer development and progression at various stages. Natural products and bioactive components are emerging as novel inhibitors of PI3K signaling and more research in this area may yield numerous drug candidates.

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.

Genome-wide scan reveals genomic regions and candidate genes underlying direct and maternal effects of preweaning calf mortality in Nellore cattle

We conducted analysis to estimate genetic parameters and to identify genomic regions and candidate genes affecting direct and maternal effects of preweaning calf mortality (PWM) in Nellore cattle. Phenotypic records of 67,196 animals, and 8443 genotypes for 410,936 SNPs were used. Analysis were performed through the weighted single-step GBLUP approach and considering a threshold animal model via Bayesian Inference. Direct and maternal heritability estimates were of 0.2143 ± 0.0348 and 0.0137 ± 0.0066, respectively. The top 10 genomic regions accounted for 13.61 and 14.23% of the direct and maternal additive genetic variances and harbored a total of 63 and 91 positional candidate genes, respectively. Two overlapping regions on BTA2 were identified for both direct and maternal effects. Candidate genes are involved in biological mechanisms i.e. embryogenesis, immune response, feto-maternal communication, circadian rhythm, hormone alterations, myometrium adaptation, and milk secretion, which are critical for the successful calf growth and survival during preweaning period.

Embryonic lethality and defective mammary gland development of activator-function impaired conditional knock-in Erbb3 V943R mice

ERBB3 is a pseudokinase domain-containing member of the ERBB family of receptor tyrosine kinases (RTKs). Following ligand binding, ERBB receptors homo- or hetero-dimerize, leading to a head-to-tail arrangement of the intracellular kinase domains, where the "receiver" kinase domain of one ERBB is activated by the "activator" domain of the other ERBB in the dimer. In ERBB3, a conserved valine at codon 943 (V943) in the kinase C-terminal domain has been shown to be important for its function as an "activator" kinase in vitro. Here we report a knock-in mouse model where we have modified the endogenous Erbb3 allele to allow for tissue-specific conditional expression of Erbb3 ^V943R (Erbb3 ^CKI-V943R ). Additionally, we generated an Erbb3 ^D850N (Erbb3 ^CKI-D850N ) conditional knock-in mouse model where the conserved aspartate in the DFG motif of the pseudokinase domain was mutated to abolish any potential residual kinase activity. While Erbb3 ^D850N/D850N animals developed normally, homozygous Erbb3 ^V943R/V943R expression during development resulted in embryonic lethality. Further, tissue specific expression of Erbb3 ^V943R/V943R in the mammary gland epithelium following its activation using MMTV-Cre resulted in delayed elongation of the ductal network during puberty. Single-cell RNA-seq analysis of Erbb3 ^V943R/V943R mammary glands showed a reduction in a specific subset of fibrinogen-producing luminal epithelial cells.

Identification and comparison of exosomal microRNAs in the milk and colostrum of two different cow breeds

Bovine milk and colostrum provide essential nutrients and immunologically active factors that are beneficial to a newborn calf. Milk-and-colostrum-derived exosomes are known as the most important for cellular communication. Exosomes also contain non-coding RNA, such as microRNA. However, there is limited information about exosomal miRNA derived from the milk and colostrum of Holstein and DAK cattle. This study aimed to identify and characterize the exosomal microRNA in the milk and colostrum of Holstein and Doğu Anadolu Kirmizisi (DAK) cows. For this purpose, total RNA isolation was carried out on the milk and colostrum samples that were collected from the Holstein and DAK cattle breeds. The RNA samples were subjected to RNA sequencing and the microRNAs were determined. Lastly, gene ontology analysis was performed for target genes. A total of 795 miRNAs that were expressed differently were identified. A total of 545 of these were known miRNAs and 260 were found to be novel miRNAs. In the functional enrichment analysis, the miRNAs expressed in Holstein milk were mostly associated with milk synthesis, and those in colostrum were mostly involved in the immunity pathways. It was also observed that the miRNAs expressed in DAK milk regulated milk fat and protein metabolism, and there were miRNAs that regulated immune pathways in the colostrum. In addition to this, many novel miRNAs were defined in DAK colostrum. When the target genes of exosomal miRNA in Holstein and DAK milk and colostrum were compared, it was suggested that the DAK breed had genes that were mostly associated with the immune system. As a result, the data obtained from this study will provide beneficial contributions to potential miRNA biomarker studies for milk yield and mastitis.

Kinase activity of ERBB3 contributes to intestinal organoids growth and intestinal tumorigenesis

As a member of the epidermal growth factor receptor (EGFR) family, ERBB3 plays an essential role in development and disease independent of inherently inactive kinase domain. Recently, ERBB3 has been found to bind to ATP and has catalytic activity in vitro. However, the biological function of ERBB3 kinase activity remains elusive in vivo. Here we have identified the physiological function of inactivated ERBB3 kinase activity by creating Erbb3‐K740M knockin mice in which ATP cannot bind to ERBB3. Unlike Erbb3 knockout mice, kinase‐inactive Erbb3^K740M homozygous mice were born in Mendelian ratios and showed normal development. After dextran sulfate sodium‐induced colitis, the kinase‐inactive Erbb3 mutant mice showed normal recovery. However, the outgrowth of ileal organoids by neuregulin‐1 treatment was more attenuated in Erbb3 mutant mice than in WT mice. Moreover, in combination with the Apc^Min mouse, the proportion of polyps less than 1 mm in diameter in mutant mice was higher than in control mice and an increase in the number of apoptotic cells was observed in polyps from mutant mice compared with polyps from control mice. Taken together, the ERBB3 kinase activity contributes to the outgrowth of ileal organoids and intestinal tumorigenesis, and the development of ERBB3 kinase inhibitors, including epidermal growth factor receptor family members, can be a potential way to target colorectal cancer.

Gas6 is dispensable for pubertal mammary gland development

Mammary gland development is a complex and dynamic process that occurs mainly postnatally. Ductal elongation and branching morphogenesis are regulated by a plethora of factors, including cytokines, hormones, growth factors and the extracellular matrix. Gas6 is a secreted gamma-carboxylated protein that binds to a family of receptors tyrosine kinase receptors known as the TAMR family (Tyro3, Axl, Mer). Gas6 function in developmental processes has been shown in nervous, reproductive and immune systems. In this study, we found that Gas6 is highly expressed in virgin adult mammary glands but declines during pregnancy and lactation. Specifically, Gas6 is highly expressed in luminal and basal mammary epithelial cells during puberty and adulthood, while TAMR expression is low. Mammary whole mount analysis revealed that Gas6 germline deletion does not impact ductal elongation, branching morphogenesis or terminal end bud formation. Masson's trichrome staining showed that collagen deposition is similar in Gas6-/- mice as compared to wildtype mice. Gas6-/- mammary glands presented an organized luminal and myoepithelial bilayer of cells, and the proportion of mammary stem cells was unchanged in Gas6-/- mammary glands as compared to wildtype. Finally, proliferation of epithelial cells and macrophage number were similar in both groups. These studies suggest that Gas6 is not essential for pubertal mammary gland development in nulliparous mice.

Chi-miR-3031 regulates beta-casein via the PI3K/AKT-mTOR signaling pathway in goat mammary epithelial cells (GMECs)

BACKGROUND

MicroRNAs can regulate gene expression at the posttranscriptional level through translational repression or target degradation. Our previous investigations examined the differential expression levels of chi-miR-3031 in caprine mammary gland tissues in colostrum and common milk stages.

RESULTS

The present study detected the role of chi-miR-3031 in the lactation mechanisms of GMECs. High-throughput sequencing was used to analyze transcriptomic landscapes of GMECs transfected with chi-miR-3031 mimics (MC) and a mimic negative control (NC). In the MC and NC groups, we acquired 39,793,503 and 36,531,517 uniquely mapped reads, respectively, accounting for 85.85 and 81.66% of total reads. In the MC group, 180 differentially expressed unigenes were downregulated, whereas 157 unigenes were upregulated. KEGG pathway analyses showed that the prolactin, TNF and ErbB signaling pathways, including TGFα, PIK3R3, IGF2, ELF5, IGFBP5 and LHβ genes, played important roles in mammary development and milk secretion. Results from transcriptome sequencing, real-time PCR and western blotting showed that chi-miR-3031 suppressed the expression of IGFBP5 mRNA and protein. The expression levels of β-casein significantly increased in the MC and siRNA-IGFBP5 groups. We observed that the down-regulation of IGFBP5 activated mTOR at the Ser2448 site in GMECs transfected with MC and siRNA-IGFBP5. Previous findings and our results showed that chi-miR-3031 activated the PI3K-AKT-mTOR pathway and increased β-casein expression by down-regulating IGFBP5.

CONCLUSIONS

These findings will afford valuable information for improving milk quality and contribute the development of potential methods for amending lactation performance.

Transcriptome sequencing to detect the potential role of long non-coding RNAs in bovine mammary gland during the dry and lactation period

BACKGROUND

It is known that long non-coding RNAs (lncRNAs) play an important role in various biological processes, including cell proliferation, differentiation and apoptosis. However, their functions and profiles in lactation cycle of dairy cows are largely unknown. In this study, lncRNA-seq technique was employed to compare the expression profiles of lncRNAs and mRNAs from Chinese Holstein mammary gland in dry and lactation period.

RESULT

Totally 3746 differentially expressed lncRNAs (DELs) and 2890 differentially expressed genes (DEGs) were identified from the dry and lactation mammary glands of Holstein cows. Functional enrichment analysis on target genes of lncRNAs indicated that these genes were involved in lactation-related signaling pathways, including cell cycle, JAK-STAT, cell adhesion, and PI3K-Akt signaling pathways. Additionally, the interaction between lncRNAs and their potential miRNAs was predicted and partly verified. The result indicated that the lactation-associated miR-221 might interact with lncRNAs TCONS_00040268, TCONS_00137654, TCONS_00071659 and TCONS_00000352, which revealed that these lncRNAs might be important regulators for lactation cycle.

CONCLUSION

This study provides a resource for lncRNA research on lactation cycle of bovine mammary gland. Besides, the interaction between lncRNAs and the specific miRNA is revealed. It expands our knowledge about lncRNA and miRNA biology as well as contributes to clarify the regulation of lactation cycle of bovine mammary gland.