Analyses of Long Non-Coding RNA and mRNA profiling using RNA sequencing during the pre-implantation phases in pig endometrium

The TCONS_00046732/chi-miR-135b-5p/PRLR regulatory axis promotes endometrial epithelial cells growth through the PI3K-Akt signaling pathway

Non-coding RNAs are considered key regulatory factors in biological and reproductive physiological processes in mammals. However, the molecular functions of long noncoding RNAs (lncRNAs) in regulating dynamic uterine development and reproductive capacity during sexual maturation remain unclear. This study analyzed the expression characteristics and molecular functions of lncRNAs in uterine tissues from 20 goats at four specific time points during sexual maturation: day 1 after birth (D1), 2 months (M2), 4 months (M4), and 6 months (M6), finding that stage-specific DE lncRNAs may regulate cell proliferation, apoptosis, hormone secretion, metabolism, and immune response through multiple action modes. Within the lncRNA-miRNA-mRNA network, a novel lncRNA, TCONS_00046732, associated with uterine development, exhibited significantly higher expression during sexual maturity compared to the prepubertal stage, correlating positively with PRLR and negatively with chi-miR-135b-5p. FISH and IF analyses revealed significant co-localization and distinct expression patterns of TCONS_00046732, chi-miR-135b-5p, and PRLR in the endometrial epithelium. Further experiments confirmed that TCONS_00046732 competitively binds to chi-miR-135b-5p to upregulate PRLR, thereby activating the PI3K-Akt signaling pathway, promoting primary endometrial epithelial cell proliferation, G1-to-S phase transition, and inhibiting apoptosis. These findings enhance our understanding of uterine molecular regulation and provide key insights into the molecular basis of goat sexual development.

LnCRNAs in the Regulation of Endometrial Receptivity for Embryo Implantation

The development of endometrial receptivity is crucial for successful embryo implantation and the initiation of pregnancy. Understanding the molecular regulatory processes that transform the endometrium into a receptive phase is essential for enhancing implantation rates in fertility treatments, such as in vitro fertilization (IVF). Long non-coding RNAs (lncRNAs) play a pivotal role as gene regulators and have been examined in the endometrium. This review offers current insights into the role of lncRNAs in regulating endometrial receptivity. Considering the significant variation in endometrial remodeling among species, we summarize the key events in the human endometrial cycle and discuss the identified lncRNAs in both humans and other species, which may play a crucial role in establishing receptivity. Notably, there are 742 lncRNAs in humans and 4438 lncRNAs that have the potential to modulate endometrial receptivity. Additionally, lncRNAs regulating matrix metalloproteinases (MMPs) and Let-7 have been observed in both species. Future investigations should explore the potential of lncRNAs as therapeutic targets and/or biomarkers for diagnosing and improving endometrial receptivity in human fertility therapy.

Dynamic intrauterine crosstalk promotes porcine embryo implantation during early pregnancy

Dynamic crosstalk between the embryo and mother is crucial during implantation. Here, we comprehensively profile the single-cell transcriptome of pig peri-implantation embryos and corresponding maternal endometrium, identifying 4 different lineages in embryos and 13 cell types in the endometrium. Cell-specific gene expression characterizes 4 distinct trophectoderm subpopulations, showing development from undifferentiated trophectoderm to polar and mural trophectoderm. Dynamic expression of genes in different types of endometrial cells illustrates their molecular response to embryos during implantation. Then, we developed a novel tool, ExtraCellTalk, generating an overall dynamic map of maternal-foetal crosstalk using uterine luminal proteins as bridges. Through cross-species comparisons, we identified a conserved RBP4/STRA6 pathway in which embryonic-derived RBP4 could target the STRA6 receptor on stromal cells to regulate the interaction with other endometrial cells. These results provide insight into the maternal-foetal crosstalk during embryo implantation and represent a valuable resource for further studies to improve embryo implantation.

Maternal Hypermethylated Genes Contribute to Intrauterine Growth Retardation of Piglets in Rongchang Pigs

The placenta is a crucial determinant of fetal survival, growth, and development. Deficiency in placental development directly causes intrauterine growth retardation (IUGR). IUGR can lead to fetal growth restriction and an increase in the mortality rate. The genetic mechanisms underlying IUGR development, however, remain unclear. In the present study, we integrated whole-genome DNA methylation and transcriptomic analyses to determine distinct gene expression patterns in various placental tissues to identify pivotal genes that are implicated with IUGR development. By performing RNA-sequencing analysis, 1487 differentially expressed genes (DEGs), with 737 upregulated and 750 downregulated genes, were identified in IUGR pigs (H_IUGR) compared with that in normal birth weight pigs (N_IUGR) (p < 0.05); furthermore, 77 miRNAs, 1331 lncRNAs, and 61 circRNAs were differentially expressed. The protein-protein interaction network analysis revealed that among these DEGs, the genes GNGT1, ANXA1, and CDC20 related to cellular developmental processes and blood vessel development were the key genes associated with the development of IUGR. A total of 495,870 differentially methylated regions were identified between the N_IUGR and H_IUGR groups, which included 25,053 differentially methylated genes (DMEs); moreover, the overall methylation level was higher in the H_IUGR group than in the N_IUGR group. Combined analysis showed an inverse correlation between methylation levels and gene expression. A total of 1375 genes involved in developmental processes, tissue development, and immune system regulation exhibited methylation differences in gene expression levels in the promoter regions and gene ontology regions. Five genes, namely, ANXA1, ADM, NRP2, SHH, and SMAD1, with high methylation levels were identified as potential contributors to IUGR development. These findings provide valuable insights that DNA methylation plays a crucial role in the epigenetic regulation of gene expression and mammalian development and that DNA-hypermethylated genes contribute to IUGR development in Rongchang pigs.

Unravelling the impact of epigenetic mechanisms on offspring growth, production, reproduction and disease susceptibility

Epigenetic mechanisms, such as DNA methylation, histone modifications and non-coding RNA molecules, play a critical role in gene expression and regulation in livestock species, influencing development, reproduction and disease resistance. DNA methylation patterns silence gene expression by blocking transcription factor binding, while histone modifications alter chromatin structure and affect DNA accessibility. Livestock-specific histone modifications contribute to gene expression and genome stability. Non-coding RNAs, including miRNAs, piRNAs, siRNAs, snoRNAs, lncRNAs and circRNAs, regulate gene expression post-transcriptionally. Transgenerational epigenetic inheritance occurs in livestock, with environmental factors impacting epigenetic modifications and phenotypic traits across generations. Epigenetic regulation revealed significant effect on gene expression profiling that can be exploited for various targeted traits like muscle hypertrophy, puberty onset, growth, metabolism, disease resistance and milk production in livestock and poultry breeds. Epigenetic regulation of imprinted genes affects cattle growth and metabolism while epigenetic modifications play a role in disease resistance and mastitis in dairy cattle, as well as milk protein gene regulation during lactation. Nutri-epigenomics research also reveals the influence of maternal nutrition on offspring's epigenetic regulation of metabolic homeostasis in cattle, sheep, goat and poultry. Integrating cyto-genomics approaches enhances understanding of epigenetic mechanisms in livestock breeding, providing insights into chromosomal structure, rearrangements and their impact on gene regulation and phenotypic traits. This review presents potential research areas to enhance production potential and deepen our understanding of epigenetic changes in livestock, offering opportunities for genetic improvement, reproductive management, disease control and milk production in diverse livestock species.

DIA-based quantitative proteomic analysis of porcine endometrium in the peri-implantation phase

The 12th day of gestation is a critical period for embryo loss and the beginning of imminent implantation in sows. Data independent acquisition (DIA) technology is one of the high-throughput, high-resolution and reproducible proteomics technologies for large-scale digital qualitative and quantitative research. The aim of this study was to identify and characterize the protein abundance landscape of Yorkshire pig endometrium on the 12th day of pregnancy (P12) and estrous cycle (C12) using DIA proteomics. A total of 1251 differentially abundant proteins (DAPs) were identified, of which 882 were up-regulated and 369 were down-regulated at P12. Functional enrichment analysis showed that the identified proteins were related to metabolism, biosynthesis and signaling pathways. Three proteins were selected for Western blot (WB) validation and the results were consistent with the DIA data. Further combined with transcriptome data, fibrinogen like 2 (FGL2) and S100 calcium binding protein A8 (S100A8) were verified to be highly abundant in the P12 endometrial epithelium. In summary, there were significantly different abundance of proteome profiles in C12 and P12 endometrium, suggesting that DAPs are associated with changes in endometrial receptivity, which laid the foundation for further research on related regulatory mechanisms. SIGNIFICANCE: The 12th day of gestation is an important point in the peri-implantation period of pigs, when the endometrium presents a receptive state under the stimulation of estrogen. DIA proteomics technology is an emerging protein identification technology in recent years, which can obtain protein information through comprehensive and unbiased scanning. In this study, DIA technology was used to characterize endometrial proteins in pigs during the peri-implantation period. The results showed that higher protein abundance was detected using the DIA technique, and some of these DAPs may be involved in regulating embryo implantation. This study will help to better reveal the related proteins involved in embryo implantation, and lay a foundation for further research on the mechanism of endometrial regulation of embryo implantation. SIGNIFICANCE OF THE STUDY: The 12th day of gestation is an important point in the peri-implantation period of pigs, when the endometrium presents a receptive state under the stimulation of estrogen. DIA proteomics technology is an emerging protein identification technology in recent years, which can obtain protein information through comprehensive and unbiased scanning. In this study, DIA technology was used to characterize endometrial proteins in pigs during the peri-implantation period. The results showed that higher protein abundance was detected using the DIA technique, and some of these DAPs may be involved in regulating embryo implantation. This study will help to better reveal the related proteins involved in embryo implantation, and lay a foundation for further research on the mechanism of endometrial regulation of embryo implantation.

