Search

Context: Ulcerative colitis has been clinically treated with Qing Hua Chang Yin (QHCY), a traditional Chinese medicine formula. However, its precise mechanisms in mitigating chronic colitis are largely uncharted. Objective: To elucidate the therapeutic efficiency of QHCY on chronic colitis and explore its underlying molecular mechanisms. Materials and methods: A total ion chromatogram fingerprint of QHCY was analysed. Chronic colitis was induced in male C57BL/6 mice using 2% dextran sodium sulphate (DSS) over 49 days. Mice were divided into control, DSS, DSS + QHCY (0.8, 1.6 and 3.2 g/kg/d dose, respectively) and DSS + mesalazine (0.2 g/kg/d) groups (n = 6). Mice were intragastrically administered QHCY or mesalazine for 49 days. The changes of disease activity index (DAI), colon length, colon histomorphology and serum pro-inflammatory factors in mice were observed. RNA sequencing was utilized to identify the differentially expressed transcripts (DETs) in colonic tissues and the associated signalling pathways. The expression of endoplasmic reticulum (ER) stress-related protein and NF-κB signalling pathway-related proteins in colonic tissues was detected by immunohistochemistry staining. Results: Forty-seven compounds were identified in QHCY. Compared with the DSS group, QHCY significantly improved symptoms of chronic colitis like DAI increase, weight loss, colon shortening and histological damage. It notably reduced serum levels of IL-6, IL-1β and TNF-α. QHCY suppressed the activation of PERK–ATF4–CHOP pathway of ER stress and NF-κB signalling pathways in colonic tissues. Discussion and conclusions: The findings in this study provide novel insights into the potential of QHCY in treating chronic colitis patients. [ABSTRACT FROM AUTHOR]

Abstract Objective Polymethylmethacrylate (PMMA) bone cement promotes the development of local thrombi. Our study found that a novel material, ES-PMMA bone cement, can reduce local thrombosis. We used a simple and reproducible animal model to confirm the reduction in local thrombosis and explored the associated molecular mechanism. Methods New Zealand rabbits, which were used to model thrombosis using extracorporeal carotid artery shunts, were divided into the following two groups, with 3 rabbits in each group: the PMMA bone cement group and the ES-PMMA bone cement group. Four hours after modelling, experimental samples, including thrombotic and vascular tissues, were collected. Thrombotic samples from the PMMA group and ES-PMMA group were subjected to lncRNA sequencing, and a lncRNA microarray was used to screen the differentially expressed lncRNAs. The expression of thrombomodulin in endothelial cells was quantified in vascular tissue samples. Differences in the lncRNA expression profiles between the thrombotic samples of the PMMA group and ES-PMMA group were assessed by base-to-base alignment in the intergenic regions of genomes. The lncRNA-miRNA-mRNA competitive endogenous RNA (ceRNA) network was established in light of ceRNA theory. Thrombosis was observed in the PMMA group and ES-PMMA group. Results The thrombotic weight was 0.00706 ± 0.00136 g/cm in the PMMA group and 0.00551 ± 0.00115 g/cm in the ES-PMMA group. Quantitative real-time polymerase chain reaction (RT–q-CR) and Western blotting revealed that the expression of CD40, which can regulate thrombosis in vascular endothelial cells, was significantly lower in the ES-PMMA group than in the PMMA group. High-throughput sequencing was used to identify 111 lncRNAs with lower expression in the ES-PMMA group than in the PMMA group. Through bioinformatics investigation, lncRNA MSTRG22719.16/ocu-miR-326-5p/CD40 binding sites were selected. Fluorescent in situ RNA hybridization (FISH) was performed to verify the lower expression of lncRNA MSTRG.22719.16 in vascular tissues from the ES-PMMA group. A dual-luciferase reporter gene assay was applied to verify that ocu-miR-326-5p binds the CD40 3ʹ-UTR and targets lncRNA MSTRG.22719.16. Conclusion Compared with PMMA bone cement, ES-PMMA bone cement can reduce thrombosis through the lncRNA MSTRG.22719.16/ocu-miR-326-5p/CD40 axis.

Long non-coding RNAs (lncRNAs) participate in acute lung injury (ALI). However, their latent biological function and molecular mechanism have not been fully understood. In the present study, the global expression profiles of lncRNAs and mRNAs between the control and lipopolysaccharide (LPS)-stimulated groups of human normal lung epithelial cells (BEAS-2B) were determined using high-throughput sequencing. Overall, a total of 433 lncRNAs and 183 mRNAs were differentially expressed. A lncRNA–mRNA co-expression network was established, and then the top 10 lncRNAs were screened using topological methods. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis results showed that the key lncRNAs targeting mRNAs were mostly enriched in the inflammatory-related biological processes. Gene set variation analysis and Pearson’s correlation coefficients confirmed the close correlation for the top 10 lncRNAs with inflammatory responses. A protein–protein interaction network analysis was conducted based on the key lncRNAs targeting mRNAs, where IL-1β, IL-6, and CXCL8 were regarded as the hub genes. A competing endogenous RNA (ceRNA) modulatory network was created with five lncRNAs, thirteen microRNAs, and twelve mRNAs. Finally, real-time quantitative reverse transcription-polymerase chain reaction was employed to verify the expression levels of several key lncRNAs in BEAS-2B cells and human serum samples.

Abstract Hemibarbus maculatus is a globally important species for its value in fisheries and aquaculture. There is relatively little information available about this species' sexual development genes. Therefore, in this study, we used the technology of high‐throughput RNA sequencing (RNA‐seq), to investigate if the ovary and testis' differently expressed genes (DEGs) during sex differentiation in H. maculatus. Our results revealed a total of 43,468,248 and 39,241,362 qualified reads were generated in the female and male libraries, and 32,127 unigenes were assembled via the annotation analysis. Moreover, 17,808 unigenes were found to be DEGs between males and females, with 12,726 up‐ and 5082 downregulated genes in the ovaries and testes, respectively. The transcriptome analysis revealed there are 23 DEGs between ovary and testes of H. maculatus, which may be involved in sex determination. The expression patterns of these genes were validated by quantitative real‐time PCR (qPCR). The results of qPCR were generally consistent with the RNA‐seq results. This work was the first to annotate genes linked with gonadal development and reproduction in H. maculatus using gonadal transcriptome analysis, which provides useful molecular data for a better understanding of sexual development, as well as functional genomics and population genetics studies of H. maculatus.

BackgroundSome bacteria enter the viable but non-culturable (VBNC) state to survive harsh environmental conditions and external stresses. This alters cell physiology and has implications for the food industry as some bacteria, such as lactobacilli, undergo similar changes during food processing.MethodsThis study aimed to investigate the transcriptomic changes of a probiotic strain, Lacticaseibacillus paracasei Zhang (L. paracasei Zhang), upon transition to the VBNC state using high throughput RNA sequencing (RNA-seq).ResultsBacteria were inoculated into the de Man, Rogosa, and Sharpe medium and maintained at low temperature and pH to induce cell transition to the VBNC state. Cells were harvested for analysis at five stages of VBNC induction: 0, 3, 30, and 180 days after induction and 210 days when the cells entered the VBNC state. Our results showed that the expression of 2,617, 2,642, 2,577, 2,829, and 2,840 genes was altered at these five different stages. The function of differentially expressed genes (DEGs, compared to healthy cells collected at day 0) and their encoded pathways were analyzed by the Gene Ontology Consortium and the Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. A total of 10 DEGs were identified in cells that entered the VBNC state: five continuously upregulated (LCAZH_0621, LCAZH_1986, LCAZH_2038, LCAZH_2040, and LCAZH_2174) and five continuously downregulated (LCAZH_0024, LCAZH_0210, LCAZH_0339, LCAZH_0621, and LCAZH_0754).ConclusionsThis study proposes a molecular model of the VBNC mechanism in L. paracasei Zhang, highlighting that changes in cell metabolism improve substrate utilization efficiency, thereby enhancing bacterial survival under adverse conditions. These data may be useful for improving the survival of probiotics in industrial food processing.

Helopeltis theivora or the tea mosquito bug (TMB) is reportedly one of the most devastating pests of tea plant (Camellia sinensis) causing threat to the beverage crop. Long non-coding RNAs (lncRNAs) constitute a group of endogenous RNAs that play gene regulatory roles in eukaryotes. In the present study, 9502 lncRNAs were identified from healthy and TMB-infested C. sinensis tissues using high-throughput strand-specific RNA sequencing, out of which 80 lncRNAs got differentially expressed in response to TMB infestation. Determination of genes that could act as potential targets of lncRNAs revealed that the identified lncRNAs could possibly target as many as 5804 genes. Differential gene expression (DGE) analysis led to the identification of 3665 differentially expressed genes (DEGs), of which, the expression of 1767 genes got upregulated and 1898 genes got downregulated during tea plant's response to TMB. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs and lncRNA-target genes have shown that TMB infestation might have triggered transcriptomic reprogramming leading to altered primary and secondary metabolism in C. sinensis. LncRNAs can act as competing endogenous RNAs (ceRNAs) to bind with common microRNA (miRNA) response elements (MREs) involving a competition between mRNAs and lncRNAs. We report 11 lncRNAs competing with 14 mRNAs to bind with 28 miRNAs forming the ceRNA network. The expression of 6 DEGs and 5 differentially expressed lncRNAs (DELs) has been validated by qRT-PCR. [ABSTRACT FROM AUTHOR]

