Integrated analysis of lncRNAs and mRNAs reveals key trans-target genes associated with ETEC-F4ac adhesion phenotype in porcine small intestine epithelial cells

Whole Genome Identification and Integrated Analysis of Long Non-Coding RNAs Responding ABA-Mediated Drought Stress in Panax ginseng C.A. Meyer

Panax ginseng C.A. Meyer is a perennial herb that is used worldwide for a number of medical purposes. Long non-coding RNAs (lncRNAs) play a crucial role in diverse biological processes but still remain poorly understood in ginseng, which has limited the application of molecular breeding in this plant. In this study, we identified 17,478 lncRNAs and 3106 novel mRNAs from ginseng by high-throughput illumine sequencing. 50 and 257 differentially expressed genes (DEGs) and DE lncRNAs (DELs) were detected under drought + ABA vs. drought conditions, respectively. The DEGs and DELs target genes main enrichment is focused on the "biosynthesis of secondary metabolites", "starch and sucrose metabolism", and "carbon metabolism" pathways under drought + ABA vs. drought conditions according to KEGG pathway enrichment analysis, suggesting that these secondary metabolites biosynthesis pathways might be crucial for ABA-mediated drought stress response in ginseng. Together, we identified drought stress response lncRNAs in ginseng for the first time and found that the target genes of these lncRNAs mainly regulate the biosynthesis of secondary metabolites pathway to response to drought stress. These findings also open up a new visual for molecular breeding in ginseng.

Transcriptomic analysis reveals role of lncRNA LOC100257036 to regulate AGAMOUS during cluster compactness of Vitis vinifera cv. sistan yaghooti

Yaghooti grape, as the earliest grape variety in Iran, is considered as more resistant to heat, drought, and salinity than other cultivars. Cluster compactness is regarded as an inappropriate feature for the productivity of Yaghooti grape as a critical commercial and nutritional product. In plants, lncRNAs play a critical role in regulating biological processes related to growth and development. However, the potential role of lncRNAs was not assessed in cluster compactness. Totally, 1549 lncRNAs were identified by RNA-Seq data analysis in three steps of cluster formation, berry formation, and final cluster size after a thorough screening process. In addition, 229 lncRNAs were differentially expressed in the cluster development steps. Based on the functional analysis, lncRNAs are related to AG and MYB, bHLH, LBD, NAC, and WRKY TFs. Further, the target genes enrichment analysis revealed a relationship between lncRNAs with grape growth and development, as well as resistance to abiotic stresses such as heat and drought, plant defense against pathogens, and early grapes ripening. The study identified four lncRNAs as precursors of miRNAs, predicting that 112 other lncRNAs could potentially be targeted by 166 miRNAs. The results provide new insights into the regulatory functions of lncRNAs in Yaghooti grape to improve overall understanding of the molecular mechanisms related to grape compactness.

ETENLNC: An end to end lncRNA identification and analysis framework to facilitate construction of known and novel lncRNA regulatory networks

Long non-coding RNAs (lncRNAs) play crucial roles in the regulation of gene expression and maintenance of genomic integrity through various interactions with DNA, RNA, and proteins. The availability of large-scale sequence data from various high-throughput platforms has opened possibilities to identify, predict, and functionally annotate lncRNAs. As a result, there is a growing demand for an integrative computational framework capable of identifying known lncRNAs, predicting novel lncRNAs, and inferring the downstream regulatory interactions of lncRNAs at the genome-scale. We present ETENLNC (End-To-End-Novel-Long-NonCoding), a user-friendly, integrative, open-source, scalable, and modular computational framework for identifying and analyzing lncRNAs from raw RNA-Seq data. ETENLNC employs six stringent filtration steps to identify novel lncRNAs, performs differential expression analysis of mRNA and lncRNA transcripts, and predicts regulatory interactions between lncRNAs, mRNAs, miRNAs, and proteins. We benchmarked ETENLNC against six existing tools and optimized it for desktop workstations and high-performance computing environments using data from three different species. ETENLNC is freely available on GitHub: https://github.com/EvolOMICS-TU/ETENLNC.

Transcriptome Analysis of Sesame (Sesamum indicum L.) Reveals the LncRNA and mRNA Regulatory Network Responding to Low Nitrogen Stress

