Integrated Transcriptome Analysis Research Articles

Integrative metabolome and transcriptome analyses reveal the differences in flavonoid and terpenoid synthesis between Glycyrrhiza uralensis (licorice) leaves and roots.

The online version contains supplementary material available at 10.1007/s10068-023-01467-y.

Integrative physiological, transcriptomic and metabolomic analysis reveals how the roots of two ornamental Hydrangea macrophylla cultivars cope with lead (Pb) toxicity

Lead (Pb) soil contamination has caused serious ecological and environmental issues. Hydrangea macrophylla is a potential Pb-contaminated soil remediation plant, however, their Pb stress defense mechanism is largely unknown. Here, the physiology, transcriptomic and metabolome of two H. macrophylla cultivars (ML, Pb-sensitive cultivar; JC, Pb-resistant cultivar) under Pb stress were investigated. The results demonstrated that JC performed superiorly, with activities of the antioxidant enzymes superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) were 1.25, 2.84, and 1.67 times higher than those of ML after Pb treatment, respectively, and the amount of soluble sugar in JC increased by 231.34 % compared with that in ML. The electrical conductivity (EC) value of the root exudates of JC was 43.71 % lower than that of ML under Pb stress. The non-targeted metabolomics analysis revealed 193 metabolites grouped into nine categories. Pb stress-induced differential expression of the 37 metabolites, among which the major metabolites up-regulated in ML were organic acids, while in JC, these were carbohydrates, fatty acids, organic acids and lipids. The transcriptomic analysis revealed that Pb exposure induced 1075 and 1314 differentially expressed genes (DEGs) in JC and ML, respectively. According to the functional annotation results, hub genes were primarily enriched in carbohydrate metabolism, root growth, and plant resistance to external stresses. A conjoint analysis of the two omics indicated that the cutin, suberine and wax biosynthesis pathway in JC played an essential role in Pb detoxification. These findings clarify the resistance mechanism of H. macrophylla to Pb stress and open up a new avenue for breeding H. macrophylla Pb-resistant cultivars.

Adipogenic characterization of immortalized CD55+ progenitor cells from human white adipose tissue

ABSTRACT Background Mature adipocytes are notoriously difficult to study ex vivo and alternative cell culture systems have therefore been developed. One of the most common models are human adipose progenitor cells (hAPCs). Unfortunately, these display replicative senescence after prolonged culture conditions, which limits their use in mechanistic studies. Methods Herein, we knocked in human telomerase reverse transcriptase (TERT) into the AAVS1 locus of CD55+ hAPCs derived from abdominal subcutaneous adipose tissue and characterized the cells before and after differentiation into adipocytes. Results Immortalized TERT-hAPCs retained proliferative and adipogenic capacities comparable to those of early-passage wild type hAPCs for > 80 passages. In line with this, our integrative transcriptomic and proteomic analyses revealed that TERT-hAPCs displayed robust adipocyte expression profiles in comparison to wild type hAPCs. This was confirmed by functional analyses of lipid turnover where TERT-hAPCs exhibited pronounced responses to insulin and pro-lipolytic stimuli such as isoprenaline, dibutyrul cAMP and tumour necrosis factor alpha. In addition, TERT-hAPCs could be readily cultured in both standard 2D and 3D-cultures and proteomic analyses revealed that the spheroid culture conditions improved adipogenesis. Conclusion Through descriptive and functional studies, we demonstrate that immortalization of human CD55+ hAPCs is feasible and results in cells with stable proliferative and adipogenic capacities over multiple passages. As these cells are cryopreservable, they provide the additional advantage over primary cells of allowing repeated studies in both 2D and 3D model systems with the same genetic background. (234/250)

Integrative whole-genome and transcriptome analysis of HER2-amplified metastatic breast cancer

