Transcriptome analysis revealed the characteristics and functions of long non-coding RNAs in the hypothalamus during sexual maturation in goats.

The hypothalamus is an essential neuroendocrine area in animals that regulates sexual development. Long non-coding RNAs (lncRNAs) are hypothesized to regulate physiological processes related to animal reproduction. However, the regulatory mechanism by which lncRNAs participate in sexual maturity in goats is poorly known, particularly from birth to sexual maturation. In this study, RNAseq analysis was conducted on the hypothalamus of four developmental stages (1day (D1, n = 5), 2 months (M2, n = 5), 4 months (M4, n = 5), and 6 months (M6, n = 5)) of Jining grey goats. The results showed that a total of 237 differentially expressed lncRNAs (DELs) were identified in the hypothalamus. Among these, 221 DELs exhibited cis-regulatory effects on 693 target genes, while 24 DELs demonstrated trans-regulatory effects on 63 target genes. The target genes of these DELs are mainly involved in biological processes related to energy metabolism, signal transduction and hormone secretion, such as sphingolipid signaling pathway, adipocytokine signaling pathway, neurotrophic signaling pathway, glutamatergic synapse, P53 signaling pathway and GnRH signaling pathway. In addition, XR_001918477.1, TCONS_00077463, XR_001918760.1, and TCONS_00029048 and their potential target genes may play a crucial role in the process of goat sexual maturation. This study advances our understanding of lncRNA in hypothalamic tissue during sexual maturation in goats and will give a theoretical foundation for improving goat reproductive features.

Comparative metabolomics reveals serum metabolites changes in goats during different developmental stages.

Goats can provide meat, milk and skins for humans and are livestock with high economic benefits. However, despite their economic significance, the comprehensive analysis of goats' serum metabolic profile and its intricate alterations throughout their developmental journey remains conspicuously absent. To investigate the stage-specificity and dynamic change characteristics of metabolites during the growth and development of goats, this study compared the alterations in serum hormone levels and serum biochemical markers across different developmental stages of female goats (1, 60, 120 and 180 days old; n = 5). Additionally, a serum untargeted LC-MS metabolomics analysis was conducted. A total of 504 DAMs were identified with age. The results indicated that PE, PC, Lyso-PE, Lyso-PC and FAFHA may play important roles in lipid metabolism in goats after birth. Weighted gene co-expression network analysis (WGCNA) identified two metabolite modules (Turquoise and Yellow) and key metabolites within these modules that were significantly associated with phenotypic features. L-carnitine may be a metabolite related to muscle development in goats. The findings of this study demonstrate notable variations in serum metabolites across distinct developmental phases in goats. Lipids and organic acids play important roles in different developmental stages of goats.

Characterization of MicroRNA expression profiles in the ovarian tissue of goats during the sexual maturity period.

BACKGROUND: The ovary is an important reproductive organ in mammals, and its development directly affects the sexual maturity and reproductive capacity of individuals. MicroRNAs (miRNAs) are recognized as regulators of reproductive physiological processes in various animals and have been shown to regulate ovarian development through typical targeting and translational repression. However, little is known about the regulatory role of miRNAs in ovarian tissue development during sexual maturity in goats. To comprehensively profile the different physiological stages of sexual maturation in goats, we performed small-RNA sequencing of ovarian tissue samples collected at four specific time points (1 day after birth (D1), 2 months old (M2), 4 months old (M4), and 6 months old (M6)). In addition, we used ELISAs to measure serum levels of reproductive hormones to study their temporal changes. RESULTS: The results showed that serum levels of gonadotropin-releasing hormone, follicle-stimulating hormone, luteinizing hormone, oestradiol, progesterone, oxytocin, and prolactin were lower in goats at the D1 stage than at the other three developmental stages (P < 0.05). The secretion patterns of these seven hormones show a similar trend, with hormone levels reaching their peaks at 4 months of age. A total of 667 miRNAs were detected in 20 libraries, and 254 differentially expressed miRNAs and 3 groups of miRNA clusters that had unique expression patterns were identified (|log2-fold change|> 1, FDR < 0.05) in the 6 comparison groups. RT‒qPCR was employed to confirm that the expression pattern of the 15 selected miRNAs was consistent with the Illumina sequencing results. Gene ontology analyses revealed significant enrichment of GO terms such as cell proliferation regulation, epithelial cell development, and amino acid transport, as well as important signaling pathways including the MAPK signaling pathway, the PI3K-Akt signaling pathway, and the oestrogen signaling pathway. Further miRNA‒mRNA regulation network analysis revealed that 8 differentially expressed miRNAs (chi-miR-1343, chi-miR-328-3p, chi-miR-877-3p, chi-miR-296-3p, chi-miR-128-5p, chi-miR-331-3p, chi-miR-342-5p and chi-miR-34a) have important regulatory roles in ovarian cell proliferation, hormone secretion and metabolism-related biological processes. CONCLUSIONS: Overall, our study investigated the changes in serum hormone and miRNA levels in the ovaries. These data provide a valuable resource for understanding the molecular regulatory mechanisms of miRNAs in ovarian tissue during the sexual maturity period in goats.

