ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article
Updated

Correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes

[version 3; peer review: 1 approved, 2 not approved]
PUBLISHED 11 Sep 2013
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Background miRNA regulation of target genes and promoter methylation are known to be the primary mechanisms underlying the epigenetic regulation of gene expression. However, how these two processes cooperatively regulate gene expression has not been extensively studied.
Methods
Gene expression and promoter methylation profiles of 270 distinct human cell lines were obtained from Gene Expression Omnibus. P-values that describe both miRNA-targeted-gene promoter methylation and miRNA regulation of target genes were computed using the MiRaGE method proposed recently by the author.
Results
Significant changes in promoter methylation were associated with miRNA targeting. It was also found that miRNA-targeted-gene promoter hypomethylation was related to differential target gene expression; the genes with miRNA-targeted-gene promoter hypomethylation were downregulated during cell senescence and upregulated during cellular differentiation. Promoter hypomethylation was especially enhanced for genes targeted by miR-548 miRNAs, which are non-conserved, primate-specific miRNAs that are typically expressed at lower levels than the frequently investigated conserved miRNAs. miRNA-targeted-gene promoter methylation may also be related to the seed region features of miRNA.
Conclusions It was found that promoter methylation was correlated to miRNA targeting. Furthermore, miRNA-targeted-gene promoter hypomethylation was especially enhanced in promoters of genes targeted by miRNAs that are not strongly expressed (e.g., miR-548 miRNAs) and was suggested to be highly related to some seed region features of miRNAs.

Updated Changes from Version 2

Update from version 2 Abstract modified: The word "affected" was removed as the purpose of this paper is simply reporting correlations, not causalities.A figure to illustrate analyses: A figure to illustrate analyses was added (Fig. 1). A figure to show data trends: In order to show data trends, the comparison of distribution of promoter methylation between miRNA-targeted genes and non-miRNA-targeted genes was added (additional file 7). A figure to visualize correlation: In order to visualize correlations between miRNA-mediated regulation of target genes and miRNA-targeted promoter methylation, scatter plots were added (additional file 8). Correlation between miRNA-targeted promoter methylation and seed region features: In order to illustrate a newly identified relation between miRNA-targeted promoter methylation and seed region features, many materials were added: Table 11, Figures 2 and 3, additional file 9, sections "Correlation between seed region features and miRNA-targeted-gene promoter methylation" and "miRNA-tageted-gene promoter methylation is correlated to miRNA seed region features".Treatment of shared miRNA seed sequences: a section "Treatment of shared miRNA seed sequences" and an Appendix "Treatment of miRNAs sharing same seed sequence for P-values computation", Figures S1 and S2, and Table S1 were added in order to explain why shared miRNA seed sequences were not considered.Fig. 1 in Ver.2 is now renumbered as Fig. 4, due to the above mentioned material additions.

See the author's detailed response to the review by Denis Dupuy
See the author's detailed response to the review by Pål Sætrom

Introduction

The epigenetic regulation of gene expression1 has recently attracted the interests of many researchers. Epigenetic modifications regulate gene expression without modifying DNA sequences. Examples include promoter methylation2, histone modification3, the binding of transcription factors to gene promoter regions4, and miRNA regulation of target genes5.

Promoter methylation and miRNA regulation of target genes are particularly important in the epigenetic regulation of gene expression. Promoter methylation is relatively stable, long term, and in some cases, heritable. It is generally believed that genes with hypermethylated promoters are repressed. In addition, there is mounting evidence that DNA methylation is involved in the development and progression of certain disease states. For example, aberrant methylation in cancer is frequently observed6, and the distinct patterns of promoter methylation between monozygotic (MZ) twin pairs have also been found to result in different health conditions7. In contrast to DNA methylation, miRNA regulation of target genes is more flexible and can change even during cellular differentiation8. miRNA expression is often tissue-specific, and similar to DNA methylation, miRNA expression has been linked to human disease9. Thus, although miRNA-directed gene regulation is thought to result in subtle changes, it is generally believed that miRNAs are involved in many important biological processes ranging from cell division to aging.

Although DNA methylation of miRNA promoters has been studied extensively (e.g., with respect to tumor formation10), the relationship between promoter methylation and miRNA regulation of target genes has not been thoroughly investigated. One likely reason for this is that the regulation of gene expression by promoter methylation is a form of pre-transcriptional control, whereas miRNA regulation of target genes is a form of post-transcriptional control, with the former taking place inside the nucleus and the latter outside the nucleus (cytoplasm). Thus, these two mechanisms are separated by both time and space, and as a result, there have not been plausible biological reasons to suspect that promoter methylation and miRNA-mediated gene regulation operate in concert.

However, Su et al.11 recently found that miRNAs have a tendency to target genes with hypomethylated promoters. To my knowledge, their study was the first report suggesting coregulation of gene expression by promoter methylation and miRNAs. In addition, Sinha et al. reported that gene promoters with high CpG content were more often targeted by miRNAs12. Saito and Sætrom also discussed the relationships between miRNA-mediated gene regulation and various features of target genes, but they did not consider the methylation status of target gene promoters13. Although the study of Su et al. represents the first evidence of a direct link between promoter methylation and miRNA regulation, the biological significance of their findings is not clear. In this study, I report that promoter methylation is associated with miRNA-targeting; that is, the amount of methylation observed at a given gene promoter is dependent on whether that gene is also a target of miRNA regulation. Furthermore, miRNA-targeted-gene promoter methylation is also related to how miRNAs regulate target gene expression. In particular, I reveal that miR-548 miRNAs target genes are associated with highly hypomethylated promoters. Finally, the data presented here indicate that miRNA-targeted-gene promoter methylation is related to the seed region features of miRNAs.

Methods

An overview of the pipeline used for data processing is presented in Figure 1.

39bce5b3-886e-4f2d-a736-8473fa282275_figure1.gif

Figure 1. Schematic representing data processing.

Data processing steps. miRNA expression, target gene expression/methylation and seed region features of miRNAs were input dataset and were processed using MiRaGE method (for expression and methylation). Obtained P-values by MiRaGE method were compared with each other, miRNA expression and seed region features and correlation coefficients were computed.

Promoter methylation profiles

In this study, I used publically available promoter methylation profiles from various resources obtained from GEO ID: GSE3065314. This included 283 human promoter methylation profiles for distinct cell lines, ranging from hESC to various somatic samples, measured using the HumanMethylation27 BeadChip (Il-lumina), which provides an efficient approach for surveying genome-wide DNA methylation profiles. The HumanMethylation27 panel targets CpG sites located within proximal promoter regions of transcription start sites (TSS). Thus, it was suitable for the purpose of this study. Promoter methylation profiles (GEO ID: GSE30653) also included data from both IMR90 and MRC5 cell lines, which were used to investigate relationships between promoter methylation and previously reported miRNA regulation and miRNA expression profile data15,16. Promoter methylation profiles in both BG02 and BG03 were also included in this study, and were compared to miRNA regulation and miRNA expression profile data (see below).

