用生物信息学方法寻找肝癌特异性表达基因转录调控模式

为了对肝癌(hepatocellularcarcinoma,HCC)的分子发病机理进行研究,首先对肝癌基因表达谱数据用t-检验算法进行了分析,找到了肝癌中特异性表达基因(characteristicgenes).然后把这些基因结合已知的肝HNF家族转录因子染色质免疫共沉淀结合DNA启动子芯片(ChIP-chip)实验数据用SAEM算法进行分析,得到了肝癌特异性表达基因的转录调控关系,并寻找到了多个HNF家族转录因子调控单基因的转录调控模式.结果表明HNF家族转录因子对大量具有重要功能的肝癌特异性表达基因进行了转录调控,并且多个HNF家族转录因子调控单基因可以形成前馈环和多输入调控等模式.     

Genomic MRI - a Public Resource for Studying Sequence Patterns within Genomic DNA

Non-coding genomic regions in complex eukaryotes, including intergenic areas, introns, and untranslated segments of exons, are profoundly non-random in their nucleotide composition and consist of a complex mosaic of sequence patterns. These patterns include so-called Mid-Range Inhomogeneity (MRI) regions -- sequences 30-10000 nucleotides in length that are enriched by a particular base or combination of bases (e.g. (G+T)-rich, purine-rich, etc.). MRI regions are associated with unusual (non-B-form) DNA structures that are often involved in regulation of gene expression, recombination, and other genetic processes (Fedorova & Fedorov 2010). The existence of a strong fixation bias within MRI regions against mutations that tend to reduce their sequence inhomogeneity additionally supports the functionality and importance of these genomic sequences (Prakash et al. 2009).Here we demonstrate a freely available Internet resource -- the Genomic MRI program package -- designed for computational analysis of genomic sequences in order to find and characterize various MRI patterns within them (Bechtel et al. 2008). This package also allows generation of randomized sequences with various properties and level of correspondence to the natural input DNA sequences. The main goal of this resource is to facilitate examination of vast regions of non-coding DNA that are still scarcely investigated and await thorough exploration and recognition.     

Isolation of Ribosome Bound Nascent Polypeptides in vitro to Identify Translational Pause Sites Along mRNA

The rate of translational elongation is non-uniform. mRNA secondary structure, codon usage and mRNA associated proteins may alter ribosome movement on the message^{for review see 1}. However, it''s now widely accepted that synonymous codon usage is the primary cause of non-uniform translational elongation rates¹. Synonymous codons are not used with identical frequency. A bias exists in the use of synonymous codons with some codons used more frequently than others². Codon bias is organism as well as tissue specific^2,3. Moreover, frequency of codon usage is directly proportional to the concentrations of cognate tRNAs⁴. Thus, a frequently used codon will have higher multitude of corresponding tRNAs, which further implies that a frequent codon will be translated faster than an infrequent one. Thus, regions on mRNA enriched in rare codons (potential pause sites) will as a rule slow down ribosome movement on the message and cause accumulation of nascent peptides of the respective sizes^5-8. These pause sites can have functional impact on the protein expression, mRNA stability and protein folding^{for review see 9}. Indeed, it was shown that alleviation of such pause sites can alter ribosome movement on mRNA and subsequently may affect the efficiency of co-translational (in vivo) protein folding^1,7,10,11. To understand the process of protein folding in vivo, in the cell, that is ultimately coupled to the process of protein synthesis it is essential to gain comprehensive insights into the impact of codon usage/tRNA content on the movement of ribosomes along mRNA during translational elongation.Here we describe a simple technique that can be used to locate major translation pause sites for a given mRNA translated in various cell-free systems^6-8. This procedure is based on isolation of nascent polypeptides accumulating on ribosomes during in vitro translation of a target mRNA. The rationale is that at low-frequency codons, the increase in the residence time of the ribosomes results in increased amounts of nascent peptides of the corresponding sizes. In vitro transcribed mRNA is used for in vitro translational reactions in the presence of radioactively labeled amino acids to allow the detection of the nascent chains. In order to isolate ribosome bound nascent polypeptide complexes the translation reaction is layered on top of 30% glycerol solution followed by centrifugation. Nascent polypeptides in polysomal pellet are further treated with ribonuclease A and resolved by SDS PAGE. This technique can be potentially used for any protein and allows analysis of ribosome movement along mRNA and the detection of the major pause sites. Additionally, this protocol can be adapted to study factors and conditions that can alter ribosome movement and thus potentially can also alter the function/conformation of the protein.     

羊FSHR基因5′端转录启动调控区生物学特性