Identification and Characterization of circRNAs in Non-Lactating Dairy Goat Mammary Glands Reveal Their Regulatory Role in Mammary Cell Involution and Remodeling.

This study conducted transcriptome sequencing of goat-mammary-gland tissue at the late lactation (LL), dry period (DP), and late gestation (LG) stages to reveal the expression characteristics and molecular functions of circRNAs during mammary involution. A total of 11,756 circRNAs were identified in this study, of which 2528 circRNAs were expressed in all three stages. The number of exonic circRNAs was the largest, and the least identified circRNAs were antisense circRNAs. circRNA source gene analysis found that 9282 circRNAs were derived from 3889 genes, and 127 circRNAs' source genes were unknown. Gene Ontology (GO) terms, such as histone modification, regulation of GTPase activity, and establishment or maintenance of cell polarity, were significantly enriched (FDR < 0.05), which indicates the functional diversity of circRNAs' source genes. A total of 218 differentially expressed circRNAs were identified during the non-lactation period. The number of specifically expressed circRNAs was the highest in the DP and the lowest in LL stages. These indicated temporal specificity of circRNA expression in mammary gland tissues at different developmental stages. In addition, this study also constructed circRNA-miRNA-mRNA competitive endogenous RNA (ceRNA) regulatory networks related to mammary development, immunity, substance metabolism, and apoptosis. These findings help understand the regulatory role of circRNAs in mammary cell involution and remodeling.

Analyses of lncRNAs, circRNAs, and the Interactions between ncRNAs and mRNAs in Goat Submandibular Glands Reveal Their Potential Function in Immune Regulation.

As part of one of the main ruminants, goat salivary glands hardly secrete digestive enzymes, but play an important role in immunity. The immune function of goat salivary glands significantly changes with age, while the expression profile and specific function of non-coding RNA during this process are unknown. In this study, transcriptome sequencing was performed on submandibular gland (SMG) tissues of 1-month-old, 12-month-old, and 24-month-old goats, revealing the expression patterns of lncRNA and circRNA at different ages. A total of 369 lncRNAs and 1699 circRNAs were found to be differentially expressed. Functional enrichment analyses showed that the lncRNA regulated target mRNAs and circRNA host genes were significantly enriched in immune-related GO terms and pathways. CeRNA network analysis showed that the key differentially expressed circRNAs and lncRNAs mainly regulate the key immune-related genes ITGB2, LCP2, PTPRC, SYK, and ZAP70 through competitive binding with miR-141-x, miR-29-y, and chi-miR-29b-3p, thereby affecting the natural killer cell-mediated cytotoxicity pathway, the T cell receptor signaling pathway, and other immune-related pathways. It should be noted that the expression of key circRNAs, lncRNAs, and key immune-related genes in goat SMGs decreased significantly with the growth of the goat. This is the first reporting of lncRNAs, circRNAs, and ceRNA network regulation in goat SMGs. Our study contributes to the knowledge of changes in the expression of non-coding RNAs during SMG development in goats and provides new insights into the relationship between non-coding RNAs and salivary gland immune function in goats.

