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

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

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

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

人与小鼠核糖体蛋白基因内含子中的转录调控位点分析

前期对酵母和果蝇核糖体蛋白(Ribosomal protein, RP)基因内含子序列中的寡核苷酸分析表明, 内含子中含有潜在的转录因子结合位点。为进一步发掘核糖体蛋白基因内含子参与转录调控的证据, 文章首先基于频率分析方法抽提出人和小鼠核糖体蛋白基因第一内含子中高频(Over-represented)出现的寡核苷酸片段 (亦称模体, Motif), 这些寡核苷酸中超过85%与已知的转录因子结合位点吻合, 是潜在的转录调控元件。对抽提出的寡核苷酸进行碱基组成分析, 发现95%以上的寡核苷酸富含碱基C和G, 而较少富含A和T。从寡核苷酸在内含子中的分布情况看, 它们相对靠近第一内含子的5′端, 即距离基因转录起始位点和上游区域较近。推测这些特征可能与基因转录调控有关。     

PKC对人mdr1基因转录调控的研究

研究了ＰＫＣα、β１和β２亚型对多药耐药基因ｍｄｒ１转录的调控作用．以ｍｄｒ１基因上游调控序列驱动的虫荧光素酶报告基因表达载体和ＰＫＣ基因表达载体共同转染ＣＯＳ７细胞后,测定虫荧光素酶报告基因表达水平．实验结果表明,ＰＫＣα、β１及β２对ｍｄｒ１基因的转录有下调作用,ＰＭＡ可进一步加强这一作用;在此转染体系中,ＰＫＣ抑制剂ｓｔａｕｒｏｓｐｏｒｉｎｅ也可抑制ｍｄｒ１基因的转录,但在只转染ｍｄｒｌｕｃ的细胞中,ｓｔａｕｒｏｓｐｏｒｉｎｅ对ｍｄｒ１的转录则没有影响,提示ｓｔａｕｒｏｓｐｏｒｉｎｅ仅在ＰＫＣ过量表达的细胞中抑制ｍｄｒ１基因的转录．     

酵母基因内含子中二聚体寡核苷酸转录调控位点的统计分析

对高效和低效转录酵母基因内含子序列中寡核苷酸的出现频率进行对照分析, 结果显示高效和低效内含子序列的结构有差异, 而且高效转录内含子序列含有较多潜在的转录因子结合位点. 观察实验获得的转录调控位点, 发现许多调控位点不是相邻接的寡核苷酸,而是由一对保守寡核苷酸构成, 这对寡核苷酸被一段长度固定的非保守区域间隔开. 于是对此形式的二聚体寡核苷酸(dyad)在高效和低效内含子序列中出现的频率进行统计比较分析,抽提出在高效内含子组出现的频率显著高于在低效内含子组出现频率的二聚体寡核苷酸, 分析这些二聚体寡核苷酸在两组内含子序列中的分布特征, 并对照实验结果, 这些二聚体寡核苷酸可能与基因转录的正调控有关.     

10个线粒体基因在甘蓝型油菜NCa不育系统不同组织中的转录调控及其与育性的关系

用10个线粒体基因为探针,对NCα不育系、保持系和可育F1的苗期叶片、幼蕾及未成熟种子的线粒体RNA进行了Northern分析。结果表明,这10个线粒体基因除atp6外,其余9个基因在同一材料的不同组织中没有表达差异,都属于组成型表达的线粒体基因。其中,off139、orf222、atp1、cox1、cox2、cob、rm5S、rm26S等8个线粒体基因在不育系、保持系和可育F1的苗期叶片、幼蕾及未成熟种子中有着相同的表达,属于表达不受核基因型影响,没有组织特异性的类型：atp9基因分别在同一材料的不同组织中的转录也基本没有变化,但是在3个不同的材料间具有表达差异：可能属于表达受核基因型影响、没有组织特异性的线粒体基因。atp6基因也在3个材料的叶、蕾和种子中都产生相同大小的转录本,但是在各个材料的不同组织中存在着信号强度的差异,可能是属于表达既受核基因型影响、又有组织特异性的线粒体基因。Orf222和off139分别在不育系和可育F1幼蕾中产生相同大小和丰度的转录本,但是在保持系幼蕾中没有检测到转录本;orf222检测到的3条转录本分别为1．1kb、0．9kb、0．6kb,而off139检测到0．8kb和0．6kb两条带。atp9探针在不育系和保持系幼蕾中都检测到1条0．6kb的转录本,而在可育F1幼蕾中检测到0．6kb和1．2kb的转录本。讨论了orf222、off139、atp9基因的表达与NCα细胞质雄性不育的可能关系。     

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

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

肺癌细胞中调控ezrin基因基本转录活性的顺式作用元件及转录因子的鉴定

