Comprehensive molecular and clinicopathological profiling of desmoid tumours

Previous studies have not clearly identified a prognostic factor for desmoid tumours (DT). Whole-exome sequencing (WES) and/or RNA sequencing (RNA-seq) were performed in 64 cases of DT to investigate the molecular profiles in combination with the clinicopathological characteristics. CTNNB1 mutations with specific hotspots were identified in 56 cases (87.5%). A copy number loss in chromosome 6 (chr6) was identified in 14 cases (21.9%). Clustering based on the mRNA expression profiles was predictive of the patients' prognoses. The risk score generated by the expression of a three-gene set (IFI6, LGMN, and CKLF) was a strong prognostic marker for recurrence-free survival (RFS) in our cohort. In risk groups stratified by the expression of IFI6, the hazard ratio for recurrence-free survival in the high-risk group relative to the low-risk group was 12.12 (95% confidence interval: 1.56-94.2; p = 8.0 × 10⁶). In conclusion, CTNNB1 mutations and a chr6 copy number loss are likely the causative mutations underlying the tumorigenesis of DT while the gene expression profiles may help to differentiate patients who would be good candidates for wait-and-see management and those who might benefit from additional systemic or radiation therapies.

Overexpression of the Bacteriophage T4 motB Gene Alters H-NS Dependent Repression of Specific Host DNA

The bacteriophage T4 early gene product MotB binds tightly but nonspecifically to DNA, copurifies with the host Nucleoid Associated Protein (NAP) H-NS in the presence of DNA and improves T4 fitness. However, the T4 transcriptome is not significantly affected by a motB knockdown. Here we have investigated the phylogeny of MotB and its predicted domains, how MotB and H-NS together interact with DNA, and how heterologous overexpression of motB impacts host gene expression. We find that motB is highly conserved among Tevenvirinae. Although the MotB sequence has no homology to proteins of known function, predicted structure homology searches suggest that MotB is composed of an N-terminal Kyprides-Onzonis-Woese (KOW) motif and a C-terminal DNA-binding domain of oligonucleotide/oligosaccharide (OB)-fold; either of which could provide MotB’s ability to bind DNA. DNase I footprinting demonstrates that MotB dramatically alters the interaction of H-NS with DNA in vitro. RNA-seq analyses indicate that expression of plasmid-borne motB up-regulates 75 host genes; no host genes are down-regulated. Approximately 1/3 of the up-regulated genes have previously been shown to be part of the H-NS regulon. Our results indicate that MotB provides a conserved function for Tevenvirinae and suggest a model in which MotB functions to alter the host transcriptome, possibly by changing the association of H-NS with the host DNA, which then leads to conditions that are more favorable for infection.

Deconstructing transcriptional variations and their effects on immunomodulatory function among human mesenchymal stromal cells

Background

Mesenchymal stromal cell (MSC)-based therapies are being actively investigated in various inflammatory disorders. However, functional variability among MSCs cultured in vitro will lead to distinct therapeutic efficacies. Until now, the mechanisms behind immunomodulatory functional variability in MSCs are still unclear.

Methods

We systemically investigated transcriptomic variations among MSC samples derived from multiple tissues to reveal their effects on immunomodulatory functions of MSCs. We then analyzed transcriptomic changes of MSCs licensed with INFγ to identify potential molecular mechanisms that result in distinct MSC samples with different immunomodulatory potency.

Results

MSCs were clustered into distinct groups showing different functional enrichment according to transcriptomic patterns. Differential expression analysis indicated that different groups of MSCs deploy common regulation networks in response to inflammatory stimulation, while expression variation of genes in the networks could lead to different immunosuppressive capability. These different responsive genes also showed high expression variability among unlicensed MSC samples. Finally, a gene panel was derived from these different responsive genes and was able to regroup unlicensed MSCs with different immunosuppressive potencies.

Conclusion

This study revealed genes with expression variation that contribute to immunomodulatory functional variability of MSCs and provided us a strategy to identify candidate markers for functional variability assessment of MSCs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13287-020-02121-8.

scReQTL: an approach to correlate SNVs to gene expression from individual scRNA-seq datasets

BACKGROUND

Recently, pioneering expression quantitative trait loci (eQTL) studies on single cell RNA sequencing (scRNA-seq) data have revealed new and cell-specific regulatory single nucleotide variants (SNVs). Here, we present an alternative QTL-related approach applicable to transcribed SNV loci from scRNA-seq data: scReQTL. ScReQTL uses Variant Allele Fraction (VAFRNA) at expressed biallelic loci, and corelates it to gene expression from the corresponding cell.

RESULTS