Identification and functional analysis of m6A in the mammary gland tissues of dairy goats at the early and peak lactation stages.

N6-methyladenosine (m6A) is the most common reversible epigenetic RNA modification in the mRNA of all higher eukaryotic organisms and plays an important role in the regulation of gene expression and cell function. In this study, m6A-modified methylated RNA immunoprecipitation sequencing (MeRIP-seq) and transcriptome sequencing (RNA-seq) were used to identify the key genes with m6A modification during mammary gland development and lactation in dairy goats. The results showed that m6A methylation occurred at 3,927 loci, which were significantly enriched in the 3' untranslated region (3'UTR) and the termination codon region. In the early stage and peak stage of lactation, m6A methylation occurred extensively in mammary tissues, and a total of 725 differentially expressed m6A-modified genes were obtained, all negatively correlated with mRNA expression. In addition, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that different methylated genes were mainly involved in the growth and apoptosis of mammary epithelial cells through signaling pathways, such as the mitogen-activated protein kinase (MAPK) and phospholipase D pathways, and then affected the development and lactation of mammary gland. All in all, we identified and analyzed the methylation events related to the development and lactation regulation of mammary gland at the early and peak lactation stages, and provided a theoretical basis to reveal the physiological regulatory system of mammary gland development and lactation in dairy goats.

Transcriptome Analysis of Goat Mammary Gland Tissue Reveals the Adaptive Strategies and Molecular Mechanisms of Lactation and Involution.

To understand how genes precisely regulate lactation physiological activity and the molecular genetic mechanisms underlying mammary gland involution, this study investigated the transcriptome characteristics of goat mammary gland tissues at the late gestation (LG), early lactation (EL), peak lactation (PL), late lactation (LL), dry period (DP), and involution (IN) stages. A total of 13,083 differentially expressed transcripts were identified by mutual comparison of mammary gland tissues at six developmental stages. Genes related to cell growth, apoptosis, immunity, nutrient transport, synthesis, and metabolism make adaptive transcriptional changes to meet the needs of mammary lactation. Notably, platelet derived growth factor receptor beta (PDGFRB) was screened as a hub gene of the mammary gland developmental network, which is highly expressed during the DP and IN. Overexpression of PDGFRB in vitro could slow down the G1/S phase arrest of goat mammary epithelial cell cycle and promote cell proliferation by regulating the PI3K/Akt signaling pathway. In addition, PDGFRB overexpression can also affect the expression of genes related to apoptosis, matrix metalloproteinase family, and vascular development, which is beneficial to the remodeling of mammary gland tissue during involution. These findings provide new insights into the molecular mechanisms involved in lactation and mammary gland involution.

Characterization of long noncoding RNA in nonlactating goat mammary glands reveals their regulatory role in mammary cell involution and remodeling.

