Metabolomics,Unknown,Transcriptomics,Genomics,Proteomics

Dataset Information

Transgenic Mouse Models Support a Protective Role of Type I IFN Response in SARS CoV-2 infection-related Immunopathology

ABSTRACT: RNA-seq analysis was performed to understand the role of type I IFN response during SARS CoV-2 infection using transgenic mice. Each sample was collected from an individual C57BL/6J mouse. The total RNA was extracted from uninfected and SARS-CoV-2 infected mice lung tissue using RNeasy mini kit (QIAGEN #74104). The quantity of RNA was determined using Qubit RNA assay kit with Qubit 4.0 and the quality of RNA was tested using agarose gel electrophoresis and High Sensitivity Tape station Kit (Agilent 2200, #5067-5576, #5067-5577 and #5067-5578). After assessing the quality of RNA, ~900 ng of total RNA was taken for library preparation using NEBNext®Ultra™ II Directional RNA Library kit for Illumina (# E7760L) and NEBNext Poly (A) mRNA Magnetic Isolation Module (# E7490L) as per manufacturer's protocol. The prepared library was quantified using Qubit dsDNA assay kit (Invitrogen, Q32851) followed by quality check (QC) and fragment size distribution using a High Sensitivity Tape station Kit (Agilent 2200, #5067-5584 and #5067-5585). The library was sequenced using the HiSeq 4000 Illumina platform. The paired-end (PE) reads quality checks for each sample were carried out using FastQC v.0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). The adapter sequence was trimmed using the BBDuk version 37.58 version 37.58 and the alignment was performed using STAR v.2.5.3a with default parameters with human hg38 genome build, gencode v21 gtf 9GRCh38) from the gencode. The duplicates were discarded using Picard-2.9.4 (https://broadinstitute.github.io/picard/) from the aligned bam files and read counts were generated using featureCount v.1.5.3 from subread-1.5.3 package (https://bioinf.wehi.edu.au/) with Q = 10 for mapping quality. The count files were used as input for downstream differential gene expression analysis with DESeq2 version 1.14.1 9. The genes with read counts of ≤ 10 in any comparison were discarded followed by count transformation and statistical analysis using DESeq “R”. The “P” value were adjusted using the Benjamini and Hochberg multiple testing correction and the differentially expressed genes were identified (fold change of ≥1.5, P-value < 0.05). A unified non-redundant gene list was made for different comparisons and subjected to gene ontology (GO) analysis using the reactome database (https://reactome.org/). The top pathways (p < 0.05) were used for generating heat maps using Complexheatmap (Version 2.0.0) through unsupervised hierarchical clustering. The expression clusters were annotated based on enriched GO terms. Normalized gene expression was used to generate the boxplots with a median depicting the trends in the expression across the different conditions using ggplot2 [version 3.3.5]. The pathways analysis was performed using Metascape database (https://metascape.org/gp/index.html#/main/step1). The top pathways (p < 0.05) were taken for constructing bubble plots using ggplot2 [version 3.3.5].

INSTRUMENT(S): Illumina HiSeq 4000

ORGANISM(S): Mus musculus

SUBMITTER: SANTOSH CHAUHAN

PROVIDER: E-MTAB-12681 | biostudies-arrayexpress |

REPOSITORIES: biostudies-arrayexpress

ACCESS DATA

Similar Datasets

Project description:Purpose: Zinc Finger MIZ-Type Containing 1 (Zmiz1) is a member of the PIAS family of protein and function as a transcriptional coactivator of Notch, Androgen Receptor (AR), p53, Estrogen Receptor (ER), and Smad3/4 . Despite Zmiz1 critical role in angiogenesis, its role in lymphatic vasculature is unknown. Here, we use HDLECs cell line to profile transcriptional changes upon Zmiz1 knockdown using siRNA. Methods: Total RNA was extracted from control and Zmiz1 siRNA transfected HDLECs using GeneJET RNA Purification Kit (Thermo Fisher Scientific, K0732) following the manufacturer instructions. RNA concentration and RNA integrity (RIN) number were determined using Qubit RNA High Sensitivity Assay Kit (Thermo Fisher Scientific, Q32852) and Bioanalyzer RNA 6000 Nano assay kit (Agilent, 5067-1511) respectively. RNA library was prepared using TruSeq RNA Library Prep Kit v2 (Illumina, RS-122-2001) according to the manufacturer instructions. mRNA library was quantified using Qubit dsDNA High Sensitivity Assay Kit (Thermo Fisher Scientific, Q32851) and verified using the Bioanalyzer DNA1000 assay kit (Agilent, 5067-1505). Verified samples were sequenced using the NextSeq1000/2000 P2 Reagents (200 Cycles) v3 (Illumina, 20046812) on a Nextseq1000/2000. Sequenced reads were aligned to the human (hg19) reference genome with RNA-Seq alignment tool (STAR aligner). The aligned reads were used to quantify mRNA expression and determine differentially expressed genes using the RNA-Seq Differential Expression tool (version 1.0.1). Both alignment and differential expression analysis were performed using the tools in the illumina BaseSpace Sequence Hub. Results: We found 2,228 differentially expressed genes of which 1,115 genes were upregulated while 1,113 genes were downregulated. Downregulated genes were enriched in biological processes such as lymph vessel development, cell migration and heart valve development. Conclusions: We identify Zmiz1 regulates Prox1 in lymphatic endothelial cells

Project description:Purpose: Next-generation sequencing (NGS) has revolutionized systems-based analysis of cellular pathways. The goals of this study are to compare adhesion gene clusters and clock genes / biomarkers at different time points throughout the day in sorted murine lympahtic endothelial cells. Methods: RNA of sorted dermal LECs harvested in TriZol-LS was purified using Direct-zol RNA mini Prep Kit (ZymoResearch) following manufacturer’s instructions. Isolated RNA was quantified using the Nanodrop (Thermofisher) and analyzed on a Bioanalyzer (Agilent) using the RNA 6000 Nano or Pico Kit (Agilent). 75 ng of eluted total RNA was digested with DNase (Thermofisher) to remove DNA contaminations. An additional purification step with RNAClean XP Beads (Agencourt) was performed. The purified, bead-bound RNA was directly used as input in the SMART-Seq v4 Ultra Low Input RNA Kit (TaKaRA Bio). Full-length cDNA was generated following manufacturer’s instructions. Full-length cDNA was quantified using the Qubit dsDNA HS Assay Kit (Invitrogen) with the Qubit fluorometer. Finally, sequencing libraries were generated using 500 pg full-length cDNA following the NexteraXT protocol (Illumina). Libraries were quantified on a Bioanalyzer (Agilent) using the DNA 1000 Kit (Agilent) and sequenced on a HiSeq1500 system (Illumina) with a readlength of 100 nt, single-end mode. Results: obtained transcriptome profiles were processed on a Galaxy web interface hosted by LAFUGA (Gene Center, Munich). After demultiplexing and trimming data was mapped against the mouse genome (mm10) using RNA-STAR mapper (Galaxy Version 2.5.2b-0). Abundant reads were counted using HTSeq-count (Galaxy Version 1.0.0). Afterwards, gene expression analysis to detect differentially expressed genes was performed using DESeq2 (Galaxy Version 2.11.40.6), setting the false discovery rate (FDR) <0.05. Comparing the expression profiles and the enrichment of gene ontology cluster adhesion revealed a highly time-of-day dependent expression profile of genes. Conclusions: Our study represents the first detailed analysis of time-of-day dependent lymphatic endothelial cell transcriptomes, with respect to adhesion and clock genes. This will provide valid input into chronotherapy, timed delivery of substances and new therapeutic avenues.

Project description:Purpose: ATP6AP2 is a single transmembrane protein that has been implicated in regulating various biological process including vascular angiogenesis. However, the exact role of ATP6AP2 during blood vessel development remains largely undefined. Here, we use an inducible endothelial cell (EC)-specific mouse to investigate the role of murine Atp6ap2 gene during both physiological and pathological angiogenesis of the postnatal retina, a well-established model system of angiogenesis. In this study we aimed to Identify what genes are differentially expressed upon loss of PRR(Atp6ap2) in ECs. Methods: ECs were isolated from lungs of both WT and PRR mutant mice at postnatal day 7 (P7). Total RNA was extracted from isolated lung ECs and quantified using Qubit RNA High Sensitivity Assay Kit. RNA integrity was determined using the Bioanalyzer RNA 6000 Nano assay kit. RNA library construction was performed with the TruSeq RNA Library Prep Kit v2 from Illumina. The resulting mRNA library was quantified using Qubit dsDNA High Sensitivity Assay Kit and verified using the Bioanalyzer DNA1000 assay kit. Verified samples were sequenced using the NextSeq 500/550 High Output Kit v2.5 (150 Cycles) on a Nextseq 550 system. Sequenced reads were aligned to the mouse (mm10) reference genome with RNA-seq alignment tool (version 2.0.1). The aligned reads were used to quantify mRNA expression and determine differentially expressed genes using the RNA-seq Differential Expression tool (version 1.0.1). Both alignment and differential expression analysis were performed using the tools in the BaseSpace Sequence Hub. Results: We found 1876 differentially expressed genes of which 823 genes were upregulated while 1053 genes were downregulated. Downregulated genes were enriched in biological processes such as angiogenesis, cell migration and extracellular matrix remodelling, while upregulated genes were enriched in DNA and cell cycle regulation. Conclusions: We identify PRR (ATP6AP2) to play an important role in regulation of the angiogenesis process.

Dataset Information

Transgenic Mouse Models Support a Protective Role of Type I IFN Response in SARS CoV-2 infection-related Immunopathology

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure