RNA-Seq Data Analysis: From Raw Data Quality Control to Differential Expression Analysis

Transcriptomics analysis of long non-coding RNAs in smooth muscle cells from patients with peripheral artery disease and diabetes mellitus

Diabetes mellitus (DM) is a significant risk factor for peripheral arterial disease (PAD), and PAD is an independent predictor of cardiovascular disorders (CVDs). Growing evidence suggests that long non-coding RNAs (lncRNAs) significantly contribute to disease development and underlying complications, particularly affecting smooth muscle cells (SMCs). So far, no study has focused on transcriptome analysis of lncRNAs in PAD patients with and without DM. Tissue samples were obtained from our Vascular Biobank. Due to the sample's heterogeneity, expression analysis of lncRNAs in whole tissue detected only ACTA2-AS1 with a 4.9-fold increase in PAD patients with DM. In contrast, transcriptomics of SMCs revealed 28 lncRNAs significantly differentially expressed between PAD with and without DM (FDR < 0.1). Sixteen lncRNAs were of unknown function, six were described in cancer, one connected with macrophages polarisation, and four were associated with CVDs, mainly with SMC function and phenotypic switch (NEAT1, MIR100HG, HIF1A-AS3, and MRI29B2CHG). The enrichment analysis detected additional lncRNAs H19, CARMN, FTX, and MEG3 linked with DM. Our study revealed several lncRNAs in diabetic PAD patients associated with the physiological function of SMCs. These lncRNAs might serve as potential therapeutic targets to improve the function of SMCs within the diseased tissue and, thus, the clinical outcome.

Effect of high-dose glucocorticoid treatment on human brown adipose tissue activity: a randomised, double-blinded, placebo-controlled cross-over trial in healthy men

BACKGROUND

Glucocorticoids (GCs) are widely applied anti-inflammatory drugs that are associated with adverse metabolic effects including insulin resistance and weight gain. Previous research indicates that GCs may negatively impact brown adipose tissue (BAT) activity in rodents and humans.

METHODS

We performed a randomised, double-blinded cross-over trial in 16 healthy men (clinicaltrials.govNCT03269747). Participants received 40 mg of prednisone per day for one week or placebo. After a washout period of four weeks, participants crossed-over to the other treatment arm. Primary endpoint was the increase in resting energy expenditure (EE) in response to a mild-cold stimulus (cold-induced thermogenesis, CIT). Secondary outcomes comprised mean 18F-FDG uptake into supraclavicular BAT (SUVmean) as determined by FDG-PET/CT, volume of the BAT depot as well as fat content determined by MRI. The plasma metabolome and the transcriptome of supraclavicular BAT and of skeletal muscle biopsies after each treatment period were analysed.

FINDINGS

Sixteen participants were recruited to the trial and completed it successfully per protocol. After prednisone treatment resting EE was higher both during warm and cold conditions. However, CIT was similar, 153 kcal/24 h (95% CI 40-266 kcal/24 h) after placebo and 186 kcal/24 h (95% CI 94-277 kcal/24 h, p = 0.38) after prednisone. SUVmean of BAT after cold exposure was not significantly affected by prednisone (3.36 g/ml, 95% CI 2.69-4.02 g/ml, vs 3.07 g/ml, 95% CI 2.52-3.62 g/ml, p = 0.28). Results of plasma metabolomics and BAT transcriptomics corroborated these findings. RNA sequencing of muscle biopsies revealed higher expression of genes involved in calcium cycling. No serious adverse events were reported and adverse events were evenly distributed between the two treatments.

INTERPRETATION

Prednisone increased EE in healthy men possibly by altering skeletal muscle calcium cycling. Cold-induced BAT activity was not affected by GC treatment, which indicates that the unfavourable metabolic effects of GCs are independent from thermogenic adipocytes.

FUNDING

Grants from Swiss National Science Foundation (PZ00P3_167823), Bangerter-Rhyner Foundation and from Nora van der Meeuwen-Häfliger Foundation to MJB. A fellowship-grant from the Swiss National Science Foundation (SNF211053) to WS. Grants from German Research Foundation (project number: 314061271-TRR 205) and Else Kröner-Fresenius (grant support 2012_A103 and 2015_A228) to MR.

The Genomic Variation and Differentially Expressed Genes on the 6P Chromosomes in Wheat-Agropyron cristatum Addition Lines 5113 and II-30-5 Confer Different Desirable Traits

