Genome-wide analysis of long non-coding RNAs at early stage of skin pigmentation in goats (Capra hircus)

Single-Cell Transcriptome Sequence Profiling on the Morphogenesis of Secondary Hair Follicles in Ordos Fine-Wool Sheep

The Ordos fine-wool sheep is a high-quality breed in China that produces superior natural textiles and raw materials such as wool and lamb meat. However, compared to the Australian Merino sheep, there is still a gap in terms of the wool fiber fineness and wool yield. The hair follicle is the main organ that controls the type of wool fiber, and the morphological changes in the secondary hair follicle are crucial in determining wool quality. However, the process and molecular mechanisms of hair follicle morphogenesis in Ordos fine-wool sheep are not yet clear. Therefore, analyzing the molecular mechanisms underlying the process of follicle formation is of great significance for improving the fiber diameter and wool production of Ordos fine-wool sheep. The differential expressed genes, APOD, POSTN, KRT5, and KRT15, which related to primary hair follicles and secondary hair follicles, were extracted from the dermal papillae. Based on pseudo-time analysis, the differentiation trajectories of dermal lineage cells and epidermal lineage cells in the Ordos fine-wool sheep were successfully constructed, providing a theoretical basis for breeding research in Ordos fine-wool sheep.

Transcriptome analysis reveals core lncRNA-mRNA networks regulating melanization and biomineralization in Patinopecten yessoensis shell-infested by Polydora

BACKGROUND

Patinopecten yessoensis, a large and old molluscan group, has been one of the most important aquaculture shellfish in Asian countries because of its high economic value. However, the aquaculture of the species has recently been seriously affected by the frequent outbreaks of Polydora disease, causing great economic losses. Long non-coding RNAs (lncRNAs) exhibit exhibit crucial effects on diverse biological processes, but still remain poorly studied in scallops, limiting our understanding of the molecular regulatory mechanism of P. yessoensis in response to Polydora infestation.

RESULTS

In this study, a high-throughput transcriptome analysis was conducted in the mantles of healthy and Polydora-infected P. yessoensis by RNA sequencing. A total of 19,133 lncRNAs with 2,203 known and 16,930 novel were identified. The genomic characterizations of lncRNAs showed shorter sequence and open reading frame (ORF) length, fewer number of exons and lower expression levels in comparison with mRNAs. There were separately 2280 and 1636 differentially expressed mRNAs and lncRNAs (DEGs and DELs) detected in diseased individuals. The target genes of DELs were determined by both co-location and co-expression analyses. Functional enrichment analysis revealed that DEGs involved in melanization and biomineralization were significantly upregulated; further, obviously increased melanin granules were observed in epithelial cells of the edge mantle in diseased scallops by histological and TEM study, indicating the crucial role of melanizaiton and biomineralization in P. yessoensis to resist against Polydora infestation. Moreover, many key genes, such as Tyrs, Frizzled, Wnts, calmodulins, Pifs, perlucin, laccase, shell matrix protein, mucins and chitins, were targeted by DELs. Finally, a core lncRNA-mRNA interactive network involved in melanization and biomineralization was constructed and validated by qRT-PCR.

CONCLUSIONS

This work provides valuable resources for studies of lncRNAs in scallops, and adds a new insight into the molecular regulatory mechanisms of P. yessoensis defending against Polydora infestation, which will contribute to Polydora disease control and breeding of disease-resistant varieties in molluscs.

Protective Effects of ζ-Carotene-like Compounds against Acute UVB-Induced Skin Damage

The previous study successfully established an expression strain of ζ-carotene-like compounds (CLC) and demonstrated its remarkable antioxidant activity, which exhibited resistance to photodamage caused by UVB radiation on the skin following gavage administration. The objective of this study was to investigate the impact and mechanism of CLC on UVB-induced skin damage through topical application. Cell viability, anti-apoptotic activity, ROS scavenging ability, the inhibition of melanin synthesis, the regulation of inflammatory factors and collagen deposition were assessed in cells and mice using qRT-PCR, WB, Elisa assays, immunohistochemistry staining and biochemical kits, etc. The experimental results demonstrated that CLC-mitigated apoptosis induced by UVB irradiation up-regulated the Keap1/Nrf2/ARE antioxidant pathway to attenuate levels of ROS and inflammatory factors (NF-κB, TNF-α, IL-6 and IL-β), and suppressed MAPK/AP-1 and CAMP/PKA/CREB signaling pathways to mitigate collagen degradation, skin aging and melanin formation. In conclusion, this study underscored the potential of CLC as a safe and efficacious source of antioxidants, positioning it as a promising ingredient in the formulation of cosmetics targeting anti-aging, skin brightening and sunburn repair.

Comprehensive transcriptional analysis of pig facial skin development

Background

Skin development is a complex process that is influenced by many factors. Pig skin is used as an ideal material for xenografts because it is more anatomically and physiologically similar to human skin. It has been shown that the skin development of different pig breeds is different, and some Chinese pig breeds have the characteristics of skin thickness and facial skin folds, but the specific regulatory mechanism of this skin development is not yet clear.

Methods

In this study, the facial skin of Chenghua sows in the four developmental stages of postnatal Day 3 (D3) , Day 90 (D90) , Day 180 (D180), and Year 3 (Y3) were used as experimental materials, and RNA sequencing (RNA-seq) analysis was used to explore the changes in RNA expression in skin development at the four developmental stages, determine the differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), microRNAs (miRNAs), and circular RNAs (circRNAs), and perform functional analysis of related genes by Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses.

Results

A pairwise comparison of the four developmental stages identified several differentially expressed genes (DEGs) and found that the number of differentially expressed RNAs (DE RNAs) increased with increasing developmental time intervals. Elastin (ELN) is an important component of the skin. Its content affects the relaxation of the epidermis and dermal connection, and its expression is continuously downregulated during the four developmental stages. The functions of DEGs at different developmental stages were examined by performing GO and KEGG analyses, and the GO terms and enrichment pathways of mRNAs, lncRNAs, miRNAs, and circRNAs highly overlapped, among which the PPAR signaling pathway, a classical pathway for skin development, was enriched by DEGs of D3 vs. D180, D90 vs. D180 and D180 vs. Y3. In addition, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA interaction networks and found genes that may be associated with skin development, but their interactions need further study.

Conclusions

We identified a number of genes associated with skin development, performed functional analyses on some important DEGs and constructed interaction networks that facilitate further studies of skin development.

The Characteristic Function of Blood-Derived Exosomes and Exosomal circRNAs Isolated from Dairy Cattle during the Dry Period and Mid-Lactation