研究发现,质膜-细胞骨架连接蛋白Ezrin在多种肿瘤细胞中异常表达,而且Ezrin的表达上调与肿瘤细胞的移动侵袭相关,但是调控ezrin基因转录的分子机制却不清楚.为了探明ezrin基因的转录调控机制,以肺癌细胞A549为材料,首先采用双荧光素酶报告基因分析系统检测ezrin基因5′侧翼嵌套缺失序列和位点突变序列的转录活性,鉴定肺癌细胞中ezrin基因的基本启动子区以及关键的顺式作用元件Sp1结合位点(-75/-69)和AP-1结合位点(-64/-58).其次,利用凝胶电泳迁移率变动分析证明,肺癌细胞核蛋白提取物能够与ezrin基因含有关键顺式作用元件的DNA序列结合,形成DNA-核蛋白复合物,而且Sp1结合位点和AP-1结合位点与重组蛋白rhSp1和rhAP-1的结合具有位点特异性.最后,利用瞬时转染实验证实,转录因子Sp1和AP-1(由c-Jun和c-Fos组成的异源二聚体)分别通过Sp1结合位点和AP-1结合位点,增强ezrin基因基本转录活性,而且,过表达转录因子Sp1、c-Jun或c-Fos上调了Ezrin蛋白表达.研究确定,肺癌细胞中调控ezrin基因基本转录活性的关键顺式作用元件是Sp1结合位点(-75/-69)和AP-1结合位点(-64/-58),与之作用的转录因子Sp1和AP-1对于ezrin基因的转录激活作用至关重要.     

人SATB1基 5'上游序列的转录激活功能分析

在对人SATB1基因进行生物信息学分析的基础上,采用PCR技术,扩增人基因组DNA中SATB1基因5'上游序列的-2955～-9片段,构建了3个分别由SATB1基因5'上游-2955～-9,-1727～-9和-760～-9序列片段驱动的报告载体-pGL3-SP2946-luc,pGL3-SP1718-luc和pGL3-SP751-luc,分别瞬时转染Jurkat T,K562,U937和HeLa细胞,通过测定荧光素酶的表达活性,观察SATB1基因5'上游序列片段3个删除突变体在不同细胞内活性的差异.结果显示,SATB1上游序列-2955/-9在4种细胞中的转录激活能力为U937>Jurkat T>K562,在HeLa细胞中基本无激活,提示SATB1的转录激活可能具有一定的细胞类型特异性.3种5'删除突变体转录激活性由大至小顺序为-760/-9>-2955/-9>-1727/-9,提示SATB1的核心启动子可能存在于其5'上游序列的-760至-9 bp区域中.     

Finite element modeling of the human kidney for probabilistic occupant models: Statistical shape analysis and mesh morphing

Statistical shape analysis was conducted on 15 pairs (left and right) of human kidneys. It was shown that the left and right kidney were significantly different in size and shape. In addition, several common modes of kidney variation were identified using statistical shape analysis. Semi-automatic mesh morphing techniques have been developed to efficiently create subject specific meshes from a template mesh with a similar geometry. Subject specific meshes as well as probabilistic kidney meshes were created from a template mesh. Mesh quality remained about the same as the template mesh while only taking a fraction of the time to create the mesh from scratch or morph with manually identified landmarks. This technique can help enhance the quality of information gathered from experimental testing with subject specific meshes as well as help to more efficiently predict injury by creating models with the mean shape as well as models at the extremes for each principal component.     

Global analysis of metabolites in rat and human urine based on gas chromatography/time-of-flight mass spectrometry

Sediment in urine may contain low-molecular-weight compounds that should be included in the analysis. To date, no systematic investigation has addressed this issue. We investigated three primary factors that influence the extraction efficiency of metabolites during preparation of urine samples for metabolomic research: centrifugation, pH, and extraction solvents. Obtained with the use of gas chromatography/time-of-flight mass spectrometry (GC/TOFMS) technique and principal component analysis (PCA), our results indicate that (1) conventional centrifugation causes an apparent loss of some metabolites, indicating that urine samples for metabolomic research should not be centrifuged before procedures are undertaken to recover the metabolites; (2) pH adjustment has a large impact on the recovery of metabolites and is therefore not encouraged; (3) with design of experiment analysis, methanol and water yield the optimal extraction efficiency. Differences between rat and human urine were observed and are discussed. Ninety-nine metabolites identified in rat and human urine are presented. An efficient protocol is proposed for the pretreatment of urine samples.     

叶用芥菜种质表型性状的遗传多样性分析