Identification cloning and functional analysis of novel natural antisense lncRNA CFL1-AS1 in cattle

Long noncoding RNAs have been identified as important regulators of gene expression and animal development. The expression of natural antisense transcripts (NATs) transcribed in the opposite direction to protein-coding genes is usually positively correlated with the expression of homologous sense genes and is the key factor for expression. Here, we identified a conserved noncoding antisense transcript, CFL1-AS1, that plays an important role in muscle growth and development. CFL1-AS1 overexpression and knockout vectors were constructed and transfected into 293T and C2C12 cells. CFL1-AS1 positively regulated CFL1 gene expression, and the expression of CFL2 was also downregulated when CFL1-AS1 was knocked down. CFL1-AS1 promoted cell proliferation, inhibited apoptosis and participated in autophagy. This study expands the research on NATs in cattle and lays a foundation for the study of the biological function of bovine CFL1 and its natural antisense chain transcript CFL1-AS1 in bovine skeletal muscle development. The discovery of this NAT can provide a reference for subsequent genetic breeding and data on the characteristics and functional mechanisms of NATs.

An Integrative ATAC-Seq and RNA-Seq Analysis of the Endometrial Tissues of Meishan and Duroc Pigs

Meishan pigs are a well-known indigenous pig breed in China characterized by a high fertility. Notably, the number of endometrial grands is significantly higher in Meishan pigs than Duroc pigs. The characteristics of the endometrial tissue are related to litter size. Therefore, we used the assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-sequencing (RNA-seq) to analyze the mechanisms underlying the differences in fecundity between the breeds. We detected the key transcription factors, including Double homeobox (Dux), Ladybird-like homeobox gene 2 (LBX2), and LIM homeobox 8 (Lhx8), with potentially pivotal roles in the regulation of the genes related to endometrial development. We identified the differentially expressed genes between the breeds, including SOX17, ANXA4, DLX3, DMRT1, FLNB, IRF6, CBFA2T2, TFCP2L1, EFNA5, SLIT2, and CYFIP2, with roles in epithelial cell differentiation, fertility, and ovulation. Interestingly, ANXA4, CBFA2T2, and TFCP2L1, which were upregulated in the Meishan pigs in the RNA-seq analysis, were identified again by the integration of the ATAC-seq and RNA-seq data. Moreover, we identified genes in the cancer or immune pathways, FoxO signaling, Wnt signaling, and phospholipase D signaling pathways. These ATAC-seq and RNA-seq analyses revealed the accessible chromatin and potential mechanisms underlying the differences in the endometrial tissues between the two types of pigs.

LncRNA-Anrel promotes the proliferation and migration of synovial fibroblasts through regulating miR-146a-mediated annexin A1 expression

OBJECTIVE

Increasing evidence demonstrates that long non-coding RNAs (lncRNAs) are closely related to multiple human autoimmune diseases, and their dysregulation is tightly linked to inflammation and disease progression. Nonetheless, little is known about the consequences of aberrant expression of lncRNAs during rheumatoid arthritis (RA) development. In this study, we screened for the expressions of lncRNAs in RA synovial fibroblasts (RA-SF) and investigated their functions in RA-SF proliferation and migration, and the relevant underlying mechanisms.

METHODS

The lncRNAs expression profiles were interrogated with microarrays. The expressions of key lncRNAs were confirmed in synovial fibroblasts from RA patients and MH7A cells using qRT-PCR. Proliferations and migrations of MH7A and HFL-1 cells were evaluated using CCK-8 assay and cell migration assay kits, respectively. The expression of inflammatory cytokines (IL-6, IL-1β, and TNF-α) and cell migration related proteins (MMP-1 and MMP-3) were evaluated using qRT-PCR and western blotting. Collagen type II-induced arthritis (CIA) in mice was used as an animal model of RA.

RESULTS

Nine lncRNAs were significantly altered in RA-SF, of which lncRNA-000239 showing the most significant upregulation. Overexpression of lncRNA-000239 significantly enhanced the proliferation and migration of human RS-SF cells (MH7A), while the opposite effect was observed with lncRNA-000239 silencing. Importantly, lncRNA-000239 enhanced annexin A1 expression by upregulating the expression of miR-146a. Moreover, locally enhanced expression of lncRNA-000239 promoted the onset of arthritis in CIA.

CONCLUSION

These data indicate that lncRNA-000239 upregulates annexin A1 expression via miR-146a and thus, promotes the proliferation and migration of RA-SF. This highlights a potential role of lncRNA-000239 as an inflammatory factor of RA.

Astragalosides inhibit proliferation of fibroblast-like synoviocytes in experimental arthritis by modulating LncRNA S56464.1/miR-152-3p/Wnt1 signaling axis

AIM

Astragalus membranaceus (Fisch.) Bunge., the dried root of the plant A. membranaceus, is widely used in the treatment of rheumatoid arthritis (RA) in many Chinese herbal remedies. Astragalosides (AST) is the primary medicinal ingredient of A. membranaceus and has a therapeutic effect on RA, but the specific mechanism of this effect has yet to be elucidated.

METHODS

In this study, MTT and flow cytometry were used to determine the effects of AST on fibroblast-like synoviocyte (FLS) proliferation and cell cycle progression. Additionally, real-time quantitative polymerase chain reaction and Western blotting were used to determine the effects of AST on the LncRNA S56464.1/miR-152-3p/Wnt1 signaling axis and on critical genes that are essential to the Wnt pathway.

RESULTS

The data showed that after the administration of AST, FLS proliferation and LncRNA S56464.1, β-catenin, C-myc, Cyclin D1, and p-GSK-3β(Ser9)/GSK-3β expression were significantly reduced, and miR-152 and SFRP4 expression was notably increased.

CONCLUSION

These results suggest that AST can inhibit FLS proliferation by modulating the LncRNA S56464.1/miR-152-3p/Wnt1 signaling axis and that AST may be a potential therapeutic drug for RA.

Identification and characterization of circRNAs in peri-implantation endometrium between Yorkshire and Erhualian pigs

BACKGROUND

One of the most critical periods for the loss of pig embryos is the 12th day of gestation when implantation begins. Recent studies have shown that non-coding RNAs (ncRNAs) play important regulatory roles during pregnancy. Circular RNAs (circRNAs) are a kind of ubiquitously expressed ncRNAs that can directly regulate the binding proteins or regulate the expression of target genes by adsorbing micro RNAs (miRNA).

RESULTS

We used the Illumina Novaseq6,000 technology to analyze the circRNA expression profile in the endometrium of three Erhualian (EH12) and three Yorkshire (YK12) pigs on day 12 of gestation. Overall, a total of 22,108 circRNAs were identified. Of these, 4051 circRNAs were specific to EH12 and 5889 circRNAs were specific to YK12, indicating a high level of breed specificity. Further analysis showed that there were 641 significant differentially expressed circRNAs (SDEcircRNAs) in EH12 compared with YK12 (FDR < 0.05). Functional enrichment of differential circRNA host genes revealed many pathways and genes associated with reproduction and regulation of embryo development. Network analysis of circRNA-miRNA interactions further supported the idea that circRNAs act as sponges for miRNAs to regulate gene expression. The prediction of differential circRNA binding proteins further explored the potential regulatory pathways of circRNAs. Analysis of SDEcircRNAs suggested a possible reason for the difference in embryo survival between the two breeds at the peri-implantation stage.

CONCLUSIONS

Together, these data suggest that circRNAs are abundantly expressed in the endometrium during the peri-implantation period in pigs and are important regulators of related genes. The results of this study will help to further understand the differences in molecular pathways between the two breeds during the critical implantation period of pregnancy, and will help to provide insight into the molecular mechanisms that contribute to the establishment of pregnancy and embryo loss in pigs.

Proteomic profiles and the function of RBP4 in endometrium during embryo implantation phases in pigs

BACKGROUND

Endometrial receptivity plays a vital role in the success of embryo implantation. However, the temporal proteomic profile of porcine endometrium during embryo implantation is still unclear.

RESULTS

In this study, the abundance of proteins in endometrium on days 9, 10, 11, 12, 13, 14, 15 and 18 of pregnancy (D9, 10, 11, 12, 13, 14, 15 and 18) was profiled via iTRAQ technology. The results showed that 25, 55, 103, 91, 100, 120, 149 proteins were up-regulated, and 24, 70, 169, 159, 164, 161, 198 proteins were down-regulated in porcine endometrium on D10, 11, 12, 13, 14, 15 and 18 compared with that on D9, respectively. Among these differentially abundance proteins (DAPs), Multiple Reaction Monitoring (MRM) results indicated that S100A9, S100A12, HRG and IFI6 were differentially abundance in endometrial during embryo implantation period. Bioinformatics analysis showed that the proteins differentially expressed in the 7 comparisons were involved in important processes and pathways related to immunization, endometrial remodeling, which have a vital effect on embryonic implantation.

CONCLUSION

Our results reveal that retinol binding protein 4 (RBP4) could regulate the cell proliferation, migration and apoptosis of endometrial epithelial cells and endometrial stromal cells to affect embryo implantation. This research also provides resources for studies of proteins in endometrium during early pregnancy.

Genome-Wide Analysis of Long Noncoding RNAs in Porcine Intestine during Weaning Stress