文章对小尾寒羊、滩羊和澳洲绵羊等繁殖性状不同的3种绵羊与排卵有关的FSHR基因5′端转录启动调控区进行了克隆和分析,通过对FSHR基因的15个转录调控元件序列进行比较,结果表明,羊不同品种FSHR基因的转录调控元件序列之间没有差异。这说明绵羊的品种与FSHR基因5′端转录启动调控区的相关性不强,排除了因转录调控元件突变而影响转录调节能力的可能性。      

Gene Expression in Bacteria Directed by Plant-specific Regulatory Sequences

The regulation of gene expression represents a specific process which has different structural and functional requirements in different groups of organisms. It is thus assumed that regulatory sequences of eucaryotes cannot be recognized in procaryotes. This assumption is of interest for risk assessments of the environmental impact of deliberate release experiments with genetically modified organisms. In order to analyse the extent of heterologous gene expression caused by the transfer of plant-specific regulatory sequences into bacteria, we constructed fusions between plant-specific regulatory sequences and the coding regions of the luxAB genes for the luciferase of the bioluminescent bacterium Vibrio harveyi, transferred the fusions into different bacterial species and measured the luminescence to quantify the expression of the luciferase genes. The regulatory sequences investigated included (a) the 35S promoter of the Cauliflower mosaic virus, (b) the B33 promoter of a class I patatin gene of potatoes, (c) the promoter of the ST-LS1 gene of potatoes and (d) the promoter of the rolC gene of Agrobacterium rhizogenes. We could show that in addition to the 35S promoter, which has already been described as being recognized in Escherichia coli, the sequences containing the B33 and the ST-LS1 promoters are recognized in bacteria. Luciferase gene expression promoted by the sequence with the ST-LS1 promoter could be observed in E. coli, Yersinia enterocolitica and Agrobacterium tumefaciens. Comparison of the luminescence caused by fusions between luxAB and different promoters on the chromosome and on an endogenous plasmid of Y. enterocolitica demonstrated that the level of the heterologous gene expression caused by the fragment with the ST-LS1 promoter was within the range of gene expression levels caused by endogenous promoters of Y. enterocolitica.     

人管家基因启动子序列中的组合转录调控元件分析

研究表明,第一内含子可能参与基因转录调控.利用统计方法提取人管家基因上游至第一内含子序列中潜在的组合转录调控模体,分析模体间的距离、区域分布等特征,探讨内含子参与基因转录调控的可能性及其参与方式.在管家基因中共获得960对潜在转录调控模体对,其中57%与实验已知的具有转录相互作用的因子对吻合,共涉及12组因子对.分析发现,绝大多数模体对(80%)偏向于上游区域及"上游-内含子"区域,进一步支持了内含子参与基因转录调控的假设,并据此推测内含子与上游序列之间具有转录协同作用,模体在基因转录起始位点(TSS)附近较为集中,模体对的两个模体之间距离较近,60%左右距离在200 bp以内,特别地,65%的模体对特征距离在100 bp以内,短距离间隔有利于转录因子间的协同作用.这些结果将有助于对人基因转录调控机制及内含子功能的深入认识.     

基于对数线性模型的酵母基因转录调控模体分析

识别真核基因的转录因子结合位点(或称模体)是后基因组时代的一项主要工作,对共表达或共调控的基因同时进行分析可以提高模体识别的准确性.本文基于2×2列联表的对数线性模型,以模体出现的基因条数计数,对酵母核糖体蛋白(RP)基因普遍使用的转录调控模体进行分析,然后用U-检验进一步筛选出相对于背景序列来说过表达的模体.这些模体为酵母RP基因潜在的转录调控元件,与实验获得的转录因子结合位点的符合率达90%.本方法的优点在于用严格的统计标准在一组基因启动子中搜索普遍使用的模体,克服了以往分析中对模体使用普遍性的模糊判断.本文的方法也可以有效地搜索共表达基因族的组合调控模体对.研究中还发现一个现象:2×2列联表中反映属性相关程度的Pearson相关系数与对数线性模型的交互效应之间存在着明显的相关性.这一结果提示,可以用对数线性模型的交互效应来评价两属性的关联情况.     

Integrating Mouse and Human Genetic Data to Move beyond GWAS and Identify Causal Genes in Cholesterol Metabolism

1. Download : Download high-res image (160KB)
2. Download : Download full-size image

     

Omics Approaches to Identify Potential Biomarkers of Inflammatory Diseases in the Focal Adhesion Complex

