Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Global population-specific variation in miRNA associated with cancer risk and clinical biomarkers

BMC Medical Genomics volume 7, Article number: 53 (2014) Cite this article

Abstract

Background

MiRNA expression profiling is being actively investigated as a clinical biomarker and diagnostic tool to detect multiple cancer types and stages as well as other complex diseases. Initial investigations, however, have not comprehensively taken into account genetic variability affecting miRNA expression and/or function in populations of different ethnic backgrounds. Therefore, more complete surveys of miRNA genetic variability are needed to assess global patterns of miRNA variation within and between diverse human populations and their effect on clinically relevant miRNA genes.

Methods

Genetic variation in 1524 miRNA genes was examined using whole genome sequencing (60x coverage) in a panel of 69 unrelated individuals from 14 global populations, including European, Asian and African populations.

Results

We identified 33 previously undescribed miRNA variants, and 31 miRNA containing variants that are globally population-differentiated in frequency between African and non-African populations (PD-miRNA). The top 1% of PD-miRNA were significantly enriched for regulation of genes involved in glucose/insulin metabolism and cell division (p < 10−7), most significantly the mitosis pathway, which is strongly linked to cancer onset. Overall, we identify 7 PD-miRNAs that are currently implicated as cancer biomarkers or diagnostics: hsa-mir-202, hsa-mir-423, hsa-mir-196a-2, hsa-mir-520h, hsa-mir-647, hsa-mir-943, and hsa-mir-1908. Notably, hsa-mir-202, a potential breast cancer biomarker, was found to show significantly high allele frequency differentiation at SNP rs12355840, which is known to affect miRNA expression levels in vivo and subsequently breast cancer mortality.

Conclusion

MiRNA expression profiles represent a promising new category of disease biomarkers. However, population specific genetic variation can affect the prevalence and baseline expression of these miRNAs in diverse populations. Consequently, miRNA genetic and expression level variation among ethnic groups may be contributing in part to health disparities observed in multiple forms of cancer, specifically breast cancer, and will be an essential consideration when assessing the utility of miRNA biomarkers for the clinic.

Peer Review reports

Background

MicroRNA (miRNA) expression profiles have been demonstrated to be unique for a wide range of human diseases, including different stages of tumor progression and metastasis [1]. MiRNA expression levels and function can also be affected by global factors, such as genomic variation due to population history, which have been less well studied. MiRNAs function mainly to inhibit protein synthesis through binding between miRNA seed sequences and complementary sequences on target messenger RNA (mRNA) genes. This binding causes degradation and/or translational repression of mRNA genes [2–6]. A single mature miRNA (21–25 base pairs) has the ability to inhibit protein synthesis of over 6,000 mRNA targets [7–13], and miRNAs are predicted to regulate the protein expression of 30–60% of all human protein-coding genes [14, 15]. Therefore, changes in miRNA expression in response to disease, the combined effects of circulating miRNA being extremely stable in blood and serum [16] and advances in miRNA detection methods, such as in situ hybridization and RT-PCR, have made miRNAs excellent candidates as diagnostic and prognostic markers in the clinic [17, 18]. As a result, miRNAs are currently under clinical investigation as biomarkers for a number of complex diseases, including breast cancer, diabetes (types 1 and 2), asthma, sepsis, lung cancer, prostate cancer, leukemia (ALL and AML), and various pediatric cancers [1, 19–27].

Most recent studies identifying potential miRNA biomarkers of disease have been performed in European or Asian populations [28–32] with only a handful of studies performed in populations of African descent [33, 34]. Nonetheless, data from these studies have demonstrated that circulating miRNA profiles were considerably different between African-Americans and European-Americans in early stage lung cancer [33] and were expressed differentially between these populations in early stage breast cancer [34]. Given the importance and ubiquitous nature of miRNA-mediated gene expression, it has been proposed that SNPs mapping within miRNA, particularly within the miRNA seed sequences base positions 2–8 of the mature miRNA, may have functional consequences resulting in expression and/or phenotypic variation [4] (Figure 1B). Therefore, genomic variation within miRNA, due to human population history, may be affecting ethnic disparities in complex diseases such as lung and breast cancer through two mechanisms of action: (1) by affecting miRNA expression patterns and (2) by disrupting miRNA/mRNA target recognition through interfering with seed sequence binding.

Figure 1
figure 1

Variants per miRNA region. (A) The different colored circles on the world map indicate the geographic location of the 14 sampled populations. Populations from Africa are shown in red (Pygmy, Hadza, Sandawe, Yoruba, Luhya, and Maasai), Asian populations in green (Gujarati Indians, Japanese, and Han Chinese), European populations in blue (Toscani Italians, Utah residents with Northern European ancestry), and recently admixed populations in stripes (African-Americans, Mexican-Americans and Puerto Ricans) (See Methods). Representative regions are enclosed by the thick black lines (B) The top panel shows the human miRNA hairpin, hsa-mir-302d, displayed in its modified form with a novel SNP found in its seed sequence circled in black. Seed sequences are defined as the 2nd -8th base of the 5’ end of mature double-stranded miRNA. The mature sequence is highlighted with the guide strand shown in red and the passenger strand shown in blue. The bottom panel depicts the summary of allele types in miRNA genes including novel and non-singleton alleles. Here, mature sequences (green) are defined as all bases within the mature miRNA excluding those within the seed sequence (red). Consistently, stem loops (blue) are defined as all bases excluded from the mature and seed sequences but included in the pre-miRNA hairpin (grey).

Full size image

To date, the genetic coverage of miRNA genes has been low (3× to 5× coverage) in large-scale resequencing projects, such as the 1000 Genomes, which can be problematic for identifying low-frequency variants with high confidence [35]. This is particularly true in African populations which have been shown to possess higher levels of genetic diversity, including low-frequency polymorphisms, compared to other populations worldwide [36–39]. However, diverse African populations are still highly underrepresented in studies of genomic variation [38]. Furthermore, recent studies of population differentiation at miRNA variants have used microarray technology, which captures only a fraction of the total number of miRNA in the genome, resulting in ascertainment bias and also reduced power to discover novel variants [3, 40, 41].

In the present study, we analyzed 1524 miRNA sequences at high coverage (60×) using whole genome sequence data from 69 individuals representing 14 worldwide populations from Europe, Asia, the Americas and Africa, (Figure 1A) including 3 African hunter-gatherer populations not included in the 1000 Genomes datasets [39]. These samples allow direct comparison of genetic variability in all annotated miRNAs without the ascertainment bias common to SNP array data in geographically and ethnically diverse populations. Based on our data, we identified 33 previously unidentified variants in miRNA genes and 31 significantly population-differentiated (PD) variants based on estimates of FST between African and non-African populations. We identified 7 PD variants within miRNA that have been experimentally linked to onset, progression, and/or metastasis of cancers with known health disparities between patients of European and African descent. Specifically, we find a T-allele at SNP rs12355840 in hsa-mir-202, that has been shown to increase miRNA expression in vivo and to be protective against breast cancer mortality [42], and to be highly PD between African and non-African populations. To our knowledge, a complete survey of genetic variation in all miRNA using high-coverage whole genome data has not previously been performed and has uncovered novel miRNA variants, and determined miRNA biomarker candidates that may differ among diverse population groups.

Results

Novel variants identified in miRNA hairpins

In a sample of 69 unrelated individuals, we identified a total of 773 polymorphisms (700 SNPs and 73 insertion-deletions) in pre-miRNA hairpins which passed strict quality control filters (Figure 1B). Of these 773 variants, 411 mutations occurred in pre-miRNA stem/loop regions, 242 in mature miRNA, and 120 in the seed sequences (Figure 1B and Additional file 1: Table S1). Among these polymorphisms, we identified 198 previously undescribed mutations that are currently not present in dbSNP v135. The number of alleles per base was calculated separately for each region of the miRNA (stem-loop, mature, and seed). In our dataset, allele frequencies were slightly higher in the stem-loop region compared to the mature miRNA and seed region (0.013, 0.011 and 0.011, respectively).

