scholarly journals Gene co-expression network analysis in zebrafish reveals chemical class specific modules

BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Prarthana Shankar ◽  
Ryan S. McClure ◽  
Katrina M. Waters ◽  
Robyn L. Tanguay

Abstract Background Zebrafish is a popular animal model used for high-throughput screening of chemical hazards, however, investigations of transcriptomic mechanisms of toxicity are still needed. Here, our goal was to identify genes and biological pathways that Aryl Hydrocarbon Receptor 2 (AHR2) Activators and flame retardant chemicals (FRCs) alter in developing zebrafish. Taking advantage of a compendium of phenotypically-anchored RNA sequencing data collected from 48-h post fertilization (hpf) zebrafish, we inferred a co-expression network that grouped genes based on their transcriptional response. Results Genes responding to the FRCs and AHR2 Activators localized to distinct regions of the network, with FRCs inducing a broader response related to neurobehavior. AHR2 Activators centered in one region related to chemical stress responses. We also discovered several highly co-expressed genes in this module, including cyp1a, and we subsequently show that these genes are definitively within the AHR2 signaling pathway. Systematic removal of the two chemical types from the data, and analysis of network changes identified neurogenesis associated with FRCs, and regulation of vascular development associated with both chemical classes. We also identified highly connected genes responding specifically to each class that are potential biomarkers of exposure. Conclusions Overall, we created the first zebrafish chemical-specific gene co-expression network illuminating how chemicals alter the transcriptome relative to each other. In addition to our conclusions regarding FRCs and AHR2 Activators, our network can be leveraged by other studies investigating chemical mechanisms of toxicity.

2021 ◽  
Author(s):  
Mathew G. Lewsey ◽  
Changyu Yi ◽  
Oliver Berkowitz ◽  
Felipe Ayora ◽  
Maurice Bernado ◽  
...  

SummaryThe Single Cell App is a cloud-based application that allows visualisation and comparison of scRNA-seq data and is scalable according to use. Users upload their own or publicly available scRNA-seq datasets after pre-processing to be visualised using a web browser. The data can be viewed in two colour modes, Cluster - representing cell identity, and Values – level of expression, and data queried using keyword or gene identification number(s). Using the app to compare four different studies we determined that some genes frequently used as cell-type markers are in fact study specific. Phosphate transporter and hormone response genes were exemplary investigated with the app. This showed that the apparent cell specific expression of PHO1;H3 differed between GFP-tagging and scRNA-seq studies. Some phosphate transporter genes were induced by protoplasting, they retained cell specificity, indicating that cell specific stress responses (i.e. protoplasting). Examination of the cell specificity of hormone response genes revealed that 132 hormone responsive genes display restricted expression and that the jasmonate response gene TIFY8 is expressed in endodermal cells which differs from previous reports. It also appears that JAZ repressors have cell-type specific functions. These differences, identified using the Single Cell App, highlight the need for resources to enable researchers to find common and different patterns of cell specific gene expression. Thus, the Single Cell App enables researchers to form new hypothesis, perform comparative studies, allows for easy re-use of data for this emerging technology to provide novel avenues to crop improvement.


2021 ◽  
Vol 9 (7) ◽  
pp. e002383
Author(s):  
Jin-Li Wei ◽  
Si-Yu Wu ◽  
Yun-Song Yang ◽  
Yi Xiao ◽  
Xi Jin ◽  
...  