Exosomes are key mediators of intercellular communication. They are secreted by most cells and contain a cargo of protein-coding genes, long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs), which modulate recipient cell behavior. Herein, we collected blood samples from Holstein cows at days 30 (mid-lactation) and 250 (dry period) of pregnancy. Prolactin, follicle-stimulating hormone, luteinizing hormone, estrogen, and progesterone levels showed an obvious increase during D250. We then extracted exosomes from bovine blood samples and found that their sizes generally ranged from 100 to 200 nm. Further, Western blotting validated that they contained CD9, CD63, and TSG101, but not calnexin. Blood-derived exosomes significantly promoted the proliferation of mammary epithelial cells, particularly from D250. This change was accompanied by increased expression levels of proliferation marker proteins PCNA, cyclin D, and cyclin E, as detected by EdU assay, cell counting kit-8 assay, and flow cytometric cell cycle analysis. Moreover, we treated mammary epithelial cells with blood-derived exosomes that were isolated from the D30 and D250 periods. And RNA-seq of two groups of cells led to the identification of 839 differentially expressed genes that were significantly enriched in KEGG signaling pathways associated with apoptosis, cell cycle and proliferation. In bovine blood-derived exosomes, we found 12,747 protein-coding genes, 31,181 lncRNAs, 9374 transcripts of uncertain coding potential (TUCP) candidates, and 460 circRNAs, and 32 protein-coding genes, 806 lncRNAs, 515 TUCP candidates, and 45 circRNAs that were differentially expressed between the D30 and D250 groups. We selected six highly expressed and four differentially expressed circRNAs to verify their head-to-tail splicing using PCR and Sanger sequencing. To summarize, our findings improve our understanding of the key roles of blood-derived exosomes and the characterization of exosomal circRNAs in mammary gland development.

Identifying Candidate Genes for Litter Size and Three Morphological Traits in Youzhou Dark Goats Based on Genome-Wide SNP Markers

This study aimed to reveal the potential genetic basis for litter size, coat colour, black middorsal stripe and skin colour by combining genome-wide association analysis (GWAS) and selection signature analysis and ROH detection within the Youzhou dark (YZD) goat population (n = 206) using the Illumina GoatSNP54 BeadChip. In the GWAS, we identified one SNP (snp54094-scaffold824-899720) on chromosome 11 for litter size, two SNPs on chromosome 26 (snp11508-scaffold142-1990450, SORCS3) and chromosome 12 (snp55048-scaffold842-324525, LOC102187779) for coat colour and one SNP on chromosome 18 (snp56013-scaffold873-22716, TCF25) for the black middorsal stripe. In contrast, no SNPs were identified for skin colour. In selection signature analysis, 295 significant iHS genomic regions with a mean |iHS| score > 2.66, containing selection signatures encompassing 232 candidate genes were detected. In particular, 43 GO terms and one KEGG pathway were significantly enriched in the selected genes, which may contribute to the excellent environmental adaptability and characteristic trait formation during the domestication of YZD goats. In ROH detection, we identified 4446 ROH segments and 282 consensus ROH regions, among which nine common genes overlapped with those detected using the iHS method. Some known candidate genes for economic traits such as reproduction (TSHR, ANGPT4, CENPF, PIBF1, DACH1, DIS3, CHST1, COL4A1, PRKD1 and DNMT3B) and development and growth (TNPO2, IFT80, UCP2, UCP3, GHRHR, SIM1, CCM2L, CTNNA3 and CTNNA1) were revealed by iHS and ROH detection. Overall, this study is limited by the small population size, which affects the results of GWAS to a certain extent. Nevertheless, our findings could provide the first overview of the genetic mechanism underlying these important traits and provide novel insights into the future conservation and utilisation of Chinese goat germplasm resources.

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.

Transcriptome and anatomical studies reveal alterations in leaf thickness under long-term drought stress in tobacco

Drought is one of the foremost environmental factors that limit the growth of plants. Leaf thickness (LT) is an important quantitative trait in plant physiology. The experiment was carried out in a growth room and the plants were divided into two groups such as well-watered and drought-stressed. This work investigated leaf growth in terms of leaf surface growth and expansion rate, leaf stomata traits, LT, anticlinal growth, and leaf cell layers. The results showed that the leaf area and leaf surface expansion rate were decreased by drought stress (DS). Similarly, LT, anticlinal expansion rate, palisade and spongy tissue thickness, and their related expansion rates were also decreased at different days' time points (DTP) of DS. However, a steady increase was observed in the aforementioned parameters after 12 DTP of DS. The stomatal density increased while stomata size decreased at 3 DTP and 12 DTP (low leaf water potential and relative leaf water content at these time points) and vice versa at 24 DTP compared with the well-watered plants indicating adaptations in these traits in response to DS, and thus the leaf water status played a role in the regulation of leaf stomata traits. The cell length decreased in the upper epidermis, palisade and spongy tissues by DS up to 12 DTP led to lower LT while an increase was observed after 12 DTP that resulted in higher LT. The increase in the LT was supported by the upregulation of starch and sucrose metabolism, glycerolipid metabolism, protein processing in endoplasmic reticulum pathways at 18 DTP along with the differentially expressed genes induced that were related to cell wall remodeling (cellulose, expansin, xyloglucans) and cell expansion (auxin response factors and aquaporin). The results explain the response of leaf thickness to drought stress and show alterations in LT and leaf stomatal traits. This study might serve as a valuable source of gene information for functional studies and provide a theoretical basis to understand leaf growth in terms of leaf anatomy and leaf stomatal traits under drought stress.

Grade follicles transcriptional profiling analysis in different laying stages in chicken

During follicular development, a series of key events such as follicular recruitment and selection are crucially governed by strict complex regulation. However, its molecular mechanisms remain obscure. To identify the dominant genes controlling chicken follicular development, the small white follicle (SWF), the small yellow follicle (SYF), and the large yellow follicle (LYF) in different laying stages (W22, W31, W51) were collected for RNA sequencing and bioinformatics analysis. There were 1866, 1211, and 1515 differentially expressed genes (DEGs) between SWF and SYF in W22, W31, and W51, respectively. 4021, 2295, and 2902 DEGs were respectively identified between SYF and LYF in W22, W31, and W51. 5618, 4016, and 4809 DEGs were respectively identified between SWF and LYF in W22, W31, and W51. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that extracellular matrix, extracellular region, extracellular region part, ECM-receptor interaction, collagen extracellular matrix, and collagen trimer were significantly enriched (P < 0.05). Protein–protein interaction analysis revealed that COL4A2, COL1A2, COL4A1, COL5A2, COL12A1, ELN, ALB, and MMP10 might be key candidate genes for follicular development in chicken. The current study identified dominant genes and pathways contributing to our understanding of chicken follicular development. 

Gene network analysis reveals candidate genes related with the hair follicle development in sheep

Background

Merino sheep are the most famous fine wool sheep in the world. They have high wool production and excellent wool quality and have attracted worldwide attention. The fleece of the Merino sheep is composed predominantly of wool fibers grown from secondary wool follicles. Therefore, it is necessary to study the development of hair follicles to understand the mechanism of wool production. The hair follicle is a complex biological system involved in a dynamic process governed by gene regulation. The hair follicle development process is very complex and poorly understood. The purpose of our research is to identify candidate genes related to hair follicle development, provide a theoretical molecular breeding basis for the cultivation of fine wool sheep, and provide a reference for the problems of hair loss and alopecia areata that affect human beings.

Results

We analyzed mRNAs data in skin tissues of 18 Merino sheep at four embryonic days (E65, E85, E105 and E135) and two postnatal days (P7 and P30). G1 to G6 represent hair follicles developmental at six stages (i.e. E65 to P30). We identified 7879 differentially expressed genes (DEGs) and 12623 novel DEGs, revealed different expression patterns of these DEGs at six stages of hair follicle development, and demonstrated their complex interactions. DEGs with stage-specific expression were significantly enriched in epidermal differentiation and development, hair follicle development and hair follicle morphogenesis and were enriched in many pathways related to hair follicle development. The key genes (LAMA5, WNT10A, KRT25, SOSTDC1, ZDHHC21, FZD1, BMP7, LRP4, TGFβ2, TMEM79, SOX10, ITGB4, KRT14, ITGA6, and GLI2) affecting hair follicle morphogenesis were identified by network analysis.