Our approach employs the advantage that, when estimated from multiple cells, VAFRNA can be used to assess effects of SNVs in a single sample or individual. In this setting scReQTL operates in the context of identical genotypes, where it is likely to capture RNA-mediated genetic interactions with cell-specific and transient effects. Applying scReQTL on scRNA-seq data generated on the 10 × Genomics Chromium platform using 26,640 mesenchymal cells derived from adipose tissue obtained from three healthy female donors, we identified 1272 unique scReQTLs. ScReQTLs common between individuals or cell types were consistent in terms of the directionality of the relationship and the effect size. Comparative assessment with eQTLs from bulk sequencing data showed that scReQTL analysis identifies a distinct set of SNV-gene correlations, that are substantially enriched in known gene-gene interactions and significant genome-wide association studies (GWAS) loci.

CONCLUSION

ScReQTL is relevant to the rapidly growing source of scRNA-seq data and can be applied to outline SNVs potentially contributing to cell type-specific and/or dynamic genetic interactions from an individual scRNA-seq dataset.

AVAILABILITY

https://github.com/HorvathLab/NGS/tree/master/scReQTL.

Transcriptomic time-series analysis of cold- and heat-shock response in psychrotrophic lactic acid bacteria

BACKGROUND

Psychrotrophic lactic acid bacteria (LAB) species are the dominant species in the microbiota of cold-stored modified-atmosphere-packaged food products and are the main cause of food spoilage. Despite the importance of psychrotrophic LAB, their response to cold or heat has not been studied. Here, we studied the transcriptome-level cold- and heat-shock response of spoilage lactic acid bacteria with time-series RNA-seq for Le. gelidum, Lc. piscium, and P. oligofermentans at 0 °C, 4 °C, 14 °C, 25 °C, and 28 °C.

RESULTS

We observed that the cold-shock protein A (cspA) gene was the main cold-shock protein gene in all three species. Our results indicated that DEAD-box RNA helicase genes (cshA, cshB) also play a critical role in cold-shock response in psychrotrophic LAB. In addition, several RNase genes were involved in cold-shock response in Lc. piscium and P. oligofermentans. Moreover, gene network inference analysis provided candidate genes involved in cold-shock response. Ribosomal proteins, tRNA modification, rRNA modification, and ABC and efflux MFS transporter genes clustered with cold-shock response genes in all three species, indicating that these genes could be part of the cold-shock response machinery. Heat-shock treatment caused upregulation of Clp protease and chaperone genes in all three species. We identified transcription binding site motifs for heat-shock response genes in Le. gelidum and Lc. piscium. Finally, we showed that food spoilage-related genes were upregulated at cold temperatures.

CONCLUSIONS

The results of this study provide new insights on the cold- and heat-shock response of psychrotrophic LAB. In addition, candidate genes involved in cold- and heat-shock response predicted using gene network inference analysis could be used as targets for future studies.

Gene expression patterns of red sea urchins (Mesocentrotus franciscanus) exposed to different combinations of temperature and pCO2 during early development

BACKGROUND

The red sea urchin Mesocentrotus franciscanus is an ecologically important kelp forest herbivore and an economically valuable wild fishery species. To examine how M. franciscanus responds to its environment on a molecular level, differences in gene expression patterns were observed in embryos raised under combinations of two temperatures (13 °C or 17 °C) and two pCO2 levels (475 μatm or 1050 μatm). These combinations mimic various present-day conditions measured during and between upwelling events in the highly dynamic California Current System with the exception of the 17 °C and 1050 μatm combination, which does not currently occur. However, as ocean warming and acidification continues, warmer temperatures and higher pCO2 conditions are expected to increase in frequency and to occur simultaneously. The transcriptomic responses of the embryos were assessed at two developmental stages (gastrula and prism) in light of previously described plasticity in body size and thermotolerance under these temperature and pCO2 treatments.

RESULTS

Although transcriptomic patterns primarily varied by developmental stage, there were pronounced differences in gene expression as a result of the treatment conditions. Temperature and pCO2 treatments led to the differential expression of genes related to the cellular stress response, transmembrane transport, metabolic processes, and the regulation of gene expression. At each developmental stage, temperature contributed significantly to the observed variance in gene expression, which was also correlated to the phenotypic attributes of the embryos. On the other hand, the transcriptomic response to pCO2 was relatively muted, particularly at the prism stage.

CONCLUSIONS

M. franciscanus exhibited transcriptomic plasticity under different temperatures, indicating their capacity for a molecular-level response that may facilitate red sea urchins facing ocean warming as climate change continues. In contrast, the lack of a robust transcriptomic response, in combination with observations of decreased body size, under elevated pCO2 levels suggest that this species may be negatively affected by ocean acidification. High present-day pCO2 conditions that occur due to coastal upwelling may already be influencing populations of M. franciscanus.