PurposeRegulatory T cells (Tregs) heavily infiltrate triple-negative breast cancer (TNBC), and their accumulation is affected by the metabolic reprogramming in cancer cells. In the present study, we sought to identify cancer cell-intrinsic metabolic modulators correlating with Tregs infiltration in TNBC.Experimental designUsing the RNA-sequencing data from our institute (n=360) and the Molecular Taxonomy of Breast Cancer International Consortium TNBC cohort (n=320), we calculated the abundance of Tregs in each sample and evaluated the correlation between gene expression levels and Tregs infiltration. Then, in vivo and in vitro experiments were performed to verify the correlation and explore the underlying mechanism.ResultsWe revealed that GTP cyclohydrolase 1 (GCH1) expression was positively correlated with Tregs infiltration and high GCH1 expression was associated with reduced overall survival in TNBC. In vivo and in vitro experiments showed that GCH1 increased Tregs infiltration, decreased apoptosis, and elevated the programmed cell death-1 (PD-1)-positive fraction. Metabolomics analysis indicated that GCH1 overexpression reprogrammed tryptophan metabolism, resulting in L-5-hydroxytryptophan (5-HTP) accumulation in the cytoplasm accompanied by kynurenine accumulation and tryptophan reduction in the supernatant. Subsequently, aryl hydrocarbon receptor, activated by 5-HTP, bound to the promoter of indoleamine 2,3-dioxygenase 1 (IDO1) and thus enhanced the transcription of IDO1. Furthermore, the inhibition of GCH1 by 2,4-diamino-6-hydroxypyrimidine (DAHP) decreased IDO1 expression, attenuated tumor growth, and enhanced the tumor response to PD-1 blockade immunotherapy.ConclusionsTumor-cell-intrinsic GCH1 induced immunosuppression through metabolic reprogramming and IDO1 upregulation in TNBC. Inhibition of GCH1 by DAHP serves as a potential immunometabolic strategy in TNBC.


Plants ◽  
2021 ◽  
Vol 10 (2) ◽  
pp. 267
Author(s):  
Axel J. Giudicatti ◽  
Ariel H. Tomassi ◽  
Pablo A. Manavella ◽  
Agustin L. Arce

MicroRNAs are small regulatory RNAs involved in several processes in plants ranging from development and stress responses to defense against pathogens. In order to accomplish their molecular functions, miRNAs are methylated and loaded into one ARGONAUTE (AGO) protein, commonly known as AGO1, to stabilize and protect the molecule and to assemble a functional RNA-induced silencing complex (RISC). A specific machinery controls miRNA turnover to ensure the silencing release of targeted-genes in given circumstances. The trimming and tailing of miRNAs are fundamental modifications related to their turnover and, hence, to their action. In order to gain a better understanding of these modifications, we analyzed Arabidopsis thaliana small RNA sequencing data from a diversity of mutants, related to miRNA biogenesis, action, and turnover, and from different cellular fractions and immunoprecipitations. Besides confirming the effects of known players in these pathways, we found increased trimming and tailing in miRNA biogenesis mutants. More importantly, our analysis allowed us to reveal the importance of ARGONAUTE 1 (AGO1) loading, slicing activity, and cellular localization in trimming and tailing of miRNAs.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Kai Zhao ◽  
Song Chen ◽  
Wenjing Yao ◽  
Zihan Cheng ◽  
Boru Zhou ◽  
...  

Abstract Background The bZIP gene family, which is widely present in plants, participates in varied biological processes including growth and development and stress responses. How do the genes regulate such biological processes? Systems biology is powerful for mechanistic understanding of gene functions. However, such studies have not yet been reported in poplar. Results In this study, we identified 86 poplar bZIP transcription factors and described their conserved domains. According to the results of phylogenetic tree, we divided these members into 12 groups with specific gene structures and motif compositions. The corresponding genes that harbor a large number of segmental duplication events are unevenly distributed on the 17 poplar chromosomes. In addition, we further examined collinearity between these genes and the related genes from six other species. Evidence from transcriptomic data indicated that the bZIP genes in poplar displayed different expression patterns in roots, stems, and leaves. Furthermore, we identified 45 bZIP genes that respond to salt stress in the three tissues. We performed co-expression analysis on the representative genes, followed by gene set enrichment analysis. The results demonstrated that tissue differentially expressed genes, especially the co-expressing genes, are mainly involved in secondary metabolic and secondary metabolite biosynthetic processes. However, salt stress responsive genes and their co-expressing genes mainly participate in the regulation of metal ion transport, and methionine biosynthetic. Conclusions Using comparative genomics and systems biology approaches, we, for the first time, systematically explore the structures and functions of the bZIP gene family in poplar. It appears that the bZIP gene family plays significant roles in regulation of poplar development and growth and salt stress responses through differential gene networks or biological processes. These findings provide the foundation for genetic breeding by engineering target regulators and corresponding gene networks into poplar lines.


Antibiotics ◽  
2020 ◽  
Vol 9 (11) ◽  
pp. 808
Author(s):  
Maurice Steenhuis ◽  
Corinne M. ten Hagen-Jongman ◽  
Peter van Ulsen ◽  
Joen Luirink

