Transcriptomic Analyses of the Biological Effects of Airborne PM2.5 Exposure on Human Bronchial Epithelial Cells

Epidemiological studies have associated high levels of airborne particulate matter (PM) with increased respiratory diseases. In order to investigate the mechanisms of air pollution-induced lung toxicity in humans, human bronchial epithelial cells (16HBE) were exposed to various concentrations of particles smaller than 2.5 μm (PM2.5) collected from Beijing, China. After observing that PM2.5 decreased cell viability in a dose-dependent manner, we first used Illumina RNA-seq to identify genes and pathways that may contribute to PM2.5-induced toxicity to 16HBE cells. A total of 539 genes, 283 up-regulated and 256 down-regulated, were identified to be significantly differentially expressed after exposure to 25 μg/cm² PM2.5. PM2.5 induced a large number of genes involved in responses to xenobtiotic stimuli, metabolic response, and inflammatory and immune response pathways such as MAPK signaling and cytokine-cytokine receptor interaction, which might contribute to PM2.5-related pulmonary diseases. We then confirmed our RNA-seq results by qPCR and by analysis of IL-6, CYP1A1, and IL-8 protein expression. Finally, ELISA assay demonstrated a significant association between exposure to PM2.5 and secretion of IL-6. This research provides a new insight into the mechanisms underlying PM2.5-induced respiratory diseases in Beijing.     

Identification of differentially expressed genes and the role of PDK4 in CD14+ monocytes of coronary artery disease

Background. Coronary artery disease (CAD) is a chronic inflammatory disease caused by development of atherosclerosis (AS), which is the leading cause of mortality and disability. Our study aimed to identify the differentially expressed genes (DEGs) in CD14⁺ monocytes from CAD patients compared with those from non-CAD controls, which might pave the way to diagnosis and treatment for CAD. Methods. The RNA-sequencing (RNA-seq) was performed by BGISEQ-500, followed by analyzing with R package to screening DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by R package. In addition, we validated the results of RNA-seq using real-time quantitative polymerase chain reaction (RT-qPCR). Furthermore, we explored the function of selected ten genes in LDL-treated CD14⁺ monocytes by RT-qPCR. Results. a total of 2897 DEGs were identified, including 753 up- and 2144 down-regulated genes in CD14⁺ monocytes from CAD patients. These DEGs were mainly enriched in plasma membrane and cell periphery of cell component, immune system process of biological process, NF-κB signaling pathway, cell adhesion molecules signaling pathway and cytokine–cytokine receptor interaction signaling pathway. In LDL-treated CD14⁺ monocytes, the mRNA expression of pyruvate dehydrogenase kinase 4 (PDK4) was significantly up-regulated. Conclusion. In the present study, we suggested that PDK4 might play a role in progression of CAD. The study will provide some pieces of evidence to investigate the role and mechanism of key genes in the pathogenesis of CAD.     

Global Analysis of Differentially Expressed Genes and Proteins in the Wheat Callus Infected by Agrobacterium tumefaciens

Agrobacterium-mediated plant transformation is an extremely complex and evolved process involving genetic determinants of both the bacteria and the host plant cells. However, the mechanism of the determinants remains obscure, especially in some cereal crops such as wheat, which is recalcitrant for Agrobacterium-mediated transformation. In this study, differentially expressed genes (DEGs) and differentially expressed proteins (DEPs) were analyzed in wheat callus cells co-cultured with Agrobacterium by using RNA sequencing (RNA-seq) and two-dimensional electrophoresis (2-DE) in conjunction with mass spectrometry (MS). A set of 4,889 DEGs and 90 DEPs were identified, respectively. Most of them are related to metabolism, chromatin assembly or disassembly and immune defense. After comparative analysis, 24 of the 90 DEPs were detected in RNA-seq and proteomics datasets simultaneously. In addition, real-time RT-PCR experiments were performed to check the differential expression of the 24 genes, and the results were consistent with the RNA-seq data. According to gene ontology (GO) analysis, we found that a big part of these differentially expressed genes were related to the process of stress or immunity response. Several putative determinants and candidate effectors responsive to Agrobacterium mediated transformation of wheat cells were discussed. We speculate that some of these genes are possibly related to Agrobacterium infection. Our results will help to understand the interaction between Agrobacterium and host cells, and may facilitate developing efficient transformation strategies in cereal crops.     

生物膜形成中间状态下副溶血弧菌的基因转录谱分析

【目的】探究生物膜形成中间状态下副溶血弧菌的差异基因表达情况,为今后研究生物膜形成调控机制提供基因信息。【方法】以非生物膜形成条件下为参照,采用Illumina HiSeq测序平台进行转录组测序（RNA sequencing,RNA-seq）研究,分析生物膜形成中间状态下副溶血弧菌的基因表达情况,并采用实时定量PCR （quantitative real-time PCR,qPCR）进行验证。【结果】本研究共获得979个差异显著性表达基因（differentially expressed gene,DEG）,其中下调基因379个,上调基因600个。基因本体（gene ontology,GO）分类分析结果显示,共有363个DEGs注释到分子功能、细胞组分和生物学过程3个一级分类和30个二级分类;京都基因和基因组百科全书（Kyoto encyclopedia of genes and genomes,KEGG）代谢途径分析结果显示,共有706个DEGs归到37个代谢通路中（Q value<0.05）,差异表达基因主要集中在细胞代谢和转运通路上;蛋白相邻类的聚簇（clusters of orthologous groups,COG）分类结果显示,有888个DEGs可归为20个类别,涉及DEGs最多的为氨基酸转运与代谢、一般功能预测基因、能量产生与转换以及未知功能基因。qPCR验证挑选的DEGs变化趋势均与RNA-seq的结果一致。此外,从转录组数据中共筛选出10个c-di-GMP代谢相关基因、1个侧生鞭毛蛋白基因（lafA）、13个极生鞭毛合成相关基因、6个荚膜多糖合成相关基因、6个胞外多糖合成相关基因、5个IV型菌毛合成相关基因、膜融合蛋白（membrane fusion protein,Mfp）基因（cpsQ-mfpABC）、14个III型分泌系统1（T3SS1）相关基因、14个Vp-PAI基因（1个tdh2和13个T3SS2基因）、3个VI型分泌系统1（T6SS1）相关基因、6个T6SS2基因、45个推定调控子基因和15个推定的外膜蛋白基因。【结论】生物膜形成引起副溶血弧菌基因表达谱发生明显变化,差异表达基因中包含生物膜形成相关基因、关键毒力基因和许多推定调控子基因等,这为后续研究生物膜形成调控机制提供重要信息。     