The R-enantiomer of ketorolac reduces ovarian cancer tumor burden in vivo

BACKGROUND

Rho-family GTPases, including Ras-related C3 botulinum toxin substrate 1 (Rac1) and cell division control protein 42 (Cdc42), are important modulators of cancer-relevant cell functions and are viewed as promising therapeutic targets. Based on high-throughput screening and cheminformatics we identified the R-enantiomer of an FDA-approved drug (ketorolac) as an inhibitor of Rac1 and Cdc42. The corresponding S-enantiomer is a non-steroidal anti-inflammatory drug (NSAID) with selective activity against cyclooxygenases. We reported previously that R-ketorolac, but not the S-enantiomer, inhibited Rac1 and Cdc42-dependent downstream signaling, growth factor stimulated actin cytoskeleton rearrangements, cell adhesion, migration and invasion in ovarian cancer cell lines and patient-derived tumor cells.

METHODS

In this study we treated mice with R-ketorolac and measured engraftment of tumor cells to the omentum, tumor burden, and target GTPase activity. In order to gain insights into the actions of R-ketorolac, we also performed global RNA-sequencing (RNA-seq) analysis on tumor samples.

RESULTS

Treatment of mice with R-ketorolac decreased omental engraftment of ovarian tumor cells at 18 h post tumor cell injection and tumor burden after 2 weeks of tumor growth. R-ketorolac treatment inhibited tumor Rac1 and Cdc42 activity with little impact on mRNA or protein expression of these GTPase targets. RNA-seq analysis revealed that R-ketorolac decreased expression of genes in the HIF-1 signaling pathway. R-ketorolac treatment also reduced expression of additional genes associated with poor prognosis in ovarian cancer.

CONCLUSION

These findings suggest that R-ketorolac may represent a novel therapeutic approach for ovarian cancer based on its pharmacologic activity as a Rac1 and Cdc42 inhibitor. R-ketorolac modulates relevant pathways and genes associated with disease progression and worse outcome.

Signature changes in the expressions of protein-coding genes, lncRNAs, and repeat elements in early and late cellular senescence

Replicative cellular senescence is the main cause of aging. It is important to note that early senescence is linked to tissue regeneration, whereas late senescence is known to trigger a chronically inflammatory phenotype. Despite the presence of various genome-wide studies, there is a lack of information on distinguishing early and late senescent phenotypes at the transcriptome level. Particularly, the changes in the noncoding RNA portion of the aging cell have not been fully elucidated. By utilising RNA sequencing data of fibroblasts, hereby, we are not only reporting changes in gene expression profiles and relevant biological processes in the early and late senescent phenotypes but also presenting significant differences in the expressions of many unravelled long noncoding RNAs (lncRNAs) and transcripts arisen from repetitive DNA. Our results indicate that, in addition to previously reported L1 elements, various LTR and DNA transposons, as well as members of the classical satellites including HSAT5 and α-satellites (ALR/Alpha), are expressed at higher levels in late senescence. Moreover, we revealed finer links between the expression levels of repeats with the genes located near them and known to be involved in cell cycle and senescence. Noncoding elements reported here provide a new perspective to be explored in further experimental studies.

Identification of critical pathways and potential therapeutic targets in poorly differentiated duodenal papilla adenocarcinoma

Background

Duodenal papilla carcinoma (DPC) is a rare malignancy of the gastrointestinal tract with high recurrence rate, and the pathogenesis of this highly malignant neoplasm is yet to be fully elucidated. This study aims to identify key genes to further understand the biology and pathogenesis underlying the molecular alterations driving DPC, which could be potential diagnostic or therapeutic targets.

Methods