Nitrogen is one of the important factors restricting the development of sesame planting and industry in China. Cultivating sesame varieties tolerant to low nitrogen is an effective way to solve the problem of crop nitrogen deficiency. To date, the mechanism of low nitrogen tolerance in sesame has not been elucidated at the transcriptional level. In this study, two sesame varieties Zhengzhi HL05 (ZZ, nitrogen efficient) and Burmese prolific (MD, nitrogen inefficient) in low nitrogen were used for RNA-sequencing. A total of 3964 DEGs (differentially expressed genes) and 221 DELs (differentially expressed lncRNAs) were identified in two sesame varieties at 3d and 9d after low nitrogen stress. Among them, 1227 genes related to low nitrogen tolerance are mainly located in amino acid metabolism, starch and sucrose metabolism and secondary metabolism, and participate in the process of transporter activity and antioxidant activity. In addition, a total of 209 pairs of lncRNA-mRNA were detected, including 21 pairs of trans and 188 cis. WGCNA (weighted gene co-expression network analysis) analysis divided the obtained genes into 29 modules; phenotypic association analysis identified three low-nitrogen response modules; through lncRNA-mRNA co-expression network, a number of hub genes and cis/trans-regulatory factors were identified in response to low-nitrogen stress including GS1-2 (glutamine synthetase 1-2), PAL (phenylalanine ammonia-lyase), CHS (chalcone synthase, CHS), CAB21 (chlorophyll a-b binding protein 21) and transcription factors MYB54, MYB88 and NAC75 and so on. As a trans regulator, lncRNA MSTRG.13854.1 affects the expression of some genes related to low nitrogen response by regulating the expression of MYB54, thus responding to low nitrogen stress. Our research is the first to provide a more comprehensive understanding of DEGs involved in the low nitrogen stress of sesame at the transcriptome level. These results may reveal insights into the molecular mechanisms of low nitrogen tolerance in sesame and provide diverse genetic resources involved in low nitrogen tolerance research.

Post-translational modifications of histone and non-histone proteins in epigenetic regulation and translational applications in alcohol-associated liver disease: Challenges and research opportunities

Epigenetic regulation is a process that takes place through adaptive cellular pathways influenced by environmental factors and metabolic changes to modulate gene activity with heritable phenotypic variations without altering the DNA sequences of many target genes. Epigenetic regulation can be facilitated by diverse mechanisms: many different types of post-translational modifications (PTMs) of histone and non-histone nuclear proteins, DNA methylation, altered levels of noncoding RNAs, incorporation of histone variants, nucleosomal positioning, chromatin remodeling, etc. These factors modulate chromatin structure and stability with or without the involvement of metabolic products, depending on the cellular context of target cells or environmental stimuli, such as intake of alcohol (ethanol) or Western-style high-fat diets. Alterations of epigenetics have been actively studied, since they are frequently associated with multiple disease states. Consequently, explorations of epigenetic regulation have recently shed light on the pathogenesis and progression of alcohol-associated disorders. In this review, we highlight the roles of various types of PTMs, including less-characterized modifications of nuclear histone and non-histone proteins, in the epigenetic regulation of alcohol-associated liver disease (ALD) and other disorders. We also describe challenges in characterizing specific PTMs and suggest future opportunities for basic and translational research to prevent or treat ALD and many other disease states.

Comprehensive analysis of circular RNAs in porcine small intestine epithelial cells associated with susceptibility to Escherichia coli F4ac diarrhea

BACKGROUND

Diarrhea is one of the most common diseases in pig industry, which seriously threatens the health of piglets and causes huge economic losses. Enterotoxigenic Escherichia coli (ETEC) F4 is regarded as the most important cause of diarrhea in piglets. Some pigs are naturally resistant to those diarrheas caused by ETEC-F4, because they have no F4 receptors (F4R) on their small intestine epithelial cells that allow F4 fimbriae adhesion. Circular RNA (circRNA) has been shown to play an important regulatory role in the pathogenesis of disease. We hypothesized that circRNAs may also regulate the adhesion of piglet small intestinal epithelial cells to ETEC F4 fimbriae. However, the circRNA expression profiles of piglets with different Enterotoxigenic Escherichia coli F4 fimbriae (ETEC-F4ac) adhesion phenotypes are still unclear, and the intermediate regulatory mechanisms need to be explored. Hence, the present study assessed the circRNA expression profiling in small intestine epithelial cells of eight male piglets with different ETEC-F4 adhesion phenotypes and ITGB5 genotypes to unravel their regulatory function in susceptibility to ETEC-F4ac diarrhea. Piglets were divided into two groups: non-adhesive group (n = 4) with CC genotype and adhesive group (n = 4) with TT genotype.