Long noncoding RNAs (lncRNAs) play crucial roles in various biological processes, and they are considered to be closely associated with the pathogenesis of intestinal diseases. However, the role and expression of lncRNAs in intestinal damage during weaning stress remain unknown. Herein, we investigated the expression profiles of jejunal tissue from weaning piglets at 4 and 7 d after weaning (groups W4 and W7, respectively) and from suckling piglets on the same days (groups S4 and S7, respectively). Genome-wide analysis of lncRNAs was also performed using RNA sequencing technology. A total of 1809 annotated lncRNAs and 1612 novel lncRNAs were obtained from the jejunum of piglets. In W4 vs. S4, a total of 331 lncRNAs showed significant differential expression, and a total of 163 significantly differentially expressed lncRNAs (DElncRNAs) was identified in W7 vs. S7. Biological analysis indicated that DElncRNAs were involved in intestinal diseases, inflammation, and immune functions, and were mainly enriched in the Jak-STAT signaling pathway, inflammatory bowel disease, T cell receptor signaling pathway, B cell receptor signaling pathway and intestinal immune network for IgA production. Moreover, we found that lnc_000884 and target gene KLF5 were significantly upregulated in the intestine of weaning piglets. The overexpression of lnc_000884 also significantly promoted the proliferation and depressed apoptosis of IPEC-J2 cells. This result suggested that lnc_000884 may contribute to repairing intestinal damage. Our study identified the characterization and expression profile of lncRNAs in the small intestine of weaning piglets and provided new insights into the molecular regulation of intestinal damage during weaning stress.

Change in Long Non-Coding RNA Expression Profile Related to the Antagonistic Effect of Clostridium perfringens Type C on Piglet Spleen

LncRNAs play important roles in resisting bacterial infection via host immune and inflammation responses. Clostridium perfringens (C. perfringens) type C is one of the main bacteria causing piglet diarrhea diseases, leading to major economic losses in the pig industry worldwide. In our previous studies, piglets resistant (SR) and susceptible (SS) to C. perfringens type C were identified based on differences in host immune capacity and total diarrhea scores. In this paper, the RNA-Seq data of the spleen were comprehensively reanalyzed to investigate antagonistic lncRNAs. Thus, 14 lncRNAs and 89 mRNAs were differentially expressed (DE) between the SR and SS groups compared to the control (SC) group. GO term enrichment, KEGG pathway enrichment and lncRNA-mRNA interactions were analyzed to identify four key lncRNA targeted genes via MAPK and NF-κB pathways to regulate cytokine genes (such as TNF-α and IL-6) against C. perfringens type C infection. The RT-qPCR results for six selected DE lncRNAs and mRNAs are consistent with the RNA-Seq data. This study analyzed the expression profiling of lncRNAs in the spleen of antagonistic and sensitive piglets and found four key lncRNAs against C. perfringens type C infection. The identification of antagonistic lncRNAs can facilitate investigations into the molecular mechanisms underlying resistance to diarrhea in piglets.

Interaction between Microbes and Host in Sow Vaginas in Early Pregnancy

Extensive research has explored the causes of embryo losses during early pregnancy by analyzing interaction mechanisms in sows' uterus, ignoring the importance of the lower reproductive tract in pregnancy development regulation. Despite recent progress in understanding the diversity of vaginal microbes under different physiological states, the dynamic of sows' vaginal microbiotas during pregnancy and the interaction between vaginal microbes and the host are poorly understood. Here, we performed a comprehensive analysis of sows' vaginal microbial communities in early pregnancy coupled with overall patterns of vaginal mucosal epithelium gene expression. The vaginal microbiota was analyzed by 16s rRNA or metagenome sequencing, and the vaginal mucosal epithelium transcriptome was analyzed by RNA sequencing, followed by integration of the data layers. We found that the sows' vaginal microbiotas in early pregnancy develop dynamically, and there is a homeostasis balance of Firmicutes and Proteobacteria. Subsequently, we identified two pregnancy-specific communities, which play diverse roles. The microbes in the vagina stimulate the epithelial cells, while vaginal epithelium changes its structure and functions in response to stimulation. These changes produce specific inflammation responses to promote pregnancy development. Our findings demonstrate the interaction between microbes and host in the sow vagina in early pregnancy to promote pregnancy development, meanwhile providing a reference data set for the study of targeted therapies of microbial homeostasis dysregulation in the female reproductive tract. IMPORTANCE This work sheds light on the dynamics of the sow vaginal microbiotas in early pregnancy and its roles in pregnancy development. Furthermore, this study provides insight into the functional mechanisms of reproductive tract microbes by outlining vaginal microbe-host interactions, which might identify new research and intervention targets for improving pregnancy development by modulating lower reproductive tract microbiota.

Screening and functional validation of lipid metabolism-related lncRNA-46546 based on the transcriptome analysis of early embryonic muscle tissue in chicken

OBJECTIVE

The study was conducted to screen differentially expressed long noncoding RNA (lncRNA) in chickens by high-throughput sequencing and explore its mechanism of action on intramuscular fat deposition.

METHODS

Herein, Rose crown and Cbb broiler chicken embryo breast and leg muscle lncRNA and mRNA expression profiles were constructed by RNA sequencing. A total of 96 and 42 differentially expressed lncRNAs were obtained in Rose crown vs Cobb broiler chicken breast and leg muscle, respectively. lncRNA-ENSGALT00000046546, with high interspecific variability and a potential regulatory role in lipid metabolism, and its predicted downstream target gene 1-acylglycerol-3-phosphate-O-acyltransferase 2 (AGPAT2), were selected for further study on the preadipocytes.

RESULTS

lncRNA-46546 overexpression in chicken preadipocyte 2 cells significantly increased (p<0.01) the expression levels of AGPAT2 and its downstream genes diacylglycerol acyltransferase 1 and diacylglycerol acyltransferase 2 and those of the fat metabolism-related genes peroxisome proliferator-activated receptor γ, CCAAT/enhancer binding protein α, fatty acid synthase, sterol regulatory element-binding transcription factor 1, and fatty acid binding protein 4. The lipid droplet concentration was higher in the overexpression group than in the control cells, and the triglyceride content in cells and medium was also significantly increased (p<0.01).

CONCLUSION

This study preliminarily concludes that lncRNA-46546 may promote intramuscular fat deposition in chickens, laying a foundation for the study of lncRNAs in chicken early embryonic development and fat deposition.

New Roles for EVs, miRNA and lncRNA in Bovine Embryo Implantation

The sine qua non of new life is fertilization. However, approximately 50% of fertilized eggs/blastocysts in cattle and up to 75% of those from human assisted reproductive procedures fail during the first 3 to 4 weeks of pregnancy, including peri-implantation periods. In these periods, blastocyst hatching and implantation to the maternal endometrium proceeds, during which physiological events such as epithelial-mesenchymal transition (EMT) and trophoblast cell fusion occur. Quite recently, extracellular vesicles (EVs) with micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs) have been found to play a pivotal role for the establishment of the proper uterine environment required for peri-implantation processes to proceed. New findings of EVs, miRNA, and lncRNAs will be described and discussed to elucidate their connections with conceptus implantation to the maternal endometrium.

Circ0001470 Acts as a miR-140-3p Sponge to Facilitate the Progression of Embryonic Development through Regulating PTGFR Expression

Embryonic implantation and development are vital in early pregnancy and assisted reproduction. Circular RNAs (circRNAs) are involved in the two physiological processes and thus regulate animal reproduction. However, their specific regulatory functions and mechanisms remain unclear. Here, a novel circ0001470, originating from the porcine GRN gene, differentially expressed on day 18 versus day 32 of gestation in Meishan and Yorkshire pigs was screened. The circularization characteristic of circ0001470 was identified based on divergent primer amplification, Sanger sequencing, RNase digestion, and RNA nuclear-cytoplasmic fractionation. Functionally, circ0001470 can promote cell proliferation and cycle progression of endometrial epithelial cells (EECs) and also inhibit apoptosis of EECs using CCK-8 assays and flow cytometry analyses. Mechanistically, bioinformatics database prediction, luciferase screening, RNA immunoprecipitation (RIP), RNA-pull down, and FISH co-localization experiments revealed that the circ0001470 acted as a competing endogenous RNA (ceRNA) through sponging miR-140-3p to regulate downstream PTGFR expression. Moreover, in vivo assays revealed that mmu_circGRN promoted embryonic development by affecting the expression of PTGFR, which can activate the MAPK reproduction pathway and facilitate pregnancy maintenance. This study enriched our understanding of circRNAs in embryo implantation and development by deciding the fate of EECs.

Integrative Analysis of lncRNA-miRNA-mRNA Regulatory Network Reveals the Key lncRNAs Implicated Potentially in the Differentiation of Adipocyte in Goats

Goats are popular in China because of their superior meat quality, delicate flesh, and unique flavor. Long noncoding RNAs (lncRNAs) play important roles in transcriptional and post-transcriptional regulation of gene expression. However, the effects of lncRNAs on adipocyte differentiation in goat has not been fully elucidated yet. In this investigation, we performed RNA-Seq analysis of intramuscular and subcutaneous adipocytes from Jianzhou Daer goat before and after differentiation, including both intramuscular preadipocytes (IMPA) vs. intramuscular adipocytes (IMA) and subcutaneous preadipocytes (SPA) vs. subcutaneous adipocytes (SA). A total of 289.49 G clean reads and 12,519 lncRNAs were obtained from 20 samples. In total, 3,733 differentially expressed RNAs (182 lncRNAs and 3,551 mRNAs) were identified by pairwise comparison. There were 135 differentially expressed lncRNAs (DELs) specific to intramuscular adipocytes, 39 DELs specific to subcutaneous adipocytes, and 8 DELs common to both adipocytes in these 182 DELs. Some well-known and novel pathways associated with preadipocyte differentiation were identified: fat acid metabolism, TGF-beta signaling pathway and PI3K-Akt signaling pathway. By integrating miRNA-seq data from another study, we also identified hub miRNAs in both types of fat cells. Our analysis revealed the unique and common lncRNA-miRNA-mRNA networks of two kinds of adipocytes. Several lncRNAs that regulate potentially goat preadipocyte differentiation were identified, such as XR_001918 647.1, XR_001917728.1, XR_001297263.2 and LNC_004191. Furthermore, our findings from the present study may contribute to a better understanding of the molecular mechanisms underlying in goat meat quality and provide a theoretical basis for further goat molecular breeding.