BackgroundIn breast cancer, the advent of anti-HER2 therapies has made HER2+ tumors a highly relevant subgroup. However, the exact characteristics which prohibit clinical response to anti-HER2 therapies and drive disease progression are not yet fully known. Integrative whole-genome and transcriptomic sequencing data from both primary and metastatic HER2-positive breast cancer will enhance our understanding of underlying biological processes.MethodsHere, we used WGS and RNA sequencing data of 700 metastatic breast tumors, of which 68 being HER2+, to search for specific genomic features of HER2+ disease and therapy resistance. Furthermore, we integrated results with transcriptomic data to associate tumors exhibiting a HER2+-specific gene expression profile with ERBB2 mutation status, prior therapy and relevant gene expression signatures.ResultsOverall genomic profiles of primary and metastatic HER2+ breast cancers were similar, and no specific acquired genomics traits connected to prior anti-HER2 treatment were observed. However, specific genomic features were predictive of progression-free survival on post-biopsy anti-HER2 treatment. Furthermore, a HER2-driven expression profile grouped HER2-amplified tumors with ERBB2-mutated cases and cases without HER2 alterations. The latter were reported as ER positive in primary disease, but the metastatic biopsy showed low ESR1 expression and upregulation of the MAPK pathway, suggesting transformation to ER independence.ConclusionsIn summary, although the quantity of variants increased throughout HER2-positive breast cancer progression, the genomic composition remained largely consistent, thus yielding no new major processes beside those already operational in primary disease. Our results suggest that integrated genomic and transcriptomic analyses may be key in establishing therapeutic options.

Integrative transcriptomic and network pharmacology analysis reveals the neuroprotective role of BYHWD through enhancing autophagy by inhibiting Ctsb in intracerebral hemorrhage mice.

In this study, we aimed to combine transcriptomic and network pharmacology to explore the crucial mRNAs and specific regulatory molecules of Buyang Huanwu Decoction (BYHWD) in intracerebral hemorrhage (ICH) treatment. C57BL/6 mice were randomly divided into three groups: sham, ICH, and BYHWD. BYHWD (43.29g/kg) was administered once a day for 7days. An equal volume of double-distilled water was used as a control. Behavioural and histopathological experiments were conducted to confirm the neuroprotective effects of BYHWD. Brain tissues were collected for transcriptomic detection. Bioinformatics analysis were performed to illustrate the target gene functions. Network pharmacology was used to predict potential targets for BYHWD. Next, transcriptomic assays were combined with network pharmacology to identify the potential differentially expressed mRNAs. Immunofluorescence staining, real-time polymerase chain reaction, western blotting, and transmission electron microscopy were performed to elucidate the underlying mechanisms. BYHWD intervention in ICH reduced neurological deficits. Network pharmacology analysis identified 203 potential therapeutic targets for ICH, whereas transcriptomic assay revealed 109 differentially expressed mRNAs post-ICH. Among these, cathepsin B, ATP binding cassette subfamily B member 1, toll-like receptor 4, chemokine (C-C motif) ligand 12, and baculoviral IAP repeat-containing 5 were identified as potential target mRNAs through the integration of transcriptomics and network pharmacology approaches. Bioinformatics analysis suggested that the beneficial effects of BYHWD in ICH may be associated with apoptosis, animal autophagy signal pathways, and PI3K-Akt and mTOR biological processes. Furthermore, BYHWD intervention decreased Ctsb expression levels and increased autophagy levels in ICH. Animal experiments in combination with bioinformatics analysis confirmed that BYHWD plays a neuroprotective role in ICH by regulating Ctsb to enhance autophagy.

Integrative analysis of non-targeted metabolome and transcriptome reveals the mechanism of volatile formation in pepper fruit.

Introduction: Aroma is a key inherent quality attributes of pepper fruit, yet the underlying mechanisms of aroma compound biosynthesis remain unclear. Methods: In this study, the volatile profile of the QH (cultivated Capsicum chinense) and WH (cultivated Capsicum annuum) pepper varieties were putatively identified during fruit development using gas chromatography-mass spectrometry (GC-MS). Results and discussion: The results identified 203 volatiles in pepper, and most of the esters, terpenes, aldehydes and alcohols were significantly down-regulated with fruit ripening. The comparison of volatile components between varieties revealed that aldehydes and alcohols were highly expressed in the WH fruit, while esters and terpenes with fruity or floral aroma were generally highly accumulated in the QH fruit, providing QH with a fruity odor. Transcriptome analysis demonstrated the close relationship between the synthesis of volatiles and the fatty acid and terpene metabolic pathways, and the high expression of the ADH, AAT and TPS genes was key in determining the accumulation of volatiles in pepper fruit. Furthermore, integrative metabolome and transcriptome analysis revealed that 208 differentially expressed genes were highly correlated with 114 volatiles, and the transcription factors of bHLH, MYB, ARF and IAA were identified as fundamental for the regulation of volatile synthesis in pepper fruit. Our results extend the understanding of the synthesis and accumulation of volatiles in pepper fruit.