Tumor samples of three DPC patients were collected and integrating RNA-seq analysis of tumor tissues and matched normal tissues were performed to discover differentially expressed genes (DEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were carried out to understand the potential bio-functions of the DPC differentially expressed genes (DEGs). Protein–protein interaction (PPI) network was constructed for functional modules analysis and identification of hub genes. qRT-PCR of clinical samples was conducted to validate the expression level of the hub genes.

Results

A total of 110 DEGs were identified from our RNA-seq data, GO and KEGG analyses showed that the DEGs were mainly enriched in multiple cancer-related functions and pathways, such as cell proliferation, IL-17signaling pathway, Jak-STAT signaling pathway, PPAR signaling pathway. The PPI network screened out five hub genes including IL-6, LCN2, FABP4, LEP and MMP1, which were identified as core genes in the network and the expression value were validated by qRT-PCR. The hub genes identified in this work were suggested to be potential therapeutic targets of DPC.

Discussion

The current study may provide new insight into the exploration of DPC pathogenesis and the screened hub genes may serve as potential diagnostic indicator and novel therapeutic target.

Transcriptome landscapes of differentially expressed genes related to fat deposits in Nandan-Yao chicken

Nandan-Yao chicken is a Chinese native chicken with lower fat deposition and better meat quality. Fat deposition is a quite complex and important economic trait. However, its molecular mechanism is still unknown in chickens. In the current study, Nandan-Yao chicken was divided into two groups based on the rate of abdominal fat at 120 days old, namely the high-fat group and low-fat group. The total RNAs were isolated and sequenced by RNA sequencing (RNA-seq). After quality control, we gained 1222, 902, 784, 624, and 736 differentially expressed genes (DEGs) in abdominal fat, back skin, liver, pectoral muscle, and leg muscle, respectively. Analysis of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that significantly enriched GO term and KEGG signaling pathway mainly involved cytosolic ribosome, growth development, PPAR signaling pathway, Wnt signaling pathway, and linoleic acid metabolism in abdominal fat, back skin, and liver. While in pectoral muscle and leg muscle, it is mainly enriched in phosphatidylinositol signaling system, adrenergic signaling in cardiomyocytes, cytosolic ribosome, and cytosolic part. Sixteen genes were differentially expressed in all five tissues. Among them, PLA2G4A and RPS4Y1 might be the key regulators for fat deposition in Nandan-Yao chicken. The protein-protein interaction (PPI) network analysis of DEGs showed that PCK1 was the most notable genes. The findings in the current study will help to understand the regulation mechanism of abdominal fat and intramuscular fat in Nandan-Yao chicken and provide a theoretical basis for Chinese local chicken breeding.

Dynamic genome-wide association analysis and identification of candidate genes involved in anaerobic germination tolerance in rice

BACKGROUND

In Asian rice production, an increasing number of countries now choose the direct seeding mode because of rising costs, labour shortages and water shortages. The ability of rice seeds to undergo anaerobic germination (AG) plays an important role in the success of direct seeding.

RESULTS

In this study, we used 2,123,725 single nucleotide polymorphism (SNP) markers based on resequencing to conduct a dynamic genome-wide association study (GWAS) of coleoptile length (CL) and coleoptile diameter (CD) in 209 natural rice populations. A total of 26 SNP loci were detected in these two phenotypes, of which 5 overlapped with previously reported loci (S1_ 39674301, S6_ 20797781, S7_ 18722403, S8_ 9946213, S11_ 19165397), and two sites were detected repeatedly at different time points (S3_ 24689629 and S5_ 27918754). We suggest that these 7 loci (-log₁₀ (P) value > 7.3271) are the key sites that affect AG tolerance. To screen the candidate genes more effectively, we sequenced the transcriptome of the flooding-tolerant variety R151 in six key stages, including anaerobic (AN) and the oxygen conversion point (AN-A), and obtained high-quality differential expression profiles. Four reliable candidate genes were identified: Os01g0911700 (OsVP1), Os05g0560900 (OsGA2ox8), Os05g0562200 (OsDi19-1) and Os06g0548200. Then qRT-PCR and LC-MS/ MS targeting metabolite detection technology were used to further verify that the up-regulated expression of these four candidate genes was closely related to AG.

CONCLUSION

The four novel candidate genes were associated with gibberellin (GA) and abscisic acid (ABA) regulation and cell wall metabolism under oxygen-deficiency conditions and promoted coleoptile elongation while avoiding adverse effects, allowing the coleoptile to obtain oxygen, escape the low-oxygen environment and germinate rapidly. The results of this study improve our understanding of the genetic basis of AG in rice seeds, which is conducive to the selection of flooding-tolerant varieties suitable for direct seeding.

Integrated analysis of the methylome and transcriptome of chickens with fatty liver hemorrhagic syndrome

Background

DNA methylation, a biochemical modification of cytosine, has an important role in lipid metabolism. Fatty liver hemorrhagic syndrome (FLHS) is a serious disease and is tightly linked to lipid homeostasis. Herein, we compared the methylome and transcriptome of chickens with and without FLHS.

Results

We found genome-wide dysregulated DNA methylation pattern in which regions up- and down-stream of gene body were hypo-methylated in chickens with FLHS. A total of 4155 differentially methylated genes and 1389 differentially expressed genes were identified. Genes were focused when a negative relationship between mRNA expression and DNA methylation in promoter and gene body were detected. Based on pathway enrichment analysis, we found expression of genes related to lipogenesis and oxygenolysis (e.g., PPAR signaling pathway, fatty acid biosynthesis, and fatty acid elongation) to be up-regulated with associated down-regulated DNA methylation. In contrast, genes related to cellular junction and communication pathways (e.g., vascular smooth muscle contraction, phosphatidylinositol signaling system, and gap junction) were inhibited and with associated up-regulation of DNA methylation.

Conclusions

In the current study, we provide a genome-wide scale landscape of DNA methylation and gene expression. The hepatic hypo-methylation feature has been identified with FLHS chickens. By integrated analysis, the results strongly suggest that increased lipid accumulation and hepatocyte rupture are central pathways that are regulated by DNA methylation in chickens with FLHS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07305-3.

Acquired mutations and transcriptional remodeling in long-term estrogen-deprived locoregional breast cancer recurrences

BACKGROUND

Endocrine therapy resistance is a hallmark of advanced estrogen receptor (ER)-positive breast cancer. In this study, we aimed to determine acquired genomic changes in endocrine-resistant disease.

METHODS

We performed DNA/RNA hybrid-capture sequencing on 12 locoregional recurrences after long-term estrogen deprivation and identified acquired genomic changes versus each tumor's matched primary.

RESULTS

Despite being up to 7 years removed from the primary lesion, most recurrences harbored similar intrinsic transcriptional and copy number profiles. Only two genes, AKAP9 and KMT2C, were found to have single nucleotide variant (SNV) enrichments in more than one recurrence. Enriched mutations in single cases included SNVs within transcriptional regulators such as ARID1A, TP53, FOXO1, BRD1, NCOA1, and NCOR2 with one local recurrence gaining three PIK3CA mutations. In contrast to DNA-level changes, we discovered recurrent outlier mRNA expression alterations were common-including outlier gains in TP63 (n = 5 cases [42%]), NTRK3 (n = 5 [42%]), NTRK2 (n = 4 [33%]), PAX3 (n = 4 [33%]), FGFR4 (n = 3 [25%]), and TERT (n = 3 [25%]). Recurrent losses involved ESR1 (n = 5 [42%]), RELN (n = 5 [42%]), SFRP4 (n = 4 [33%]), and FOSB (n = 4 [33%]). ESR1-depleted recurrences harbored shared transcriptional remodeling events including upregulation of PROM1 and other basal cancer markers.

CONCLUSIONS

Taken together, this study defines acquired genomic changes in long-term, estrogen-deprived disease; highlights the importance of longitudinal RNA profiling; and identifies a common ESR1-depleted endocrine-resistant breast cancer subtype with basal-like transcriptional reprogramming.

[Related genes and characteristic analysis of trophoblast cells during early embryo developmental cessation]

Trophoblast cells play essential roles in the maintenance of normal embryo implantation, growth and development. The study of abnormal gene changes in trophoblastic cells from arrested embryos is helpful to understand the developmental mechanism of embryo developmental cessation or adverse pregnancy outcomes. In this study, we sequenced and analyzed the transcriptomes of the villi from ten women who have undergone abortion with either normal pregnancy or embryo development cessation. We found that there were 436 differentially expressed genes, of which 406 mRNA were significantly up-regulated and 32 mRNA were significantly down-regulated. Gene enrichment analysis showed that these genes were significantly enriched in immune-related functions and intercellular adhesion, such as lymphocyte activation, myeloid cell activation, extracellular matrix and collagen junction. And their potential regulatory pathways were enriched in terms of complement and coagulation cascade, extracellular matrix degradation. In addition, in this study the co-expression analysis of WGCNA was used to obtain the lncRNA with co-expression relationship with the differential genes. According to the different functions of the modules, two network diagrams were drawn, and four key genes were obtained, namely VSIG4, C1QC, CD36 and SPP1. These differential genes obtained in this study can be used as key molecules with potential effects on embryo development cessation. The enriched entries can provide a theoretical basis and new direction for further understanding of the etiology and mechanism of embryo development cessation or adverse pregnancy outcomes.

Identification of Translational microRNA Biomarker Candidates for Ketoconazole-Induced Liver Injury Using Next-Generation Sequencing

Drug-induced liver injury (DILI) is a leading cause of acute liver failure. Reliable and translational biomarkers are needed for early detection of DILI. microRNAs (miRNAs) have received wide attention as a novel class of potential DILI biomarkers. However, it is unclear how DILI drugs other than acetaminophen may influence miRNA expression or which miRNAs could serve as useful biomarkers in humans. We selected ketoconazole (KCZ), a classic hepatotoxin, to study miRNA biomarkers for DILI as a proof of concept for a workflow that integrated in vivo, in vitro, and bioinformatics analyses. We examined hepatic miRNA expression in KCZ-treated rats at multiple doses and durations using miRNA-sequencing and correlated our results with conventional DILI biomarkers such as liver histology. Significant dysregulation of rno-miR-34a-5p, rno-miR-331-3p, rno-miR-15b-3p, and rno-miR-676 was associated with cytoplasmic vacuolization, a phenotype in rat livers with KCZ-induced injury, which preceded the elevation of serum liver transaminases (ALT and AST). Between rats and humans, miR-34a-5p, miR-331-3p, and miR-15b-3p were evolutionarily conserved with identical sequences, whereas miR-676 showed 73% sequence similarity. Using quantitative PCR, we found that the levels of hsa-miR-34a-5p, hsa-miR-331-3p, and hsa-miR-15b-3p were significantly elevated in the culture media of HepaRG cells treated with 100 µM KCZ (a concentration that induced cytotoxicity). Additionally, we computationally characterized the miRNA candidates for their gene targeting, target functions, and miRNA/target evolutionary conservation. In conclusion, we identified miR-34a-5p, miR-331-3p, and miR-15b-3p as translational biomarker candidates for early detection of KCZ-induced liver injury with a workflow applicable to computational toxicology studies.

Sequencing, de novo assembly and annotation of Digitalis ferruginea subsp. schischkinii transcriptome

There is an increasing demand for elucidating the biosynthetic pathway of medicinal plants, which are capable of producing several metabolites with great potentials for industrial drug production. Digitalis species are important medicinal plants for the production of cardenolide compounds. Advancement on culture techniques is strictly related to our understanding of the genomic background of species. There are a limited number of genomic studies on Digitalis species. The goal of this study is to contribute to the genomic data of Digitalis ferruginea subsp. schischkinii by presenting transcriptome annotation. Digitalis ferruginea subsp. schischkinii has a limited distribution in Turkey and Transcaucasia, and has a high level of lanatoside C, an important cardenolide. In the study, we sequenced the cDNA library prepared from RNA pools of D. ferruginea subsp. schischkinii tissues treated with various stress conditions. Comprehensive bioinformatics approaches were used for de novo assembly and functional annotation of D. ferruginea subsp. schischkinii transcriptome sequence data along with TF families predictions and phylogenetic analysis. In the study, 58,369 unigenes were predicted and unigenes were annotated by analyzing the sequence data in the non-redundant (NR) protein database, the non-redundant nucleotide (NT) database, Gene Orthology (GO), EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), SwissProt, and InterPro databases. This study is the first transcriptome data for D. ferruginea subsp. schischkinii.

