The role of exosomal circular RNAs (circRNAs) in Hepatocellular carcinoma (HCC) cells with high metastatic potential has been little studied.

Alternative splicing (AS) is an important biological process for regulating the expression of various isoforms from a single gene and thus to promote proteome diversity. In this study, RNA-seq data from 15 pairs of matched esophageal squamous cell carcinoma (ESCC) and normal tissue samples as well as two cell lines were analyzed. AS events with significant differences were identified between ESCC and matched normal tissues, which were re-annotated to find protein coding genes or non-coding RNAs. A total of 45,439 AS events were found. Of these, 6019 (13.25%) significant differentially AS events were identified. Exon skipping (SE) events occupied the largest proportion of abnormal splicing events. Fifteen differential splicing events with the same trends of ΔΨ values in ESCC tissues, as well in the two cell lines were found. Four pathways and 20 biological processes related to pro-metastasis cell junction and migration were significantly enriched for the differentially spliced genes. The upregulated splicing factor SF3B4, which regulates 92 gene splicing events, could be a potential prognostic factor of ESCC. Differentially spliced genes, including HNRNPC, VCL, ZNF207, KIAA1217, TPM1 and CALD1 are shown with a sashimi plot. These results suggest that cell junction- and migration-related biological processes are influenced by AS abnormalities, and aberrant splicing events can be affected by splicing factor expression changes. The involved splicing factor SF3B4 was found to be a survival-related gene in ESCC and is presumed to regulate AS in multiple cancers. In summary, we identified significant differentially expressed AS events which may be related to the development of ESCC.

Accurately reconstructing gene co-expression network is of great importance for uncovering the genetic architecture underlying complex and various phenotypes. The recent availability of high-throughput RNA-seq sequencing has made genome-wide detecting and quantifying of the novel, rare and low-abundance transcripts practical. However, its potential merits in reconstructing gene co-expression network have still not been well explored. Using massive-scale RNA-seq samples, we have designed an ensemble pipeline, called NetMiner, for building genome-scale and high-quality Gene Co-expression Network (GCN) by integrating three frequently used inference algorithms. We constructed a RNA-seq-based GCN in one species of monocot rice. The quality of network obtained by our method was verified and evaluated by the curated gene functional association data sets, which obviously outperformed each single method. In addition, the powerful capability of network for associating genes with functions and agronomic traits was shown by enrichment analysis and case studies. In particular, we demonstrated the potential value of our proposed method to predict the biological roles of unknown protein-coding genes, long non-coding RNA (lncRNA) genes and circular RNA (circRNA) genes. Our results provided a valuable and highly reliable data source to select key candidate genes for subsequent experimental validation. To facilitate identification of novel genes regulating important biological processes and phenotypes in other plants or animals, we have published the source code of NetMiner, making it freely available at https://github.com/czllab/NetMiner.

Combination targeted therapy is commonly used to treat acute myeloid leukemia (AML) patients, particularly in refractory/relapse (RR) population. However, concerns have been raised regarding the safety and patient tolerance of combination chemotherapy. It is critical to choose the appropriate treatment for precision therapy. We performed genome-wide RNA profiling using RNA-Seq to compare the RR group and the complete remission (CR) group (a total of 42 adult AML patients). The Hedgehog (Hh) and PI3K/AKT pathways were upregulated in the RR population, which was further confirmed by western blot and/or qPCR. Overexpression of GLI1 in AML cells led to increased AKT phosphorylation and decreased drug sensitivity, which was attenuated by GLI1 inhibition. By contrast, neither the expression of GLI1 nor apoptosis in response to Ara-C treatment of AML cells was significantly affected by PI3K inhibition. Furthermore, co-inhibition of GLI1 and PI3K induced apoptosis of hematopoietic stem/progenitor cells (HSPCs), which raised serious concerns about the side effects of this treatment. These results indicated that GLI1 inhibition alone, but not combined inhibition, is sufficient to enhance AML drug sensitivity, which provides a novel therapeutic strategy for AML treatment.