Integrative transcriptomic and metabolomic analyses reveal the phenylpropanoid and flavonoid biosynthesis of Prunus mume.

Prunus mume is an important medicinal plant with ornamental and edible value. Its flowers contain phenylpropanoids, flavonoids and other active components, that have important medicinal and edible value, yet their molecular regulatory mechanisms in P. mume remain unclear. In this study, the content of total flavonoid and total phenylpropanoid of P. mume at different developmental periods was measured first, and the results showed that the content of total flavonoid and total phenylpropanoid gradually decreased in three developmental periods. Then, an integrated analysis of transcriptome and metabolome was conducted on three developmental periods of P. mume to investigate the law of synthetic accumulation for P. mume metabolites, and the key enzyme genes for the biosynthesis of phenylpropanoids and flavonoids were screened out according to the differentially expressed genes (DEGs). A total of 14,332 DEGs and 38 differentially accumulate metabolites (DAMs) were obtained by transcriptomics and metabolomics analysis. The key enzyme genes and metabolites in the bud (HL) were significantly different from those in the half-opening (BK) and full-opening (QK) periods. In the phenylpropanoid and flavonoid biosynthesis pathway, the ion abundance of chlorogenic acid, naringenin, kaempferol, isoquercitrin, rutin and other metabolites decreased with the development of flowers, while the ion abundance of cinnamic acid increased. Key enzyme genes such as HCT, CCR, COMT, CHS, F3H, and FLS positively regulate the downstream metabolites, while PAL, C4H, and 4CL negatively regulate the downstream metabolites. Moreover, the key genes FLS (CL4312-2, CL4312-3, CL4312-4, CL4312-5, CL4312-6) regulating the synthesis of flavonols are highly expressed in bud samples. The dynamic changes of these metabolites were validated by determining the content of 14 phenylpropanoids and flavonoids in P. mume at different developmental periods, and the transcription expression levels of these genes were validated by real-time PCR. Our study provides new insights into the molecular mechanism of phenylpropanoid and flavonoid accumulation in P. mume.

Genetic and immunological insights into COVID-19 with acute myocardial infarction: integrated analysis of mendelian randomization, transcriptomics, and clinical samples.

Globally, most deaths result from cardiovascular diseases, particularly ischemic heart disease. COVID-19 affects the heart, worsening existing heart conditions and causing myocardial injury. The mechanistic link between COVID-19 and acute myocardial infarction (AMI) is still being investigated to elucidate the underlying molecular perspectives. Genetic risk assessment was conducted using two-sample Mendelian randomization (TSMR) to determine the causality between COVID-19 and AMI. Weighted gene co-expression network analysis (WGCNA) and machine learning were used to discover and validate shared hub genes for the two diseases using bulk RNA sequencing (RNA-seq) datasets. Additionally, gene set enrichment analysis (GSEA) and single-cell RNA-seq (scRNA-seq) analyses were performed to characterize immune cell infiltration, communication, and immune correlation of the hub genes. To validate the findings, the expression patterns of hub genes were confirmed in clinical blood samples collected from COVID-19 patients with AMI. TSMR did not find evidence supporting a causal association between COVID-19 or severe COVID-19 and AMI. In the bulk RNA-seq discovery cohorts for both COVID-19 and AMI, WGCNA's intersection analysis and machine learning identified TLR4 and ABCA1 as significant hub genes, demonstrating high diagnostic and predictive value in the RNA-seq validation cohort. Single-gene GSEA and single-sample GSEA (ssGSEA) revealed immune and inflammatory roles for TLR4 and ABCA1, linked to various immune cell infiltrations. Furthermore, scRNA-seq analysis unveiled significant immune dysregulation in COVID-19 patients, characterized by altered immune cell proportions, phenotypic shifts, enhanced cell-cell communication, and elevated TLR4 and ABCA1 in CD16 monocytes. Lastly, the increased expression of TLR4, but not ABCA1, was validated in clinical blood samples from COVID-19 patients with AMI. No genetic causal link between COVID-19 and AMI and dysregulated TLR4 and ABCA1 may be responsible for the development of immune and inflammatory responses in COVID-19 patients with AMI.

