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

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.

Sox2基因3'非翻译区保守元件对基因表达的调控作用

转录因子Sox2是脊椎动物早期发育中最早表达的神经系统特异性基因之一,同时在干细胞的维持中也起着关键作用.通过生物信息学分析,作者发现在脊椎动物Sox2 mRNA 3'非翻译区中存在4段非常保守的富含AU的区域.将这些片段按照不同的组合克隆到GFP和荧光素酶两种报告基因载体中,在非洲爪蟾胚胎和培养细胞中检测了这些片段对报告基因表达的影响.结果显示,Sox2的3'UTR可影响报告基因的表达水平,特别是其中的保守片段2可显著提高报告基因的表达水平,表明Sox2 3'非翻译区有可能参与Sox2表达的转录后调控.

Differential expression of the Brunol/CELF family genes during Xenopus laevis early development

The BRUNOL/CELF family of RNA-binding proteins plays important roles in post-transcriptional regulation and has been implicated in several developmental processes. In this study, we describe the cloning and expression patterns of five Brunol genes in Xenopus laevis. Among them, only Brunol2 is maternally expressed and the zygotic expression of the other four Brunol genes starts at different developmental stages. During Xenopus development, Brunol1, 4-5 are exclusively expressed in the nervous system including domains in the brain, spinal cord, optic and otic vesicles. Brunol2 and 3 are expressed in both the somatic mesoderm and the nervous system. Brunol2 is also extensively expressed in the lens. In transfected Hela cells, BRUNOL1, 2 and 3 proteins are localized in both the cytoplasm and the nucleus, while BRUNOL4 and 5 are only present in the cytoplasm, indicating their different functions.

非洲爪蟾早期发育中的转录后调控机制及眼睛发育相关基因的表达研究

近来的研究表明，转录后调控对于调节脊椎动物发育过程中的细胞分化，细胞分裂及基因区域特异性表达都具有重要作用。转录后调控包括对mRNA稳定性、翻译效率、细胞内定位及poly(A)水平的调控等。Sox2基因是脊椎动物早期发育中最早表达的神经系统特异性基因之一，是脊椎动物早期神经系统发育的重要调节因子。通过生物信息学分析，我们发现，在脊椎动物Sox2 mRNA 3’非翻译区中存在4段非常保守的富含AU的区域，通过报告基因分析等手段研究发现Sox2 3’非翻译区中的部分元件可显著提高报告基因表达，提示我们Sox2的表达可能受到转录后调控。 我们通过对爪蟾Xfhl3基因序列分析时发现其3’非翻译区存在一段保守的只在两栖类具有的序列，我们克隆并检测了该基因的表达图式，并采用报告基因分析等手段研究了Xfhl3基因 3’非翻译区对报告基因表达的影响。结果发现其3’UTR可抑制报告基因表达水平。由于Xfhl3基因3’UTR中这段序列只在爪蟾基因中高度保守，而在在进化过程中两栖动物最独特的便是变态现象，这提示我们去探索这段爪蟾特有的保守序列是否与两栖类变态发育密切相关。由于甲状腺激素在两栖类的变态中的重要作用，因此我们设想Xfhl3基因的3’UTR中的保守序列可能与甲状腺激素相互作用共同调节爪蟾的变态过程。我们的初步结果表明，在爪蟾胚胎中，甲状腺激素对于正常报告基因表达没有明显的作用，但是在插入Xfhl3基因3’UTR中保守序列后，甲状腺激素处理可显著提高报告基因的表达，表明甲状腺素可能直接或间接通过与该段保守序列参与基因的表达调控。 脊椎动物的眼是一个功能非常特殊的器官，受到复杂的调控网络的调节，众多对神经发生重要的基因在眼中表达并参与了这一调节过程。我们克隆了非洲爪蟾的Sox1基因并研究了它在非洲爪蟾早期发育过程中的时空表达图式，比较了Sox1-3基因在发育的脑和眼中的表达图式，进一步阐明SoxB1基因家族在脊椎动物神经系统发生过程中的作用。此外，我们还克隆了非洲爪蟾MGC85160基因并利用RT-PCR和胚胎整体原位杂交技术探测它在不同胚胎阶段的时空表达图式。结果表明母源性表达的MGC85160基因早期主要在动物极表达；从神经板期开始在发育的中枢神经系统和眼中表达，石蜡切片显示它主要在视网膜和晶状体中表达，说明该基因在爪蟾早期外胚层的模式化以及中枢神经系统的发育过程中可能起到重要作用。 此外，我们还研究了鱇浪白鱼的早期发育分期和眼睛特异基因的表达图式。鱇浪白鱼（Anabarilius grahami ）是云南抚仙湖的特有鱼种。我们首次完成了鱇浪白鱼早期发育的完整分期，主要包括合子期，卵裂期，囊胚期，原肠期，体节期和孵化期六个主要的时期。为了理解鱇浪白鱼眼睛的发育，我们克隆并检测了在眼发育早期起关键作用的基因Sox2, Pax6a, Six3a 和 Rx2的表达图式。结果表明这四个基因全部在尾芽期的前端神经板中表达，随后在视网膜原基细胞中表达明显。在晚期阶段，除Rx2外其它三个基因也在晶状体中表达。其表达模式与斑马鱼中同源基因的表达很相似，说明涉及眼发育的分子网络在鱇浪白鱼中也是高度保守的。