The structural integrity of the Gram-negative cell envelope is guarded by several stress responses, such as the σE, Cpx and Rcs systems. Here, we report on assays that monitor these responses in E. coli upon addition of antibacterial compounds. Interestingly, compromised peptidoglycan synthesis, outer membrane biogenesis and LPS integrity predominantly activated the Rcs response, which we developed into a robust HTS (high-throughput screening) assay that is suited for phenotypic compound screening. Furthermore, by interrogating all three cell envelope stress reporters, and a reporter for the cytosolic heat-shock response as control, we found that inhibitors of specific envelope targets induce stress reporter profiles that are distinct in quality, amplitude and kinetics. Finally, we show that by using a host strain with a more permeable outer membrane, large-scaffold antibiotics can also be identified by the reporter assays. Together, the data suggest that stress profiling is a useful first filter for HTS aimed at inhibitors of cell envelope processes.


2021 ◽  
Vol 36 (Supplement_1) ◽  
Author(s):  
Floranne Boulogne ◽  
Laura Claus ◽  
Henry Wiersma ◽  
Roy Oelen ◽  
Floor Schukking ◽  
...  

Abstract Background and Aims Genetic testing in patients with suspected hereditary kidney disease does not always reveal the genetic cause for the patient's disorder. Potentially pathogenic variants can reside in genes that are not known to be involved in kidney disease, which makes it difficult to prioritize and interpret the relevance of these variants. As such, there is a clear need for methods that predict the phenotypic consequences of gene expression in a way that is as unbiased as possible. To help identify candidate genes we have developed KidneyNetwork, in which tissue-specific expression is utilized to predict kidney-specific gene functions. Method We combined gene co-expression in 878 publicly available kidney RNA-sequencing samples with the co-expression of a multi-tissue RNA-sequencing dataset of 31,499 samples to build KidneyNetwork. The expression patterns were used to predict which genes have a kidney-related function, and which (disease) phenotypes might be caused when these genes are mutated. By integrating the information from the HPO database, in which known phenotypic consequences of disease genes are annotated, with the gene co-expression network we obtained prediction scores for each gene per HPO term. As proof of principle, we applied KidneyNetwork to prioritize variants in exome-sequencing data from 13 kidney disease patients without a genetic diagnosis. Results We assessed the prediction performance of KidneyNetwork by comparing it to GeneNetwork, a multi-tissue co-expression network we previously developed. In KidneyNetwork, we observe a significantly improved prediction accuracy of kidney-related HPO-terms, as well as an increase in the total number of significantly predicted kidney-related HPO-terms (figure 1). To examine its clinical utility, we applied KidneyNetwork to 13 patients with a suspected hereditary kidney disease without a genetic diagnosis. Based on the HPO terms “Renal cyst” and “Hepatic cysts”, combined with a list of potentially damaging variants in one of the undiagnosed patients with mild ADPKD/PCLD, we identified ALG6 as a new candidate gene. ALG6 bears a high resemblance to other genes implicated in this phenotype in recent years. Through the 100,000 Genomes Project and collaborators we identified three additional patients with kidney and/or liver cysts carrying a suspected deleterious variant in ALG6. Conclusion We present KidneyNetwork, a kidney specific co-expression network that accurately predicts what genes have kidney-specific functions and may result in kidney disease. Gene-phenotype associations of genes unknown for kidney-related phenotypes can be predicted by KidneyNetwork. We show the added value of KidneyNetwork by applying it to exome sequencing data of kidney disease patients without a molecular diagnosis and consequently we propose ALG6 as a promising candidate gene. KidneyNetwork can be applied to clinically unsolved kidney disease cases, but it can also be used by researchers to gain insight into individual genes to better understand kidney physiology and pathophysiology. Acknowledgments This research was made possible through access to the data and findings generated by the 100,000 Genomes Project; http://www.genomicsengland.co.uk.


2020 ◽  
Author(s):  
Albert Maibam ◽  
Sunil Nigombam ◽  
Harinder Vishwakarma ◽  
Showkat Ahmad Lone ◽  
Kishor Gaikwad ◽  
...  