Effect of MiR-100-5p on proliferation and apoptosis of goat endometrial stromal cell in vitro and embryo implantation in vivo

The growth of endometrial stromal cells (ESCs) at implantation sites may be a potential factor affecting the success rate of embryo implantation. Incremental proofs demonstrated that ncRNAs (e.g. miRNAs, lncRNAs and circRNAs) were involved in various biological procedures, including proliferation and apoptosis. In this study, the role of miR‐100‐5p on proliferation and apoptosis of goat ESCs in vitro and embryo implantation in vivo was determined. The mRNA expression of miR‐100‐5p was significantly inhibited in the receptive phase (RE) rather than in the pre‐receptive phase (PE). Overexpression of miR‐100‐5p suppressed ESCs proliferation and induced apoptosis. The molecular target of MiR‐100‐5p, HOXA1, was confirmed by 3′‐UTR assays. Meanwhile, the product of HOXA1 mRNA RT‐PCR increased in the RE more than that in the PE. The HOXA1‐siRNA exerted significant negative effects on growth arrest. Instead, incubation of ESCs with miR‐100‐5p inhibitor or overexpressed HOXA1 promoted the cell proliferation. In addition, Circ‐9110 which acted as a sponge for miR‐100‐5p reversed the relevant biological effects of miR‐100‐5p. The intrinsic apoptosis pathway was suppressed in ESCs, revealing a crosstalk between Circ‐9110/miR‐100‐5p/HOXA1 axis, PI3K/AKT/mTOR, and ERK1/2 pathways. To further evaluate the progress in study on embryo implantation regulating mechanism of miR‐100‐5p in vivo, the pinopodes of two phases were observed and analysed, suggesting that, as similar as in situ, miR‐100‐5p was involved in significantly regulating embryo implantation in vivo. Mechanistically, miR‐100‐5p performed its embryo implantation function through regulation of PI3K/AKT/mTOR and ERK1/2 pathways by targeting Circ‐9110/miR‐100‐5p/HOXA1 axis in vivo.

Heat Stress Impairs Maternal Endometrial Integrity and Results in Embryo Implantation Failure by Regulating Transport-Related Gene Expression in Tongcheng Pigs

Heat stress (HS) poses a significant threat to production and survival in the global swine industry. However, the molecular regulatory effects of heat stress on maternal endometrial cells are poorly understood in pigs during early embryo implantation. In this study, we systematically examined morphological changes in the endometrium and the corresponding regulation mechanism in response to HS by combining scanning electron microscopy (SEM), hematoxylin/eosin (H&E) staining, western blot, and RNA-seq analyses. Our results showed that HS led to porcine endometrium damage and endometrial thinness during embryo implantation. The expression levels of cell adhesion-related proteins, including N-cadherin and E-cadherin, in the uterus were significantly lower in the heat stress group (39 ± 1 °C, n = 3) than in the control group (28 ± 1 °C, n = 3). A total of 338 up-regulated genes and 378 down-regulated genes were identified in porcine endometrium under HS. The down-regulated genes were found to be mainly enriched in the pathways related to the microtubule complex, immune system process, and metalloendopeptidase activity, whereas the up-regulated genes were mainly involved in calcium ion binding, the extracellular region, and molecular function regulation. S100A9 was found to be one of the most significant differentially expressed genes (DEGs) in the endometrium under HS, and this gene could promote proliferation of endometrial cells and inhibit their apoptosis. Meanwhile, HS caused endometrial epithelial cell (EEC) damage and inhibited its proliferation. Overall, our results demonstrated that HS induced uterine morphological change and tissue damage by regulating the expression of genes associated with calcium ions and amino acid transport. These findings may provide novel molecular insights into endometrial damage under HS during embryo implantation.

Characterization of lncRNA functioning in ovine conceptuses and endometria during the peri-implantation period

In ruminants, RNA-sequence analyses have revealed many characteristics of transcripts expressed in conceptuses (embryo and extraembryonic membrane) during peri-implantation periods; however, lncRNA profiles are yet characterized. In this study, we aimed to characterize the lncRNA expression profile in conceptuses during peri-implantation periods in sheep. We analyzed the RNA-sequence data of ovine conceptuses and endometria obtained from pregnant animals on days 15, 17, 19 and 21 (day 0 = day of estrus, n = 3 or 4/day). We predicted the protein coding ability of the assembled transcripts to identify the lncRNA candidates. This analysis identified 8808 lncRNAs, 3423 of which were novel lncRNAs. Gene ontology analysis revealed that lncRNA target genes were enriched for biological processes involved in the respiratory electron transport chain (RETC). qPCR analysis demonstrated that the expression levels on transcripts encoding RETC such as mitochondrially encoded cytochrome c oxidase II (MTCO2) and mitochondria DNA copy number in conceptuses were not increased on P21, although western blotting analysis and immunohistochemistry demonstrated that MTCO2 protein in conceptuses was increased on P21. NAD/NADH assay revealed that NADH level in conceptuses was increased on P21. These results indicate that lncRNAs could regulate the RETC through post-transcriptional levels in the conceptuses. Therefore, lncRNA is a potential new regulator in ovine conceptus development during peri-implantation periods.

Transcriptomic analysis of human endometrial stromal cells during early embryo invasion

PURPOSE

During early embryo invasion (48 h after embryo attachment), what functional changes accompany dynamic gene expression alterations in human endometrial stromal cells?

METHOD

In the present study, primary human endometrial stromal cells (phESCs) were cultured. After in vitro decidualization, primary human endometrial stromal cells (phESCs) were cultured with blastocysts for 48 h. During this process, blastocysts attached and invaded the phESCs (embryo-invaded primary human endometrial stromal cells, ehESCs). We performed comprehensive transcriptomic profiling of phESCs (two replicates) and ehESCs (five replicates) and analyzed the differentially expressed gene (DEGs) sets for gene ontology (GO) terms and Kyoto encyclopaedia of genes and genomes (KEGG) pathway enrichment. To analyse potential connectivity patterns between the transcripts in these DEG sets, a protein-protein interaction (PPI) network was constructed using the STRING database.

RESULTS

A total of 592 DEGs were identified between phESCs and ehESCs after embryo invasion. Primary human endometrial stromal cells underwent significant transcriptomic changes that occur in a stepwise fashion. Oxidative phosphorylation, mitochondrial organization, and P53 signalling pathways were significantly altered in phESCs after embryo invasion. EP300 may play a key role in regulating transcription via chromatin remodelling to facilitate the adaptive gene expression changes that occur during embryo invasion.

CONCLUSIONS

Our data identify dynamic transcriptome changes that occur in endometrial stromal cells within 48 h after embryo invasion. The pathways that we found to be enriched in phESCs after embryo invasion (oxidative phosphorylation, mitochondrial organization, and P53 signalling) may represent novel mechanisms underlying embryo implantation, and may illuminate the reasons that some women experience reproductive failure.Key messagesHuman endometrial stromal cells have undergone changes in gene expression regulation and signalling pathways during the embryo invasion.Mitochondrial-oxidative phosphorylation changes in human stromal cells manifested as down-regulation of gene expression in the electron transport chain.TP53 signalling pathway and transcriptional regulator EP300 assist stromal cells to get adaptive changes during embryo invasion phase.

S100A8 expression in oviduct mucosal epithelial cells is regulated by estrogen and affects mucosal immune homeostasis

Chronic inflammation can cause oviduct mucosal damage and immune dysfunction, leading to infertility, early pregnancy loss, ectopic pregnancy, tumors, and a decrease in reproductive capacities in female animals. Estrogen can suppress immune responses in different tissues and oviducts, and regulate the oviduct immune balance; however, the underlying mechanisms remain unclear. The objective of this study was to explore the mechanism of estrogen-regulated oviduct mucosal immunity and discover new estrogen targets for regulating oviduct mucosal immune homeostasis. Sheep oviduct epithelial cells (SOECs) were treated with 17-β estradiol (E2). Transcriptome sequencing and analysis showed differentially expressed S100 calcium-binding protein A (S100A) genes that may participate in the oviduct mucosa immunoregulation of estrogen. Quantitative polymerase chain reaction and immunocytochemistry analysis showed that S100A8 expression changed dynamically in E2-treated SOECs and peaked after 7 h of treatment. Estrogen nuclear receptors and G protein-coupled membrane receptors promoted E2-dependent S100A8 upregulation. The S100A8 gene was disrupted using the clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 method. Levels of inflammatory factors interleukin (IL)-1β and IL-4 were significantly upregulated in S100A8-knockdown SOECs, whereas those of the anti-inflammatory factor IL-10 was downregulated. Following S100A8 knockdown in SOECs treated with E2 for 7 h, IL-10 levels increased significantly. Estrogen affected oviduct mucosa immune function and dynamically regulated S100A8 in SOECs. S100A8 knockdown caused an excessive immune response, indicating that S100A8 is beneficial for maintaining immune homeostasis in the oviduct mucosa. Moreover, estrogen can compensate for the effect of S100A8 knockdown by upregulating IL-10.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

