Imputation of missing covariate values in epigenome-wide analysis of DNA methylation data

DNA methylation is a widely studied epigenetic mechanism and alterations in methylation patterns may be involved in the development of common diseases. Unlike inherited changes in genetic sequence, variation in site-specific methylation varies by tissue, developmental stage, and disease status, and may be impacted by aging and exposure to environmental factors, such as diet or smoking. These non-genetic factors are typically included in epigenome-wide association studies (EWAS) because they may be confounding factors to the association between methylation and disease. However, missing values in these variables can lead to reduced sample size and decrease the statistical power of EWAS. We propose a site selection and multiple imputation (MI) method to impute missing covariate values and to perform association tests in EWAS. Then, we compare this method to an alternative projection-based method. Through simulations, we show that the MI-based method is slightly conservative, but provides consistent estimates for effect size. We also illustrate these methods with data from the Atherosclerosis Risk in Communities (ARIC) study to carry out an EWAS between methylation levels and smoking status, in which missing cell type compositions and white blood cell counts are imputed.     

Cytochromes P450 in phenolic metabolism

Three independent cytochrome P450 enzyme families catalyze the three rate-limiting hydroxylation steps in the phenylpropanoid pathway leading to the biosynthesis of lignin and numerous other phenolic compounds in plants. Their characterization at the molecular and enzymatic level has revealed an unexpected complexity of phenolic metabolism as the major route involves shikimate/quinate esters and alcohol/aldehyde intermediates. Engineering expression of CYP73s (encoding cinnamate 4-hydroxylase), CYP98s (encoding 4-coumaroylshikimate 3′-hydroxylase) or CYP84s (encoding coniferaldehyde 5-hydroxylase) leads to modified lignin and seed phenolic composition. In particular CYP73s and CYP98s also play essential roles in plant growth and development, while CYP84 constitutes a check-point for the synthesis of syringyl lignin and sinapate esters. Although recent data shed new light on the main path for lignin synthesis, they also raised new questions. Mutants and engineered plants revealed the existence of (an) alternative pathway(s), which most likely involve(s) different precursors and oxygenases. On the other hand, phylogenetic analysis of plant genomes show the existence of P450 gene duplications in each family, which may have led to the acquisition of novel or additional physiological functions in planta. In addition to the main lignin pathway, P450s contribute to the biosynthesis of many bioactive phenolic derivatives, with potential applications in medicine and plant defense, including lignans, phenylethanoids, benzoic acids, xanthones or quinoid compounds. A very small proportion of these P450s have been characterized so far, and rarely at a molecular level. The possible involvement of P450s in salicylic acid is discussed.     

Phanerochaete chrysosporium NADPH-cytochrome P450 reductase kinetic mechanism

The recently completed genome of the basidiomycete, Phanerochaete chrysosporium, revealed the presence of one NADPH-cytochrome P450 oxidoreductase (CPR; EC 1.6.2.4) gene and >123 cytochrome P450 (CYP) genes. How a single CPR can drive many CYPs is an important area of study. We have investigated this CPR to gain insight into the mechanistic and structural biodiversity of the cytochrome P450 catalytic system. Native CPR and a NH(2)-terminally truncated derivative lacking 23 amino acids have been overexpressed in Escherichia coli and purified to electrophoretic homogeneity. Steady-state kinetics of cytochrome c reductase activity revealed a random sequential bireactant kinetic mechanism in which both products form dead-end complexes reflecting differences in CPR kinetic mechanisms even within a single kingdom of life. Removal of the N-terminal anchor of P. chrysosporium CPR did not alter the kinetic properties displayed by the enzyme in vitro, indicating it was a useful modification for structural studies.     

定向进化细胞色素P450BM-3的研究

以来自于巨大芽孢杆菌的细胞色素P450BM-3为研究对象,采用随机突变和饱和定点突变定向进化技术对P450BM-3进行改造,通过突变体催化靛蓝显色的特性采用活性琼脂平板分析和96微孔板相结合的高通量筛选成功获得了几个具有更高催化性能的突变体。     

CYP79B1 from Sinapis alba converts tryptophan to indole-3-acetaldoxime

The cytochrome P450 CYP79B1 from Sinapis alba has been heterologously expressed in Escherichia coli and shown to catalyze the conversion of tryptophan to indole-3-acetaldoxime. Three expression constructs were made, one expressing the native protein and two expressing proteins with different N-terminal modifications. The native construct gave the highest yield as estimated by enzymatic activity per liter of culture. Spheroplasts of E. coli expressing CYP79B1 were reconstituted with the Arabidopsis thaliana NADPH:cytochrome P450 reductase ATR1 heterologously expressed in E. coli to obtain enzymatic activity. This indicates that the E. coli electron-donating system, flavodoxin/flavodoxin reductase, does not support CYP79B1 activity. Recombinant CYP79B1 has a K(m) for tryptophan of 29+/-2 microM and a V(max) of 36.5+/-0.7nmolh(-1)(mlculture)(-1). The identity at the amino acid level of CYP79B1 is, respectively, 93 and 84% to CYP79B2 and CYP79B3 from A. thaliana, and 96% to CYP79B5 (Accession No. AF453287) from Brassica napus. The CYP79B subfamily of cytochromes P450 is likely to constitute a group of orthologous genes in the biosynthesis of indole glucosinolates.     

