Comparative transcriptomics reveals candidate transcription factors involved in costunolide biosynthesis in medicinal plant-Saussurea lappa

Effect of strontium on transcription factors identified by transcriptome analyses of bovine ruminal epithelial cells

BACKGROUND

Strontium (Sr) has similar physicochemical properties as calcium (Ca) and is often used to evaluate the absorption of this mineral. Because the major route of Ca absorption in the bovine occurs in the rumen, it is essential to understand whether Sr impacts the ruminal epithelial cells and to what extent.

RESULTS

In the present study, RNA sequencing and assembled transcriptome assembly were used to identify transcription factors (TFs), screening and bioinformatics analysis in bovine ruminal epithelial cells treated with Sr. A total of 1405 TFs were identified and classified into 64 families based on an alignment of conserved domains. A total of 174 differently expressed TFs (DE-TFs) were increased and 52 DE-TFs were decreased; the biological process-epithelial cell differentiation was inhibited according to the GSEA-GO analysis of TFs; The GO analysis of DE-TFs was enriched in the DNA binding. Protein-protein interaction network (PPI) found 12 hubs, including SMAD4, SMAD2, SMAD3, SP1, GATA2, NR3C1, PPARG, FOXO1, MEF2A, NCOA2, LEF1, and ETS1, which verified genes expression levels by real-time PCR.

CONCLUSIONS

In this study, SMAD2, PPARG, LEF1, ETS1, GATA2, MEF2A, and NCOA2 are potential candidates that could be targeted by Sr to mediate cell proliferation and differentiation, as well as lipid metabolism. Hence, these results enhance the comprehension of Sr in the regulation of transcription factors and provide new insight into the study of Sr biological function in ruminant animals.

Integrated analysis of smRNAome, transcriptome, and degradome data to decipher microRNAs regulating costunolide biosynthesis in Saussurea lappa

Saussurea lappa (S. lappa) has been known to synthesize medicinally important, costunolide. Due to its immense therapeutic importance, understanding of regulatory mechanism associated with its biosynthesis is crucial. The identification of genes and transcription factors (TFs) in S. lappa, created a clear picture of costunolide biosynthesis pathways. Further to understand the regulation of costunolide biosynthesis by miRNAs, an integrated study of transcriptome, miRNAs, and degradome was performed. Identified candidate miRNAs and associated feed-forward loops (FFLs) illustrates their regulatory role in secondary metabolite biosynthesis. Small RNA and degradome sequencing were performed for leaf and root tissues to determine miRNAs-targets pairs. A total of 711 and 525 such targets were obtained for novel and known miRNAs respectively. This data was used to generate costunolide-specific miRNA-TF-gene interactome to perform systematic analyses through graph theoretical approach. Interestingly, miR171c.1 and sla-miR121 were identified as key regulators to connect and co-regulate both mevalonate and sesquiterpenoid pathways to bio-synthesize costunolide. Tissue-specific FFLs were identified to be involved in costunolide biosynthesis which further suggests the evolutionary co-relation of root-specific networks in synthesis of secondary metabolites in addition to leaf-specific networks. This integrative approach allowed us to determine candidate miRNAs and associated tissue-specific motifs involved in the diversification of secondary metabolites. MiRNAs identified in present study can provide alternatives for bioengineering tool to enhance the synthesis of costunolide and other secondary metabolites in S. lappa.

Metabolomics-based exploration the response mechanisms of Saussurea involucrata leaves under different levels of low temperature stress

BACKGROUND

Saussurea involucrata (Sik.) is alpine plant that have developed special adaptive mechanisms to resist adverse environmental conditions such as low temperature chilling during long-term adaptation and evolution. Exploring the changes of its metabolites under different temperature stresses is helpful to gain insight into its cold stress tolerance.

METHODS

Ultra-performance liquid chromatography and tandem mass spectrometry were used to analyze the metabolites in the leaves of Sik. under low different temperature stress conditions.

