Co-treatment with melatonin and ortho-topolin riboside reduces cell viability by altering metabolic profiles in non-small cell lung cancer cells

Lung cancer is a highly prevalent and lethal malignancy worldwide, with non-small cell lung cancer (NSCLC) accounting for 85% of cancer-related deaths. In this study, the effects of co-treatment with melatonin and ortho-topolin riboside (oTR) on the cell viability and alteration of metabolites and transcripts were investigated in NSCLC cells using gas chromatography-mass spectrometry (GC-MS) and next-generation sequencing (NGS). The co-treatment of melatonin and oTR exhibited synergistic effects on the reduction of cell viability and alteration of metabolic and transcriptomic profiles in NSCLC cells. We observed that the co-treatment inhibited glycolytic function and mitochondria respiration, and downregulated glycine, serine and threonine metabolism alongside tyrosine metabolism in NSCLC cells. In the glycine, serine and threonine metabolism pathway, the co-treatment resulted in a significant 8.4-fold reduction in the expression level of the SDS gene, which encodes the enzyme responsible for the breakdown of serine to pyruvate. Moreover, co-treatment decreased the gene expression of TH, DDC, and CYP1A1 in tyrosine metabolism. Additionally, we observed that the co-treatment resulted in a significant 146.9-fold reduction in the expression of the DISC1 gene. The alteration in metabolites and transcript expressions might provide information to explain the cytotoxicity of co-treatment of melatonin and oTR in NSCLC cells. Our study presents insights into the synergistic anticancer effect of the co-treatment of melatonin and oTR, which could be a potential future therapeutic strategy for the treatment of NSCLC patients.

Identification of MCM4 and PRKDC as new regulators of osteosarcoma cell dormancy based on 3D cell cultures

Dormancy is a potential way for tumors to develop drug resistance and escape treatment. However, the mechanisms involved in cancer dormancy remain poorly understood. This is mainly because there is no in vitro culture model making it possible to spontaneously induce dormancy. In this context, the present work proposes the use of three-dimensional (3D) spheroids developed from osteosarcoma cell lines as a relevant model for studying cancer dormancy. MNNG-HOS, SaOS-2, 143B, MG-63, U2OS and SJSA-1 cell lines were cultured in 3D using the Liquid Overlay Technique (LOT). Dormancy was studied by staining cancer cells with a lipophilic dye (DiD), and long-term DiD+ cells were considered as dormant cancer cells. The role of the extracellular matrix in inducing dormancy was investigated by embedding cells into methylcellulose or Geltrex™. Gene expression of DiD+ cells was assessed with a Nanostring™ approach and the role of the genes detected in dormancy was validated by a transient down-expression model using siRNA treatment. Proliferation was measured using fluorescence microscopy and the xCELLigence technology. We observed that MNNG-HOS, 143B and MG-G3 cell lines had a reduced proliferation rate in 3D compared to 2D. U2OS cells had an increased proliferation rate when they were cultured in Geltrex™ compared to other 3D culture methods. Using 3D cultures, a transcriptomic signature of dormancy was obtained and showed a decreased expression of 18 genes including ETV4, HELLS, ITGA6, MCM4, PRKDC, RAD21 and UBE2T. The treatment with siRNA targeting these genes showed that cancer cell proliferation was reduced when the expression of ETV4 and MCM4 were decreased, whereas proliferation was increased when the expression of RAD21 was decreased. 3D culture facilitates the maintenance of dormant cancer cells characterized by a reduced proliferation and less differential gene expression as compared to proliferative cells. Further studies of the genes involved has enabled us to envisage their role in regulating cell proliferation.

Is ketogenic diet a 'precision medicine'? Recent developments and future challenges

Recently, precision medicine has attracted much attention in the management of epilepsies, but it remains unclear if the increasingly utilized ketogenic diet approaches can truly be considered precision medicine in all epilepsy treatment. Currently, it is the standard treatment for patients with GLUT1 deficiency and the latest NICE guidelines highlight ketogenic diet as a therapeutic option for multi-drug resistant epilepsy patients. Ketogenic diet is presumed to be a precision medicine tool when applied to the treatment of seizures secondary to GLUT1 transporter deficiency. In contrast, the genetic and epigenetic mechanisms modulated by ketogenic diet and underlying its efficacy in other epilepsy types can only be hypothesized to relate to mechanisms of neuroprotection, neuromodulation, and reduction of neuroinflammation. Early ketogenic diet initiation in well-selected patients, would allow immediate action in the direction of neuroprotection and modulation of neuroinflammation, ensuring higher success rates and lower "cost" to the patient in terms of quality of life and comorbidities. These considerations have fueled an increasing interest in investigating the efficacy, side effects, and adherence to long-term use of the ketogenic diet in epilepsy treatment in large contemporary cohorts, available within the scope of multicentric collaborations, such as the European Network for Therapy in Rare Epilepsies (NETRE). Future directions should involve the use of precision medicine, applied to each patient with the help of "omics", whose use should be expanded and inclusive.

Identification of shared disease marker genes and underlying mechanisms between major depression and rheumatoid arthritis

Both depression and rheumatoid arthritis (RA) have a very high comorbidity rate. A bilateral association is estimated to increase the mutual risk and the common denominator is inflammation being observed in both diseases. Previous studies have mainly focused on assessing peripheral blood's inflammatory and pro-inflammatory cytokines levels. We aimed to extend insights into the molecular mechanisms of depression based on hub RA genes. To do so, we prioritized RA-related genes using in-silico tools. We then investigated whether RA-related genes undergo altered expression in patients with unipolar and bipolar depression without a concurrent RA diagnosis and any exponents of active inflammation. In addition, we selected a homogeneous group of patients treated with lithium (Li), which has immunomodulatory properties. The study was performed on patients with bipolar depression (BD, n = 45; Li, n = 20), unipolar depression (UD, n = 27), and healthy controls (HC, n = 22) of both sexes. To identify DEGs in peripheral blood mononuclear cells (PBMCs), we used the SurePrint G3 Microarray and GeneSpring software. We selected a list of 180 hub genes whose altered expression we analyzed using the expression microarray results. In the entire study group, we identified altered expression of 93 of the 180 genes, including 35 down-regulated (OPRM1 gene with highest FC > 3) and 58 up-regulated (TLR4 gene with highest FC > 3). In UD patients, we observed maximally up-regulated expression of the TEK gene (FC > 3), and in BD of the CXCL8 gene (FC > 5). On the other hand, in lithium-treated patients, the gene with the most reduced expression was the TRPV1 gene. The study proved that depression and RA are produced by a partially shared "inflammatory interactome" in which the opioid and angiogenesis pathways are important.