热带假丝酵母多倍体变种的细胞色素P450和尿素在烷烃代谢中的作用

建立了可行的细胞色素 P450测定方法。考察了以烷烃为单一碳源的酵母细胞色素P450的一氧化碳差示光谱,峰值约为455nm。观察了烷烃培养的酵母细胞色素 P450在生长期中的消长。比较了十四醇、十四醇添加苯巴比妥、以及正十四烷等三种不同培养条件下酵母细胞色素P450的含量和发酵产物的成份。结果表明,细胞色素 P450为烷烃转化成二元酸所必需。烷烃为单一碳源培养酵母时,培养基中过量尿素(0.2%以上)促进烷烃利用和酵母生长,降低细胞色素 P450生成和二元酸的积累。根据上述实验结果和本研究室以前报道,提出了烷烃代谢调控模式。     

热带假丝酵母1230细胞色素P450基因CYPA14与CYPA16的克隆和表达

细胞色素P450(CYP)是一种单加氧酶,在热带假丝酵母(Candidatropicalis)ω-氧化过程中发挥关键作用。通过对来源不同的P450基因进行同源性分析,首先克隆到热带假丝酵母1230中P450基因的部分序列,再利用基因组步行法克隆其未知序列,结果分别获得了两个P450同工酶基因CYPA14和CYPA16的完整序列。经PCR方法证实,二者在染色体上的位置相邻,其读码框分别编码522和540个氨基酸残基的肽链。经NCBIBLAST搜索比较后发现,二者与热带假丝酵母ATCC20336中的P450成员CYP52A14和CYP52A16分别编码的序列几乎完全一致,与热带假丝酵母ATCC750中的P450成员CYP52A2和CYP52A1也具有较高的相似性。同时,对经诱变后的几株二元酸生产菌株的CYPA14与CYPA16也进行了克隆和序列比较,发现部分序列中的个别氨基酸残基发生了突变。CYPA14和CYPA16均在酿酒酵母中获得了有效表达,其中CYPA16的P450表达含量高于CYPA14,后者有部分表达产物发生了变性。     

家蚕细胞色素P450基因CYP6AE21的克隆、表达分析及亚细胞定位

信号失活是嗅觉动态过程的一个重要步骤, 这一过程涉及多样的气味降解酶类。本研究利用RT-PCR方法从家蚕Bombyx mori雄蛾的触角中克隆了一个细胞色素P450基因CYP6AE21。该基因含有一个1 572 bp的开放阅读框（open reading frame, ORF）, 编码523个氨基酸, 预测分子量为60.5 kD, 等电点为8.4, 含有细胞色素P450的特征序列血红素结合位点区域。CYP6AE21和CYP6AE2基因一样在相同位置含有1个内含子序列, 且相应的2个外显子大小相同。两者的核苷酸序列相似性达到94.5%, 且在基因组上以头尾相连的方式成簇排列, 中间由约7.6 kb核苷酸序列隔开。CYP6AE21在幼虫的头部和脂肪体, 以及雄蛾和雌蛾的触角中表达量较高; 在幼虫的其他组织和蛾的多个组织中也有一定的表达。P450酶系的重要组分之一--NADPH细胞色素P450还原酶（cytochrome P450 reductase, CPR）基因也在雌蛾和雄蛾触角中高水平表达, 而在其他组织中表达量相对较低。亚细胞定位分析表明CYP6AE21表达产物定位于细胞质中。推测CYP6AE21和CYP6AE2是由其中一个基因加倍复制形成的; 此P450的功能之一可能是参与内化进细胞的气味分子的降解清除。     

Induction of Two Forms of Eel Cytochrome P450 1A Genes by 3-Methylcholanthrene

In eel (Anguilla japonica), exposure to polyaromatic hydrocarbons such as 3-methylcholanthrene leads to induction of two CYP1A enzymes, CYP1A1 and CYP1A6. We studied the time course and tissue specificity of induction of messenger RNAs for CYP1A1 and CYP1A6 in eel by administering 3-methylcholanthrene intraperitoneally. In both cases, the drug induced a rapid increase of mRNAs and biphasic expression. In the liver, mRNA levels of CYP1A1 and CYP1A6 increased 22-fold at 3 hours and 27-fold at 6 hours after the administration, respectively, showing initial peaks in the induction. After the initial inductions, mRNA levels decreased unexpectedly. Following these temporary decreases, the mRNA levels again increased and reached levels that were 35 and 41 times the basal levels at 24 hours after administration, respectively. CYP1A1 and CYP1A6 resembled each other also in the tissue specificity of gene expression; the expression levels were liver ≫ gill > intestine > kidney. The rapid induction, the biphasic expression, and the tissue-specific expression were common features of gene expression in CYP1A1 and CYP1A6 and may come from common structures of the regulatory regions of the two genes. Received December 7, 1998; accepted February 15, 1999