PAIP1 regulates expression of immune and inflammatory response associated genes at transcript level in liver cancer cell.

Poly(A) binding protein interacting protein 1 (PAIP1) is a translation regulator and also regulate the decay of mRNA. PAIP1 has also been reported to be a marker of increased invasive potential of liver cancer. However, the roles and underlying molecular mechanism of PAIP1 in liver cancer is still unclear. Here, cell viability and the gene expression profile of liver cancer line HepG2 transfected with PAIP1 siRNA was compared with cells transfected with non-targeting control siRNA. The results showed that PAIP1 knockdown inhibited cell viability, and extensively affects expression of 893 genes at transcriptional level in HepG2 cells. Gene function analysis showed that a large number of PAIP1 up-regulated genes were enriched in term of DNA-dependent transcription and the down-regulated genes were enriched in some pathways including immune response and inflammatory response. qPCR confirmed that PAIP1 knockdown positively regulated the expression of selected immune and inflammatory factor genes in HepG2 cells. Expression analysis of TCGA revealed that PAIP1 had positive correlations with two immune associated genes IL1R2 and PTAFR in liver tumor tissue. Taken together, our results demonstrated that PAIP1 was not only a translation regulator, but also a transcription regulator in liver cancer. Moreover, PAIP1 could function as a regulatory factor of immune and inflammatory genes in liver cancer. Thus, our study provides important cues for further study on the regulatory mechanism of PAIP1 in liver cancer.

Attenuated Duck Hepatitis A Virus Infection Is Associated With High mRNA Maintenance in Duckling Liver via m6A Modification.

Host translation is generally modulated by viral infection, including duck hepatitis A virus (DHAV) infection. Previously, we reported that cellular protein synthesis in a cell model of duck embryo fibroblasts is significantly inhibited by DHAV infection but not viral proteins, suggesting that an important viral-host interaction occurs at the translational level. In this study, we aim to further understand the impact of DHAV virulence on cellular N6-methyladenosine (m6A) modification, which is essential to a wide variety of RNA biological processes, such as mRNA stabilization and translation. Using m6A antibody-based immunoprecipitation, m6A-seq, and LC-MS/MS, we observed that m6A-modified mRNA exists in both virulent and attenuated DHAV-infected duckling livers. Importantly, m6A levels in mRNA were much higher in attenuated DHAV-infected livers compared with virulent DHAV-infected livers, suggesting virulence-dependent regulation of m6A modification. Analysis of modification motifs indicated that GAAGAAG is the most enriched motif. Combined m6A-seq and RNA-seq data analysis indicated a generally positive correlation between m6A and mRNA expression levels in DHAV-infected duckling livers. GO analysis of genes with decreased or increased m6A levels showed that these genes were enriched in various terms, including oxidation-reduction processes and antiviral immune responses. Collectively, our work reveals DHAV virulence-dependent coordination between m6A modification and mRNA expression in duckling livers.

Transcriptome analysis reveals the potential biological function of FSCN1 in HeLa cervical cancer cells.

Fascin actin-bundling protein 1 (FSCN1), an actin-bundling protein associated with cell migration and invasion, is highly expressed in various tumor tissues. FSCN1 has also been reported to be a marker of increased invasive potential in cervical cancers. However, the functions of FSCN1 are still not fully understood in cervical cancers. Here, the gene expression profile of HeLa cells transfected with FSCN1 shRNA (shFSCN1) was compared with that of cells transfected with empty vector (shCtrl). The results showed that shFSCN1 extensively affected the transcription level of 5,043 genes in HeLa cells. In particular, Gene Ontology (GO) analysis showed that a large number of upregulated genes were annotated with terms including transcription regulation and DNA binding. The downregulated genes were enriched in some cancer pathways, including angiogenesis and cell adhesion. qPCR validation confirmed that FSCN1 knockdown significantly affected the expression of selected genes in HeLa cells either negatively or positively. Expression analysis in TCGA (The Cancer Genome Atlas) revealed that FSCN1 had negative correlations with several transcription factors and a positive correlation with an angiogenic factor (angiopoietin like 4, ANGPTL4) in cervical tumor tissue. In particular, validation by Western blotting showed that FSCN1 knockdown decreased the protein level of ANGPTL4. Our results demonstrated that FSCN1 is not only an actin-binding protein but also a transcriptional regulator and an angiogenic factor in cervical cancer. Thus, our study provides important insights for further study on the regulatory mechanism of FSCN1 in cervical cancer.

