Evolution of conserved secondary structures and their function in transcriptional regulation networks

Background: Many conserved secondary structures have been identified within conserved elements in the human genome, but only a small fraction of them are known to be functional RNAs. The evolutionary variations of these conserved secondary structures in h

Single molecule dynamics and statistical fluctuations of gene regulatory networks: A repressilator

The authors developed a time dependent method to study the single molecule dynamics of a simple gene regulatory network: a repressilator with three genes mutually repressing each other. They quantitatively characterize the time evolution dynamics of the repressilator. Furthermore, they study purely dynamical issues such as statistical fluctuations and noise evolution. They illustrated some important features of the biological network such as monostability, spirals, and limit cycle oscillation. Explicit time dependent Fano factors which describe noise evolution and show statistical fluctuations out of equilibrium can be significant and far from the Poisson distribution. They explore the phase space and the interrelationships among fluctuations, order, amplitude, and period of oscillations of the repressilators. The authors found that repressilators follow ordered limit cycle orbits and are more likely to appear in the lower fluctuating regions. The amplitude of the repressilators increases as the suppressing of the genes decreases and production of proteins increases. The oscillation period of the repressilators decreases as the suppressing of the genes decreases and production of proteins increases.

Detection of potential positive regulatory motifs of transcription in yeast introns by comparative analysis of oligonucleotide frequencies

We conducted a comparative statistical analysis of tetra- through hexanucleotide frequencies in two sets of introns of yeast genes. The first set consisted of introns of genes that have transcription rates higher than 30 mRNAs/h while the second set contained introns of genes whose transcription rates were lower than or equal to 10 mRNAs/h. Some oligonucleotides whose occurrence frequencies in the first set of introns are significantly higher than those in the second set of introns were detected. The frequencies of occurrence of most of these detected oligonucleotides are also significantly higher than those in the exons flanking the introns of the first set. Interestingly some of these detected oligonucleotides are the same as well known "signature" sequences of transcriptional regulatory elements. This could imply the existence of potential positive regulatory motifs of transcription in yeast introns. (C) 2003 Elsevier Ltd. All rights reserved.

Gene transfer into genomes of human cells by the sleeping beauty transposon system

The Sleeping Beauty (SB) transposon system, derived from teleost fish sequences, is extremely effective at delivering DNA to vertebrate genomes, including those of humans. We have examined several parameters of the SB system to improve it as a potential, nonviral vector for gene therapy. Our investigation centered on three features: the carrying capacity of the transposon for efficient integration into chromosomes of HeLa cells, the effects of overexpression of the SB transposase gene on transposition rates, and improvements in the activity of SB transposase to increase insertion rates of transgenes into cellular chromosomes. We found that SB transposons of about 6 kb retained 50% of the maximal efficiency of transposition, which is sufficient to deliver 70-80% of identified human cDNAs with appropriate transcriptional regulatory sequences. Overexpression inhibition studies revealed that there are optimal ratios of SB transposase to transposon for maximal rates of transposition, suggesting that conditions of delivery of the two-part transposon system are important for the best gene-transfer efficiencies. We further refined the SB transposase to incorporate several amino acid substitutions, the result of which led to an improved transposase called SB11. With SB11 we are able to achieve transposition rates that are about 100-fold above those achieved with plasmids that insert into chromosomes by random recombination. With the recently described improvements to the transposon itself, the SB system appears to be a potential gene-transfer tool for human gene therapy.

人类基因组中保守二级结构的纯净化选择及其在转录调控网络中的作用