Conclusion

This study provides a new reference for the molecular basis of hair follicle development and lays a foundation for further improving sheep hair follicle breeding. Candidate genes related to hair follicular development were found, which provided a theoretical basis for molecular breeding for the culture of fine wool sheep. These results are a valuable resource for biological investigations of fleece evolution in animals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08552-2.

Effects of ischemic postconditioning and long non-coding RNAs in ischemic stroke

Stroke is a main cause of disability and death among adults in China, and acute ischemic stroke accounts for 80% of cases. The key to ischemic stroke treatment is to recanalize the blocked blood vessels. However, more than 90% of patients cannot receive effective treatment within an appropriate time, and delayed recanalization of blood vessels causes reperfusion injury. Recent research has revealed that ischemic postconditioning has a neuroprotective effect on the brain, but the mechanism has not been fully clarified. Long non-coding RNAs (lncRNAs) have previously been associated with ischemic reperfusion injury in ischemic stroke. LncRNAs regulate important cellular and molecular events through a variety of mechanisms, but a comprehensive analysis of potential lncRNAs involved in the brain protection produced by ischemic postconditioning has not been conducted. In this review, we summarize the common mechanisms of cerebral injury in ischemic stroke and the effect of ischemic postconditioning, and we describe the potential mechanisms of some lncRNAs associated with ischemic stroke.

Integrated analysis of lncRNAs and mRNAs by RNA-Seq in secondary hair follicle development and cycling (anagen, catagen and telogen) of Jiangnan cashmere goat (Capra hircus)

BACKGROUND

Among the world's finest natural fiber composites is derived from the secondary hair follicles (SHFs) of cashmere goats yield one of the world's best natural fibres. Their development and cycling are characterized by photoperiodism with diverse, well-orchestrated stimulatory and inhibitory signals. Long non-coding RNA (lncRNAs) and mRNAs play important roles in hair follicle (HF) development. However, not many studies have explored their specific functions in cashmere development and cycling. This study detected mRNAs and lncRNAs with their candidate genes and related pathways in SHF development and cycling of cashmere goat. We utilized RNA sequencing (RNA-Seq) and bioinformatics analysis on lncRNA and mRNA expressions in goat hair follicles to discover candidate genes and metabolic pathways that could affect development and cycling (anagen, catagen, and telogen).

RESULTS

We identified 228 differentially expressed (DE) mRNAs and 256 DE lncRNA. For mRNAs, catagen and anagen had 16 upregulated and 35 downregulated DEGs, catagen and telogen had 18 upregulated and 9 downregulated DEGs and telogen and anagen had 52 upregulated and 98 downregulated DEGs. LncRNA witnessed 22 upregulated and 39 downregulated DEGs for catagen and anagen, 36 upregulated and 29 downregulated DEGs for catagen and telogen as well as 66 upregulated and 97 downregulated DEGs for telogen and anagen. Several key genes, including MSTRG.5451.2, MSTRG.45465.3, MSTRG.11609.2, CHST1, SH3BP4, CDKN1A, GAREM1, GSK-3β, DEFB103A KRTAP9-2, YAP1, S100A7A, FA2H, LOC102190037, LOC102179090, LOC102173866, KRT2, KRT39, FAM167A, FAT4 and EGFL6 were shown to be potentially important in hair follicle development and cycling. They were related to, WNT/β-catenin, mTORC1, ERK/MAPK, Hedgehog, TGFβ, NFkB/p38MAPK, caspase-1, and interleukin (IL)-1a signaling pathways.

CONCLUSION

This work adds to existing understanding of the regulation of HF development and cycling in cashmere goats via lncRNAs and mRNAs. It also serves as theoretical foundation for future SHF research in cashmere goats.

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.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

Differential expression of long non-coding RNAs under Peste des petits ruminants virus (PPRV) infection in goats

Peste des petits ruminants (PPR) characterized by fever, sore mouth, conjunctivitis, gastroenteritis, and pneumonia, is an acute, highly contagious viral disease of sheep and goats. The role of long non-coding RNAs (lncRNAs) in PPRV infection has not been explored to date. In this study, the transcriptome profiles of virulent Peste des petits ruminants virus (PPRV) infected goat tissues – lung and spleen were analyzed to identify the role of lncRNAs in PPRV infection. A total of 13,928 lncRNA transcripts were identified, out of which 170 were known lncRNAs. Intergenic lncRNAs (7625) formed the major chunk of the novel lncRNA transcripts. Differential expression analysis revealed that 15 lncRNAs (11 downregulated and 4 upregulated) in the PPRV infected spleen samples and 16 lncRNAs (13 downregulated and 3 upregulated) in PPRV infected lung samples were differentially expressed as compared to control. The differentially expressed lncRNAs (DElncRNAs) possibly regulate various immunological processes related to natural killer cell activation, antigen processing and presentation, and B cell activity, by regulating the expression of mRNAs through the cis- or trans-regulatory mechanism. Functional enrichment analysis of differentially expressed mRNAs (DEmRNAs) revealed enrichment of immune pathways and biological processes in concordance with the pathways in which correlated lncRNA-neighboring genes were enriched. The results suggest that a coordinated immune response is raised in both lung and spleen tissues of the goat through mRNA-lncRNA crosstalk.

Identification of Potential Blind-Side Hypermelanosis-Related lncRNA–miRNA–mRNA Regulatory Network in a Flatfish Species, Chinese Tongue Sole (Cynoglossus semilaevis)

Blind-side hypermelanosis has emerged as a major concern in commercial rearing environments of the flatfish aquaculture industry. To date, the underlying molecular mechanisms are not well understood. To fill this gap, in this study, whole transcriptomic sequencing and analyses were performed using normal skins and hypermelanic skins of the blind side of Chinese tongue sole (Cynoglossus semilaevis). Differentially expressed long non-coding RNAs (DElncRNAs), miRNAs (DEmiRNAs), and differentially expressed genes as well as their competing endogenous RNA (ceRNA) networks were identified. A total of 34 DElncRNAs, 226 DEmiRNAs, and 610 DEGs were identified. Finally, lncRNA–miRNA–mRNA regulatory networks (involving 29 DElncRNAs, 106 DEmiRNAs, and 162 DEGs) associated with blind-side hypermelanosis were constructed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of 162 DEGs in ceRNA networks identified DEGs (e.g., oca2, mc1r, and ihhb) in pigmentation-related biological processes and DEGs (e.g., ca4, glul, and fut9) in nitrogen metabolism, glycosphingolipid biosynthesis, and folate biosynthesis pathways, as well as their corresponding DElncRNAs and DEmiRNAs to potentially play key regulatory roles in blind-side hypermelanosis. In conclusion, this is the first study on the ceRNA regulatory network associated with blind-side hypermelanosis in flatfish. These new findings expand the spectrum of non-coding regulatory mechanisms underpinning blind-side hypermelanosis, which facilitates the further exploration of molecular regulatory mechanisms of malpigmentation in flatfish.

Characterization of XR_311113.2 as a MicroRNA Sponge for Pre-ovulatory Ovarian Follicles of Goats via Long Noncoding RNA Profile and Bioinformatics Analysis