Long noncoding RNA (lncRNA) can regulate mammary gland development and lactation physiological activities. However, the molecular genetic mechanisms of lncRNA in mammary gland involution and cell remodeling remain unclear. This work analyzed the expression characteristics and molecular functions of lncRNA in goat mammary gland tissue at the late lactation (LL), dry period (DP), and late gestation (LG) stages. Sequencing results showed that 3074 lncRNAs were identified in non-lactating goat mammary gland tissue. Statistical analysis of lncRNA length characteristics and exon number found that goat lncRNAs were shorter in length, had fewer exons, and significantly lower expression levels than those of protein-coding genes. 331 differentially expressed lncRNAs were identified in the three comparison groups (LLvsDP, DPvsLG, and LLvsLG), which indicated that the lncRNAs expression at the transcriptional level were changed during mammary involution. Interestingly, lncRNAs were more actively expressed during the dry period compared to lactation, suggesting that lncRNAs in mammary glands are developmentally specific. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses revealed that lncRNAs could regulate immune function, cell proliferation, apoptosis, hormones, substance metabolism, transport, and intercellular communication in the mammary gland through various action modes. Among them, cis-acting lncRNAs enhanced the protection of mammary gland health during the dry period and late gestation. The above reflects the particular mechanisms of lncRNA to adapt to the developmental needs of mammary involution and remodeling. Furthermore, in the lncRNA-miRNA-mRNA network associated with mammary gland development, the expression of LOC102168552 was higher in late gestation than in the dry period and late lactation. Its expression was positively correlated with PRLR and negatively correlated with chi-miR-324-3p. Overexpression of LOC102168552 in goat mammary epithelial cells cultured in vitro could up-regulate PRLR to activate the prolactin signaling pathway by competitively binding to chi-miR-324-3p, promoting cell proliferation, reducing cell cycle arrest in the G1 / S phase, and inhibiting apoptosis. However, overexpression of LOC102168552 alone did not affect mammary cell growth status and the prolactin signaling pathway. This indicates that LOC102168552 must rely on chi-miR-324-3p to inhibit mammary cell apoptosis. In conclusion, the above analysis revealed that lncRNAs in goat mammary tissue are differentially expressed at different stages of involution. As expected, lncRNAs adaptively regulate various physiological activities during mammary gland involution through multiple modes of action, in preparation for a new round of lactation. These findings provide a reference and help further understand the regulatory role of lncRNAs in mammary cell involution and remodeling.

High-Throughput Sequencing Reveals Transcriptome Signature of Early Liver Development in Goat Kids.

As a vital metabolic and immune organ in animals, the liver plays an important role in protein synthesis, detoxification, metabolism, and immune defense. The primary research purpose of this study was to reveal the effect of breast-feeding, weaning transition, and weaning on the gene expression profile in the goat kid liver and to elucidate the transcriptome-level signatures associated with liver metabolic adaptation. Therefore, transcriptome sequencing was performed on liver tissues, which was collected at 1 day (D1), 2 weeks (W2), 4 weeks (W4), 8 weeks (W8), and 12 weeks (W12) after birth in Laiwu black goats at five different time-points, with five goats at each time point. From 25 libraries, a total of 37497 mRNAs were found to be expressed in goat kid livers, and 1271 genes were differentially expressed between at least two of the five time points. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that these genes were annotated as an extracellular region fraction, exhibiting monooxygenase activity, positive regulation of T cell activation, mitotic spindle mid-region assembly, cytokinesis, cytoskeleton-dependent cytokinesis, regulation of cytokinesis, regulation of lymphocyte proliferation, and so on. In addition, it mainly deals with metabolism, endocrine, cell proliferation and apoptosis, and immune processes. Finally, a gene regulatory network was constructed, and a total of 14 key genes were screened, which were mainly enriched for cell growth and development, endocrine, immune, and signal transduction-related pathways. Our results provide new information on the molecular mechanisms and pathways involved in liver development, metabolism, and immunity of goats.

Transcriptome profiling of the nonlactating mammary glands of dairy goats reveals the molecular genetic mechanism of mammary cell remodeling.