Candida haemulonii, a close relative of Candida auris, is an emerging pathogen which frequently shows multidrug resistance especially to triazoles, the most used antifungal drugs. The mechanisms of drug resistance in C. haemulonii, however, are largely unknown. Here, we sequenced and annotated the genomes of two reference strains from the C. haemulonii complex, compared the phenotypes, genomes, and transcriptomes of a triazole-susceptible and two triazole-resistant C. haemulonii strains, and identified triazole susceptibility, morphology, fitness, and the major genetic and gene expression differences between the strains. A multidrug efflux gene, CDR1, was recurrently found to be upregulated for expression in triazole-resistant strains. Blocking the activity of Cdr1 increased the susceptibility to triazoles strikingly. Comparative transcriptome analysis also demonstrated impaired cell wall integrity, filamentous growth, and iron homeostasis in triazole-resistant strains. Finally, we also identified a zinc-binding MHR family transcription regulator gene that mutated in triazole-resistant strains spontaneously, contributing to the changes of morphology and, possibly, cell wall integrity between the strains. The study provided important clues for future studies exploring the mechanisms of multidrug resistance and related phenotypic differences of C. haemulonii strains.IMPORTANCE A comprehensive, multi-omic survey was performed to disclose the genetic backgrounds and differences between multidrug-resistant and -susceptible C. haemulonii strains. Genes were identified with mutations or significant expression differences in multidrug-resistant compared to multidrug-susceptible strains, which were mainly involved in multidrug resistance, stress fitness, and morphology. The Cdr1-encoding gene, C. haemulonii 2486 (CH_2486), was expressed at a significantly increased level in multidrug-resistant strains. Functional inhibition experiments further implicated potential roles of CH_2486 in drug resistance. A gene spontaneously mutated in resistant strains, CH_4347, was experimentally validated to influence the morphology of spores, possibly by controlling cell wall integrity.

Differentiated somatic cells can be reprogrammed into totipotent embryos through somatic cell nuclear transfer. However, most cloned embryos arrest at early stages and the underlying molecular mechanism remains largely unexplored. Here, we first developed a somatic cell nuclear transfer embryo biopsy system at two- or four-cell stage, which allows us to trace the developmental fate of the biopsied embryos precisely. Then, through single-cell transcriptome sequencing of somatic cell nuclear transfer embryos with different developmental fates, we identified that inactivation of Kdm4b, a histone H3 lysine 9 trimethylation demethylase, functions as a barrier for two-cell arrest of cloned embryos. Moreover, we discovered that inactivation of another histone demethylase Kdm5b accounts for the arrest of cloned embryos at the four-cell stage through single-cell analysis. Co-injection of Kdm4b and Kdm5b can restore transcriptional profiles of somatic cell nuclear transfer embryos and greatly improve the blastocyst development (over 95%) as well as the production of cloned mice. Our study therefore provides an effective approach to identify key factors responsible for the developmental arrest of somatic cell cloned embryos.

Glucocorticoids are important therapeutic compounds for acute lymphoblastic leukemia (ALL). AKT1 or the protein kinase B is frequently activated in ALL, and contributes to the development of glucocorticoid resistance. We examined impact of AKT1 on glucocorticoid receptor (GR)-induced transcriptional activity in cooperation with phospho-serine/threonine-binding protein 14-3-3. AKT1 has two distinct actions on GR transcriptional activity, one through segregation of GR in the cytoplasm by phosphorylating GR at Ser-134 and subsequent association of 14-3-3, and the other through direct modulation of GR transcriptional activity in the nucleus. For the latter, AKT1 and 14-3-3 are attracted to DNA-bound GR, accompanied by AKT1-dependent p300 phosphorylation, H3S10 phosphorylation and H3K14 acetylation at the DNA site. These two actions of AKT1 regulate distinct sets of glucocorticoid-responsive genes. Our results suggest that specific inhibition of the AKT1/14-3-3 activity on the cytoplasmic retention of GR may be a promising target for treating glucocorticoid resistance observed in ALL.