Identification of RNA reads encoding different channels in isolated rat ventricular myocytes and the effect of cell stretching on L-type Ca2+current

BACKGROUND

The study aimed to identify transcripts of specific ion channels in rat ventricular cardiomyocytes and determine their potential role in the regulation of ionic currents in response to mechanical stimulation. The gene expression levels of various ion channels in freshly isolated rat ventricular cardiomyocytes were investigated using the RNA-seq technique. We also measured changes in current through CaV1.2 channels under cell stretching using the whole-cell patch-clamp method.

RESULTS

Among channels that showed mechanosensitivity, significant amounts of TRPM7, TRPC1, and TRPM4 transcripts were found. We suppose that the recorded L-type Ca2+ current is probably expressed through CaV1.2. Furthermore, stretching cells by 6, 8, and 10 μm, which increases ISAC through the TRPM7, TRPC1, and TRPM4 channels, also decreased ICa,L through the CaV1.2 channels in K+ in/K+ out, Cs+ in/K+ out, K+ in/Cs+ out, and Cs+ in/Cs+ out solutions. The application of a nonspecific ISAC blocker, Gd3+, during cell stretching eliminated ISAC through nonselective cation channels and ICa,L through CaV1.2 channels. Since the response to Gd3+ was maintained in Cs+ in/Cs+ out solutions, we suggest that voltage-gated CaV1.2 channels in the ventricular myocytes of adult rats also exhibit mechanosensitive properties.

CONCLUSIONS

Our findings suggest that TRPM7, TRPC1, and TRPM4 channels represent stretch-activated nonselective cation channels in rat ventricular myocytes. Probably the CaV1.2 channels in these cells exhibit mechanosensitive properties. Our results provide insight into the molecular mechanisms underlying stretch-induced responses in rat ventricular myocytes, which may have implications for understanding cardiac physiology and pathophysiology.

Effects of lactic acid bacteria-fermented formula milk supplementation on ileal microbiota, transcriptomic profile, and mucosal immunity in weaned piglets

BACKGROUND

Lactic acid bacteria (LAB) participating in milk fermentation naturally release and enrich the fermented dairy product with a broad range of bioactive metabolites, which has numerous roles in the intestinal health-promoting of the consumer. However, information is lacking regarding the application prospect of LAB fermented milk in the animal industry. This study investigated the effects of lactic acid bacteria-fermented formula milk (LFM) on the growth performance, intestinal immunity, microbiota composition, and transcriptomic responses in weaned piglets. A total of 24 male weaned piglets were randomly divided into the control (CON) and LFM groups. Each group consisted of 6 replicates (cages) with 2 piglets per cage. Each piglet in the LFM group were supplemented with 80 mL LFM three times a day, while the CON group was treated with the same amount of drinking water.

RESULTS

LFM significantly increased the average daily gain of piglets over the entire 14 d (P < 0.01) and the average daily feed intake from 7 to 14 d (P < 0.05). Compared to the CON group, ileal goblet cell count, villus-crypt ratio, sIgA, and lactate concentrations in the LFM group were significantly increased (P < 0.05). Transcriptomic analysis of ileal mucosa identified 487 differentially expressed genes (DEGs) between two groups. Especially, DEGs involved in the intestinal immune network for IgA production pathways, such as polymeric immunoglobulin receptor (PIGR), were significantly up-regulated (P < 0.01) by LFM supplementation. Moreover, trefoil factor 2 (TFF2) in the LFM group, one of the DEGs involved in the secretory function of goblet cells, was also significantly up-regulated (P < 0.01). Sequencing of the 16S rRNA gene of microbiota demonstrated that LFM led to selective enrichment of lactate-producing and short-chain fatty acid (SCFA)-producing bacteria in the ileum, such as an increase in the relative abundance of Enterococcus (P = 0.09) and Acetitomaculum (P < 0.05).

CONCLUSIONS

LFM can improve intestinal health and immune tolerance, thus enhancing the growth performance of weaned piglets. The changes in microbiota and metabolites induced by LFM might mediate the regulation of the secretory function of goblet cells.

Ectopic expression of a combination of 5 genes detects high risk forms of T-cell acute lymphoblastic leukemia

BACKGROUND

T cell acute lymphoblastic leukemia (T-ALL) defines a group of hematological malignancies with heterogeneous aggressiveness and highly variable outcome, making therapeutic decisions a challenging task. We tried to discover new predictive model for T-ALL before treatment by using a specific pipeline designed to discover aberrantly active gene.

RESULTS

The expression of 18 genes was significantly associated with shorter survival, including ACTRT2, GOT1L1, SPATA45, TOPAZ1 and ZPBP (5-GEC), which were used as a basis to design a prognostic classifier for T-ALL patients. The molecular characterization of the 5-GEC positive T-ALL unveiled specific characteristics inherent to the most aggressive T leukemic cells, including a drastic shut-down of genes located on the mitochondrial genome and an upregulation of histone genes, the latter characterizing high risk forms in adult patients. These cases fail to respond to the induction treatment, since 5-GEC either predicted positive minimal residual disease (MRD) or a short-term relapse in MRD negative patients.

CONCLUSION

Overall, our investigations led to the discovery of a homogenous group of leukemic cells with profound alterations of their biology. It also resulted in an accurate predictive tool that could significantly improve the management of T-ALL patients.

Determining oncogenic patterns and cancer predisposition through the transcriptomic profile in Mitchell-Riley syndrome with heterotopic gastric mucosa and duodenal atresia: a case report

Background

Homozygous mutations in the transcription factor RFX6 are the cause of the Mitchell–Riley syndrome (MRS) associating neonatal diabetes, congenital digestive system, such as biliary atresia, pancreatic hypoplasia, duodenal and/or jejunal atresia, intestinal malrotation, gallbladder aplasia, cholestasis. A constitutive inactivation of RFX6 leads also to gastric heterotopia. Application of RNA-seq in human diseases may help to better understand pathogenic mechanism of diseases and to predict the risk of developing chronic disorders and personalizing their prevention and treatment. We evaluated oncogenic patterns and cancer predisposition using the transcriptomic profile in a case of MRS with neonatal diabetes, duodenal atresia, and extensive intestinal tract gastric heterotopia.