Wild relatives of wheat are essential gene pools for broadening the genetic basis of wheat. Chromosome rearrangements and genomic variation in alien chromosomes are widespread. Knowledge of the genetic variation between alien homologous chromosomes is valuable for discovering and utilizing alien genes. In this study, we found that 5113 and II-30-5, two wheat-A. cristatum 6P addition lines, exhibited considerable differences in heading date, grain number per spike, and grain weight. Genome resequencing and transcriptome analysis revealed significant differences in the 6P chromosomes of the two addition lines, including 143,511 single-nucleotide polymorphisms, 62,103 insertion/deletion polymorphisms, and 757 differentially expressed genes. Intriguingly, genomic variations were mainly distributed in the middle of the chromosome arms and the proximal centromere region. GO and KEGG analyses of the variant genes and differentially expressed genes showed the enrichment of genes involved in the circadian rhythm, carbon metabolism, carbon fixation, and lipid metabolism, suggesting that the differential genes on the 6P chromosome are closely related to the phenotypic differences. For example, the photosynthesis-related genes PsbA, PsbT, and YCF48 were upregulated in II-30-5 compared with 5113. ACS and FabG are related to carbon fixation and fatty acid biosynthesis, respectively, and both carried modification variations and were upregulated in 5113 relative to II-30-5. Therefore, this study provides important guidance for cloning desirable genes from alien homologous chromosomes and for their effective utilization in wheat improvement.

Transcriptomic profiling of canine decidualization and effects of antigestagens on decidualized dog uterine stromal cells

Maternal-stroma derived decidual cells, the only cell population in the canine placenta expressing the nuclear progesterone (P4) receptor (PGR), are crucial for the maintenance of canine pregnancy. Decreased circulating progesterone (P4) levels, or blockage of PGR function with antigestagens, terminate canine pregnancy. As an in vitro model for canine decidualization, dog uterine stromal (DUS) cells can be decidualized in vitro with cAMP. The antigestagens aglepristone and mifepristone ablate the expression of decidualization markers in DUS cells (e.g., PGR, PRLR, IGF1 or PTGES). Here, the transcriptome profile of DUS cells was investigated to acquire deeper insights into decidualization-associated changes. Additionally, effects mediated by antigestagens (competitive PGR blockers) in decidualized cells were assessed. Decidualization led to the upregulation of 1841 differentially expressed genes (DEGs, P and FDR < 0.01) involved in cellular proliferation and adhesion, mesenchymal-epithelial transition, extracellular matrix organization, and vaso- and immunomodulation. The 1475 DEGs downregulated after decidualization were mostly associated with apoptosis and cell migration. In decidualized DUS cells, aglepristone modulated 1400 DEGs and mifepristone 1558 DEGs. Interestingly, around half of the identified DEGs were modulated by only one of the antigestagens. In all cases, however, PGR-blockage was mainly associated with an inversion of several decidualization-induced effects. Comparison between antigestagen-mediated effects and transcriptional changes in the canine placenta at term allowed the identification of 191 DEGs associated with diminished cell proliferation and adhesion, and vascular and immune modulation. This study emphasizes the importance of P4/PGR signaling for decidual cell function, providing new insights into the maintenance of canine pregnancy.

Effects of ACTH-Induced Long-Term Hypercortisolism on the Transcriptome of Canine Visceral Adipose Tissue

Cushing’s syndrome, or hypercortisolism (HC), a common endocrinopathy in adult dogs, is caused by chronic hypercortisolemia. Among different metabolic disorders, this syndrome is associated with enhanced subcutaneous lipolysis and visceral adiposity. However, effects of HC in adipose tissue, especially regarding visceral adipose tissue (VAT), are still poorly understood. Herein, the transcriptomic effects of chronic HC on VAT of dogs were evaluated. For this, subcutaneously implanted ACTH-releasing pumps were used, followed by deep RNA sequencing of the canine VAT. Prolonged HC seems to affect a plethora of regulatory mechanisms in VAT of treated dogs, with 1190 differentially expressed genes (DEGs, p and FDR < 0.01) being found. The 691 downregulated DEGs were mostly associated with functional terms like cell adhesion and migration, intracellular signaling, immune response, extracellular matrix and angiogenesis. Treatment also appeared to modulate local glucocorticoid and insulin signaling and hormonal sensitivity, and several factors, e.g., TIMP4, FGF1, CCR2, CXCR4 and HSD11B1/2, were identified as possible important players in the glucocorticoid-related expansion of VAT. Modulation of their function during chronic HC might present interesting targets for further clinical studies. Similarities in the effects of chronic HC on VAT of dogs and humans are highlighted.