保守序列是一种跨物种保守的基因组序列，而且绝大多数为非蛋白编码序 列。保守序列在人类遗传疾病中发挥着重要作用。其中，一部分保守序列能够 折叠形成二级结构。已鉴定的一些保守二级结构编码一些RNA 分子，如 microRNA、RNA 编辑序列和组蛋白mRNA 3’端非翻译区茎环结构等。但是，对 于绝大部分的保守二级结构，它们的生物学功能以及作用于它们上面的进化作 用力依然是未知的。 群体的SNP 数据在分析序列上的进化作用力时非常有效。SNP 在群体中的 频率会因为受到不同的进化作用力而表现出差异，而与其是否位于基因组中的 突变热点无关。对于受纯净化选择作用的SNP，它们的频率一般会比中性SNP 具有低的新生型等位基因频率（DAF）。我们运用生物信息学的方法，在人类基 因组保守二级结构中找到746 个SNP。这746 个SNP 与基因组其它区段的SNP 在突变模式上并不存在显著差异，在保守二级结构内同样存在突变热点。通过 与侧翼序列SNP 的分布比较发现，保守二级结构上SNP 密度约为其侧翼序列的 2/3。相比于侧翼序列SNP，有更高比例的保守二级结构SNP 具有低的DAF 值。 这些结果提示，有很多保守二级结构上的SNP 因为受到纯净化选择作用而在现 代人群中被剔除了。保守二级结构与侧翼序列在SNP 密度和DAF 上的差异要高 于保守序列与非保守序列之间的差异，提示保守二级结构是受到纯净化选择作 用最为严格的一类保守序列。我们发现，在保守二级结构内部，纯净化选择作 用的强度也有差异。茎区比环区具有更低的SNP 密度，而且有更高比例的茎区 SNP 具有低的DAF 值。这个结果提示，保守二级结构上的纯净化选择力主要作 用于茎区上的位点。我们推测，这可能是茎区上的突变往往比环区的突变对二级结构的造成更大的影响导致的。 我们通过寻找保守二级结构与转录因子SOX2、OCT4、NANOG、SUZ12 和C-MYC 结合位点之间的重叠，还分析了保守二级结构在转录调控网络中的作用。结果 显示，很多保守二级结构是作为转录因子的结合位点调控了许多与发育相关的 转录因子编码基因的表达。转录因子与保守二级结构之间的结合模式非常复杂， 可以有多个转录因子结合到同一个保守二级结构上，也可以是一个转录因子结 合到自身编码基因相关的保守二级结构上。不同的转录因子和保守二级结构结 合可以主导靶基因的特异模式，当绝大多数相关的保守二级结构与SUZ12 结合 时，基因表达受到抑制，而当绝大多数相关的保守二级结构不与SUZ12 结合时， 基因表达受到激活。在转录调控网络中，约有30%的保守二级结构是作为启动 子来调控基因的表达。因为转录因子SOX2、OCT4、NANOG、SUZ12 和C-MYC 仅仅 只结合到很小一部分保守二级结构上，提示可能还有更多的转录因子会结合到 保守二级结构上。因此，保守二级结构介导的转录调控网络要比目前已知的复 杂得多。

Statistical analysis of sequence features of introns with positive transcriptional regulation in yeast genes

A great deal of experimental studies have shown that many introns of eukaryotic genes function as regulators of transcription. However, comprehensive studies of this problem have not yet been conducted. After checking the transcription frequencies of some Saccharomyces cerevisiae (yeast), genes and their introns, a remarkable phenomenon was discovered that generally the introns of the genes with higher transcription frequencies are longer, and the introns of the genes with lower transcription frequencies are shorter. This suggests that the longer introns of genes with higher transcription frequencies may contain some characteristic sequence structures, which could enhance the transcription of genes. Therefore, two sets of introns of yeast genes were chosen for further study. The transcription frequencies of the first set of genes are higher (>30), and those of the second set of genes are lower (less than or equal to10). Some oligonucleotides are detected by statistically comparative analyses of the occurrence frequencies of oligonucleotides (mainly tetranucleotides and pentanucleotides), whose occurrence frequencies in the first set of introns; are significantly higher than those in the second set of introns, and are also significantly higher than those in the exons flanking the introns of the first set. Some of these extracted oligonucleotides are the same as the regulatory elements of transcription revealed by experimental analyses. Besides, the distributions of these extracted oligonucleotides in the two sets of introns and the exons show that the sequence structures of the first set of introns are favorable for transcription of genes.

Potential energy landscape and robustness of a gene regulatory network: Toggle switch

Finding a multidimensional potential landscape is the key for addressing important global issues, such as the robustness of cellular networks. We have uncovered the underlying potential energy landscape of a simple gene regulatory network: a toggle switch. This was realized by explicitly constructing the steady state probability of the gene switch in the protein concentration space in the presence of the intrinsic statistical fluctuations due to the small number of proteins in the cell. We explored the global phase space for the system. We found that the protein synthesis rate and the unbinding rate of proteins to the gene were small relative to the protein degradation rate; the gene switch is monostable with only one stable basin of attraction. When both the protein synthesis rate and the unbinding rate of proteins to the gene are large compared with the protein degradation rate, two global basins of attraction emerge for a toggle switch. These basins correspond to the biologically stable functional states. The potential energy barrier between the two basins determines the time scale of conversion from one to the other. We found as the protein synthesis rate and protein unbinding rate to the gene relative to the protein degradation rate became larger, the potential energy barrier became larger. This also corresponded to systems with less noise or the fluctuations on the protein numbers.