Analysis of Bulk RNA Sequencing Data Reveals Novel Transcription Factors Associated With Immune Infiltration Among Multiple Cancers.

Tumor-infiltrating immune cells shape the tumor microenvironment and are closely related to clinical outcomes. Several transcription factors (TFs) have also been reported to regulate the antitumor activity and immune cell infiltration. This study aimed to quantify the populations of different immune cells infiltrated in tumor samples based on the bulk RNA sequencing data obtained from 50 cancer patients using the CIBERSORT and the EPIC algorithm. Weighted gene coexpression network analysis (WGCNA) identified eigengene modules strongly associated with tumorigenesis and the activation of CD4+ memory T cells, dendritic cells, and macrophages. TF genes FOXM1, MYBL2, TAL1, and ERG are central in the subnetworks of the eigengene modules associated with immune-related genes. The analysis of The Cancer Genome Atlas (TCGA) cancer data confirmed these findings and further showed that the expression of these potential TF genes regulating immune infiltration, and the immune-related genes that they regulated, was associated with the survival of patients within multiple cancers. Exome-seq was performed on 24 paired samples that also had RNA-seq data. The expression quantitative trait loci (eQTL) analysis showed that mutations were significantly more frequent in the regions flanking the TF genes compared with those of non-TF genes, suggesting a driver role of these TF genes regulating immune infiltration. Taken together, this study presented a practical method for identifying genes that regulate immune infiltration. These genes could be potential biomarkers for cancer prognosis and possible therapeutic targets.

SUVA: splicing site usage variation analysis from RNA-seq data reveals highly conserved complex splicing biomarkers in liver cancer.

Most of the current alternative splicing (AS) analysis tools are powerless to analyse complex splicing. To address this, we developed SUVA (Splice sites Usage Variation Analysis) that decomposes complex splicing events into five types of splice junction pairs. By analysing real and simulated data, SUVA showed higher sensitivity and accuracy in detecting AS events than the compared methods. Notably, SUVA detected extensive complex AS events and screened out 69 highly conserved and dominant AS events associated with cancer. The cancer-associated complex AS events in FN1 and the co-regulated RNA-binding proteins were significantly correlated with patient survival.

An integrated analysis of lncRNA and mRNA expression profiles in the kidneys of mice with lupus nephritis.

Long noncoding RNAs (lncRNAs) are persistently expressed and have been described as potential biomarkers and therapeutic targets in various diseases. However, there is limited information regarding lncRNA expression in the tissue of kidney exhibiting lupus nephritis (LN)a serious complication of systemic lupus erythematosus (SLE). In this study, RNA sequencing (RNA-seq) was performed to characterize the lncRNA and mRNA expression in kidney tissues from LN (MRL/lpr) and control mice. We identified 12,979 novel lncRNAs in mouse. The expression profiles of both mRNAs and lncRNAs were differed significantly between LN and control mice. In particular, there were more upregulated lncRNAs and mRNAs than downregulated ones in the kidney tissues of LN mice. However, GO analysis showed that more downregulated genes were enriched in immune and inflammatory response-associated pathways. KEGG analysis showed that both downregulated and upregulated genes were enriched in a number of pathways, including the SLE pathway, and approximately half of these SLE-associated genes encoded inflammatory factors. Moreover, we observed that 2,181 DElncRNAs may have targeted and regulated the expression of 778 mRNAs in LN kidney tissues. The results of this study showed that 11 DElncRNAs targeted and were co-expressed with six immune and SLE-associated genes. qPCR analysis confirmed that lncRNA Gm20513 positively regulated the expression of the SLE-associated gene H2-Aa. In conclusion, the results of our study demonstrates that lncRNAs influence the progression of LN and provide some cues for further study of lncRNAs in LN. These results regarding the lncRNA-mRNAregulatory network may have important value in LN diagnosis and therapy.

