scholarly journals Transcriptome sequencing and analysis reveals the molecular response to selenium stimuli in Pueraria lobata (willd.) Ohwi

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e8768 ◽  
Author(s):  
Kunyuan Guo ◽  
Yiwei Yao ◽  
Meng Yang ◽  
Yanni Li ◽  
Bin Wu ◽  
...  

Pueraria lobata (willd.) Ohwi is a consumable selenium-enriched plant used for medicinal purposes. The molecular response to selenium (Se) stimuli in P. lobata is currently unknown. We used RNA-Seq to identify potential genes involved in selenite metabolism and analyzed their expression profiles. We obtained a total of 150,567 unigenes, of which 90,961 were annotated, including 16 structural genes, 14 sulfate transporters, and 13 phosphate transporters that may be involved in Se metabolism, and 33 candidate structural genes involved in isoflavone biosynthesis. The genes with a —foldchange— >2 and q value <0.05 after sodium selenite treatment were identified as differentially expressed genes (DEGs). We obtained a total of 4,246 DEGs, which were enriched in GO terms that included “response to stimulus”, “response to stress”, “signal transduction”, “response to abiotic stimulus”, and “response to chemical”. Of the 4,246 DEGs, one sulfate transporter and five phosphate transporter genes involved Se metabolism, and nine structural genes involved in isoflavone biosynthesis were up-regulated. The expression patterns of 10 DEGs were selected randomly and validated using qRT-PCR. The Pearson Correlation Coefficient (r) was 0.86, indicating the reliability of RNA-Seq results. 22 Reactive Oxygen Species (ROS) scavenging DEGs were found, 11 of which were up-regulated. 436, 624 transcription factors (TFs) correlated with structural genes were identified that may be involved in Se and isoflavone biosynthesis, respectively, using r (r > 0.7 or r <  − 0.7). 556 TFs were related to at least one sulfate and phosphate transporter. Our results provided a comprehensive description of gene expression and regulation in response to Se stimuli in P. lobata.

Genes ◽  
2021 ◽  
Vol 12 (10) ◽  
pp. 1610
Author(s):  
Mohammad Vatanparast ◽  
Youngjin Park

Solenopsis japonica, as a fire ant species, shows some predatory behavior towards earthworms and woodlice, and preys on the larvae of other ant species by tunneling into a neighboring colony’s brood chamber. This study focused on the molecular response process and gene expression profiles of S. japonica to low (9 °C)-temperature stress in comparison with normal temperature (25 °C) conditions. A total of 89,657 unigenes (the clustered non-redundant transcripts that are filtered from the longest assembled contigs) were obtained, of which 32,782 were annotated in the NR (nonredundant protein) database with gene ontology (GO) terms, gene descriptions, and metabolic pathways. The results were 81 GO subgroups and 18 EggNOG (evolutionary genealogy of genes: Non-supervised Orthologous Groups) keywords. Differentially expressed genes (DEGs) with log2fold change (FC) > 1 and log2FC < −1 with p-value ≤ 0.05 were screened for cold stress temperature. We found 215 unigenes up-regulated and 115 unigenes down-regulated. Comparing transcriptome profiles for differential gene expression resulted in various DE proteins and genes, including fatty acid synthases and lipid metabolism, which have previously been reported to be involved in cold resistance. We verified the RNA-seq data by qPCR on 20 up- and down-regulated DEGs. These findings facilitate the basis for the future understanding of the adaptation mechanisms of S. japonica and the molecular mechanisms underlying the response to low temperatures.


BMC Genomics ◽  
2019 ◽  
Vol 20 (1) ◽  
Author(s):  
Meng-qi Wang ◽  
Zhen-yu Bai ◽  
Ya-fang Xiao ◽  
Yan Li ◽  
Qing-lin Liu ◽  
...  