Integrated analysis of miRNA, transcriptome, and degradome sequencing provides new insights into lipid metabolism in perilla seed

MicroRNAs (miRNA) are small noncoding RNAs that play a crucial as molecular regulators in lipid metabolism in various oil crops. Perilla (Perilla frutescens) is a specific oil crop known for its high alpha-linolenic acid (C18:3n3, ALA) content (>65 %) in their seed oils. In view of the regulatory mechanism of miRNAs in perilla remains unclear, we conducted miRNAs and transcriptome sequencing in two cultivars with distinct lipid compositions. A total of 525 unique miRNAs, including 142 differentially expressed miRNAs was identified in perilla seeds. The 318 miRNAs targeted 7,761 genes. Furthermore, we identified 112 regulated miRNAs and their 610 target genes involved in lipid metabolism. MiR159b and miR167a as the core nodes to regulate the expression of genes in oil biosynthesis (e.g., KAS, FATB, GPAT, FAD, DGK, LPAAT) and key regulatory TFs (e.g., MYB, ARF, DOF, SPL, NAC, TCP, and bHLH). The 1,219 miRNA-mRNA regulation modules were confirmed through degradome sequencing. Notably, pf-miR159b-MYBs and pf-miR167a-ARFs regulation modules were confirmed. They exhibited significantly different expression levels in two cultivars and believed to play important roles in oil biosynthesis in perilla seeds. This provides valuable insights into the functional analysis of miRNA-regulated lipid metabolism in perilla seeds.

Integrative transcriptome analysis reveals alternative polyadenylation potentially contributes to GCRV early infection.

Grass carp reovirus (GCRV), a member of the Aquareovirus genus in the Reoviridae family, is considered to be the most pathogenic aquareovirus. Productive viral infection requires extensive interactions between viruses and host cells. However, the molecular mechanisms underlying GCRV early infection remains elusive. In this study we performed transcriptome and DNA methylome analyses with Ctenopharyngodon idellus kidney (CIK) cells infected with GCRV at 0, 4, and 8 h post infection (hpi), respectively. We found that at early infection stage the differentially expressed genes related to defense response and immune response in CIK cells are activated. Although DNA methylation pattern of CIK cells 8 hpi is similar to mock-infected cells, we identified a considerable number of genes that selectively utilize alternative polyadenylation sites. Particularly, we found that biological processes of cytoskeleton organization and regulation of microtubule polymerization are statistically enriched in the genes with altered 3'UTRs. Our results suggest that alternative polyadenylation potentially contributes to GCRV early infection.

Integrated analysis of transcriptome and metabolome reveals insights for low-temperature germination in hybrid rapeseeds (Brassica napus L.)

Rapeseed (Brassica napus L.) is an important oil-producing crop in China. However, cold stress in winter can adversely affect rapeseed germination and subsequently result in poor seed yield at the mature stage. Studies of differences in the transcriptional and metabolic levels of rapeseed under cold stress can improve our understanding of low-temperature germination (LTG). The current study aimed to identify the cold stress-responsive genes, metabolites, and metabolic pathways based on a combined transcriptome and metabolome analysis to understand the difference of LTG and tolerance mechanisms in the cold-tolerant (Yueyou1301, YY1301) and cold-normal (Fengyou737, FY737) rapeseed varieties. Compared to FY737, YY1301 had a higher germination rate, indole acetic acid (IAA) and gibberellic acid (GA)/(abscisic acid) ABA levels at 7.5 °C. A total of 951 differentially expressed genes (DEGs) and 86 differentially accumulated metabolites (DAMs) were identified in two rapeseed varieties. Conjoint analysis revealed 12 DAMs and 5 DEGs that were strongly correlated in inducing rapeseed LTG, which were mainly related to carbohydrate and amino acid metabolism, specifically the pathway of glutathione metabolism and starch and sucrose metabolism. These results suggest that the DAMs and DEGs involved in crucial biological pathways may regulate the LTG of rapeseed. It increases the understanding of the molecular mechanisms underlying the adaptation of rapeseed to LTG.