The mammary gland redevelops to the prepregnancy state during involution, which shows the mammary cells have the characteristics of remodeling. The rapidity and degree of mammary gland involution vary across species (e.g., between model organism mice and dairy livestock). However, the molecular genetic mechanism of involution and remodeling of goat mammary gland has not yet been clarified. This work investigated the structural changes and transcriptome characteristics of the mammary gland tissue of nonlactating dairy goats during the late lactation (LL), the dry period (DP), and late gestation (LG). Terminal deoxynucleotidyl transferase dUTP nick-end labeling (TUNEL) staining revealed significant changes in the structure of the nonlactating goat mammary gland, and obvious cell apoptosis occurred at LL and DP. Sequencing identified 1,381 genes that are differentially expressed in mammary gland tissue at the 3 developmental stages. Genes related to cell growth, apoptosis, immunity, nutrient transport, synthesis, and metabolism exhibited adaptive transcriptional changes to meet the needs of a new set of mammary gland lactation functions. The significant enrichment of Gene Ontology terms such as humoral immune response, complement activation, and neutrophil-mediated immunity indicates that the innate immune system plays an important role in maintaining the health of degenerative mammary glands and eliminating apoptotic cells. The peroxisome proliferator-activated receptor signaling pathway plays an important regulatory role in lipid metabolism, especially the adaptive changes in expression of genes encoded lipid transport and enzymes, which promote the formation of milk fat during the lactation. The mammary gland development gene module revealed that pregnancy hormone receptors, cell growth factors and their receptors, and genes encoding insulin-like growth factor binding proteins regulate the physiological process of mammary gland involution through adaptive transcriptional changes. Interestingly, ERBB4 was identified as the hub gene of the network that regulates mammary gland growth and development. Overexpression of ERBB4 in mammary epithelial cells cultured in vitro can reduce cell cycle arrest in G1/S phase and apoptosis by regulating the PI3K/Akt signaling pathway and promote the proliferation of mammary epithelial cells. The gene ERBB4 also affects the expression of genes that initiate mammary gland involution and promote mammary gland remodeling. These findings contribute to an in-depth understanding of the molecular mechanisms involved in mammary gland involution and remodeling.

Differentially Expressed MiRNAs of Goat Submandibular Glands Among Three Developmental Stages Are Involved in Immune Functions.

Submandibular glands (SMGs) are one of the primary components of salivary glands in goats. The proteins and biologically active substances secreted by the SMGs change with growth and development. Our previous studies showed that most of the differentially expressed genes in the SMGs of goats at different developmental stages are involved in immune-related signaling pathways, but the miRNA expression patterns in the same tissues are unknown. The aim of this study was to reveal the expression profile of miRNAs at three different developmental stages, detect differentially expressed miRNAs (DE miRNAs) and predict disease-related DE miRNAs. SMG tissue samples were collected from groups of 1-month-old kids, 12-month-old maiden goats and 24-month-old adult goats (three samples from each group), and high-throughout transcriptome sequencing was conducted. A total of 178, 241 and 7 DE miRNAs were discovered between 1-month-old kids and 12-month-old maiden goats, between 1-month-old kids and 24-month-old adult goats, and between 12-month-old maiden goats and 24-month-old adult goats, respectively. Among these DE miRNAs, 88 DE miRNAs with medium or high expression levels (TPM ≥50) were classified into five expression pattern clusters. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that some of the predicted target genes of the DE miRNAs in the five clusters were enriched in disease-related GO terms and pathways. MiRNA target genes in significant pathways were significantly enriched in Hepatitis B (FDR = 9.03E-10) and Pathways in cancer (FDR = 4.2E-10). Further analysis was performed with a PPI network, and 10 miRNAs were predicted to play an important role in the occurrence and prevention of diseases during the growth and development of goats.

Characterization of the microRNA Expression Profiles in the Goat Kid Liver.

The liver is the largest digestive gland in goats with an important role in early metabolic function development. MicroRNAs (miRNA) are crucial for regulating the development and metabolism in the goat liver. In the study, we sequenced the miRNAs in the liver tissues of the goat kid to further research their regulation roles in early liver development. The liver tissues were procured at 5-time points from the Laiwu black goats of 1 day (D1), 2 weeks (W2), 4 weeks (W4), 8 weeks (W8), and 12 weeks (W12) after birth, respectively with five goats per time point, for a total of 25 goats. Our study identified 214 differential expression miRNAs, and the expression patterns of 15 randomly selected miRNAs were examined among all five age groups. The Gene ontology annotation results showed that differential expression miRNA (DE miRNA) target genes were significantly enriched in the fatty acid synthase activity, toxin metabolic process, cell surface, and antibiotic metabolic process. The KEGG analysis result was significantly enriched in steroid hormone synthesis and retinol metabolism pathways. Further miRNA-mRNA regulation network analysis reveals 9 differently expressed miRNA with important regulation roles. Overall, the DE miRNAs were mainly involved in liver development, lipid metabolism, toxin related metabolism-related biological process, and pathways. Our results provide new information about the molecular mechanisms and pathways in the goat kid liver development.