RESULTS

The RNA-seq data analysis identified 13,199 circRNAs from eight samples, most of which were exon-derived. In the small intestine epithelial cells, 305 were differentially expressed (DE) circRNAs between the adhesive and non-adhesive groups; of which 46 circRNAs were upregulated, and 259 were downregulated. Gene ontology and KEGG enrichment analysis revealed that most significantly enriched DE circRNAs' host genes were linked to cytoskeletal components, protein phosphorylation, cell adhesion, ion transport and pathways (such as adherens junction, gap junction) associated with ETEC diarrhea. The circRNA-miRNA-mRNA interaction network was also constructed to elucidate their underlying regulatory relationships. Our results identified several candidate circRNAs that affects susceptibility to ETEC diarrhea. Among them, circ-SORBS1 can adsorb ssc-miR-345-3p to regulate the expression of its host gene SORBS1, thus improving cell adhesion.

CONCLUSION

Our results provided insights into the regulation function of circRNAs in susceptibility to ETEC diarrhea of piglets, and enhanced our understanding of the role of circRNAs in regulating ETEC diarrhea, and reveal the great potential of circRNA as a diagnostic marker for susceptibility of ETEC diarrhea in piglets.

Exosomal long non-coding RNA TRAFD1-4:1 derived from fibroblast-like synoviocytes suppresses chondrocyte proliferation and migration by degrading cartilage extracellular matrix in rheumatoid arthritis

Rheumatoid arthritis (RA) is a chronic, autoimmune and systemic inflammatory disease affecting 1% of the population worldwide. Immune suppression of the activity and progress of RA is vital to reduce the disability and mortality rate as well as improve the quality of life of RA patients. However, the immune molecular mechanism of RA has not been clarified yet. Our results indicated that exosomes derived from TNFα-stimulated RA fibroblast-like synoviocytes (RA-FLSs) suppressed chondrocyte proliferation and migration through modulating cartilage extracellular matrix (CECM) determining by MTS assay, cell cycle analysis, Transwell assay and Western blot (WB). Besides, RNA sequencing and verification by qRT-PCR revealed that exosomal long non-coding RNA (lncRNA) tumor necrosis factor-associated factor 1 (TRAF1)-4:1 derived from RA-FLSs treated with TNFα was a candidate lncRNA, which also inhibited chondrocyte proliferation and migration through degrading CECM. Moreover, RNA sequencing and bioinformatics analysis identified that C-X-C motif chemokine ligand 1 (CXCL1) was a target mRNA of miR-27a-3p while miR-27a-3p was a target miRNA of lnc-TRAF1-4:1 in chondrocytes. Mechanistically, lnc-TRAF1-4:1 upregulated CXCL1 expression through sponging miR-27a-3p as a competing endogenous RNA (ceRNA) in chondrocytes identifying by Dual-luciferase reporter gene assay. Summarily, exosomal lncRNA TRAFD1-4:1 derived from RA-FLSs suppressed chondrocyte proliferation and migration through degrading CECM by upregulating CXCL1 as a sponge of miR-27a-3p. This study uncovered a novel RA-related lncRNA and investigated the roles of RA-FLS-derived exosomes and exosomal lnc-TRAF1-4:1 in articular cartilage impairment, which might provide novel therapeutic targets for RA.

RNA-seq reveals a novel porcine lncRNA MPHOSPH9-OT1 induces CXCL8/IL-8 expression in ETEC infected IPEC-J2 cells