Effect of emodin on long non-coding RNA-mRNA networks in rats with severe acute pancreatitis-induced acute lung injury.

Long non-coding RNAs (lncRNAs) contribute to disease pathogenesis and drug treatment effects. Both emodin and dexamethasone (DEX) have been used for treating severe acute pancreatitis-associated acute lung injury (SAP-ALI). However, lncRNA regulation networks related to SAP-ALI pathogenesis and drug treatment are unreported. In this study, lncRNAs and mRNAs in the lung tissue of SAP-ALI and control rats, with or without drug treatment (emodin or DEX), were assessed by RNA sequencing. Results showed both emodin and DEX were therapeutic for SAP-ALI and that mRNA and lncRNA levels differed between untreated and treated SAP-ALI rats. Gene expression profile relationships for emodin-treated and control rats were higher than DEX-treated and -untreated animals. By comparison of control and SAP-ALI animals, more up-regulated than down-regulated mRNAs and lncRNAs were observed with emodin treatment. For DEX treatment, more down-regulated than up-regulated mRNAs and lncRNAs were observed. Functional analysis demonstrated both up-regulated mRNA and co-expressed genes with up-regulated lncRNAs were enriched in inflammatory and immune response pathways. Further, emodin-associated lncRNAs and mRNAs co-expressed modules were different from those associated with DEX. Quantitative polymerase chain reaction demonstrates selected lncRNA and mRNA co-expressed modules were different in the lung tissue of emodin- and DEX-treated rats. Also, emodin had different effects compared with DEX on co-expression network of lncRNAs Rn60_7_1164.1 and AABR07062477.2 for the blue lncRNA module and Nrp1 for the green mRNA module. In conclusion, this study provides evidence that emodin may be a suitable alternative or complementary medicine for treating SAP-ALI.

Comparative transcriptomics of the pheromone glands provides new insights into the differentiation of sex pheromone between two host populations of Chilo suppressalis.

Reproductive isolation between different host populations is often based on intraspecific sex pheromone differences. The mechanisms underlying these differences have not been thoroughly elucidated to date. Previous studies suggested that Chilo suppressalis has differentiated into rice and water-oat host populations, and these two populations manifest clear differences in sex pheromone titer and mating rhythm. Hence, this moth is an ideal model to investigate the endogenous mechanisms of intraspecific reproductive isolation. Here, we identified a series of putative genes associated with sex pheromone biosynthesis based on the C. suppressalis pheromone gland transcriptome data. Transcripts of most genes were at higher level in the rice population. Then we obtained 11 pivotal differentially expressed genes (DEGs). The expression levels of these DEGs exhibited a distinct increase in the rice population. Moreover, we also observed the expression rhythm of these DEGs is discrepant between two host populations. Our study offers a new understanding to elucidate the mechanisms of intraspecific reproductive isolation.

RBM4 modulates the proliferation and expression of inflammatory factors via the alternative splicing of regulatory factors in HeLa cells.