RESULTS

A total of 753 metabolites were identified, and 360 different metabolites were identified according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) involved in the biosynthesis of secondary metabolites and amino acids and sugars. Sucrose and trehalose synthesis, glycolysis, tricarboxylic acid cycle, pentose phosphate pathway, glutamic acid-mediated proline biosynthesis, purine metabolism, amino acid metabolism, phenylpropane synthesis pathway metabolites all respond to low temperature stress. Under cold stress conditions, carbohydrates in Sik. leaves accumulate first than under freezing conditions, and the lower the temperature under freezing conditions, the less amino acids accumulate, while the phenolic substances increase. The expression of various substances in LPE and LPC increased more than 10-fold after low temperature stress compared with the control, but the content of LPE and LPC substances decreased after cold adaptation. In addition, purines and phenolics decreased and amino acids accumulated significantly under freezing conditions.

CONCLUSION

The metabolic network of Sik. leaves under different low temperature stress conditions was proposed, which provided a reference for further exploration of the metabolic mechanism related to low temperature stress tolerance of Sik.

Transcriptome-wide analysis reveals core transcriptional regulators associated with culm development and variation in Dendrocalamus sinicus, the strongest woody bamboo in the world

Transcription factors (TFs) play indispensable roles in plant development and stress responses. As the largest woody bamboo species in the world, Dendrocalamus sinicus is endemic to Yunnan Province, China, and possesses two natural variants characterized by culm shape, namely straight or bent culms. Understanding the transcriptional regulation network of D. sinicus provides a unique opportunity to clarify the growth and development characteristics of woody bamboos. In this study, 10,236 TF transcripts belonging to 57 families were identified from transcriptome data of two variants at different developmental stages, from which we constructed a transcriptional regulatory network and unigene-coding protein-TFs interactive network of culm development for this attractive species. Gene function enrichment analysis revealed that hormone signaling and MAPK signaling pathways were two most enriched pathways in TF-regulated network. Based on PPI analysis, 50 genes interacting with nine TFs were screened as the core regulation components related to culm development. Of them, 18 synergistic genes of seven TFs, including nuclear cap-binding protein subunit 1, transcription factor GTE9-like, and ATP-dependent DNA helicase DDX11 isoform X1, involved in culm-shape variation. Most of these genes would interact with MYB, C3H, and ARF transcription factors. Six members with two each from ARF, C3H, and MYB transcription factor families and six key interacting genes (IAA3, IAA19, leucine-tRNA ligase, nuclear cap-binding protein subunit 1, elongation factor 2, and coiled-coil domain-containing protein 94) cooperate with these transcription factors were differentially expressed at development stage of young culms, and were validated by quantitative PCR. Our results represent a crucial step towards understanding the regulatory mechanisms of TFs involved in culm development and variation of D. sinicus.

An integrated network analysis approach to identify potential key genes, transcription factors, and microRNAs regulating human hematopoietic stem cell aging

Hematopoietic stem cells (HSCs) undergo functional deterioration with increasing age that causes loss of their self-renewal and regenerative potential. Despite various efforts, significant success in identifying molecular regulators of HSC aging has not been achieved, one prime reason being the non-availability of appropriate human HSC samples. To demonstrate the scope of integrating and re-analyzing the HSC transcriptomics data available, we used existing tools and databases to structure a sequential data analysis pipeline to predict potential candidate genes, transcription factors, and microRNAs simultaneously. This sequential approach comprises (i) collecting matched young and aged mice HSC sample datasets, (ii) identifying differentially expressed genes, (iii) identifying human homologs of differentially expressed genes, (iv) inferring gene co-expression network modules, and (v) inferring the microRNA-transcription factor-gene regulatory network. Systems-level analyses of HSC interaction networks provided various insights based on which several candidates were predicted. For example, 16 HSC aging-related candidate genes were predicted (e.g., CD38, BRCA1, AGTR1, GSTM1, etc.) from GCN analysis. Following this, the shortest path distance-based analyses of the regulatory network predicted several novel candidate miRNAs and TFs. Among these, miR-124-3p was a common regulator in candidate gene modules, while TFs MYC and SP1 were identified to regulate various candidate genes. Based on the regulatory interactions among candidate genes, TFs, and miRNAs, a potential regulation model of biological processes in each of the candidate modules was predicted, which provided systems-level insights into the molecular complexity of each module to regulate HSC aging.

