Powerful identification of cis-regulatory SNPs in human primary monocytes using allele-specific gene expression.

A large number of genome-wide association studies have been performed during the past five years to identify associations between SNPs and human complex diseases and traits. The assignment of a functional role for the identified disease-associated SNP is not straight-forward. Genome-wide expression quantitative trait locus (eQTL) analysis is frequently used as the initial step to define a function while allele-specific gene expression (ASE) analysis has not yet gained a wide-spread use in disease mapping studies. We compared the power to identify cis-acting regulatory SNPs (cis-rSNPs) by genome-wide allele-specific gene expression (ASE) analysis with that of traditional expression quantitative trait locus (eQTL) mapping. Our study included 395 healthy blood donors for whom global gene expression profiles in circulating monocytes were determined by Illumina BeadArrays. ASE was assessed in a subset of these monocytes from 188 donors by quantitative genotyping of mRNA using a genome-wide panel of SNP markers. The performance of the two methods for detecting cis-rSNPs was evaluated by comparing associations between SNP genotypes and gene expression levels in sample sets of varying size. We found that up to 8-fold more samples are required for eQTL mapping to reach the same statistical power as that obtained by ASE analysis for the same rSNPs. The performance of ASE is insensitive to SNPs with low minor allele frequencies and detects a larger number of significantly associated rSNPs using the same sample size as eQTL mapping. An unequivocal conclusion from our comparison is that ASE analysis is more sensitive for detecting cis-rSNPs than standard eQTL mapping. Our study shows the potential of ASE mapping in tissue samples and primary cells which are difficult to obtain in large numbers.

Identification of epigenetically regulated genes that predict patient outcome in neuroblastoma.

BACKGROUND: Epigenetic mechanisms such as DNA methylation and histone modifications are important regulators of gene expression and are frequently involved in silencing tumor suppressor genes. METHODS: In order to identify genes that are epigenetically regulated in neuroblastoma tumors, we treated four neuroblastoma cell lines with the demethylating agent 5-Aza-2'-deoxycytidine (5-Aza-dC) either separately or in conjunction with the histone deacetylase inhibitor trichostatin A (TSA). Expression was analyzed using whole-genome expression arrays to identify genes activated by the treatment. These data were then combined with data from genome-wide DNA methylation arrays to identify candidate genes silenced in neuroblastoma due to DNA methylation. RESULTS: We present eight genes (KRT19, PRKCDBP, SCNN1A, POU2F2, TGFBI, COL1A2, DHRS3 and DUSP23) that are methylated in neuroblastoma, most of them not previously reported as such, some of which also distinguish between biological subsets of neuroblastoma tumors. Differential methylation was observed for the genes SCNN1A (p < 0.001), PRKCDBP (p < 0.001) and KRT19 (p < 0.01). Among these, the mRNA expression of KRT19 and PRKCDBP was significantly lower in patients that have died from the disease compared with patients with no evidence of disease (fold change -8.3, p = 0.01 for KRT19 and fold change -2.4, p = 0.04 for PRKCDBP). CONCLUSIONS: In our study, a low methylation frequency of SCNN1A, PRKCDBP and KRT19 is significantly associated with favorable outcome in neuroblastoma. It is likely that analysis of specific DNA methylation will be one of several methods in future patient therapy stratification protocols for treatment of childhood neuroblastomas.

Differential genome-wide array-based methylation profiles in prognostic subsets of chronic lymphocytic leukemia.

Global hypomethylation and regional hypermethylation are well-known epigenetic features of cancer; however, in chronic lymphocytic leukemia (CLL), studies on genome-wide epigenetic modifications are limited. Here, we analyzed the global methylation profiles in CLL, by applying high-resolution methylation microarrays (27,578 CpG sites) to 23 CLL samples, belonging to the immunoglobulin heavy-chain variable (IGHV) mutated (favorable) and IGHV unmutated/IGHV3-21 (poor-prognostic) subsets. Overall, results demonstrated significant differences in methylation patterns between these subgroups. Specifically, in IGHV unmutated CLL, we identified methylation of 7 known or candidate tumor suppressor genes (eg, VHL, ABI3, and IGSF4) as well as 8 unmethylated genes involved in cell proliferation and tumor progression (eg, ADORA3 and PRF1 enhancing the nuclear factor-kappaB and mitogen-activated protein kinase pathways, respectively). In contrast, these latter genes were silenced by methylation in IGHV mutated patients. The array data were validated for selected genes using methylation-specific polymerase chain reaction, quantitative reverse transcriptase-polymerase chain reaction, and bisulfite sequencing. Finally, the significance of DNA methylation in regulating gene promoters was shown by reinducing 4 methylated tumor suppressor genes (eg, VHL and ABI3) in IGHV unmutated samples using the methyl-inhibitor 5-aza-2'-deoxycytidine. Taken together, our data for the first time reveal differences in global methylation profiles between prognostic subsets of CLL, which may unfold epigenetic silencing mechanisms involved in CLL pathogenesis.