Objective: Polymethylmethacrylate (PMMA) bone cement promotes the development of local thrombi. Our study found that a novel material, ES-PMMA bone cement, can reduce local thrombosis. We used a simple and reproducible animal model to confirm the reduction in local thrombosis and explored the associated molecular mechanism. Methods: New Zealand rabbits, which were used to model thrombosis using extracorporeal carotid artery shunts, were divided into the following two groups, with 3 rabbits in each group: the PMMA bone cement group and the ES-PMMA bone cement group. Four hours after modelling, experimental samples, including thrombotic and vascular tissues, were collected. Thrombotic samples from the PMMA group and ES-PMMA group were subjected to lncRNA sequencing, and a lncRNA microarray was used to screen the differentially expressed lncRNAs. The expression of thrombomodulin in endothelial cells was quantified in vascular tissue samples. Differences in the lncRNA expression profiles between the thrombotic samples of the PMMA group and ES-PMMA group were assessed by base-to-base alignment in the intergenic regions of genomes. The lncRNA-miRNA-mRNA competitive endogenous RNA (ceRNA) network was established in light of ceRNA theory. Thrombosis was observed in the PMMA group and ES-PMMA group. Results: The thrombotic weight was 0.00706 ± 0.00136 g/cm in the PMMA group and 0.00551 ± 0.00115 g/cm in the ES-PMMA group. Quantitative real-time polymerase chain reaction (RT–q-CR) and Western blotting revealed that the expression of CD40, which can regulate thrombosis in vascular endothelial cells, was significantly lower in the ES-PMMA group than in the PMMA group. High-throughput sequencing was used to identify 111 lncRNAs with lower expression in the ES-PMMA group than in the PMMA group. Through bioinformatics investigation, lncRNA MSTRG22719.16/ocu-miR-326-5p/CD40 binding sites were selected. Fluorescent in situ RNA hybridization (FISH) was performed to verify the lower expression of lncRNA MSTRG.22719.16 in vascular tissues from the ES-PMMA group. A dual-luciferase reporter gene assay was applied to verify that ocu-miR-326-5p binds the CD40 3ʹ-UTR and targets lncRNA MSTRG.22719.16. Conclusion: Compared with PMMA bone cement, ES-PMMA bone cement can reduce thrombosis through the lncRNA MSTRG.22719.16/ocu-miR-326-5p/CD40 axis. [ABSTRACT FROM AUTHOR]

Background: Highly pathogenic avian influenza virus (HPAIV) is considered a global threat to both human health and the poultry industry. MicroRNAs (miRNA) can modulate the immune system by affecting gene expression patterns in HPAIV-infected chickens. Objectives: To gain further insights into the role of miRNAs in immune responses against H5N1 infection, as well as the development of strategies for breeding disease-resistant chickens, we characterized miRNA expression patterns in tracheal tissues from H5N1-infected Ri chickens. Methods: miRNAs expression was analyzed from two H5N1-infected Ri chicken lines using small RNA sequencing. The target genes of differentially expressed (DE) miRNAs were predicted using miRDB. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis were then conducted. Furthermore, using quantitative real-time polymerase chain reaction, we validated the expression levels of DE miRNAs (miR-22-3p, miR-146b-3p, miR-27b-3p, miR-128-3p, miR-2188-5p, miR-451, miR-205a, miR-203a, miR-21-3p, and miR-200a-3p) from all comparisons and their immune-related target genes. Results: A total of 53 miRNAs were significantly expressed in the infection samples of the resistant compared to the susceptible line. Network analyses between the DE miRNAs and target genes revealed that DE miRNAs may regulate the expression of target genes involved in the transforming growth factor-beta, mitogen-activated protein kinase, and Toll-like receptor signaling pathways, all of which are related to influenza A virus progression. Conclusions: Collectively, our results provided novel insights into the miRNA expression patterns of tracheal tissues from H5N1-infected Ri chickens. More importantly, our findings offer insights into the relationship between miRNA and immune-related target genes and the role of miRNA in HPAIV infections in chickens. [ABSTRACT FROM AUTHOR]

Hemibarbus maculatus is a globally important species for its value in fisheries and aquaculture. There is relatively little information available about this species' sexual development genes. Therefore, in this study, we used the technology of high‐throughput RNA sequencing (RNA‐seq), to investigate if the ovary and testis' differently expressed genes (DEGs) during sex differentiation in H. maculatus. Our results revealed a total of 43,468,248 and 39,241,362 qualified reads were generated in the female and male libraries, and 32,127 unigenes were assembled via the annotation analysis. Moreover, 17,808 unigenes were found to be DEGs between males and females, with 12,726 up‐ and 5082 downregulated genes in the ovaries and testes, respectively. The transcriptome analysis revealed there are 23 DEGs between ovary and testes of H. maculatus, which may be involved in sex determination. The expression patterns of these genes were validated by quantitative real‐time PCR (qPCR). The results of qPCR were generally consistent with the RNA‐seq results. This work was the first to annotate genes linked with gonadal development and reproduction in H. maculatus using gonadal transcriptome analysis, which provides useful molecular data for a better understanding of sexual development, as well as functional genomics and population genetics studies of H. maculatus. [ABSTRACT FROM AUTHOR]

Abstract Background Long non-coding RNA (lncRNA) plays crucial role in osteogenic differentiation of bone marrow mesenchymal stem cells (BMMSCs), involving in regulation of competing endogenous RNA (ceRNA) mechanisms and conduction of signaling pathways. However, its mechanisms are poorly understood. This study aimed to investigate lncRNAs, miRNAs and mRNAs expression profiles in rat BMMSCs (rBMMSCs) osteogenic differentiation, screen the potential key lncRNA-miRNA-mRNA networks, explore the putative functions and identify the key molecules, as the basis of studying potential mechanism of rBMMSCs osteogenic differentiation driven by lncRNA, providing molecular targets for the management of bone defect. Methods High-throughput RNA sequencing (RNA-seq) was used to determine lncRNAs, miRNAs, and mRNAs expression profiles at 14-day rBMMSCs osteogenesis. The pivotal lncRNA-miRNA and miRNA-mRNA networks were predicted from sequencing data and bioinformatic analysis, and the results were exported by Cytoscape 3.9.0 software. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used for functional exploration. Real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was performed to validate lncRNAs, miRNAs and mRNAs. Results rBMMSCs were identified, and the osteogenic and adipogenic differentiation ability were detected. A total of 8634 lncRNAs were detected by RNA-seq, and 1524 differential expressed lncRNAs, of which 812 up-regulated and 712 down-regulated in osteo-inductive groups compared with control groups. 30 up-regulated and 61 down-regulated miRNAs, 91 miRNAs were differentially expressed in total. 2453 differentially expressed mRNAs including 1272 up-expressed and 1181 down-expressed were detected. 10 up-regulated lncRNAs were chosen to predict 21 down-regulated miRNAs and 650 up-regulated mRNAs. 49 lncRNA-miRNA and 1515 miRNA–mRNA interactive networks were constructed. GO analysis showed the most important enrichment in cell component and molecular function were “cytoplasm” and “protein binding”, respectively. Biological process related to osteogenic differentiation such as “cell proliferation”, “wound healing”, “cell migration”, “osteoblast differentiation”, “extracellular matrix organization” and “response to hypoxia” were enriched. KEGG analysis showed differentially expressed genes were mainly enriched in “PI3K-Akt signaling pathway”, “Signaling pathway regulating pluripotency of stem cells”, “cGMP-PKG signaling pathway”, “Axon guidance” and “Calcium signaling pathway”. qRT-PCR verified that lncRNA Tug1, lncRNA AABR07011996.1, rno-miR-93-5p, rno-miR-322-5p, Sgk1 and Fzd4 were consistent with the sequencing results, and 4 lncRNA-miRNA-mRNA networks based on validations were constructed, and enrichment pathways were closely related to “PI3K-Akt signaling pathway”, “Signaling pathway regulating pluripotency of stem cells” and “Wnt signaling pathway”. Conclusions lncRNAs, miRNAs and mRNAs expression profiles provide clues for future studies on their roles for BMMSCs osteogenic differentiation. Furthermore, lncRNA–miRNA–mRNA networks give more information on potential new mechanisms and targets for management on bone defect.