Identification of Alternative Splicing Events Associated with Paratuberculosis in Dairy Cattle Using Multi-Tissue RNA Sequencing Data

Paratuberculosis is a major endemic disease caused by Mycobacterium avium subspecies paratuberculosis (MAP) infection and leads to huge economic loss in the dairy sector worldwide. Alternative splicing (AS) events, playing indispensable regulatory roles in many protein functions and biological pathways, are shown to be associated with complex traits and diseases. In this study, by integrating the RNA sequencing (RNA-seq) data of 24 samples from three tissues (peripheral blood, jejunum and salivary gland) of Holstein cows, we obtained 2,706,541,696 uniquely mapped reads in total that represented 12,870 expressed genes, and detected 4285 differentially expressed genes (DEGs) between MAP-infected and healthy cows (p < 0.05). Of them, 92 differentially expressed splicing factors (DESFs) were included. Further, 119, 150 and 68 differential alternative splicing (DAS) events between MAP-infected and healthy cows were identified in peripheral blood, jejunum and salivary glands, respectively. Of note, six DAS events were highly and significantly correlated with the DESFs (R2 > 0.9; p < 0.01), and their corresponding genes (COPI coat complex subunit gamma 2gene (COPG2), kinesin family member 2C gene (KIF2C), exocyst complex component 7 (EXOC7), Rab9 effector protein with kelch motifs gene (RABEPK), deoxyribonuclease 1 gene (DNASE1) and early endosome antigen 1gene (EEA1)) were significantly enriched in immune response such as vesicle-mediated transport, regulation of acute inflammatory response and tuberculosis through gene ontology (GO) and KEGG analysis. KS test showed that the DAS events in the EXOC7 and KIF2C genes indeed displayed differences between MAP-infected cows and healthy cows. The DAS in EXOC7 might produce a new protein sequence with lack of 23 amino acids, and the DAS in KIF2C induced a stop codon of premature occurrence and resulted in a lack of functional domain. In summary, this study identified the DAS events and corresponding genes related to MAP-infection base on the RNA-seq data from multiple tissues of Holstein cows, providing novel insights into the regulatory mechanisms underpinning paratuberculosis in dairy cattle.

Trends and Applications of Omics Technologies to Functional Characterisation of Enzymes and Protein Metabolites Produced by Fungi

Identifying and adopting industrial applications for proteins and enzymes derived from fungi strains have been at the focal point of several studies in recent times. To facilitate such studies, it is necessary that advancements and innovation in mycological and molecular characterisation are concomitant. This review aims to provide a detailed overview of the necessary steps employed in both qualitative and quantitative research using the omics technologies that are pertinent to fungi characterisation. This stems from the understanding that data provided from the functional characterisation of fungi and their metabolites is important towards the techno-economic feasibility of large-scale production of biological products. The review further describes how the functional gaps left by genomics, internal transcribe spacer (ITS) regions are addressed by transcriptomics and the various techniques and platforms utilised, including quantitive reverse transcription polymerase chain reaction (RT-qPCR), hybridisation techniques, and RNA-seq, and the insights such data provide on the effect of environmental changes on fungal enzyme production from an expressional standpoint. The review also offers information on the many available bioinformatics tools of analysis necessary for the analysis of the overwhelming data synonymous with the omics approach to fungal characterisation.

Transcriptomic Analysis Exploring the Molecular Mechanisms of Hanchuan Zupa Granules in Alleviating Asthma in Rat

Objective

To investigate the molecular mechanisms of HCZP treatment of asthma.

Materials and Methods

Thirty Sprague Dawley (SD) rats were divided into normal, asthma, and HCZP groups (n = 10). The asthma model was sensitized by 1 mg ovalbumin (OVA)/aluminum hydroxide Al(OH)₃mixture and then challenged with 1% aerosolized OVA for four weeks. Rats in the HCZP group received 10.08 g/kg/d HCZP for four weeks during OVA challenge. Then, lung tissues of rats in each group were collected for RNA sequencing. Moreover, the expression level of some core genes was detected by using western blotting and immunohistochemistry.

Results