非洲爪蟾Brunol和MCPH基因在早期胚胎发育过程中表达的研究及随机定向酵母双杂交cDNA文库的构建

BRUNOL/CELF家族RNA结合蛋白在转录后调控（post-transcriptional regulation）中起着至关重要的作用，参与多种组织的发育过程。本研究中，我们描述了非洲爪蟾5个Brunol 基因的克隆与表达。其中只有Brunol2是母源性以及合子表达的，其它的4个Brunol基因都是合子表达的，但起始表达的发育时期有所不同。爪蟾发育过程中，Brunol1、4和5基因特异性地在神经系统中表达，包括脑、脊髓、眼泡和耳泡。Brunol2和3基因在体节中胚层与神经系统表达。Brunol2也在晶状体中有非常高的表达。在转染的Hela细胞中，BRUNOL1、2和3蛋白定位于细胞质和细胞核中，BRUNOL4和5只是定位于细胞质中，显示它们具有不同的功能。 人的microcephalin1（MCPH1）基因的遗传突变产生原发性小头症，而在人类的进化过程中，这个基因的变化可能对人脑体积的增加和认知能力的增强也起到重要的作用。但是对于MCPH基因在其它物种中功能的研究才刚刚开始。我们克隆了非洲爪蟾MCPH基因，发现非洲爪蟾具有A，B两个同源基因，其功能域的保守性较高，暗示非洲爪蟾MCPH基因仍然执行一些保守的功能，但MCPHB由于突变只编码一种截短的蛋白，目前尚不清楚它是否是有功能的。胚胎原位杂交的结果显示MCPHA，B基因在胚胎发育中的表达图式相似，但MCPHB的表达水平较低。在神经胚期，二者均表达于头部基板区，在尾芽期主要表达于咽鳃区，而在脑区的表达并不显著，与小鼠中的表达模式不同，提示在爪蟾中MCPH基因可能主要参与咽鳃区而不是脑的发育。 为了进一步筛选这些蛋白可能的结合因子，我们构建了非洲爪蟾双杂交cDNA文库。其中利用修饰的随机引物和特别设计的连接头在合成双链cDNA时，在下游合成一个SalI限制性酶切位点，可以将cDNA定向插入载体。通过对于空克隆率和插入片段长度等一系列参数的分析，表明这个定向cDNA文库的构建是成功的。

Interaction between Thymidylate Synthase and Its Cognate mRNA in Zebrafish Embryos

Thymidylate synthase (TS), which catalyzes the de novo synthesis of dUMP, is an important target for cancer therapy. In this report, the effects of 5-fluorouracil (5-FU) and ZD1694 on the regulation of TS gene expression were evaluated in zebrafish embryos. Our results revealed that the expression of TS was increased by about six-fold when embryos were treated with 1.0 mu M 5-FU and there was a greater than 10-fold increase in the TS protein level after treatment with 0.4 mu M ZD1694. Northern blot analysis confirmed that expression of TS mRNA was identical in treated or untreated embryos. Gel shift and immunoprecipitation assays revealed that zebrafish TS was specifically bound with its cognate mRNA in vitro and in vivo. We identified a 20 nt RNA sequence, TS:N20, localized to the 5'-UTR of TS mRNA, which corresponded to nt 13-32; TS:N20 bound to the TS protein with an affinity similar to that of the full-length TS mRNA. The MFold program predicted that TS:N20 formed a stable stem-loop structure similar to that of the cis-acting element found in human TS mRNA. Variant RNAs with either a deletion or mutation in the core motif of TS:N20 were unable to bind to the TS protein. In vitro translation experiments, using the rabbit lysate system, confirmed that zebrafish TS mRNA translation was significantly repressed when an excess amount of TS protein was included in the system. Additionally, a TS stability experiment confirmed that treatment of zebrafish embryos with 5-FU could increase the TS stability significantly, and the half life of TS protein was about 2.7 times longer than in untreated embryos. Our study revealed a structural requirement for the interaction of TS RNA with TS protein. These findings also demonstrated that the increase in TS protein induced by 5-FU occurs at the post-transcriptional level and that increased stability and translation efficiency both contributed to the increase in TS protein levels induced by TS inhibitors.