Transcriptome profiling unveils GAP43 regulates ABC transporters and EIF2 signaling in colorectal cancer cells

Background

The growth- and plasticity-associated protein-43 (GAP43) is biasedly expressed in indigestive system and nervous system. Recent study has shown that GAP43 is responsible for the development of neuronal growth and axonal regeneration in normal nervous tissue, while serves as a specific biomarker of relapsed or refractory neuroblastoma. However, its expression pattern and function in digestive system cancer remains to be clarified.

Methods

In this study, we examined the GAP43 status with qRT-PCR and bisulfite genomic sequencing in colorectal cancer (CRC). We investigated the effect of overexpressed GAP43 in CRC cells with RNA-seq. The RNA-seq data was analyzed with DAVID and IPA.

Results

GAP43 was downregulated in CRC compared to the adjacent tissues. DNA methylase inhibitor 5-Aza-CdR treatment could significantly induce GAP43, indicated that the silencing of GAP43 gene in CRC is closely related to DNA methylation. Bisulfite genomic sequencing confirmed the promoter methylation of GAP43 in CRC. To explore the transcriptional alterations by overexpressed GAP43 in CRC, we performed RNA-seq and found that upregulated genes were significantly enriched in the signaling pathways of ABC transporters and ECM-receptor interaction, while downregulated genes were significantly enriched in Ribosome signaling pathway. Further Ingenuity Pathway Analysis (IPA) showed that EIF2 signaling pathway was significantly repressed by overexpression of GAP43.