Long noncoding RNAs (lncRNAs) were identified recently as a large class of noncoding RNAs (ncRNAs) with a length ≥200 base pairs (bp). The function and mechanism of lncRNAs have been reported in a growing number of species and tissues. In contrast, the regulatory mechanism of lncRNAs in the goat reproductive system has rarely been reported. In the present study, we sequenced and analyzed the lncRNAs using bioinformatics to identify their expression profiles. As a result, 895 lncRNAs were predicted in the pre-ovulatory ovarian follicles of goats. Eighty-eight lncRNAs were differentially expressed in the Macheng black goat when compared with Boer goat. In addition, the lncRNA XR_311113.2 acted as a sponge of chi-miR-424-5p, as assessed via a luciferase activity assay. Taken together, our findings demonstrate that lncRNAs have potential effects in the ovarian follicles of goats and may represent a promising new research field to understand follicular development.

Neuroprotective effects of long noncoding RNAs involved in ischemic postconditioning after ischemic stroke

During acute reperfusion, the expression profiles of long noncoding RNAs in adult rats with focal cerebral ischemia undergo broad changes. However, whether long noncoding RNAs are involved in neuroprotective effects following focal ischemic stroke in rats remains unclear. In this study, RNA isolation and library preparation was performed for long noncoding RNA sequencing, followed by determining the coding potential of identified long noncoding RNAs and target gene prediction. Differential expression analysis, long noncoding RNA functional enrichment analysis, and co-expression network analysis were performed comparing ischemic rats with and without ischemic postconditioning rats. Rats were subjected to ischemic postconditioning via the brief and repeated occlusion of the middle cerebral artery or femoral artery. Quantitative real-time reverse transcription-polymerase chain reaction was used to detect the expression levels of differentially expressed long noncoding RNAs after ischemic postconditioning in a rat model of ischemic stroke. The results showed that ischemic postconditioning greatly affected the expression profile of long noncoding RNAs and mRNAs in the brains of rats that underwent ischemic stroke. The predicted target genes of some of the identified long noncoding RNAs (cis targets) were related to the cellular response to ischemia and stress, cytokine signal transduction, inflammation, and apoptosis signal transduction pathways. In addition, 15 significantly differentially expressed long noncoding RNAs were identified in the brains of rats subjected to ischemic postconditioning. Nine candidate long noncoding RNAs that may be related to ischemic postconditioning were identified by a long noncoding RNA expression profile and long noncoding RNA-mRNA co-expression network analysis. Expression levels were verified by quantitative real-time reverse transcription-polymerase chain reaction. These results suggested that the identified long noncoding RNAs may be involved in the neuroprotective effects associated with ischemic postconditioning following ischemic stroke. The experimental animal procedures were approved by the Animal Experiment Ethics Committee of Kunming Medical University (approval No. KMMU2018018) in January 2018.

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.

A Single-cell Transcriptome Atlas of Cashmere Goat Hair Follicle Morphogenesis

Cashmere, also known as soft gold, is produced from the secondary hair follicles (SHFs) of cashmere goats. The number of SHFs determines the yield and quality of cashmere; therefore, it is of interest to investigate the transcriptional profiles present during cashmere goat hair follicle development. However, mechanisms underlying this development process remain largely unexplored, and studies regarding hair follicle development mostly use a murine research model. In this study, to provide a comprehensive understanding of cellular heterogeneity and cell fate decisions, single-cell RNA sequencing was performed on 19,705 single cells of the dorsal skin from cashmere goat fetuses at induction (embryonic day 60; E60), organogenesis (E90), and cytodifferentiation (E120) stages. For the first time, unsupervised clustering analysis identified 16 cell clusters, and their corresponding cell types were also characterized. Based on lineage inference, a detailed molecular landscape was revealed along the dermal and epidermal cell lineage developmental pathways. Notably, our current data also confirmed the heterogeneity of dermal papillae from different hair follicle types, which was further validated by immunofluorescence analysis. The current study identifies different biomarkers during cashmere goat hair follicle development and has implications for cashmere goat breeding in the future.

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.

Exosomes in Immune Regulation

Exosomes, small extracellular vesicles mediate intercellular communication by transferring their cargo including DNA, RNA, proteins and lipids from cell to cell. Notably, in the immune system, they have protective functions. However in cancer, exosomes acquire new, immunosuppressive properties that cause the dysregulation of immune cells and immune escape of tumor cells supporting cancer progression and metastasis. Therefore, current investigations focus on the regulation of exosome levels for immunotherapeutic interventions. In this review, we discuss the role of exosomes in immunomodulation of lymphoid and myeloid cells, and their use as immune stimulatory agents to elicit specific cytotoxic responses against the tumor.

Transcriptome analysis of long noncoding RNAs ribonucleic acids from the livers of Hu sheep with different residual feed intake

Long noncoding RNAs (LncRNAs), as key regulators, have vital functions in various biological activities. However, in sheep, little has been reported concerning the genetic mechanism of LncRNA regulation of feed efficiency. In the present study, we explored the genome-wide expression of LncRNAs and transcripts of uncertain coding potential (TUCPs) in the livers of sheep with extreme residual feed intake (RFI) using RNA sequencing. We identified 1 523 TUCPs and 1 996 LncRNAs, among which 10 LncRNAs and 16 TUCPs were identified as being differentially expressed between the High-RFI and Low-RFI groups. Co-expression and co-localization methods were used to search for LncRNA and TUCP target genes, which identified 970/1 538 and 23/27 genes, respectively. Ontology and pathways analysis revealed that the LncRNAs/TUCPs that were highly expressed in the Low-RFI group are mostly concentrated in energy metabolism pathways. For example, LNC_000890 and TUCP_000582 might regulate liver tissue metabolic efficiency. The LncRNAs/TUCPs that were highly expressed in the High-RFI group are mostly enriched in immune function pathways. For example, TUCP_000832 might regulate animal health, thereby affecting feed efficiency. Subsequently, a co-expression network was established by applying the expression information of both the differentially expressed LncRNAs and TUCPs and their target mRNAs. The network indicated that differentially expressed genes targeted by the upregulated LncRNAs and TUCPs were mainly related to energy metabolism, while those genes targeted by the downregulated LncRNAs and TUCPs were mainly related to immune response. These results provide the basis for further study of LncRNA/TUCP-mediated regulation of feed efficiency.

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.

Hepatic transcriptomic adaptation from prepartum to postpartum in dairy cows

The transition from pregnancy to lactation is the most challenging period for high-producing dairy cows. The liver plays a key role in biological adaptation during the peripartum. Prior works have demonstrated that hepatic glucose synthesis, cholesterol metabolism, lipogenesis, and inflammatory response are increased or activated during the peripartum in dairy cows; however, those works were limited by a low number of animals used or by the use of microarray technology, or both. To overcome such limitations, an RNA sequencing analysis was performed on liver biopsies from 20 Holstein cows at 7 ± 5d before (Pre-P) and 16 ± 2d after calving (Post-P). We found 1,475 upregulated and 1,199 downregulated differently expressed genes (DEG) with a false discovery rate adjusted P-value < 0.01 between Pre-P and Post-P. Bioinformatic analysis revealed an activation of the metabolism, especially lipid, glucose, and amino acid metabolism, with increased importance of the mitochondria and a key role of several signaling pathways, chiefly peroxisome proliferators-activated receptor (PPAR) and adipocytokines signaling. Fatty acid oxidation and gluconeogenesis, with a likely increase in amino acid utilization to produce glucose, were among the most important functions revealed by the transcriptomic adaptation to lactation in the liver. Although gluconeogenesis was induced, data indicated decrease in expression of glucose transporters. The analysis also revealed high activation of cell proliferation but inhibition of xenobiotic metabolism, likely due to the liver response to inflammatory-like conditions. Co-expression network analysis disclosed a tight connection and coordination among genes driving biological processes associated with protein synthesis, energy and lipid metabolism, and cell proliferation. Our data confirmed the importance of metabolic adaptation to lipid and glucose metabolism in the liver of early Post-P cows, with a pivotal role of PPAR and adipocytokines.

