scholarly journals Transcriptomic Analysis Reveals Hormonal Control of Shoot Branching in Salix matsudana

Forests ◽  
2020 ◽  
Vol 11 (3) ◽  
pp. 287
Author(s):  
Juanjuan Liu ◽  
Bingbing Ni ◽  
Yanfei Zeng ◽  
Caiyun He ◽  
Jianguo Zhang

Shoot branching is regulated by axillary bud activities, which subsequently grow into branches. Phytohormones play a central role in shoot branching control, particularly with regard to auxin, cytokinins (CKs), strigolactones (SLs), and gibberellins (GAs). To further study the molecular basis for the shoot branching in Salix matsudana, how shoot branching responds to hormones and regulatory pathways was investigated, and potential genes involved in the regulation of shoot branching were identified. However, how these positive and inhibitory processes work on the molecular level remains unknown. RNA-Seq transcriptome expression analysis was used to elucidate the mechanisms underlying shoot branching. In total, 102 genes related to auxin, CKs, SLs, and GAs were differentially expressed in willow development. A majority of the potential genes associated with branching were differentially expressed at the time of shoot branching in S. matsudana, which have more number of branching. These findings are consistent with the growth and physiological results. A regulatory network model was proposed to explain the interaction between the four hormones that control shoot branching. Collectively, the results presented here contribute to a more comprehensive understanding of the hormonal effects on shoot branching in S. matsudana. In the future, these findings will help uncover the interactions among auxin, SLs, CKs, and GAs that control shoot branching in willow, which could help improve plant structures through the implementation of molecular techniques in targeted breeding.

Genes ◽  
2021 ◽  
Vol 13 (1) ◽  
pp. 21
Author(s):  
Vladimir Babenko ◽  
Olga Redina ◽  
Dmitry Smagin ◽  
Irina Kovalenko ◽  
Anna Galyamina ◽  
...  

Both aggressive and aggression-deprived (AD) species represent pathologic cases intensely addressed in psychiatry and substance abuse disciplines. Previously, we reported that AD mice displayed a higher aggressive behavior score than the aggressive group, implying the manifestation of a withdrawal effect. We employed an animal model of chronic social conflicts, curated in our lab for more than 30 years. In the study, we pursued the task of evaluating key events in the dorsal striatum transcriptome of aggression experienced mice and AD species compared to controls using RNA-Seq profiling. Aggressive species were subjected to repeated social conflict encounters (fights) with regular positive (winners) experience in the course of 20 consecutive days (A20 group). This led to a profoundly shifted transcriptome expression profile relative to the control group, outlined by more than 1000 differentially expressed genes (DEGs). RNA-Seq cluster analysis revealed that elevated cyclic AMP (cAMP) signaling cascade and associated genes comprising 170 differentially expressed genes (DEGs) in aggressive (A20) species were accompanied by a downturn in the majority of other metabolic/signaling gene networks (839 DEGs) via the activation of transcriptional repressor DEGs. Fourteen days of a consecutive fighting deprivation period (AD group) featured the basic restoration of the normal (control) transcriptome expression profile yielding only 62 DEGs against the control. Notably, we observed a network of 12 coordinated DEG Transcription Factor (TF) activators from 62 DEGs in total that were distinctly altered in AD compared to control group, underlining the distinct transcription programs featuring AD group, partly retained from the aggressive encounters and not restored to normal in 14 days. We found circadian clock TFs among them, reported previously as a withdrawal effect factor. We conclude that the aggressive phenotype selection with positive reward effect (winning) manifests an addiction model featuring a distinct opioid-related withdrawal effect in AD group. Along with reporting profound transcriptome alteration in A20 group and gaining some insight on its specifics, we outline specific TF activator gene networks associated with transcriptional repression in affected species compared to controls, outlining Nr1d1 as a primary candidate, thus offering putative therapeutic targets in opioid-induced withdrawal treatment.


Genes ◽  
2020 ◽  
Vol 11 (8) ◽  
pp. 935 ◽  
Author(s):  
Chong Wang ◽  
Miao Zhao ◽  
Wen-Qian Zhang ◽  
Ming-Yu Huang ◽  
Can Zhu ◽  
...  