A genotyping method for VKORC1 1173C > T by Pyrosequencing technology.

Vitamin K epoxide reductase complex subunit 1 (VKORC1) is the site of inhibition by warfarin and other antivitamin K drugs during oral anticoagulant therapy. The SNP rs9934438 in intron 1 of VKORC1 (c.173+1000C > T or 1173C > T) discriminating the VKORC1*2 haplotype is associated with low warfarin dose requirement and unstable prothrombin time - international normalized ratio. To genotype this SNP, we have developed a rapid method using Pyrosequencing technology. The proposed method takes a post-PCR sample preparation of less than 1 h and a DNA sequencing time of less than 15 min to genotype 96 samples. The current method was compared with a dHPLC method that we reported previously. Genotype frequencies at VKORC1 1173C>T for our Swedish population were 38% wild-type, 40% heterozygote and 22% homozygote. The frequency of the T-allele was 0.42, which exactly matches the frequency previously reported for Germans. The current method can be used to determine whether patients initiating warfarin therapy are carriers of SNP 1173 C>T that is strongly associated with low warfarin dose requirement.

Platelet PAR1 receptor density--correlation to platelet activation response and changes in exposure after platelet activation.

INTRODUCTION: A polymorphism (-14 A/T) affecting PAR1 expression on the platelet surface has recently been identified. A two-fold variation in receptor density, which correlated with the platelet response to PAR1-activating peptide (PAR1-AP), has been reported. MATERIALS AND METHODS: We used flow cytometry to measure the correlation between the number of PAR1 receptors and platelet activation. We also measured the changes in receptor exposure after platelet activation with PAR1-AP, ADP, PAR4-AP or a collagen-related peptide (CRP). RESULTS: In our study, the PAR1 receptor number varied almost four-fold, from 547 to 2063 copies/platelet (mean+/-S.D. 1276+/-320, n=70). The number of PAR1 receptors on resting platelets correlated to platelet fibrinogen binding and P-selectin expression following platelet activation with PAR1-AP (r(2)=0.30, p<0.01 and r(2)=0.15, p<0.05, respectively, n=36). The correlation was not improved by exclusion of the ADP-component from the PAR1-AP-induced response. We found a trend, but no statistically significant differences in PAR1 receptor number and platelet reactivity between A/A individuals and T/A or T/T individuals. Ex vivo activation with PAR1-AP decreased PAR1 surface exposure to 71+/-19% of the exposure on resting platelets (mean+/-S.D., p<0.01, n=19), while activation by ADP, PAR4-AP or CRP significantly increased the exposure, to 151+/-27%, 120+/-21% and 138+/-25%, respectively (n=11, 11 and 10). CONCLUSIONS: This study shows a large variation in PAR1 receptor number in healthy individuals, a variation correlated to the platelet activation response. We found a significant reduction in PAR1 surface exposure after adding PAR1-AP, while activation with ADP, PAR4-AP or CRP increased the exposure.

Plasma S/R ratio of warfarin co-varies with VKORC1 haplotype.

We recently reported that the low-dose VKORC1*2 haplotype is an important genetic determinant for warfarin dose requirement and is associated with difficulties to attain stable therapeutic prothrombin time--International Normalized Ratio in patients undergoing anticoagulation therapy. The aim of this study was to investigate whether patients with VKORC1*2 compared with patients carrying high-dose haplotypes VKORC1*3 or VKORC1*4 had different warfarin S/R ratios in their plasma, and whether that was related to CYP2C9 variants CYP2C9*2 and CYP2C9*3 or other factors. Samples from patients previously haplotyped for VKORC1 and measured for plasma warfarin concentration were genotyped for the CYP2C9 variants CYP2C9*2 and CYP2C9*3. Nonparametric statistical analysis was performed to elucidate whether there was any significant difference in the warfarin S/R ratio between the two patient groups. Our result shows that there is a significant difference (P<0.01) in warfarin S/R ratios between VKORC1*2 and VKORC1*3 or VKORC1*4 patients. This difference did not originate from CYP2C9 variants CYP2C9*2 and CYP2C9*3. We speculate that VKORC1 haplotypes possibly are linked to some unidentified factors involved in the metabolic clearance of warfarin enantiomers. Dose-dependent variations in (S)-warfarin and (R)-warfarin clearance in these patients can also be a probable explanation for the difference in warfarin S/R ratios.