Next Generation Sequencing Facilitates Quantitative Analysis of Wild Type and Sirt3K223R Bone-Marrow derived Macrophages Transcriptomes
Ontology highlight
ABSTRACT: Purpose: Next-generation sequencing (NGS) has revolutionized systems-based analysis of cellular pathways. The goals of this study are to compare Transcriptomes of Sirt3 WT and KR macrophages with high-throughput data analysis Methods: Macrophages' mRNA profiles of 8 weeks-old wild-type (WT) and Sirt3 K223R mice were generated by deep sequencing, in quadruplicate, using Total RNA from each sample was quantified using a NanoDrop ND-1000 instrument. 1μg total RNA was used to prepare the sequencing library in the following steps: 1. Total RNA is enriched by oligo (dT) magnetic beads (rRNA removed); 2. RNA-seq library preparation using KAPA Stranded RNA-Seq Library Prep Kit (Illumina), which incorporates dUTP into the second cDNA strand and renders the RNA-seq library strand-specific. The completed libraries were qualified with Agilent 2100 Bioanalyzer and quantified by absolute quantification qPCR method. To sequence the libraries on the Illumina HiSeq 4000 instrument, the barcoded libraries were mixed, denatured to single stranded DNA in NaOH, captured on Illumina flow cell, amplified in situ, and subsequently sequenced for 150 cycles for both ends on Illumina HiSeq instrument. Results: Image analysis and base calling were performed using Solexa pipeline v1.8 (Off-Line Base Caller software, v1.8). Sequence quality was examined using the FastQC software. The trimmed reads (trimmed 5’, 3’-adaptor bases using cutadapt) were aligned to reference genome using Hisat2 software. The transcript abundances for each sample was estimated with StringTie , and the FPKM value for gene and transcript level were calculated with R package Ballgown. The differentially expressed genes and transcripts were filtered using R package Ballgown. The novel genes and transcripts were predicted from assembled results by comparing to the reference annotation using StringTie and Ballgown, then use CPAT to assess the coding potential of those sequences. Then use rMATS to detecting alternative splicing events and plots. Principle Component Analysis (PCA) and correlation analysis were based on gene expression level, Hierarchical Clustering, Gene Ontology, Pathway analysis, Gene Ontology, Pathway analysis, scatter plots and volcano plots were performed with the differentially expressed genes in R, Python or shell environment for statistical computing and graphics. Conclusions: Our study represents the first detailed analysis of retinal transcriptomes, with 4 biologic replicates, generated by RNA-seq technology. The optimized data analysis workflows reported here should provide a framework for comparative investigations of expression profiles. Our results offer a comprehensive and more accurate quantitative and qualitative evaluation of mRNA content within a cell or tissue. We conclude that RNA-seq based transcriptome characterization would expedite genetic network analyses and permit the dissection of complex biologic functions.
Project description:Purpose: The goals of this study are to compare type 2 diabetes lacrimal gland transcriptome profiling (RNA-seq) with the counterparts. Methods: Lacrimal gland mRNA profiles of 24 weeks T2DM mice (C57 mice with high fat diet + STZ intraperitoneal injection) and counterparts were enriched by NEBNext® Poly(A) mRNA Magnetic Isolation Module (New England Biolabs), and the produced RNA was used for construction library, via KAPA Stranded RNA-Seq Library Prep Kit (Illumina). The prepared RNA-seq libraries were qualified using Agilent 2100 Bioanalyzer and quantified by qPCR absolute quantification method. The sequencing was then performed using Illumina NovaSeq 6000. After quality control, the fragments were 5’, 3’-adaptor trimmed and filtered <=16 bp reads with cutadapt software. The trimmed reads were aligned to a reference genome with Hisat2 software. The transcript abundance for each sample was estimated with StringTie (v1.2.3), and the FPKM values for gene level were calculated with R package Ballgown (v2.6.0). The differentially expressed genes (DEGs) analysis was also performed with Ballgown. Fold change (cutoff 1.5), p-value (cutoff 0.05) and FPKM (≥0.5 mean in one group) were used for filtering differentially expressed genes and transcripts. Conclusions: Our study represents the first detailed analysis of T2DM lacrimal gland transcriptomes, with biologic replicates, generated by RNA-seq technology. The optimized data analysis workflows reported here should provide a framework for comparative investigations of expression profiles. Our results showed that the vascular related genes were altered in the T2DM.
Project description:Purpose: The goal of this study is to analyze microRNA affected by wuho and mei-p26 mutations in Drosophila ovary Ovarian microRNA prolife of 7 day old control flies and flies bearing wuho and mei-p26 mutations, in duplicate.Sample libraries were prepared using the QIAseq miRNA Library Kit (QIAGEN) according to the manufacturer's protocols. Adaptors are ligated sequentially to the 3' and 5' ends of miRNAs. Subsequently, universal cDNA synthesis with UMI assignment, cDNA cleanup, library amplification and library cleanup are performed. Libraries were sequenced on an Illumina instrument (75-cycle single-end read, 75SE). Sequencing data was processed using the Illumina software BCL2FASTQ v2.20.
Project description:RNA-sequencing analysis from whole heart ribosome depleted RNA from the 2-week old WT and Lmna-/- mice (N=5) . Strand specific RNA seq libraries where prepared form ribosome-depleted cardiac RNA samples using the Illumina TruSeq stranded total RNA library preparation kit. The samples weresequenced on the Illumina HiSeq 4000 instrument using the paired-end sequencing reagents to generate100 base paired end reads.
Project description:Mouse BM extracellular vesicles were isolated from Ocn-GFP Topaz mice using the Exoeasy and miRneasy (Qiagen). Small RNA libraries were constructed using NEBNext Multiplex Small RNA Library Prep Set for Illumina (New England biolabs). Libraries were sequenced using Hiseq2500 instrument
Project description:We report the application of MeRIP sequencing technology for high-throughput profiling of RNA m6A modifications in wide-type and knock-down METTL3 or WTAP Human Umbilical Vein Endothelial Cells (HUVECs). Both the input samples without immunoprecipitation and the m6A IP samples were used for RNA-seq library generation with NEBNext® Ultra II Directional RNA Library Prep Kit (New England Biolabs, Inc., USA). Library sequencing was performed on an illumina Hiseq instrument with 150bp paired-end reads. Clean reads of all libraries were aligned to the reference genome (HG19) by Hisat2 software (v2.0.4). Methylated sites on RNAs (peaks) were identified by MACS software. Differentially methylated sites were identified by diffReps. And, guided by the Ensembl gtf gene annotation file, cuffdiff software (part of cufflinks) was used to get the gene level FPKM as the expression profiles of mRNA, and fold change and p-value were calculated based on FPKM, differentially expressed mRNA were identified. qRT-PCR validation was performed using SYBR Green assays. Finally, we find that METTL3/WTAP can regulate the expression level of target genes through m6A modification in HUVECs. This study provides a framework for applying MeRIP sequencing profiles to characterize vascular endothelial cells.
Project description:Purpose: The aim of this study was to profile the transcriptome of differentiating embryonic stem cells (ES) with nonsilencing shRNA mediated knockdown and Setdb1 geneTrap heterozygous cells with Setdb1 shRNA mediated knockdown. Methods: Female Setdb1+/+ Xist∆A/+ ES cells were induced to differentiate and RNA sampled on day 0, day 3 and day 5 of differentiation. RNAseq libraries were prepared with TruSeq RNA sample preparation v2 kit. Libraries were pooled and sequenced on the Illumina HiSeq 2000 platform for 100 bp single-end reads. Image analysis was performed in real time by the HiSeq Control Software (HCS) v1.4.8 and Real Time Analysis (RTA) v1.12.4.2, running on the instrument computer. Real-time base calling on the HiSeq instrument computer was performed with the RTA software. Illumina CASAVA1.8 pipeline was used to generate the sequence data.
Project description:Hep3B and Huh7 cells pre-treated with XL413 for 10 days to induce senescence prior to sertraline treatment for 24 hours. For RNA sequencing, the library was prepared using TruSeq RNA sample prep kit according to the manufacturer’s protocol (Illumina). Gene set enrichment analysis was performed using gene set enrichment analysis software.
Project description:Purpose: to analyze gene expression in XL413 treated liver cancer cells. Methods: PLC/PRF/5 cells are treated with XL413 (10uM) for 96 hours .For RNA sequencing, the library was prepared using TruSeq RNA sample prep kit according to the manufacturer’s protocol (Illumina). Gene set enrichment analysis was performed using gene set enrichment analysis software.
Project description:The aim is to analyze the transcriptome profiling (RNA-seq) of Mettl3 knockout mesenchymal stem cells (MSCs) and identify m6A modified mRNAs (m6A-seq). For RNA-seq, total RNAs from MSCs of Prx1-Cre;Mettl3fl/fl mice and their control littermates were extracted and purified using poly-T oligo-attached magnetic beads. Sequencing libraries were generated using NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA) following manufacturer’s recommendations, and were then subjected to Illumina HiSeq 2500. We used FastQC (v0.11.5) and FASTX toolkit (0.0.13) to control the quality of RNA-seq data and mapped them to Mus musculus reference genomes using HISAT2 (v.2.0.4). Ballgown software (v.3.4.0) was performed to identify differentially expressed genes and transcripts. Genes were considered significantly differentially expressed if showing ≥ 2.0 fold change and P value < 0.05. For m6A-seq, total RNAs were isolated from MSCs and then mRNAs were purified using poly-T oligo-attached magnetic beads. The purified mRNAs were fragmentated and used for immunoprecipitation with anti-m6A antibody, then the pulled down fragments were purified and used for libraries construction with the Illumina TruSeq Stranded mRNA Sample Prep Kit following the manufacturer's instructions. m6A libraries were sequenced with Illumina NextSeq500. Sequencing reads were aligned to mm9 using tophat and peak calling was performed with MACS2.
Project description:Purpose: The aim of this study was to determine the DNA methylation state of wildtype female mouse embryonic fibroblasts with nonsilencing shRNA mediated knockdown and Setdb1 geneTrap heterozygous cells with Setdb1 shRNA mediated knockdown. Methods: Enhanced Reduced Representation Bisulfite Libraries (eRRBS) were produced as previously descirbed (Akalin et al. 2012). Libraries were pooled and sequenced on the Illumina HiSeq 2000 platform for 100 bp single-end reads with dark cycle parameters (Boyle et al. 2012). Image analysis was performed in real time by the HiSeq Control Software (HCS) v1.4.8 and Real Time Analysis (RTA) v1.12.4.2, running on the instrument computer. Real-time base calling on the HiSeq instrument computer was performed with the RTA software. Illumina CASAVA1.8 pipeline was used to generate the sequence data.