Transcriptome analysis reveals potential immune function-related regulatory genes/pathways of female Lubo goat submandibular glands at different developmental stages.

BACKGROUND: The submandibular glands, as major salivary glands, participate in rumen digestion in goats. Sialic acid, lysozyme, immunoglobulin A (IgA), lactoferrin and other biologically active substances secreted in the submandibular glands were reported in succession, which suggests that the submandibular gland may have immune functions in addition to participating in digestion. The aim of this study was to map the expression profile of differentially expressed genes (DEGs) at three different stages by transcriptome sequencing, screen immune-related genes and pathways by bioinformatics methods, and predict the immune function of submandibular glands at different developmental stages. METHODS: Nine submandibular gland tissue samples were collected from groups of 1-month-old kids, 12-month-old adolescent goats and 24-month-old adult goats (3 samples from each group), and high-throughput transcriptome sequencing was conducted on these samples. The DEGs among the three stages were screened and analysed. Key genes and signalling pathways were selected via protein-protein interaction (PPI) network analysis. RESULTS: The results revealed 2,706, 2,525 and 52 DEGs between 1-month-old and 12-month-old goats, between 1-month-old and 24-month-old goats, and between 12-month-old and 24-month-old goats, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses indicated that most of the DEGs were enriched in immune- related GO terms and pathways. Based on functional enrichment analysis and network analysis, 10 genes (PTPRC, CD28, SELL, LCP2, MYC, LCK, ZAP70, ITGB2, SYK and CCR7), two signalling pathways (the T cell receptor signalling pathway and the NF-κß signalling pathway) and eight GO terms (T cell receptor signalling pathway, neutrophil mediated immunity, B cell mediated immunity, regulation of alpha-beta T cell activation, positive regulation of T cell proliferation, regulation of leukocyte differentiation, positive regulation of antigen receptor-mediated signalling pathway, positive regulation of lymphocyte proliferation) that may play key roles in the immune functions of the goat submandibular glands at different developmental stages were identified. Moreover, we found that eight antibacterial peptide-encoding genes were downregulated in the tuberculosis and salivary secretion pathways, while all immunoglobulins were upregulated in 10 immune system pathways. These findings indicate that the submandibular glands may be important immunological organs during the growth process of goats and that the immune function of these glands gradually weakens with age up to 12 months but remains relatively stable after 12 months of age. Overall, this study will improve our understanding of transcriptional regulation related to goat submandibular gland immune function.

Characterization of microRNA profiles in the mammary gland tissue of dairy goats at the late lactation, dry period and late gestation stages.

MicroRNAs (miRNAs) play an important role in regulating mammary gland development and lactation. We previously analyzed miRNA expression profiles in Laoshan dairy goat mammary glands at the early (20 d postpartum), peak (90 d postpartum) and late lactation (210 d postpartum) stages. To further enrich and clarify the miRNA expression profiles during the lactation physiological cycle, we sequenced miRNAs in the mammary gland tissues of Laoshan dairy goats at three newly selected stages: the late lactation (240 d postpartum), dry period (300 d postpartum) and late gestation (140 d after mating) stages. We obtained 4038 miRNAs and 385 important miRNA families, including mir-10, let-7 and mir-9. We also identified 754 differentially expressed miRNAs in the mammary gland tissue at the 3 different stages and 6 groups of miRNA clusters that had unique expression patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that GO terms such as mammary gland development (GO:0030879) and mammary gland morphogenesis (GO:0060443) and important signaling pathways, including the insulin signaling pathway (chx04910), hippo signaling pathway (chx04390) and estrogen signaling pathway (chx04915), were enriched. We screened miRNAs and potential target genes that may be involved in the regulation of lactation, mammary gland growth and differentiation, cell apoptosis, and substance transport and synthesis and detected the expression patterns of important genes at the three stages. These miRNAs and critical target genes may be important factors for mammary gland development and lactation regulation and potentially valuable molecular markers, which may provide a theoretical reference for further investigation of mammary gland physiology.