Enterotoxigenic Escherichia coli (ETEC) is a major cause of bacterial diarrhea in piglets, leading to economic losses in the pig industry. In past decades, long non-coding RNAs (lncRNAs) have shown to be widely involved in the regulation of host immunity in porcine infection diseases. In this study, we explored the lncRNAs associated with ETEC F41 infection in IPEC-J2 cells by high-throughput sequencing and bioinformatic analysis. A total of 10150 novel porcine lncRNAs were identified. There were 161 differentially expressed (DE) lncRNAs associated with ETEC F41 infection, of which 65 DE lncRNAs were up-regulated and 96 DE lncRNAs were down-regulated. Functional and KEGG enrichment analysis of predicted target genes of DE lncRNAs indicated they are enriched in cell growth and inflammation-related pathways, such as endocytosis, focal adhesion, TGF-β signaling pathway, and adherens junctions. We revealed a novel candidate lncRNA MPHOSPH9-OT1 that was up-regulated after ETEC infection. The qRT-PCR validation and ELISA assessment showed the knockdown and overexpression of MPHOSPH9-OT1 resulted in significantly down- and up-regulation of cellular mRNA levels and secreted cytokine levels of CXCL8/IL-8, respectively. Meanwhile, MPHOSPH9-OT1 equilibrium is important to maintain the transepithelial electric resistance value and tight junction protein expression of IPEC-J2 cells. This study provides insights into the functionality of novel porcine lncRNAs in host immune responses to ETEC infection.

Identification and integrated analysis of lncRNAs and miRNAs in IPEC-J2 cells provide novel insight into the regulation of the innate immune response by PDCoV infection

Background

Noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), are pivotal regulators involved in the pathogenic mechanism of multiple coronaviruses. Porcine deltacoronavirus (PDCoV) has evolved multiple strategies to escape the innate immune response of host cells, but whether ncRNAs are involved in this process during PDCoV infection is still unknown.

Results

In this study, the expression profiles of miRNAs, lncRNAs and mRNAs in IPEC-J2 cells infected with PDCoV at 0, 12 and 24 hours postinfection (hpi) were identified through small RNA and RNA sequencing. The differentially expressed miRNAs (DEmiRNAs), lncRNAs (DElncRNAs) and mRNAs (DEmRNAs) were screened from the comparison group of IPEC-J2 cells at 0 and 12 hpi as well as the comparison group of IPEC-J2 cells at 12 and 24 hpi. The target genes of these DEncRNAs were predicted. The bioinformatics analysis of the target genes revealed multiple significantly enriched functions and pathways. Among them, the genes that were associated with innate immunity were specifically screened. The expression of innate immunity-related ncRNAs and mRNAs was validated by RT–qPCR. Competing endogenous RNA (ceRNA) regulatory networks among innate immunity-related ncRNAs and their target mRNAs were established. Moreover, we found that the replication of PDCoV was significantly inhibited by two innate immunity-related miRNAs, ssc-miR-30c-3p and ssc-miR-374b-3p, in IPEC-J2 cells.

Conclusions

This study provides a data platform to conduct studies of the pathogenic mechanism of PDCoV from a new perspective and will be helpful for further elucidation of the functional role of ncRNAs involved in PDCoV escaping the innate immune response.

Supplementary Information

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

Insights Into Long Non-Coding RNA and mRNA Expression in the Jejunum of Lambs Challenged With Escherichia coli F17

It has long been recognized that enterotoxigenic Escherichia coli (ETEC) is the major pathogen responsible for vomiting and diarrhea. E. coli F17, a main subtype of ETEC, is characterized by high morbidity and mortality in young livestock. However, the transcriptomic basis underlying E. coli F17 infection has not been fully understood. In the present study, RNA sequencing was conducted to explore the expression profiles of mRNAs and long non-coding RNAs (lncRNAs) in the jejunum of lambs who were identified as resistant or sensitive to E. coli F17 that was obtained in a challenge experiment. A total of 772 differentially expressed (DE) mRNAs and 190 DE lncRNAs were detected between the E. coli F17—resistance and E. coli F17-sensitive lambs (i.e., TFF2, LOC105606142, OLFM4, LYPD8, REG4, APOA4, TCONS_00223467, and TCONS_00241897). Then, a two-step machine learning approach (RX) combination Random Forest and Extreme Gradient Boosting were performed, which identified 16 mRNAs and 17 lncRNAs as potential biomarkers, within which PPP2R3A and TCONS_00182693 were prioritized as key biomarkers involved in E. coli F17 infection. Furthermore, functional enrichment analysis showed that peroxisome proliferator-activated receptor (PPAR) pathway was significantly enriched in response to E. coli F17 infection. Our finding will help to improve the knowledge of the mechanisms underlying E. coli F17 infection and may provide novel targets for future treatment of E. coli F17 infection.

Integrated analysis of long non-coding RNAs and mRNAs reveals the regulatory network of maize seedling root responding to salt stress

Background