The accurate generation of an appropriate number of different neuronal and glial subtypes is fundamental to normal brain functions and requires tightly orchestrated spatial and temporal developmental programmes to maintain the balance between the proliferation and the differentiation of neural progenitor cells. However, the molecular mechanism governing this process has not been fully elucidated. Here, we found that miR-214-3p was highly expressed in neural progenitor cells and dynamically regulated during neocortical development. Moreover, our in vivo and in vitro studies showed that miR-214 inhibited self-renewal of neural progenitor cells and promoted neurogenesis. In addition, after target screening, we identified miR-214 targets including Quaking (Qki) by binding the 3'- untranslated region (3'-UTR) of the Qki mRNA, which was specifically expressed in the progenitor cells of the proliferative ventricular zone as 3 Qki isoforms. Furthermore, overexpression and knockdown of Qki showed that the different isoforms of Qki had different functions in the regulation of neural progenitor cells differentiation. Moreover, overexpression of Qki could counteract the function of miR-214 in neurogenesis. Our results revealed that miR-214 maintains the balance between neural progenitor/stem cell proliferation and differentiation together with Quaking, its target gene.

D-2-hydroxyglutaric aciduria (D2HGA) type II is a rare neurometabolic disorder caused by germline gain-of-function mutations in isocitrate dehydrogenase 2 (IDH2), resulting in accumulation of D-2-hydroxyglutarate (D2HG). Patients exhibit a wide spectrum of symptoms including cardiomyopathy, epilepsy, developmental delay and limited life span. Currently, there are no effective therapeutic interventions. We generated a D2HGA type II mouse model by introducing the Idh2R140Q mutation at the native chromosomal locus. Idh2R140Q mice displayed significantly elevated 2HG levels and recapitulated multiple defects seen in patients. AGI-026, a potent, selective inhibitor of the human IDH2R140Q-mutant enzyme, suppressed 2HG production, rescued cardiomyopathy, and provided a survival benefit in Idh2R140Q mice; treatment withdrawal resulted in deterioration of cardiac function. We observed differential expression of multiple genes and metabolites that are associated with cardiomyopathy, which were largely reversed by AGI-026. These findings demonstrate the potential therapeutic benefit of an IDH2R140Q inhibitor in patients with D2HGA type II.

The heterogeneous nature of acute myeloid leukemia (AML) and its poor prognosis necessitate therapeutic improvement. Current advances in AML research yield important insights regarding both AML genetics and epigenetics. MicroRNAs (miRNAs) play important roles in cell proliferation, differentiation, and survival and may be useful for AML diagnosis and prognosis. In this study, a novel miRNA, hsa-miR-12462, was identified in bone marrow (BM) samples from AML patients at diagnosis by small RNA sequencing. A significant higher level of hsa-miR-12462 was found in patients who achieve complete remission (AML-CR) after induction therapy compared with those who suffer relapse/refractory (AML-RR). FosB was predicted to be the target of hsa-miR-12462 through RNA sequencing, bioinformatics analysis, and protein-protein interaction (PPI) network analysis and then verified by luciferase activity assay. T-5224, the inhibitor of FosB, was administered to AML cell lines, which could inhibit cell proliferation, promote apoptosis, and restore the sensitivity of AML cells to cytarabine (Ara-C). In summary, a higher level of hsa-miR-12462 in AML cells is associated with increased sensitivity to Ara-C via targeting FosB.

As a class of transcription regulators, numerous miRNAs have been verified to participate in regulating ovary follicular development in chickens (Gallus gallus). Previously we showed that gga-miR-135a-5p has significant differential expression between high and low-yield chicken ovaries, and the abundance of gga-miR-135a-5p is significantly higher in follicular theca cells than in granulosa cells. However, the exact role of gga-miR-135a-5p in chicken follicular theca cells is unclear. In this study, primary chicken follicular theca cells were isolated and then transfected with gga-miR-135a-5p inhibitor. Transcriptome sequencing was performed in chicken follicular theca cells with or without transfection. Differentially expressed genes (DEGs) were analyzed using bioinformatics. A dual-luciferase reporter assay was used to verify the target relationship between gga-miR-135a-5p and predicted targets within the DEGs. Compared with the normal chicken follicle theca cells, 953 up-regulated and 1060 down-regulated genes were detected in cells with gga-miR-135a-5p inhibited. The up-regulated genes were significantly enriched in Gene Ontology terms and pathways involved in cell proliferation and differentiation. In chicken follicular theca cells, Krüppel-like factor 4 (KLF4), ATPase phospholipid transporting 8A1 (ATP8A1), and Complexin-1 (CPLX1) were significantly up-regulated when the expression of gga-miR-135a-5p was inhibited. In addition, KLF4, ATP8A1, and CPLX1 confirmed as targets of gga-miR-135a-5p by using a dual-luciferase assay in vitro The results suggest that gga-mir-135a-5p may involve in proliferation and differentiation in chicken ovarian follicular theca cells by targeting KLF4, ATP8A1, and CPLX1.