Regulatory factors function by modulating a variety of cascade mechanisms in cells. RBM4 is a multifunctional RNA-binding protein in post-transcriptional gene regulation. Cytoplasmic RBM4 interacts with Ago2 to regulate inflammatory responses by affecting mRNA decay and cap-dependent translation. However, it is unclear whether RBM4 functions in inflammation regulation by its splicing factor role. Here, the cell biology, gene expression profile and alternative splicing pattern of HeLa cells with RBM4 overexpression (RBM-OE) were compared with the control. The results showed that RBM4-OE inhibited proliferation. RBM4-OE extensively affects the transcriptional level of genes involved in cell surface receptor signalling pathway, inflammatory responses and the response to lipopolysaccharide. RBM4 broadly regulated the alternative splicing of hundreds of genes with functions of protein binding, helicase activity, DNA binding and transcription co-activator. RBM4-regulated splicing of these genes plays an important role in apoptotic process and gene transcription regulation. As an example, exon inclusion of TNIP1 mediated by RBM4 affects the expression of its targets in inflammatory pathways. These results indicated that RBM4 can mediate the inflammatory response via splicing regulation, which adds to the understanding of the critical role of RBM4 in cancer complicated by inflammation. In conclusion, this study indicated a mechanism in which the dysregulation of alternative splicing can influence cellular biology and lead to various immune-related diseases.

circRNAs Are Involved in the Rice-Magnaporthe oryzae Interaction.

Circular RNAs (circRNAs) play roles in various biological processes, but their functions in the rice (Oryza sativa) response to Magnaporthe oryzae remain elusive. Here, we demonstrate that circRNAs are involved in the rice-M. oryzae interaction using comparative circRNA-sequencing and transgenic approaches. We identified 2932 high-confidence circRNAs from young leaves of the blast-resistant accession International Rice Blast Line Pyricularia-Kanto51-m-Tsuyuake (IR25) and the blast-susceptible accession Lijiangxin Tuan Heigu (LTH) under M oryzae-infected or uninfected conditions; 636 were detected specifically upon M oryzae infection. The circRNAs in IR25 were significantly more diverse than those in LTH, especially under M. oryzae infection. Particularly, the number of circRNAs generated per parent gene was much higher in IR25 than in LTH and increased in IR25 but decreased in LTH upon M. oryzae infection. The higher diversity of circRNAs in IR25 was further associated with more frequent 3' and 5' alternative back-splicing and usage of complex splice sites. Moreover, a subset of circRNAs was differentially responsive to M oryzae in IR25 and LTH. We further confirmed that circR5g05160 promotes rice immunity against M oryzae Therefore, our data indicate that circRNA diversity is associated with different responses to M oryzae infection in rice and provide a starting point to investigate a new layer of regulation in the rice-M oryzae interaction.

A draft genome assembly of halophyte Suaeda aralocaspica, a plant that performs C4 photosynthesis within individual cells.

BACKGROUND: The halophyte Suaeda aralocaspica performs complete C4 photosynthesis within individual cells (SCC4), which is distinct from typical C4 plants, which require the collaboration of 2 types of photosynthetic cells. However, despite SCC4 plants having features that are valuable in engineering higher photosynthetic efficiencies in agriculturally important C3 species such as rice, there are no reported sequenced SCC4 plant genomes, limiting our understanding of the mechanisms involved in, and evolution of, SCC4 photosynthesis. FINDINGS: Using Illumina and Pacific Biosciences sequencing platforms, we generated ∼202 Gb of clean genomic DNA sequences having a 433-fold coverage based on the 467 Mb estimated genome size of S. aralocaspica. The final genome assembly was 452 Mb, consisting of 4,033 scaffolds, with a scaffold N50 length of 1.83 Mb. We annotated 29,604 protein-coding genes using Evidence Modeler based on the gene information from ab initio predictions, homology levels with known genes, and RNA sequencing-based transcriptome evidence. We also annotated noncoding genes, including 1,651 long noncoding RNAs, 21 microRNAs, 382 transfer RNAs, 88 small nuclear RNAs, and 325 ribosomal RNAs. A complete (circular with no gaps) chloroplast genome of S. aralocaspica 146,654 bp in length was also assembled. CONCLUSIONS: We have presented the genome sequence of the SCC4 plant S. aralocaspica. Knowledge of the genome of S. aralocaspica should increase our understanding of the evolution of SCC4 photosynthesis and contribute to the engineering of C4 photosynthesis into economically important C3 crops.

The dynamics of FTO binding and demethylation from the m6A motifs.