Genome­wide integrated analysis demonstrates widespread functions of lncRNAs in mammary gland development and lactation in dairy goats.

BACKGROUND: The mammary gland is a unique organ for milk synthesis, secretion and storage, and it undergoes cyclical processes of development, differentiation, lactation and degeneration. At different developmental periods, the biological processes governing mammary gland physiology and internal environmental homeostasis depend on a complex network of genes and regulatory factors. Emerging evidence indicates that lncRNAs have arbitrarily critical functions in regulating gene expression in many organisms; however, the systematic characteristics, expression, and regulatory roles of lncRNAs in the mammary gland tissues of dairy goats have not been determined. RESULT: In the present study, we profiled long noncoding RNA (lncRNA) expression in the mammary gland tissues of Laoshan dairy goats (Capra hircus) from different lactation periods at the whole-genome level, to identify, characterize and explore the regulatory functions of lncRNAs. A total of 37,249 transcripts were obtained, of which 2381 lncRNAs and 37,249 mRNAs were identified, 22,488 transcripts, including 800 noncoding transcripts and 21,688 coding transcripts, differed significantly (p ≤ 0.01) among the different lactation stages. The results of lncRNA-RNA interaction analysis showed that six known lncRNAs belonging to four families were identified as the precursors of 67 known microRNAs; 1478 and 573 mRNAs were predicted as hypothetical cis-regulation elements and antisense mRNAs, respectively. GO annotation and KEGG analysis indicated that the coexpressed mRNAs were largely enriched in biological processes related to such activities as metabolism, immune activation, and stress,., and most genes were involved in pathways related to such phenomena as inflammation, cancer, signal transduction, and metabolism. CONCLUSIONS: Our results clearly indicated that lncRNAs involved in responses to stimuli, multiorganism processes, development, reproductive processes and growth, are closely related to mammary gland development and lactation.

chi-miR-143-3p Promotes Apoptosis of Mammary Gland Epithelial Cells from Dairy Goats by Targeting Ndfip1.

The mammary gland is an important organ for lactation in dairy goats. Mammary gland development and lactation functions are primarily regulated by natural hormones and certain crucial regulatory factors. Nedd4 family-interacting protein 1 (Ndfip1) can specifically bind to neural precursor cell-expressed, developmentally downregulated protein 4 (Nedd4) family members to participate in ubiquitination, which in turn regulates a range of biological processes in the body. However, the effects of Ndfip1 expression regulation at the post-transcriptional level on the development of mammary gland cells have not been previously reported. To study the regulation of Ndfip1 at post-transcriptional level, the overexpression and interference vectors of Ndfip1 were constructed, and co-transfected into the primary mammary gland epithelial cells cultured in vitro with miR-143 mimics and inhibitor. Dual luciferase reporter gene system, real-time quantitative polymerase chain reaction, western blotting, cholecystokinin octapeptide assays, and flow cytometry were used to identify their regulation and function. As a result, Ndfip1 was targeted and regulated by miR-143, which influences the development of mammary gland epithelial cells in dairy goats cultured in vitro. This study will lay an experimental foundation for further understanding the functions of Ndfip1 and miR-143.

Transcriptome Analysis of Dairy Goat Mammary Gland Tissues from Different Lactation Stages.

The mammary gland is a unique organ involved in lactation in dairy livestock, and its development and lactation ability are affected by both genetic and environmental factors. To explore the molecular regulatory mechanisms of these factors, this study used high-throughput sequencing technology and bioinformatics methods to systematically analyze the transcriptome of Laoshan dairy goat mammary gland tissues from different lactation stages. From three libraries, 36,336,892, 36,469,596, and 35,759,380, pure sequences were obtained, with 25,292, 23,665 and 27,220 expressed genes, respectively, resulting in a total of 14,892 nonredundant differentially expressed genes. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that these genes were annotated to biological regulation, cellular processes, metabolic processes, cells, organelles, binding, catalytic activity and transcriptional activity, and so on. The genes were mainly involved in cellular processes, environmental information processes, genetic information processes, human diseases, metabolism, and organ systems. Finally, 20 genes related to mammary gland development and lactation were screened to construct a gene regulatory network. These findings support the involvement of the constructed network in the regulation of mammary gland development and lactation, and they not only lay the foundation to further investigate and screen the main genes or molecular genetic markers controlling the development of the mammary gland and lactation ability of dairy goats but also increase the in-depth understanding of the lactation physiology of dairy goats.