To control for somatic mutations, undescribed mutations found in a single individual (ie. singletons) were removed from analyses and 33 novel variants present in multiple individuals (non-singletons) were analyzed (Figure 1B). Among the non-singleton mutations identified in pre-miRNA hairpins, 5 novel mutations were in highly conserved miRNA seed sequences (Figure 1B and Additional file 1: Table S2). The first novel seed variant, a “C/T” SNP at chr3 128081086, was located in the 3’ strand of human miRNA hsa-mir-1280 and was present in one Hadza and one Sandawe, two hunter-gatherer populations from Tanzania. The second, a “C/T” SNP at chr1 62544469, was located in the 3’ strand of hsa-mir-942 and was present in one Hadza and two Sandawe individuals. Two novel seed sequence indels found at low frequency were an “ACA” deletion in miRNA hsa-mir-4483 at chr10 115537763-115537766, found in 2 Yoruban individuals, and a “T”-allele insertion in hsa-mir-3940 at chr19 6416444 in 2 individuals with northern European ancestry from Utah. Lastly, a “CT” deletion located on chromosome 6 in the 3’ strand of hsa-mir-4640 was found in 9 individuals from 7 global populations (namely, 1 Pygmy, 2 Sandawe, 1 Yoruba, 1 Maasai, 2 African-Americans, 1 Mexican-American and 1 Gujarati Indian individual). Using miRNA target prediction software, we determined that in the absence of the “CT” mutation in miRNA hsa-mir-4640, the 3’ strand was predicted to individually target 316 binding sites in 79 genes (Additional file 1: Table S3). With the “CT” deletion, however, the number of predicted targets dropped to 11 binding sites in 3 genes, where gene targets did not overlap between the original and modified miRNA.

MicroRNA conservation and frequency of population-specific miRNA alleles

We compared levels of nucleotide diversity in miRNA to other genomic sequence classes (for examples, exons and introns) in each population group using Watterson’s estimator of theta (θW) (Figure 2A). We found in our high-coverage sequence data that miRNAs were among the most conserved sequences in the genome, at the same level as exons (Figure 2A), consistent with results from prior lower-coverage sequencing studies [41]. Overall, however, African populations had the highest level of nucleotide diversity across all sequence classes (θW = 5.0 ± 0.7 × 104) compared to European and Asian populations (θW = 3.6 ± 0.2 × 104) (Figure 2A).

Figure 2
figure 2

Population specific SNPs. (A) Genomic nucleotide diversity, as estimated by Watterson’s (θ), for each population. Within populations, nucleotide diversity was measured for different genomic sequence classes including miRNA genes (black) which were conserved to a level comparable to exons (blue), in all populations. (B) The number of population specific miRNA alleles (PSMA) divided by the number of genomes collected for each population (shown in black). The right axis displays the number of high frequency PSMAs, with an allele frequency greater than 50% in a single population (red).

Full size image

In addition, we identified 319 population-specific miRNA alleles (PSMAs), defined as variants present exclusively in 1 of the 14 globally sampled populations. About two-thirds (66.8%) of PSMAs were present in African populations and the proportion of population-specific alleles (ie. population-specific density) was highest in Africans relative to non-Africans, consistent with prior analyses of human genetic variation [41]. In particular, Pygmy hunter-gatherers had the highest population-specific density with 11 PSMAs per genome (Figure 2B). Among the Hadza hunter-gatherers, we found two high frequency PSMAs with allele frequency ≥ 50% in the Hadza only, a novel “A/C” SNP in the stem-loop of hsa-mir-1291 and an “A/G” SNP (rs111566161) in the 3’ mature sequence of miRNA hsa-mir-4711. This “A/G” SNP was found to be exclusively shared by two Hadza and a Southern Kalahari San individual, a population thought to share an ancient genetic ancestry with the Hadza and other African click speaking populations [43, 44].

Population differentiation of human miRNA

To measure population differentiation, we calculated pairwise FST at variants in miRNA genes (See Methods). FST ranges from 0 to 1 with an estimate of 0 indicating no population differentiation and an estimate of 1 indicating complete differentiation. Estimates of FST at miRNA variants were calculated hierarchically: (a) between individual populations (Figure 3A-B), (b) between major geographic regions (for example, Europe-Africa, Europe-Asia, Africa-Asia) (Additional file 2: Figure S1) and (c) between pooled African and pooled non-African populations (Figure 3C). We then compared the individual estimates of FST to empirical distributions of miRNA FST. The pairwise FST values that were outliers (ie. within the top 5% of the empirical distribution or above the 95th percentile) were classified as population-differentiated (PD) and variants with these extreme values were inferred to be enriched for SNPs under recent selection [45]. Our data showed that among major geographic regions, African populations had the highest average FST values among populations, with 34 PD-miRNA alleles between Africa and Asia, 33 PD-miRNA alleles between Europe and Africa, and 18 PD-miRNA alleles between Europe and Asia (Figure 3B; Additional file 2: Figure S1).

Figure 3
figure 3

Pairwise F ST for each miRNA allele. (A) Boxplots of pairwise FST between African populations (green), between non-African populations (red), and between African and non-African populations (blue) are grouped hierarchically by major geographic region: Europe-Asia, Africa-Europe, Africa-Asia, and between African and non-African samples (grey). Black circles are outliers above 1.5 times the interquartile range. The red circles and the black solid lines are the mean and median of the distributions, respectively. (B) Zoom in of boxplots in A showing the median (black lines) and average FST of each pairwise group (red circles). (C) Pairwise FST for all miRNA variants between African (n = 32) and non-African (n = 25) samples. Population differentiated miRNAs (PD-miRNAs) contain variants above the 95th percentile (FST > 0.186) shown as a solid blue line. Highly population differentiated miRNAs (HPD-miRNAs) contain variants above the 99th percentile (FST > 0.366) shown as a dotted red line and are annotated with miRNA affiliation and rs numbers.

Full size image

In a comparison of pooled African and pooled non-African populations, we identified 31 PD-miRNAs containing variants with outlier FST values above the 95th percentile of the distribution of miRNA FST values (FST ≥ 0.186 with p < 0.05) (Figure 3C and S2; Additional file 1: Table S4), with 4 PD-miRNA variants in seed sequences (Additional file 2: Figure S3). Furthermore, when we apply a more stringent criteria for population-differentiation among miRNA variants, we found 8 highly population-differentiated miRNAs (HPD-miRNAs) with variants above the 99th percentile representing the top 1% of all pairwise FST estimates (FST ≥ 0.366; p < 0.05) (Figure 3C and Additional file 1: Table S4).

Messenger RNA target and functional enrichment of HPD-miRNAs

Experimentally-validated mRNA targets of the 8 HPD-miRNAs were identified by querying 45 publically available deep sequencing and microarray datasets (See Methods). We found that these 8 HPD-miRNAs experimentally down-regulated the expression of 2,139 unique human mRNA target genes in at least two datasets (Figure 4A). Target gene enrichment analysis revealed that 72 of the 2,139 mRNA targets were significantly over-regulated (p < 0.05) by these 8 HPD-miRNAs compared to a null or random set of 8 miRNA (See Methods) (Figure 4A-B and Additional file 1: Table S5). The genes of these over-regulated mRNA targets were involved in immune response, metabolism, developmental processes, cell communication, transport and response to stress (Additional file 2: Figure S4). Additionally, 14 of the 72 enriched gene targets were reported to be candidate loci in genome-wide association studies (GWAS) (Additional file 1: Table S6).

Figure 4
figure 4

Analysis flowchart. (A) Blue boxes represent allele counts, with 773 alleles found in miRNA genes and 31 population-differentiated alleles within miRNA (PD-miRNA). Red boxes represent messenger RNA targets regulated by the 8 highly population-differentiated miRNAs (HPD-miRNAs). These targets were found to be involved in over 8000 gene ontology functions (GO terms) (orange boxes) and over 5,000 biological processes (green boxes). Seven PD-miRNAs were also found to be differentially expressed in 37 diseases, with cancer being the top hit. (B) Description of Target enrichment, Genomic Functional Enrichment, and MiRNA functional enrichment that was performed and displayed in A. Output of statistical function enrichments were clustered to determine overrepresented functions of HPD-miRNAs. Sugar and insulin metabolism were among the top hits in genomic functional enrichment (p < 0.001), including regulation of the glucose/insulin response pathway. Mitosis was the top hit in miRNA functional enrichment (p <0.001), including regulation of mitosis transition and checkpoint genes by HPD-miRNA.

Full size image