Inflammatory cell infiltration and pathological damage of the lungs improved in the HCZP group. Compared with the asthma group (0.049 ± 0.002 mm²/mm; 0.036 ± 0.006 mm²/mm; and 0.014 ± 0.001 mm²/mm), total wall thickness (0.042 ± 0.001 mm²/mm), inner wall thickness (0.013 ± 0.001 mm²/mm), and smooth muscle layer thickness (0.012 ± 0.001 mm²/mm) significantly decreased in the HCZP group. Bioinformatics analysis showed that hub genes such as bradykinin receptor B2 (Bdkrb2) and CD4 molecule (Cd4) had different expression patterns between model and HCZP groups. Two transcription factors, forkhead box Q1 (Foxq1) and nuclear factor of activated T cells 2 (Nfatc2), served important regulatory roles in asthma. Compared with the model group, Bdkrb2 protein expression increased and Nfatc2 protein expression decreased in the HCZP group. Discussion and Conclusion. HCZP could alleviate asthma via regulating the expression of several hub genes, which might serve as therapeutic targets for asthma. However, the mechanism of these genes will be studied in the future.

Global transcriptome analysis implicates cholesterol availability in the regulation of canine cyclic luteal function

Considering the key role of the corpus luteum in the regulation of the canine diestrus, the present study aimed to investigate changes in the luteal transcriptome of pseudopregnant dogs (n = 18) from days (D) 10, 20, 30, 40, 50 and 60 post-ovulation. After RNAsequencing was performed, data was analyzed by resorting to several informatic tools. A total of 3300 genes were differently expressed among all samples (FDR < 0.01). By comparing different time points, enriched biological processes as response to estradiol and lipids (D20 vs D10) and intracellular cholesterol transport (D40 vs D60) were observed. Moreover, LXR/RXR (liver X receptor- retinoid X receptor) signaling appeared as an overrepresented pathway in all comparisons. Thus, the expression of 19 genes involved in intracellular cholesterol availability was further evaluated; most were affected by time (P < 0.05). Adding to the deep transcriptomic analysis, presented data implies the importance of cholesterol regulation in luteal physiology of pseudopregnant dogs.

Curcumin suppresses LGR5(+) colorectal cancer stem cells by inducing autophagy and via repressing TFAP2A-mediated ECM pathway

Abstract

Colorectal cancer stem cells (CSCs) have the potential for self-renewal, proliferation, and differentiation. And LGR5 is a stem cell marker gene of colorectal cancer. Curcumin can suppress oncogenicity of many cancer cells, yet the effect and mechanism of curcumin in LGR5(+) colorectal cancer stem cells (CSCs) have not been studied. In this study, we studied the effect of curcumin on LGR5(+) colorectal CSCs using the experiments of tumorsphere formation, cell viability and cell apoptosis. Then autophagy analysis, RNA-Seq, and real-time PCR were used to identify the mechanism responsible for the inhibition of LGR5(+) colorectal CSCs. Our results showed that curcumin inhibited tumorsphere formation, decreased cell viability in a dose-dependent manner, and also promoted apoptosis of LGR5(+) colorectal CSCs. Next, we found curcumin induced autophagy of LGR5(+) colorectal CSCs. When LGR5(+) colorectal CSCs were co-treated with curcumin and the autophagy inhibitor (hydroxychloroquine), curcumin-induced cell proliferation inhibition decreased. In addition, we also found that curcumin inhibited the extracellular matrix (ECM)-receptor interaction pathway via the downregulation of the following genes: GP1BB, COL9A3, COMP, AGRN, ITGB4, LAMA5, COL2A1, ITGB6, ITGA1, and TNC. Further, these genes were transcriptionally regulated by TFAP2A, and the high expression of TFAP2A was associated with poor prognosis in colorectal cancer. In conclusion, curcumin suppressed LGR5(+) colorectal CSCs, potentially by inducing autophagy and repressing the oncogenic TFAP2A-mediated ECM pathway.

Graphic abstract

Downregulation of DUSP9 Promotes Tumor Progression and Contributes to Poor Prognosis in Human Colorectal Cancer

Background

Dual-specificity phosphatase 9 (DUSP9) belongs to the dual-specificity protein phosphatase subfamily. Recently, increasing attention has been paid on the role of DUSP9 in a variety of cancers. However, its functional role in tumor development is still unclear, especially in colorectal cancer (CRC).

Methods

The functional role of DUSP9 in inhibiting the progression of CRC was verified using colony formation assay, wound healing assay, nude mice xenograft model, etc. RNA-seq was performed to assess the gene expression profiling in SW480 cells with DUSP9 stable knockdown and shControl cells. Bisulfite sequencing (BSE) was performed to reveal the methylation status of CpG island in the promoter of DUSP9.