Sheep skeletal muscle transcriptome analysis reveals muscle growth regulatory lncRNAs.

As widely distributed domestic animals, sheep are an important species and the source of mutton. In this study, we aimed to evaluate the regulatory lncRNAs associated with muscle growth and development between high production mutton sheep (Dorper sheep and Qianhua Mutton Merino sheep) and low production mutton sheep (Small-tailed Han sheep). In total, 39 lncRNAs were found to be differentially expressed. Using co-expression analysis and functional annotation, 1,206 co-expression interactions were found between 32 lncRNAs and 369 genes, and 29 of these lncRNAs were found to be associated with muscle development, metabolism, cell proliferation and apoptosis. lncRNA-mRNA interactions revealed 6 lncRNAs as hub lncRNAs. Moreover, three lncRNAs and their associated co-expressed genes were demonstrated by cis-regulatory gene analyses, and we also found a potential regulatory relationship between the pseudogene lncRNA LOC101121401 and its parent gene FTH1. This study provides a genome-wide resolution of lncRNA and mRNA regulation in muscles from mutton sheep.

Identification and characterization of microRNAs in the intestinal tissues of sheep (Ovis aries).

Sheep are small ruminants, and their long intestines exhibit high digestive and absorptive capacity in many different rearing conditions; however, the genetic bases of this characteristic remains unclear. MicroRNAs (miRNAs) play a major role in maintaining both intestinal morphological structure as well as in regulating the physiological functions of this organ. However, no study has reported on the miRNA expression profile in the intestinal tissues of sheep. Here, we analyzed and identified the miRNA expression profile of three different intestinal tissues (i.e., duodenum, cecum, and colon) of sheep (Ovis aries) using high-throughput sequencing and bioinformatic methods. In total, 106 known miRNAs were identified, 458 conserved miRNAs were detected, 192 unannotated novel miRNAs were predicted, and 195 differentially expressed miRNAs were found between the different tissues. Additionally, 3,437 candidate target genes were predicted, and 17 non-redundant significantly enriched GO terms were identified using enrichment analysis. A total of 99 candidate target genes were found to significantly enriched in 4 KEGG biological pathways. A combined regulatory network was constructed based on 92 metabolism-related candidate target genes and 65 differentially expressed miRNAs, among which 7 miRNAs were identified as hub miRNAs. Via these mechanisms, miRNAs may play a role in maintaining intestinal homeostasis and metabolism. This study helps to further explain the mechanisms that underlie differences in tissue morphology and function in three intestinal segments of sheep.

Screening of miRNA profiles and construction of regulation networks in early and late lactation of dairy goat mammary glands.

In recent years, studies related to the expression profiles of miRNAs in the dairy goat mammary gland were performed, but regulatory mechanisms in the physiological environment and the dynamic homeostasis of mammary gland development and lactation are not clear. In the present study, sequencing data analysis of early and late lactation uncovered a total of 1,487 unique miRNAs, including 45 novel miRNA candidates and 1,442 known and conserved miRNAs, of which 758 miRNAs were co-expressed and 378 differentially expressed with P < 0.05. Moreover, 76 non-redundant target genes were annotated in 347 GO consortiums, with 3,143 candidate target genes grouped into 33 pathways. Additionally, 18 predicted target genes of 214 miRNAs were directly annotated in mammary gland development and used to construct regulatory networks based on GO annotation and the KEGG pathway. The expression levels of seven known miRNAs and three novel miRNAs were examined using quantitative real-time PCR. The results showed that miRNAs might play important roles in early and late lactation during dairy goat mammary gland development, which will be helpful to obtain a better understanding of the genetic control of mammary gland lactation and development.