Abstract Background Cadmium (Cd) is a serious heavy metal (HM) soil pollutant. To alleviate or even eliminate HM pollution in soil, environmental-friendly methods are applied. One is that special plants are cultivated to absorb the HM in the contaminated soil. As an excellent economical plant with ornamental value and sound adaptability, V. bonariensis could be adapted to this very situation. In our study, the Cd tolerance in V. bonariensis was analyzed as well as an overall analysis of transcriptome. Results In this study, the tolerance of V. bonariensis to Cd stress was investigated in four aspects: germination, development, physiological changes, and molecular alterations. The results showed that as a non-hyperaccumulator, V. bonariensis did possess the Cd tolerance and the capability to concentration Cd. Under Cd stress, all 237, 866 transcripts and 191, 370 unigenes were constructed in the transcriptome data of V. bonariensis roots. The enrichment analysis of gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway revealed that differentially expressed genes (DEGs) under Cd stress were predominately related to cell structure, reactive oxygen species (ROS) scavenging system, chelating reaction and secondary metabolites, transpiration and photosynthesis. DEGs encoding lignin synthesis, chalcone synthase (CHS) and anthocyanidin synthase (ANS) were prominent in V. bonariensis under Cd stress. The expression patterns of 10 DEGs, validated by quantitative real-time polymerase chain reaction (qRT-PCR), were in highly accordance with the RNA-Sequence (RNA-Seq) results. The novel strategies brought by our study was not only benefit for further studies on the tolerance of Cd and functional genomics in V. bonariensis, but also for the improvement molecular breeding and phytoremediation.


Blood ◽  
2006 ◽  
Vol 108 (11) ◽  
pp. 2355-2355
Author(s):  
Steven M. Kornblau ◽  
David McCue ◽  
Kang L. Lu ◽  
Wenjing Chen ◽  
Kevin R. Coombes

Abstract Protein expression and activation determines the pathophysiology of leukemic cells in Myelodysplasia (MDS) and Acute Myelogenous Leukemia (AML) and is dependent on endogenous changes (e.g mutation, methylation) and exogenous signals from stromal interactions, cytokines (CTKN) and chemokines. We have previously performed proteomics on primary AML sample (using reverse phase protein arrays) and wanted to assess how cytokines affect protein expression and phosphorylation. Prior studies of CTKN expression in AML and MDS have generally measured individual CTKNs, but not provided an overall assessment of CTKN expression. We measured the level of 26 CTKN (IL-1RA, 1B, 2, 4 5, 6, 7 , 8 , 9, 10,12, 13, 15, 17, Eotaxin, FGFB, G-CSF, GM-CSF, IFNγ, IP10, MCP1, MIP1α, MIP1β, PDGF, TNFα and VEGF) using multiplex cytometry (Bioplex™, Biorad) in serum samples from 176 AML (138 untreated (New), 38 relapsed (REL)) and 114 MDS patients (97 New, 10 post biological therapy, 7 REL) and 19 normal (NL) subjects. Individual CTKN expression was not correlated with clinical features (e.g. age, gender, cytogenetics, FAB, HB, WBC, platelet etc). The levels of IL -1β, 4, 5, 6, 7,10,12, 13, 17, IFNγ, FGFB and MIP1α were significantly lower and IL-8 and 15 higher in AML/MDS compared to NL. The expression profiles of AML and MDS were statistically indistinguishable whether analyzed individually or by unsupervised hierarchical clustering, except for IL-8 and 13 (higher in AML) and VEGF (higher in MDS). When CTKN were evaluated individually in new AML cases higher levels of IL4, 5 and 10 correlated significantly with remission attainment, and higher levels of IL8, Il1Ra, IP-10, Mip1β, VEGF and ILR, correlated significantly with shorter survival. No CTKN predicted remission attainment or survival in MDS. Unsupervised hierarchical bootstrap clustering using Pearson correlation and average linkage of CTKN expression relative to other CTKN expression, where high levels of one CTKN correlated with high expression of the other, revealed 6 highly reproducible expression patterns: IL-1β 4, 7, 10, 12, 13, G-CSF, IFNγ, MIP1α and PDGF IL 1ra, 6, 8 Eotaxin, IP-10, MIP1β and VEGF, IL2, 9, 15 and GMCSF, IL5 IL-7, FGF-Basic, TNFα and MCP1. Similar unsupervised clustering of the samples based on CTKN expression using average linkage also revealed 5 disease clusters and a NL sample cluster (containing all 19 NL samples). Average expression levels of each CTKN in these 5 clusters show diminished expression of most CTKN that had high expression in the NL samples, with each group showing increase in expression in a distinct subset of CTKN relative to NL. Remission attainment was strongly associated with cytokine signature (P=0.005). Additional CTKN are being studied (SCF, TGFβ, IL3). Comparison of CTKN expression patterns with proteomic profiling of expression and phosphorylation status is pending. In summary, this is the largest sample set studied for multiple CTKN expression in AML and MDS and the first assessment of many of these CTKN in these diseases. Most CTKNs showed different expression in AML and MDS compared to NL. Interestingly, CTKN expression in AML and MDS were similar. Many CTKN are predictive of outcome individually. CTKN signatures distinguish groups of patients and are predictive of outcome. Correlation with proteomic profiling may suggest CTKN to target in combination with other targeted therapies to maximally affect activated pathways.