A normalization strategy for comparing tag count data

Background

High-throughput sequencing, such as ribonucleic acid sequencing (RNA-seq) and chromatin immunoprecipitation sequencing (ChIP-seq) analyses, enables various features of organisms to be compared through tag counts. Recent studies have demonstrated that the normalization step for RNA-seq data is critical for a more accurate subsequent analysis of differential gene expression. Development of a more robust normalization method is desirable for identifying the true difference in tag count data.

Results

We describe a strategy for normalizing tag count data, focusing on RNA-seq. The key concept is to remove data assigned as potential differentially expressed genes (DEGs) before calculating the normalization factor. Several R packages for identifying DEGs are currently available, and each package uses its own normalization method and gene ranking algorithm. We compared a total of eight package combinations: four R packages (edgeR, DESeq, baySeq, and NBPSeq) with their default normalization settings and with our normalization strategy. Many synthetic datasets under various scenarios were evaluated on the basis of the area under the curve (AUC) as a measure for both sensitivity and specificity. We found that packages using our strategy in the data normalization step overall performed well. This result was also observed for a real experimental dataset.

Conclusion

Our results showed that the elimination of potential DEGs is essential for more accurate normalization of RNA-seq data. The concept of this normalization strategy can widely be applied to other types of tag count data and to microarray data.     

Reduced CYFIP1 in Human Neural Progenitors Results in Dysregulation of Schizophrenia and Epilepsy Gene Networks

Deletions encompassing the BP1-2 region at 15q11.2 increase schizophrenia and epilepsy risk, but only some carriers have either disorder. To investigate the role of CYFIP1, a gene within the region, we performed knockdown experiments in human neural progenitors derived from donors with 2 copies of each gene at the BP1-2 locus. RNA-seq and cellular assays determined that knockdown of CYFIP1 compromised cytoskeletal remodeling. FMRP targets and postsynaptic density genes, each implicated in schizophrenia, were significantly overrepresented among differentially expressed genes (DEGs). Schizophrenia and/or epilepsy genes, but not those associated with randomly selected disorders, were likewise significantly overrepresented. Mirroring the variable expressivity seen in deletion carriers, marked between-line differences were observed for dysregulation of disease genes. Finally, a subset of DEGs showed a striking similarity to known epilepsy genes and represents novel disease candidates. Results support a role for CYFIP1 in disease and demonstrate that disease-related biological signatures are apparent prior to neuronal differentiation.     

Comprehensive analysis of coding sequence architecture features and gene expression in Arachis duranensis

Coding sequence (CDS) architecture affects gene expression levels in organisms. Codon optimization can increase the gene expression level. Therefore, understanding codon usage patterns has important implications for research on genetic engineering and exogenous gene expression. To date, the codon usage patterns of many model plants have been analyzed. However, the relationship between CDS architecture and gene expression in Arachis duranensis remains poorly understood. According to the results of genome sequencing, A. duranensis has many resistant genes that can be used to improve the cultivated peanut. In this study, bioinformatic approaches were used to estimate A. duranensis CDS architectures, including frequency of the optimal codon (Fop), polypeptide length and GC contents at the first (GC1), second (GC2) and third (GC3) codon positions. In addition, Arachis RNA-seq datasets were downloaded from PeanutBase. The relationships between gene expression and CDS architecture were assessed both under normal growth as well as nematode and drought stress conditions. A total of 26 codons with high frequency were identified, which preferentially ended with A or T in A. duranensis CDSs under the above-mentioned three conditions. A similar CDS architecture was found in differentially expressed genes (DEGs) under nematode and drought stresses. The GC1 content differed between DEGs and non-differentially expressed genes (NDEGs) under both drought and nematode stresses. The expression levels of DEGs were affected by different CDS architectures compared with NDEGs under drought stress. In addition, no correlation was found between differential gene expression and CDS architecture neither under nematode nor under drought stress. These results aid the understanding of gene expression in A. duranensis.     

Exploring the correlation between deltamethrin stress and Keap1-Nrf2-ARE pathway from Drosophila melanogaster RNASeq data

Deltamethrin (DM) is widely used in a variety of pest control, resulting in serious drug resistance. Keap1-Nrf2-ARE is the antioxidant stress pathway. We identified 268 genes differentially expressed (DEGs) in Drosophila Kc cells treated with DM, including up-regulated 180 genes and down-regulated 88 genes compared with the control group (fold-change≥2, qValue≤0.001) by RNA-seq, they are mainly linked to metabolic process, stimulation response, immune system process. When the cells are treated with DM in the case of overexpression of the Keap1 gene, the cytochrome P450 family genes were significantly down-regulated, and some diseases-related genes and non-coding genes also changed. Our data shown that Keap1-Nrf2-ARE pathway may play an important role in DM stress, which will provide a new direction for studying the mechanism of insect resistance.