Exploration of Long Non-coding RNAs and Circular RNAs in Porcine Milk Exosomes

RNA in milk exosomes can be absorbed in the mammalian intestinal tract and function in gene expression regulations. Our previous work demonstrated that porcine milk exosomes (PME) contain large amounts of miRNAs and mRNAs. Increasing evidence suggests that long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are of particular interest, given their key role in diverse biological processes of animals. However, the expression profiles and the potential functions of lncRNAs and circRNAs in PME are still unknown. In the present study, we isolated PME by ultracentrifugation and performed a comprehensive analysis of lncRNA and circRNA in PME by using RNA sequencing. As a result, 2,466 novel lncRNAs, 809 annotated lncRNAs, and 61 circRNAs were identified in PME. The lncRNAs shared similar characteristics with other mammals in terms of length, exon number, and open reading frames. However, lncRNAs showed a higher level compared with mRNAs. Eight lncRNAs and five circRNAs in PME were selected for PCR identification. A functional enrichment analysis on the target genes of lncRNAs indicated that these genes were involved in cellular macromolecule metabolic, RNA metabolic, and immune processes. The circRNAs host genes were enriched in nucleic acid binding and adherence junction. We also evaluated the potential interaction targets between miRNAs and PME lncRNAs or circRNAs, and the results showed that the PME lncRNAs and the circRNAs have a high density of miRNA target sites. The top 20 highly expressed lncRNAs were found to interact with the proliferation-related miRNAs, and the circRNAs potentially targeted many miRNAs that are associated with the intestinal barrier. This study is the first to provide a resource for lncRNA and circRNA research of porcine milk. Moreover, the potential interaction between lncRNA/circRNA and miRNA is revealed. The present study expands our knowledge of non-coding RNAs in milk, and additional research is necessary to confirm their exactly physiological functions.

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.

Clinical Ketosis-Associated Alteration of Gene Expression in Holstein Cows

Ketosis is one of the most prevalent transition metabolic disorders in dairy cows, and has been intrinsically influenced by both genetic and nutritional factors. However, altered gene expression with respective to dairy cow ketosis has not been addressed yet, especially at the genome-wide level. In this study, we recruited nine Holsteins diagnosed with clinical ketosis and ten healthy controls, for which whole blood samples were collected at both prepartum and postpartum. Four groups of blood samples were defined: from cows with ketosis at prepartum (PCK, N = 9) and postpartum (CK, N = 9), respectively, and controls at prepartum (PHC, N = 10) and postpartum (HC, N = 10). RNA-Seq approach was used for investigating gene expression, by which a total of 27,233 genes were quantified with four billion high-quality reads. Subsequently, we revealed 75 and four differentially expressed genes (DEGs) between sick and control cows at postpartum and prepartum, respectively, which indicated that sick and control cows had similar gene expression patterns at prepartum. Meanwhile, there were 95 DEGs between postpartum and prepartum for sick cows, which showed depressed changes of gene expression during this transition period in comparison with healthy cows (428 DEGs). Functional analyses revealed the associated DEGs with ketosis were mainly involved in biological stress response, ion homeostasis, AA metabolism, energy signaling, and disease related pathways. Finally, we proposed that the expression level of STX1A would be potentially used as a new biomarker because it was the only gene that was highly expressed in sick cows at both prepartum and postpartum. These results could significantly help us to understand the underlying molecular mechanisms for incidence and progression of ketosis in dairy cows.

Genome-wide identification and characterization of novel long non-coding RNA in Ruminal tissue affected with sub-acute Ruminal acidosis from Holstein cattle

Sub-acute ruminal acidosis is a type of metabolic disorder in which affected cattle show a considerable depression of rumen pH. This leads to a dramatic decline in productivity and consequent loss of income for many dairy farms. The objective of the present study is to identify and characterize novel long non-coding RNAs (lncRNAs) in Holstein cattle affected by sub-acute ruminal acidosis. Two replicates from six animals were sequenced that bioinformatically analyzed. Results showed 6679 novel lncRNAs among which 12 intergenic lncRNAs showed differential expression (p value ≤0.05). GO and KEGG analysis revealed that calcium signaling and G protein couple-receptor pathways may be involved in regulating metabolic processes during sub-acute ruminal acidosis. Furthermore, other biological processes including transmembrane transport, adult behavior, neuroactive ligand-receptor interaction, GABAergic synapse, cholinergic synapse were significantly enriched. The present data suggest that these differentially expressed lncRNAs may play regulatory roles in modulating biological processes associated with sub-acute ruminal acidosis in cattle rumen.

The comprehensive detection of miRNA, lncRNA, and circRNA in regulation of mouse melanocyte and skin development

BACKGROUND

Pigmentation development, is a complex process regulated by many transcription factors during development. With the development of the RNA sequencing (RNA-seq), non-coding RNAs, such as miRNAs, lncRNAs, and circRNAs, are found to play an important role in the function detection of related regulation factors. In this study, we provided the expression profiles and development of ncRNAs related to melanocyte and skin development in mice with black coat color skin and mice with white coat color skin during embryonic day 15 (E15) and postnatal day 7 (P7). The expression profiles of different ncRNAs were detected via RNA-seq and also confirmed by the quantitative real-time PCR (qRT-PCR) method. GO and KEGG used to analyze the function the related target genes.

RESULTS

We identified an extensive catalogue of 206 and 183 differently expressed miRNAs, 600 and 800 differently expressed lncRNAs, and 50 and 54 differently expressed circRNAs, respectively. GO terms and pathway analysis showed the target genes of differentially expressed miRNA and lncRNA. The host genes of circRNA were mainly enriched in cellular process, single organism process. The target genes of miRNAs were mainly enriched in chromatin binding and calcium ion binding in the nucleus. The function of genes related to lncRNAs are post translation modification. The competing endogenous RNA (ceRNA) network of lncRNAs and circRNAs displays a complex interaction between ncRNA and mRNA related to skin development, such as Tcf4, Gnas, and Gpnms related to melanocyte development.

CONCLUSIONS

The ceRNA network of lncRNA and circRNA displays a complex interaction between ncRNA and mRNA related to skin development and melanocyte development. The embryonic and postnatal development of skin provide a reference for further studies on the development mechanisms of ncRNA during pigmentation.

Genome-wide analysis of mRNAs and lncRNAs in Mycoplasma bovis infected and non-infected bovine mammary gland tissues