2021 ◽  
Author(s):  
Jakub Jankowski ◽  
Hye Kyung Lee ◽  
Julia Wilflingseder ◽  
Lothar Hennighausen

SummaryRecently, a short, interferon-inducible isoform of Angiotensin-Converting Enzyme 2 (ACE2), dACE2 was identified. ACE2 is a SARS-Cov-2 receptor and changes in its renal expression have been linked to several human nephropathies. These changes were never analyzed in context of dACE2, as its expression was not investigated in the kidney. We used Human Primary Proximal Tubule (HPPT) cells to show genome-wide gene expression patterns after cytokine stimulation, with emphasis on the ACE2/dACE2 locus. Putative regulatory elements controlling dACE2 expression were identified using ChIP-seq and RNA-seq. qRT-PCR differentiating between ACE2 and dACE2 revealed 300- and 600-fold upregulation of dACE2 by IFNα and IFNβ, respectively, while full length ACE2 expression was almost unchanged. JAK inhibitor ruxolitinib ablated STAT1 and dACE2 expression after interferon treatment. Finally, with RNA-seq, we identified a set of genes, largely immune-related, induced by cytokine treatment. These gene expression profiles provide new insights into cytokine response of proximal tubule cells.


2021 ◽  
Vol 12 ◽  
Author(s):  
Yani Dong ◽  
Likang Lyu ◽  
Haishen Wen ◽  
Bao Shi

Long noncoding RNAs (lncRNAs) have been identified to be involved in half-smooth tongue sole (Cynoglossus semilaevis) reproduction. However, studies of their roles in reproduction have focused mainly on the ovary, and their expression patterns and potential roles in the brain and pituitary are unclear. Thus, to explore the mRNAs and lncRNAs that are closely associated with reproduction in the brain and pituitary, we collected tongue sole brain and pituitary tissues at three stages for RNA sequencing (RNA-seq), the 5,135 and 5,630 differentially expressed (DE) mRNAs and 378 and 532 DE lncRNAs were identified in the brain and pituitary, respectively. The RNA-seq results were verified by RT-qPCR. Moreover, enrichment analyses were performed to analyze the functions of DE mRNAs and lncRNAs. Interestingly, their involvement in pathways related to metabolism, signal transduction and endocrine signaling was revealed. LncRNA-target gene interaction networks were constructed based on antisense, cis and trans regulatory mechanisms. Moreover, we constructed competing endogenous RNA (ceRNA) networks. In summary, this study provides mRNA and lncRNA expression profiles in the brain and pituitary to understand the molecular mechanisms regulating tongue sole reproduction.


2021 ◽  
Vol 12 ◽  
Author(s):  
Pengpeng Zhang ◽  
Mingxuan Sheng ◽  
Chunyu Du ◽  
Zhe Chao ◽  
Haixia Xu ◽  
...  

Brown adipose tissue (BAT) is specialized for energy expenditure, thus a better understanding of the regulators influencing BAT development could provide novel strategies to defense obesity. Many protein-coding genes, miRNAs, and lncRNAs have been investigated in BAT development, however, the expression patterns and functions of circRNA in brown adipogenesis have not been reported yet. This study determined the circRNA expression profiles across brown adipogenesis (proliferation, early differentiated, and fully differentiated stages) by RNA-seq. We identified 3,869 circRNAs and 36.9% of them were novel. We found the biogenesis of circRNA was significantly related to linear mRNA transcription, meanwhile, almost 70% of circRNAs were generated by alternative back-splicing. Next, we examined the cell-specific and differentiation stage-specific expression of circRNAs. Compared to white adipocytes, nearly 30% of them were specifically expressed in brown adipocytes. Further, time-series expression analysis showed circRNAs were dynamically expressed, and 117 differential expression circRNAs (DECs) in brown adipogenesis were identified, with 77 upregulated and 40 downregulated. Experimental validation showed the identified circRNAs could be successfully amplified and the expression levels detected by RNA-seq were reliable. For the potential functions of the circRNAs, GO analysis suggested that the decreased circRNAs were enriched in cell proliferation terms, while the increased circRNAs were enriched in development and thermogenic terms. Bioinformatics predictions showed that DECs contained numerous binding sites of functional miRNAs. More interestingly, most of the circRNAs contained multiple binding sites for the same miRNA, indicating that they may facilitate functions by acting as microRNA sponges. Collectively, we characterized the circRNA expression profiles during brown adipogenesis and provide numerous novel circRNAs candidates for future brown adipogenesis regulating studies.