Background. Numerous studies have shown circular RNA (circRNA) dysregulation is associated with the pathogenesis of cervical cancer, particularly in individual carcinoma variants. The aim of this study is to investigate and contrastively analyze the expression pattern of circRNAs in cervical squamous carcinoma and adenocarcinoma mediated by human papillomavirus type 16 (HPV-16). Methods. The expression of circRNAs in cervical squamous carcinoma (SCC), adeno-carcinoma (ADC) and adenosquamous carcinoma (ASC) tissues, together with the adjacent normal tissues (ANT), was profiled by high-throughput RNA sequencing (RNA-seq). Bioinformatics analysis and quantitative real time polymerase chain reaction (qRT-PCR) validation of the sequencing data were performed. A network of circRNA-miRNA (microRNA)-mRNA was then constructed according to predicted targets and function of candidate circRNAs. Results. A total of 11,685 annotated circRNAs were identified in six cervical samples. There were 42 up-regulated and 98 down-regulated circRNAs. 215 circRNAs were up-regulated in SCC but down-regulated circRNAs in ADC, while 50 circRNAs displayed the opposite trend. Function enrichment analysis based on different expressions of circRNAs found that the most enriched pathway in all the three pathologic variants of cervical cancer was the "ubiquitin mediated proteolysis" pathway. Eight key candidate circRNAs derived from this pathway were further validated, and we noticed that several target miRNAs of candidate circRNAs could target the source genes. Based on this we constructed a related competing endogenous RNA (ceRNA) network. Conclusion. Through a comprehensive interpretation of differentially expressed circRNAs in different pathologic variants of cervical cancer, this study provides new insights into the process of tumor differentiation mediated by HPV. Our results may help to complement the molecular typing and stem cell theory of cervical cancer. [ABSTRACT FROM AUTHOR]

Postoperative atrial fibrillation (POAF) is a frequent complication associated with increased periprocedural mortality and morbidity after cardiac surgery. Our study aimed to identify the difference in exosomal miRNA and further explore its role in the diagnosis of POAF. First, the differentially expressed miRNAs (DEMs) were obtained by high-throughput RNA sequencing. Second, the DEMs target genes were put into gene ontology (GO) and KEGG pathway analysis. Third, real-time quantification PCR (RT-qPCR) was used to verify the DEMs. Finally, we revealed 23 DEMs in POAF patients. Furthermore, analysis of gene function revealed that DEMs may affect atrial structure through many signaling pathways. We also found that miR-122-5p was up-regulated in POAF patients, but there are no significant changes in miR-191-5p, miR-181a-5p, miR-155-5p and miR-151a-5p. Our study revealed that exosomal miRNAs exert enormous potential in evaluating the severity or prognostic of POAF. [ABSTRACT FROM AUTHOR]

The holocrine process of sebum, derived from the sebocytes, begins with the proliferation of this undifferentiated and mitotically active population of peripheral, immature cells [[10]]; furthermore, the cells undergo nuclear degeneration and move more centrally, which can lead to the induction of apoptosis and the bursting of cells [[12]]. TdT-Mediated dUTP Nick-End Labeling (TUNEL) Assay for Cell Apoptosis The cell apoptotic analysis was performed on three slides randomly selected to be technical replicates from the E36, E41, and NH1 using commercial TUNEL assay kit (C1089, Beyotime, China). As the most important regulators of cell apoptosis, the key regulator of B-cell lymphoma 2 (Bcl-2) shows expression patterns in the uropygial gland; this indicates that the terminal differentiation of sebocytes during holocrine secretion in DGs is realized mainly through cell apoptosis [[13]]. For the cell apoptosis pathway, five genes including MCP1, -tubulin 3, CRLM2, CSF2RB, and BIRC5, were down-regulated in the E41 compared with those within the E36 ( I p i < 0.01), whereas no significant difference was detected between E41 and NH1 ( I p i > 0.05). [Extracted from the article]

Copyright of Journal of Henan Agricultural Sciences is the property of Editorial Board of Journal of Henan Agricultural Sciences and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Toxoplasmosis is an important zoonotic parasitic disease caused by Toxoplasma gondii (T. gondii). However, the functions of circRNAs and miRNAs in response to T. gondii infection in the livers of mice at acute and chronic stages remain unknown. Here, high-throughput RNA sequencing was performed for detecting the expression of circRNAs and miRNAs in livers of mice infected with 20 T. gondii cysts at the acute and chronic stages, in order to understand the potential molecular mechanisms underlying hepatic toxoplasmosis. Overall, 265 and 97 differentially expressed (DE) circRNAs were found in livers at the acute and chronic infection stages in comparison with controls, respectively. In addition, 171 and 77 DEmiRNAs were found in livers at the acute and chronic infection stages, respectively. Functional annotation showed that some immunity-related Gene ontology terms, such as “positive regulation of cytokine production”, “regulation of T cell activation”, and “immune receptor activity”, were enriched at the two infection stages. Moreover, the pathways “Valine, leucine, and isoleucine degradation”, “Fatty acid metabolism”, and “Glycine, serine, and threonine metabolism” were involved in liver disease. Remarkably, DEcircRNA 6:124519352|124575359 was significantly correlated with DEmiRNAs mmu-miR-146a-5p and mmu-miR-150-5p in the network that was associated with liver immunity and pathogenesis of disease. This study revealed that the expression profiling of circRNAs in the livers was changed after T. gondii infection, and improved our understanding of the transcriptomic landscape of hepatic toxoplasmosis in mice.