Functional enrichment was performed by two methods for the 2,139 mRNA gene targets of HPD-miRNA. First, genomic functional enrichment was used to determine if particular biological functions regulated by the 8 HPD-miRNAs were overrepresented compared to the null expectation across the genome (Figure 4B) [46]. Based on this analysis, we identified 319 statistically over-represented functions (p < 0.05; n = 319) (Figure 4). The top hits included sugar (mannose, fructose, and glucose) metabolism and the regulation of insulin (p < 10−7) (Additional file 1: Table S7). Secondly, we performed a bootstrapping analysis of miRNA functional enrichment to account for the non-random distribution of miRNA targets throughout the genome and identified biological pathways overrepresented by HPD-miRNAs as compared to other non population-differentiated miRNA sets (Figure 4B). In particular, the 2,139 experimentally-validated gene targets of the 8 HPD-miRNAs were found to function in 5,475 annotated biological processes based on the ENSEMBL gene ontology. Each of the 5,475 functions were assigned a probability of being regulated by a random null set of 8 miRNAs (See Methods) and this probability was compared to the observed regulation by the 8 HPD-miRNAs. We found 55 biological processes that were significantly enriched for regulation by our set of 8 HPD-miRNAs (p < 0.05) (Additional file 1: Table S8), and identified 120 HPD-miRNA gene targets involved in enriched biological processes; of particular interest, these targets include DICER1 and BRCA1 among other genes (Additional file 1: Table S8).

In addition, clustering analysis of the 55 significant biological processes, based on gene ontology similarity, revealed nine major pathways affected specifically by the 8 HPD-miRNAs: the mitosis pathway, axons and morphogenesis, response to toxins, signaling, acid regulation, transcription, ion sequestration, muscle movement and immune response (Additional file 2: Figure S5). The mitosis pathway contained the highest number of significantly over-represented processes, specifically mitotic transitions (G1 and metaphase/anaphase transition), cell cycle checkpoints, and cytokinesis, representing over 20% of the 55 biological processes. Interestingly, the altered expression of mitotic genes, which can occur through miRNA regulation, is one of the known characteristics of the onset and progression of cancer.

Disease association analysis reveals links to multiple cancers

Disease associations were also identified for all 31 PD-miRNAs through the miRNA disease database MiRGator [47] (Figure 5A). Of the 31 PD-miRNAs, 7 miRNAs (hsa-mir-202, hsa-mir-196a-2, hsa-mir-423, hsa-mir-943, hsa-mir-520h, hsa-mir-1908 and hsa-mir-647) were identified in prior studies as being differentially expressed in human diseases representing a significant increase above expectation (1.1 out of 31) for disease-associated miRNAs (95% C.I. = 1.08 – 1.11; p = 2.2 × 10−16) (Figure 5B). The 7 PD-miRNAs were associated with 37 diseases, most extensively cancer, with all 7 PD-miRNAs being implicated in cancer risk or as biomarkers for one or multiple cancers, specifically, breast cancer, prostate cancer, testicular cancer, ovarian cancer, lung cancer, renal cell carcinoma, gastric cancer, pancreatic cancer, head and neck cancer, colon cancer, esophageal cancer, brain tumors, thyroid cancer, glioblastoma, glioma, and germ cell tumors (Figure 5C). The hsa-mir-196a-2 T-allele at SNP rs11614913 has been significantly associated with increased risk for esophageal cancer in non-smoking European males [48] but decreased risk for breast, lung and gastric cancers in Chinese populations [48]. In the present study, we observed a significantly lower frequency of the hsa-mir-196a-2 T-allele at SNP rs11614913 in Africans compared to non-African populations (FST = 0.41; p < 0.001) (Figure 5D).

Figure 5
figure 5

Disease network analysis. (A) The number of miRNA disease associations plotted against the pairwise FST between Africans and non-Africans for each miRNA allele. The 99th percentile (FST > 0.366) is shown with a dotted red line and the 95th percentile (FST > 0.186) is shown with a solid blue line. (B) Table of FST values for the 7 PD-miRNAs associated with human diseases and containing alleles with FST > 95th percentile and p < 0.05. MRNA targets were obtained using 45 public deep sequencing and microarray datasets (See Methods) (C) Disease network graph depicting the interactions between the 7 PD-miRNA (shown in black) and human disease, where up-regulation is shown in red, down-regulation (blue), unknown regulation (grey), and miRNA with reports of both up and down regulation in the same disease state shown in green. Labels are scaled by number of disease associations where miRNA or diseases with larger numbers of connecting lines are indicated by larger dot sizes within the network. (D) Population allele frequency differentiation of the breast, lung, gastric and esophageal cancer associated T-allele at SNP rs11614913 within miRNA hsa-mir-196a-2. African populations have decreased frequency of the T allele (green) as compared to non-African populations in blue (p < 0.001).

Full size image