Results

We signalled the interactors of RFX6 with other up and downregulated genes, that may be interested in severity of diabetic condition, in multi-organs impairment and cancer predisposition. Furthermore, several dysregulated genes are involved in biological processes that can lead to promote cancer including “Evading apoptosis” (BAD, BBC3, EGF, FGFR2, FLT3LG, HMOX1, HRAS, IFNAR2, IGF1R, IL12RB1, IL13RA1, IL15, IL2RB, IL2RG, IL6R, KEAP1, MGST1, PDGFA, PDGFRB, PIK3R3, RALB, RALGDS, RASSF1, SOS1, TGFA, TXNRD3), “Proliferation” (APC, BRAF, CCND2, CCND3, CCNE2, FGFR2, FLT3LG, FZD1, FZD6, HMOX1, HRAS, IGF1R, KEAP1, LRP6, MAPK3, MGST1, PDGFA, PDGFB, PDGFRB, RB1, SOS1, TGFA, TXNRD3, WNT10B), “Sustained angiogenesis” (BRAF, FGFR2, FLT3LG, HRAS, IGF1R, JAG1, MAPK3, NOTCH2, PDGFA, PDGFB, PDGFRB, SOS1, TGFA, TGFB1), “Genomic instability” (BAD, BBC3) and “Insensitivity to anti-growth signals” (SMAD2, TGFB1). We also inspected the signalings and their related genes in cancer, such as “PI3K signaling”, “ERK signaling”, “JAK-STAT signaling”, “Calcium signaling”, “Other RAS signaling”, “WNT signaling”.

Conclusions

In our MRS patient, we signaled the interactors of RFX6 with other up- and downregulated genes that may be related to severe diabetic condition, multi-organ impairment, and cancer predisposition. Notably, many dysregulated genes may lead to triggering carcinogenesis. The possibility of the patient developing cancer degeneration in heterotopic gastric mucosa and/or additional long-term tumoral sequelae is not excluded. Personalized prevention and treatment strategies should be proposed.

Comparative analysis of the transcriptional responses of five Leishmania species to trivalent antimony

BACKGROUND

Leishmaniasis is a neglected tropical disease caused by several species of Leishmania. The resistance phenotype of these parasites depends on the characteristics of each species, which contributes to increased therapeutic failures. Understanding the mechanism used by the parasite to survive under treatment pressure in order to identify potential common and specific therapeutic targets is essential for the control of leishmaniasis. The aim of this study was to investigate the expression profiles and potential shared and specific resistance markers of the main Leishmania species of medical importance [subgenus L. (Leishmania): L. donovani, L. infantum and L. amazonensis; subgenus L. (Viannia): L. panamensis and L. braziliensis)] resistant and sensitive to trivalent stibogluconate (SbIII).

METHODS

We conducted comparative analysis of the transcriptomic profiles (only coding sequences) of lines with experimentally induced resistance to SbIII from biological replicates of five Leishmania species available in the databases of four articles based on ortholog attribution. Simultaneously, we carried out functional analysis of ontology and reconstruction of metabolic pathways of the resulting differentially expressed genes (DEGs).

RESULTS

Resistant lines for each species had differential responses in metabolic processes, compound binding, and membrane components concerning their sensitive counterpart. One hundred and thirty-nine metabolic pathways were found, with the three main pathways comprising cysteine and methionine metabolism, glycolysis, and the ribosome. Differentially expressed orthologous genes assigned to species-specific responses predominated, with 899 self-genes. No differentially expressed genes were found in common among the five species. Two common upregulated orthologous genes were found among four species (L. donovani, L. braziliensis, L. amazonensis, and L. panamensis) related to an RNA-binding protein and the NAD(P)H cytochrome-B5-oxidoreductase complex, associated with transcriptional control and de novo synthesis of linoleic acid, critical mechanisms in resistance to antimonials.

CONCLUSION

Herein, we identified potential species-specific genes related to resistance to SbIII. Therefore, we suggest that future studies consider a treatment scheme that is species-specific. Despite the limitations of our study, this is the first approach toward unraveling the pan-genus genetic mechanisms of resistance in leishmaniasis.

Transcriptomic Analysis of Peripheral Monocytes upon Fingolimod Treatment in Relapsing Remitting Multiple Sclerosis Patients

Fingolimod (FTY), a second-line oral drug approved for relapsing remitting Multiple Sclerosis (RRMS) acts in preventing lymphocyte migration outside lymph nodes; moreover, several lines of evidence suggest that it also inhibits myeloid cell activation. In this study, we investigated the transcriptional changes induced by FTY in monocytes in order to better elucidate its mechanism of action. CD14⁺ monocytes were collected from 24 RRMS patients sampled at baseline and after 6 months of treatment and RNA profiles were obtained through next-generation sequencing. We conducted pathway and sub-paths analysis, followed by centrality analysis of cell-specific interactomes on differentially expressed genes (DEGs). We investigated also the predictive role of baseline monocyte transcription profile in influencing the response to FTY therapy. We observed a marked down-regulation effect (60 down-regulated vs. 0 up-regulated genes). Most of the down-regulated DEGs resulted related with monocyte activation and migration like IL7R, CCR7 and the Wnt signaling mediators LEF1 and TCF7. The involvement of Wnt signaling was also confirmed by subpaths analyses. Furthermore, pathway and network analyses showed an involvement of processes related to immune function and cell migration. Baseline transcriptional profile of the HLA class II gene HLA-DQA1 and HLA-DPA1 were associated with evidence of disease activity after 2 years of treatment. Our data support the evidence that FTY induces major transcriptional changes in monocytes, mainly regarding genes involved in cell trafficking and immune cell activation. The baseline transcriptional levels of genes associated with antigen presenting function were associated with disease activity after 2 years of FTY treatment.

B cells acquire a unique and differential transcriptomic profile during pregnancy