Long non-coding RNAs (lncRNAs) play important roles in response to abiotic stresses in plants, by acting as cis- or trans-acting regulators of protein-coding genes. As a widely cultivated crop worldwide, maize is sensitive to salt stress particularly at the seedling stage. However, it is unclear how the expressions of protein-coding genes are affected by non-coding RNAs in maize responding to salt tolerance.

Results

The whole transcriptome sequencing was employed to investigate the differential lncRNAs and target transcripts responding to salt stress between two maize inbred lines with contrasting salt tolerance. We developed a flexible, user-friendly, and modular RNA analysis workflow, which facilitated the identification of lncRNAs and novel mRNAs from whole transcriptome data. Using the workflow, 12,817 lncRNAs and 8,320 novel mRNAs in maize seedling roots were identified and characterized. A total of 742 lncRNAs and 7,835 mRNAs were identified as salt stress-responsive transcripts. Moreover, we obtained 41 cis- and 81 trans-target mRNA for 88 of the lncRNAs. Among these target transcripts, 11 belonged to 7 transcription factor (TF) families including bHLH, C2H2, Hap3/NF-YB, HAS, MYB, WD40, and WRKY. The above 8,577 salt stress-responsive transcripts were further classified into 28 modules by weighted gene co-expression network analysis. In the salt-tolerant module, we constructed an interaction network containing 79 nodes and 3081 edges, which included 5 lncRNAs, 18 TFs and 56 functional transcripts (FTs). As a trans-acting regulator, the lncRNA MSTRG.8888.1 affected the expressions of some salt tolerance-relative FTs, including protein-serine/threonine phosphatase 2C and galactinol synthase 1, by regulating the expression of the bHLH TF.

Conclusions

The contrasting genetic backgrounds of the two inbred lines generated considerable variations in the expression abundance of lncRNAs and protein-coding transcripts. In the co-expression networks responding to salt stress, some TFs were targeted by the lncRNAs, which further regulated the salt tolerance-related functional transcripts. We constructed a regulatory pathway of maize seedlings to salt stress, which was mediated by the hub lncRNA MSTRG.8888.1 and participated by the bHLH TF and its downstream target transcripts. Future work will be focused on the functional revelation of the regulatory pathway.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08286-7.

N6-Methyladenosine Methylation Analysis of Long Noncoding RNAs and mRNAs in IPEC-J2 Cells Treated With Clostridium perfringens beta2 Toxin

Background

The n6-methyladenosine (m6A) modification is present widely in mRNAs and long non-coding RNAs (lncRNAs), and is related to the occurrence and development of certain diseases. However, the role of m6A methylation in Clostridium perfringens type C infectious diarrhea remains unclear.

Methods

Here, we treated intestinal porcine jejunum epithelial cells (IPEC-J2 cells) with Clostridium perfringens beta2 (CPB2) toxin to construct an in vitro model of Clostridium perfringens type C (C. perfringens type C) infectious diarrhea, and then used methylated RNA immunoprecipitation sequencing (MeRIP-seq) and RNA sequencing (RNA-seq) to identify the methylation profiles of mRNAs and lncRNAs in IPEC-J2 cells.

Results

We identified 6,413 peaks, representing 5,825 m6A-modified mRNAs and 433 modified lncRNAs, of which 4,356 m6A modified mRNAs and 221 m6A modified lncRNAs were significantly differential expressed between the control group and CPB2 group. The motif GGACU was enriched significantly in both the control group and the CPB2 group. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation analysis showed that the differentially methylated modified mRNAs were mainly enriched in Hippo signaling pathway and Wnt signaling pathway. In addition, the target genes of the differentially m6A modified lncRNAs were related to defense response to virus and immune response. For example, ENSSSCG00000042575, ENSSSCG00000048701 and ENSSSCG00000048785 might regulate the defense response to virus, immune and inflammatory response to resist the harmful effects of viruses on cells.

Conclusion

In summary, this study established the m6A transcription profile of mRNAs and lncRNAs in IPEC-J2 cells treated by CPB2 toxin. Further analysis showed that m6A-modified RNAs were related to defense against viruses and immune response after CPB2 toxin treatment of the cells. Threem6A-modified lncRNAs, ENSSSCG00000042575, ENSSSCG00000048785 and ENSSSCG00000048701, were most likely to play a key role in CPB2 toxin-treated IPEC-J2 cells. The results provide a theoretical basis for further research on the role of m6A modification in piglet diarrhea.