The mouse is widely used to study decidualization and there are three well-established mouse models of decidualization, namely natural pregnancy decidualization (NPD), artificial decidualization (AD), and in vitro decidualization (IVD). However, the extent of similarity and difference between these models at the molecular level remains largely unknown. Here, we performed a comparative analysis using the RNA-seq approach. In the NPD model, which is thought to be the golden standard of mouse decidualization, we found a total of 5277 differentially expressed genes, with 3158 genes being up-regulated and 2119 genes being down-regulated. A total of 4294 differentially expressed genes were identified in the AD model: 1127 up-regulated genes and 3167 down-regulated genes. In comparison to NPD, 1977 genes were consistently expressed, whereas only 217 genes were inconsistently expressed, indicating that AD is a reliable model for mouse decidualization. In the IVD model, RNA-seq analysis revealed that 513 genes were up-regulated and 988 genes were down-regulated. Compared to NPD, 310 genes were consistently expressed, whereas 456 genes were inconsistently expressed. Moreover, although the decidualization marker Prl8a2 (prolactin family 8 subfamily a member 2) was up-regulated, the widely-used marker Alpl (alkaline phosphatase liver/bone/kidney) was down-regulated in the IVD model. Therefore, we suggest that the IVD model should be optimized to mimic NPD at the transcriptomic level. Our study contributes to an increase in the knowledge about mouse models of decidualization.


2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 12.2-12
Author(s):  
I. Muller ◽  
M. Verhoeven ◽  
H. Gosselt ◽  
M. Lin ◽  
T. De Jong ◽  
...  

Background:Tocilizumab (TCZ) is a monoclonal antibody that binds to the interleukin 6 receptor (IL-6R), inhibiting IL-6R signal transduction to downstream inflammatory mediators. TCZ has shown to be effective as monotherapy in early rheumatoid arthritis (RA) patients (1). However, approximately one third of patients inadequately respond to therapy and the biological mechanisms underlying lack of efficacy for TCZ remain elusive (1). Here we report gene expression differences, in both whole blood and peripheral blood mononuclear cells (PBMC) RNA samples between early RA patients, categorized by clinical TCZ response (reaching DAS28 < 3.2 at 6 months). These findings could lead to identification of predictive biomarkers for TCZ response and improve RA treatment strategies.Objectives:To identify potential baseline gene expression markers for TCZ response in early RA patients using an RNA-sequencing approach.Methods:Two cohorts of RA patients were included and blood was collected at baseline, before initiating TCZ treatment (8 mg/kg every 4 weeks, intravenously). DAS28-ESR scores were calculated at baseline and clinical response to TCZ was defined as DAS28 < 3.2 at 6 months of treatment. In the first cohort (n=21 patients, previously treated with DMARDs), RNA-sequencing (RNA-seq) was performed on baseline whole blood PAXgene RNA (Illumina TruSeq mRNA Stranded) and differential gene expression (DGE) profiles were measured between responders (n=14) and non-responders (n=7). For external replication, in a second cohort (n=95 therapy-naïve patients receiving TCZ monotherapy), RNA-seq was conducted on baseline PBMC RNA (SMARTer Stranded Total RNA-Seq Kit, Takara Bio) from the 2-year, multicenter, double-blind, placebo-controlled, randomized U-Act-Early trial (ClinicalTrials.gov identifier: NCT01034137) and DGE was analyzed between 84 responders and 11 non-responders.Results:Whole blood DGE analysis showed two significantly higher expressed genes in TCZ non-responders (False Discovery Rate, FDR < 0.05): urotensin 2 (UTS2) and caveolin-1 (CAV1). Subsequent analysis of U-Act-Early PBMC DGE showed nine differentially expressed genes (FDR < 0.05) of which expression in clinical TCZ non-responders was significantly higher for eight genes (MTCOP12, ZNF774, UTS2, SLC4A1, FECH, IFIT1B, AHSP, and SPTB) and significantly lower for one gene (TND2P28M). Both analyses were corrected for baseline DAS28-ESR, age and gender. Expression of UTS2, with a proposed function in regulatory T-cells (2), was significantly higher in TCZ non-responders in both cohorts. Furthermore, gene ontology enrichment analysis revealed no distinct gene ontology or IL-6 related pathway(s) that were significantly different between TCZ-responders and non-responders.Conclusion:Several genes are differentially expressed at baseline between responders and non-responders to TCZ therapy at 6 months. Most notably, UTS2 expression is significantly higher in TCZ non-responders in both whole blood as well as PBMC cohorts. UTS2 could be a promising target for further analyses as a potential predictive biomarker for TCZ response in RA patients in combination with clinical parameters (3).References:[1]Bijlsma JWJ, Welsing PMJ, Woodworth TG, et al. Early rheumatoid arthritis treated with tocilizumab, methotrexate, or their combination (U-Act-Early): a multicentre, randomised, double-blind, double-dummy, strategy trial. Lancet. 2016;388(10042):343-55.[2]Bhairavabhotla R, Kim YC, Glass DD, et al. Transcriptome profiling of human FoxP3+ regulatory T cells. Human Immunology. 2016;77(2):201-13.[3]Gosselt HR, Verhoeven MMA, Bulatovic-Calasan M, et al. Complex machine-learning algorithms and multivariable logistic regression on par in the prediction of insufficient clinical response to methotrexate in rheumatoid arthritis. Journal of Personalized Medicine. 2021;11(1).Disclosure of Interests:None declared