Comprehensive analysis of lncRNAs involved in skeletal muscle development in ZBED6-knockout Bama pigs

BACKGROUND

The mutation of insulin-like growth factor 2 (IGF2 mutation) that a single-nucleotide substitution (G→A) in the third intron of IGF2 abrogates the interaction with zinc finger BED-type containing 6 (ZBED6) and leads to increased muscle mass in pigs. IGF2 mutation knock-in (IGF2 KI) and ZBED6 knockout (ZBED6 KO) lead to changes in IGF2 expression and increase muscle mass in mice and pigs. Long noncoding RNAs (lncRNAs) may participate in numerous biological processes, including skeletal muscle development. However, the role of the ZBED6-lncRNA axis in skeletal muscle development is poorly characterized.

RESULTS

In this study, we assembled transcriptomes using RNA-seq data published in previous studies by our group and identified 11,408 known lncRNAs and 2269 potential lncRNAs in seven tissues, heart, longissimus dorsi, gastrocnemius muscle, liver, spleen, lung and kidney, of ZBED6 KO (lean mass model) and WT Bama pigs. ZBED6 affected the expression of 1570 lncRNAs (differentially expressed lncRNAs [DE-lncRNAs]; log2-fold change ≥ 1, nominal p-value ≤ 0.05) in the seven examined tissues. The expressed lncRNAs (FPKM > 0.1) exhibited tissue-specific patterns in WT pigs. Specifically, 3410 lncRNAs were expressed exclusively in only one tissue. Potential functions of lncRNAs were indirectly predicted by searching their target cis- and trans-regulated protein-coding genes. LncRNAs with tissue-specific expression influence numerous genes related to tissue functions. Weighted gene coexpression network analysis (WGCNA) of 1570 DE-lncRNAs between WT and ZBED6 KO pigs was used to define the following six lncRNA modules specific to different tissues: skeletal muscle, heart, lung, spleen, kidney and liver modules. Furthermore, by conjoint analysis of longissimus dorsi data (tissue-specific expression, muscle module and DE-lncRNAs) and ChIP-PCR revealed NONSUSG002145.1 (adjusted p-values = 0.044), which is coexpressed with the IGF2 gene and binding with ZBED6, may play important roles in ZBED6 KO pig skeletal muscle development.

CONCLUSIONS

These findings indicate that the identified lncRNAs may play essential roles in tissue function and regulate the mechanism of ZBED6 action in skeletal muscle development in pigs. To our knowledge, this is the first study describing lncRNAs in ZBED6 KO pigs. These results may open new research directions leading to a better understanding of the global functions of ZBED6 and of lncRNA functions in skeletal muscle development in pigs.

Long Non-coding RNA Signatures Associated With Liver Aging in Senescence-Accelerated Mouse Prone 8 Model

The liver is sensitive to aging because the risk of hepatopathy, including fatty liver, hepatitis, fibrosis, cirrhosis, and hepatocellular carcinoma, increases dramatically with age. Long non-coding RNAs (lncRNAs) are >200 nucleotides long and affect many pathological and physiological processes. A potential link was recently discovered between lncRNAs and liver aging; however, comprehensive and systematic research on this topic is still limited. In this study, the mouse liver genome-wide lncRNA profiles of 8-month-old SAMP8 and SAMR1 models were explored through deep RNA sequencing. A total of 605,801,688 clean reads were generated. Among the 2,182 identified lncRNAs, 28 were differentially expressed between SAMP8 and SAMR1 mice. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) surveys showed that these substantially dysregulated lncRNAs participated in liver aging from different aspects, such as lipid catabolic (GO: 0016042) and metabolic pathways. Further assessment was conducted on lncRNAs that are most likely to be involved in liver aging and related diseases, such as LNC_000027, LNC_000204E, NSMUST00000144661.1, and ENSMUST00000181906.1 acted on Ces1g. This study provided the first comprehensive dissection of lncRNA landscape in SAMP8 mouse liver. These lncRNAs could be exploited as potential targets for the molecular-based diagnosis and therapy of age-related liver diseases.

20-Hydroxy-3-Oxolupan-28-Oic Acid, a Minor Component From Mahonia bealei (Fort.) Carr. Leaves Alleviates Lipopolysaccharide-Induced Inflammatory in Murine Macrophages

20-Hydroxyl-3-oxolupan-28-oic acid (HOA), a minor component from Mahonia bealei (Fort.) Carr. leaves, has been found to attenuate inflammatory responses. However, the underlying molecular mechanism is still unclear. In this study, we performed a comprehensive transcriptional study to investigate genetic changes. We used RNA sequencing technology to analyses the transcriptional changes in RAW 264.7 cells in a control group, lipopolysaccharide (LPS)-induced group, and HOA-treated group. We identified 1,313 and 388 differentially expressed genes (DEGs) in the control/LPS group and LPS/HOA group, respectively. Gene Ontology (GO) classification revealed that the DEGs were mainly enriched in a series of inflammatory and immune-related processes. The results of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that the DEGs were mainly enriched in inflammatory-related pathways such as the nuclear factor-kappa B (NF-κB) signaling pathway, cytokine-cytokine receptor interaction, chemokine signaling pathway, mitogen-activated protein kinase (MAPK) pathway, and Janus kinase-signal transducer and activator of transcription proteins (JAK-STAT) signaling pathway. The results of qPCR validation revealed that dynamic changes in immune-related mRNAs such as Saa3, Bcl2l1, Mapkapk2, Ccl9, Sdc4, Ddx3x, Socs3, Prdx5, Tlr4, Lif, IL15, Tnfaip3, Tet2, Tgf-β1, and Ccl20, which were significantly upregulated in the LPS group and downregulated in the HOA group. Taken together, these results suggest that HOA may be used as a source of anti-inflammatory agents as well as a dietary complement for health promotion.

Analysis of lncRNA, miRNA, and mRNA Expression Profiling in Type I IFN and Type II IFN Overexpressed in Porcine Alveolar Macrophages

Current data is scarce regarding the function of noncoding RNAs (ncRNAs) such as microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) in the interferon- (IFN-) mediated immune response. This is a comprehensive study that analyzes the lncRNA and miRNA expression profiles of the type I IFN and type II IFN in porcine alveolar macrophages using RNA sequencing. There was a total of 152 overexpressed differentially expressed (DE) lncRNAs and 21 DE miRNAs across type I IFN and type II IFN in porcine alveolar macrophages. Subsequent lncRNA-miRNA-mRNA network construction revealed the involvement of 36 DE lncRNAs and 12 DE miRNAs. LncRNAs such as the XLOC_211306, XLOC_100516, XLOC_00695, XLOC_149196, and XLOC_014459 were expressed at a higher degree in the type I IFN group, while XLOC_222640, XLOC_047290, XLOC_147777, XLOC_162298, XLOC_220210, and XLOC_165237 were expressed at a higher degree in the type II IFN group. These lncRNAs were found to act as "sponges" for miRNAs such as miR-34a, miR-328, miR-885-3p, miR-149, miR-30c-3p, miR-30b-5p, miR-708-5p, miR-193a-5p, miR-365-5p, and miR-7. Their target genes FADS2, RPS6KA1, PIM1, and NOD1 were found to be associated with several immune-related signaling pathways including the NOD-like receptor, Jak-STAT, mTOR, and PPAR signaling pathways. These experiments provide a comprehensive profile of overexpressed noncoding RNAs in porcine alveolar macrophages, providing new insights regarding the IFN-mediated immune response.

Profiling and integrated analysis of whole-transcriptome changes in uterine caruncles of pregnant and non-pregnant buffaloes

Improved reproductive performance in buffaloes can be achieved by understanding the basic mechanism governing the embryonic attachment and feto-maternal communication. Considering this, trascriptomic profiling and integrative analysis of long intergenic non-coding RNAs were carried out in the uterine caruncles of pregnant and non-pregnant buffaloes. Transcriptome data of pregnant and non-pregnant uterine caruncles after quality control was used to perform the analysis. Total of 86 novel lincRNAs expressed in uterine caruncular tissues were identified and characterized. Differential expression analysis revealed that 447 mRNAs and 185 mRNAs were up- and down- regulated, respectively. The number of up- and down- regulated lincRNAs were 114 and 13, respectively. Of the identified 86 novel lincRNAs, six novel lincRNAs were up-regulated in the pregnant uterine caruncles. GO terms (biological process) and PANTHER pathways associated with reproduction and embryogenesis were over-represented in differentially expressed genes. Through miRNA interaction analysis, interactions of 16 differentially expressed lincRNAs with mi-RNAs involved in reproduction were identified. This study has provided a catalogue of differentially expressed genes and novel regions previously unknown to play a significant role in buffalo reproduction. The results from the current study extends the buffalo uterine lncRNAs database and provides candidate regulators for future molecular genetic studies on buffalo uterine physiology to improve the embryo implantation and successful completion of pregnancy.

Profiling Non-Coding RNA Changes Associated with 16 Different Engineered Nanomaterials in a Mouse Airway Exposure Model