Abstract Background Pennisetum glaucum (L.) R. Br. is mainly grown in arid and semi-arid regions. Being naturally tolerant to various adverse condtitions, it is a good biological resource for deciphering the molecular basis of abiotic stresses such as heat stress in plants but limited studies have been carried out till date to this effect. Here, we performed RNA-sequencing from the leaf of two contrasting genotypes of pearl millet (841-B and PPMI-69) subjected to heat stress (42 °C for 6 h). Results Over 274 million high quality reads with an average length of 150 nt were generated. Assembly was carried out using trinity, obtaining 47,310 unigenes having an average length of 1254 nucleotides, N50 length of 1853 nucleotides and GC content of 53.11%. Blastx resulted in annotation of 35,628 unigenes and functional classification showed 15,950 unigenes designated to 51 Gene Ontology terms, 13,786 unigenes allocated to 23 Clusters of Orthologous Groups and 4,255 unigenes distributed into 132 functional KEGG pathways. 12,976 simple sequence repeats were identified from 10,294 unigenes for the development of functional markers. A total of 3,05,759 SNPs were observed in the transcriptome data. Out of 2,301 differentially expressed genes, 10 potential candidates genes were selected based on log2 fold change and adjusted p-value parameters for their differential gene expression by qRT-PCR. Conclusions The dynamic expression changes in two genotypes of P. glaucum reflect transcriptome regulation of signaling pathways in heat stress response. In order to develop genetic markers, 12,976 simple sequence repeats (SSRs) were identified. The sequencing data generated in this study shall serve as an important resource for further research in the area of crop biotechnology.


Author(s):  
Roman Farooq Alvi ◽  
Bilal Aslam ◽  
Muhammad Hidayat Rasool ◽  
Saima Muzammil ◽  
Abu Baker Siddique ◽  
...  

Background. The term “persisters” refers to a small bacterial population that persists during treatment with high antibiotic concentration or dose in the absence of genetic resistance. The present study was designed to investigate the transcriptional response in indigenous Klebsiella pneumoniae under the ciprofloxacin stress. Methods. Isolation and identification of K. pneumoniae were carried out through standard microbiological protocols. The characterization of quinolone resistance was performed by estimating the quinolone susceptibility testing, MIC estimation, and detecting the QRDR and PMQR. Transcriptional response of the isolates to ciprofloxacin was determined using qPCR. Results. Among 34 isolates, 23 (67%) were resistant to ciprofloxacin. Both QRDR (gyrA and gyrB) and PMQR (qnrA, qnrB, and qnrS) were detected in the isolates, and all were found resistant to ciprofloxacin. The mRNA levels of both mutS and euTu under the influence of ciprofloxacin were significantly increased. On ciprofloxacin exposure, the mRNA levels of the DNA damage response element (mutS) were raised in a time-dependent fashion. K. pneumoniae showed high-level resistance to ciprofloxacin in the presence of mutations in QRDR and PMQR genes. Conclusion. The transcriptional response revealed the upregulation of DNA repair and protein folding elements (mutS and euTu) in ciprofloxacin stress and delayed cell division. The ciprofloxacin was found to trigger various stress responses in a time- and concentration-dependent manner.


2021 ◽  
Author(s):  
Wen Feng ◽  
Lei Zhou ◽  
Pengju Zhao ◽  
Heng Du ◽  
Chenguang Diao ◽  
...  

As warthog (Phacochoerus africanus) has innate immunity against African swine fever (ASF), it is critical to understanding the evolutionary novelty of warthog to explain its specific ASF resistance. Here, we present two completed new genomes of one warthog and one Kenyan domestic pig, as the fundamental genomic references to decode the genetic mechanism on ASF tolerance. Our results indicated, multiple genomic variations, including gene losses, independent contraction and expansion of specific gene families, likely moulded warthog's genome to adapt the environment. Importantly, the analysis of presence and absence of genomic sequences revealed that, the warthog genome had a DNA sequence absence of the lactate dehydrogenase B (LDHB) gene on chromosome 2 compared to the reference genome. The overexpression and siRNA of LDHB indicated that its inhibition on the replication of ASFV. The Combining with large scale sequencing data of 123 pigs from all over world, contraction and expansion of TRIM genes families revealed that TRIM family genes in the warthog genome were potentially responsible for its tolerance to ASF. Our results will help further improve the understanding of genetic resistance ASF in pigs.


Sign in / Sign up

Export Citation Format

Share Document