2020 ◽  
Vol 79 (Suppl 1) ◽  
pp. 178-179
Author(s):  
S. Alehashemi ◽  
M. Garg ◽  
B. Sellers ◽  
A. De Jesus ◽  
A. Biancotto ◽  
...  

Background:Systemic Autoinflammatory diseases present with sterile inflammation. NOMID (Neonatal-Onset Multisystem Inflammatory Disease) is caused by gain-of-function mutations inNLRP3and excess IL-1 production, presents with fever, neutrophilic dermatosis, aseptic meningitis, hearing loss and eye inflammation; CANDLE (Chronic Atypical Neutrophilic Dermatosis, Lipodystrophy and Elevated Temperature) is caused by loss-of-function mutations in proteasome genes that lead to type-1 interferon signaling, characterized by fever, panniculitis, lipodystrophy, cytopenia, systemic and pulmonary hypertension and basal ganglia calcification. IL-1 blockers are approved for NOMID and JAK-inhibitors show efficacy in CANDLE treatment.Objectives:We used proteomic analysis to compare differentially expressed proteins in active NOMID and CANDLE compared to healthy controls before and after treatment, and whole blood bulk RNA seq to identify the immune cell signatures.Methods:Serum samples from active NOMID (n=12) and CANDLE (n=7) before and after treatment (table 1) and age matched healthy controls (HC) (n=7) were profiled using the SomaLogic platform (n=1125 proteins). Differentially expressed proteins in NOMID and CANDLE were ranked after non-parametric tests for unpaired (NOMIDp<0.05, CANDLE,p<0.1) and paired (p<0.05) analysis and assessed by enriched Gene Ontology pathways and network visualization. Whole blood RNA seq was performed (NOMID=7, CANDLE=7, Controls =5) and RPKM values were used to assess immune cells signatures.Table 1.Patient’s characteristicsNOMIDN=12, Male =6CANDLEN=7, Male =6AgeMedian (range)12 (2, 28)16 (3, 20)Ethnicity%White (Hispanic)80 (20)100 (30)GeneticsNLRP3mutation(2 Somatic, 10 Germline)mutations in proteasome component genes(1 digenic, 6 Homozygous/compound Heterozygous)Before treatmentAfter treatmentBefore treatmentAfter treatmentCRPMedian (range) mg/L52 (16-110)5 (0-23)5 (0-101)1 (0-4)IFN scoremedian (range)0NA328 (211-1135)3 (0-548)Results:Compared to control, 205 proteins (127 upregulated, 78 downregulated) were significantly different at baseline in NOMID, compared to 163 proteins (101 upregulated, and 62 downregulated) in CANDLE. 134 dysregulated proteins (85 upregulated, 49 downregulated) overlapped in NOMID and CANDLE (Figure 1). Pathway analysis identified neutrophil and monocyte chemotaxis signature in both NOMID and CANDLE. NOMID patients had neutrophilia and active neutrophils. CANDLE patients exhibited active neutrophils in whole blood RNA. Endothelial cell activation was the most prominent non-hematopoietic signature and suggest distinct endothelial cell dysregulation in NOMID and CANDLE. In NOMID, the signature included neutrophil transmigration (SELE) endothelial cell motility in response to angiogenesis (HGF, VEGF), while in CANDLE the endothelial signatures included extracellular matrix protein deposition (COL8A) suggesting increased vascular stiffness. CANDLE patients had higher expression of Renin, 4 out of 7 had hypertension, NOMID patients did not have hypertension. Treatment with anakinra and baricitinib normalized 143 and 142 of dysregulated proteins in NOMID and CANDLE respectively.Conclusion:Differentially expressed proteins in NOMID and CANDLE are consistent with innate immune cell activation. Distinct endothelial cell signatures in NOMID and CANDLE may provide mechanistic insight into differences in vascular phenotypes. Treatment with anakinra and Baricitinib in NOMID and CANDLE leaves 30% and 13% of the dysregulated proteins unchanged.Acknowledgments:This work was supported by Intramural Research atNational Institute of Allergy Immunology and Infectious Diseases of National Institutes of Health, Bethesda, Maryland, the Center of Human Immunology and was approved by the IRB.Disclosure of Interests:None declared