Integrative Genomic and Transcriptomic Analysis Reveals Targetable Vulnerabilities in Angioimmunoblastic T-Cell Lymphoma

Follicular helper T-cell lymphoma of the angioimmunoblastic type (AITL) is associated with dismal prognosis. We performed functional genomic approaches including whole-exome sequencing (WES; n=119), transcriptomic (n=78) and methylation (n=40) analysis. We identified recurrent mutations in known epigenetic drivers ( TET2, DNMT3A, IDH2 R172), and also identified novel ones (TET3, KMT2D). Somatic mutation of all three epigenetic drivers ( TET2, IDH2, and DNMT3A) was associated with poor prognosis (p<.001). Mutations in genes regulating T-cell receptor (TCR) signaling ( CD28 VAV1, FYN, PLCG1) or activation ( RHOA G17V), and regulators of the PI3K pathway (PIK(3)C members, PTEN,PHLPP-1/-2) were also found. Genome-wide DNA-methylation analysis integrated with mRNA expression profiling also revealed epigenetic alterations in genes regulating TCR-RHOA/B/C or PI3K-signaling. TET2 loss was noted in 85% AITLs and was significantly associated with RHOA G17V, CD28 and IDH2 R172mutations. AITLs lacking RHOA G17V tended to have mutations regulating the JAK-STAT pathway ( JAK2, JAK3, STAT1, STAT3, SOCS1). RNA-seq analysis identified fusion transcripts in genes regulating TCR activation (8%), revealed a restricted TCR repertoire in the majority of cases (a=87%, b=72%), and showed the presence of Epstein-Barr virus transcriptome (73%). GEP demonstrated association of B-cells in the tumor-milieu with better prognosis (p=.006), while dendritic cells were associated with worse prognosis (p=.001), which was further validated by immunohistochemistry using CD20, CD68, and CD163 antibodies. RNA-seq and corresponding WES analysis of 12 AITL patient-derived-xenografts (PDX) showed that bi-allelic TET2 mutations, DNMT3A mutations or sub-clonal mutations ( PLCG1 PHLPP2) werepropagated in sequential passages. Gene signatures related to T FH (follicular helper) and T CM (central memory) were also well-maintained in secondary passages in PDX models. Gene signatures of late PDX passages (3 rd-5 th) were enriched with genes related to proliferation and metabolic reprogramming, and in an independent cohort of AITLs, high expression of T3/T5 related signatures was associated with worse outcome (p=0.02/p=0.009). Low mRNA expression of PHLPP2 predicted poor prognosis (p=.03) and engineered PHLPP2 loss showed enhanced PI(3)K activation and FOXO1 inactivation in CD4+ T-cells in-vitro. Thus, we defined the genomic landscape for AITL, which is largely characterized by epigenetic alterations, TCR signaling and PI3K/AKT dysregulation, which may be amenable for therapeutic targeting.

Identification of a Novel Target Implicated in Chronic Obstructive Sleep Apnea-Related Atrial Fibrillation by Integrative Analysis of Transcriptome and Proteome.

This study aimed to identify a newly identified target involved in atrial fibrillation (AF) linked to chronic obstructive sleep apnea (COSA) through an integrative analysis of transcriptome and proteome. Fifteen beagle canines were randomly assigned to three groups: control (CON), obstructive sleep apnea (OSA), and OSA with superior left ganglionated plexi ablation (OSA+GP). A COSA model was established by intermittently obstructing the endotracheal cannula during exhalation for 12 weeks. Left parasternal thoracotomy through the fourth intercostal space allowed for superior left ganglionated plexi (SLGP) ablation. In vivo open-chest electrophysiological programmed stimulation was performed to assess AF inducibility. Histological, transcriptomic, and proteomic analyses were conducted on atrial samples. After 12 weeks, the OSA group exhibited increased AF inducibility and longer AF durations compared to the CON group. Integrated transcriptomic and proteomic analyses identified 2422 differentially expressed genes (DEGs) and 1194 differentially expressed proteins (DEPs) between OSA and CON groups, as well as between OSA+GP and OSA groups (1850 DEGs and 1418 DEPs). The analysis revealed that differentially regulated DEGs were primarily enriched in mitochondrial biological processes in the CON-vs.-OSA and OSA-vs.-GP comparisons. Notably, the key regulatory molecule GSTZ1 was activated in OSA and inhibited by GP ablation. These findings suggest that GSTZ1 may play a pivotal role in mitochondrial damage, triggering AF substrate formation, and increasing susceptibility to AF in the context of COSA.