Mycoplasma bovis (M. bovis) causes diseases such as arthritis, pneumonia, abortion, and mastitis, leading to great losses in the bovine dairy industries. RNA types such as messenger RNAs (mRNAs) and long non-coding RNAs (lncRNAs) play significant roles in regulating the immune responses triggered by bacteria. The expression profiles of mRNA and lncRNA as they occur in bovine mammary gland tissues infected with M. bovis are still not well understood. To illuminate this issue, transcription analysis of mRNA and LncRNAs were conducted on the mammary gland tissues belonging to Holstein cattle infected and not infected with M. bovis. The analysis revealed 1310 differentially expressed mRNAs and 57 differentially expressed lncRNAs in the bovine mammary gland tissues infected and not infected with M. bovis. In addition, 392 novel lncRNAs were detected, 19 of which were differentially expressed. Gene ontology analysis reveals that differentially expressed mRNAs and lncRNAs play significant roles in such vital biological pathways as metabolic pathways, T-cell receptor signaling, TGF-beta signaling, pathways in cancer, PI3K-Akt signaling, NF-kappa B signaling, mTOR signaling, and apoptosis, including in the immune response to cancer. Based on our literature review, this study is the first genome-wide lncRNA research conducted on bovine mammary gland tissues infected with M. bovis. Our results provide bovine mammary gland lncRNA and mRNA resources to understand their roles in the regulation of the immune response against the agent M. bovis in bovine mammary gland tissues.

Transcriptional Differences of Coding and Non-Coding Genes Related to the Absence of Melanocyte in Skins of Bama Pig

Skin is the body’s largest organ, and the main function of skin is to protect underlying organs from possible external damage. Melanocytes play an important role in skin pigmentation. The Bama pig has a “two-end-black” phenotype with different coat colors across skin regions, e.g., white skin (without melanocytes) and black skin (with melanocytes), which could be a model to investigate skin-related disorders, specifically loss of melanocytes. Here, we generated expression profiles of mRNAs and long noncoding RNAs in Bama pig skins with different coat colors. In total, 14,900 mRNAs and 7549 lncRNAs were expressed. Overall, 2338 mRNAs/113 lncRNAs with FDR-adjusted p-value ≤ 0.05 were considered to be differentially expressed (DE) mRNAs/lncRNAs, with 1305 down-regulated mRNAs and 1033 up-regulated mRNAs in white skin with|log₂(fold change)| > 1. The genes down-regulated in white skin were associated with pigmentation, melanocyte–keratinocyte interaction, and keratin, while up-regulated ones were mainly associated with cellular energy metabolisms. Furthermore, those DE lncRNAs were predicted to be implicated in pigmentation, keratin synthesis and cellular energy metabolism. In general, this study provides insight into the transcriptional difference involved in melanocyte-loss-induced keratinocyte changes and promotes the Bama pig as a biomedical model in skin research.

Transcriptional profiling of long noncoding RNAs associated with leaf-color mutation in Ginkgo biloba L

BACKGROUND

Long noncoding RNAs (lncRNAs) play an important role in diverse biological processes and have been widely studied in recent years. However, the roles of lncRNAs in leaf pigment formation in ginkgo (Ginkgo biloba L.) remain poorly understood.

RESULTS

In this study, lncRNA libraries for mutant yellow-leaf and normal green-leaf ginkgo trees were constructed via high-throughput sequencing. A total of 2044 lncRNAs were obtained with an average length of 702 nt and typically harbored 2 exons. We identified 238 differentially expressed lncRNAs (DELs), 32 DELs and 49 differentially expressed mRNAs (DEGs) that constituted coexpression networks. We also found that 48 cis-acting DELs regulated 72 target genes, and 31 trans-acting DELs regulated 31 different target genes, which provides a new perspective for the regulation of the leaf-color mutation. Due to the crucial regulatory roles of lncRNAs in a wide range of biological processes, we conducted in-depth studies on the DELs and their targets and found that the chloroplast thylakoid membrane subcategory and the photosynthesis pathways (ko00195) were most enriched, suggesting their potential roles in leaf coloration mechanisms. In addition, our correlation analysis indicates that eight DELs and 68 transcription factors (TFs) might be involved in interaction networks.

CONCLUSIONS

This study has enriched the knowledge concerning lncRNAs and provides new insights into the function of lncRNAs in leaf-color mutations, which will benefit future selective breeding of ginkgo.

Genome-wide identification and characterization of long non-coding RNAs during differentiation of visceral preadipocytes in rabbit

Emerging evidence suggests that long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes, including adipogenesis. Despite being considered an ideal animal model for studying adipogenesis, little is known about the roles of lncRNAs in the regulation of rabbit preadipocyte differentiation. In the present study, visceral preadipocytes isolated from newborn rabbits were cultured in vitro and induced for differentiation, and global lncRNA expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. A total of 2066 lncRNAs were identified from nine RNA-seq libraries. Compared to protein-coding transcripts, lncRNA transcripts exhibited characteristics of a longer length and lower expression level. Furthermore, 486 and 357 differentially expressed (DE) lncRNAs were identified when comparing day 3 vs. day 0 and day 9 vs. day 3, respectively. Target genes of DE lncRNAs were predicted by the cis-regulating approach. Prediction of functions revealed that DE lncRNAs when comparing day 3 vs. day 0 were involved in gene ontology (GO) terms of developmental growth, growth, developmental cell growth, and stem cell proliferation, and involved in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of PI3K-Akt signaling pathway, fatty acid biosynthesis, and the insulin signaling pathway. The DE lncRNAs when comparing day 9 vs. day 3 were involved in GO terms that associated with epigenetic modification and were involved in the KEGG pathway of cAMP signaling pathway. This study provides further insight into the regulatory function of lncRNAs in rabbit visceral adipose and facilitates a better understanding of different stages of preadipocyte differentiation.

Circulating non-coding RNAs in recurrent and metastatic ovarian cancer

Ovarian cancer has a poor outcome because it is usually detected at advanced tumor stages, and the majority of the patients develop disease relapse as a result of chemotherapy resistance. This most lethal gynecological malignancy metastasizes within the peritoneal fluid or ascites to pelvic and distal organs. In ovarian cancer progression and metastasis, small non-coding RNAs (ncRNAs), including long noncoding RNAs and microRNAs have been recognized as important regulators. Their dysregulation modulates gene expression and cellular signal pathways and can be detected in liquid biopsies. In this review, we provide an overview on circulating plasma and serum ncRNAs participating in tumor cell migration and invasion, and contributing to recurrence and metastasis of ovarian cancer. We will also discuss the development of potential, novel therapies using ncRNAs as target molecules or tumor markers for ovarian cancer.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

Heat stress negatively affects the transcriptome related to overall metabolism and milk protein synthesis in mammary tissue of midlactating dairy cows

Inadequate dry matter intake only partially accounts for the decrease in milk protein synthesis during heat stress (HS) in dairy cows. Our hypothesis is that reduced milk protein synthesis during HS in dairy cows is also caused by biological changes within the mammary gland. The objective of this study was to assess the hypothesis via RNA-Seq analysis of mammary tissue. Herein, four dairy cows were used in a crossover design where HS was induced for 9 days in environmental chambers. There was a 30-day washout between periods. Mammary tissue was collected via biopsy at the end of each environmental period (HS or pair-fed and thermal neutral) for transcriptomic analysis. RNA-Seq analysis revealed HS affected >2,777 genes (false discovery rate-adjusted P value < 0.05) in mammary tissue. Expression of main milk protein-encoding genes and several key genes related to regulation of protein synthesis and amino acid and glucose transport were downregulated by HS. Bioinformatics analysis revealed an overall decrease of mammary tissue metabolic activity by HS (especially carbohydrate and lipid metabolism) and an increase in immune activation and inflammation. Network analysis revealed a major role of TNF, IFNG, S100A8, S100A9, and IGF-1 in inducing/controlling the inflammatory response, with a central role of NF-κB in the process of immunoactivation. The same analysis indicated an overall inhibition of PPARγ. Collectively, these data suggest HS directly controls milk protein synthesis via reducing the transcription of metabolic-related genes and increasing inflammation-related genes.