Notably, one of our PD-miRNAs, hsa-mir-202, contains a T-allele with known effect on miRNA expression and a protective effect on breast cancer mortality [42]. In our dataset, we found that the frequency of the T-allele at hsa-mir-202 was lower in African populations (26%) compared to non-African populations (65%), on average (Figure 6A and Additional file 1: Table S9). We observed similar frequencies of the T-allele in African and non-African populations in the 1000 Genomes Project data, that sampled only one African population (YRI) (16%) and 2 non-African populations (CEU, CHB + JPT) (83.3% and 92%, respectively) (See Methods). In addition, we also observed variability in the frequency of the T-allele within Africa; specifically, the T-allele occurred at lower frequency in the non-hunter-gatherer Yoruba, Maasai, and Luhya populations ((0.08%) compared to hunter-gather Pygmy, Hadza, and Sandawe populations (43%), on average (Figure 6A and Additional file 1: Table S9). We also found the T-allele at moderate frequency in African-American populations (Figure 6B).

Figure 6
figure 6

The role of miRNA hsa-mir-202 SNP in gene expression and breast cancer. (A) Differences in frequency of the T-allele at mir-202 across populations. African populations are highlighted in green and non-African populations in blue. (B) Average frequency of the T-allele at mir-202 in Africa, Europe and Asia (double circles) are compared to a limited sample of US population groups.

Full size image

Discussion

Population differentiation and functional enrichment of miRNA

Based on analysis of high-coverage whole genome sequence data in 69 individuals from 14 worldwide populations, we identified 33 novel polymorphisms in miRNA genes and 31 variants with high levels of genetic population differentiation between Africans and non-Africans (PD-variants). Five novel variants (2 SNPs and 3 indels) were located in miRNA seed sequences where they could have a large effect on the number, strength, and specificity of each miRNA/mRNA target interaction. Specifically, a novel “CT” deletion found in 7 of the 14 globally diverse populations altered the predicted mRNA targets of hsa-mir-4640 to include 3 additional targets and removed all of its 79 original predicted targets from regulation, indicative of the disruptive potential of seed sequencing indels on miRNA function. In addition, among the 31 PD-variants we found 8 HPD-miRNAs with variants that lie in the top 1% of pairwise estimates of FST, indicating extensive population differentiation at these loci between African and non-African populations consistent with a model of local adaptation (Figure 4) [5, 45, 49].

From functional enrichment analysis we observe that HPD-miRNAs are significantly enriched for regulation of genes involved in the glucose/insulin metabolism pathway (p < 10−7), and cellular division (p < 0.001), notably the regulation of genes involved in mitosis (specifically, mitotic checkpoints, transitions, and cytokinesis) (Figures 4 and 5; Additional file 1: Tables S7 and S8). Aberrant cell division during mitosis or aberrant gene expression levels during chromosome segregation often result in chromosomal instability, which is a key diagnostic feature of most cancers [50]. The disruption of cell division by altering gene expression levels in response to mitotic instability has been strongly correlated with tumor development and progression; both in vitro and in vivo evidence have demonstrated that in the absence of other cell cycle and DNA repair defects, mitotic disruption can transform cells and predispose them toward cancer [50]. Given that mis-regulation of cellular division is the hallmark of cancer, it is striking that miRNAs with highly population-differentiated alleles are observed to be significantly enriched for regulation of mitotic pathway genes including genes such as S100A8 and P2RX3 whose expression profiles are currently used as biomarkers for multiple cancers.

The Role of miRNA in ethnic disparities in cancer susceptibility

In addition, we identified 7 population-differentiated miRNAs where expression level differences of these PD-miRNAs have been correlated with cancer and other disease phenotypes (Figure 5B-C). Among the identified cancers, higher mortality rates have been reported for breast, ovarian, gastric, prostate and testicular-germ cell cancers in individuals of recent African ancestry compared to individuals of either European or Asian descent (p < 0.001) [51–54]. Of particular interest is hsa-mir-202 which contained one of the most highly population-differentiated variants in our dataset and is one of two miRNAs currently under investigation as a circulating blood-based marker for the detection of non-Hodgkin lymphoma and early stage breast cancer [42, 55]. Recent in vitro functional data demonstrated that the T-allele was protective against breast cancer mortality by first increasing mature hsa-mir-202 expression levels, leading to subsequent down-regulation of its gene targets, including cancer related genes CRYBB2, DICER1, SART1, S100A8, P2RX3, and BRCA1[42]. Diminished expression of mature hsa-mir-202 in individuals harboring at least one non T-allele resulted in a significantly elevated risk of non-Hodgkin lymphoma (OR = 1.83, 95% CI: 1.17–2.85; P = 0.008) [42]. Our data showed that African and African American populations had a lower frequency of the T-allele compared to European and Asian populations, suggesting decreased baseline expression levels of mature hsa-mir-202 in African populations. In the context of tissue specific expression, hsa-mir-202 expression varies considerably by tissue type and is more highly expressed in prostate cells compared to breast or ovarian tissues (Additional file 2: Figure S7). Therefore, lower frequencies of the T-allele in African populations have the potential to reduce gene expression levels down below critical thresholds in the breast and ovarian tissues. This reduction may contribute to population differences seen in the onset of breast and ovarian cancer in women of African and European descent (Additional file 2: Figure S7).

Similarly, the other 6 PD-miRNAs identified (hsa-mir-423, hsa-mir-196a-2, hsa-mir-520h, hsa-mir-1908, hsa-mir-647 and hsa-mir-943) have also been implicated in cancer susceptibility in human populations [56]. First, the hsa-mir-423 SNP rs6505162 has been shown to confer reduced risk for breast cancer in women of European decent in GWAS [57]. Second, the hsa-mir-196a-2 SNP rs11614913 CC genotype has been significantly associated with increased risk for breast, lung and gastric cancers in Chinese populations; conversely, the homozygous TT genotype has been significantly associated with esophageal cancer in non-smoking European males [48]. Based on these results, recent studies have called for more detailed analysis of the frequency of this allele in different ethnic groups [58]. We observed a significantly higher frequency of the hsa-mir-196a-2 C-allele at SNP rs11614913 in African compared to non-African populations (FST = 0.41; p < 0.001) (Figure 5D). Third, hsa-mir-520h expression was determined to be significantly associated with E1A-mediated tumor suppression and cell migration during cancer metastasis and inhibition of hsa-mir-520h significantly decreased the downstream ability of cancer cells to migrate and invade other areas of the body [59]. This pattern has also been observed consistently in different types of cancer including pancreatic, breast and ovarian cancer [59–61]. Also, up-regulation of hsa-mir-520h was shown to increase the effects of the anticancer drug resveratrol in slowing lung cancer tumor mobility [62]. Finally, multiple studies have linked miRNAs hsa-mir-1908, hsa-mir-647 and hsa-mir-943 expression to various cancers known to have ethnic specific disparities [63–65]. Overall, these studies demonstrate that genetic variability within miRNA has the potential to vary miRNA expression and/or mRNA target binding which can be strongly correlated with the onset of multiple cancers, the progression of cancer metastasis and the response to drug therapies. We demonstrate that the frequency of clinically important genomic miRNA variants varies significantly among ethnic populations, particularly between African and non-African groups. Thus, we suggest that population-differentiated variation in miRNA may contribute to ethnic disparities seen in certain forms of cancer.

Conclusions

Here, we identified several miRNA genetic variants that are highly differentiated among human populations and uncovered a set of HPD-miRNAs that play a role in the suppression, susceptibility, and metastasis of cancer cells. We also found that some of the HPD-miRNA variants are in regions of strong linkage disequilibrium (D' = 1) with markers in a commonly used genotyping array (Illumina 1 M Duo) in African populations that could be included in genome-wide association studies of disease (Additional file 2: Figure S6). Finally, although we focused on population-differentiated miRNAs known to be associated with disease, we also identified an additional 24 PD-miRNAs that represent interesting candidate loci for further study of differential disease risk in ethnically diverse populations. Further investigation is needed in order to understand the patterns of variation at miRNA and their role in phenotypic variation and human adaptation, particularly in African populations which are greatly underrepresented in genomic studies. Additional RNA-sequencing studies, together with eQTL mapping, will be needed in order to assess the effect of PD-variants on gene expression in global populations. Furthermore, future follow-up studies could integrate SNPs within downstream miRNA target sites [66] and upstream miRNA-regulomes (i.e. transcription factors that regulate miRNA genes) [67] with our findings to examine population differentiation or disease association in all phases of the miRNA cycle.

Methods

Whole genome sequencing and sample collection

High quality whole genome sequencing (~60× coverage) was obtained for 69 globally diverse individuals from publically-available datasets. Fifteen African hunter-gathers were obtained from Lachance et al. 2012 [39], including 5 Pygmy (Py) (three Baka, one Bakola, and one Bedzan), 5 Hadza (Hz) (plus two technical replicates), and 5 Sandawe (Sw) using the Complete Genomics sequencing platform [68, 69].

Additionally, 54 unrelated individuals were obtained directly from Complete Genomics including 9 individuals of Northern European ancestry (combined as CEPH - including 4 CEPH and 5 CEU individuals), 9 individuals of Yoruban ancestry (YRI), 5 individuals of Mexican ancestry (MEX), 5 African-Americans living in Dallas (ASW), and 2 individuals of Puerto Rican ancestry (PUR), and 4 individuals each of Toscani Italians (TSI), Japanese ancestry (JPT), Han Chinese (CHB), Gujarati Indian (GIH), Maasai Kenyan ancestry (MKK), Luhya Kenyan ancestry (LWK) [70]. Variants were defined as autosomal alleles that differ from the human reference genome build (GRCh37/hg19), and novel variants are defined as variants that are absent from dbSNP (db135). Complete Genomic Data is available through their public site: http://www.completegenomics.com/public-data/69-Genomes/.

Nucleotide diversity and population differentiation

Pre-miRNA locations and mature miRNA sequences were downloaded from the database mirbase (v18) and included 1524 known pre-miRNAs. We filtered known pre-miRNAs into high and low confidence miRNA, based on experimental validation of their function as target suppressors [71]. High confidence miRNAs were classified as those with at least one experimentally-validated mRNA target, using 45 public datasets (See Messenger RNA Target Enrichment and Analysis). Among novel variants, 165 mutations were identified within high confidence miRNA, 83 were found in miRNA stem-loops, 82 in mature miRNA, and 28 in seed sequences. Seed sequence locations were computed in R software and defined as the 2nd – 8th base pair of the mature miRNA [72]. Novel miRNA target prediction was done in DIANA software [73] for hsa-mir-4640 “CT” deletion. Predicted targets were also filtered for those with experimental validation, based on 45 public datasets (See Messenger RNA Target Enrichment and Analysis). Nucleotide diversity for each sequence class was measured using Watterson’s estimator of θ (θW) [74].

Whole genome sequence data were annotated using the UCSC genome browser for sequence class determination. Pairwise population FST values were calculated using Weir and Cockerham’s weighted equations adjusting for small sample size using the R statistical software [75, 76]. MiRNAs were determined to be population differentiated (PD-miRNA) if they contained an allele with a FST value in the top 95th percentile of the empirical distribution of FST values (FST ≥ top 5% with p < 0.05). Highly population differentiated miRNA (HPD-miRNA) were defined as miRNA containing variants with FST values above the 99th percentile between African and non-Africans (FST ≥ top 1% with p < 0.05). P-values for FST’s were determined by testing the allele frequency difference at each allele between African and non-African samples using a Welsh two-sided t-test (Additional file 2: Figure S2).

Messenger RNA target enrichment and analysis

For HPD-miRNA, messenger RNA (mRNA) targets were identified from the consensus of 45 paired miRNA/mRNA experimental datasets, and subject to the following quality control filters: (1) the mRNA targets have to be correlated with miRNA expression with a correlation coefficient of r < −0.5 and (2) that this level of correlation was observed in at least 2 publicly available experimental deep sequencing or microarray datasets. Experimental datasets consisted of: 6 ENCODE (Encyclopedia of DNA elements) cell line RNA sequencing datasets (GM1287- a lymphoblastoid cell line produced from the blood of a female donor with northern and western European ancestry, H1_hESC- a human embryonic stem cell line, Hela_S3- an immortalized cell line from an African-American female patient with cervical cancer, K562- an immortalized cell line from a female patient with chronic myelogenous leukemia (CML), HepG2- a cell line produced from a male patient with liver carcinoma, and NHEK- a epidermal keratinocyte cell line), 13 cancer deep sequencing datasets from the Cancer Genome Atlas (TCGA) (BLCA - Bladder Urothelial Carcinoma, BRCA - Breast invasive carcinoma, COAD - Colon adenocarcinoma, HNSC - Head and Neck squamous cell carcinoma, KIRC - Kidney renal clear cell carcinoma, KIRP - Kidney renal papillary cell carcinoma, LAML - Acute Myeloid Leukemia, LIHC - Liver hepatocellular carcinoma, LUAD - Lung adenocarcinoma, LUSC - Lung squamous cell carcinoma, READ - Rectum adenocarcinoma, STAD - Stomach adenocarcinoma, and UCEC - Uterine Corpus Endometrioid Carcinoma), 2 deep sequencing Gene Expression Omnibus (GEO) datasets (GSE31999 and GSE37765), and 24 microarray datasets from GEO (GSE2564, GSE9234, GSE11255, GSE12250, GSE14224, GSE14473, GSE14794, GSE14834, GSE15387, GSE15745, GSE16558, GSE16654, GSE16759, GSE17306, GSE17491, GSE17498, GSE18155, GSE18693, GSE18899, GSE19350, GSE20692, GSE21032_1, GSE21032_2, GSE21321). Dataset correlation was done through MiRGator v3 software [77].

Two HPD-miRNA (hsa-mir-5007 and hsa-mir-4634) had no targets that passed quality control filters and were excluded from downstream analyses. The number of targets that passed quality control filters for the remaining HPD-miRNA were hsa-mir-202 (1145 targets), hsa-mir-1304 (710 targets), hsa-mir-1269a (488 targets), hsa-mir-4482-1 (144 targets), hsa-mir-449c (69 targets), and hsa-mir-4707 (1 target). Each unique mRNA target was individually analyzed using bootstrapping analysis to test for significant enrichment in HPD-miRNA. Specifically, targets were identified for (n = 2000) randomly chosen sets of 8 miRNA. The distribution of random targets was then compared to the 2,139 targets regulated by the 8 HPD-miRNA, and 72 targets were significantly enriched for regulation by HPD-miRNA (p < 0.05) (Figure 4 and Additional file 1: Table S5). Enriched targets were then annotated with the PANTHER pathway software (Additional file 2: Figure S4).

MiRNA target prediction was done with and without the novel “CT” deletion observed in hsa-mir-4640 using DIANA prediction software [73]. Canonical hsa-mir-4640 targets were compared among 6 different target prediction algorithms (TargetScan, miRNAorg, Microcosm Targets, PITA, PICTAR and miRDB) using MiRGator v3 software [77].

Experimentally validated gene targets were analyzed for reported genome wide significance in a GWAS study using the NHGRI GWAS catalog (Additional file 1: Table S6) [78]. We determine the extent of linkage disequilibrium (LD) in the regions surrounding the 8 HPD-miRNA SNPs, by using common tagging SNPs from the Illumina 1 M-Duo array in African samples (Additional file 2: Figure S6). The program PHASE v.2.1, which implements a Bayesian statistical method [79], was used to reconstruct multi-site haplotypes from genotype data for 15 SNPs from the Illumina 1 M-Duo SNP array flanking the miRNAs of interest on chromosomes 4, 5, 10,11,13 and 14 in 697 African individuals (Additional file 2: Figure S6). Haploview [80] was then used to calculate pairwise measures of LD among SNP loci, creating a graphical representation of the LD relationships among these loci (Additional file 2: Figure S6).

Functional enrichment

Genomic functional enrichment was analyzed with GOEAST software to determine statistically overrepresented GO terms within our set of 2,139 mRNA gene targets [46]. GO terms include annotated biological processes, molecular functions, and cellular components. The GOEAST algorithm assumes genes should be evenly distributed across the genome. MiRNA, as a class, may favor targeting a specific range of biological processes. To address this fact we perform miRNA functional enrichment.

MiRNA functional enrichment identifies biological pathways overrepresented by HPD-miRNA as compared to a set of 8 randomly sampled miRNA. We found the 2,139 mRNA gene targets of HPD-miRNA to be involved in 5,475 annotated biological processes using the ENSEMBL Gene Ontology database (Figure 4) [81]. For miRNA functional enrichment, the null distribution was created by, (1) randomly resampling 8 miRNA from the mirbase database, (2) identifying all mRNA targets (meeting the quality control filters described above in Messenger RNA Target Enrichment and Analysis), (3) identifying all biological processes associated with mRNA targets (using the ENSEMBL Gene Ontology database). Steps 1–3 are replicated n = 2000 times. The frequency of obtaining each of the 5,475 biological processes seen in HPD-miRNA is calculated. Finally, the actual frequency of each of the 5,475 biological processes associated with HPD-miRNAs is then compared to the null distribution above to ascertain enrichment of biological processes in our sample. Biological processes overrepresented in HPD-miRNA, having p-values < 0.05, were considered significantly enriched in HPD-miRNA and further clustered by gene ontology similarities using Revigo software [82].

Disease association

Disease associations were assessed through the MiRGator v2 disease database [47]. Disease associations are defined as differentially expressed miRNA in a published gene expression profile [47]. The significance of disease association results was tested by bootstrapping 31 miRNAs from mirbase and testing for disease association, n = 10,000 times. Network analysis was done using R with disease network visualization performed with Gephi software [83].

Abbreviations

miRNA:

microRNAs

PD-miRNA:

Population-differentiated miRNA

HPD-miRNA:

Highly population-differentiated miRNA

mRNA:

Messenger RNA

SNP:

Single nucleotide polymorphism

indels:

Insertions and deletions.

References

  1. Jansson MD, Lund AH: MicroRNA and cancer. Mol Oncol. 2012, 6 (6): 590-610. 10.1016/j.molonc.2012.09.006.

    Article  CAS  PubMed  Google Scholar 

  2. Friedman RC, Farh KK, Burge CB, Bartel DP: Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19 (1): 92-105.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Huang RS, Gamazon ER, Ziliak D, Wen Y, Im HK, Zhang W, Wing C, Duan S, Bleibel WK, Cox NJ, Dolan ME: Population differences in microRNA expression and biological implications. RNA Biol. 2011, 8 (4): 692-701. 10.4161/rna.8.4.16029.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Lewis BP, Burge CB, Bartel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.

    Article  CAS  PubMed  Google Scholar 

  5. Lu J, Clark AG: Impact of microRNA regulation on variation in human gene expression. Genome Res. 2012, 22 (7): 1243-1254. 10.1101/gr.132514.111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Xie Z, Allen E, Fahlgren N, Calamar A, Givan SA, Carrington JC: Expression of Arabidopsis MIRNA genes. Plant Physiol. 2005, 138 (4): 2145-2154. 10.1104/pp.105.062943.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Muniategui A, Pey J, Planes FJ, Rubio A: Joint analysis of miRNA and mRNA expression data. Brief Bioinform. 2013, 14 (3): 263-278. 10.1093/bib/bbs028.

    Article  CAS  PubMed  Google Scholar 

  8. Huntzinger E, Izaurralde E: Gene silencing by microRNAs: contributions of translational repression and mRNA decay. Nat Rev Genet. 2011, 12 (2): 99-110. 10.1038/nrg2936.

    Article  CAS  PubMed  Google Scholar 

  9. Chekulaeva M, Filipowicz W: Mechanisms of miRNA-mediated post-transcriptional regulation in animal cells. Curr Opin Cell Biol. 2009, 21 (3): 452-460. 10.1016/j.ceb.2009.04.009.

    Article  CAS  PubMed  Google Scholar 

  10. Flynt AS, Lai EC: Biological principles of microRNA-mediated regulation: shared themes amid diversity. Nat Rev Genet. 2008, 9 (11): 831-842. 10.1038/nrg2455.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Hobert O: Gene regulation by transcription factors and microRNAs. Science. 2008, 319 (5871): 1785-1786. 10.1126/science.1151651.

    Article  CAS  PubMed  Google Scholar 

  12. Filipowicz W, Bhattacharyya SN, Sonenberg N: Mechanisms of post-transcriptional regulation by microRNAs: are the answers in sight?. Nat Rev Genet. 2008, 9 (2): 102-114.

    Article  CAS  PubMed  Google Scholar 

  13. Eulalio A, Huntzinger E, Izaurralde E: Getting to the root of miRNA-mediated gene silencing. Cell. 2008, 132 (1): 9-14. 10.1016/j.cell.2007.12.024.

    Article  CAS  PubMed  Google Scholar 

  14. Berezikov E, Guryev V, van de Belt J, Wienholds E, Plasterk RH, Cuppen E: Phylogenetic shadowing and computational identification of human microRNA genes. Cell. 2005, 120 (1): 21-24. 10.1016/j.cell.2004.12.031.

    Article  CAS  PubMed  Google Scholar 

  15. Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36 (Database issue): D154-D158.

    CAS  PubMed  Google Scholar 

  16. Chugh PE, Sin SH, Ozgur S, Henry DH, Menezes P, Griffith J, Eron JJ, Damania B, Dittmer DP: Systemically Circulating Viral and Tumor-Derived MicroRNAs in KSHV-Associated Malignancies. PLoS Pathog. 2013, 9 (7): e1003484-10.1371/journal.ppat.1003484.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Farazi TA, Hoell JI, Morozov P, Tuschl T: MicroRNAs in human cancer. Adv Exp Med Biol. 2013, 774: 1-20. 10.1007/978-94-007-5590-1_1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Nana-Sinkam SP, Croce CM: Clinical applications for microRNAs in cancer. Clin Pharmacol Ther. 2013, 93 (1): 98-104. 10.1038/clpt.2012.192.

    Article  CAS  PubMed  Google Scholar 

  19. Chicago ARHLCsHo: Circulating microRNAs as Disease Markers in Pediatric Cancers. ClinicalTrialsgov [Internet]. 2012, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT01541800:NLM Identifier: NCT01541800

    Google Scholar 

  20. Florida UoS: Identification of Plasma miRNAs as Potential Biomarkers in Asthma Exacerbation. ClinicalTrialsgov [Internet]. 2012, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT01631760:NLM Identifier: NCT01631760

    Google Scholar 

  21. Group CsO: Studying Biomarkers in Cell Samples From Patients With Acute Myeloid Leukemia. ClinicalTrialsgov [Internet]. 2010, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT01057199:NLM Identifier: NCT01057199

    Google Scholar 

  22. Hospital C: Circulating microRNAs as Biomarkers of Sepsis. ClinicalTrialsgov [Internet]. 2000, Bethesda (MD): National Library of Medicine (US), [cited 2002 Feb 27] 2009, Available from: http://clinicaltrials.gov/show/NCT00862290:NLM Identifier: NCT00862290

    Google Scholar 

  23. Hospital NU: Profiling of Original Cellular and Humoral Biomarkers of Type 1 Diabetes (Lymphoscreen). ClinicalTrialsgov [Internet]. 2012, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT01042301:NLM Identifier:NCT01042301

    Google Scholar 

  24. Hospital WU: Micro-RNA Expression Profiles in High Risk Prostate Cancer. ClinicalTrialsgov [Internet]. 2010, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT01220427:NLM Identifier: NCT01220427

    Google Scholar 

  25. IRCCS CSdS: Gene Expression Profiles and Metformin Efficacy in Type 2 Diabetes. ClinicalTrialsgov [Internet]. 2012, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT01334684:NLM Identifier:NCT01334684

    Google Scholar 

  26. Party ELCW: Biological Factors Predicting Response to Chemotherapy in Advanced Non Small Cell Lung Cancer. ClinicalTrialsgov [Internet]. 2009, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT00864266:NLM Identifier: NCT00864266

    Google Scholar 

  27. Regaud IC: Circulating miRNAs as Biomarkers of Hormone Sensitivity in Breast Cancer (MIRHO). ClinicalTrialsgov [Internet]. 2012, Bethesda (MD): National Library of Medicine (US), Available from: http://clinicaltrials.gov/show/NCT01612871:NLM Identifier: NCT01612871

    Google Scholar 

  28. Zhang ZC, Huang Y, Wang XJ, Wang M, Ma LL: Expression of circulating microRNAs in patients with bladder urothelial carcinoma. Beijing Da Xue Xue Bao. 2013, 45 (4): 532-536.

    CAS  PubMed  Google Scholar 

  29. Markou A, Sourvinou I, Vorkas PA, Yousef GM, Lianidou E: Clinical evaluation of microRNA expression profiling in non small cell lung cancer. Lung Cancer. 2013, 81 (3): 388-396. 10.1016/j.lungcan.2013.05.007.

    Article  CAS  PubMed  Google Scholar 

  30. Wei C, Henderson H, Spradley C, Li L, Kim IK, Kumar S, Hong N, Arroliga AC, Gupta S: Circulating miRNAs as potential marker for pulmonary hypertension. PLoS One. 2013, 8 (5): e64396-10.1371/journal.pone.0064396.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Zhi F, Cao X, Xie X, Wang B, Dong W, Gu W, Ling Y, Wang R, Yang Y, Liu Y: Identification of circulating microRNAs as potential biomarkers for detecting acute myeloid leukemia. PLoS One. 2013, 8 (2): e56718-10.1371/journal.pone.0056718.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Zhang H, Li QY, Guo ZZ, Guan Y, Du J, Lu YY, Hu YY, Liu P, Huang S, Su SB: Serum levels of microRNAs can specifically predict liver injury of chronic hepatitis B. World J Gastroenterol. 2012, 18 (37): 5188-5196.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Heegaard NH, Schetter AJ, Welsh JA, Yoneda M, Bowman ED, Harris CC: Circulating micro-RNA expression profiles in early stage nonsmall cell lung cancer. Int J Cancer. 2012, 130 (6): 1378-1386. 10.1002/ijc.26153.

    Article  CAS  PubMed  Google Scholar 

  34. Zhao H, Shen J, Medico L, Wang D, Ambrosone CB, Liu S: A pilot study of circulating miRNAs as potential biomarkers of early stage breast cancer. PLoS One. 2010, 5 (10): e13735-10.1371/journal.pone.0013735.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Ramsay M: Africa: continent of genome contrasts with implications for biomedical research and health. FEBS Lett. 2012, 586 (18): 2813-2819. 10.1016/j.febslet.2012.07.061.

    Article  CAS  PubMed  Google Scholar 

  36. Campbell MC, Tishkoff SA: African genetic diversity: implications for human demographic history, modern human origins, and complex disease mapping. Annu Rev Genomics Hum Genet. 2008, 9: 403-433. 10.1146/annurev.genom.9.081307.164258.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Lambert CA, Tishkoff SA: Genetic structure in African populations: implications for human demographic history. Cold Spring Harb Symp Quant Biol. 2009, 74: 395-402. 10.1101/sqb.2009.74.053.

    Article  CAS  PubMed  Google Scholar 

  38. Tishkoff SA, Reed FA, Friedlaender FR, Ehret C, Ranciaro A, Froment A, Hirbo JB, Awomoyi AA, Bodo JM, Doumbo O, Ibrahim M, Juma AT, Kotze MJ, Lema G, Moore JH, Mortensen H, Nyambo TB, Omar SA, Powell K, Pretorius GS, Smith MW, Thera MA, Wambebe C, Weber JL, Williams SM: The genetic structure and history of Africans and African Americans. Science. 2009, 324 (5930): 1035-1044. 10.1126/science.1172257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Lachance J, Vernot B, Elbers CC, Ferwerda B, Froment A, Bodo JM, Lema G, Fu W, Nyambo TB, Rebbeck TR, Zhang K, Akey JM, Tishkoff SA: Evolutionary history and adaptation from high-coverage whole-genome sequences of diverse African hunter-gatherers. Cell. 2012, 150 (3): 457-469. 10.1016/j.cell.2012.07.009.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Mu XJ, Lu ZJ, Kong Y, Lam HY, Gerstein MB: Analysis of genomic variation in non-coding elements using population-scale sequencing data from the 1000 Genomes Project. Nucleic Acids Res. 2011, 39 (16): 7058-7076. 10.1093/nar/gkr342.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Quach H, Barreiro LB, Laval G, Zidane N, Patin E, Kidd KK, Kidd JR, Bouchier C, Veuille M, Antoniewski C, Quintana-Murci L: Signatures of purifying and local positive selection in human miRNAs. Am J Hum Genet. 2009, 84 (3): 316-327. 10.1016/j.ajhg.2009.01.022.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Hoffman AE, Liu R, Fu A, Zheng T, Slack F, Zhu Y: Targetome profiling, pathway analysis and genetic association study implicate miR-202 in lymphomagenesis. Cancer Epidemiol Biomarkers Prev. 2013, 22 (3): 327-336. 10.1158/1055-9965.EPI-12-1131-T.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Schuster SC, Miller W, Ratan A, Tomsho LP, Giardine B, Kasson LR, Harris RS, Petersen DC, Zhao F, Qi J, Alkan C, Kidd JM, Sun Y, Drautz DI, Bouffard P, Muzny DM, Reid JG, Nazareth LV, Wang Q, Burhans R, Riemer C, Wittekindt NE, Moorjani P, Tindall EA, Danko CG, Teo WS, Buboltz AM, Zhang Z, Ma Q, Oosthuysen A, et al: Complete Khoisan and Bantu genomes from southern Africa. Nature. 2010, 463 (7283): 943-947. 10.1038/nature08795.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Sherry STWM, Kholodov M, Baker J, Phan L, Smigielski EM, Sirotkin K: dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 2001, 29 (1): 308-311. 10.1093/nar/29.1.308.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Akey JM, Zhang G, Zhang K, Jin L, Shriver MD: Interrogating a high-density SNP map for signatures of natural selection. Genome Res. 2002, 12 (12): 1805-1814. 10.1101/gr.631202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Zheng Q, Wang XJ: GOEAST: a web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 2008, 36 (Web Server issue): W358-W363.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Cho S, Jun Y, Lee S, Choi HS, Jung S, Jang Y, Park C, Kim S, Kim W: miRGator v2.0: an integrated system for functional investigation of microRNAs. Nucleic Acids Res. 2011, 39 (Database issue): D158-D162.

    Article  CAS  PubMed  Google Scholar 

  48. Ryan BM, Robles AI, Harris CC: Genetic variation in microRNA networks: the implications for cancer research. Nat Rev Cancer. 2010, 10 (6): 389-402. 10.1038/nrc2867.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Li J, Liu Y, Xin X, Kim TS, Cabeza EA, Ren J, Nielsen R, Wrana JL, Zhang Z: Evidence for positive selection on a number of MicroRNA regulatory interactions during recent human evolution. PLoS Genet. 2012, 8 (3): e1002578-10.1371/journal.pgen.1002578.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Olson JE, Wang X, Goode EL, Pankratz VS, Fredericksen ZS, Vierkant RA, Pharoah PD, Cerhan JR, Couch FJ: Variation in genes required for normal mitosis and risk of breast cancer. Breast Cancer Res Treat. 2010, 119 (2): 423-430. 10.1007/s10549-009-0386-1.

    Article  CAS  PubMed  Google Scholar 

  51. (CDC) CfDCaP: Vital signs: racial disparities in breast cancer severity--United States, 2005-2009. MMWR Morb Mortal Wkly Rep. 2012, 61 (45): 922-926.

    Google Scholar 

  52. Al-Refaie WB, Tseng JF, Gay G, Patel-Parekh L, Mansfield PF, Pisters PW, Yao JC, Feig BW: The impact of ethnicity on the presentation and prognosis of patients with gastric adenocarcinoma. Results from the National Cancer Data Base. Cancer. 2008, 113 (3): 461-469.

    PubMed  Google Scholar 

  53. Chornokur G, Amankwah EK, Schildkraut JM, Phelan CM: Global ovarian cancer health disparities. Gynecol Oncol. 2012, 129 (1): 258-264.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Sun M, Abdollah F, Liberman D, Abdo A, Thuret R, Tian Z, Shariat SF, Montorsi F, Perrotte P, Karakiewicz PI: Racial disparities and socioeconomic status in men diagnosed with testicular germ cell tumors: a survival analysis. Cancer. 2011, 117 (18): 4277-4285. 10.1002/cncr.25969.

    Article  PubMed  Google Scholar 

  55. Schrauder MG, Strick R, Schulz-Wendtland R, Strissel PL, Kahmann L, Loehberg CR, Lux MP, Jud SM, Hartmann A, Hein A, Bayer CM, Bani MR, Richter S, Adamietz BR, Wenkel E, Rauh C, Beckmann MW, Fasching PA: Circulating micro-RNAs as potential blood-based markers for early stage breast cancer detection. PLoS One. 2012, 7 (1): e29770-10.1371/journal.pone.0029770.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. He B, Pan Y, Cho WC, Xu Y, Gu L, Nie Z, Chen L, Song G, Gao T, Li R, Wang S: The Association between Four Genetic Variants in MicroRNAs (rs11614913, rs2910164, rs3746444, rs2292832) and Cancer Risk: Evidence from Published Studies. PLoS One. 2012, 7 (11): e49032-10.1371/journal.pone.0049032.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Smith RA, Jedlinski DJ, Gabrovska PN, Weinstein SR, Haupt L, Griffiths LR: A genetic variant located in miR-423 is associated with reduced breast cancer risk. Cancer Genomics Proteomics. 2012, 9 (3): 115-118.

    CAS  PubMed  Google Scholar 

  58. Jedlinski DJ, Gabrovska PN, Weinstein SR, Smith RA, Griffiths LR: Single nucleotide polymorphism in hsa-mir-196a-2 and breast cancer risk: a case control study. Twin Res Hum Genet. 2011, 14 (5): 417-421. 10.1375/twin.14.5.417.

    Article  PubMed  Google Scholar 

  59. Su JL, Chen PB, Chen YH, Chen SC, Chang YW, Jan YH, Cheng X, Hsiao M, Hung MC: Downregulation of microRNA miR-520h by E1A contributes to anticancer activity. Cancer Res. 2010, 70 (12): 5096-5108. 10.1158/0008-5472.CAN-09-4148.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  60. Li X, Pan YZ, Seigel GM, Hu ZH, Huang M, Yu AM: Breast cancer resistance protein BCRP/ABCG2 regulatory microRNAs (hsa-miR-328, -519c and -520h) and their differential expression in stem-like ABCG2+ cancer cells. Biochem Pharmacol. 2011, 81 (6): 783-792. 10.1016/j.bcp.2010.12.018.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Wang F, Xue X, Wei J, An Y, Yao J, Cai H, Wu J, Dai C, Qian Z, Xu Z, Miao Y: hsa-miR-520h downregulates ABCG2 in pancreatic cancer cells to inhibit migration, invasion, and side populations. Br J Cancer. 2010, 103 (4): 567-574. 10.1038/sj.bjc.6605724.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Yu YH, Chen HA, Chen PS, Cheng YJ, Hsu WH, Chang YW, Chen YH, Jan Y, Hsiao M, Chang TY, Liu YH, Jeng YM, Wu CH, Huang MT, Su YH, Hung MC, Chien MH, Chen CY, Kuo ML, Su JL: MiR-520h-mediated FOXC2 regulation is critical for inhibition of lung cancer progression by resveratrol. Oncogene. 2012, 32 (4): 431-443.

    Article  PubMed  Google Scholar 

  63. Kim YW, Kim EY, Jeon D, Liu JL, Kim HS, Choi JW, Ahn WS: Differential microRNA expression signatures and cell type-specific association with Taxol resistance in ovarian cancer cells. Drug Des Devel Ther. 2014, 8: 293-314.

    CAS  PubMed  PubMed Central  Google Scholar 

  64. Pencheva N, Tran H, Buss C, Huh D, Drobnjak M, Busam K, Tavazoie SF: Convergent multi-miRNA targeting of ApoE drives LRP1/LRP8-dependent melanoma metastasis and angiogenesis. Cell. 2012, 151 (5): 1068-1082. 10.1016/j.cell.2012.10.028.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Long Q, Johnson BA, Osunkoya AO, Lai YH, Zhou W, Abramovitz M, Xia M, Bouzyk MB, Nam RK, Sugar L, Stanimirovic A, Williams DJ, Leyland-Jones BR, Seth AK, Petros JA, Moreno CS: Protein-coding and microRNA biomarkers of recurrence of prostate cancer following radical prostatectomy. Am J Pathol. 2011, 179 (1): 46-54. 10.1016/j.ajpath.2011.03.008.

    Article  PubMed  PubMed Central  Google Scholar 

  66. Richardson K, Lai CQ, Parnell LD, Lee YC, Ordovas JM: A genome-wide survey for SNPs altering microRNA seed sites identifies functional candidates in GWAS. BMC Genomics. 2011, 12: 504-10.1186/1471-2164-12-504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Bulik-Sullivan B, Selitsky S, Sethupathy P: Prioritization of genetic variants in the microRNA regulome as functional candidates in genome-wide association studies. Hum Mutat. 2013, 34 (8): 1049-1056. 10.1002/humu.22337.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Drmanac R, Sparks AB, Callow MJ, Halpern AL, Burns NL, Kermani BG, Carnevali P, Nazarenko I, Nilsen GB, Yeung G, Dahl F, Fernandez A, Staker B, Pant KP, Baccash J, Borcherding AP, Brownley A, Cedeno R, Chen L, Chernikoff D, Cheung A, Chirita R, Curson B, Ebert JC, Hacker CR, Hartlage R, Hauser B, Huang S, Jiang Y, Karpinchyk V, et al: Human genome sequencing using unchained base reads on self-assembling DNA nanoarrays. Science. 2010, 327 (5961): 78-81. 10.1126/science.1181498.

    Article  CAS  PubMed  Google Scholar 

  69. Roach JC, Glusman G, Smit AF, Huff CD, Hubley R, Shannon PT, Rowen L, Pant KP, Goodman N, Bamshad M, Shendure J, Drmanac R, Jorde LB, Hood L, Galas DJ: Analysis of genetic inheritance in a family quartet by whole-genome sequencing. Science. 2010, 328 (5978): 636-639. 10.1126/science.1186802.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, Belmont JW, Boudreau A, Hardenbol P, Leal SM, Pasternak S, Wheeler DA, Willis TD, Yu F, Yang H, Zeng C, Gao Y, Hu H, Hu W, Li C, Lin W, Liu S, Pan H, Tang X, Wang J, Wang W, Yu J, Zhang B, Zhang Q, Zhao H, et al: A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007, 449 (7164): 851-861. 10.1038/nature06258.

    Article  CAS  PubMed  Google Scholar 

  71. Ritchie W, Gao D, Rasko JE: Defining and providing robust controls for microRNA prediction. Bioinformatics. 2012, 28 (8): 1058-1061. 10.1093/bioinformatics/bts114.

    Article  CAS  PubMed  Google Scholar 

  72. Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136 (2): 215-233. 10.1016/j.cell.2009.01.002.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Maragkakis M, Alexiou P, Papadopoulos GL, Reczko M, Dalamagas T, Giannopoulos G, Goumas G, Koukis E, Kourtis K, Simossis VA, Sethupathy P, Vergoulis T, Koziris N, Sellis T, Tsanakas P, Hatzigeorgiou AG: Accurate microRNA target prediction correlates with protein repression levels. BMC Bioinformatics. 2009, 10: 295-10.1186/1471-2105-10-295.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Watterson GA: On the number of segregating sites in genetical models without recombination. Theor Popul Biol. 1975, 7 (2): 256-276. 10.1016/0040-5809(75)90020-9.

    Article  CAS  PubMed  Google Scholar 

  75. Cockerham CC, Weir BS: Covariances of relatives stemming from a population undergoing mixed self and random mating. Biometrics. 1984, 40 (1): 157-164. 10.2307/2530754.

    Article  CAS  PubMed  Google Scholar 

  76. Holsinger KE, Weir BS: Genetics in geographically structured populations: defining, estimating and interpreting F(ST). Nat Rev Genet. 2009, 10 (9): 639-650. 10.1038/nrg2611.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Cho S, Jang I, Jun Y, Yoon S, Ko M, Kwon Y, Choi I, Chang H, Ryu D, Lee B, Kim VN, Kim W, Lee S: MiRGator v3.0: a microRNA portal for deep sequencing, expression profiling and mRNA targeting. Nucleic Acids Res. 2013, 41 (Database issue): D252-D257.

    Article  CAS  PubMed  Google Scholar 

  78. Hindorff LA, MacArthur J, Morales J, Junkins HA, Hall PN, Klemm AK, Manolio TA: A Catalog of Published Genome-Wide Association Studies. Available at: http://www.genome.gov/gwastudies Accessed July 2013

  79. Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001, 68 (4): 978-989. 10.1086/319501.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  80. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.

    Article  CAS  PubMed  Google Scholar 

  81. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29.

    CAS  PubMed  Google Scholar 

  82. Supek F, Bosnjak M, Skunca N, Smuc T: REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011, 6 (7): e21800-10.1371/journal.pone.0021800.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  83. Bastian M, Heymann S, Jacomy M: Gephi: an open source software for exploring and manipulating networks. International AAAI Conference on Weblogs and Social Media. 2009, Available at: https://gephi.org/

    Google Scholar 

Pre-publication history

Download references

Acknowledgements

We would like to thank Dr. Joseph LaChance and Dr. Alessia Ranciaro for useful discussions and input on this manuscript. We would especially like to thank all of the participants who generously donated their blood for this work. This study was supported by a NIEHS Pioneer Award (DP1ES022577) to S.A.T and a NIH postdoctoral fellowship (K12GM081259) to R.A.R.-G.

Author information

Authors and Affiliations

  1. Department of Genetics, University of Pennsylvania, Philadelphia, PA, 19104, USA

    Renata A Rawlings-Goss, Michael C Campbell & Sarah A Tishkoff

  2. Department of Biology, University of Pennsylvania, Philadelphia, PA, 19104, USA

    Sarah A Tishkoff

  3. Department of Biostatistics, Yale University, New Haven, CT, 06520, USA

    Michael C Campbell

Authors
  1. Renata A Rawlings-Goss
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Michael C Campbell
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Sarah A Tishkoff
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Renata A Rawlings-Goss or Sarah A Tishkoff.

Additional information

Competing interests

The authors declare that they have no competing interests. Whole genome variant data has been made available through dbSNP, dbGAP and Complete Genomics.

Authors’ contributions

RR conceived of the study, performed analysis, and wrote the manuscript. MC performed linkage disequilibrium analysis and contributed to manuscript revisions. ST supervised the project and contributed to the writing of the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

12920_2014_509_MOESM1_ESM.xls

Additional file 1: Table S1: All miRNA variants. Table S2. Novel variants within miRNA seed sequences. Table S3. Predicted changes in gene targets of hsa-mi-4640-3p with novel "CT" deletion. Table S4. All PD-miRNA with FST above the 95th percentile between African and non-Africans. Table S5. Significantly enriched targets of HPD-miRNA. Table S6. Significantly enriched genes that were reported in genome-wide association studies. Table S7. GOEAST: genetic functional enrichment of the 2,139 mRNA gene targets. Table S8. MiRNA functional enrichment of the 2,139 mRNA gene targets. Table S9. Allele Frequency for hsa-mir-202 T-allele. T-allele has known effect on miRNA expression and a protective effect on breast cancer mortality [42]. N is the number of chromosomes in each sample group. (XLS 680 KB)

12920_2014_509_MOESM2_ESM.docx

Additional file 2: Figure S1: Pairwise FST for populations grouped by continent. Figure S2. Distribution of p-values for pairwise FST’s measured between African and non-African populations. Figure S3. Allele frequencies for the 4 PD-variants in miRNA seed sequences. Figure S4. Biological functions of the 72 significant gene targets of HPD-miRNA. Figure S5. Linkage disequilibrium plots for HPD-SNPs in African populations based on the Illumina 1M Duo. Figure S6. Tissue specific expression of hsa-mir-202. (DOCX 7 MB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Rawlings-Goss, R.A., Campbell, M.C. & Tishkoff, S.A. Global population-specific variation in miRNA associated with cancer risk and clinical biomarkers. BMC Med Genomics 7, 53 (2014). https://doi.org/10.1186/1755-8794-7-53

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1755-8794-7-53

Keywords

BMC Medical Genomics

ISSN: 1755-8794

Contact us