Pregnancy alters B cell development and function. B cell activation is initiated by antigens binding to the BCR leading to B cell survival, proliferation, antigen presentation and antibody production. We performed a genome-wide transcriptome profiling of splenic B cells from pregnant (P) and non-pregnant (NP) mice and identified 1136 genes exhibiting differential expression in B cells from P mice (625 up- and 511 down-regulated) compared to NP animals. In silico analysis showed that B cell activation through BCR seems to be lowered during pregnancy. RT-qPCR analysis confirmed these data. Additionally, B cells from pregnant women stimulated in vitro through BCR produced lower levels of inflammatory cytokines compared to non-pregnant women. Our results suggest that B cells acquire a state of hypo-responsiveness during gestation, probably as part of the maternal immune strategy for fetal tolerance but also open new avenues to understand why pregnant women are at highest risk for infections.

Characterization of terminal-ileal and colonic Crohn's disease in treatment-naïve paediatric patients based on transcriptomic profile using logistic regression

BACKGROUND

Inflammatory bowel disease (IBD) is a chronic and idiopathic inflammatory disorder of the gastrointestinal tract and comprises ulcerative colitis (UC) and Crohn's disease (CD). Crohn's disease can affect any part of the gastrointestinal tract, but mainly the terminal ileum and colon. In the present study, we aimed to characterize terminal-ileal CD (ICD) and colonic CD (CCD) at the molecular level, which might enable a more optimized approach for the clinical care and scientific research of CD.

METHODS

We analyzed differentially expressed genes in samples from 23 treatment-naïve paediatric patients with CD and 25 non-IBD controls, and compared the data with previously published RNA-Seq data using multi-statistical tests and confidence intervals. We implemented functional profiling and proposed statistical methods for feature selection using a logistic regression model to identify genes that are highly associated in ICD or CCD. We also validated our final candidate genes in independent paediatric and adult cohorts.

RESULTS

We identified 550 genes specifically expressed in patients with CD compared with those in healthy controls (p < 0.05). Among these DEGs, 240 from patients with CCD were mainly involved in mitochondrial dysfunction, whereas 310 from patients with ICD were enriched in the ileum functions such as digestion, absorption, and metabolism. To choose the most effective gene set, we selected the most powerful genes (p-value ≤ 0.05, accuracy ≥ 0.8, and AUC ≥ 0.8) using logistic regression. Consequently, 33 genes were identified as useful for discriminating CD location; the accuracy and AUC were 0.86 and 0.83, respectively. We then validated the 33 genes with data from another independent paediatric cohort (accuracy = 0.93, AUC = 0.92) and adult cohort (accuracy = 0.88, AUC = 0.72).

CONCLUSIONS

In summary, we identified DEGs that are specifically expressed in CCD and ICD compared with those in healthy controls and patients with UC. Based on the feature selection analysis, 33 genes were identified as useful for discriminating CCD and ICD with high accuracy and AUC, for not only paediatric patients but also independent cohorts. We propose that our approach and the final gene set are useful for the molecular classification of patients with CD, and it could be beneficial in treatments based on disease location.

Alterations in 'inflammatory' pathways in the rat prefrontal cortex as early biological predictors of the long-term negative consequences of exposure to stress early in life

Early life stress, especially when experienced during the first period of life, affects the brain developmental trajectories leading to an enhanced vulnerability for stress-related psychiatric disorders later in life. Although both clinical and preclinical studies clearly support this association, the biological pathways deregulated by such exposure, and the effects in shaping the neurodevelopmental trajectories, have so far been poorly investigated. By using the prenatal stress (PNS) model, a well-established rat model of early life stress, we performed transcriptomic analyses in the prefrontal cortex of rats exposed or not to PNS and sacrificed at different postnatal days (PNDs 21, 40, 62). We first investigated the long-lasting mechanisms and pathways affected in the PFC. We have decided to focus on the prefrontal cortex because we have previously shown that this brain region is highly sensitive to PNS exposure. We found that adult animals exposed to PNS show alterations in 389 genes, mainly involved in stress and inflammatory signalling. We then wanted to establish whether PNS exposure could also affect the neurodevelopmental trajectories in order to identify the most critical temporal window. We found that PNS rats show the most significant changes during adolescence (between PND 40 versus PND 21), with alterations of several pathways related to stress, inflammation and metabolism, which were maintained until adulthood.

Influence of Vitamin A deficiency on the transcriptomic profile of rat meibomian glands

PURPOSE

Vitamin A deficiency (VAD) is associated with chalazion in young children. However, the underlying molecular mechanism remains unclear. In the present study, transcriptome data from rat meibomian glands (MGs) were analyzed to reveal specific molecular responses to VAD.

METHODS

Total RNA was extracted and purified for library preparation and transcriptome sequencing. Differentially expressed genes (DEGs) between vitamin A normal (VAN) and VAD rats were analyzed using DESeq software. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of DEGs were performed using the GO seq R package and KOBAS software. Real-time quantitative reverse transcription polymerase chain reaction was used to validate the RNA sequencing results.

RESULTS

The number of DEGs in the VAD group compared to the VAN group was 3129 (1531 upregulated and 1598 downregulated) in the rat MGs. VAD upregulated a large number of lipid metabolism-related genes. GO analysis showed that the most enriched and meaningful terms were related to lipid metabolism (e.g., "oxidation-reduction process, GO: 0,055,114," "lipid metabolic process, GO: 000,662"). KEGG pathway analysis showed that most of the enriched signaling pathways were involved in lipid metabolism, including the PPAR signaling pathway associated with retinoic acid (RA)-mediated nuclear receptors.

CONCLUSION

These findings demonstrate that VAD regulates the expression of numerous genes in the rat MG and that many of these genes are involved in lipid metabolic pathways.

Global Transcriptome Profiling of Enterobacter Strain NRS-1 in Response to Hydrogen Peroxide Stress Treatment