Expression variation plays an important role in plant adaptation, but little is known about the factors impacting the expression variation when population adapts to changing environment. We used RNA-seq data from 80 individuals in 14 Miscanthus lutarioriparius populations, which were transplanted into a harsh environment from native habitat, to investigate the expression level, expression diversity and genetic diversity for genes expressed in both environments. The expression level of genes with lower expression level or without SNP tended to be more changeable in new environment, which suggested highly expressed genes experienced stronger purifying selection than those at lower level. Low proportion of genes with population effect confirmed the weak population structure and frequent gene flow in these populations. Meanwhile, the number of genes with environment effect was the most frequent compared with that with population effect. Our results showed that environment and genetic diversity were the main factors determining gene expression variation in population. This study could facilitate understanding the mechanisms of global gene expression variation when plant population adapts to changing environment.

Seawater aspiration‑induced acute lung injury (ALI) is a syndrome associated with a high mortality rate, which is characterized by severe hypoxemia, pulmonary edema and inflammation. The present study is the first, to the best of our knowledge, to analyze gene expression profiles from a rat model of seawater aspiration‑induced ALI. Adult male Sprague‑Dawley rats were instilled with seawater (4 ml/kg) in the seawater aspiration‑induced ALI group (S group) or with distilled water (4 ml/kg) in the distilled water negative control group (D group). In the blank control group (C group) the rats' tracheae were exposed without instillation. Subsequently, lung samples were examined by histopathology; total protein concentration was detected in bronchoalveolar lavage fluid (BALF); lung wet/dry weight ratios were determined; and transcript expression was detected by gene sequencing analysis. The results demonstrated that histopathological alterations, pulmonary edema and total protein concentrations in BALF were increased in the S group compared with in the D group. Analysis of differential gene expression identified up and downregulated genes in the S group compared with in the D and C groups. A gene ontology analysis of the differential gene expression revealed enrichment of genes in the functional pathways associated with neutrophil chemotaxis, immune and defense responses, and cytokine activity. Kyoto Encyclopedia of Genes and Genomes analysis revealed that the cytokine‑cytokine receptor interaction pathway was one of the most important pathways involved in seawater aspiration‑induced ALI. In conclusion, activation of the cytokine‑cytokine receptor interaction pathway may have an essential role in the progression of seawater aspiration‑induced ALI, and the downregulation of tumor necrosis factor superfamily member 10 may enhance inflammation. Furthermore, IL‑6 may be considered a biomarker in seawater aspiration‑induced ALI.

Fiber quality is an important economic trait of cotton, and its improvement is a major goal of cotton breeding. To better understand the genetic mechanisms responsible for fiber quality traits, we conducted a genome-wide association study to identify and mine fiber-quality-related quantitative trait loci (QTLs) and genes.

Overexposure to manganese (Mn) is neurotoxic. Our previous research has demonstrated that the interaction of endoplasmic reticulum (ER) stress and autophagy participates in the early stage of Mn-mediated neurotoxicity in mouse. However, the mechanisms of ER stress signalling pathways in the initiation of autophagy remain confused. In the current study, we first validated that ER stress-mediated cell apoptosis is accompanied by autophagy in SH-SY5Y cells. Then, we found that inhibiting ER stress with 4-phenylbutyrate (4-PBA) decreased ER stress-related protein expression and reduced cell apoptosis, whereas blocking autophagy with 3-methyladenine (3-MA) increased cell apoptosis. These data indicate that protective autophagy was activated to alleviate ER stress-mediated apoptosis. Knockdown of the protein kinase RNA-like ER kinase (PERK) gene inhibited Mn-induced autophagy and weakened the interaction between ATF4 and the LC3 promoter. Our results reveal a novel molecular mechanism in which ER stress may regulate autophagy via the PERK/eIF2α/ATF4 signalling pathway. Additionally, Mn may activate protective autophagy to alleviate ER stress-mediated apoptosis via the PERK/eIF2α/ATF4 signalling pathway in SH-SY5Y cells.