Integrated DNA Methylation and Transcriptomics Analyses of Lacrimal Glands Identify the Potential Genes Implicated in the Development of Sjögren's Syndrome-Related Dry Eye.

Sjögren's syndrome-related dry eye (SS-related dry eye) is an intractable autoimmune disease characterized by chronic inflammation of lacrimal glands (LGs), where epigenetic factors are proven to play a crucial role in the pathogenesis of this disease. However, the alteration of DNA methylation in LGs and its role in the pathogenesis of SS-related dry eye is still unknown. Here, we performed an integrated analysis of DNA methylation and RNA-Seq data in LGs to identify novel DNA methylation-regulated differentially expressed genes (MeDEGs) in the pathogenesis of SS-related dry eye. The DNA methylation and transcription profiles of LGs in NOD mice at different stages of SS-related dry eye (4-, 8-, 12- and 16 weeks old) were generated by reduced representation bisulfite sequencing (RRBS) and RNA-Seq. The differentially methylated genes (DMGs) and differentially expressed genes (DEGs) were analyzed by MethylKit R package and edgeR. Correlation analysis between methylation level and mRNA expression was conducted with R software. The functional correlation of DMGs and DEGs was analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Finally, LG tissues from another litter of NOD mice were collected for methylation-specific polymerase chain reaction (MSP) and quantitative real-time PCR (qRT-PCR) to validate the methylation and expression levels of key genes. CD4+ cell infiltration of LGs was detected by immunofluorescence staining. Hypermethylation of LGs was identified in NOD mice with the progression of SS-related dry eye and the DMGs were mainly enriched in the GTPases activation and Ras signaling pathway. RNA-seq analysis revealed 1321, 2549, and 3712 DEGs in the 8-, 12- and 16-week-old NOD mice compared with 4-week-old normal control mice. For GO analysis, the DEGs were mainly enriched in T cell immune responses. Further, a total of 140 MeDEGs were obtained by integrated analysis of methylome and transcriptome, which were primarily enriched in T cell activation, proliferation and differentiation. Based on the main GO terms and KEGG pathways of MeDEGs, 8 genes were screened out. The expression levels of these key genes, especially Itgal, Vav1, Irf4 and Icosl, were verified to elevate after the onset of SS-related dry eye in NOD mice and positively correlated with the extent of inflammatory cell infiltration in LGs. Immunofluorescence assay revealed that CD4+ cell infiltration dramatically increased in LGs of SS-related dry eye mice compared with the control mice. And the expression levels of four genes showed significantly positive correlation with the extent of CD4+ cell infiltration in LGs. MSP showed the hypomethylation of the Irf4 and Itgal promoters in NOD mice with SS-related dry eye compared to control group. Our study revealed the critical role of epigenetic regulation of T cell immunity-related genes in the progression of SS-related dry eye and reminded us that DNA methylation-regulated genes such as Itgal, Vav1, Irf4 and Icosl may be used as new targets for SS-related dry eye therapy.

Integrative metabolome and transcriptome analyses reveal the mechanism by which white light increases apple phenolics