Results

DUSP9 was significantly downregulated in tumor tissues compared with peritumor tissues. Mechanistically, the high methylation status of CpG island in the promoter of DUSP9 may lead to the downregulation of DUSP9 in CRC. Clinically, low DUSP9 expression in CRC was closely associated with depth of invasion, metastasis (TNM) stage, and poor survival, indicating that DUSP9 may be involved in the progression of CRC. Functional study revealed that DUSP9 inhibited proliferation, migration, invasion, and epithelial–mesenchymal transition of CRC cells both in vitro and in vivo. Transcriptome profiling studies revealed that Erk signaling was involved in the tumor progression mediated by DUSP9 silencing, which is confirmed by cell experiments and clinical tissue sample staining analysis.

Conclusion

Our findings demonstrate that DUSP9 plays a critical role in the progression of CRC, and therapeutic intervention to increase the expression or activity of DUSP9 may be a potential target for CRC treatment in the future.

Gene expression profiling of the canine placenta during normal and antigestagen-induced luteolysis

The domestic dog is the only domestic animal species that does not produce steroids in the placenta and instead relies on luteal steroids throughout pregnancy. Nevertheless, the canine placenta is highly responsive to steroids, and withdrawal of progesterone (P4) affects the feto-maternal unit, initializing the parturition cascade. Similar effects can be observed during antigestagen-induced abortion. Here, aiming to provide new insights into mechanisms involved in the termination of canine pregnancy, next generation sequencing (NGS, RNA-seq) was applied. Placental transcriptomes derived from natural prepartum and antigestagen-induced abortions were analyzed and compared with fully developed mid-gestation placentas. The contrast "prepartum luteolysis over mid-gestation" revealed 1973 differentially expressed genes (DEG). Terms associated with apoptosis, impairment of vascular function and activation of signaling of several cytokines (e.g., IL-8, IL-3, TGF-β) were overrepresented at natural luteolysis. When compared with mid-term, antigestagen treatment revealed 135 highly regulated DEG that were involved in the induced luteolysis and showed similar associations with functional terms and expression patterns as during natural luteolysis. The contrast "antigestagen-induced luteolysis over prepartum luteolysis" revealed that, although similar changes occur in both conditions, they are more pronounced during natural prepartum. Among P4-regulated DEG were those related to immune system and cortisol metabolism. It appears that, besides inducing placental PGF2α output, both natural and induced P4 withdrawal is associated with disruption of the feto-maternal interface, leading to impaired vascular functions, apoptosis and controlled modulation of the immune response. The time-related maturation of the feto-maternal interface needs to be considered because it may be clinically relevant.

Global Transcriptomic Analysis of the Canine corpus luteum (CL) During the First Half of Diestrus and Changes Induced by in vivo Inhibition of Prostaglandin Synthase 2 (PTGS2/COX2)

The canine luteal phase exhibits several peculiarities compared with other species. In early diestrus, the corpus luteum (CL) is, at least in part, independent of gonadotropins, and prostaglandins (PGs) appear to be among its main regulators. This was also observed with the inhibition in vivo of COX2, when also transcriptional capacity, vascularization and immune-related factors were affected. Here, we aimed to further investigate the potential effects of PGs withdrawal on the CL transcriptome by performing deep RNA sequencing (RNA-Seq). Samples from a previous in vivo study were used; bitches were treated for 5, 10, 20, or 30 days after ovulation with firocoxib (Previcox®), a PTGS2/COX2 inhibitor, or a placebo. Analysis of results was performed with SUSHI (framework from FGCZ) and with pathways and functional networks analyzers. Time-dependent effects were also investigated and used for quality control. More highly represented differentially expressed genes (DEGs, P < 0.01, FDR < 0.1) in the early CL (days 5 and 10) referred to proliferation and immune system, while in the mature CL (days 20 and 30) they were related with steroidogenesis. The absence of genes concomitantly affected by the treatment at all time-points suggested stage-dependency in the observed effects. Little effect was observed on days 5 and 10. Day 20 had the highest number of DEGs (n = 1,741), related with increased immune response. On day 30, DEGs found (n = 552) referred to decreased steroidogenesis and vascularization. Our results suggest the presence of strong compensatory effects in the early CL and multidirectional effects toward gonadotropin-dependency of the CL after COX2 inhibition.