Additional promoter methylation profiles in IMR90 cell lines were obtained from GEO IDs, GSM86800814, GSM73994017, and GSM37544218. They were compared to IMR90 promoter methylation profiles (GEO ID: GSM760387 within GSE30653). GSM868008 was included in the GEO ID, GSE31848, which were generated using the Illumina HumanMethylation450 BeadChip. This BeadChip allowed us to interrogate > 485000 methylation sites per sample at single-nucleotide resolution. In addition, because this array also includes CpG sites outside of promoter regions, I restricted probes to a subset labeled as either TSS200 or TSS1500. Data from GSM739940 includes IMR90 promoter methylation profiles measured by the Illumina HumanMethylation27 BeadChip. However, because these data were generated by a different research group than that of GSE30653, I tested profiles from this dataset to confirm that obtained results were not research group dependent. Finally, I also used methylation profile data from GSM375442, which were generated using next generation sequencing (NGS). CpG methylation profiles from promoter regions were extracted using Bismark Software (Ver. 0.7.4)19 (see below); promoter regions were defined as nucleotide positions between -200 and +1200 basepairs (bp) from transcription start sites (TSSs).

mRNA and miRNA expression profiles

In order to compare miRNA-targeted-gene promoter methylation with target gene miRNA regulation and miRNA expression profile data from BG02 and BG03 cell lines, both miRNA and mRNA profiles were obtained from GEO ID: GSE1447320. Gene (mRNA) expression profiles of undifferentiated and differentiated BG02 and BG03 cell lines were obtained from the GEO IDs, GSM551204 and GSM551206, and GSM551216 and GSM551218, respectively. Corresponding miRNA expression profiles of these two cell lines were obtained from the GEO IDs, GSM361147 and GSM361271 (BG02) and GSM361288 and GSM361289 (BG03). Raw data files were downloaded for further analysis and were normalized so as to have a mean of 0 and a variance of 1.

Investigation of miRNA-targeted-gene promoter methylation

In order to infer miRNA-targeted-gene promoter methylation, I employed the MiRaGE method21 (see below). The MiRaGE method, which was implemented on a public domain MiRaGE server and Bioconductor MiRaGE package, was first used to infer the contribution of miRNA to the measured expression levels. This software was designed to accept expression profiles of the target genes in question; however, in this study, I used this method to infer miRNA-targeted-gene promoter methylation by substituting gene expression profiles with the promoter methylation profiles of each gene.

I first prepared a control dataset (pseudo) in which the expression level of all genes was assigned a value of 1. Then, the amount of methylation at each gene was used in place of the values of treatment data set. Although, the ratio of the number of methylated sites to the total number of methylated and non-methylated sites is typically used to describe promoter methylation levels, I employed a method in which total methylation values were used. I used this approach because I found that P-values computed when using methylation data were more strongly correlated to the P-values calculated from target gene miRNA regulation data (see below), which is likely due to the fact that the frequency of CpGs is also related to miRNA targeting12; i.e., genes with promoters that contain more CpGs were more often targeted by miRNAs as mentioned above. Using this procedure, I attributed two P-values to each miRNA, one expressing the degree of promoter hypermethylation, and the other expressing the degree of promoter hypomethylation. The approach used to compute P-values representing promoter methylation are described below for each of the different methodologies and/or datasets used.

GSM868008. Promoter methylation profiles used to replace “gene expression” values within the MiRaGE method were

xcontrol,i=1,(pseudo) control set

xtreat,i=exp(M0i),treatment set

where M0i represented the scaled values of signal_B (intensity estimated of methylated DNA), which were expressed as the amount of promoter methylation of ith gene,

            M0i  ≡  Mi − MiσMi   〈Mi〉   ≡  1NiMi             σMi  ≡  1Ni(Mi − Mi)2,

where N was the total number of genes considered and Mi was the raw value of signal_B. This signified that the amount of promoter methylation was scaled so as to have a mean 〈M0i〉 of 0 and a standard deviation σM0i of 1. exp was applied in this instance because I wanted to consider the amount of methylation rather than the ratio of methylation. Because P-values were computed after the pair of input values were transformed to a logarithmic ratio, substituting 1 in the control dataset and an exponential value in the treatment dataset based on raw values resulted in the usage of raw values of differential expression/promoter methylation (see below).

GSE30653. Promoter methylation profiles used to replace “gene expression” values within the MiRaGE method were

xcontrol,i=1,(pseudo) control set

xtreat,i=Mi+Ci,treatment set

where Ci took on the value of 1 only when Mi = 0; otherwise, it took on the value 0, so as to avoid infinite values after transformation to the logarithmic ratio.

GSM739940. Promoter methylation profiles used to replace “gene expression” values within the MiRaGE method were

xcontrol,i=1,(pseudo) control set

xtreat,i=exp(βi),treatment set

where βi was the ratio of methylated sites to unmethylated sites,

βiMi+CMi+Ui+C,

where Ui was the signal from unmethylated sites (signal_A) and C was the regulation constant, which typically took the value of 100. Since only β values were deposited in the public datasets used, their use could not be avoided; however, as a result, the correlation with target gene miRNA regulation was substantially decreased. An explanation for this is noted above.

GSM375442. Promoter methylation profiles used to replace “gene expression” values within the MiRaGE method were

xcontrol,i  =   1, (pseudo) control setxtreat,i     =   exp (Mimax(Mi)) ,  treatment set

where max(Mi) was the maximum value of Mi and Mi was computed in this case as follows:

Mijpromoter of gene iyj,

where yj, 0 ≤ yj ≤ 100 was the percentage of methylation at site j, which was computed using the Bismark Software19 (see below). The summation was taken over the length of the promoter region as defined above (i.e., between -200 bp and +1200 bp from the TSS).

Methylation computation by Bismark Software

The following command line inputs were used to generate methylation values of CpG sites within the Bismark Software package19.

% bismark_genome_preparation \

– – path_to_bowtie bowtie_dir \

– – verbose ./hg19/ &

% R

>x <– scan("GSM375442_CpgMIP–IMR90.seq.txt",

sep="\n", what=character(0))

>write.table(file="sequence.fa",

paste(paste(">p", 1:length(x), sep=""),

x, sep="\n"), sep="\n", row.names=F,

quote=F, col.names=F)

>q()

% bismark ./hg19/ \

– – path_to_bowtie bowtie_dir \

– – bowtie2 –f sequence.fa

% methylation_extractor –s – – comprehensive \

sequence.fa_bt2_bismark.sam

% genome_methylation_bismark2bedGraph_v3.pl \

CpG_content_sequence.fa_bt2_bismark.sam.txt \

> sequence.fa_bt2_bismark.sam.bed

where bowtie_dir is the directory where bowtie2 is installed. In this study R22 was used, but scripts executed by R can also be performed using other alternative languages. GSM375442_CpgMIP-IMR90.seq.txt is a file downloaded from GEO and sequence.fa_bt2_bismark.sam.bed includes methylation percentages of each CpG site, yj, which was explained above.

Inference of miRNA-targeted-gene promoter methylation/miRNA regulation of target genes using the MiRaGE method

The inference of miRNA-targeted-gene promoter methylation/miRNA regulation of target genes was carried out using the MiRaGE method, which has been described previously21.

Although the MiRaGE method is typically used for datasets with two experimental conditions, each of which contains more than one replicate, I used this method for instances in which each condition consisted of only a single replicate. In each case, I based the analysis on the premise that for a given gene i, there were a pair of gene expression datasets or promoter methylation profile datasets, which were measured under a control condition xcontrol,i and a treatment condition xtreat,i. From this, I computed the logarithmic ratio

Δxilog(xtreat,ixcontrol,i)

In such an example, when a difference between raw values is favorable, exponential values exp(xi) can be used instead of xi, in which case, I got