为有效利用叶用芥菜种质资源,对筛选出的24份叶用芥菜品系资源的11个表型性状,应用相关系数、聚类分析和主成分分析等多元统计方法进行遗传多样性分析。结果表明:叶用芥菜种质资源具有丰富的遗传多样性;11个表型性状平均变异系数为44.74%,叶片数的变异系数最大,为150.80%;净菜率变异系数最小,仅为9.54%。通过系统聚类,将参试的24份品系分为3个类群,第Ⅰ类群有14份材料,第Ⅱ类群有5份材料,第Ⅲ类群有5份材料,各类群性状之间的差异较明显,明确了品系类群间存在的亲缘关系。在主成分分析中,可选取方差累计贡献率为86.62%的前5个主成分来评价24份叶用芥菜品系资源。本研究揭示了叶用芥菜不同品系的表型特异性和遗传多样性,筛选出一些特异品系资源,为高产优质叶用芥菜新品种(系)的选育提供重要的科学依据。     

基于高频系数小波分析的高维蛋白质波谱数据特征提取

高维蛋白质波谱数据分析过程中,对于数据的特征提取一直是许多学者专注解决的问题。本文提出了一种基于高频系数的小波分析和主成份分析技术(Principal component analysis,PCA)的特征提取方法,首先采用小波分析技术对数据进行降噪,提取高频系数作为特征,之后用主成份分析技术进行降维。实验显示:本论文中提出的方法在8-7-02、4/3/02数据集上的实验识别率分别可以达到100%和99.45%,可以有效提高分类识别率。     

柠檬明串珠菌及相近种部分持家基因的系统发育分析

[目的]利用16S rRNA、dnaA、murC和pyrG基因分子标记研究Leuconostoc citreum(Leu.citreum)及相近种间的种系发育关系,并比较这些基因序列对Leu.citreum及相近种的区分能力.[方法]以分离自酸面团中的7株Leu.citreum为研究对象,以dnaA、murC和pyrG基因片段为标记,通过PCR扩增、测序,结合已公布的近缘种及亚种相应序列,计算遗传距离,构建系统发育树,并与16S rRNA基因进行比较.[结果]研究发现Leu.citreum及相近种间的dnaA、murC和pyrG基因构建的系统发育树拓扑结构与16S rRNA基因基本一致,区别在于相似性的不同,其分别为75.5％-97.2％、50.2％-99.7％、65.0％-99.8％和98.5％-100％.[结论]在Leu.citreum及相近种间的种系发育关系中,dnaA、murC和pyrG基因与16S rRNA基因系统进化关系都具有很好的一致性,但这三个持家基因的遗传距离显著高于16S rRNA基因.因此,采用dnaA、murC和pyrG基因可以用于Leu.citreum及相近种的分类鉴定.     

Independent component analysis of Raman spectra: Application on paraffin-embedded skin biopsies

Raman spectra provide wealthy but complex information about the chemical constituents of biological samples. Digital processing techniques are usually needed to extract the spectra of chemical constituents and their associated concentration profiles. However, spectral signatures may admit transformations from those recorded on pure constituents and these techniques require a priori knowledge of spectra to be estimated. We propose in this study to analyse paraffin-embedded skin biopsies of malignant and benign tumors dedicated to oncology researches by Raman spectroscopy and advanced signal processing methods. We show that the commonly used principal component analysis (PCA) does not give physically interpretable estimators of spectra and associated concentration profiles. Based on a linear model and taking into account the statistical properties of spectra, independent component analysis (ICA) is used to better estimate the spectra of chemical constituents. The estimators of associated concentration profiles are no longer orthogonal and have only positive values, contrary to PCA. ICA allows to model the paraffin by three Raman spectra and provides good estimators of underlying spectra of the human skin, which is of great interest in oncology since the retrieval of spectral features of different types of skin tumors is sufficient for their discrimination.     

Application of principal component analysis in clinical gait research: Identification of systematic differences between healthy and medial knee-osteoarthritic gait

For a successful completion of a movement task the motor control system has to observe a multitude of internal constraints that govern the coordination of its segments. The purpose of this study was to apply principal component (PC) analysis to detect differences in the segmental coordination between healthy subjects and patients with medial knee osteoarthritis (OA). It was hypothesized that (1) systematic differences in patterns of whole body movement would be identifiable with this method even in small sample sized groups and that (2) these differences will include compensatory movements in the OA patients in both the lower and upper body segments. Marker positions and ground reaction forces of three gait trials of 5 healthy and 5 OA participants with full body marker sets were analyzed using a principal component analysis. Group differences in the PC-scores were determined for the first 10 PC-vectors and a linear combination of those PC-vectors where differences were found defined a discriminant vector. Projecting the original trials onto this discriminant vector yielded significant group differences (t(d=8)=3.011; p=0.017) with greater upper body movement in patients with knee OA that was correlated with the medial–lateral ground reaction force. These results help to characterize the adaptation of whole-body gait patterns to knee OA in a relatively small population and may provide an improved basis for the development of interventions to modify knee load. The PC-based motion analysis offered a highly sensitive approach to identify characteristic whole body patterns of movement associated with pathological gait.     