N6-methyladenosine (m6A) is considered as a reversible RNA modification occurring more frequently on the GAC than AAC context in vivo, which regulates post-transcriptional gene expression in mammalian cells. m6A 'writers' METTL3 and METTL14 demonstrate a strong preference for binding AC-containing motifs in living cells. However, this evidence is currently lacking for m6A erasers, leaving the dynamics of the internal m6A modification under debate recently. We analysed three recently published FTO CLIP-seq data sets and two generated in this study, one of the two known m6A 'erasers'. FTO binding peaks from all cell lines contain RRACH motifs. Only those from K562, 3T3-L1and HeLa cells were enriched in AC-containing motifs, while those from HEK293 were not. The exogenously overexpressed FTO effectively binds to m6A motif-containing RNA sites. FTO overexpression specifically removed m6A modification from GGACU and RRACU motifs in a concentration-dependent manner. These findings underline the dynamics of FTO in target selection, which is predicted to contribute to both the m6A dynamics and the FTO plasticity in biological functions and diseases.

In vitro expansion affects the response of human bone marrow stromal cells to irradiation.

BACKGROUND: Bone marrow stromal cells (BMSCs) are extensively used in regeneration therapy and cytology experiments simulate how BMSCs respond to radiation. Due to the small number and the heterogeneity of primary isolated BMSCs, extensive in vitro expansion is usually required before application, which affects the cellular characteristics and gene expression of BMSCs. However, whether the radiation response of BMSCs changes during in vitro expansion is unclear. METHODS: In this study, BMSCs were passaged in vitro and irradiated at passage 6 (P6) and passage 10 (P10). Then, apoptosis, the cell cycle, senescence, the cytokine secretion and the gene expression profile were analysed for the P6, P10, and non-irradiated (control) BMSCs at different post-irradiation time points. RESULTS: The P6 BMSCs had a lower percentage of apoptotic cells than the P10 BMSCs at 24 and 48 h post-irradiation but not compared to that of the controls at 2 and 8 h post-irradiation. The P6 BMSCs had a lower percentage of cells in S phase and a higher percentage in G1 phase than the P10 BMSCs at 2 and 8 h post-irradiation. The radiation had similar effects on the senescent cell level and impaired immunomodulation capacity of the P6 and P10 BMSCs. Regardless of whether they were irradiated, the P6 and P10 BMSCs always expressed a distinctive set of genes. The upregulated genes were enriched in pathways including the cell cycle, DNA replication and oocyte meiosis. Then, a subset of conserved irradiation response genes across the BMSCs was identified, comprising 12 differentially upregulated genes and 5 differentially downregulated genes. These genes were especially associated with the p53 signaling pathway, DNA damage and DNA repair. Furthermore, validation experiments revealed that the mRNA and protein levels of these conserved genes were different between the P6 and P10 BMSCs after irradiation. Weighted gene co-expression network analysis supported these findings and further revealed the effects of cell passage on the irradiation response in BMSCs. CONCLUSION: The results indicated that cell passage in vitro affected the irradiation response of BMSCs via molecular mechanisms that mediated differences in apoptosis, the cell cycle, senescence and the cytokine secretion. Thus, accurate cell passage information is not only important for transplantation therapy but also for future studies on the radiation response in BMSCs.

Divergence in larval diapause induction between the rice and water-oat populations of the striped stem borer, Chilo suppressalis (Walker) (Lepidoptera: Crambidae).