Conclusion

Our findings provide a novel mechanistic insight of GAP43 in CRC. Transcriptome profiling of overexpressed GAP43 in CRC uncovered the functional roles of GAP43 in the development of human CRC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-020-07728-x.

Effect of the AM Fungus Sieverdingia tortuosa on Common Vetch Responses to an Anthracnose Pathogen

Colletotrichum lentis Damm causes anthracnose in Vicia sativa L, otherwise known as common vetch. It was first reported in China in 2019. This study evaluates the effects of the arbuscular mycorrhizal (AM) fungus Sieverdingia tortuosa (N.C. Schenck & G.S. Sm.) Błaszk., Niezgoda, & B.T. Goto on growth and disease severity in common vetch. Our main finding is that the AM fungus increased root biomass and reduced anthracnose severity of common vetch. Responses correlated with defense, such as chitinase activity, polyphenol oxidase (PPO) activity, the concentrations of jasmonic acid and proline, and the expression of resistance-related genes (e.g., upregulated "signal transduction," "MAPK signaling pathway," "chitinase activity," "response to stress," and the KEGG pathways "phenylpropanoid biosynthesis," "MAPK signaling pathways," and "plant-pathogen interactions"), were also affected These findings provide insight into the mechanism by which this AM fungus regulates the defense response of common vetch to C. lentis.