Perturbations in cellular molecular events and their associated biological processes provide opportunities for hazard assessment based on toxicogenomic profiling. Long non-coding RNAs (lncRNAs) are transcribed from DNA but are typically not translated into full-length proteins. Via epigenetic regulation, they play important roles in organismal response to environmental stress. The effects of nanoparticles on this important part of the epigenome are understudied. In this study, we investigated changes in lncRNA associated with hazardous inhalatory exposure of mice to 16 engineered nanomaterials (ENM)–4 ENM (copper oxide, multi-walled carbon nanotubes, spherical titanium dioxide, and rod-like titanium dioxide particles) with 4 different surface chemistries (pristine, COOH, NH₂, and PEG). Mice were exposed to 10 µg of ENM by oropharyngeal aspiration for 4 consecutive days, followed by cytological analyses and transcriptomic characterization of whole lung tissues. The number of significantly altered non-coding RNA transcripts, suggestive of their degrees of toxicity, was different for each ENM type. Particle surface chemistry and shape also had varying effects on lncRNA expression. NH₂ and PEG caused the strongest and weakest responses, respectively. Via correlational analyses to mRNA expression from the same samples, we could deduce that significantly altered lncRNAs are potential regulators of genes involved in mitotic cell division and DNA damage response. This study sheds more light on epigenetic mechanisms of ENM toxicity and also emphasizes the importance of the lncRNA superfamily as toxicogenomic markers of adverse ENM exposure.

Deep Insight Into Long Non-coding RNA and mRNA Transcriptome Profiling in HepG2 Cells Expressing Genotype IV Swine Hepatitis E Virus ORF3

Swine hepatitis E (swine HE) is a new type of zoonotic infectious disease caused by the swine hepatitis E virus (swine HEV). Open reading frame 3 (ORF3) is an important virulent protein of swine HEV, but its function still is mainly unclear. In this study, we generated adenoviruses ADV4-ORF3 and ADV4 negative control (ADV4-NC), which successfully mediated overexpression of enhanced green fluorescent protein (EGFP)-ORF3 and EGFP, respectively, in HepG2 cells. High-throughput sequencing was used to screen for differentially expressed long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs). The cis-target genes of lncRNAs were predicted, functional enrichment (Gene Ontology [GO] and Kyoto Encyclopedia of Genes and Genomes [KEGG]) was performed, and 12 lncRNAs with statistically significant different expressions (p ≤ 0.05 and q ≤ 1) were selected for further quantitative real-time reverse transcription (qRT-PCR) validation. In HepG2 cells, we identified 62 significantly differentially expressed genes (DEGs) (6,564 transcripts) and 319 lncRNAs (124 known lncRNAs and 195 novel lncRNAs) that were affected by ORF3, which were involved in systemic lupus erythematosus, Staphylococcus aureus infection, signaling pathways pluripotency regulation of stem cells, the peroxisome proliferator-activated receptor (PPAR) signaling pathway, and platinum drug resistance pathways. Cis-target gene prediction identified 45 lncRNAs corresponding to candidate mRNAs, among which eight were validated by qRT-PCR: LINC02476 (two transcripts), RAP2C-AS1, AC016526, AL139099, and ZNF337-AS1 (3 transcripts). Our results revealed that the lncRNA profile in host cells affected by ORF3, swine HEV ORF3, might affect the pentose and glucuronate interconversions and mediate the formation of obstructive jaundice by influencing bile secretion, which will help to determine the function of ORF3 and the infection mechanism and treatment of swine HE.

A systematic evaluation of bioinformatics tools for identification of long noncoding RNAs

High-throughput RNA sequencing unveiled the complexity of transcriptome and significantly increased the records of long noncoding RNAs (lncRNAs), which were reported to participate in a variety of biological processes. Identification of lncRNAs is a key step in lncRNA analysis, and a bunch of bioinformatics tools have been developed for this purpose in recent years. While these tools allow us to identify lncRNA more efficiently and accurately, they may produce inconsistent results, making selection a confusing issue. We compared the performance of 41 analysis models based on 14 software packages and different data sets, including high-quality data and low-quality data from 33 species. In addition, computational efficiency, robustness, and joint prediction of the models were explored. As a practical guidance, key points for lncRNA identification under different situations were summarized. In this investigation, no one of these models could be superior to others under all test conditions. The performance of a model relied to a great extent on the source of transcripts and the quality of assemblies. As general references, FEELnc_all_cl, CPC, and CPAT_mouse work well in most species while COME, CNCI, and lncScore are good choices for model organisms. Since these tools are sensitive to different factors such as the species involved and the quality of assembly, researchers must carefully select the appropriate tool based on the actual data. Alternatively, our test suggests that joint prediction could behave better than any single model if proper models were chosen. All scripts/data used in this research can be accessed at http://bioinfo.ihb.ac.cn/elit.

Molecular mechanism of Chinese alligator (Alligator sinensis) adapting to hibernation

Hibernation is a physiological state for Chinese alligators to cope with cold weather. In mammals, gene expression changes during hibernation and their regulatory mechanisms have been extensively studied, however, these studies in reptiles are still rare. Here, integrated analysis of messenger RNA (mRNA), microRNA (miRNA), and long noncoding RNA (lncRNA) reveals the molecular mechanisms of the hypothalamus, liver, and skeletal muscle in hibernating and active individuals. During hibernation, the number of genes increased in the hypothalamus, liver, and skeletal muscle was 585, 282, and 297, while the number of genes decreased was 215, 561, and 627, respectively, as compared with active individuals. Through Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis, the differential expressed genes were mainly enriched in DNA damage repair, biological rhythm, energy metabolism, myoprotein degradation, and other related items and pathways. Besides, 4740 miRNAs were identified in three tissues. Through the comprehensive analysis of miRNA and mRNA abundance profiles, 12,291, 6997, and 8232 miRNA-mRNA pairs all showed a negative correlation in the hypothalamus, liver, and skeletal muscle, respectively. Some miRNA target genes were related tobiological rhythm and energy metabolism, suggesting that miRNA may play an important role in the physiological metabolism of the hibernating adaptability of Chinese alligators. Moreover, 402, 230, and 130 differentially expressed lncRNAs were identified in the hypothalamus, liver, and skeletal muscle, respectively. The targeting relationship of four lncRNA-mRNA pairs were predicted, with the main function of target genes involved in the amino acid transportation. These results are helpful to further understand the molecular regulatory basis of the hibernation adaptation in Chinese alligators.

Analysis of long non-coding RNAs in epididymis of cattleyak associated with male infertility

Cattleyak (CY), is a cross breed between cattle and yak (YK), which display equal adaptability to the harsh environment as YK and much higher performances than YK. However, the CY is female fertile and male sterile. Previous studies were conducted on testes tissues to investigate the mechanism of male infertility in CY. There is no systematic research on genes, especially lncRNAs between CY and YK epididymis. In this study, Illumina Hiseq was performed to profile the epididymis transcriptome (lncRNA and mRNA) of CY and YK. In total 18859 lncRNAs were identified, from which lincRNAs 12458, antisense lncRNAs 2345, intronic lncRNAs 3101, and sense lncRNAs 955 respectively. We have identified 345 DE lncRNAs and 3008 DE mRNAs between YK and CY epididymis. Thirteen DEGs were validated by quantitative real-time PCR. Combing with DEG, 14 couples of lncRNAs and their target genes were both DE, and 6 of them including CCDC39, KCNJ16, NECTIN2, MRPL20, PSMC4, and DEFB112 show their potential infertility-related terms such as cellular motility, sperm maturation, sperm storage, cellular junction, folate metabolism, and capacitation. On the other hand, several down-regulated genes such as DEFB124, DEFB126, DEFB125, DEFB127, DEFB129, CES5A, TKDP1, CST3, RNASE9 and CD52 in CY compared to YK were involved in the immune response and sperm maturation. Therefore, comprehensive analysis for lncRNAs and their target genes may enhance our understanding of the molecular mechanisms underlying the process of sperm maturation in CY and may provide important resources for further research.

Analysis and Screening of Reproductive Long Non-coding RNAs Through Genome-Wide Analyses of Goat Endometrium During the Pre-attachment Phase

Reproduction in goat is highly impeded by implantation failure. Of concern, the underlying mechanism leading to embryo implantation remains unclear. In this study, deep sequencing was employed through strand-specific Ribo-Zero RNA-Seq to characterize transcriptome changes in the endometrium during the maternal recognition of pregnancy. A total of 996 differential transcripts (115 lncRNAs and 881 mRNAs) existing between the pregnant and non-pregnant endometrium were revealed through bioinformatics analysis. The screening was performed on lncRNAs (XR_001918173.1, LNC_002760, and LNC_000599) and LNC_009053, to determine their potential role in regulating the synthesis of retinol and endometrium remolding through the proteasome pathway, respectively. The hypothesis of whether certain lncRNAs, namely, LNC_007223, LNC_005256, and LNC_010092 could play important roles in embryo implantation was tested. These novel findings are of paramount relevance to further elucidate the molecular mechanisms of embryo implantation and uncover new targets to improve goat reproduction.

Genome-Wide Identification of Long Non-coding RNAs in the Gravid Ectoparasite Varroa destructor

Long non-coding RNAs (lncRNAs) emerge as critical regulators across a wide variety of biological functions in living organisms. However, to date, no systematic characterization of lncRNAs has been investigated in the ectoparasitic mite Varroa destructor, the most severe biotic threat to honey bees worldwide. Here, we performed an initial genome-wide identification of lncRNAs in V. destructor via high-throughput sequencing technology and reported, for the first time, the transcriptomic landscape of lncRNAs in the devastating parasite. By means of a lncRNA identification pipeline, 6,645 novel lncRNA transcripts, encoded by 3,897 gene loci, were identified, including 2,066 sense lncRNAs, 2,772 lincRNAs, and 1,807 lncNATs. Compared with protein-coding mRNAs, V. destructor lncRNAs are shorter in terms of full length, as well as of the ORF length, contain less exons, and express at lower level. GO term and KEGG pathway enrichment analyses of the lncRNA target genes demonstrated that these predicted lncRNAs may be potentially responsible for the regulatory functions of cellular and biological progresses in the reproductive phase of V. destructor. To our knowledge, this is the first catalog of lncRNA profile in the parasitiformes species, providing a valuable resource for genetic and genomic studies. Understanding the characteristics and features of lncRNAs in V. destructor would promote sustainable parasite control.