Microbes are often subjected to oxidative stress in nature that badly affects their growth rate and viability. Although the response of microbes against oxidative stress has been characterized at the chemical, physiological, and molecular levels, the mechanism of gene-regulation network adaptations of bacteria in response to oxidative stress remains largely unknown. In this study, transcriptomic profiling of glyphosate-tolerant Enterobacter strain NRS-1 was analyzed under 9 mM H2O2 stress using RNA-seq and qRT-PCR. The lag period in the growth of NRS-1 was very short compared with wild-type strain under H2O2 treatment. A total of 113 genes are identified as differentially expressed genes (DEGs) under H2O2 that include 38 upregulated and 75 downregulated transcripts. But not any genes regulated by major oxidative regulons, viz., oxyR, soxR, rpoS, perR, ohrR, and σв, have been reported in DEGs, hence potentially reflecting that specific changes have occurred in NRS-1 for adaptation to oxidative stress. Based on the functions of the DEGs, six elements namely formate dehydrogenase, processes associated with iron ions, repair programs, multidrug resistance, antioxidant defense, and energy generation (mqo, sdhC) might have contributed for stress tolerance in NRS-1. These elements are proposed to form a molecular network explaining gene response of NRS-1 to stress, and ensure global cell protection and growth recovery of NRS-1. These findings enrich the view of gene regulation in bacteria in response to H2O2 oxidative stress.

A high fat diet with a high C18:0/C16:0 ratio induced worse metabolic and transcriptomic profiles in C57BL/6 mice

BACKGROUND

Differential effects of individual saturated fatty acids (SFAs), particularly stearic acid (C18:0), relative to the shorter-chain SFAs have drawn interest for more accurate nutritional guidelines. However, specific biologic and pathologic functions that can be assigned to particular SFAs are very limited. The present study was designed to compare changes in metabolic and transcriptomic profiles in mice caused by a high C18:0 diet and high palmitic acid (C16:0) diet.

METHODS

Male C57BL/6 mice were assigned to a normal fat diet (NFD), a high fat diet with high C18:0/C16:0 ratio (HSF) or an isocaloric high fat diet with a low C18:0/C16:0 ratio (LSF) for 10 weeks. An oral glucose tolerance test, 72-h energy expenditure measurement and CT scan of body fat were done before sacrifice. Fasting glucose and lipids were determined by an autobiochemical analyzer. Blood insulin, tumor necrosis factor-α (TNF-α), and interleukin-6 (IL-6) levels were measured by enzyme-linked immunosorbent assay methods. Free fatty acids (FFAs) profiles in blood and liver were determined by using gas chromatography-mass spectrometry. Microarray analysis was applied to investigate changes in transcriptomic profiles in the liver. Pathway analysis and gene ontology analysis were applied to describe the roles of differentially expressed mRNAs.

RESULTS

Compared with the NFD group, body weight, body fat ratio, fasting blood glucose, insulin, homeostasis model assessment of insulin resistance (HOMA-IR), triglyceride, IL-6, serum and liver FFAs including total FFAs, C16:0 and C18:0 were increased in both high fat diet groups and were much higher in the HSF group than those in the LSF group. Both HSF and LSF mice exhibited distinguishable long non-coding RNA (lncRNA), microRNA and mRNA expression profiles when compared with those of NFD mice. Additionally, more differentially expressed lncRNAs and mRNAs were observed in the HSF group than in the LSF group. Some biological functions and pathways, other than energy metabolism regulation, were identified as differentially expressed mRNAs between the HSF group and the LSF group.

CONCLUSION

The high fat diet with a high C18:0/C16:0 ratio induced more severe glucose and lipid metabolic disorders and inflammation and affected expression of more lncRNAs and mRNAs than an isocaloric low C18:0/C16:0 ratio diet in mice. These results provide new insights into the differences in biological functions and related mechanisms, other than glucose and lipid metabolism, between C16:0 and C18:0.

Transcriptomic profiling disclosed the role of DNA methylation and histone modifications in tumor-infiltrating myeloid-derived suppressor cell subsets in colorectal cancer

Background

Increased numbers of myeloid-derived suppressor cells (MDSCs) are positively correlated with poor prognosis and reduced survivals of cancer patients. They play central roles in tumor immune evasion and tumor metastasis. However, limited data are available on phenotypic/transcriptomic characteristics of the different MDSCs subsets in cancer. These cells include immature (I-MDSCs), monocytic (M-MDSCs), and polymorphonuclear/granulocytic (PMN-MDSCs).

Methods

Phenotypic characterization of myeloid subsets from 27 colorectal cancer (CRC) patients was assessed by flow cytometric analyses. RNA-sequencing of sorted I-MDSCs, PMN-MDSCs, and antigen-presenting cells (APCs) was also performed.

Results

We found that the levels of I-MDSCs and PMN-MDSCs were increased in tumor tissues (TT), compared with normal tissues (NT) in colorectal cancer. Our functional annotation analyses showed that genes associated with histone deacetylase (HDAC) activation- and DNA methylation-mediated transcriptional silencing were upregulated, and histone acetyl transferase (HAT)-related genes were downregulated in tumor-infiltrating I-MDSCs. Moreover, pathways implicated in cell trafficking and immune suppression, including Wnt, interleukin-6 (IL-6), and mitogen-activated protein kinase (MAPK) signaling, were upregulated in I-MDSCs. Notably, PMN-MDSCs showed downregulation in genes related to DNA methylation and HDAC binding. Using an ex vivo model, we found that inhibition of HDAC activation or neutralization of IL-6 in CRC tumor tissues downregulates the expression of genes associated with immunosuppression and myeloid cell chemotaxis, confirming the importance of HDAC activation and IL-6 signaling pathway in MDSC function and chemotaxis.

Conclusions

This study provides novel insights into the epigenetic regulations and other molecular pathways in different myeloid cell subsets within the CRC tumor microenvironment (TME), giving opportunities to potential targets for therapeutic benefits.

Predictive early gene signature during mouse Bhas 42 cell transformation induced by synthetic amorphous silica nanoparticles

Synthetic amorphous silica nanoparticles (SAS) are used widely in industrial applications. These nanoparticles are not classified for their carcinogenicity in humans. However, some data still demonstrate a potential carcinogenic risk of these compounds in humans. The Bhas 42 cell line was developed to screen chemicals, as tumor-initiators or -promoters according to their ability to trigger cell-to-cell transformation, in a cell transformation assay. In the present study, we performed unsupervised transcriptomic analysis after exposure of Bhas 42 cells to NM-203 SAS as well as to positive (Min-U-Sil 5® crystalline silica microparticles, and 12-O-tetradecanoylphorbol-13-acetate) and negative (diatomaceous earth) control compounds. We identified a common gene signature for 21 genes involved in the early stage of the SAS- Min-U-Sil 5®- or TPA-induced cell transformation. These genes were related to cell proliferation (over expression) and cell adhesion (under expression). Among them, 12 were selected on the basis of their potential impact on cell transformation. RT-qPCR and western blotting were used to confirm the transcriptomic data. Moreover, similar gene alterations were found when Bhas 42 cells were treated with two other transforming SAS. In conclusion, the results obtained in the current study highlight a 12-gene signature that could be considered as a potential early "bio-marker" of cell transformation induced by SAS and perhaps other chemicals.