2020 ◽  
Author(s):  
Dan Huang ◽  
Jian Liu ◽  
Lei Wan ◽  
Yanyan Fang ◽  
Yan Long ◽  
...  

Abstract Background Ankylosing spondylitis (AS) is a chronic autoimmune disease affecting the sacroiliac joint. To date, few studies have examined the association between long non-coding RNAs (lncRNAs) and AS pathogenesis. As such, we herein sought to characterize patterns of AS-related lncRNA expression and to evaluate the potential role played by these lncRNAs in this complex autoimmune context. Methods We conducted an RNA-seq analysis of peripheral blood mononuclear cell samples isolated from five AS patients and corresponding controls. These data were then leveraged to characterize AS-related lncRNA expression patterns. We further conducted GO and KEGG enrichment analyses of the parental genes encoding these lncRNAs, and we confirmed the validity of our RNA-seq data by assessing the expression of six lncRNAs via qRT-PCR in 15 AS and control patient samples. Pearson correlation analyses were additionally employed to examine the associations between the expression levels of these six lncRNAs and patient clinical index values. Results We detected 56575 total lncRNAs in AS and control patient samples during our initial RNA-seq analysis, of which 200 and 70 were found to be up- and down-regulated (FC > 2 or < 0.05; P < 0.05), respectively, in AS samples relative to controls. In qRT-PCR validation assays, we confirmed the significant upregulation of NONHSAT118801.2, ENST00000444046, and NONHSAT183847.1 and the significant downregulation of NONHSAT205110.1, NONHSAT105444.2, and NONHSAT051856.2 in AS patient samples. We further found the expression of NONHSAT118801.2 and NONHSAT183847.1 to be positively correlated with disease severity. Conclusion Overall, our findings highlight several lncRNAs that are specifically expressed in the context of AS, indicating that they may play key functions in the pathogenesis of this autoimmune disease. Specifically, we determined that NONHSAT118801.2 and NONHSAT183847.1 may be valuable biomarkers of AS.


2021 ◽  
Vol 12 ◽  
Author(s):  
Alexandre Hild Aono ◽  
Ricardo José Gonzaga Pimenta ◽  
Ana Letycia Basso Garcia ◽  
Fernando Henrique Correr ◽  
Guilherme Kenichi Hosaka ◽  
...  

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and the core regulator of cellular signaling. Despite this substantial importance, the kinomes of sugarcane and sorghum have not been profiled. Here, we identified and profiled the complete kinomes of the polyploid Saccharum spontaneum (Ssp) and Sorghum bicolor (Sbi), a close diploid relative. The Sbi kinome was composed of 1,210 PKs; for Ssp, we identified 2,919 PKs when disregarding duplications and allelic copies, and these were related to 1,345 representative gene models. The Ssp and Sbi PKs were grouped into 20 groups and 120 subfamilies and exhibited high compositional similarities and evolutionary divergences. By utilizing the collinearity between the species, this study offers insights into Sbi and Ssp speciation, PK differentiation and selection. We assessed the PK subfamily expression profiles via RNA-Seq and identified significant similarities between Sbi and Ssp. Moreover, coexpression networks allowed inference of a core structure of kinase interactions with specific key elements. This study provides the first categorization of the allelic specificity of a kinome and offers a wide reservoir of molecular and genetic information, thereby enhancing the understanding of Sbi and Ssp PK evolutionary history.


Blood ◽  
2012 ◽  
Vol 120 (21) ◽  
pp. 3298-3298 ◽  
Author(s):  
Eric R. Londin ◽  
Eleftheria Hatzimichael ◽  
Phillipe Loher ◽  
Yue Zhao ◽  
Yi Jing ◽  
...  