Inflammatory diseases such as inflammatory bowel disease(IBD) require recurrent invasive tests, including blood tests, radiology, and endoscopic evaluation both to diagnose and assess disease activity, and to determine optimal therapeutic strategies. Simple ‘bedside' biomarkers could be used in all phases of patient management to avoid unnecessary investigation and guide further management. The focal adhesion complex(FAC) has been implicated in the pathogenesis of multiple inflammatory diseases, including IBD, rheumatoid arthritis, and multiple sclerosis. Utilizing omics technologies has proven to be an efficient approach to identify biomarkers from within the FAC in the field of cancer medicine. Predictive biomarkers are paving the way for the success of precision medicine for cancer patients, but inflammatory diseases have lagged behind in this respect. This review explores the current status of biomarker prediction for inflammatory diseases from within the FAC using omics technologies and highlights the benefits of future potential biomarker identification approaches.     

MotViz: a tool for sequence motif prediction in parallel to structural visualization and analyses

Linking similar proteins structurally is a challenging task that may help in finding the novel members of a protein family. In this respect, identification of conserved sequence can facilitate understanding and classifying the exact role of proteins. However, the exact role of these conserved elements cannot be elucidated without structural and physiochemical information. In this work, we present a novel desktop application MotViz designed for searching and analyzing the conserved sequence segments within protein structure. With MotViz, the user can extract a complete list of sequence motifs from loaded 3D structures, annotate the motifs structurally and analyze their physiochemical properties. The conservation value calculated for an individual motif can be visualized graphically. To check the efficiency, predicted motifs from the data sets of 9 protein families were analyzed and MotViz algorithm was more efficient in comparison to other online motif prediction tools. Furthermore, a database was also integrated for storing, retrieving and performing the detailed functional annotation studies. In summary, MotViz effectively predicts motifs with high sensitivity and simultaneously visualizes them into 3D strucures. Moreover, MotViz is user-friendly with optimized graphical parameters and better processing speed due to the inclusion of a database at the back end. MotViz is available at http://www.fi-pk.com/motviz.html.     

Detection of Histone Modifications in Plant Leaves

Chromatin structure is important for the regulation of gene expression in eukaryotes. In this process, chromatin remodeling, DNA methylation, and covalent modifications on the amino-terminal tails of histones H3 and H4 play essential roles^1-2. H3 and H4 histone modifications include methylation of lysine and arginine, acetylation of lysine, and phosphorylation of serine residues^1-2. These modifications are associated either with gene activation, repression, or a primed state of gene that supports more rapid and robust activation of expression after perception of appropriate signals (microbe-associated molecular patterns, light, hormones, etc.)^3-7. Here, we present a method for the reliable and sensitive detection of specific chromatin modifications on selected plant genes. The technique is based on the crosslinking of (modified) histones and DNA with formaldehyde^8,9, extraction and sonication of chromatin, chromatin immunoprecipitation (ChIP) with modification-specific antibodies^9,10, de-crosslinking of histone-DNA complexes, and gene-specific real-time quantitative PCR. The approach has proven useful for detecting specific histone modifications associated with C₄ photosynthesis in maize^5,11 and systemic immunity in Arabidopsis³.     

Using Human Microarrays to Identify Differentially Expressed Genes Associated with Increased Steroidogenesis in Boars

Human microarrays are readily available, and it would be advantageous if they could be used to study gene expression in other species, such as pigs. The objectives of this research were to validate the use of human microarrays in the analysis of porcine gene expression, to assess the variability of the data generated, and to compare gene expression in boars with different levels of steroidogenesis. Cytochrome b5 (CYB5) expression was used to assess array detection sensitivity. Samples having high or low CYB5 RNA levels were hybridized to microarrays to determine if the known expression difference could be detected. Six hybridizations were conducted using human microarrays containing 3840 total spots representing 1718 characterized human ESTs. To analyze gene expression in boars with different levels of steroidogenesis, testis RNA from four boars with high levels of plasma estrone sulphate was hybridized to testis RNA from four boars with lower levels. Eight microarray hybridizations were conducted including fluor-flips. Self-self hybridizations were also conducted to assess the variability of array experiments. The Cy5 and Cy3 intensity values for each array were normalized using a locally weighted linear regression (LOESS). Statistical significance was assessed using a Student's t-test followed by the Benjamini and Hochberg multiple testing correction procedure. Quantitative real-time PCR (Q-RT-PCR) was used to verify select gene expression differences. The results show that CYB5 was significantly overexpressed in the high CYB5 sample by 1.8 fold (P < 0.05), verifying the known expression difference. The average log2 ratio of the majority of genes (1643) falls within one standard deviation of the mean, indicating the data were reproducible. In the high versus low steroidogenesis experiment, seven genes were significantly overexpressed in the high group (P < 0.05). Quantitative real-time PCR was used to validate five genes with the highest fold change, and the results corroborated those found by the microarray experiments. The results of the self-self hybridizations showed that no genes were significantly differentially expressed following the application of the Benjamini and Hochberg multiple testing correction procedure. The results presented in this report show that human arrays can be used for gene expression analysis in pigs.     