Physiological and transcriptomic analysis reveals the roles of 1-MCP in the ripening and fruit aroma quality of banana fruit (Fenjiao)

Fenjiao (Musa ABB Pisang Awak) is a popular banana cultivar due to its good taste and stress resistance, but it has a short shelf-life and deteriorates rapidly post-harvest. The effects of 1-methylcyclopropene (1-MCP) treatment on fruit physiology and quality and transcriptomic profiles are investigated in this study. The results showed that 1-MCP significantly delayed fruit ripening by repressing fruit softening and inhibiting the respiratory rate and ethylene production. The 1-MCP treatment delayed sugar accumulation and influenced the content of the precursors of the biosynthesis of aroma volatiles. 1-MCP reduced the production of flavor-contributing volatile esters isoamyl isobutyrate, isoamyl acetate and trans-2-hexenal and hexanal, but dramatically increased the hexyl acetate production at the full-ripening stage. The transcriptomic analysis showed that 1-MCP dramatically affected the transcript profiles during fruit ripening, especially the KEGG pathways involved in amino acid metabolism, biosynthesis of other secondary metabolites, carbohydrate metabolism, lipid metabolism, signal transduction, and translation classes. The key genes and the corresponding enzyme activities involved in the volatile and ethylene synthesis were severely repressed due to the 1-MCP treatment. The 1-MCP treatment effectively delayed Fenjiao fruit ripening, but affected volatile production by reducing the precursor production and expression level of genes involved in the metabolism pathways of ethylene, auxin and volatiles.

A transcriptomic profile of topping responsive non-coding RNAs in tobacco roots (Nicotiana tabacum)

BACKGROUND

Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long ncRNAs (lncRNAs) and circular RNAs (circRNAs), accomplish remarkable variety of biological functions. However, the composition of ncRNAs and their interactions with coding RNAs in modulating and controlling of cellular process in plants is largely unknown. Using a diverse group of high-throughput sequencing strategies, the mRNA, miRNA, lncRNA and circRNA compositions of tobacco (Nicotiana tabacum) roots determined and their alteration and potential biological functions in response to topping treatment analyzed.

RESULTS

A total of 688 miRNAs, 7423 non-redundant lncRNAs and 12,414 circRNAs were identified, among which, some selected differentially expressed RNAs were verified by quantitative real-time PCR. Using the differentially expressed RNAs, a co-expression network was established that included all four types of RNAs. The number of circRNAs identified were higher than that of miRNAs and lncRNAs, but only two circRNAs were present in the co-expression network. LncRNAs appear to be the most active ncRNAs based on their numbers presented in the co-expression network, but none of them seems to be an eTM (endogenous Target Mimicry) of miRNAs. Integrated with analyses of sequence interaction, several mRNA-circRNA-miRNA interaction networks with a potential role in the regulation of nicotine biosynthesis were uncovered, including a QS-circQS-miR6024 interaction network. In this network miR6024 was significantly down-regulated, while the expression levels of its two targets, circQS and its host gene QS, were sharply increased following the topping treatment.

CONCLUSIONS

These results illustrated the transcriptomic profiles of tobacco roots, the organ responsible for nicotine biosynthesis. mRNAs always play the most important roles, while ncRNAs are also expressed extensively for topping treatment response, especially circRNAs are the most activated in the ncRNA pool. These studies also provided insights on the coordinated regulation module of coding and non-coding RNAs in a single plant biological sample. The findings reported here indicate that ncRNAs appear to form interaction complex for the regulation of stress response forming regulation networks with transcripts involved in nicotine biosynthesis in tobacco.

Oral administration of oat beta-glucan preparations of different molecular weight results in regulation of genes connected with immune response in peripheral blood of rats with LPS-induced enteritis

Purpose

Beta-glucans are biologically active polysaccharides having antioxidant, immunomodulatory, and antiinflammatory properties. This study investigated the transcriptomic profile in peripheral blood of rats with LPS-induced enteritis, which were fed a diet supplemented with high- (G1) and low- (G2) molecular-weight oat beta-glucans.

Methods

Two-color rat gene expression microarrays were applied and the analysis was performed using a common reference design to provide easy means of comparing samples from various experimental conditions against one another. Common reference sample was labeled with cyanine 3 (Cy3) and investigated samples from each experimental group: C-G0 (control group fed semi-synthetic diet), LPS-G0 (LPS-challenged group fed semi-synthetic diet), LPS-G1 (LPS-challenged group fed G1 beta-glucan enriched diet), and LPS-G2 (LPS-challenged group fed G2 beta-glucan enriched diet) were labeled with cyanine 5 (Cy5). Each microarray was performed in quadruplicate. Statistical analysis was performed using one-way ANOVA and Tukey’s HSD post-hoc test (p < 0.05). A multiple testing correction was performed using Benjamini and Hochberg False Discovery Rate < 5%. A quantitative real-time RT-PCR was performed to verify the expression of chosen transcripts.

Results

The microarray analyses revealed differentially expressed transcripts between: the LPS-G0 and the control groups: C-G0 (138 genes), the LPS-G1 and LPS-G0 groups (533 genes), and the LPS-G2 and LPS-G0 groups (97 genes). Several differentially expressed genes in the beta-glucan-supplemented groups encoded proteins belonging to TLR and NLR signaling pathways, as well as prostaglandin synthesis and regulation pathways. Both beta-glucans up-regulated the expression of Atg10, which belongs to the family of autophagy-related genes, suggesting a possible link between autophagy induction and beta-glucan supplementation.

Conclusion

The changes in gene expression observed in the peripheral blood indicate that oat beta-glucans exerted a protective effect in rats with an induced inflammatory state caused by LPS challenge. The greater number of differentially expressed genes was observed in group supplemented with G1 beta-glucan, pointing at the differences in the mode of action of high- and low-molecular-weight beta-glucans in the organism.