Interaction between Thymidylate Synthase and Its Cognate mRNA in Zebrafish Embryos

Thymidylate synthase (TS), which catalyzes the de novo synthesis of dUMP, is an important target for cancer therapy. In this report, the effects of 5-fluorouracil (5-FU) and ZD1694 on the regulation of TS gene expression were evaluated in zebrafish embryos. Our results revealed that the expression of TS was increased by about six-fold when embryos were treated with 1.0 mu M 5-FU and there was a greater than 10-fold increase in the TS protein level after treatment with 0.4 mu M ZD1694. Northern blot analysis confirmed that expression of TS mRNA was identical in treated or untreated embryos. Gel shift and immunoprecipitation assays revealed that zebrafish TS was specifically bound with its cognate mRNA in vitro and in vivo. We identified a 20 nt RNA sequence, TS:N20, localized to the 5'-UTR of TS mRNA, which corresponded to nt 13-32; TS:N20 bound to the TS protein with an affinity similar to that of the full-length TS mRNA. The MFold program predicted that TS:N20 formed a stable stem-loop structure similar to that of the cis-acting element found in human TS mRNA. Variant RNAs with either a deletion or mutation in the core motif of TS:N20 were unable to bind to the TS protein. In vitro translation experiments, using the rabbit lysate system, confirmed that zebrafish TS mRNA translation was significantly repressed when an excess amount of TS protein was included in the system. Additionally, a TS stability experiment confirmed that treatment of zebrafish embryos with 5-FU could increase the TS stability significantly, and the half life of TS protein was about 2.7 times longer than in untreated embryos. Our study revealed a structural requirement for the interaction of TS RNA with TS protein. These findings also demonstrated that the increase in TS protein induced by 5-FU occurs at the post-transcriptional level and that increased stability and translation efficiency both contributed to the increase in TS protein levels induced by TS inhibitors.

Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome

Understanding the dynamics of eukaryotic transcriptome is essential for studying the complexity of transcriptional regulation and its impact on phenotype. However, comprehensive studies of transcriptomes at single base resolution are rare, even for modern organisms, and lacking for rice. Here, we present the first transcriptome atlas for eight organs of cultivated rice. Using high-throughput paired-end RNA-seq, we unambiguously detected transcripts expressing at an extremely low level, as well as a substantial number of novel transcripts, exons, and untranslated regions. An analysis of alternative splicing in the rice transcriptome revealed that alternative cis-splicing occurred in similar to 33% of all rice genes. This is far more than previously reported. In addition, we also identified 234 putative chimeric transcripts that seem to be produced by trans-splicing, indicating that transcript fusion events are more common than expected. In-depth analysis revealed a multitude of fusion transcripts that might be by-products of alternative splicing. Validation and chimeric transcript structural analysis provided evidence that some of these transcripts are likely to be functional in the cell. Taken together, our data provide extensive evidence that transcriptional regulation in rice is vastly more complex than previously believed.

Transcription rates of yeast genes are influenced by the distribution of introns

A comparative analysis on the intron sequence oligonucleotide usages in two sets of yeast genes with higher and lower transcription frequencies, respectively, has shown that the intron sequence structures of the two sets of genes are different. There are more potential binding sites for transcription factors in the introns of the genes with high transcription frequencies. So it is speculated that introns regulate the transcription of genes. But more evidences are needed to favor this speculation. The detailed comparative analyses on the distribution ( length and position) of introns and exons in the two sets of gene sequences also show that there is an obvious boundary between the lengths of the two sets of introns. There is no boundary between the lengths of the two sets of exons, although the means of their lengths are of discrepancy. The situation of the gene lengths ( length of intron and exon) is similar to exon lengths. As far as the relative position, the introns in two sets of genes all have a bias toward the 5' ends of genes. But as the actual position is considered, more introns in high transcription genes have a tendency to be located toward the 5' ends of genes, some even located at 5'-UTR. These results suggest that the gene transcription rates are related to the length of intron, but not to the lengths of exons and genes sequences. The positions of introns may also influence the transcription rates. The transcriptional regulation of introns may be correlative with the transcriptional regulation of the upstream of genes, or be its continuous action.

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