Systematic identification of intergenic long-noncoding RNAs in mouse retinas using full-length isoform sequencing

Background

A great mass of long noncoding RNAs (lncRNAs) have been identified in mouse genome and increasing evidences in the last decades have revealed their crucial roles in diverse biological processes. Nevertheless, the biological roles of lncRNAs in the mouse retina remains largely unknown due to the lack of a comprehensive annotation of lncRNAs expressed in the retina.

Results

In this study, we applied the long-reads sequencing strategy to unravel the transcriptomes of developing mouse retinas and identified a total of 940 intergenic lncRNAs (lincRNAs) in embryonic and neonatal retinas, including about 13% of them were transcribed from unannotated gene loci. Subsequent analysis revealed that function of lincRNAs expressed in mouse retinas were closely related to the physiological roles of this tissue, including 90 lincRNAs that were differentially expressed after the functional loss of key regulators of retinal ganglion cell (RGC) differentiation. In situ hybridization results demonstrated the enrichment of three class IV POU-homeobox genes adjacent lincRNAs (linc-3a, linc-3b and linc-3c) in ganglion cell layer and indicated they were potentially RGC-specific.

Conclusions

In summary, this study systematically annotated the lincRNAs expressed in embryonic and neonatal mouse retinas and implied their crucial regulatory roles in retinal development such as RGC differentiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5903-y) contains supplementary material, which is available to authorized users.

Integrated analysis of long non-coding RNA and mRNA expression in different colored skin of koi carp

BACKGROUND

Long non-coding RNAs (lncRNAs) perform crucial roles in biological process involving complex mechanisms. However, information regarding their abundance, characteristics and potential functions linked to fish skin color is limited. Herein, Illumina sequencing and bioinformatics were conducted on black, white, and red skin of Koi carp (Cyprinus carpio L.).

RESULTS

A total of 590,415,050 clean reads, 446,614 putative transcripts, 4252 known and 72,907 novel lncRNAs were simultaneously obtained, including 92 significant differentially expressed lncRNAs and 722 mRNAs. Ccr_lnc5622441 and Ccr_lnc765201 were up-regulated in black and red skin, Ccr_lnc14074601 and Ccr_lnc2382951 were up-regulated in white skin, and premelanosome protein a (Pmela), Pmelb and tyrosinase (Tyr) were up-regulated in black skin. The expression patterns of 18 randomly selected differentially expressed genes were validated using the quantitative real-time PCR method. Moreover, 70 lncRNAs acting on 107 target mRNAs in cis and 79 lncRNAs acting on 41,625 target mRNAs in trans were investigated. The resulting co-expression networks revealed that a single lncRNA can connect with numerous mRNAs, and vice versa. To further reveal their potential functions, Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed, and membrane, pigment cell development, cAMP signaling, melanogenesis and tyrosine metabolism appear to affect skin pigmentation. Additionally, three lncRNAs (Ccr_lnc142711, Ccr_lnc17214525 and Ccr_lnc14830101) and three mRNAs (Asip, Mitf and Tyr) involved in the melanogenesis pathway were investigated in terms of potential functions in embryogenesis and different tissues.

CONCLUSIONS

The findings broaden our understanding of lncRNAs and skin color genetics, and provide new insight into the mechanisms underlying lncRNA-mediated pigmentation and differentiation in Koi carp.

Synchronous profiling and analysis of mRNAs and ncRNAs in the dermal papilla cells from cashmere goats

Background

Dermal papilla cells (DPCs), the “signaling center” of hair follicle (HF), delicately master continual growth of hair in mammals including cashmere, the fine fiber annually produced by secondary HF embedded in cashmere goat skins. Such unparalleled capacity bases on their exquisite character in instructing the cellular activity of hair-forming keratinocytes via secreting numerous molecular signals. Past studies suggested microRNA (miRNAs) and long non-coding RNAs (lncRNAs) play essential roles in a wide variety of biological process, including HF cycling. However, their roles and related molecular mechanisms in modulating DPCs secretory activities are still poorly understood.

Results

Here, we separately cultivated DPCs and their functionally and morphologically distinct dermal fibroblasts (DFs) from cashmere goat skins at anagen. With the advantage of high throughput RNA-seq, we synchronously identified 2540 lncRNAs and 536 miRNAs from two types of cellular samples at 4th passages. Compared with DFs, 1286 mRNAs, 18 lncRNAs, and 42 miRNAs were upregulated, while 1254 mRNAs, 53 lncRNAs and 44 miRNAs were downregulated in DPCs. Through overlapping with mice data, we ultimately defined 25 core signatures of DPCs, including HOXC8 and RSPO1, two crucial activators for hair follicle stem cells (HFSCs). Subsequently, we emphatically investigated the impacts of miRNAs and lncRNAs (cis- and trans- acting) on the genes, indicating that ncRNAs extensively exert negative and positive effects on their expressions. Furthermore, we screened lncRNAs acting as competing endogenous RNAs (ceRNAs) to sponge miRNAs and relief their repressive effects on targeted genes, and constructed related lncRNAs-miRNAs-HOXC8/RSPO1 interactive lines using bioinformatic tools. As a result, XR_310320.3-chi-miR-144-5p-HOXC8, XR_311077.2-novel_624-RSPO1 and others lines appeared, displaying that lncRNAs might serve as ceRNAs to indirectly adjust HFSCs status in hair growth.

Conclusion

The present study provides an unprecedented inventory of lncRNAs, miRNAs and mRNAs in goat DPCs and DFs. We also exhibit some miRNAs and lncRNAs potentially participate in the modulation of HFSCs activation via delicately adjusting core signatures of DPCs. Our report shines new light on the latent roles and underlying molecular mechanisms of ncRNAs on hair growth.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5861-4) contains supplementary material, which is available to authorized users.

Screening and evaluating of long non-coding RNAs in prenatal and postnatal pituitary gland of sheep

Long non-coding RNAs are transcribed into RNA molecules that are >200 nucleotides in length. However, the expression and function analysis of lncRNAs in the sheep pituitary gland are still lacking. In this study, we identified 1755 lncRNAs (545 annotated lncRNAs and 1210 novel lncRNAs) from RNA-seq data in the pituitary gland of embryonic and adult sheep. A total of 235 lncRNAs were differentially expressed between embryonic and adult group. We verified the presence of some lncRNAs using RT-PCR and DNA sequencing, and identified some differentially expressed lncRNAs using qPCR. We also investigated the role of cis-acting lncRNAs on target genes. GO and KEGG enrichment analysis revealed that the target genes of lncRNAs were involved in the regulation of hormones secretion and some signaling pathways in the sheep pituitary gland. Our study provides comprehensive expression profiles of lncRNAs and valuable resource for understanding their function in the pituitary gland.

Comprehensive Analysis of lncRNAs and circRNAs Reveals the Metabolic Specialization in Oxidative and Glycolytic Skeletal Muscles