2021 ◽  
Vol 15 (1) ◽  
Author(s):  
Weitong Cui ◽  
Huaru Xue ◽  
Lei Wei ◽  
Jinghua Jin ◽  
Xuewen Tian ◽  
...  

Abstract Background RNA sequencing (RNA-Seq) has been widely applied in oncology for monitoring transcriptome changes. However, the emerging problem that high variation of gene expression levels caused by tumor heterogeneity may affect the reproducibility of differential expression (DE) results has rarely been studied. Here, we investigated the reproducibility of DE results for any given number of biological replicates between 3 and 24 and explored why a great many differentially expressed genes (DEGs) were not reproducible. Results Our findings demonstrate that poor reproducibility of DE results exists not only for small sample sizes, but also for relatively large sample sizes. Quite a few of the DEGs detected are specific to the samples in use, rather than genuinely differentially expressed under different conditions. Poor reproducibility of DE results is mainly caused by high variation of gene expression levels for the same gene in different samples. Even though biological variation may account for much of the high variation of gene expression levels, the effect of outlier count data also needs to be treated seriously, as outlier data severely interfere with DE analysis. Conclusions High heterogeneity exists not only in tumor tissue samples of each cancer type studied, but also in normal samples. High heterogeneity leads to poor reproducibility of DEGs, undermining generalization of differential expression results. Therefore, it is necessary to use large sample sizes (at least 10 if possible) in RNA-Seq experimental designs to reduce the impact of biological variability and DE results should be interpreted cautiously unless soundly validated.


Pathogens ◽  
2021 ◽  
Vol 10 (4) ◽  
pp. 405
Author(s):  
Anna Matysiak ◽  
Michal Kabza ◽  
Justyna A. Karolak ◽  
Marcelina M. Jaworska ◽  
Malgorzata Rydzanicz ◽  
...  

The ocular microbiome composition has only been partially characterized. Here, we used RNA-sequencing (RNA-Seq) data to assess microbial diversity in human corneal tissue. Additionally, conjunctival swab samples were examined to characterize ocular surface microbiota. Short RNA-Seq reads, obtained from a previous transcriptome study of 50 corneal tissues, were mapped to the human reference genome GRCh38 to remove sequences of human origin. The unmapped reads were then used for taxonomic classification by comparing them with known bacterial, archaeal, and viral sequences from public databases. The components of microbial communities were identified and characterized using both conventional microbiology and polymerase chain reaction (PCR) techniques in 36 conjunctival swabs. The majority of ocular samples examined by conventional and molecular techniques showed very similar microbial taxonomic profiles, with most of the microorganisms being classified into Proteobacteria, Firmicutes, and Actinobacteria phyla. Only 50% of conjunctival samples exhibited bacterial growth. The PCR detection provided a broader overview of positive results for conjunctival materials. The RNA-Seq assessment revealed significant variability of the corneal microbial communities, including fastidious bacteria and viruses. The use of the combined techniques allowed for a comprehensive characterization of the eye microbiome’s elements, especially in aspects of microbiota diversity.


Animals ◽  
2021 ◽  
Vol 11 (5) ◽  
pp. 1268
Author(s):  
Shengchao Zhang ◽  
Sibtain Ahmad ◽  
Yuxia Zhang ◽  
Guohua Hua ◽  
Jianming Yi