Understanding the genetic basis of water use efficiency (WUE) and its roles in plant adaptation to a drought environment is essential for the production of second-generation energy crops in water-deficit marginal land. In this study, RNA-Seq and WUE measurements were performed for 78 individuals of Miscanthus lutarioriparius grown in two common gardens, one located in warm and wet Central China near the native habitats of the species and the other located in the semiarid Loess Plateau, the domestication site of the energy crop. The field measurements showed that WUE of M. lutarioriparius in the semiarid location was significantly higher than that in the wet location. A matrix correlation analysis was conducted between gene expression levels and WUE to identify candidate genes involved in the improvement of WUE from the native to the domestication site. A total of 48 candidate genes were identified and assigned to functional categories, including photosynthesis, stomatal regulation, protein metabolism, and abiotic stress responses. Of these genes, nearly 73% were up-regulated in the semiarid site. It was also found that the relatively high expression variation of the WUE-related genes was affected to a larger extent by environment than by genetic variation. The study demonstrates that transcriptome-wide correlation between physiological phenotypes and expression levels offers an effective means for identifying candidate genes involved in the adaptation to environmental changes.

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

Long non-coding RNA (lncRNA) is a class of non-coding RNA with important regulatory roles in biological process of organisms. The systematic comparison of lncRNAs with protein coding mRNAs in population expression and their response to environmental change are still poorly understood. Here we identified 17,610 lncRNAs and calculated their expression levels based on RNA-seq of 80 individuals of Miscanthus lutarioriparius from two environments, the nearly native habitats and transplanted field, respectively.

Peach (Prunus persica L.) generally exhibits self-pollination, however, they can also be pollinated by other varieties of pollen. Here we found two varieties that are different from other peaches: 'Daifei' and 'Liuyefeitao'. 'Daifei' produces less pollen, which needs artificial pollination, honeybee pollination, and the fruit setting depends on other varieties of peach pollen. 'Liuyefeitao' exhibits strictly self-pollination, hence pollen from other species is rejected. To explore the mechanism of this phenomenon, we performed a high-throughput sequencing of the stigma (including style) of 'Daifei' and 'Liuyefeitao' to explain the rejection mechanism of other varieties of pollen of 'Liuyefeitao' peach. In our study, we found one S gene, and lots of non-S-locus factors such as: F-box proteins, Ub/26S, MAPKs, RLK, and transcription factor were differential expressed between 'Daifei' and 'Liuyefeitao'. We supposed that the strictly self-compatible of 'Liuyefeitao' may result from the synthesis of these factors.

PCIF1 (phosphorylated CTD interacting factor 1) is the first reported RNA N6,2'-O-dimethyladenosine (m6Am) methyltransferase. However, the pathological significance of PCIF1 and m6Am modification remains unknown. Here we find that both PCIF1 expression and m6Am modification are significantly elevated in gastric cancer tissues. Increased PCIF1 is associated with gastric cancer progression, and predicts poor prognosis. Silence of PCIF1 inhibits the proliferation and invasion of gastric cancer cells, and suppresses tumor growth and metastasis in mouse model. m6Am-seq analysis reveals TM9SF1 (transmembrane 9 superfamily member 1) as a target of PCIF1. PCIF1 modifies TM9SF1 mRNA with m6Am leading to decreased TM9SF1 translation. TM9SF1 reverses the effects of PCIF1 on gastric cancer cell aggressiveness. Collectively, our work uncovers an oncogenic function of PCIF1, providing insights into the critical role of m6Am modification in cancer progression.