Electronic supplementary material

The online version of this article (10.1007/s00394-018-1838-3) contains supplementary material, which is available to authorized users.

Micropterus salmoides rhabdovirus (MSRV) infection induced apoptosis and activated interferon signaling pathway in largemouth bass skin cells

Largemouth bass (Micropterus salmoides) rhabdovirus (MSRV) was isolated from infected juveniles of largemouth bass, and the infected fish exhibited corkscrew, irregular swimming, and crooked body. To our knowledge, the potential molecular mechanisms underlying the pathogenesis of MSRV infection remain largely unknown. In the current study, we found that MSRV infection in largemouth bass skin (LBS) cells induced typical apoptosis, evidenced by the presence of apoptotic bodies and caspase-3 activation. To further analyze the host factors involved in MSRV infection in LBS cells, the transcriptomic profiles during MSRV infection were uncovered using deep RNA sequencing technique, and several differentially expressed genes (DEGs) were validated by quantitative PCR. Our results showed that a total of 124483 unigenes were assembled. Among them, 34465 and 27273 had significant hits to those in the NR and SwissProt databases. After MSRV infection, a total of 2432 and 2480 genes which involved in multiples pathways including TNF signaling, NF-κB signaling, Toll-like receptor signaling and RIG-I signaling pathway were differentially expressed in MSRV infected LBS cells compared to mock-infected cells at 12 h, respectively. Furthermore, quantitative PCR showed that the expression levels of 9 differentially expressed genes (DEGs) related to apoptosis and interferon signaling pathway was consistent with that from transcriptomic profiles. Together, our results not only demonstrated that interferon signaling pathway and apoptosis pathway might exerted crucial roles during MSRV infection, but also provided a useful resource for subsequent investigation of other immune-related genes related to virus infection.

Transcriptome analysis of response to Plasmodiophora brassicae infection in the Arabidopsis shoot and root

BACKGROUND

Clubroot is an important disease caused by the obligate parasite Plasmodiophora brassicae that infects the Brassicaceae. As a soil-borne pathogen, P. brassicae induces the generation of abnormal tissue in the root, resulting in the formation of galls. Root infection negatively affects the uptake of water and nutrients in host plants, severely reducing their growth and productivity. Many studies have emphasized the molecular and physiological effects of the clubroot disease on root tissues. The aim of the present study is to better understand the effect of P. brassicae on the transcriptome of both shoot and root tissues of Arabidopsis thaliana.

RESULTS

Transcriptome profiling using RNA-seq was performed on both shoot and root tissues at 17, 20 and 24 days post inoculation (dpi) of A. thaliana, a model plant host for P. brassicae. The number of differentially expressed genes (DEGs) between infected and uninfected samples was larger in shoot than in root. In both shoot and root, more genes were differentially regulated at 24 dpi than the two earlier time points. Genes that were highly regulated in response to infection in both shoot and root primarily were involved in the metabolism of cell wall compounds, lipids, and shikimate pathway metabolites. Among hormone-related pathways, several jasmonic acid biosynthesis genes were upregulated in both shoot and root tissue. Genes encoding enzymes involved in cell wall modification, biosynthesis of sucrose and starch, and several classes of transcription factors were generally differently regulated in shoot and root.

CONCLUSIONS

These results highlight the similarities and differences in the transcriptomic response of above- and below-ground tissues of the model host Arabidopsis following P. brassicae infection. The main transcriptomic changes in root metabolism during clubroot disease progression were identified. An overview of DEGs in the shoot underlined the physiological changes in above-ground tissues following pathogen establishment and disease progression. This study provides insights into host tissue-specific molecular responses to clubroot development and may have applications in the development of clubroot markers for more effective breeding strategies.

RNA-Sequencing for profiling goat milk transcriptome in colostrum and mature milk

BACKGROUND

In this work we aimed at sequencing and assembling the goat milk transcriptome corresponding at colostrum and 120 days of lactation. To reconstruct transcripts we used both the genome as reference, and a de novo assembly approach. Additionally, we aimed at identifying the differentially expressed genes (DEGs) between the two lactation stages and at analyzing the expression of genes involved in oligosaccharides metabolism.

RESULTS

A total of 44,635 different transcripts, organized in 33,757 tentative genes, were obtained using the goat genome as reference. A significant sequence similarity match was found for 40,353 transcripts (90%) against the NCBI NT and for 35,701 (80%) against the NR databases. 68% and 69% of the de novo assembled transcripts, in colostrum and 120 days of lactation samples respectively, have a significant match with the merged transcriptome obtained using Cufflinks/Cuffmerge. CSN2, PAEP, CSN1S2, CSN3, LALBA, TPT1, FTH1, M-SAA3, SPP1, GLYCAM1, EEF1A1, CTSD, FASN, RPS29, CSN1S1, KRT19 and CHEK1 were found between the top fifteen highly expressed genes. 418 loci were differentially expressed between lactation stages, among which 207 and 122 were significantly up- and down-regulated in colostrum, respectively. Functional annotation and pathway enrichment analysis showed that in goat colostrum somatic cells predominate biological processes involved in glycolysis, carbohydrate metabolism, defense response, cytokine activity, regulation of cell proliferation and cell death, vasculature development, while in mature milk, biological process associated with positive regulation of lymphocyte activation and anatomical structure morphogenesis are enriched. The analysis of 144 different oligosaccharide metabolism-related genes showed that most of these (64%) were more expressed in colostrum than in mature milk, with eight expressed at very high levels (SLCA3, GMSD, NME2, SLC2A1, B4GALT1, B3GNT2, NANS, HEXB).

CONCLUSIONS

To our knowledge, this is the first study comparing goat transcriptome of two lactation stages: colostrum and 120 days. Our findings suggest putative differences of expression between stages and can be envisioned as a base for further research in the topic. Moreover because a higher expression of genes involved in immune defense response, carbohydrate metabolism and related to oligosaccharide metabolism was identified in colostrum we here corroborate the potential of goat milk as a natural source of lactose-derived oligosaccharides and for the development of functional foods.

Characterisation of equine satellite cell transcriptomic profile response to β-hydroxy-β-methylbutyrate (HMB)