Abstract Background Air pollution may be associated with increased airway responsiveness to allergens in allergic rhinitis (AR). Ozone-aged environmental black carbon (O3BC) is an important constituent of atmospheric particulate matter (PM), for which the mechanisms underlying its effects have not been fully elucidated in AR. The objective of the present study was to determine the O3BC and pollen-induced alterations in the transcriptome in human nasal epithelial cells (hNECs) in vitro. Methods hNECs from nasal epithelial mucosal samples of healthy individuals undergoing nasal surgery (turbinoplasty or septoplasty) were established as air–liquid interface (ALI) cultures and exposed to O3BC, pollen, or a combination of O3BC+ pollen. Changes in cell viability were analyzed by fluorescence and changes in the transcriptome by high-throughput RNA sequencing (RNA-seq). Several differentially expressed genes were verified by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Enrichment analysis, based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database, was performed to determine major biological functions and pathways involved. Results Exposure to ≥ 50 μg/ml O3BC or 25 μg/ml O3BC+ 200 μg /ml pollen significantly decreased cell viability of the hNECs compared to control (p

Background: Long non-coding RNA (lncRNA) plays crucial role in osteogenic differentiation of bone marrow mesenchymal stem cells (BMMSCs), involving in regulation of competing endogenous RNA (ceRNA) mechanisms and conduction of signaling pathways. However, its mechanisms are poorly understood. This study aimed to investigate lncRNAs, miRNAs and mRNAs expression profiles in rat BMMSCs (rBMMSCs) osteogenic differentiation, screen the potential key lncRNA-miRNA-mRNA networks, explore the putative functions and identify the key molecules, as the basis of studying potential mechanism of rBMMSCs osteogenic differentiation driven by lncRNA, providing molecular targets for the management of bone defect. Methods: High-throughput RNA sequencing (RNA-seq) was used to determine lncRNAs, miRNAs, and mRNAs expression profiles at 14-day rBMMSCs osteogenesis. The pivotal lncRNA-miRNA and miRNA-mRNA networks were predicted from sequencing data and bioinformatic analysis, and the results were exported by Cytoscape 3.9.0 software. Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were used for functional exploration. Real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was performed to validate lncRNAs, miRNAs and mRNAs. Results: rBMMSCs were identified, and the osteogenic and adipogenic differentiation ability were detected. A total of 8634 lncRNAs were detected by RNA-seq, and 1524 differential expressed lncRNAs, of which 812 up-regulated and 712 down-regulated in osteo-inductive groups compared with control groups. 30 up-regulated and 61 down-regulated miRNAs, 91 miRNAs were differentially expressed in total. 2453 differentially expressed mRNAs including 1272 up-expressed and 1181 down-expressed were detected. 10 up-regulated lncRNAs were chosen to predict 21 down-regulated miRNAs and 650 up-regulated mRNAs. 49 lncRNA-miRNA and 1515 miRNA–mRNA interactive networks were constructed. GO analysis showed the most important enrichment in cell component and molecular function were "cytoplasm" and "protein binding", respectively. Biological process related to osteogenic differentiation such as "cell proliferation", "wound healing", "cell migration", "osteoblast differentiation", "extracellular matrix organization" and "response to hypoxia" were enriched. KEGG analysis showed differentially expressed genes were mainly enriched in "PI3K-Akt signaling pathway", "Signaling pathway regulating pluripotency of stem cells", "cGMP-PKG signaling pathway", "Axon guidance" and "Calcium signaling pathway". qRT-PCR verified that lncRNA Tug1, lncRNA AABR07011996.1, rno-miR-93-5p, rno-miR-322-5p, Sgk1 and Fzd4 were consistent with the sequencing results, and 4 lncRNA-miRNA-mRNA networks based on validations were constructed, and enrichment pathways were closely related to "PI3K-Akt signaling pathway", "Signaling pathway regulating pluripotency of stem cells" and "Wnt signaling pathway". Conclusions: lncRNAs, miRNAs and mRNAs expression profiles provide clues for future studies on their roles for BMMSCs osteogenic differentiation. Furthermore, lncRNA–miRNA–mRNA networks give more information on potential new mechanisms and targets for management on bone defect. [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) play an important role in growth, development, stress resilience, and epigenetic modifications of plants. However, the effect of calcium (Ca2+) deficiency on miRNA expression in the orphan crop tef (Eragrostis tef) remains unknown. In this study, we analyzed expression of miRNAs in roots and shoots of tef in response to Ca2+ treatment. miRNA‐seq followed by bioinformatic analysis allowed us to identify a large number of small RNAs (sRNAs) ranging from 17 to 35 nt in length. A total of 1380 miRNAs were identified in tef experiencing long‐term Ca2+ deficiency while 1495 miRNAs were detected in control plants. Among the miRNAs identified in this study, 161 miRNAs were similar with those previously characterized in other plant species and 348 miRNAs were novel, while the remaining miRNAs were uncharacterized. Putative target genes and their functions were predicted for all the known and novel miRNAs that we identified. Based on gene ontology (GO) analysis, the predicted target genes are known to have various biological and molecular functions including calcium uptake and transport. Pairwise comparison of differentially expressed miRNAs revealed that some miRNAs were specifically enriched in roots or shoots of low Ca2+‐treated plants. Further characterization of the miRNAs and their targets identified in this study may help in understanding Ca2+ deficiency responses in tef and related orphan crops. [ABSTRACT FROM AUTHOR]

Abstract MicroRNAs (miRNAs) play an important role in growth, development, stress resilience, and epigenetic modifications of plants. However, the effect of calcium (Ca2+) deficiency on miRNA expression in the orphan crop tef (Eragrostis tef) remains unknown. In this study, we analyzed expression of miRNAs in roots and shoots of tef in response to Ca2+ treatment. miRNA‐seq followed by bioinformatic analysis allowed us to identify a large number of small RNAs (sRNAs) ranging from 17 to 35 nt in length. A total of 1380 miRNAs were identified in tef experiencing long‐term Ca2+ deficiency while 1495 miRNAs were detected in control plants. Among the miRNAs identified in this study, 161 miRNAs were similar with those previously characterized in other plant species and 348 miRNAs were novel, while the remaining miRNAs were uncharacterized. Putative target genes and their functions were predicted for all the known and novel miRNAs that we identified. Based on gene ontology (GO) analysis, the predicted target genes are known to have various biological and molecular functions including calcium uptake and transport. Pairwise comparison of differentially expressed miRNAs revealed that some miRNAs were specifically enriched in roots or shoots of low Ca2+‐treated plants. Further characterization of the miRNAs and their targets identified in this study may help in understanding Ca2+ deficiency responses in tef and related orphan crops.

Abstract Growing evidence has indicated that circular RNAs (circRNAs) play crucial roles in multiple biological processes. However, alterations in circRNA profiles during bladder cancer progression and the clinical significance thereof remain unclear. Therefore, high‐throughput RNA sequencing was conducted to identify circRNA and mRNA profiles in five pairs of bladder cancer tissues and adjacent noncancerous tissues. A total of 87 differentially expressed circRNAs and 2756 mRNAs were detected in above bladder cancer samples compared with paired noncancerous samples. Functional enrichment analyses, circRNA‐microRNA‐mRNA, and protein‐protein interaction networks revealed that these dysregulated circRNAs were potentially involved in carcinogenesis and evolution of bladder cancer. Subsequently, the differential expression of eight circRNAs was detected by real‐time qPCR. Hsa_circ_0003141 and hsa_circ_0008039 were significantly upregulated as well as hsa_circ_0026782, hsa_circ_0077837, hsa_circ_0004826, and hsa_circ_0001946 were significantly downregulated among validation of 70 matched bladder cancer tissues (≥75%). Moreover, hsa_circ_0077837 and hsa_circ_0004826 were also verified as markedly downregulated in four bladder cancer cells (100%). Naturally, hsa_circ_0077837 and hsa_circ_0004826 were also demonstrated using RNase‐R+ resistance experiments. In addition, Fisherʹs exact test, Kaplan‐Meier plots, Cox regression analyses, and receiver operating characteristic curve was performed to assess their clinical value. Downregulation of hsa_circ_0077837 and hsa_circ_0004826 all was significantly correlated with worse clinicopathological features and poor prognosis of bladder cancer patients. The area under the receiver operating characteristic curve of them was 0.775 (P

Abstract Bone cancer pain (BCP) is one of the most common types of chronic cancer pain and its pathogenesis has not been fully understood. Long non-coding RNAs (lncRNAs) are new promising targets in the field of pain research, however, their involvements in BCP have not been reported. In the present study, we established the BCP model by implantation of Walker 256 carcinoma cells into rats’ tibial medullary cavity and performed transcriptome sequencing of the ipsilateral lumbar spinal cord to explore changes in expression profiles of lncRNA and mRNA. We identified 1220 differently expressed mRNAs (DEmRNAs) (1171 up-regulated and 49 down-regulated) and 323 differently expressed lncRNAs (DElncRNAs) (246 up-regulated and 77 down-regulated) in BCP model, among which 10 DEmRNAs (5 up-regulated and 5 down-regulated) and 10 DElncRNAs (5 up-regulated and 5 down-regulated) were validated the expression by RT-qPCR. Then, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis on the expression of DEmRNAs and DElncRNAs, showing that they were mainly enriched in inflammatory and immunologic processes/pathways. Finally, we constructed a co-expression network and a ceRNA network of DEmRNAs and DElncRNAs to exhibit a potential regulatory mechanism of DElncRNAs, directly regulating protein coding gene expression in cis or in trans and indirectly regulating protein coding gene expression by sponging miRNA. In conclusion, our study provided a landscape of dysregulated lncRNA and mRNA in spinal cord of bone cancer pain and detected novel potential targets for treatment in the future.

Non-coding RNAs have remarkable roles in acute lung injury (ALI) initiation. Nevertheless, the significance of long non-coding RNAs (lncRNAs) in ALI is still unknown. Herein, we purposed to identify potential key genes in ALI and create a competitive endogenous RNA (ceRNA) modulatory network to uncover possible molecular mechanisms that affect lung injury. We generated a lipopolysaccharide-triggered ALI mouse model, whose lung tissue was subjected to RNA sequencing, and then we conducted bioinformatics analysis to select genes showing differential expression (DE) and to build a lncRNA-miRNA (microRNA)- mRNA (messenger RNA) modulatory network. Besides, GO along with KEGG assessments were conducted to identify major biological processes and pathways, respectively, involved in ALI. Then, RT-qPCR assay was employed to verify levels of major RNAs. A protein-protein interaction (PPI) network was created using the Search Tool for the Retrieval of Interacting Genes (STRING) database, and the hub genes were obtained with the Molecular Complex Detection plugin. Finally, a key ceRNA subnetwork was built from major genes and their docking sites. Overall, a total of 8,610 lncRNAs were identified in the normal and LPS groups. Based on the 308 DE lncRNAs [ p -value < 0.05, |log2 (fold change) | > 1] and 3,357 DE mRNAs [ p -value < 0.05, |log2 (fold change) | > 1], lncRNA-miRNA and miRNA-mRNA pairs were predicted using miRanda. The lncRNA-miRNA-mRNA network was created from 175 lncRNAs, 22 miRNAs, and 209 mRNAs in ALI. The RT-qPCR data keep in step with the RNA sequencing data. GO along with KEGG analyses illustrated that DE mRNAs in this network were mainly bound up with the inflammatory response, developmental process, cell differentiation, cell proliferation, apoptosis, and the NF-kappa B, PI3K-Akt, HIF-1, MAPK, Jak-STAT, and Notch signaling pathways. A PPI network on the basis of the 209 genes was established, and three hub genes (Nkx2-1, Tbx2, and Atf5) were obtained from the network. Additionally, a lncRNA-miRNA-hub gene subnetwork was built from 15 lncRNAs, 3 miRNAs, and 3 mRNAs. Herein, novel ideas are presented to expand our knowledge on the regulation mechanisms of lncRNA-related ceRNAs in the pathogenesis of ALI. [ABSTRACT FROM AUTHOR]

RationaleTreatment of sinus node disease with regenerative or cell-based therapies will require a detailed understanding of gene regulatory networks in cardiac pacemaker cells (PCs).ObjectiveTo characterize the transcriptome of PCs using RNA sequencing and to identify transcriptional networks responsible for PC gene expression.Methods and resultsWe used laser capture microdissection on a sinus node reporter mouse line to isolate RNA from PCs for RNA sequencing. Differential expression and network analysis identified novel sinoatrial node-enriched genes and predicted that the transcription factor Islet-1 is active in developing PCs. RNA sequencing on sinoatrial node tissue lacking Islet-1 established that Islet-1 is an important transcriptional regulator within the developing sinoatrial node.Conclusions(1) The PC transcriptome diverges sharply from other cardiomyocytes; (2) Islet-1 is a positive transcriptional regulator of the PC gene expression program.

Non-coding RNAs have remarkable roles in acute lung injury (ALI) initiation. Nevertheless, the significance of long non-coding RNAs (lncRNAs) in ALI is still unknown. Herein, we purposed to identify potential key genes in ALI and create a competitive endogenous RNA (ceRNA) modulatory network to uncover possible molecular mechanisms that affect lung injury. We generated a lipopolysaccharide-triggered ALI mouse model, whose lung tissue was subjected to RNA sequencing, and then we conducted bioinformatics analysis to select genes showing differential expression (DE) and to build a lncRNA-miRNA (microRNA)- mRNA (messenger RNA) modulatory network. Besides, GO along with KEGG assessments were conducted to identify major biological processes and pathways, respectively, involved in ALI. Then, RT-qPCR assay was employed to verify levels of major RNAs. A protein-protein interaction (PPI) network was created using the Search Tool for the Retrieval of Interacting Genes (STRING) database, and the hub genes were obtained with the Molecular Complex Detection plugin. Finally, a key ceRNA subnetwork was built from major genes and their docking sites. Overall, a total of 8,610 lncRNAs were identified in the normal and LPS groups. Based on the 308 DE lncRNAs [p-value < 0.05, |log2 (fold change) | > 1] and 3,357 DE mRNAs [p-value < 0.05, |log2 (fold change) | > 1], lncRNA-miRNA and miRNA-mRNA pairs were predicted using miRanda. The lncRNA-miRNA-mRNA network was created from 175 lncRNAs, 22 miRNAs, and 209 mRNAs in ALI. The RT-qPCR data keep in step with the RNA sequencing data. GO along with KEGG analyses illustrated that DE mRNAs in this network were mainly bound up with the inflammatory response, developmental process, cell differentiation, cell proliferation, apoptosis, and the NF-kappa B, PI3K-Akt, HIF-1, MAPK, Jak-STAT, and Notch signaling pathways. A PPI network on the basis of the 209 genes was established, and three hub genes (Nkx2-1, Tbx2, and Atf5) were obtained from the network. Additionally, a lncRNA-miRNA-hub gene subnetwork was built from 15 lncRNAs, 3 miRNAs, and 3 mRNAs. Herein, novel ideas are presented to expand our knowledge on the regulation mechanisms of lncRNA-related ceRNAs in the pathogenesis of ALI.

Brown adipose tissue (BAT) is specialized for energy expenditure, thus a better understanding of the regulators influencing BAT development could provide novel strategies to defense obesity. Many protein-coding genes, miRNAs, and lncRNAs have been investigated in BAT development, however, the expression patterns and functions of circRNA in brown adipogenesis have not been reported yet. This study determined the circRNA expression profiles across brown adipogenesis (proliferation, early differentiated, and fully differentiated stages) by RNA-seq. We identified 3,869 circRNAs and 36.9% of them were novel. We found the biogenesis of circRNA was significantly related to linear mRNA transcription, meanwhile, almost 70% of circRNAs were generated by alternative back-splicing. Next, we examined the cell-specific and differentiation stage-specific expression of circRNAs. Compared to white adipocytes, nearly 30% of them were specifically expressed in brown adipocytes. Further, time-series expression analysis showed circRNAs were dynamically expressed, and 117 differential expression circRNAs (DECs) in brown adipogenesis were identified, with 77 upregulated and 40 downregulated. Experimental validation showed the identified circRNAs could be successfully amplified and the expression levels detected by RNA-seq were reliable. For the potential functions of the circRNAs, GO analysis suggested that the decreased circRNAs were enriched in cell proliferation terms, while the increased circRNAs were enriched in development and thermogenic terms. Bioinformatics predictions showed that DECs contained numerous binding sites of functional miRNAs. More interestingly, most of the circRNAs contained multiple binding sites for the same miRNA, indicating that they may facilitate functions by acting as microRNA sponges. Collectively, we characterized the circRNA expression profiles during brown adipogenesis and provide numerous novel circRNAs candidates for future brown adipogenesis regulating studies.

Brown adipose tissue (BAT) is specialized for energy expenditure, thus a better understanding of the regulators influencing BAT development could provide novel strategies to defense obesity. Many protein-coding genes, miRNAs, and lncRNAs have been investigated in BAT development, however, the expression patterns and functions of circRNA in brown adipogenesis have not been reported yet. This study determined the circRNA expression profiles across brown adipogenesis (proliferation, early differentiated, and fully differentiated stages) by RNA-seq. We identified 3,869 circRNAs and 36.9% of them were novel. We found the biogenesis of circRNA was significantly related to linear mRNA transcription, meanwhile, almost 70% of circRNAs were generated by alternative back-splicing. Next, we examined the cell-specific and differentiation stage-specific expression of circRNAs. Compared to white adipocytes, nearly 30% of them were specifically expressed in brown adipocytes. Further, time-series expression analysis showed circRNAs were dynamically expressed, and 117 differential expression circRNAs (DECs) in brown adipogenesis were identified, with 77 upregulated and 40 downregulated. Experimental validation showed the identified circRNAs could be successfully amplified and the expression levels detected by RNA-seq were reliable. For the potential functions of the circRNAs, GO analysis suggested that the decreased circRNAs were enriched in cell proliferation terms, while the increased circRNAs were enriched in development and thermogenic terms. Bioinformatics predictions showed that DECs contained numerous binding sites of functional miRNAs. More interestingly, most of the circRNAs contained multiple binding sites for the same miRNA, indicating that they may facilitate functions by acting as microRNA sponges. Collectively, we characterized the circRNA expression profiles during brown adipogenesis and provide numerous novel circRNAs candidates for future brown adipogenesis regulating studies. [ABSTRACT FROM AUTHOR]

Wheat streak mosaic (WSM), a viral disease affecting cereals and grasses, causes substantial losses in crop yields. Wheat streak mosaic virus (WSMV) is the main causal agent of the complex, but mixed infections with Triticum mosaic virus (TriMV) and High plains wheat mosaic emaravirus (HPWMoV) were reported as well. Although resistant varieties are effective for the disease control, a WSMV resistance-breaking isolate and several potential resistance-breaking isolates have been reported, suggesting that viral populations are genetically diverse. Previous phylogenetic studies of WSMV were conducted by focusing only on the virus coat protein (CP) sequence, while there is no such study for either TriMV or HPWMoV. Here, we studied the genetic variation and evolutionary mechanisms of natural populations of WSM-associated viruses mainly in Kansas fields and fields in some other parts of the Great Plains using high-throughput RNA sequencing. In total, 28 historic and field samples were used for total RNA sequencing to obtain full genome sequences of WSM-associated viruses. Field survey results showed WSMV as the predominant virus followed by mixed infections of WSMV + TriMV. Phylogenetic analyses of the full genome sequences demonstrated that WSMV Kansas isolates are widely distributed in sub-clades. In contrast, phylogenetic analyses for TriMV isolates showed no significant diversity. Recombination was identified as the major evolutionary force of WSMV and TriMV variation in KS fields, and positive selection was detected in some encoding genomic regions in the genome of both viruses. Furthermore, the full genome sequence of a second Kansas HPWMoV isolate was reported. Here, we also identified previously unknown WSMV isolates in the Great Plains sharing clades and high nucleotide sequence similarities with Central Europe isolates. The findings of this study will provide more insights into the genetic structure of WSM-associated viruses and, in turn, help in improving strategies for disease management.

Wheat streak mosaic (WSM), a viral disease affecting cereals and grasses, causes substantial losses in crop yields. Wheat streak mosaic virus (WSMV) is the main causal agent of the complex, but mixed infections with Triticum mosaic virus (TriMV) and High plains wheat mosaic emaravirus (HPWMoV) were reported as well. Although resistant varieties are effective for the disease control, a WSMV resistance-breaking isolate and several potential resistance-breaking isolates have been reported, suggesting that viral populations are genetically diverse. Previous phylogenetic studies of WSMV were conducted by focusing only on the virus coat protein (CP) sequence, while there is no such study for either TriMV or HPWMoV. Here, we studied the genetic variation and evolutionary mechanisms of natural populations of WSM-associated viruses mainly in Kansas fields and fields in some other parts of the Great Plains using high-throughput RNA sequencing. In total, 28 historic and field samples were used for total RNA sequencing to obtain full genome sequences of WSM-associated viruses. Field survey results showed WSMV as the predominant virus followed by mixed infections of WSMV + TriMV. Phylogenetic analyses of the full genome sequences demonstrated that WSMV Kansas isolates are widely distributed in sub-clades. In contrast, phylogenetic analyses for TriMV isolates showed no significant diversity. Recombination was identified as the major evolutionary force of WSMV and TriMV variation in KS fields, and positive selection was detected in some encoding genomic regions in the genome of both viruses. Furthermore, the full genome sequence of a second Kansas HPWMoV isolate was reported. Here, we also identified previously unknown WSMV isolates in the Great Plains sharing clades and high nucleotide sequence similarities with Central Europe isolates. The findings of this study will provide more insights into the genetic structure of WSM-associated viruses and, in turn, help in improving strategies for disease management. [ABSTRACT FROM AUTHOR]

Homozygous familial hypercholesterolemia (HoFH) is a rare, life-threatening genetic disorder characterized by an extremely elevated serum level of low-density lipoprotein cholesterol (LDL-C) and accelerated premature atherosclerotic cardiovascular diseases (ASCVD). However, the detailed mechanism of how the pathogenic mutations of HoFH trigger the acceleration of ASCVD is not well understood. Therefore, we performed high-throughput RNA and small RNA sequencing on the peripheral blood RNA samples of six HoFH patients and three healthy controls. The gene and miRNA expression differences were analyzed, and seven miRNAs and six corresponding genes were screened out through regulatory network analysis. Validation through quantitative PCR of genes and miRNAs from 52 HoFH patients and 20 healthy controls revealed that the expression levels of hsa-miR-486-3p, hsa-miR-941, and BIRC5 were significantly upregulated in HoFH, while ID1 , PLA2G4C , and CACNA2D2 were downregulated. Spearman correlation analysis found that the levels of ID1 , hsa-miR-941, and hsa-miR-486-3p were significantly correlated with additional ASCVD risk factors in HoFH patients. This study represents the first integrated analysis of transcriptome and miRNA expression profiles in patients with HoFH, a rare disease, and as a result, six differentially expressed miRNAs/genes that may be related to atherosclerosis in HoFH are reported. The miRNA–mRNA regulatory network may be the critical regulation mechanism by which ASCVD is accelerated in HoFH. • Integrated analysis of transcriptome and miRNA expression profiles in patients with HoFH was performed. • The significantly differently expression of hsa-miR-486-3p, hsa-miR-941, and other four genes were validated in HoFH. • The levels of hsa-miR-941, and hsa-miR-486-3p were significantly correlated with additional ASCVD risk factors in HoFH. [ABSTRACT FROM AUTHOR]

Objective: There is an urgent clinical need for identifying blood-based diagnostic biomarkers for Dementia with Lewy Bodies (DLB). Transcriptomic studies have reported unique RNA changes in postmortem DLB brains. Small extracellular vesicles (SEV) that transport RNA between brain and peripheral circulation enable identifying molecular changes in living human brain. Hence, we aimed to identify differentially expressed RNA in serum SEVs from people with DLB.Methods: We investigated serum SEV total RNA profiles in people with DLB (n = 10) and age and gender matched comparisons (n = 10) using next-generation RNA-sequencing. SEVs were separated by ultracentrifugation with density gradient and were characterized by nanoparticle analysis and western blotting. We verified the differential expression levels of identified differentially expressed genes (DEG) using high-throughput qPCR. Functional implications of identified DEG were evaluated using Ingenuity pathway analyses.Results: We identified 846 nominally significant DEG including 30 miRNAs in DLB serum SEVs. We identified significant downregulation of proinflammatory genes, IL1B, CXCL8, and IKBKB. Previously reported postmortem DLB brain DEGs were significantly enriched (χ2=4.99; df=1; p = 0.03) among the identified DEGs, and the differential expression of 40 postmortem DLB brain DEGs could be detected in serum SEVs of people living with DLB. Functional pathway and network analyses highlighted the importance of immunosenescence, ubiquitin proteasome system (UPS) dysfunction, DNA repair, and RNA post-transcriptional modification deficits in DLB pathology.Conclusion: Identified DEGs, especially reduced expression levels of inflammation, and UPS-associated RNA, may aid diagnosing DLB, and their biomarker potential warrants further investigation in larger clinical cohorts. Our findings corroborate the absence of chronic neuroinflammation in DLB. [ABSTRACT FROM AUTHOR]

Background: Air pollution may be associated with increased airway responsiveness to allergens in allergic rhinitis (AR). Ozone-aged environmental black carbon (O3BC) is an important constituent of atmospheric particulate matter (PM), for which the mechanisms underlying its effects have not been fully elucidated in AR. The objective of the present study was to determine the O3BC and pollen-induced alterations in the transcriptome in human nasal epithelial cells (hNECs) in vitro. Methods: hNECs from nasal epithelial mucosal samples of healthy individuals undergoing nasal surgery (turbinoplasty or septoplasty) were established as air–liquid interface (ALI) cultures and exposed to O3BC, pollen, or a combination of O3BC+ pollen. Changes in cell viability were analyzed by fluorescence and changes in the transcriptome by high-throughput RNA sequencing (RNA-seq). Several differentially expressed genes were verified by reverse transcription-quantitative polymerase chain reaction (RT-qPCR). Enrichment analysis, based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) database, was performed to determine major biological functions and pathways involved. Results: Exposure to ≥ 50 μg/ml O3BC or 25 μg/ml O3BC+ 200 μg /ml pollen significantly decreased cell viability of the hNECs compared to control (p < 0.05) or 25 μg/ml O3BC alone (p < 0.05); whereas exposure to pollen alone did not alter cell viability at any concentration investigated. High-throughput RNA sequencing analysis indicated that there was significant difference in gene expression between pollen or O3BC alone and O3BC+ pollen exposed cells. Exposure to 200 μg/ml O3BC was associated with hypoxia stress response GO terms, whereas exposure to 25 μg/ml O3BC+ 200 μg/ml pollen was associated with inflammatory response GO terms; including regulation of neutrophil migration and chemotaxis, macrophage differentiation and chemotaxis, mast cell activation, and phagocytosis. KEGG pathway analysis indicated the top 10 upstream regulators to be IL1B, CSF1, CCL2, TLR2, LPL, IGF8, SPP1, CXCL8, FCER1G and IL1RN; of which expressions of inflammation-related genes IL1B, CSF1 and FCER1G were significantly increased. Conclusion: O3BC and pollen allergen combined exposure may induce innate immune and allergic inflammation in hNECs, and therefore potentially exacerbate the symptoms of AR in affected individuals. [ABSTRACT FROM AUTHOR]

Anorectal malformations (ARMs) are among the most common congenital terminal digestive tract malformations. Circular RNAs (circRNAs), a novel type of endogenous non-coding RNAs, play roles in the development of the digestive system; however, their contributions to the pathogenesis of ARMs are not well-established. In this study, we explored the mechanism underlying ethylenethiourea (ETU)-induced ARMs by profiling circRNA expression via RNA-seq and constructing a regulatory circRNA-miRNA-mRNA network. Nine pregnant rats were gavage-fed a single dose of 125 mg/kg 1% ETU (ARM group) on gestational day 10 (GD10), and another 9 pregnant rats received a similar dose of saline (normal group) as a control. Embryos were obtained by cesarean section on the key time-points of anorectal development (GD14, GD15, and GD16). Hindgut samples isolated from the fetuses were evaluated by high-throughput sequencing and differentially expressed circRNAs were validated by reverse transcription-quantitative polymerase chain reaction, agarose gel electrophoresis, and Sanger cloning and sequencing. A total of 18295 circRNAs were identified in the normal and ARM groups. Based on the 425 differentially expressed circRNAs (|Fc| > 2, p < 0.05), circRNA-miRNA and miRNA-mRNA pairs were predicted using miREAP, miRanda, and TargetScan. A total of 55 circRNAs (14 up- and 41 downregulated in the ARM group compared to the normal group) were predicted to bind to 195 miRNAs and 947 mRNAs. Competing endogenous RNA networks and a Kyoto Encyclopedia of Genes and Genomes analysis revealed that novel_circ_001042 had the greatest connectivity and was closely related to ARM-associated signaling pathways, such as the Wingless Type MMTV integration site family, mitogen-activated protein kinase, and transforming growth factor-β pathways. These results provide original insight into the roles of circRNAs in ARMs and provide a valuable resource for further analyses of molecular mechanisms and signaling networks.

Anorectal malformations (ARMs) are among the most common congenital terminal digestive tract malformations. Circular RNAs (circRNAs), a novel type of endogenous non-coding RNAs, play roles in the development of the digestive system; however, their contributions to the pathogenesis of ARMs are not well-established. In this study, we explored the mechanism underlying ethylenethiourea (ETU)-induced ARMs by profiling circRNA expression via RNA-seq and constructing a regulatory circRNA-miRNA-mRNA network. Nine pregnant rats were gavage-fed a single dose of 125 mg/kg 1% ETU (ARM group) on gestational day 10 (GD10), and another 9 pregnant rats received a similar dose of saline (normal group) as a control. Embryos were obtained by cesarean section on the key time-points of anorectal development (GD14, GD15, and GD16). Hindgut samples isolated from the fetuses were evaluated by high-throughput sequencing and differentially expressed circRNAs were validated by reverse transcription-quantitative polymerase chain reaction, agarose gel electrophoresis, and Sanger cloning and sequencing. A total of 18295 circRNAs were identified in the normal and ARM groups. Based on the 425 differentially expressed circRNAs (|Fc| > 2, p < 0.05), circRNA-miRNA and miRNA-mRNA pairs were predicted using miREAP, miRanda, and TargetScan. A total of 55 circRNAs (14 up- and 41 downregulated in the ARM group compared to the normal group) were predicted to bind to 195 miRNAs and 947 mRNAs. Competing endogenous RNA networks and a Kyoto Encyclopedia of Genes and Genomes analysis revealed that novel_circ_001042 had the greatest connectivity and was closely related to ARM-associated signaling pathways, such as the Wingless Type MMTV integration site family, mitogen-activated protein kinase, and transforming growth factor-β pathways. These results provide original insight into the roles of circRNAs in ARMs and provide a valuable resource for further analyses of molecular mechanisms and signaling networks. [ABSTRACT FROM AUTHOR]

Abstract Alternative polyadenylation (APA) affects most mammalian genes. The genome-wide investigation of APA has been hampered by an inability to reliably profile it using conventional RNA-seq. We describe ‘Quantification of APA’ (QAPA), a method that infers APA from conventional RNA-seq data. QAPA is faster and more sensitive than other methods. Application of QAPA reveals discrete, temporally coordinated APA programs during neurogenesis and that there is little overlap between genes regulated by alternative splicing and those by APA. Modeling of these data uncovers an APA sequence code. QAPA thus enables the discovery and characterization of programs of regulated APA using conventional RNA-seq.

Growing evidence has indicated that circular RNAs (circRNAs) play crucial roles in multiple biological processes. However, alterations in circRNA profiles during bladder cancer progression and the clinical significance thereof remain unclear. Therefore, high‐throughput RNA sequencing was conducted to identify circRNA and mRNA profiles in five pairs of bladder cancer tissues and adjacent noncancerous tissues. A total of 87 differentially expressed circRNAs and 2756 mRNAs were detected in above bladder cancer samples compared with paired noncancerous samples. Functional enrichment analyses, circRNA‐microRNA‐mRNA, and protein‐protein interaction networks revealed that these dysregulated circRNAs were potentially involved in carcinogenesis and evolution of bladder cancer. Subsequently, the differential expression of eight circRNAs was detected by real‐time qPCR. Hsa_circ_0003141 and hsa_circ_0008039 were significantly upregulated as well as hsa_circ_0026782, hsa_circ_0077837, hsa_circ_0004826, and hsa_circ_0001946 were significantly downregulated among validation of 70 matched bladder cancer tissues (≥75%). Moreover, hsa_circ_0077837 and hsa_circ_0004826 were also verified as markedly downregulated in four bladder cancer cells (100%). Naturally, hsa_circ_0077837 and hsa_circ_0004826 were also demonstrated using RNase‐R+ resistance experiments. In addition, Fisherʹs exact test, Kaplan‐Meier plots, Cox regression analyses, and receiver operating characteristic curve was performed to assess their clinical value. Downregulation of hsa_circ_0077837 and hsa_circ_0004826 all was significantly correlated with worse clinicopathological features and poor prognosis of bladder cancer patients. The area under the receiver operating characteristic curve of them was 0.775 (P <.0001) and 0.790 (P <.0001), respectively. Not surprisingly, in vitro functional experiments also demonstrated that the overexpression of hsa_circ_0077837 and hsa_circ_0004826 significantly weakened the proliferation, migration, and invasion of bladder cancer cells. Overall, hsa_circ_0077837 and hsa_circ_0004826 might act as tumor suppressors in the bladder cancer progression and serve as a potential biomarker for the diagnosis, prognosis, and therapy of bladder cancer. [ABSTRACT FROM AUTHOR]

Bone cancer pain (BCP) is one of the most common types of chronic cancer pain and its pathogenesis has not been fully understood. Long non-coding RNAs (lncRNAs) are new promising targets in the field of pain research, however, their involvements in BCP have not been reported. In the present study, we established the BCP model by implantation of Walker 256 carcinoma cells into rats' tibial medullary cavity and performed transcriptome sequencing of the ipsilateral lumbar spinal cord to explore changes in expression profiles of lncRNA and mRNA. We identified 1220 differently expressed mRNAs (DEmRNAs) (1171 up-regulated and 49 down-regulated) and 323 differently expressed lncRNAs (DElncRNAs) (246 up-regulated and 77 down-regulated) in BCP model, among which 10 DEmRNAs (5 up-regulated and 5 down-regulated) and 10 DElncRNAs (5 up-regulated and 5 down-regulated) were validated the expression by RT-qPCR. Then, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis on the expression of DEmRNAs and DElncRNAs, showing that they were mainly enriched in inflammatory and immunologic processes/pathways. Finally, we constructed a co-expression network and a ceRNA network of DEmRNAs and DElncRNAs to exhibit a potential regulatory mechanism of DElncRNAs, directly regulating protein coding gene expression in cis or in trans and indirectly regulating protein coding gene expression by sponging miRNA. In conclusion, our study provided a landscape of dysregulated lncRNA and mRNA in spinal cord of bone cancer pain and detected novel potential targets for treatment in the future. [ABSTRACT FROM AUTHOR]

Bone cancer pain (BCP)–depression comorbidity has become a complex clinical problem during cancer treatment; however, its underlying molecular mechanisms have not been clarified. Several long noncoding RNAs (lncRNAs) have been demonstrated to be promising therapeutic targets in depression, but research on the role of lncRNAs in BCP–depression comorbidity has been limited. Therefore, high-throughput RNA sequencing was performed to detect differentially expressed profiles in the amygdala of a BCP–depression rat model in this study. We detected 330 differentially expressed mRNAs (DEmRNAs) and 78 differentially expressed lncRNAs (DElncRNAs) in the BCP–depression comorbidity model and then verified the expression of six DEmRNAs and six DElncRNAs with the greatest degrees of difference by RT-qPCR. Furthermore, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that differentially expressed genes were strongly enriched in inflammatory and immunologic systemic responses. Then the nuclear factor kappa B (NF-κB) signaling pathway and the Th17 differentiation pathway showed significant differences, as determined by Western blot analysis. Finally, we constructed a protein–protein interaction (PPI) network to explore the potential regulatory mechanism of DEmRNAs. In conclusion, our study reveals a new resource for the understanding of dysregulated lncRNAs and mRNAs in BCP–depression comorbidity and provides novel potential therapeutic targets for further approaches.

Objective: Prevalence of Lewy body dementias (LBD) is second only to Alzheimer's disease (AD) among people with neurodegenerative dementia. LBD cause earlier mortality, more intense neuropsychiatric symptoms, more caregivers' burden, and higher costs than AD. The molecular mechanisms underlying LBD are largely unknown. As advancing molecular level mechanistic understanding is essential for identifying reliable peripheral biomarkers and novel therapeutic targets for LBD, the authors aimed to identify differentially expressed genes (DEG), and dysfunctional molecular networks in postmortem LBD brains.Methods: The authors investigated the transcriptomics of postmortem anterior cingulate and dorsolateral prefrontal cortices of people with pathology-verified LBD using next-generation RNA-sequencing. The authors verified the identified DEG using high-throughput quantitative polymerase chain reactions. Functional implications of identified DEG and the consequent metabolic reprogramming were evaluated by Ingenuity pathway analyses, genome-scale metabolic modeling, reporter metabolite analyses, and in silico gene silencing.Results: The authors identified and verified 12 novel DEGs (MPO, SELE, CTSG, ALPI, ABCA13, GALNT6, SST, RBM3, CSF3, SLC4A1, OXTR, and RAB44) in LBD brains with genome-wide statistical significance. The authors documented statistically significant down-regulation of several cytokine genes. Identified dysfunctional molecular networks highlighted the contributions of mitochondrial dysfunction, oxidative stress, and immunosenescence toward neurodegeneration in LBD.Conclusion: Our findings support that chronic microglial activation and neuroinflammation, well-documented in AD, are notably absent in LBD. The lack of neuroinflammation in LBD brains was corroborated by statistically significant down-regulation of several inflammatory markers. Identified DEGs, especially down-regulated inflammatory markers, may aid distinguishing LBD from AD, and their biomarker potential warrant further investigation. [ABSTRACT FROM AUTHOR]

Acute lung injury (ALI) is a serious illness that develops suddenly, progresses rapidly, has a poor treatment response and a high mortality rate. Studies have found that circular RNAs (circRNA) play a critical role in several diseases, but their role in ALI remains unclear. The aim of this study was to identify circRNAs that are associated with ALI and investigate their potential molecular mechanisms. A comparison of lung circRNA and microRNA expression profiles in mice with ALI and controls was performed by RNA-sequencing. A bioinformatic analysis was conducted to identify differentially expressed (DE) RNAs, to construct competitive endogenous RNA (ceRNA) networks, and to analyze their function and pathways. Then, a protein-protein interaction (PPI) network was generated by the Search Tool for the Retrieval of Interacting Genes database, and hub genes were identified using Cytoscape. Furthermore, a key ceRNA subnetwork was constructed based on these hub genes. Overall, we found 239 DE circRNAs and 42 DE microRNAs in ALI mice compared to controls. Additionally, the molecular mechanism of ALI was further understood by building ceRNA networks based on these DE genes. ALI-induced circRNAs are mostly function in the inflammatory response and metabolic processes. Moreover, DE circRNAs are primarily involved in the nuclear factor (NF)-kappa B and PI3K-Akt signaling pathways. Seven hub genes were derived from the PPI network of 191 genes, followed by the construction of circRNA-miRNA-hub gene subnetworks. In this study, circRNA profiles are remarkably changed in mice with LPS-triggered ALI, and their potential contribution to the disease is revealed. • Differential expression of circRNAs involved in ALI and their interactions with miRNAs and mRNAs. • Bioinformatics pipeline to study circRNA-miRNA-mRNA regulation network. • 7 hub genes (Ptpn6, Itga7, Bcl2l1, Tlr4, Vwf, Nfkb2, and Kdr). • 24 circRNA-6 miRNA- 7 hub genes interaction. [ABSTRACT FROM AUTHOR]

RNA is a polymeric molecule implicated in various biological processes, such as the coding, decoding, regulation, and expression of genes. Numerous studies have examined RNA features using whole transcriptome sequencing (RNA-seq) approaches. RNA-seq is a powerful technique for characterizing and quantifying the transcriptome and accelerates the development of bioinformatics software. In this review, we introduce routine RNA-seq workflow together with related software, focusing particularly on transcriptome reconstruction and expression quantification.

To better understand gene expression in the intestine after Shewanella algae infection and provide insights into its immune roles in the tongue sole, Cynoglossus semilaevis , sequencing-based high-throughput RNA analysis (RNA-Seq) for the intestines between the control group and 12 h post-injection group was performed. After assembly, there was an average of 23,957,159 raw sequencing reads, and 23,943,491 clean reads were obtained after filtering out low-quality reads. Then, 383 differentially expressed genes (DEGs) in the intestines in response to S. algae infection were identified. Subsequently, gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of the DEGs were conducted to further explore their functions. Among all of the pathways involved, sixteen pathways were related to the immune system, among which the complement and coagulation cascades pathway was the most prominent for immunity-related DEGs, followed by the leukocyte transendothelial migration pathway. Furthermore, the expression levels of twelve selected DEGs in the immune-related pathways were identified by quantitative real-time polymerase chain reaction, substantiating the reliability and reproducibility of the RNA-Seq results. In summary, this study represents an important genomic resource for understanding the potential immune role of the tongue sole intestine from the perspective of gene expression. [ABSTRACT FROM AUTHOR]

Downy mildew caused by Hyaloperonospora brassicae is one of the most damaging diseases of Chinese cabbage. In this study, we analyzed genome-wide gene expression profiles of a resistant line 318R and a susceptible line 361S of Chinese cabbage in response to downy mildew infection by high throughput RNA sequencing technology. More transcripts were identified in H. brassicae infected samples in both lines, especially in the susceptible line 361S. A total of 2297 and 3338 differentially expressed genes were identified in 318R and 361S, respectively, which are mainly involved in the plant-pathogen interaction pathway, plant hormone signal transduction pathway, and biosynthesis of secondary metabolites pathway. A total of 1284 genes were found to be differentially expressed in both lines, among which, 68 genes showed higher variation in the resistant line than in the susceptible line. These 68 genes, along with the 1012 genes that were only differentially expressed in the resistant line, were annotated to defense response genes, transcriptional factors, and cell wall related genes, which function in the plant-pathogen interaction pathway, protein processing in endoplasmic reticulum pathway, or metabolism pathways. These genes are important candidates for understanding the molecular mechanisms of host response to downy mildew. [ABSTRACT FROM AUTHOR]

Myoblast differentiation is a highly complex process that is regulated by proteins as well as by non-coding RNAs. Circular RNAs have been identified as an emerging new class of non-coding RNA in the modulation of skeletal muscle development, whereas their expression profiles and functional regulation in myoblast differentiation remain unknown. In the present study, we performed deep RNA-sequencing of C2C12 myoblasts during cell differentiation and uncovered 37,751 unique circular RNAs derived from 6943 hosting genes. The ensuing qRT-PCR and RNA fluorescence in situ hybridization verification were carried out to confirm the RNA-sequencing results. An unbiased analysis demonstrated dynamic circular RNA expression changes in the process of myoblast differentiation, and the circular RNA abundances were independent from their cognate linear RNAs. Gene ontology analysis showed that many down-regulated circular RNAs were exclusive to cell division and the cell cycle, whereas up-regulated circular RNAs were related to the cell development process. Furthermore, interaction networks of circular RNA-microRNA were constructed. Several microRNAs well-known for myoblast regulation, such as miR-133, miR-24 and miR-23a, were in this network. In summary, this study showed that circular RNA expression dynamics changed during myoblast differentiation. Circular RNAs play a role in regulating the myoblast cell cycle and development by acting as microRNA binding sites to facilitate their regulation of gene expression during myoblast differentiation. These findings open a new avenue for future investigation of this emerging RNA class in skeletal muscle growth and development. [ABSTRACT FROM AUTHOR]

Leishmania (L.) tropica is the main causative agent of anthroponotic cutaneous leishmaniasis (CL) in Iran. Defining the host inflammatory response in the L. tropica lesions are crucial for the development of new treatment modalities. High-throughput RNA sequencing provides a powerful method for characterization of the human gene expression profile in L. tropica lesions. Comparing the transcription profile of the L. tropica skin lesions with normal skin identified over 5000 differentially regulated genes. Gene set enrichment analysis indicated significant activation of key immunological pathways related to antigen processing and presentation. In addition, we observed a substantial upregulation of immunoglobulin genes in lesion samples, highlighting the remarkable involvement of B cells in the infection site. To our knowledge, this study is the first report to build a comprehensive picture of transcriptome changes in acute human skin lesions during infection by L. tropica . [ABSTRACT FROM AUTHOR]

MicroRNAs (miRNAs) exhibit aberrant expression in the initiation and progression of a variety of human cancers, including colorectal cancer (CRC). However, the exact mechanisms are not well defined. miRNA expression profiles were characterized by microarrays in CRC samples, and miRNA 18b (miR-18b) was increased significantly in tumor tissues. The expression of miR-18b was confirmed in the CRC cell lines SW480 and HCT116 and 44 clinical specimens by quantitative real-time PCR (qRT-PCR). Multiple linear regression analysis showed a strong correlation of miR-18b expression with lymph node and distant metastasis. Overexpression of miR-18b promoted cell proliferation by facilitating cell cycle progression, and knockdown of miR-18b significantly suppressed migration in CRC cells. CDKN2B was identified as a target of miR-18b by highthroughput RNA sequencing and bioinformatics. After transfection with a miR-18b mimic, expression of CDKN2B was reduced significantly in CRC cells, and the effect was restored when a miR-18b inhibitor was transfected. A luciferase assay indicated miR-18b directly binds to the 3= untranslated region (UTR) of CDKN2B. Expression of CDKN2B was downregulated in patient cancer tissues and negatively correlated with miR-18b. In a model of ectopic expression of miR-18b and CDKN2B, CDKN2B overexpression antagonized the effects of miR-18b in vitro and in vivo. The data show that miR-18b is involved in CRC carcinogenesis through targeting CDKN2B. [ABSTRACT FROM AUTHOR]