Integrated analysis of lncRNA, miRNA and mRNA reveals novel insights into the fertility regulation of large white sows

Background

Improving sow fertility is extremely important as it can lead to increased reproductive efficiency and thus profitability for swine producers. There are considerable differences in fertility rates among individual animals, but the underlying molecular mechanisms remain unclear. In this study, by using different types of RNA libraries, we investigated the complete transcriptome of ovarian tissue during the luteal (L) and follicular (F) phases of the estrous cycle in Large White pigs with high (H) and low (L) fecundity, and performed a comprehensive analysis of long noncoding RNAs (lncRNAs), mRNAs and micro RNAs (miRNAs) from 16 samples by combining RNA sequencing (RNA-seq) with bioinformatics.

Results

In total, 24,447 lncRNAs, 27,370 mRNAs, and 216 known miRNAs were identified in ovarian tissues. The genomic features of lncRNAs, such as length distribution and number of exons, were further analyzed. We selected a threshold of P < 0.05 and |log₂ (fold change)| ≥ 1 to obtain the differentially expressed lncRNAs, miRNAs and mRNAs by pairwise comparison (LH vs. LL, FH vs. FL). Bioinformatics analysis of these differentially expressed RNAs revealed multiple significantly enriched pathways (P < 0.05) that were closely involved in the reproductive process, such as ovarian steroidogenesis, lysosome, steroid biosynthesis, and the estrogen and GnRH signaling pathways. Moreover, bioinformatics screening of differentially expressed miRNAs that share common miRNA response elements (MREs) with lncRNAs and their downstream mRNA targets were performed. Finally, we constructed lncRNA–miRNA–mRNA regulation networks. The key genes in these networks were verified by Reverse Transcription Real-time Quantitative PCR (RT-qRCR), which were consistent with the results from RNA-Seq data.

Conclusions

These results provide further insights into the fertility of pigs andcan contribute to further experimental investigation of the functions of these genes.

Long non-coding RNA366.2 controls endometrial epithelial cell proliferation and migration by upregulating WNT6 as a ceRNA of miR-1576 in sheep uterus

Long non-coding RNAs (lncRNAs) play an important regulatory role in mammalian fecundity. Currently, most studies are primarily concentrated on ovarian lncRNAs, ignoring the influence of uterine lncRNAs on the fecundity of female sheep. In this study, we found a higher density of uterine glands and endometrial microvessel density (MVD) in high prolificacy group of Hu sheep compared to low prolificacy groups (p < 0.05) as well as an increased level of serum placental growth factor (PLGF). Hundreds of differentially expressed (DE) lncRNAs were identified in Hu sheep with different fecundity by RNA sequencing (RNA-seq), and their targets were enriched in some signaling pathways involved in endometrial functions, such as the estrogen signaling pathway, nuclear factor kappa B (NF-κB) signaling pathway, oxytocin signaling pathway, and Wnt signaling pathway. Furthermore, the underlying mechanisms of competitive endogenous RNA (ceRNA) of lncRNA366.2-miR-1576- WNT6 were determined by bioinformatics analysis. Functionally, our results indicated that lncRNA366.2 promoted endometrial epithelial cell (EEC) proliferation, migration, and growth factor expression by sponging miR-1576 to upregulate WNT6 expression and activate the Wnt/β-catenin pathway. Taken together, our research indicated the regulatory mechanism of the lncRNA366.2-miR-1576-WNT6 in EEC proliferation and migration. Furthermore, this study provides a new theoretical reference for the identification of candidate genes related to fecundity.

Identification of Differentially Expressed Gene Transcripts in Porcine Endometrium during Early Stages of Pregnancy

During the early stages of pregnancy, the uterine endometrium undergoes dramatic morphologic and functional changes accompanied with dynamic variation in gene expression. Pregnancy-stage specific differentially expressed gene (DEG)-transcript-probes were investigated and identified by comparing endometrium transcriptome at 9th day (9D), 12th day (12D) and 16th day (16D) of early pregnancy in Polish large-white (PLW) gilts. Endometrium comparisons between 9D-vs-12D, 9D-vs-16D and 12D-vs-16D of early pregnancy identified 6049, 374 and 6034 highly significant DEG-transcript-probes (p < 0.001; >2 FC). GO term enrichment analysis identified commonly shared upregulated endometrial DEG-transcript-probes (p < 0.001; >2 FC), that were regulating the gene functions of anatomic structure development and transport (TG), DNA-binding and methyltransferase activity (ZBTB2), ion-binding and kinase activity (CKM), cell proliferation and apoptosis activity (IL1B). Downregulated DEG-transcript-probes (p < 0.001; >2 FC) were involved in regulating the gene functions of phosphatase activity (PTPN11), TC616413 gene-transcript and Sus-scrofa LOC100525539. Moreover, blastn comparison of microarray-probes sequences against sus-scrofa11 assembly identified commonly shared upregulated endometrial DEG-transcript-probes (E < 0.06; >2 FC), that were regulating the gene functions of reproduction and growth (SELENOP), cytoskeleton organization and kinase activity (CDC42BPA), phosphatase activity (MINPP1), enzyme-binding and cell-population proliferation (VAV3), cancer-susceptibility candidate gene (CASC4), cytoskeletal protein-binding (COBLL1), ion-binding, enzyme regulator activity (ACAP2) Downregulated endometrial DEG-transcript-probes (E < 0.06; >2FC) were involved in regulating the gene functions of signal-transduction (TMEM33), catabolic and metabolic processes (KLHL15). Microarray validation experiment on selected candidate genes showed complementarity to significant endometrial DEG-transcript-probes responsible for the regulation of immune response (IL1B, S100A11), lipid metabolism (FABP3, PPARG), cell-adhesion (ITGAV), angiogenesis (IL1B), intercellular transmission (NMB), cell-adhesion (OPN) and response to stimuli (RBP4) was confirmed by RT-PCR. This study provides a clue that identified pregnancy-stage specific microarray transcript probes could be considered as candidate genes for recognition and establishment of early pregnancy in the pig.

Review: Long non-coding RNA in livestock

Less than 2% of mammalian genomes code for proteins, but 'the majority of its bases can be found in primary transcripts' - a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species - Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals' genomes is critical in bridging the gap between genotype and phenotype.

The Interaction of lncRNA XLOC-2222497, AKR1C1, and Progesterone in Porcine Endometrium and Pregnancy

The endometrium is an important tissue for pregnancy and plays an important role in reproduction. In this study, high-throughput transcriptome sequencing was performed in endometrium samples of Meishan and Yorkshire pigs on days 18 and 32 of pregnancy. Aldo-keto reductase family 1 member C1 (AKR1C1) was found to be a differentially expressed gene, and was identified by quantitative real-time PCR (qRT-PCR) and Western blot. Immunohistochemistry results revealed the cellular localization of the AKR1C1 protein in the endometrium. Luciferase activity assay demonstrated that the AKR1C1 core promoter region was located in the region from -706 to -564, containing two nuclear factor erythroid 2-related factor 2 (NRF2) binding sites (antioxidant response elements, AREs). XLOC-2222497 was identified as a nuclear long non-coding RNA (lncRNA) highly expressed in the endometrium. XLOC-2222497 overexpression and knockdown have an effect on the expression of AKR1C1. Endocrinologic measurement showed the difference in progesterone levels between Meishan and Yorkshire pigs. Progesterone treatment upregulated AKR1C1 and XLOC-2222497 expression in porcine endometrial epithelial cells. In conclusion, transcriptome analysis revealed differentially expressed transcripts during the early pregnancy process. Further experiments demonstrated the interaction of XLOC-2222497/AKR1C1/progesterone in the endometrium and provided new potential targets for pregnancy maintenance and its control.

Long noncoding RNA and mRNA profiling of hypothalamic-pituitary-mammary gland axis in lactating sows under heat stress

Heat stress (HS) seriously affects sow lactation performance and Long non-coding RNAs (lncRNAs) play vital roles in the regulation of transcription and post transcription. However, the mechanism of lncRNAs expression affecting lactation performance on the hypothalamus-pituitary-mammary axis of sows is still unclear. In this study, we performed RNA sequencing and bioinformatics analysis of the hypothalamus, pituitary, and mammary gland tissues of lactating sows under HS and thermal comfort. In total, the analysis identified 658, 6021, and 6745 differently expressed (DE) mRNAs, 26, 126, and 169 DE lncRNAs between comparison groups in the hypothalamus, pituitary, and mammary glands, respectively. The hormone genes and most DE mRNAs encoding heat shock protein were differently expressed in the HS group. In addition, 2, 60, and 86 pairs of DE lncRNAs and mRNAs correlation were observed in those tissues, respectively. Some lncRNAs may be involved in the regulation of lactation performance in the HS sows.

Exploring the roles of fecundity-related long non-coding RNAs and mRNAs in the adrenal glands of small-tailed Han Sheep

BACKGROUND

Long non-coding RNAs (lncRNAs) can play important roles in uterine and ovarian functions. However, little researches have been done on the role of lncRNAs in the adrenal gland of sheep. Herein, RNA sequencing was used to compare and analyze gene expressions in adrenal tissues between follicular phases and luteal phases in FecBBB (MM) and FecB++ (WW) sheep, respectively, and differentially expressed lncRNAs and genes associated with reproduction were identified.

RESULTS