Enhanced plane of nutrition at pre-weaning stage can promote the development of mammary gland especially heifer calves. Although several genes are involved in this process, long intergenic non-coding RNAs (lincRNAs) are regarded as key regulators in the regulated network and are still largely unknown. We identified and characterized 534 putative lincRNAs based on the published RNA-seq data, including heifer calves in two groups: fed enhanced milk replacer (EH, 1.13 kg/day, including 28% crude protein, 25% fat) group and fed restricted milk replacer (R, 0.45 kg/day, including 20% crude protein, 20% fat) group. Sub-samples from the mammary parenchyma (PAR) and mammary fat pad (MFP) were harvested from heifer calves. According to the information of these lincRNAs’ quantitative trait loci (QTLs), the neighboring and co-expression genes were used to predict their function. By comparing EH vs R, 79 lincRNAs (61 upregulated, 18 downregulated) and 86 lincRNAs (54 upregulated, 32 downregulated) were differentially expressed in MFP and PAR, respectively. In MFP, some differentially expressed lincRNAs (DELs) are involved in lipid metabolism pathways, while, in PAR, among of DELs are involved in cell proliferation pathways. Taken together, this study explored the potential regulatory mechanism of lincRNAs in the mammary gland development of calves under different planes of nutrition.


Author(s):  
Mayukh Banerjee ◽  
Ana Ferragut Cardoso ◽  
Laila Al-Eryani ◽  
Jianmin Pan ◽  
Theodore S. Kalbfleisch ◽  
...  

AbstractChronic arsenic exposure causes skin cancer, although the underlying molecular mechanisms are not well defined. Altered microRNA and mRNA expression likely play a pivotal role in carcinogenesis. Changes in genome-wide differential expression of miRNA and mRNA at 3 strategic time points upon chronic sodium arsenite (As3+) exposure were investigated in a well-validated HaCaT cell line model of arsenic-induced cutaneous squamous cell carcinoma (cSCC). Quadruplicate independent HaCaT cell cultures were exposed to 0 or 100 nM As3+ for up to 28-weeks (wk). Cell growth was monitored throughout the course of exposure and epithelial-mesenchymal transition (EMT) was examined employing immunoblot. Differentially expressed miRNA and mRNA profiles were generated at 7, 19, and 28-wk by RNA-seq, followed by identification of differentially expressed mRNA targets of differentially expressed miRNAs through expression pairing at each time point. Pathway analyses were performed for total differentially expressed mRNAs and for the miRNA targeted mRNAs at each time point. RNA-seq predictions were validated by immunoblot of selected target proteins. While the As3+-exposed cells grew slower initially, growth was equal to that of unexposed cells by 19-wk (transformation initiation), and exposed cells subsequently grew faster than passage-matched unexposed cells. As3+-exposed cells had undergone EMT at 28-wk. Pathway analyses demonstrate dysregulation of carcinogenesis-related pathways and networks in a complex coordinated manner at each time point. Immunoblot data largely corroborate RNA-seq predictions in the endoplasmic reticulum stress (ER stress) pathway. This study provides a detailed molecular picture of changes occurring during the arsenic-induced transformation of human keratinocytes.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Zefang Sun ◽  
Jia Tan ◽  
Minqiong Zhao ◽  
Qiyao Peng ◽  
Mingqing Zhou ◽  
...  

AbstracttRNAs and tRNA-derived RNA fragments (tRFs) play various roles in many cellular processes outside of protein synthesis. However, comprehensive investigations of tRNA/tRF regulation are rare. In this study, we used new algorithms to extensively analyze the publicly available data from 1332 ChIP-Seq and 42 small-RNA-Seq experiments in human cell lines and tissues to investigate the transcriptional and posttranscriptional regulatory mechanisms of tRNAs. We found that histone acetylation, cAMP, and pluripotency pathways play important roles in the regulation of the tRNA gene transcription in a cell-specific manner. Analysis of RNA-Seq data identified 950 high-confidence tRFs, and the results suggested that tRNA pools are dramatically distinct across the samples in terms of expression profiles and tRF composition. The mismatch analysis identified new potential modification sites and specific modification patterns in tRNA families. The results also show that RNA library preparation technologies have a considerable impact on tRNA profiling and need to be optimized in the future.


Sign in / Sign up

Export Citation Format

Share Document