Bioinformatic identification of hub genes and related transcription factors in low shear stress treated endothelial cells

BACKGROUND

Recent evidences indicated that shear stress is critical in orchestrating gene expression in cardiovascular disease. It is necessary to identify the mechanism of shear stress influencing gene expression in physiology and pathophysiology conditions. This paper aimed to identify candidate hub genes and its transcription factors with bioinformatics.

METHODS

We analyzed microarray expression profile of GSE16706 to identify differentially expressed genes (DEGs) in low shear stress (1 dyne/cm2) treated human umbilical vein endothelial cells (HUVECs) compared with static condition for 24 h.

RESULTS

652 DEGs, including 333 up-regulated and 319 down-regulated DEGs, were screen out. Functional enrichment analysis indicated enrichment items mainly included cytokine-cytokine receptor interaction and cell cycle. Five hub genes (CDC20, CCNA2, KIF11, KIF2C and PLK1) and one significant module (score = 17.39) were identified through protein-protein interaction (PPI) analysis. Key transcriptional factor FOXC1 displayed close interaction with all the hub genes via gene-transcriptional factor network. Single-gene GSEA analysis indicated that CDC20 was linked to the G2M_CHECKPOINT pathway and cell cycle pathway.

CONCLUSIONS

By using integrated bioinformatic analysis, a new transcriptional factor and hub-genes network related to HUVECs treated with low shear stress were identified. The new regulation mechanism we discovered may be a promising potential therapeutic target for cardiovascular disease.

Identification and characterization of SlbHLH, SlDof and SlWRKY transcription factors interacting with SlDPD gene involved in costunolide biosynthesis in Saussurea lappa

The genes involved in costunolide biosynthesis in Saussurea lappa have been identified recently by our lab. However, the study of transcriptional regulators of these genes was lacking for better opportunities for engineering the pharmacologically important biosynthetic pathway. Therefore, we cloned the promoter region of diphosphomevalonate decarboxylase gene (DPD) and analyzed its cis-acting regulatory elements to reveal the potential transcription factor (TF) binding sites for Dof, bHLH and WRKY family proteins in the gene promoter. The transcriptome study approach followed by the hidden Markov model based search, digital gene expression, co-expression network analysis, conserved domain properties and evolutionary analyses were carried out to screen out seven putative TFs for the DPD-TF interaction studies. Yeast one-hybrid assays were performed and three TFs were reported, namely, SlDOF2, SlbHLH3 and SlWRKY2 from Dof, bHLH and WRKY families, respectively that interacted positively with the DPD gene of the costunolide biosynthetic pathway. The tissue specific relative gene expression studies also supported the linked co-expression of the gene and its interacting TFs The present report will improve the understanding of transcriptional regulation pattern of costunolide biosynthetic pathway.

MYB Transcription Factors as Regulators of Secondary Metabolism in Plants

MYB transcription factors (TFs), as one of the largest gene families in plants, play important roles in multiple biological processes, such as plant growth and development, cell morphology and pattern building, physiological activity metabolism, primary and secondary metabolic reactions, and responses to environmental stresses. The function of MYB TFs in crops has been widely studied, but few studies have been done on medicinal plants. In this review, we summarized the MYB TFs that play important roles in secondary metabolism and emphasized the possible mechanisms underlying how MYB TFs are regulated at the protein, posttranscriptional, and transcriptional levels, as well as how they regulate the downstream target gene networks related to secondary metabolism in plants, especially in medicinal plants.