The biochemical and functional differences between oxidative and glycolytic muscles could affect human muscle health and animal meat quality. However, present understanding of the epigenetic regulation with respect to lncRNAs and circRNAs is rudimentary. Here, porcine oxidative and glycolytic skeletal muscles, which were at the growth curve inflection point, were sampled to survey variant global expression of lncRNAs and circRNAs using RNA-seq. A total of 4046 lncRNAs were identified, including 911 differentially expressed lncRNAs (p < 0.05). The cis-regulatory analysis identified target genes that were enriched for specific GO terms and pathways (p < 0.05), including the oxidation-reduction process, glycolytic process, and fatty acid metabolic. All these were closely related to different phenotypes between oxidative and glycolytic muscles. Additionally, 810 circRNAs were identified, of which 137 were differentially expressed (p < 0.05). Interestingly, some circRNA-miRNA-mRNA networks were found, which were closely linked to muscle fiber-type switching and mitochondria biogenesis in muscles. Furthermore, 44.69%, 39.19%, and 54.01% of differentially expressed mRNAs, lncRNAs, and circRNAs respectively were significantly enriched in pig quantitative trait loci (QTL) regions for growth and meat quality traits. This study reveals a mass of candidate lncRNAs and circRNAs involved in muscle physiological functions, which may improve understanding of muscle metabolism and development from an epigenetic perspective.

Genome-wide identification and characterization of long non-coding RNAs expressed during sheep fetal and postnatal hair follicle development

Long non-coding RNAs (lncRNAs), >200 nt in length, are transcribed from mammalian genomes. They play important regulatory roles in various biological processes; However, the function and expression profile of lncRNAs involved in the development of hair follicles in the fetus, have been relatively under-explored area. To investigate the specific role of lncRNAs and mRNAs that regulate hair follicle development, we herein performed a comprehensive study on the lncRNA and mRNA expression profiles of sheep at multiple embryonic days (E65, E85, E105, and E135) and six lambs aged one week (D7) and one month (D30) using RNA-seq technology. The number of genes (471 lncRNAs and 12,812 mRNAs) differentially expressed and potential targets of differentially expressed lncRNAs were predicted. Differentially expressed lncRNAs were grouped into 10 clusters based on their expression pattern by K-means clustering. Moreover, Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that some differentially expressed mRNAs, such as DKK1, DSG4, FOXE1, Hoxc13, SFRP1, SFRP2, and Wnt10A overlapped with lncRNAs targets, and enriched in important hair follicle developmental pathways, including Wnt, TNF, and MAPK signaling pathways. In addition, 9 differentially expressed lncRNAs and 4 differentially expressed mRNAs were validated using quantitative real-time PCR (qRT-PCR). This study helps enrich the Ovis lncRNA databases and provides a comprehensive lncRNA transcriptome profile of fetal and postnatal skin of sheep. Additionally, it provides a foundation for further experiments on the role of lncRNAs in the regulation of hair growth in sheep.

Exploration of long noncoding RNA in bovine milk exosomes and their stability during digestion in vitro

Previous studies have demonstrated that bovine milk contains mRNA and microRNA that are largely encapsulated in milk-derived exosomes. However, little information is available about long noncoding RNAs (lncRNA) in bovine milk. Increasing evidence suggests that lncRNA are of particular interest given their key role in gene expression and development. We performed a comprehensive analysis of lncRNA in bovine milk exosomes by RNA sequencing. We used a validated human in vitro digestion model to investigate the stability of lncRNA encapsulated in bovine milk exosomes during the digestion process. We identified 3,475 novel lncRNA and 6 annotated lncRNA. The lncRNA shared characteristics with those of other mammals in terms of length, exon number, and open reading frames. However, lncRNA showed higher expression than mRNAs. We selected 12 lncRNA of high-expression abundance and identified them by PCR. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses showed that lncRNA regulate immune function, osteoblastogenesis, neurodevelopment, reproduction, cell proliferation, and cell-cell communication. We also investigated the 12 lncRNA using quantitative real-time PCR to reveal their expression profiles in milk exosomes during different stages of lactation (colostrum 2 d, 30 d, 150 d, and 270 d); their resulting expression levels in milk exosomes showed variations across the stages. A digestion experiment showed that bovine milk exosome lncRNA was resistant to in vitro digestion with different digestive juices, including saliva, gastric juice, pancreatic juice, and bile juice. Taken together, these results show for the first time that cow milk contains lncRNA, and that their abundance varied at different stages of lactation. As expected, bovine milk exosomal lncRNA were stable during in vitro digestion. These findings provide a basis for further understanding of the physiological role of milk lncRNA.

An Integrated Analysis of Cashmere Fineness lncRNAs in Cashmere Goats

Animal growth and development are regulated by long non-coding RNAs (lncRNAs). However, the functions of lncRNAs in regulating cashmere fineness are poorly understood. To identify the key lncRNAs that are related to cashmere fineness in skin, we have collected skin samples of Liaoning cashmere goats (LCG) and Inner Mongolia cashmere goats (MCG) in the anagen phase, and have performed RNA sequencing (RNA-seq) approach on these samples. The high-throughput sequencing and bioinformatics analyses identified 437 novel lncRNAs, including 93 differentially expressed lncRNAs. We also identified 3084 differentially expressed messenger RNAs (mRNAs) out of 27,947 mRNAs. Gene ontology (GO) analyses of lncRNAs and target genes in cis show a predominant enrichment of targets that are related to intermediate filament and intermediate filament cytoskeleton. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, sphingolipid metabolism is a significant pathway for lncRNA targets. In addition, this is the first report to reveal the possible lncRNA–mRNA regulatory network for cashmere fineness in cashmere goats. We also found that lncRNA XLOC_008679 and its target gene, KRT35, may be related to cashmere fineness in the anagen phase. The characterization and expression analyses of lncRNAs will facilitate future studies on the potential value of fiber development in LCG.

Analysis of long non-coding RNA expression profiles in neonatal rats with hypoxic-ischemic brain damage

Hypoxic-ischemic brain damage (HIBD) which is a common cause of acute mortality and neurological dysfunction in neonates still lacks effective therapeutic methods. Long non-coding RNAs (lncRNAs) were demonstrated to play a crucial role in many diseases. To give a foundation for subsequent functional studies of lncRNAs in HIBD, we investigated the profiling of lncRNAs and messenger RNAs (mRNAs) using neonatal HIBD rat model. Six neonatal rats were divided into sham-operated group (n = 3) and HIBD group (n = 3) randomly. Deep RNA sequencing was implemented to find out the meaningful lncRNAs and mRNAs. Quantitative real-time PCR was used to validate expressions of lncRNAs and mRNAs. The Gene Ontology (GO) and kyoto encyclopedia of genes a genomes (KEGG) database were used to predict functions of lncRNAs. A total of 328 differentially expressed lncRNAs (177 down-regulated vs 151 up-regulated) and 7157 differentially expressed mRNAs (2552 down-regulated vs 4605 up-regulated) were identified. The Quantitative real-time PCR results showed significant differential expressions of five lncRNAs and five mRNAs which were consistent with the RNA-Seq data. Gene ontology and KEGG analysis showed these lncRNAs and their expression-correlated mRNAs were closely related to the Janus tyrosine kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway, NF-kappa B signaling pathway, Toll-like receptor signaling pathway, calcium signaling pathway, Notch signaling pathway, mitogen activated protein kinase signaling pathway, neuroactive ligand-receptor interaction pathway and more. The results of our study identified the characterization and expression profiles of lncRNAs in neonatal HIBD and may be a basis for further therapeutic research. OPEN SCIENCE BADGES: This article has received a badge for *Open Materials* and *Open Data* because it provided all relevant information to reproduce the study in the manuscript and because it made the data publicly available. The data can be accessed at https://osf.io/yf3da/. The complete Open Science Disclosure form for this article can be found at the end of the article. More information about the Open Practices badges can be found at https://cos.io/our-services/open-science-badges/.