Apple is a photophilic fruit and the nutrients content such as phenolics loss is a major problem during fruit storage. Light can affect the yield, appearance and quality of apples. However, the molecular mechanism behind how white light induces phenolics accumulation to improve quality in fruits remains unclear. In this study, we discovered that white light with an intensity of 4 w/m2 can enhance the total phenolic and flavonoid content in postharvest apple fruits. To investigate the specific phenolics that undergo changes during this process, we analyzed the metabolomes of apples exposed to white light for 5 days after harvesting. We identified a total of 97 differentially accumulated metabolites (DAMs), including 15 phenolics. Isorhamnetin-O-hexos, tangereti, phloretin, cis-p-coumaroyltartaric acid and cis-p-coumaroyltartaric acid contents decreased in white light. While myricetin 3-O-galactoside, astilbin, sieboldin 2′-hydroxygenistein, syringaresinol, pinoresino, maleoyl-caffeoylquinic acid, rosmarinyl glucosid, coniferyl alcohol, p-coumaryl alcoho, and feruloylferuloyltartaric acid contents increased in white light. Additionally, RNA-seq analysis allowed us to identify 1988 differentially expressed genes (DEGs), which included three key enzymes (phenylalanine ammonia-lyase [MdPAL-1], chalcone isomerase [MdCHI-like], and 4-coumarate: -CoA ligase [Md4CL-1]) involved in the phenolic synthesis pathway, as well as a basic leucine zipper domain (bZIP) transcription factor known as MdbZIP67. Further investigations revealed that white light can regulate the transcription of MdPAL-1, Md4CL-1, and MdCHI-like genes through MdbZIP67, resulting in an increase in phenolic content during fruit storage. Here, we have identified a new method to enhance the phenolic content in apple fruits with characteristics of easy to use and low cost. And we have elucidated the mechanism by which white light affects phenolic content during the fruit's shelf life.

Integrative transcriptome and whole-genome bisulfite sequencing analyses of a temperature-sensitive albino tea plant cultivar.

Green tea made from albino buds and leaves has a strong umami taste and aroma. The cultivar 'Zhonghuang 2' (ZH2, Camellia sinensis) is a natural mutant with young shoots that are yellow in spring and green or yellow-green in summer. However, the mechanism of leaf color change remains unclear. Here, we found that young shoots of ZH2 were yellow at low temperature (LT) and green at high temperature (HT), indicating that ZH2 is a temperature-sensitive cultivar. Transmission electron microscopy analysis showed that the grana in the chloroplasts of young shoots grown at LT were poorly stacked, which caused a lack of photoreactions and chlorophyll. RNA-seq results showed 1279 genes differentially expressed in the young shoots grown at LT compared with those at HT, including genes related to cytochrome synthesis, chloroplast development, photosynthesis, and DNA methylation. A whole-genome bisulfite sequencing assay revealed that the dynamics of DNA methylation levels in the CG, CHG, and CHH contexts decreased under LT, and the change was most obvious in the CHH context. Furthermore, 72 genes showed significant changes in both expression and DNA methylation levels, and most of them were related to cytochrome synthesis, chloroplast development, photosynthesis, transcription factors, and signaling pathways. These results demonstrate that DNA methylation is involved in the LT-regulated albino processes of ZH2. Changes in DNA methylation levels were associated with changes in gene expression levels, affecting the structure and function of chloroplasts, which may have a phenotypic impact on shoot and leaf color.

Integrative transcriptomic and proteomics profiling analyses reveals PLK1 as a target of curcumol in hepatocellular carcinoma cells

Curcumol is a natural drug isolated from Chinese medicinal plant Rhizoma Curcumae. Though many studies have revealed its anti-tumor mechanism, there are still few studies reveal its comprehensive anti-cancer effect used multi-omics method. In this study, we aim to comprehensively analyze its anti-cancer mechanism and find its potential target in human hepatocellular carcinoma cells (HCC). In total, 1330 differential genes and 1381 differential proteins were significantly altered by curcumol in HCC. Transcriptomics and proteomics combined analysis revealed that 47 downregulated of the 97 common molecules regulated by curcumol were mainly enriched in the cell cycle pathway, including Cell Division Cycle Associated 2 (CDCA2), CDCA3, CDCA8, CDC25C, S-phase kinase-associated protein 2 (SKP2) and Polo-like kinase 1(PLK1). Furthermore, silenced PLK1 inhibited the cell growth and reversed curcumol’s effect on HCC. This paper provided a comprehensive understanding of the effect of curcumol and identified PLK1 was its potential target in HCC cells.

Integrated Transcriptome Analysis of miRNAs and mRNAs in the Skeletal Muscle of Wuranke Sheep.

