Exploring genome-wide DNA methylation profiles altered in hepatocellular carcinoma using Infinium HumanMethylation 450 BeadChips

Hepatocellular carcinoma (HCC) incidence has increased in the US and also has one of the fastest growing death rates of any cancer. The purpose of the current study was to discover novel genome-wide aberrant DNA methylation patterns in HCC tumors that are predominantly HCV-related. Infinium HumanMethylation 450K BeadChip arrays were used to examine genome-wide DNA methylation profiles in 66 pairs of HCC tumor and adjacent non-tumor tissues. After Bonferroni adjustment, a total of 130,512 CpG sites significantly differed in methylation level in tumor compared with non-tumor tissues, with 28,017 CpG sites hypermethylated and 102,495 hypomethylated in tumor tissues. Absolute tumor/non-tumor methylation differences ≥ 20% were found in 24.9% of the hypermethylated and 43.1% of the hypomethylated CpG sites; almost 10,000 CpG sites have ≥ 30% DNA methylation differences. Most (60.1%) significantly hypermethylated CpG sites are located in CpG islands, with 21.6% in CpG shores and 3.6% in shelves. In contrast, only a small proportion (8.2%) of significantly hypomethylated CpG sites are situated in islands, while most are found in open sea (60.2%), shore (17.3%) or shelf (14.3%) regions. A total of 2,568 significant CpG sites (2,441 hypermethylated and 127 hypomethylated) covering 589 genes are located within 684 differentially methylated regions defined as regions with at least two significant CpG sites displaying > 20% methylation differences in the same direction within 250-bp. The top 500 significant CpG sites can significantly distinguish HCC tumor from adjacent tissues with one misclassification. Within adjacent non-tumor tissues, we also identified 75 CpG sites significantly associated with gender, 228 with HCV infection, 17,207 with cirrhosis, and 56 with both HCV infection and cirrhosis after multiple comparisons adjustment. Aberrant DNA methylation profiles across the genome were identified in tumor tissues from US HCC cases that are predominantly related to HCV infection. These results demonstrate the significance of aberrant DNA methylation in HCC tumorigenesis.     

The length of CpG islands is associated with the distribution of Alu and L1 retroelements

Alu and L1 retroelements have been suggested to initiate the spread of CpG methylation. In this study, the spread of CpG methylation was estimated based on the distance between the CpG islands and the nearest retroelements. All human genes (23,116) were examined and the correlations between the length of the CpG islands and the distance and density of the confronting retroelements were examined using nonoverlapping 5-kb windows. There was a linear relationship between the length of the CpG islands and the density of the Alu elements and an inverse relationship between the CpG islands and the L1 elements located more distantly, suggesting a suppressive effect of the Alu's on the spread of L1 methylation. Methylation analysis of the transitional CpG sites between the CpG islands and the nearest retroelements upstream of 16 genes was then carried out using DNA preparations from 11 different human tissues. Methylation-variable transitional CpGs were observed for the selected genes and the different tissues.     

Genetic variation in the developmental regulation of cortical avpr1a among prairie voles

Early experiences can have enduring impacts on brain and behavior, but the strength of these effects can be influenced by genetic variation. In principle, polymorphic CpGs (polyCpGs) may contribute to gene‐by‐environment interactions (G × E) by altering DNA methylation. In this study, we investigate the influence of polyCpGs on the development of vasopressin receptor 1a abundance in the retrosplenial cortex (RSC‐V1aR) of prairie voles (Microtus ochrogaster). Two alternative alleles (‘HI’/‘LO’) predict RSC avpr1a expression, V1aR abundance and sexual fidelity in adulthood; these alleles differ in the frequency of CpG sites and in methylation at a putative intron enhancer. We hypothesized that the elevated CpG abundance in the LO allele would make homozygous LO/LO voles more sensitive to developmental perturbations. We found that genotype differences in RSC‐V1aR abundance emerged early in ontogeny and were accompanied by differences in methylation of the putative enhancer. As predicted, postnatal treatment with an oxytocin receptor antagonist (OTA) reduced RSC‐V1aR abundance in LO/LO adults but not their HI/HI siblings. Similarly, methylation inhibition by zebularine increased RSC‐V1aR in LO/LO adults, but not in HI/HI siblings. These data show a gene‐by‐environment interaction in RSC‐V1aR. Surprisingly, however, neither OTA nor zebularine altered adult methylation of the intronic enhancer, suggesting that differences in sensitivity could not be explained by CpG density at the enhancer alone. Methylated DNA immunoprecipiation‐sequencing showed additional differentially methylated regions between HI/HI and LO/LO voles. Future research should examine the role of these regions and other regulatory elements in the ontogeny of RSC‐V1aR and its developmentally induced changes.