β-Hydroxy-β-methylbutyrate (HMB) is a popular ergogenic aid used by human athletes and as a supplement to sport horses, because of its ability to aid muscle recovery, improve performance and body composition. Recent findings suggest that HMB may stimulate satellite cells and affect expressions of genes regulating skeletal muscle cell growth. Despite the scientific data showing benefits of HMB supplementation in horses, no previous study has explained the mechanism of action of HMB in this species. The aim of this study was to reveal the molecular background of HMB action on equine skeletal muscle by investigating the transcriptomic profile changes induced by HMB in equine satellite cells in vitro. Upon isolation from the semitendinosus muscle, equine satellite cells were cultured until the 2nd day of differentiation. Differentiating cells were incubated with HMB for 24 h. Total cellular RNA was isolated, amplified, labelled and hybridised to microarray slides. Microarray data validation was performed with real-time quantitative PCR. HMB induced differential expressions of 361 genes. Functional analysis revealed that the main biological processes influenced by HMB in equine satellite cells were related to muscle organ development, protein metabolism, energy homoeostasis and lipid metabolism. In conclusion, this study demonstrated for the first time that HMB has the potential to influence equine satellite cells by controlling global gene expression. Genes and biological processes targeted by HMB in equine satellite cells may support HMB utility in improving growth and regeneration of equine skeletal muscle; however, the overall role of HMB in horses remains equivocal and requires further proteomic, biochemical and pharmacokinetic studies.

Transcriptome analysis in oak uncovers a strong impact of endogenous rhythmic growth on the interaction with plant-parasitic nematodes

BACKGROUND

Pedunculate oak (Quercus robur L.), an important forest tree in temperate ecosystems, displays an endogenous rhythmic growth pattern, characterized by alternating shoot and root growth flushes paralleled by oscillations in carbon allocation to below- and aboveground tissues. However, these common plant traits so far have largely been neglected as a determining factor for the outcome of plant biotic interactions. This study investigates the response of oak to migratory root-parasitic nematodes in relation to rhythmic growth, and how this plant-nematode interaction is modulated by an ectomycorrhizal symbiont. Oaks roots were inoculated with the nematode Pratylenchus penetrans solely and in combination with the fungus Piloderma croceum, and the systemic impact on oak plants was assessed by RNA transcriptomic profiles in leaves.

RESULTS

The response of oaks to the plant-parasitic nematode was strongest during shoot flush, with a 16-fold increase in the number of differentially expressed genes as compared to root flush. Multi-layered defence mechanisms were induced at shoot flush, comprising upregulation of reactive oxygen species formation, hormone signalling (e.g. jasmonic acid synthesis), and proteins involved in the shikimate pathway. In contrast during root flush production of glycerolipids involved in signalling cascades was repressed, suggesting that P. penetrans actively suppressed host defence. With the presence of the mycorrhizal symbiont, the gene expression pattern was vice versa with a distinctly stronger effect of P. penetrans at root flush, including attenuated defence, cell and carbon metabolism, likely a response to the enhanced carbon sink strength in roots induced by the presence of both, nematode and fungus. Meanwhile at shoot flush, when nutrients are retained in aboveground tissue, oak defence reactions, such as altered photosynthesis and sugar pathways, diminished.

CONCLUSIONS

The results highlight that gene response patterns of plants to biotic interactions, both negative (i.e. plant-parasitic nematodes) and beneficial (i.e. mycorrhiza), are largely modulated by endogenous rhythmic growth, and that such plant traits should be considered as an important driver of these relationships in future studies.

Astaxanthin supplementation enhances adult hippocampal neurogenesis and spatial memory in mice

SCOPE

There is a growing necessity for efficacious natural supplements with antioxidant effects on the brain, in particular, hippocampal function. One such compound, which also has a neuroprotective effect, is the carotenoid astaxanthin (ASX). Despite ASX's potential benefit to the brain, very little is known about its effect on hippocampal plasticity and cognition. Thus, we investigated the effect of ASX on adult hippocampal neurogenesis (AHN) and spatial memory using a mouse model.

METHODS AND RESULTS

Dose-response was examined in mice fed ASX-supplemented diets (0, 0.02, 0.1, and 0.5%) to define the effect of ASX on AHN. In conjunction with AHN results, hippocampus-dependent cognitive function was assessed. We delineated molecular mechanisms associated with ASX-enhanced AHN using DNA microarray analysis. Results revealed that ASX enhanced cell proliferation and survival at 0.1% and 0.5% doses. Newborn mature neurons were higher only with 0.5% ASX, which also enhanced spatial memory. Transcriptomic profiling revealed potential AHN-associated molecules (Prl, Itga4, and Il4) that were ASX induced. Their downstream factors, identified through Ingenuity Pathway Analysis, were positively correlated with ASX-induced increases in spatial memory.

CONCLUSION

ASX supplementation enhanced AHN and spatial memory, and a DNA microarray approach provided, for the first time, novel molecular insights into ASX action.

Transcriptomic analysis of the porcine endometrium during early pregnancy and the estrous cycle

The goal of this study was to describe the alterations in the transcriptome of the endometrium in pigs during the beginning of implantation (days 15-16 of pregnancy) compared to cyclic pigs during the onset of luteolysis (days 15-16 of the estrous cycle). The global expression of genes in porcine gravid and non-gravid endometria was investigated using the Porcine (V2) Two-color gene expression microarray, 4 × 44 (Agilent, USA). Analysis of the microarray data showed that, of 589 accurately annotated genes, the expression of 266 genes was up-regulated and expression of 323 was down-regulated in the endometrium harvested during early pregnancy compared with the endometrium during the estrous cycle. In pregnant pigs, genes with the most significantly altered expression were involved in the following biological processes: the metabolic process, cellular process, cell communication, immune system process, developmental process, cell adhesion, antigen processing and presentation, antigen processing and presentation of peptide or polysaccharide antigen via major histocompatibility complex (MHC) class II, immune response, and the polysaccharide metabolic process. In the pregnant endometrium, cell adhesion molecules and steroid hormone biosynthesis pathways were the most significantly enriched biological pathways. Analysis of the interaction network among selected genes showed that androgen receptor (AR) encoding genes interact with genes involved in important processes occurring during early pregnancy. The bioinformatic analysis revealed information about the meaning of differentially expressed genes. The data provided new insight into the dynamic changes of the endometrial gene expression profile during days 15-16 of pregnancy.