Differences in diapause traits can result in the seasonal reproductive isolation of host plant-associated insect populations and thereby facilitate the population divergence. The striped stem borer, Chilo suppressalis, has two host plant-associated populations: rice population and water-oat population. Several studies have found evidence that seasonal reproductive isolation between these populations is at least partially due to interpopulation differences in diapause. However, there still lack unambiguous evidence comparing characteristics of diapause induction for both populations. We compared the photoperiodic response and the age of peak photoperiod sensitivity of these populations and used RNA-Seq to compare the molecular response of diapause induction between populations. The photoperiodic response of the two populations differed at 25 °C; the critical night length of larvae from the rice population was 11 h and 20 min, whereas no obvious critical night length was in those from the water-oat population. In rice population, larvae were most sensitive to photoperiod at 9-12 days of age, whereas in water-oat population, larvae were the most sensitive to photoperiod at 9-10 days of age. The RNA-Seq results indicated that there were several differences in the molecular response of diapause induction and small overlap in differentially expressed genes (DEGs) between populations. Furthermore, GO analysis indicated that both rice and water-oat population's DEGs were significantly enriched in heme and iron binding. Besides, water-oat population's DEGs were significantly enriched in metabolizing nutrients but rice population's DEGs do not. Thus, our results described differences in diapause induction between rice and water-oat populations of C. suppressalis which could affect the timing of diapause and thereby contribute to the seasonal reproductive isolation of these host plant-associated populations. In conclusion, this work suggests that difference in diapause induction could promote the population divergence in insects associated with different host plants.

Difference in diel mating time contributes to assortative mating between host plant-associated populations of Chilo suppressalis.

Behavioral isolation in animals can be mediated by inherent mating preferences and assortative traits, such as divergence in the diel timing of mating activity. Although divergence in the diel mating time could, in principle, promote the reproductive isolation of sympatric, conspecific populations, there is currently no unequivocal evidence of this. We conducted different mate-choice experiments to investigate the contribution of differences in diel mating activity to the reproductive isolation of the rice and water-oat populations of Chilo suppressalis. The results show that inter-population difference in diel mating activity contributes to assortative mating in these populations. In the rice population, most mating activity occurred during the first half of the scotophase, whereas in the water-oat population virtually all mating activity was confined to the second half of the scotophase. However, when the photoperiod of individuals from the water-oat population was altered to more closely align their mating activity with that of the rice population, mate choice was random. We conclude that inter-population differences in diel mating time contribute to assortative mating, and thereby the partial reproductive isolation, of these host-associated populations of C. suppressalis.

Thermal effects on development and adult longevity of endoparasitoid Chelonus murakatae Munakata (Hymenoptera: Braconidae).

The effects of temperature on the development duration and longevity of adult of Chelonus murakatae were studied under five constant temperatures including 22.5, 25, 27.5, 30 and 32.5 ± 0.5 °C under laboratory conditions. It was observed that the development time was inversely proportional to the temperature within the range of 22.5 to 32.5 °C. The results indicated that the optimum temperature for development ranged from 25 to 30 °C. Thermal threshold was estimated by a linear model which was recorded as 15.5 and 18.5 °C for males and females, respectively. Number of degree days required to complete the development from egg to adult were 439.6 degree days in males and 336.8 degree days in females. Adult longevity also decreased with increase in temperature. This information can be used for optimizing mass culturing and field release for an efficient biological control of Chilo suppressalis in this specie.

Do differences in life-history traits and the timing of peak mating activity between host-associated populations of Chilo suppressalis have a genetic basis?

The development of host races, genetically distinct populations of the same species with different hosts, is considered to be the initial stage of ecological speciation. Ecological and biological differences consistent with host race formation have been reported between water-oat and rice-associated populations of Chilo suppressalis. In order to confirm whether these differences have a genetic basis, we conducted experiments to determine the extent to which various life-history traits and the time of peak mating activity of these populations were influenced by the species of host plant larvae were raised on. Individuals from each population were reared for three consecutive generations on either water-oat fruit pulp or rice seedlings. Descendants of both populations had higher larval survival rates, shorter larval developmental periods, higher pupal weight, and longer adult forewings, when reared on water-oats than when reared on rice. The time of peak of mating activity differed between the descendants of each population, irrespective of whether they were raised on water-oats or rice. These results indicate that although some life-history traits of host-associated populations of C. suppressalis are influenced by the host plant larvae are raised on, time of peak mating activity is not. Because it is a stable, objective, phenotypic trait, further research on difference in the time of peak mating activity between host-associated populations of C. suppressalis should be conducted to clarify the mechanism responsible for host race formation in this species.