Abstract Abstract 3298 The anucleate platelets play a critical role in the formation of thrombi and prevention of bleeding. While the repertoire of platelet transcripts is a reflection of the megakaryocyte at the time of platelet differentiation, post-transcriptional events are known to occur. Furthermore, a strong correlation between the expressed mRNAs and proteome has been identified. Having a complete understanding of the platelet transcriptome is important for generating insights into the genetic basis of platelet disease traits. To capture the complexity of the platelet transcriptome, we performed RNA sequencing (RNA-seq) in leukocyte-depleted platelets from 10 males, with median age of 24.5 yrs and unremarkable medical history. Their short and long RNA platelet transcriptomes were analyzed on the SOLiD 5500xl sequencing platform. We generated ∼3.5 billion sequence reads ∼40% of which could be mapped uniquely to the human genome. Our analysis revealed that ∼9,000 distinct protein-coding mRNAs and ∼800 microRNAs (miRNAs) were present in the transcriptome of each of the 10 sequenced individuals. Comparison of the levels of mRNA expression across the 10 individuals showed an exceptional level of consistency with pair-wise Pearson correlation values ≥0.98. The miRNA expression profiles across the 10 individuals showed a similar consistency with pair-wise Pearson correlation values ≥0.98. Surprisingly, we found that these mRNAs and miRNAs accounted for a little over 1/2 of all of the uniquely mapped sequence reads suggesting the abundant presence of additional non-protein coding RNA (ncRNA) transcripts. Using the annotated entries of the latest release of the ENSEMBL database, we investigated the genetic make-up of these other transcripts. We found that ∼25% of each individual's uniquely mapped reads corresponded to non-protein coding transcripts from mRNA-coding loci. These reads accounted for more than 10,000 distinct such transcripts. In addition, each of the individuals in our cohort expressed an average of ∼1,500 pseudogenes and ∼200 long intergenic non-coding RNAs (lincRNAs). The short RNA profiles of the ten individuals revealed an abundance of diverse categories of ncRNAs including the signal recognition particle RNA (srpRNA), small nuclear RNA (snRNA) and small cytoplasmic RNAs (scRNA). These ncRNAs are involved in the processing of pre-mRNAs and their presence and prevalence in the anucleate platetet suggests the existence of a complex network of mRNA processing that persists after the megakaryocyte fragmentation. We also investigated the RNA-omes of the ten individuals for evidence of transcription of the pyknon category of ncRNAs. Pyknons are of particular interest because each has numerous intergenic and intronic copies whereas nearly all known human protein-coding genes contain one or more pyknons in their mRNA. Recent experimental work has shown that intergenic instances of the pyknons are transcribed in a tissue- and cell-state specific manner. An average of ∼100,000 pyknons are transcribed in each of the 10 sequenced individuals suggesting the possibility of a far-reaching network of interactions that link exonic space to distant non-exonic regions and are active in platelets. Lastly, we found that a large variety of distinct repeat element categories are expressed in the RNA-omes (both short and long) of these individuals. Among the most abundantly represented categories of repeat elements were DNA transposons, long terminal repeat (LTR) retrotransposons, and non-LTR retrotransposons such as long interspersed elements (LINEs) and short interspersed elements (SINEs). In summary, our RNA-seq analyses have revealed a spectrum of platelet transcripts that transcends protein-coding genes and miRNAs. Indeed, the transcripts that have their source in genomic features not previously discussed or analyzed in the platelet context represent a very significant portion of all platelet transcripts. This in turn suggests an unanticipated richness, and presumably commensurate complexity, for the platelet transcriptome. While the role of these novel non-protein coding RNAs is currently unknown it is expected that at least some of them may be of functional significance which will in turn permit a better understanding of the molecular mechanisms that regulate platelet physiology and may contribute to processes beyond thrombosis and hemostasis. Disclosures: No relevant conflicts of interest to declare.


2019 ◽  
Author(s):  
Xiaomu Liu ◽  
Wei You ◽  
Xianglun Zhang ◽  
Qing Jin ◽  
Xiuwen Tan ◽  
...  

Abstract Background:Previously, we found that mevalonic acid stimulates HMGR expression in bovine intramuscular adipocytes, and influences adipocyte differentiation. However, it remains unclear whether there is any direct link between HMGR, steroidogenic genes, and cholesterol content. RNA-Seq was conducted to determine the differences between the gene expression profiles of bovine adipocytes containing different HMGR expression constructs. Results:In total, 10 234 differentially expressed genes (DEGs) were found. Of these, 35 and 6 DEGs between the control and the overexpression groups were functionally related to lipid and energy metabolisms, respectively. Additionally, 43 and 8 DEGs between the control and the HMGR inhibition groups were related to lipid and energy metabolism. Additionally, several DEGs related to lipid and energy metabolism were identified between the HMGR overexpression group and the HMGR interference group. Several DEGs correlated positively or negatively with overexpression or inhibition of HMGR. We also found that, following activation or inhibition of the HMGR gene, AMPK and SIRT1 had opposite expression patterns in bovine intramuscular adipocytes. Interestingly, the HMGR gene was downregulated when HMGR was overexpressed, and upregulated when HMGR was inhibited. Conclusion:Our findings establish a theoretical understanding of signaling pathways involved in cholesterol synthesis by elucidating the relationships between key genes.


Sign in / Sign up

Export Citation Format

Share Document