核糖体蛋白基因上游转录调控元件协同作用的分析

在对酵母基因上游区调控位点的研究中发现,具有高转录水平的基因几乎都是编码核糖体蛋白的,在低转录水平的基因中却没有发现核糖体蛋白基因,这一现象似乎提示着核糖体蛋白基因上游序列中含有较多潜在的有利于基因转录的转录因子结合位点.为了能对核糖体蛋白基因上游序列有一个系统全面的认识,本文采用多种统计方法从不同角度探测了核糖体蛋白基因上游序列中潜在的转录正调控位点,并对它们之间可能存在的协同作用模式进行了分析.     

应用绿色荧光蛋白研究外源基因在造血细胞的表达调控

以增强型绿色荧光蛋白为报道分子,研究双基因逆转录病毒载体介导小鼠骨髓细胞的基因转移中,SV(simianvirus)40启动子和内部核糖体进入位点(internalribosomeentrysite,IRES)对双基因共表达的影响。构建双基因中间序列分别为SV40启动子和IRES的载体pLESN和pLEIN。经包装获得较高滴度的病毒上清,以共培养的方式感染5-氟尿嘧啶预刺激的小鼠骨髓细胞。流式细胞仪检测表明转染效率约25%,PCR证明EGFP基因整合至骨髓细胞基因组。半固体培养转基因细胞7d,LEIN组获得具有G418抗性的集落中98%表达绿色荧光,而LESN组54%。结果表明:在双基因逆转录病毒载体介导的小鼠骨髓细胞的基因转导中,IRES与内部SV40启动子相比,更能保证双基因的共同表达。     

Using RNA-mediated Interference Feeding Strategy to Screen for Genes Involved in Body Size Regulation in the Nematode C. elegans

Double-strand RNA-mediated interference (RNAi) is an effective strategy to knock down target gene expression^1-3. It has been applied to many model systems including plants, invertebrates and vertebrates. There are various methods to achieve RNAi in vivo^4,5. For example, the target gene may be transformed into an RNAi vector, and then either permanently or transiently transformed into cell lines or primary cells to achieve gene knockdown effects; alternatively synthesized double-strand oligonucleotides from specific target genes (RNAi oligos) may be transiently transformed into cell lines or primary cells to silence target genes; or synthesized double-strand RNA molecules may be microinjected into an organism. Since the nematode C. elegans uses bacteria as a food source, feeding the animals with bacteria expressing double-strand RNA against target genes provides a viable strategy⁶. Here we present an RNAi feeding method to score body size phenotype. Body size in C. elegans is regulated primarily by the TGF- β - like ligand DBL-1, so this assay is appropriate for identification of TGF-β signaling components⁷. We used different strains including two RNAi hypersensitive strains to repeat the RNAi feeding experiments. Our results showed that rrf-3 strain gave us the best expected RNAi phenotype. The method is easy to perform, reproducible, and easily quantified. Furthermore, our protocol minimizes the use of specialized equipment, so it is suitable for smaller laboratories or those at predominantly undergraduate institutions.     

Deciphering Ascorbic Acid Regulatory Pathways in Ripening Tomato Fruit Using a Weighted Gene Correlation Network Analysis Approach

Genotype is generally determined by the co-expression of diverse genes and multiple regulatory pathways in plants. Gene co-expression analysis combining with physiological trait data provides very important information about the gene function and regulatory mechanism. L-Ascorbic acid （AsA）, which is an essential nutrient component for human health and plant metabolism, plays key roles in diverse biological processes such as cell cycle, cell expansion, stress resistance, hormone synthesis, and signaling. Here, we applied a weighted gene correlation network analysis approach based on gene expression values and AsA content data in ripening tomato （Solanum lycopersicum L.） fruit with different AsA content levels, which leads to identification of AsA relevant modules and vital genes in AsA regulatory pathways. Twenty- four modules were compartmentalized according to gene expression profiling. Among these modules, one negatively related module containing genes involved in redox processes and one positively related module enriched with genes involved in AsA biosynthetic and recycling pathways were further analyzed. The present work herein indicates that redox pathways as well as hormone-signal pathways are closely correlated with AsA accumulation in ripening tomato fruit, and allowed us to prioritize candidate genes for follow-up studies to dissect this interplay at the biochemical and molecular level.