保守序列是一种跨物种保守的基因组序列，而且绝大多数为非蛋白编码序 列。保守序列在人类遗传疾病中发挥着重要作用。其中，一部分保守序列能够 折叠形成二级结构。已鉴定的一些保守二级结构编码一些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 仅仅 只结合到很小一部分保守二级结构上，提示可能还有更多的转录因子会结合到 保守二级结构上。因此，保守二级结构介导的转录调控网络要比目前已知的复 杂得多。

A key gene of the RNA interference pathway in the black tiger shrimp, Penaeus monodon: Identification and functional characterisation of Dicer-1

RNA interference (RNAi) is an evolutionarily conserved mechanism by which double-stranded RNA (dsRNA) initiates post-transcriptional silencing of homologous genes. Here we report the amplification and characterisation of a full length cDNA from black tiger shrimp (Penaeus monodon) that encodes the bidentate RNAase III Dicer, a key component of the RNAi pathway. The full length of the shrimp Dicer (Pm Dcr1) cDNA is 7629 bp in length, including a 51 untranslated region (UTR) of 130 bp, a 3' UTR of 77 bp, and an open reading frame of 7422 bp encoding a polypeptide of 2473 amino acids with an estimated molecular mass of 277.895 kDa and a predicted isoelectric point of 4.86. Analysis of the deduced amino acid sequence indicated that the mature peptide contains all the seven recognised functional domains and is most similar to the mosquito (Aedes aegypti) Dicer-1 sequence with a similarity of 34.6%. Quantitative RT-PCR analysis showed that Pm Dcr1 mRNA is most highly expressed in haemolymph and lymphoid organ tissues (P 0.05). However, there was no correlation between Pm Dcr1 mRNA levels in lymphoid organ and the viral genetic loads in shrimp naturally infected with gill-associated virus (GAV) and Mourilyan virus (P > 0.05). Treatment with synthetic dsRNA corresponding to Pm Dcr1 sequence resulted in knock-down of Pm Dcr1 mRNA expression in both uninfected shrimp and shrimp infected experimentally with GAV. Knock-down of Pm Dcr1 expression resulted in more rapid mortalities and higher viral loads. These data demonstrated that Dicer is involved in antiviral defence in shrimp. (c) 2007 Elsevier Ltd. All rights reserved.

Molecular cloning, characterization and evolutionary analysis of phytoene desaturase (PDS) gene from Haematococcus pluvialis

Phytoene desaturase is one of the most important enzymes necessary for the biosynthesis of carotenoids in some cyanobacteria, green algae and plants. In this study, genomic DNA and cDNA of pds were cloned from unicellular green alga Haematococcus pluvialis strain323 using PCR and RT-PCR methods. The cDNA was cloned into plasmid pET-28a and efficiently expressed in Escherichia coli BL21. The complete genomic PDS gene of H. pluvialis, 3.3 kb in size, included eight exons and seven introns. To locate transcriptional regulation elements, an approximate 1 kb of 5'-flanking region was isolated by genome-walking method. Results of bioinformatic analysis showed several putative cis-elements e.g. the ABRE motif (abscisic acid responsive element), the C-repeat/DRE (dehydration responsive element) motif and the GCN4 motif were located in 5'-flanking region of pds. Results of phylogenetic analyses reveal that different sources of PDS genes form a separate clade, respectively, with 100% bootstrap support. Moreover, a maximum likelihood approach was employed to detect evidence of positive selection in the evolution of PDS genes. Results of branch-site model analysis suggest that 7.9% of sites along the green algal branch are under positive selection, and the PDS gene in green algae exhibits a different evolutionary pattern from its counterparts in cyanobacteria and plants.

Reversed-phase high-performance liquid chromatographic investigation of urinary normal and modified nucleosides of cancer patients

Post-transcriptional modifications in RNA give rise to free modified ribonucleosides circulating in the blood stream and excreted in urine. Due to their abnormal levels in conjunction with several tumor diseases, they have been suggested as possible tumor markers. The developed RP-HPLC method has been applied to analyze the urinary nucleosides in 34 urinary samples from 15 kinds of cancer patients. The statistical analyses showed the urinary nucleoside excretion, especially modified nucleoside levels, in cancer patients were significantly higher than those in normal healthy volunteers. Factor analysis was used to classify the patients with cancer and normal healthy humans. It was found that using 15 urinary nucleoside levels or only five modified nucleoside levels as data vectors the factor analysis plot displayed two almost separate clusters representing each group. (C) 1999 Elsevier Science B.V. All rights reserved.