Δxi=log{exp(xtreat,i)exp(xcontrol,i)}=xtreat,i-xcontrol,i

When computing P-values that rejected the null hypothesis, by using the alternative hypothesis that ∆xis of the target genes of the miRNA m are less (greater) than those that are off-target but form a target of other miRNAs, I computed

Pm<(>)P[{Δxi|iGm}<(>){Δxi|iGm'}]

where P[A < (>) B] represented P-values computed by statistical tests when sets A and B were compared. The tests implemented within the MiRaGE Server/package are the t-test, Wilcoxon rank sum test, and Kolmogorov-Smironov test.

Thus, determining whether A < (>)B depends on the selected statistical test used. I used Gm to represent the set of genes targeted by miRNA m and m as the intersection of the set of off-target genes of miRNA m and the set of target genes of all other miRNAs. It should be noted that genes that were not the targets of any miRNAs were totally excluded from the analysis; however, all miRNAs were considered and no exclusions based on conservation were applied. When inferring promoter methylation, xcontrol,i = 1 and xtreat,i were used to represent the amount of promoter methylation. When inferring miRNA regulation of target genes during cell senescence in IMR90 and MRC5 cell lines, xcontrol,i was used to represent gene expression of young cell lines and xtreat,i was used to represent gene expression in senescent cell lines. When inferring regulation of target genes during differentiation in BG02 and BG03 cell lines, xcontrol,i was used to represent gene expression in undifferentiated cell lines and xtreat,i for gene expression in differentiated cell lines.

Correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes

I used two types of P-values: P<(>)m:methy, which corresponded to the miRNA-targeted-gene promoter methylation; and P<(>)m:regul, which corresponded to the miRNA regulation of target genes of miRNA m. When P<(>)m:methy[regul] was small enough, the target genes of miRNA m were significantly hypermethylated or hypomethylated and thus downregulated or upregulated. In order to see if these two types of P-values were correlated, I computed various correlation coefficients:

ρPearson     ≡  ρ‌[Pm:methy,Pm:regul]ρlogPearson     ≡  ρ‌[log(Pm:methy), log(Pm:regul)]ρSpearman  ≡  ρ‌[rank(Pm:methy), rank(Pm:regul)]

and accompanied P-values to reject null hypothesis that ρ = 0 by using the alternative hypothesis that ρ ≠ 0. I used ρ[a, b] as the Pearson's correlation coefficients between a and b and rank(xm) as the rank order of xm among {xm}. Pm:z where z ∈ {methyl, regul} was either P<m, P>m, 1 – P<m or 1 – P>m. Thus, there were 4 × 4 = 16 possible combinations of Pm:methy and Pm:regul. ρs of the Kolmogorov-Smirnov test and ρlogPearsons for all tests changed when P<(>) was replaced with 1 – P>(<) because P<(>) ≠ 1 – P>(<) for the Kolmogorov-Smirnov test and log(P<(>)) ≠ log(1 – P>(<)) for all tests. Thus, optimal combinations of the maximum absolute correlation coefficients were employed.

The reciprocal relationship between miRNA regulation of target genes and miRNA expression during differentiation in BG02 and BG03

In contrast to the cell senescence study15 in which miRNA expression was investigated by NGS, only microarray measurements were available for differentiation in BG02 and BG03 cell lines. Due to issues related to the accuracy and quality of microarray data and the relatively low levels of miRNA expression, very few miRNAs were found to be differentially expressed between undifferentiated and differentiated cell lines. Thus, miRNAs that were differentially expressed between undifferentiated and differentiated cell lines were selected based on two criteria before being subjected to further analyses:

  • Absolute differential expression | xdiff,ixundiff,i |> ∆xc, where xdiff,i and xundiff,i were the normalized expression of gene i of differentiated and undifferentiated cell lines, respectively. ∆xc was set as the threshold value that could be used to select genes associated with significant differential expression during differentiation.

  • For this method, genes were considered to be significantly upregulated or downregulated based on an adjusted P-value threshold (P < 0.05) by using the BH criterion23. Here, P-values were computed using a t-test between two sets of probe values attributed to gene i.

After ∆xc was suitably selected, the correlation coefficient between differential expression xdiff,ixundiff,i and log (P<m:regul) were computed. Positive values were taken to indicate a reciprocal relationship between miRNA expression and miRNA-mediated regulation of target genes; smaller P<m:regul were assumed to signify that target genes were upregulated (and vice versa). Thus, a reciprocal relationship required that miRNA were downregulated, i.e., xdiff,ixundiff,i < (>)0, which should result in a positive correlation between xdiff,ixundiff,i and log (P<m:regul). Since miRNA names in GSE14473 were old (miRBase, release 9.1), they were converted to miRNA names used in the present version of the MiRaGE software package (miRBase release 18) using the miRConverter implemented in miRSystem24.

Ranking of miRNAs having target genes with miRNA-targeted-hypomethylated promoters

miRNAs were ranked based on P>m in each of the 283 samples in GSE30653, after excluding 12 control samples. Each miRNA was ranked within each sample, and then, the cumulative rank was determined for all 270 samples.

Ranking of miRNA regulation of target genes

miRNAs were ranked based on P-values for either downregulation during cell senescence (IMR90 and MRC5) or upregulation during differentiation (BG02 and BG03) for each statistical test. Then, each miRNA was ranked based on the order summed up over three statistical tests.

Comparison of promoter methylation between miRNA-target genes and miRNA-non-target genes

I compared promoter methylation between miRNA-target genes and miRNA-non-target genes for each of the 270 samples in GSE30653, excluding 12 control samples. Statistical tests used were the two-sided t-test, Wilcoxon rank sum test, and Kolmogorov-Smirnov test, resulting in 270 P-values for each test.

P-value computation for miRNA-targeted-gene promoter methylation after random permutation

In order to see if P-values for miRNA-targeted-gene promoter methylation, P>m, changed after random gene permutation, gene IDs were randomly permuted for each of the 270 samples. Then, the number of significant miRNAs whose target gene promoter was significantly hypomethylated (P-values computed with t-test and corrected using the BH criterion23 were less than 0.05) was counted. The numbers were averaged over one hundred independent random permutations.

Significance of correlation between miRNA-targeted-gene promoter methylation and miRNA-mediated regulation of target genes after random gene permutation

In order to see if the correlation between miRNA-targeted-gene promoter methylation and miRNA-mediated regulation of target genes changed after random gene permutation while conserving the correlation between promoter methylation and gene expression, promoter methylation and gene expression profiles needed to be permuted simultaneously. However, because gene expression and promoter methylation were measured by distinct microarrays, the number of probes attributed to each gene in some cases differed between the microarray used for gene expression measurement and promoter methylation. Also, multiple CpG sites may be attributed to each gene, while gene expression is unique to each gene. In order to overcome this difficulty, multiple probes or CpG sites were first grouped based on overlapping RefSeq mRNA IDs. RefSeq gene IDs were then shuffled and P-values were computed using MiRaGE, which employs a t-test for both miRNA-targeted-gene promoter methylation and miRNA-mediated regulation of target genes, before and after random permutation. The methylation data sets used were GSM868008 (for IMR90 cell line) and GSE30653 (for MRC5, BG02 and BG03). ρPearson were computed with Pm:methy = P<m and Pm:regul = P<m. Gene expression profiles reported previously for both IMR90 cell lines and MRC5 cell lines were used15. The maximum and minimum values of ρPearsons among 100 independent random permutations are presented. The geometric average for over 100 random permutations of P-values attributed to ρPearson was also computed. In addition, P-values associated with ρPearson before random permutation was estimated based on 100 random permutation ensembles assuming a normal distribution. The normality of the distribution of correlation coefficients from 100 random permutations was checked using both the Shapiro-Wilk normality test and the Kolmogorov-Smirnov test (no significant differences from a normal distribution were detected).

Correlation between seed region features and miRNA-targeted-gene promoter methylation

miRNA seed region features were downloaded from miRmap25. According to the target gene table of miRNA26, seed region features attributed to each miRNA were obtained as averages over target genes. Correlations were computed in each cell line between seed region features and P>m:methy computed by t-test. The distributions over all cell lines are presented.

P-values that represent the significance of differences in seed region features between sets of miRNAs (i.e., miRNA-548 miRNAs or let-7 miRNAs) were computed using the Wilcoxon rank sum test. Corresponding binary logarithmic ratios were also computed.

Results and discussion

Significant promoter hypomethylation of miRNA target genes

Based on the inference of data produced using the MiRaGE method21 (see Methods), the promoters of genes that are targets of 70–90% of human miRNAs were significantly hypomethylated, dependent on the statistical tests used, and the definition of promoter methylation levels: the β-value or the amount of methylation (see Table 1 and Additional files 1 to 6). This finding was consistent with conclusions made by Su et al.11, who stated that miRNAs had a tendency to target genes with hypomethylated promoters. Given that promoter methylation patterns do not change drastically between certain cell types, the amount of miRNA-targeted-gene promoter methylation is also highly cell-type-independent. Mean correlation coefficients of P>m ranged between 0.8 and 0.9, again depending on the statistical test used and the definition of promoter methylation levels (Table 2).

Table 1. miRNA-targeted-gene promoter hypomethylation.

Amount of methylation
Statistical tests Number of miRNAs
Not significant Significant
t-test13605438453774%
Wilcoxon rank sum test14457237601972%
Kolmogorov-Smirnov test12999839059375%
β -values
Statistical tests Number of genes
Not significant Significant
t-test3939048120192%
Wilcoxon rank sum test4518847540391%
Kolmogorov-Smirnov test4195247863992%

The summation of the number of miRNAs that target genes with significantly hypomethylated promoters in each cell line. The total numbers of miRNAs and cell lines were 1921 and 270, respectively (the number of cell lines was out of 283 in GSE30653, excluding 12 control samples). Thus, the total number of combinations of an miRNA and a cell line is 1921 × 271 = 520591. Promoter methylation levels were defined as either the amount of methylation or the β-values. Observed P-values based on the BH criterion23, P>m of < 0.05 were used to signify the promoters of the genes targeted by miRNAs that significantly hypomethylated. P-values were adjusted using the p.adjust function in R22. A complete list of raw (not adjusted) P-values is available in Additional files 1 to 6. Cell line names can be also identified as column names of the Additional files 1 to 6.

Table 2. Mean correlation coefficients between P-values of miRNA-targeted-gene promoter hypomethylation of samples.

Mean Pearson’s correlation coefficients of P>m between pairs of cell lines.

Statistical test Averaged correlation coefficients
Amount of methylation
t-test0.9112043
Wilcoxon rank sum test0.8981256
Kolmogorov-Smirnov test0.8263153
β-values
t-test0.8919359
Wilcoxon rank sum test0.8832241
Kolmogorov-Smirnov test0.8257556

In order to confirm that the observed miRNA-targeted-gene promoter methylation was really correlated to miRNAs targeting, I compared promoter methylation of miRNA target genes with that of miRNA non-target genes. In order to check this point, I have computed P-values that compare promoter methylation between miRNA target genes and miRNA non-target genes (Table 3). Histograms that show promoter methylation in miRNA target genes and miRNA non-target genes for each cell line are available as Additional file 7. I found that the difference was strongly significant. P-values were clearly small enough to be significant even after correction for multiple comparisons. These significant differences corroborate the findings of Su et al.11 and support the argument that observed miRNA-targeted-gene promoter methylation was correlated to miRNA targeting.

Table 3. P-values that represent the difference of promoter methylation between miRNA target genes and miRNA non-target genes.

t-testWilcoxon rank sum testKolmogorov-Smirnov test
MinimumMaximumMinimumMaximumMinimumMaximum
1.38 × 10–60 4.66 × 10–25 4.26 × 10–105 1.56 × 10–50 00

The minimum and maximum values are among 270 cell line samples. 0s for Kolmogorov-Smirnov test mean that they cannot be distinct from 0 within numerical accuracy.

Because promoter hypermethylation and hypomethylation can have contrasting functional effects, these modifications are often investigated as separate data classes. However, in the context of miRNA-mediated gene regulation, there are many known cases in which the functional outcomes of miRNA targeting can vary depending on the presence of additional factors. For example, some miRNAs can remove Ago proteins that bind other miRNAs27. In such cases, target genes of these miRNAs are assumed to be upregulated, as miRNAs are generally believed to suppress the expression of their target genes. miRNA sponges can also inhibit miRNA suppression of target genes28. miRNA sponges absorb miRNAs that bind to the RNA-induced silencing complex (RISC). Because of this effect, the expression of miRNA targets tends to increase with increasing concentrations of miRNA sponges. The biological mechanisms that mediate the interactions between miRNA and methylated promoters, and the role of other factors are unknown, thus, it is not always necessary to distinguish the mechanisms underlying the hyper- and hypomethylation of miRNA target genes.

Another possible issue is that miRNA-targeted-gene promoter methylation may represent an artifact caused by the complex algorithm employed by MiRaGE. In order to exclude this possibility, I computed miRNA-targeted-gene promoter methylation after random permutation (Table 4), from which I determined that there was no significant miRNA-targeted-gene promoter methylation. Thus, it is apparent that random permutation destroys miRNA-targeted-gene promoter methylation completely, demonstrating that these data provide robust evidence for the occurrence of miRNA-targeted-gene promoter methylation.

Table 4. miRNA-targeted-gene promoter hypomethylation after random gene permutation.

Amount of methylation
Statistical testsNumber of miRNAs
Not significantSignificant
t-test52058561.1 × 10–3%
Wilcoxon rank sum test52058835.8 × 10–4%
Kolmogorov-Smirnov test52058561.1 × 10–3%

The summation of the number of miRNAs that target genes with significantly hypomethylated promoters in each cell line after random permutation. See Table 1 for other notations. The values are averaged over 100 random permutation ensembles.

Significant correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation

In order to see if miRNA-targeted-gene promoter methylation was genuinely related to the miRNA regulation of target genes, I compared P-values of miRNA-targeted-gene promoter methylation to P-values of miRNA regulation of target genes during the senescence of IMR90 and MRC5 cell lines15,16 and during the differentiation of BG02 and BG03 cell lines20. It was clear that promoter methylation and the regulation of target genes by miRNA were significantly correlated during both cell senescence and differentiation (Table 5), despite the fact that correlation coefficients exhibited opposite directionalities in cell senescence and differentiation. This means that genes with miRNA-targeted-gene promoter hypomethylation are downregulated during cell senescence, but upregulated during differentiation.

Table 5. Correlation coefficients between P-values of miRNA-targeted-gene promoter hypomethylation and miRNA regulation of target genes.

Statistical testρPearsonρPearsonlogρSpearman
Cell senescence Pm:methy = P>m vs. Pm:regul = P>m
IMR90
t-test0.330.400.41
P-values***
Wilcoxon rank sum test0.280.350.33
P-values***
Kolmogorov-Smirnov test0.310.720.61
P-values***
MRC5
t-test0.220.230.25
P-values***
Wilcoxon rank sum test0.530.760.73
P-values***
Kolmogorov-Smirnov test0.450.820.67
P-values***
Differentiation Pm:methy = P>m vs. Pm:regul = 1 – P<m
BG02
t-test-0.25-0.68-0.67
P-values***
Wilcoxon rank sum test-0.14-0.53-0.53
P-values2.1 × 10–9 **
Kolmogorov-Smirnov test-0.13-0.71-0.43
P-values5.5 × 10–12 **
BG03
t-test-0.31-0.57-0.54
P-values***
Wilcoxon rank sum test-0.15-0.23-0.24
P-values5.4 × 10–11 **
Kolmogorov-Smirnov test-0.13-0.22-0.20
P-values3.5 × 10–10 **

Various correlation coefficients between P-values of promoter methylation and regulation of target genes in IMR90 and MRC5 during cell senescence and in BG02 and BG03 during differentiation; P-values associated with the correlation coefficients are indicated. * indicates < 2.2 × 10–16. P>ms were employed for cell senescence while P>ms of promoter methylation and 1 – P>m of regulation of target genes were employed for differentiation as shown, as these combinations exhibited the most significant correlations. Genes with miRNA-targeted-gene promoter hypomethylation were recognized to be downregulated during cell senescence based on positive correlations and upregulated during differentiation based on negative correlations (for details, see Methods).

To confirm that the significant correlation noted above between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes was not a false-positive finding, I also tested for this association in IMR90 cell lines using a different microarray and NGS approach, as described previously18 (Table 6). Again, I observed a significant correlation between promoter methylation status and miRNA regulation, suggesting that the relationship found between these two mechanisms was independent of experimental condition. Scatter plots generated from data in Table 5 and Table 6 are presented in Additional file 8.

Table 6. Evaluation of the correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes.

Statistical testρPearsonρPearsonlogρSpearman
GSM868008Pm:methy = P>m vs. Pm:regul = P>m
t-test0.280.370.37
P-values***
Wilcoxon rank sum test0.180.320.28
P-values4.4 × 10–16**
Kolmogorov-Smirnov test0.160.700.54
P-values1.9 × 10–12**
GSM739940Pm:methy = P>m vs. Pm:regul = P>m
t-test0.110.400.34
P-values3.0 × 10–6**
Wilcoxon rank sum test0.070.440.27
P-values2.5 × 10–3**
Kolmogorov-Smirnov test0.070.790.51
P-values2.7 × 10–3**
GSM375442Pm:methy = P>m vs. Pm:regul = P>m
t-test0.190.170.26
P-values*1.5 × 10–14*
Wilcoxon rank sum test0.190.270.24
P-values2.2 × 10–16**
Kolmogorov-Smirnov test0.280.570.50
P-values***

Correlation coefficients between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes during cell senescence shown in Table 5 were evaluated by a comparison with the results produced by different array designs (GSM868008), distinct research groups, alternative measures of promoter methylation (GSM739940), and NGS (GSM375442). Independent of these factors, miRNA-targeted-gene promoter methylation was found to always be positively correlated to miRNA regulation of target genes (i.e., genes with miRNA-targeted-gene promoter hypomethylation were downregulated during cell senescence).

However, these findings may not represent a true correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation. For example, if the observed correlation is genuine, then the correlation between gene expression and promoter methylation should differ between miRNA target genes and non-target genes. Alternatively, a correlation between miRNA-mediated gene regulation and miRNA-targeted-gene promoter methylation could simply be a bi-product of the well-documented direct correlation between gene expression and promoter methylation. These issues arise from two assumptions: (1) the correlation between gene expression and promoter methylation is strong enough to mediate a correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation; and (2) the correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation should be strongly related to, or even directly caused by, the direct correlation between gene expression and promoter methylation. However, these two assumptions do not necessarily need to be true in all cases. For example, the correlation between gene expression and promoter methylation is not always strong and the correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation could be highly significant in instances in which the relationship between gene expression and promoter methylation is weak.

In order to test the first assumption, I computed the correlation between gene expression and promoter methylation, and revealed that the correlation between gene expression and promoter methylation is much weaker than the correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation (see the column entitled “gene expression vs promoter methylation” in Table 7). This suggests that the correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation is not likely the outcome of the correlation between gene expression and promoter methylation, but possibly vice versa. Next, I considered the second assumption, which also need not always be true. From a mathematical point of view, a correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation does not have to be directly related to the correlation between gene expression and promoter methylation. The reason is as follows. It is always possible to shuffle gene order while conserving the correlation between gene expression and promoter methylation, because the correlation does not change between pre- and post-permutation if pairs of promoter methylation and gene expression can be conserved. However, permutation inevitably alters both miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation, because permutation changes the set of genes targeted by each miRNA.

Table 7. The correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes after random gene permutation.

miRNA regulation vs miRNA-targeted-gene promoter methylationGene expression vs
promoter methylation
(P-value)
Before random permutation
(P-value)
After random permutation
MinimumMaximum(Mean P-value)
IMR900.39 (2.5 × 10–13)-0.1270.133(0.03)-0.016 (2.27 × 10–5)
MRC50.22 (1.7 × 10–7)-0.1050.128(0.07)-0.003 (0.77)
BG02-0.27 (2.0 × 10–7)-0.1410.078(0.01)-0.048 (9.03 × 10–8)
BG03-0.30 (2.3 × 10–14)-0.0860.148(0.02)0.021 (0.02)

ρPearson with Pm:methy = P<m and Pm:regul = P<m where P-values were computed by t-test. Before random permutation: correlations between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes after random permutation and associated P-values computed from 100 random permutations ensemble assuming the normal distribution. After random permutation: minimum and maximum correlations among 100 random permutation ensembles. P-values are computed via geometric mean over 100 random permutation ensembles. Gene expression vs promoter methylation: correlation between gene expression and promoter methylation. Please note that genes are permuted such that correlation between gene expression and methylation is conserved. P-values are computed by t-test.

Suppose we have N genes and M miRNAs. Even if we permute gene id g in both microarrays that measure gene expression and promoter methylation simultaneously, the correlation between gene expression and promoter methylation does not change, because pairs of gene expression and promoter methylation are conserved. However, Gm and G'm are always altered, as permutation inevitably removes some genes from Gm, and genes that were not included in Gm before permutation are added to Gm after permutation. Since both miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation are computed based on the comparison between Gm and G'm, miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation must be altered by the permutation.

In spite of the above discussion, one may still argue that a correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation may not be altered or may be altered only slightly by the permutation. This is due to the fact that the modification of Gm is necessary for the computation of both miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation, even though miRNA-targeted-gene promoter methylation itself is not likely to be significant because of permutation, as demonstrated in the previous section (see Table 4).

In order to refute this idea, I calculated the correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation after random permutation (Table 7) and determined that random permutation drastically reduces the correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation, although it does not alter the correlation between gene expression and promoter methylation, as discussed above. This suggests that the correlation between gene expression and promoter methylation can hardly be a cause of the correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation, because correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation can take either smaller or larger values while the correlation between gene expression and promoter methylation is unaltered.

All of the above findings and discussions, together with the findings and discussions in the previous section, strongly suggest that the significant correlation between miRNA-mediated regulation of target genes and miRNA-targeted-gene promoter methylation is neither an artifact generated by the complicated algorithm employed in MiRaGE, nor a bi-product of a direct correlation between gene expression and promoter methylation, but is itself a genuine biologically significant phenomenon.

Significant reciprocal relationship between miRNA expression and target gene expression

In addition to these analyses, I also confirmed a significant reciprocal relationship between miRNA expression and target gene expression. Because it is usually believed that miRNAs downregulate target gene expression, downregulated genes should be targeted by upregulated miRNAs and vice versa. In fact, the reciprocal relationship during the senescence process of IMR90 cells has been reported15,16. In this study, I observed a similar reciprocal relationship during the differentiation of BG02 and BG03 cells20 (Table 8).

Table 8. The reciprocal relationship between target gene regulation by miRNA and miRNA expression.

Cell linesΔxcCorrelation coefficientsP-value
BG020.000.181.9 × 10–4
BG030.040.510.04

The reciprocal relationship between target gene regulation by miRNA and miRNA expression was evaluated in BG02 and BG03 cell lines during differentiation. Given that only a limited number of miRNAs are differentially expressed during differentiation, miRNAs were screened as described in the Methods section before computing a correlation. Then the correlation coefficients between differential expression of miRNA and logarithmic value of P-values of miRNA regulation of target genes were computed. A positive correlation implied the reciprocal relationship, that is, if the target genes are expressive, then, the miRNA itself should be suppressive (for details, see Methods).

Thus, I conclude that miRNA-targeted-gene promoter methylation is not an artifact, but a genuine biological process. One may doubt this conclusion given the finding that miRNA-targeted-gene promoter methylation is cell-line-independent and miRNAs are believed to regulate target gene expression in a cell-line-dependent manner. Although this implies an apparent discrepancy, the correlation coefficients obtained, which were at most 0.7 to 0.8, indicate that miRNA-targeted-gene promoter methylation reflects, at most, 50% to 60% of miRNA regulation of target genes. Thus, the cell line specific miRNA regulation of target genes can act within the remaining 50% to 40%. In actuality, as can be seen in Table 5, the correlations between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes have opposite signs dependent on the biological processes being considered. Therefore, the fact that cell-line-independent and miRNA-targeted-gene promoter methylation partially governs cell line dependent miRNA regulation of target genes is not a discrepancy. At present, I have no hypotheses that explain why miRNA-mediated regulation of target genes is the opposite between cell senescence and cell differentiation despite the fact that miRNA-targeted-gene promoters are similarly hypomethylated. However, as discussed above, the miRNA-mediated regulation of target genes is often reversed because of the existence of unobserved factors. The underlying biological processes, which are unknown at present, will be elucidated in the future.

One should understand that my objective is not to insist that miRNA up/downregulation is associated with hyper/hypomethylation of the promoter. Nor can I confirm whether miRNA mediates promoter hyper-methylation or hypomethylation. To my knowledge, there are no experimental observations that illustrate how miRNA methylate target gene promoters in animal cells. However, in plants, it has been reported that mature miRNAs in nuclei can mediate DNA methylation29. Mature miRNAs can also be transferred to and reside within the nucleus30, although an explicit role in mediating promoter methylation has not observed. Details in animal cells are not yet known. In the present study, I observed correlations between gene expression and promoter methylation from an miRNA-centric viewpoint. The biological mechanism mediating this correlation, however, is not clear. It is plausible that an unknown biological factor regulates both gene expression and promoter methylation, and that the correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes may simply be an outcome of indirect regulation and promoter methylation. I cannot deny this possibility, but if the correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes represents a real phenomena and can be found in a wide range of biological processes, there must be a biological/mechanistic underpinning of this effect.

Significant promoter hypomethylation of genes targeted by miR-548 miRNAs

In order to better understand the biological significance of the findings regarding miRNA-targeted-gene promoter methylation, I sought to identify specific miRNAs that target genes with hypomethylated promoters. I found that most gene promoters targeted by miR-548 miRNAs were highly hypomethylated (Table 9; here, I define miRNAs called “miR-548” as miR-548 miRNAs). At present, little is known about miR-548 miRNAs, which likely results from the fact that miR-548 miRNAs are primate-specific miRNAs (ps-miRNAs). Because of this, ps-miRNAs were previously not expected to play critical roles in fundamental biological processes like cell differentiation and development, as it has been widely assumed that such basic biological processes are conserved among mammals. This assumption has resulted in the exclusion of ps-miRNAs from consideration as candidates for basic and important biological processes. Thus, it is thought that the roles of these ps-miRNAs may be limited to primate-specific biological processes. Another reason that miR-548 miRNAs have not been extensively investigated is that they are generally expressed at low levels; miRNAs that play critical roles are thought to be highly expressed. For example, although miRBase stores the number of short reads detected NGS for various miRNAs, miR-548 miRNAs have at most 100 short reads; this is in comparison to more abundant miRNAs like those from the let-7 family31, which are typically represented by several million reads. This makes it difficult to detect miR-548 miRNAs using microarray technology or sequencing. Actually, although there have been a few reports noting that miR-548 family members, including miR-548 miRNAs, play critical roles, these data come exclusively from quantitative PCR (qPCR) experiments32,33.

Table 9. Evaluation of miR-548 miRNAs in the context of miRNA-targeted-gene promoter hypomethylation.

Statistical test Number of miR-548 miRNAs
Within top 100 ranked Total number Mean rank
t-test4368216.2
Wilcoxon rank sum test4368214.0
Kolmogorov-Smirnov test4368215.0

68 miR-548 miRNAs were ranked as described in the Methods (promoter methylation was defined as the amount of methylation). Independent of the statistical methods used for computation of P>m, 43 out of 68 miR-548 miRNAs were ranked within the top-ranked 100 miRNAs. The mean rank of all miR-548 miRNAs is ~200 independent of the statistical method used. Given that the total number of miRNAs considered was 1921, these values are highly significant, P = 2 × 10–35.

In this study, I found that promoters of genes targeted by miR-548 miRNAs were hypomethylated. Because it is generally assumed that genes with hypomethylated promoters are actively expressed, the fact that promoters of target genes of miR-548 miRNAs are hypomethylated may mean that, compared to genes that are targets of other miRNAs, target genes of miR-548 miRNAs are more sensitive to changes in miRNA expression. This could explain why promoters of target genes of non-expressed miR-548 miRNAs are specifically hypomethylated.

miR-548 miRNAs represent a large set of miRNAs, with 68 family members; thus, although these miRNAs exhibit low levels of expression, the fact that there are many family members may imply that they have important biological roles. Indeed, it has been reported that potential target genes of miR-548 miRNAs play critical roles in various biological processes34,35. Since it is also suggested that they originate from transposable elements (TEs), miR-548 miRNAs exhibit high sequence similarities while maintaining significant sequence diversity34.

During senescence processes of IMR90 and MRC5 cell lines, the target genes of miR-548 miRNAs were significantly downregulated (Table 10), whereas during differentiation of the BG02 cell line, the target genes of miR-548 miRNAs were significantly upregulated (Table 10). The fact that I observed significant changes in expression of target genes without significant changes in expression of miRNAs that target these genes implicates a role of miRNA-targeted-gene promoter hypomethylation.

Table 10. Ranking of miRNA regulation of target genes.

Cell linesTarget genesNumber of miR-548 miRNAs
Within top 100 rankedTotal numberMean rankP-values
IMR90Downregulated2468438.29 × 10–14
MRC5Downregulated1668382.65 × 10–7
BG02Upregulated2468427.89 × 10–14
BG03Upregulated368762.47 × 10–1

miRNAs were ranked as described in the Methods. Target genes of miR-548 miRNAs were significantly downregulated during cell senescence (IMR90 and MRC5 cell lines) and upregulated during differentiation (BG02). Insignificant upregulation of target genes of miR-548 miRNAs in the BG03 cell line was possibly because of the weak correlation between miRNA expression and miRNA regulation of target genes (Table 8). Given that the total number of miRNAs considered was 1921, the values for IMR90, MRC5, and BG02 cell lines are highly significant.

Of course, the possibility remains that there are unknown biological factors that both mediate promoter methylaion of miR-548 target genes and regulate gene expression targeted by miR-548 miRNAs, and all findings here are just outcomes of this hidden biological processes. Although I cannot deny this possibility, even if there are some hidden biological processes, the fact is that miRNA-548 target genes are both methylated in promoters and their regulation is robust. Future work investigating the potential underlying mechanisms is warranted.

miRNA-targeted-gene promoter methylation is correlated to miRNA seed region features

Previous studies have reported that miRNA activity, as defined by target gene expression, is often not strongly correlated to miRNA expression. For example, Liang et al.36 tried to relate the expression levels of miRNA to their regulatory activities, but failed. Furthermore, Fu et al.37 found that miRNA expression and target gene expression are correlated for only a limited number of miRNAs. Thus, attempting to relate miRNA expression to miRNA-targeted-gene promoter methylation may also be unavailing. In addition to this, miRNA-targeted-gene promoter methylation is not cell line specific (Table 2), but miRNA expression itself is usually believed to be cell-type specific. Finally, correlations between miRNA and mRNA expression38 are weaker than that observed for miRNA-targeted-gene promoter methylation and miRNA-mediated regulation of target genes (Table 5 and Table 6), further suggesting that the relationship between miRNA activity and miRNA expression may not be straightforward.

Actually, Fu et al.37 found that miRNA activity is more reproducible than miRNA expression across different datasets. This may suggest that miRNA activity is functionally more important than miRNA expression. Following Fu et al.37, who found that the efficiency of miRNA-mRNA binding is positively correlated with miRNA activity, I attempted to find something other than miRNA expression that is related to miRNA-targeted-gene promoter methylation.

Briefly, seed region features were downloaded from miRmap25 and were attributed to each miRNA (for detail, see Method). Correlations between seed region features and miRNA-targeted-gene promoter methylation in each of the 270 cell lines were then computed (Figure 2). A full set of histograms are available as Additional file 9. Particular seed region features, mainly those related to binding affinity, were highly correlated to miRNA-targeted-gene promoter methylation (absolute values of correlation coefficients are 0.5). Based on this, I conjectured that miRNA-targeted-gene promoter methylation is not mediated by miRNA expression, but instead by the strength of the binding affinity of miRNA seed regions in target genes.

39bce5b3-886e-4f2d-a736-8473fa282275_figure2.gif

Figure 2. Boxplot of correlation coefficients between miRNA-targeted-gene promoter methylation and seed region features.

Boxplot of correlation coefficients between miRNA-targeted-gene promoter methylation (P-values computed by t-test) and seed region features25 across all cell lines. Upper (lower) whisker represents 1.5 times extension from the upper (lower) bound of box, i.e., the third (first) quartile.

Interestingly, miR-548 miRNAs had many distinct seed region features from other miRNAs (Table 11). The significance was even stronger than that of let-7 miRNAs, which are highly abundant and known to be related to various biological processes31 (Table 11; here, I define miRNAs called “let-7” as let-7 miRNAs). Not only significance values, but also the degree of differences between miR-548 miRNAs were generally larger than those observed for let-7 miRNAs (Figure 3). Thus, it is plausible that promoters of genes targeted by miRNA-548 miRNAs are highly hypomethylated, even though miR-548 miRNAs are not highly expressed. However, it is not clear at present how binding affinity is biologically related to miRNA-targeted-gene promoter methylation.

Table 11. Comparison of seed region features between miR-548 or let-7 miRNAs and other miRNAs.

let-7 miRNAs
MeanP-values
Featurelet-7 miRNAsOthersP<P>q-value
1site_nb1.371.463.36e-029.66e-015.19e-02
2seed6_nb0.991.075.69e-029.43e-018.06e-02
3seed7_nb0.380.404.17e-015.83e-014.17e-01
4tgs_au0.610.551.00e+003.45e-04 1.17e-03
5tgs_position419.36395.699.27e-017.27e-029.51e-02
6tgs_pairing3p2.192.292.08e-039.98e-01 5.89e-03
7dg_duplex-12.09-14.509.92e-017.91e-03 1.92e-02
8dg_binding-11.93-13.929.87e-011.26e-02 2.68e-02
9dg_duplex_seed-6.05-7.309.69e-013.14e-025.19e-02
10dg_binding_seed-6.30-7.639.85e-011.46e-02 2.76e-02
11dg_open18.1020.681.83e-041.00e+00 7.78e-04
12dg_total6.416.711.98e-018.02e-012.25e-01
13prob_exact0.170.251.12e-081.00e+00 9.52e-08
14prob_binomial0.210.378.97e-091.00e+00 9.52e-08
15cons_bls0.450.473.37e-016.63e-013.58e-01
16selec_phylop0.570.611.97e-051.00e+00 1.12e-04
17mirmap_score-0.07-0.061.29e-018.71e-011.57e-01
miR-548 miRNAs
MeanP-values
FeaturemiR-548 miRNAsOthersP<P>q-values
1site_nb1.801.451.00e+002.43e-21 4.14e-21
2seed6_nb1.301.061.00e+004.81e-13 6.29e-13
3seed7_nb0.500.391.00e+001.50e-11 1.82e-11
4tgs_au0.650.551.00e+005.52e-30 3.13e-29
5tgs_position334.06398.159.03e-191.00e+00 1.40e-18
6tgs_pairing3p2.362.291.00e+003.69e-04 3.69e-04
7dg_duplex-9.45-14.671.00e+004.00e-31 6.80e-30
8dg_binding-9.22-14.081.00e+001.74e-30 1.48e-29
9dg_duplex_seed-3.78-7.421.00e+001.55e-29 6.59e-29
10dg_binding_seed-4.39-7.741.00e+003.25e-28 1.11e-27
11dg_open16.7520.806.08e-271.00e+00 1.72e-26
12dg_total7.976.661.00e+002.29e-08 2.44e-08
13prob_exact0.310.251.00e+006.35e-22 1.20e-21
14prob_binomial0.480.371.00e+001.61e-15 2.28e-15
15cons_bls0.530.471.00e+004.85e-09 5.49e-09
16selec_phylop0.540.621.15e-221.00e+00 2.44e-22
17mirmap_score-0.01-0.061.00e+001.17e-25 2.84e-25

Features included are those defined in miRmap25. Mean values were computed for miR-548 or let-7 miRNAs and other miRNAs. P< (P>) indicates that mean values in miR-548 or let-7 miRNAs are lower (greater) than those in other miRNAs. P-values between P< and P> in each feature were compared with each other and only the smaller of the two was considered. Selected P-values were adjusted based on the BH criterion23 and were presented as q-values. q-values < 0.05, which were regarded to be significant, are bold faced.

39bce5b3-886e-4f2d-a736-8473fa282275_figure3.gif

Figure 3. Scatter plot of the logarithmic ratio of seed region features for miR-548/let-7 miRNAs compared to all other miRNAs.

Scatter plot of absolute binary logarithmic ratios of seed region features between miR-548 (horizontal) or let-7 (vertical) miRNAs and other miRNAs. Solid line indicates that logarithmic ratios are equal between horizontal and vertical axes.

Treatment of shared miRNA seed sequences

Some miRNAs were known to share seed sequences. These seed-sharing-miRNAs consequently also share sets of target genes. Because P-values were attributed to the set of miRNA target genes, I attributed the same P-value to multiple seed-sharing-miRNAs without merging them. This approach can be rationalized by considering the fact that each seed-sharing-miRNA may still target a given gene, even though seed sequences are shared between them. Furthermore, it is important to note that P-values attributed to a seed sequence cannot be directly compared with miRNA expression as I have done above (see Table 8) and in my previous work15. There are additional reasons for not merging seed-sharing-miRNAs into one, which are discussed in the Appendix.

Conclusions

Figure 4 is a schematic summary of the principal findings of this study. To the best of my knowledge, this is the first time that miRNA-targeted-gene promoter methylation in many types of cell lines have been observed (Table 1 and Table 2). miRNA-targeted-gene hypomethylation was correlated (Table 5 and Table 6) with the regulation of miRNA target genes, which had a reciprocal relationship with miRNA expression; genes with miRNA-targeted-gene hypomethylated promoters were downregulated during cell senescence15,16 and upregulated during differentiation (Table 8). I also found that genes with miRNA-targeted-gene hypomethylated promoters were specific targets of miR-548 miRNAs (Table 9 and Table 10). The regulation of target genes by miR-548 miRNAs, which were typically expressed at low levels, is likely the result of miRNA-targeted-gene promoter hypomethylation. However, all findings here were correlative and thus lacked the consideration of an underlying biological mechanism. However, some seed region features were suggested to have significant correlations with miRNA-targeted-gene promoter methylation (Figure 2 and Figure 3 and Table 11).

39bce5b3-886e-4f2d-a736-8473fa282275_figure4.gif

Figure 4. Schematic summarizing the principal findings of the present study, as well as those reported previously.

Promoters are methylated dependent upon miRNAs that target each gene (Table 1 and Table 2). Genes with miRNA-targeted-gene hypomethylated promoters are downregulated during cell senescence and upregulated during cellular differentiation (Table 5 and Table 6), whereas miRNAs that target genes with hypomethylated promoters were upregulated during cell senescence15,16 and downregulated during differentiation (Table 8). miR-548 miRNAs, which are expressed at low levels, were suggested to regulate target genes with the assistance of miRNA-targeted-gene promoter hypomethylation (Table 9 and Table 10). Color allows just below “MiRaGE” box indicated statistical tests: t-test (blue), Wilcoxon rank sum test (red), and Kolmogorov-Smirnov test (green).

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 23 Jan 2013
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Taguchi Yh. Correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes [version 3; peer review: 1 approved, 2 not approved]. F1000Research 2013, 2:21 (https://doi.org/10.12688/f1000research.2-21.v3)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 2
VERSION 2
PUBLISHED 27 Mar 2013
Views
49
Cite
Reviewer Report 18 Jun 2013
Pål Sætrom, Department of Cancer Research and Molecular Medicine, Norwegian University of Science and Technology, Trondheim, Norway 
Not Approved
VIEWS 49
Although this article’s topic is interesting, I found that its contents were confusing and inaccessible. I also have some reservations regarding the methods. 
  • The abstract is misleading.
The author states in the abstract: “It was found that promoter methylation was affected by ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Sætrom P. Reviewer Report For: Correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes [version 3; peer review: 1 approved, 2 not approved]. F1000Research 2013, 2:21 (https://doi.org/10.5256/f1000research.1268.r1002)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 09 Sep 2013
    Y-h. Taguchi, Department of Physics, Chuo University, Tokyo, Japan
    09 Sep 2013
    Author Response
    Thank you very much for your efforts to review my article. In the following, I would like to reply to your valuable comments point by point.

    Although this article’s topic is interesting, ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 09 Sep 2013
    Y-h. Taguchi, Department of Physics, Chuo University, Tokyo, Japan
    09 Sep 2013
    Author Response
    Thank you very much for your efforts to review my article. In the following, I would like to reply to your valuable comments point by point.

    Although this article’s topic is interesting, ... Continue reading
Views
53
Cite
Reviewer Report 14 May 2013
Denis Dupuy, Institut Europeen de Chimie et Biologie, Pessac, France 
Not Approved
VIEWS 53
I had initially understood that the title of this study could be: "Correlation between promoter methylation and miRNA regulation of target genes." Which in simpler terms could mean that miRNA targeted genes tend to have more hypomethylated promoters (or hypermethylated maybe) ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Dupuy D. Reviewer Report For: Correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes [version 3; peer review: 1 approved, 2 not approved]. F1000Research 2013, 2:21 (https://doi.org/10.5256/f1000research.1268.r949)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 16 May 2013
    Y-h. Taguchi, Department of Physics, Chuo University, Tokyo, Japan
    16 May 2013
    Author Response
    Dear Denis, 


    Thank you very much for your efforts to review my paper. Here I respond to you point by point. If these responses resolve disagreements between us and then if ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 16 May 2013
    Y-h. Taguchi, Department of Physics, Chuo University, Tokyo, Japan
    16 May 2013
    Author Response
    Dear Denis, 


    Thank you very much for your efforts to review my paper. Here I respond to you point by point. If these responses resolve disagreements between us and then if ... Continue reading
Version 1
VERSION 1
PUBLISHED 23 Jan 2013
Views
30
Cite
Reviewer Report 06 Feb 2013
Denis Dupuy, Institut Europeen de Chimie et Biologie, Pessac, France 
Approved with Reservations
VIEWS 30
I find that the often-used term of ‘miRNA-targeting-specific promoter methylation’ (sic) should be replaced by a less confusing expression, e.g. ‘miRNA-targeted-gene promoter methylation’, as ‘targeting’ and ‘specific’ are redundant and confuse the point. Specifically, the grammatical construction would seem to
... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Dupuy D. Reviewer Report For: Correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes [version 3; peer review: 1 approved, 2 not approved]. F1000Research 2013, 2:21 (https://doi.org/10.5256/f1000research.1183.r751)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
33
Cite
Reviewer Report 04 Feb 2013
Wenqian Hu, Center for RNA Molecular Biology, Whitehead Institute for Biomedical Research, Cambridge, MA, USA 
Approved
VIEWS 33
This paper reports interesting correlations between the regulation of genes by miRNAs and the methylation status of the miRNA target genes. Although the two processes occur in two different compartments of the cell, these correlations suggest that there may be ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Hu W. Reviewer Report For: Correlation between miRNA-targeted-gene promoter methylation and miRNA regulation of target genes [version 3; peer review: 1 approved, 2 not approved]. F1000Research 2013, 2:21 (https://doi.org/10.5256/f1000research.1183.r748)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 23 Jan 2013
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.