In MM sheep, 38 lncRNAs and 545 mRNAs were differentially expressed in the adrenal gland between the luteal and follicular phases; In WW sheep, 513 differentially expressed lncRNAs and 2481 mRNAs were identified. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses indicated that differentially expressed lncRNAs and their target genes are mainly involved in the circadian rhythm, the mitogen activated protein kinase, thyroid, ovarian steroidogenesis and transforming growth factor beta signaling pathways. Differentially expressed lncRNAs can regulate reproduction by modulating genes involved in these signaling pathways and biological processes. Specifically, XLOC_254761, XLOC_357966, 105,614,839 and XLOC_212877 targeting CREB1, PER3, SMAD1 and TGFBR2, respectively, appear to play key regulatory roles.

CONCLUSION

These results broaden our understanding of lncRNAs in adrenal gland of sheep and provide new insights into the molecular mechanisms underlying sheep reproduction.

The regulation mechanism of lncRNAs and mRNAs in sea cucumbers under global climate changes: Defense against thermal and hypoxic stresses

The aquatic environment can be greatly impacted by thermal and hypoxic stresses, particularly caused by intensified global warming. Hence, there is an urgency to understand the response mechanisms of marine organisms to adverse environment. Although long non-coding RNAs (lncRNAs) are involved in many biological processes, their roles in stress responses still remain unclear. Here, differentially expressed (DE) lncRNAs and mRNAs were identified as responses to environmental stresses in the economically important sea cucumber, Apostichopus japonicus, and their potential roles were explored. Based on a total of 159, 355 and 495 significantly upregulated genes and 230, 518 and 647 significantly downregulated genes identified in the thermal, hypoxic and combination thermal + hypoxic stress treatments, respectively, we constructed DE-lncRNA-mRNA coexpression networks. Among the networks, eight shared pairs were identified from the three treatments, and based on the connectivity degree, MSTRG.27265, MSTRG.19729 and MSTRG.95524 were shown to be crucial lncRNAs. Among all the significantly changed lncRNAs identified by RT-qPCR and sequencing data, binding sites were found in four other lncRNAs (MSTRG.34610, MSTRG.10941, MSTRG.81281 and MSTRG.93731) with Aja-miR-2013-3p, a key miRNA that responds to hypoxia in sea cucumbers. The hypoxia-inducible factor (HIF-1α) was also shown as the possible targeted mRNA of Aja-miR-2013-3p. As indicated by a dual-luciferase reporter assay system, "HIF-1α gene/Aja-miR-2013-3p/MSTRG.34610" network and the "HIF-1α gene/Aja-miR-2013-3p/MSTRG.10941" network may play important roles in sea cucumbers under environmental stresses. Moreover, environmental stress altered the expression of multiple lncRNAs and mRNAs, thus affecting various biological processes in A. japonicus, including immunity, energy metabolism and the cell cycle. At the molecular level, more comprehensive responses were elicited by the combined thermal/hypoxic stress treatment than by individual stresses alone in sea cucumbers. This study lays the groundwork for future research on molecular mechanisms of echinoderm responses to thermal and hypoxic stress in the context of global climate changes.

Growth hormone promotes human endometrial glandular cells proliferation and motion through the GHR-STAT3/5 pathway

Background

This study aims at investigating the effect of growth hormone (GH) on the growth of human endometrial glandular cells (hEGCs) and preliminary exploring its mechanism.

Methods

HEGCs were isolated from the endometrial biopsies and exposed to different dose of GH (0, 50, 100, and 200 ng/mL). Cell proliferation and cell cycle assay, migration assay was performed to investigate the growth and motivation of hEGCs, respectively. Reverse transcription-polymerase chain reaction (RT-PCR), immunocytochemistry (ICC), and western blot (WB) were processed to investigate its related gene or protein expression.

Results

The results revealed that GH administration promoted the proliferation, cell cycle, migration, and growth hormone receptors (GHRs) expression of the hEGC. We further inhibited GHRs with AG490, and the inhibitor reversed the effects of GH on cell growth, motion, and the activation of GHR and STAT3/5.

Conclusions

GH promoted hEGCs proliferation and motion, which is GHR-JAK-STAT3/5 signaling pathway-dependent. These findings reveal the essential roles of GH in the hEGCs growth and provide evidence for potential GH therapy in intrauterine adhesion (IUA) treatment.

Revealing stage-specific expression patterns of long noncoding RNAs along mouse spermatogenesis

The discovery of a large number of long noncoding RNAs (lncRNAs), and the finding that they may play key roles in different biological processes, have started to provide a new perspective in the understanding of gene regulation. It has been shown that the testes express the highest amount of lncRNAs among different vertebrate tissues. However, although some studies have addressed the characterization of lncRNAs along spermatogenesis, an exhaustive analysis of the differential expression of lncRNAs at its different stages is still lacking. Here, we present the results for lncRNA transcriptome profiling along mouse spermatogenesis, employing highly pure flow sorted spermatogenic stage-specific cell populations, strand-specific RNAseq, and a combination of up-to-date bioinformatic pipelines for analysis. We found that the vast majority of testicular lncRNA genes are expressed at post-meiotic stages (i.e. spermiogenesis), which are characterized by extensive post-transcriptional regulation. LncRNAs at different spermatogenic stages shared common traits in terms of transcript length, exon number, and biotypes. Most lncRNAs were lincRNAs, followed by a high representation of antisense (AS) lncRNAs. Co-expression analyses showed a high correlation along the different spermatogenic stage transitions between the expression patterns of AS lncRNAs and their overlapping protein-coding genes, raising possible clues about lncRNA-related regulatory mechanisms. Interestingly, we observed the co-localization of an AS lncRNA and its host sense mRNA in the chromatoid body, a round spermatids-specific organelle that has been proposed as a reservoir of RNA-related regulatory machinery. An additional, intriguing observation is the almost complete lack of detectable expression for Y-linked testicular lncRNAs, despite that a high number of lncRNA genes are annotated for this chromosome.

Neonatal lactocrine deficiency affects the adult porcine endometrial transcriptome at pregnancy day 13

Reproductive performance of female pigs that do not receive sufficient colostrum from birth is permanently impaired. Whether lactocrine deficiency, reflected by low serum immunoglobulin immunocrit (iCrit), affects patterns of endometrial gene expression during the periattachment period of early pregnancy is unknown. Here, objectives were to determine effects of low iCrit at birth on the adult endometrial transcriptome on pregnancy day (PxD) 13. On the first day of postnatal life, gilts were assigned to high or low iCrit groups. Adult high (n = 8) and low (n = 7) iCrit gilts were bred (PxD 0), and humanely slaughtered on PxD 13 when tissues and fluids were collected. The endometrial transcriptome was defined for each group using mRNAseq and microRNAseq. Reads were mapped to the Sus scrofa 11.1 genome build. Mature microRNAs were annotated using miRBase 21. Differential expression was defined based on fold change (≥ ±1.5). Lactocrine deficiency did not affect corpora lutea number, uterine horn length, uterine wet weight, conceptus recovery, or uterine luminal fluid estrogen content on PxD 13. However, mRNAseq revealed 1157 differentially expressed endometrial mRNAs in high versus low iCrit gilts. Differentially expressed genes had functions related to solute transport, endometrial receptivity, and immune response. Six differentially expressed endometrial microRNAs included five predicted to target 62 differentially expressed mRNAs, affecting similar biological processes. Thus, lactocrine deficiency on the first day of postnatal life can alter uterine developmental trajectory with lasting effects on endometrial responses to pregnancy as reflected at the level of the transcriptome on PxD 13.

Identification, functional prediction, and key lncRNA verification of cold stress-related lncRNAs in rats liver

Cold stimulation reduces the quality of animal products and increases animal mortality, causing huge losses to the livestock industry in cold regions. Long non-coding RNAs (lncRNAs) take part in many biological processes through transcriptional regulation, intracellular material transport, and chromosome remodeling. Although cold stress-related lncRNAs have been reported in plants, no research is available on the characteristic and functional analysis of lncRNAs after cold stress in rats. Here, we built a cold stress animal model firstly. Six SPF male Wistar rats were randomly divided to the acute cold stress group (4 °C, 12 h) and the normal group (24 °C, 12 h). lncRNA libraries were constructed by high-throughput sequencing (HTS) using rat livers. 2,120 new lncRNAs and 273 differentially expressed (DE) lncRNAs were identified in low temperature environments. The target genes of DElncRNA were predicted by cis and trans, and then functional and pathway analysis were performed to them. GO and KEGG analysis revealed that lncRNA targets were mainly participated in the regulation of nucleic acid binding, cold stimulation reaction, metabolic process, immune system processes, PI3K-Akt signaling pathway and pathways in cancer. Next, a interaction network between lncRNA and its targets was constructed. To further reveal the mechanism of cold stress, DElncRNA and DEmRNA were extracted to reconstruct a co-expression sub-network. We found the key lncRNA MSTRG.80946.2 in sub-network. Functional analysis of key lncRNA targets showed that targets were significantly enriched in fatty acid metabolism, the PI3K-Akt signaling pathway and pathways in cancer under cold stress. qRT-PCR confirmed the sequencing results. Finally, hub lncRNA MSTRG.80946.2 was characterized, and verified its relationship with related mRNAs by antisense oligonucleotide (ASO) interference and qRT-PCR. Results confirmed the accuracy of our analysis. To sum up, our work was the first to perform detailed characterization and functional analysis of cold stress-related lncRNAs in rats liver. lncRNAs played crucial roles in energy metabolism, growth and development, immunity and reproductive performance in cold stressed rats. The MSTRG.80946.2 was verified by network and experiments to be a key functional lncRNA under cold stress, regulating ACP1, TSPY1 and Tsn.