RNA-Seq Provides New Insights into the Molecular Events Involved in "Ball-Skin versus Bladder Effect" on Fruit Cracking in Litchi

Fruit cracking is a disorder of fruit development in response to internal or external cues, which causes a loss in the economic value of fruit. Therefore, exploring the mechanism underlying fruit cracking is of great significance to increase the economic yield of fruit trees. However, the molecular mechanism underlying fruit cracking is still poorly understood. Litchi, as an important tropical and subtropical fruit crop, contributes significantly to the gross agricultural product in Southeast Asia. One important agricultural concern in the litchi industry is that some famous varieties with high economic value such as ‘Nuomici’ are susceptible to fruit cracking. Here, the cracking-susceptible cultivar ‘Nuomici’ and cracking-resistant cultivar ‘Huaizhi’ were selected, and the samples including pericarp and aril during fruit development and cracking were collected for RNA-Seq analysis. Based on weighted gene co-expression network analysis (WGCNA) and the “ball-skin versus bladder effect” theory (fruit cracking occurs upon the aril expanding pressure exceeds the pericarp strength), it was found that seven co-expression modules genes (1733 candidate genes) were closely associated with fruit cracking in ‘Nuomici’. Importantly, we propose that the low expression level of genes related to plant hormones (Auxin, Gibberellins, Ethylene), transcription factors, calcium transport and signaling, and lipid synthesis might decrease the mechanical strength of pericarp in ‘Nuomici’, while high expression level of genes associated with plant hormones (Auxin and abscisic acid), transcription factors, starch/sucrose metabolism, and sugar/water transport might increase the aril expanding pressure, thereby resulting in fruit cracking in ‘Nuomici’. In conclusion, our results provide comprehensive molecular events involved in the “ball-skin versus bladder effect” on fruit cracking in litchi.

Comparative RNA-Seq analysis reveals genes associated with masculinization in female Cannabis sativa

Using RNA profiling, we identified several silver thiosulfate-induced genes that potentially control the masculinization of female Cannabis sativa plants. Genetically female Cannabis sativa plants normally bear female flowers, but can develop male flowers in response to environmental and developmental cues. In an attempt to elucidate the molecular elements responsible for sex expression in C. sativa plants, we developed genetically female lines producing both female and chemically-induced male flowers. Furthermore, we carried out RNA-Seq assays aimed at identifying differentially expressed genes responsible for male flower development in female plants. The results revealed over 10,500 differentially expressed genes, of which around 200 potentially control masculinization of female cannabis plants. These genes include transcription factors and other genes involved in male organ (i.e., anther and pollen) development, as well as genes involved in phytohormone signalling and male-biased phenotypes. The expressions of 15 of these genes were further validated by qPCR assay confirming similar expression patterns to that of RNA-Seq data. These genes would be useful for understanding predisposed plants producing flowers of both sex types in the same plant, and help breeders to regulate the masculinization of female plants through targeted breeding and plant biotechnology.

Comparative transcriptomic and metabolomic analyses of carotenoid biosynthesis reveal the basis of white petal color in Brassica napus

The molecular mechanism underlying white petal color in Brassica napus was revealed by transcriptomic and metabolomic analyses. Rapeseed (Brassica napus L.) is one of the most important oilseed crops worldwide, but the mechanisms underlying flower color in this crop are known less. Here, we performed metabolomic and transcriptomic analyses of the yellow-flowered rapeseed cultivar 'Zhongshuang 11' (ZS11) and the white-flowered inbred line 'White Petal' (WP). The total carotenoid contents were 1.778-fold and 1.969-fold higher in ZS11 vs. WP petals at stages S2 and S4, respectively. Our findings suggest that white petal color in WP flowers is primarily due to decreased lutein and zeaxanthin contents. Transcriptome analysis revealed 10,116 differentially expressed genes with a fourfold or greater change in expression (P-value less than 0.001) in WP vs. ZS11 petals, including 1,209 genes that were differentially expressed at four different stages and 20 genes in the carotenoid metabolism pathway. BnNCED4b, encoding a protein involved in carotenoid degradation, was expressed at abnormally high levels in WP petals, suggesting it might play a key role in white petal formation. The results of qRT-PCR were consistent with the transcriptome data. The results of this study provide important insights into the molecular mechanisms of the carotenoid metabolic pathway in rapeseed petals, and the candidate genes identified in this study provide a resource for the creation of new B. napus germplasms with different petal colors.