生信分析筛选心衰氧化应激相关关键[]基因并预测候选药物

为探究氧化应激相关基因在心力衰竭发生发展中的作用,并发现核心基因进行靶基因药物预测。从GEO数据库下载GSE120895基因表达图谱,通过GEO2R筛选差异表达基因,将差异表达基因与GeneCard数据库中筛选的氧化应激相关基因取交集,得到心力衰竭氧化应激相关差异表达基因,利用R软件对差异表达基因进行GO及KEGG分析,利用Cytoscape进行PPI网络的模块以及关键基因的筛选。之后在GSE17800基因表达图谱中验证关键基因的表达,并针对关键基因进行相互作用药物预测。差异表达基因与氧化应激相关基因取交集后,共筛选出52个上调的氧化应激相关差异表达基因,在此基础上,筛选出ACTB,STAT3,FN1,EDN1,CAT共5个关键基因,在GSE17800基因表达图谱中验证后,针对4个关键基因预测了19个靶基因潜在药物。总之,本研究通过生物信息学方法鉴定关键基因,并预测潜在治疗药物,从而为了解心力衰竭的分子机制及其诊治方法提供新的见解。     

粪肠、屎肠球菌及相近种部分持家基因的系统发育分析

【目的】利用16S rRNA、clpX和recA基因分子标记研究Enterococcus faecalis、Enterococcus faecium及相近种间的种系发育关系,并比较这些基因序列对E.faecalis、E.faecium及相近种的区分能力。【方法】以分离自传统乳制品中的9株E.faecium和1株E.durans分离株为研究对象,以clpX和recA基因片段为标记,通过PCR扩增、测序,结合已公布的近缘种相应序列构建系统发育树并与16S rRNA基因进行比较。【结果】在基于clpX和recA基因的进化树中,10株试验菌株与E.faecalis始终处于同一分支。与该物种这两个基因的平均相似性为99.6%和98.6%,与另一分支的Faecium-group(E.durans和E.faecium)的平均相似性仅为61.5%和33.5%。相近种E.durans和E.hirae间这两个基因的差异性为20.3%和39.0%;在基于16S rRNA基因的进化树中,试验菌株与Faecium-group(E.lactis、E.faecium、E.durans、E.hirae)处于同一分支。与这些成员间该基因的相似性大于99.6%,与E.faecalis基因的平均相似性可达98.4%。相近种间该基因相似性无明显差异。【结论】按照10株试验菌株clpX和recA基因的分析结果可将由传统生理生化和16S rRNA基因序列鉴定的9株E.faecium和1株E.durans归类为E.faecalis,clpX和recA基因可用于部分相近种的分类鉴定。     

Genetic diversity and phylogeny analysis of Anabaena azollae based on RFLPs detected in Azolla-Anabaena azollae DNA complexes using nif gene probes

The cyanobacterium Anabaena has both symbiotic and free-living forms. The genetic diversity of Anabaena strains symbiotically associated with the aquatic fern Azolla and the evolutionary relationships among these symbionts were evaluated by means of RFLP (restriction fragment length polymorphism) experiments. Three DNA fragments corresponding to nif genes were cloned from the free-living cyanobacterium Anabaena PCC 7120 and used as probes. A mixture of Azolla, Anabaena and bacterial DNA was extracted from Azolla fronds and digested with two restriction enzymes. Single-copy RFLP signals were detected with two of the probes in all Azolla Anabaena examined. Multiple-copy RFLP signals were obtained from the third probe which corresponded to a part of the nif N gene. A total of 46 probe/enzyme combinations were scored as present or absent and used to calculate pairwise Nei's genetic distances among symbiotic Anaebaena strains. Phylogenetic trees summarizing phenetic and cladistic relationships among strains were generated according to three different evolutionary scenarios: parsimony, UPGMA and neighbour joining. All trees revealed identical phylogenetic relationships. Principal component analysis was also used to evaluate genetic similarities and revealed three groups: group one contains the cyanobacteria associated with plants from the Azolla section, group two contains those associated with plants from the pinnata species and group three contains those associated with plants from the nilotica species. The same groups had already been identified earlier in a random amplified polymorphic DNA (RAPD) analysis of Azolla-Anbaena DNA complexes, suggesting that the present Azolla taxonomy should be revised. We now suggest a taxonomy of Anabaena azollae that is parallel to such a revised Azolla taxonomy. An Azolla chloroplast DNA sequence derived from Oryza sativa was also used as an RFLP probe on Azolla DNA to confirm the presence of plant DNA in the total genomic DNA extracted from ferns with or without the symbiont. Our results also suggest that total DNA extracted from the Azolla-Anabaena complexes includes both plant and symbiont DNA and can be used equally well for RFLP analysis of host plant or symbiotic cyanobacteria.