MicroRNAs (miRNAs) are regarded as important regulators in skeletal muscle development. To reveal the regulatory roles of miRNAs and their target mRNAs underlying the skeletal muscle development of Wuranke sheep, we investigated the miRNA and mRNA expression profiles in the biceps femoris of these sheep at the fetal (3 months of gestation) and 3- and 15-month-old postnatal stages. Consequently, a total of 1195 miRNAs and 24,959 genes were identified. Furthermore, 474, 461, and 54 differentially expressed miRNAs (DEMs) and 6783, 7407, and 78 differentially expressed genes (DEGs) were detected among three comparative groups. Functional analysis demonstrated that the target mRNAs of the DEMs were enriched in multiple pathways related to muscle development. Moreover, the interactions among several predicted miRNA-mRNA pairs (oar-miR-133-HDAC1, oar-miR-1185-5p-MYH1/HADHA/OXCT1, and PC-5p-3703_578-INSR/ACTG1) that potentially affect skeletal muscle development were verified using dual-luciferase reporter assays. In this study, we identified the miRNA and mRNA differences in the skeletal muscle of Wuranke sheep at different developmental stages and revealed that a series of candidate miRNA-mRNA pairs may act as modulators of muscle development. These results will contribute to future studies on the function of miRNAs and their target mRNAs during skeletal muscle development in Wuranke sheep.

A 15-Gene-Based Risk Signature for Predicting Overall Survival in SCLC Patients Who Have Undergone Surgical Resection.

Small cell lung cancer (SCLC) is a malignancy with a poor prognosis whose treatment has not progressed for decades. The survival benefit of surgery and the selection of surgical candidates are still controversial in SCLC. This study is the first report to identify transcriptomic alterations associated with prognosis and propose a gene expression-based risk signature that can be used to predict overall survival (OS) in SCLC patients who have undergone potentially curative surgery. An integrative transcriptome analysis of three gene expression datasets (GSE30219, GSE43346, and GSE149507) revealed 1734 up-regulated and 2907 down-regulated genes. Cox-Mantel test, Cox regression, and Lasso regression analyses were used to identify genes to be included in the risk signature. EGAD00001001244 and GSE60052-cohorts were used for internal and external validation, respectively. Overall survival was significantly poorer in patients with high-risk scores compared to the low-risk group. The discriminatory performance of the risk signature was superior to other parameters. Multivariate analysis showed that the risk signature has the potential to be an independent predictor of prognosis. The prognostic genes were enriched in pathways including regulation of transcription, cell cycle, cell metabolism, and angiogenesis. Determining the roles of the identified prognostic genes in the pathogenesis of SCLC may contribute to the development of new treatment strategies. The risk signature needs to be validated in a larger cohort of patients to test its usefulness in clinical decision-making.

Integrative analysis of transcriptome and metabolome provide new insights into mechanisms of Capilliposide A against cisplatin-induced nephrotoxicity

ObjectiveCisplatin (CDDP) has been widely used for chemotherapy against tumours. However，the nephrotoxicity has limited its clinical use. Here, we reported a novel compound, Capilliposide A (CPS-A), to exhibit therapeutic effects on CDDP-induced acute kidney injury (AKI) and explored its potential mechanisms via transcriptome and metabolome. Materials and methodsHK-2 cells were treated with CPS-A, after which cell viability, apoptosis and inflammation were investigated. A mouse model of AKI was constructed by single injection of CDDP in vivo. The renal function and morphology and mitochondrial function were assessed by pathological section and transmission electron microscope (TEM). Transcriptomics and metabolomics are used to explore possible mechanisms which was later verified in vitro. ResultsCPS-A administration improved the survival rates of HK-2 cells with a significant decrease in the expression of KIM-1, NGAL, IL-6, IL-8 and IL-1β. In vivo results also suggested that CPS-A attenuates CDDP-induced kidney injury by reducing serum creatinine (Cr) and blood urea nitrogen (BUN) levels. Furthermore, TEM also showed the improvement of mitochondrial ultrastructure both in vivo and vitro. Transcriptomics analysis of the mice's renal cortex indicated the expression of ATF4 and CHOP were upregulated, which was further validated by qPCR and Western blotting in vitro. Integrative analysis of transcriptome and metabolome indicated that L-Leucine enriched in Valine, leucine and isoleucine degradation might be potential targets. ConclusionsCPS-A can effectively regulate endogenous metabolites associated with amino acid metabolism and ameliorate apoptosis and oxidative stress in CDDP-induced AKI by reducing endoplasmic reticulum stress.