Identification of a novel therapeutic candidate, NRK, in primary cancer-associated fibroblasts of lung adenocarcinoma microenvironment

PURPOSE

Lung adenocarcinoma (LUAD) accounts for approximately half of patients in lung cancer. Cancer-associated fibroblasts (CAFs) are the major component in the tumor microenvironment (TME). Targeting CAFs is a promising therapeutic strategy for cancer treatment. However, therapeutic targets of CAFs in LUAD remains largely unclear.

METHODS

Seven CAFs and nine normal fibroblasts (NFs) were isolated from tumor and paratumor tissues of LUAD patients undergoing surgery, respectively. RNA-seq and bioinformatics analysis were performed to identify the differentially expressed genes (DEGs) and their functions in CAFs compared with NFs. DEGs of ten overlaying were obtained from RNA-seq, our previously reported lncRNA microarray and public datasets (E-MTAB-6149, E-MTAB-6653) and validated by RT-qPCR. Nik-related kinase (NRK) was further validated by RT-qPCR, immunofluorescence (IF), Western Blot (WB) in vitro, and in Cancer Cell Line Encyclopedia (CCLE) database. Survival analysis was performed on Kaplan-Meier plotter.

RESULTS

A total of 1799 DEGs were identified, including 650 upregulated DEGs and 1149 downregulated DEGs. The upregulated and downregulated DEGs were mostly enriched in extracellular matrix (ECM) functions and in glycolysis/gluconeogenesis pathways. Interestingly, NRK was the most significantly upregulated overlaying DEGs which was rarely associated with CAFs before. NRK was predominantly expressed in CAFs, but weakly expressed in NFs, normal lung bronchial epithelial cell line BEAS-2B, LUAD cell lines A549 and H1299, as well as in the majority of 191 lung cancer cell lines including LUAD. Moreover, elevated NRK predicted poor survival in LUAD patients.

CONCLUSION

Here, we first report that NRK is significantly elevated in LUAD-associated CAFs and may function as a promising therapeutic target for cancer combination treatment. Besides, modulation of ECM and glycolysis/gluconeogenesis pathways may be an efficient approach to alter CAFs functionality in LUAD.

Deciphering the Dynamic Landscape of Transcription-Associated mRNP Quality Control Components Over the Whole Yeast Genome

In eukaryotic cells, aberrant mRNPs with processing and packaging defects are targeted co-transcriptionally by a surveillance system that triggers their nuclear retention and ultimately the degradation of their mRNA component by the 3'-5' activity of the exosome-associated exonuclease Rrp6. This mRNP quality control process is stimulated by the NNS complex (Nrd1-Nab3-Sen1), which otherwise mediates termination, processing, and decay of ncRNAs. The process involves also the exosome co-activator TRAMP complex (Trf4-Air2-Mtr4). Here, we describe a genome-wide approach to visualize the dynamic movement and coordination of these quality control components over the yeast chromosomes upon perturbation of mRNP biogenesis. The method provides valuable information on how the surveillance system is precisely coordinated both physically and functionally with the transcription machinery to detect the faulty events during perturbation of mRNP biogenesis. The overview shows also that the gathering of the quality control components over affected mRNA genes takes place at the expense of their commitment to be recruited at ncRNA genomic features, provoking termination and processing defects of ncRNAs.

Engineering Sensors for Gene Expression Burden

RNA-seq enables the analysis of gene expression profiles across different conditions and organisms. Gene expression burden slows down growth, which results in poor predictability of gene constructs and product yields. Here, we describe how we applied RNA-seq to study the transcriptional profiles of Escherichia coli when burden is elicited during heterologous gene expression. We then present how we selected early responsive promoters from our RNA-seq results to design sensors for gene expression burden. Finally, we describe how we used one of these sensors to develop a burden-driven feedback regulator to improve cellular fitness in engineered E. coli.

Characterizing Genetic Parts and Devices Using RNA Sequencing

Synthetic genetic circuits are composed of many parts that must interact and function together to produce a desired pattern of gene expression. A challenge when assembling circuits is that genetic parts often behave differently within a circuit, potentially impacting the desired functionality. Existing debugging methods based on fluorescent reporter proteins allow for only a few internal states to be monitored simultaneously, making diagnosis of the root cause impossible for large systems. Here, we present a tool called the Genetic Analyzer which uses RNA sequencing data to simultaneously characterize all transcriptional parts (e.g., promoters and terminators) and devices (e.g., sensors and logic gates) in complex genetic circuits. This provides a complete picture of the inner workings of a genetic circuit enabling faults to be easily identified and fixed. We construct a complete workflow to coordinate the execution of the various data processing and analysis steps and explain the options available when adapting these for the characterization of new systems.