ORIGINAL RESEARCH article

Front. Physiol., 07 June 2018

Sec. Integrative Physiology

Volume 9 - 2018 | https://doi.org/10.3389/fphys.2018.00705

No Association Between MicroRNA-608 rs4919510 G>C Polymorphism and Digestive System Cancers Susceptibility: A Meta-Analysis Based on 10,836 Individuals

  • 1. Department of Endocrinology, Taihe Hospital, Hubei University of Medicine, Shiyan, China

  • 2. Department of Oral and Maxillary Surgery, Gui Zhou Provincial People's Hospital, Guiyang, China

  • 3. Trade Union, Taihe Hospital, Hubei University of Medicine, Shiyan, China

  • 4. Department of Stomatology, Center for Evidence-Based Medicine, Taihe Hospital, Hubei University of Medicine, Shiyan, China

Abstract

Previous epidemiologic studies have revealed a possible association between microRNA-608 rs4919510 G>C polymorphism and digestive system cancers (DSCs) risk, but the results were not consistent. We therefore performed an updated meta-analysis to explore the association between microRNA-608 rs4919510 G>C polymorphism and DSCs risk. Crude odds ratios (ORs) with 95% confidence intervals (CIs) were calculated to assess the relationship between the microRNA-608 rs4919510 G>C polymorphism and DSCs risk. Heterogeneity, cumulative analyses, sensitivity analyses, and publication bias were also conducted to examine the statistical power. Eight published articles with nine independent case-control studies involving 10,836 individuals were included in this meta-analysis. Overall, no significant association was found between microRNA-608 rs4919510 G>C polymorphism and DSCs risk in general populations. But some significant protective effects were observed in the subgroup of Caucasian population group in three genetic models (C vs. G: OR = 0.82, 95% CI, 0.68–0.99, P = 0.03, I2 = 0%; CC vs. GG: OR = 0.59, 95% CI = 0.36–0.97, P = 0.04, I2 = 0%; GC+CC vs. GG: OR = 0.61, 95% CI = 0.37–0.99, P = 0.05, I2 = 0%). In summary, current evidence indicates that the microRNA-608 rs4919510 G>C polymorphism maybe an important factor of DSCs susceptibility, especially in Caucasian population.

Introduction

Digestive system cancers (DSCs), which comprise esophageal cancer (EC), gastric cancer (GC), hepatocellular cancer (HC), pancreatic cancer (PC) colorectal cancer (CRC) and other solid carcinoma, is one of the most common malignancies with increasing incidence and mortality worldwide (Torre et al., 2015). According to the recent published cancer statistics, there were approximately 310,440 new DSCs cases and 157,700 DSCs-related deaths in United States in 2017 (Siegel et al., 2017). The treatment and care of patients with DSCs result in a heavy economic and psychological burden to both society and the patient's family (Hsu et al., 2017; Jinjuvadia et al., 2017). Surgery, radiotherapy, chemotherapy and other treatments may cause extensive damage to patients' tissues and organs and may result in various complications that lead to local dysfunction and seriously decrease the quality of life of patients (Nederlof et al., 2016; Bosch et al., 2017; Motoyama et al., 2017). Regrettably, no clear explanation of the mechanism underlying DSCs development and the susceptibility of different patients exist. Some researches indicate that unhealthy life styles, cigarette and alcohol use, viral infection, local inflammation and stress can trigger DSCs development (Erren et al., 2016; Gao et al., 2017; Jarzynski et al., 2017; Jayasekara et al., 2018).

Evidences suggested that abnormal gene expression and an alteration of amino acid structure can lead to abnormal biologic activities that result in an imbalance of normal physiologic function and induce the formation of tumors. MicroRNAs were discovered in 1993 and were defined as small noncoding RNA molecules with approximately 21–25 nucleotides in length and a characteristic double-stranded structure (Lee et al., 1993; Ambros, 2004). MicroRNAs always originate from endogenous transcripts and bind to imperfect complementary sequences of the 3′-untranslated regions (3'-UTR) of target mRNAs to regulate their post-transcriptional repression (Bartel, 2004). Some studies indicated that abnormal microRNA expression could result in tumor occurrence (Valeri et al., 2010; Tessitore et al., 2014). MicroRNA-608 is a newly discovered microRNA, as reported by Wang et al., who indicated that the expression of microRNA-608 was significantly reduced in hepatocellular cancer and its expression levels were associated with tumor size, differentiation, clinical stage, and overall prognosis (Wang K. et al., 2016). Further, other researchers found that the expression of microRNA-608 was significantly down-regulated in glioma stem cells (GSCs). And the up-regulation of microRNA-608 expression would inhibit the proliferation, migration, and invasion of GSCs and promotes their apoptosis (Wang Z. et al., 2016).

Single nucleotide polymorphisms (SNPs) are the most common type of genetic variation in people. SNPs alter gene function and/or expression, consequently affecting downstream biologic pathways and increasing cancer risk (Martini et al., 2013; Niu et al., 2015). In microRNA-608, rs4919510 G>C is the most common locus and has attracted increasing attention. Ye et al. conducted the first case-control study in 2008 and did not find any significant association between the microRNA-608 rs4919510 G>C polymorphism and esophageal cancer (Ye et al., 2008). Recently, a published meta-analysis addressed the association between microRNA-608 rs4919510 G>C polymorphism and cancer risk (Liu et al., 2016), but the results were not consistent with subsequent studies, especially in DSCs. Therefore, we conducted the present meta-analysis to more precisely and comprehensively assess of the association between microRNA-608 rs4919510 G>C polymorphism and DSCs risk.

Materials and methods

This present meta-analysis followed the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) statement (Moher et al., 2009). All included data were collected from published studies, and no ethical issues were involved.

Search strategy

Five electronic English databases (PubMed, Embase, Web of Science, CNKI and Wanfang) were searched for relevant studies that focused on the association between microRNA-608 rs4919510 G>C polymorphism and DSCs risk from inception up to January 1, 2018. Only studies that were written in English and Chinese were included. Moreover, the bibliographies of the collected studies and relevant reviews were retrospected to identify additional articles. The following search terms and strategy was used (e.g., in PubMed):

#1 microRNA 608

#2 microRNA-608

#3 mir 608

#4 mir-608

#5 rs4919510

#6 #1 OR #2 OR #3 OR #4 OR #5

#7 polymorphism

#8 variant

#9 mutation

#10 #7 OR #8 OR #9

#11 cancer

#12 tumor

#13 neoplasm

#14 #11 OR #12 OR #13

#15 #6 AND #10 AND #14

Eligibility criteria

Studies were included based on the following criteria: (1) only case-control studies that focusing on the association between the microRNA-608 rs4919510 G>C polymorphism and DSCs risk; and (2) studies had to provided sufficient frequency data on the genotype distribution to evaluate the crude odds ratios (ORs) and 95% confidence intervals (CIs); and; (3) studies had to be published only in English and Chinese; and (4) only publications with the largest or most recently updated sample data were included when there were some overlapping or duplicate publications on the same theme.

Data extraction and quality evaluation

Two researchers (Li and Song) independently reviewed and extracted relevant information from all included studies, including name of the first author, publishing date, country or region where the study was conducted, race, control design, sample sizes of the cases and controls, frequency data of genotype distribution, genotyping method, Hardy-Weinberg equilibrium (HWE) assessment in controls, and cancer type. Quality evaluation of the included studies was performed by the two researchers using the modified Newcastle-Ottawa scale (NOS). The scores ranged from 0 points (worst) to 11 points (best) (Table 4); studies with a score of 9 or higher were classified as high quality.

Statistical analysis

We calculated ORs and 95% CIs to assess the association between microRNA-608 rs4919510 G>C polymorphism and DSCs risk. All pooled genetic models were examined, including the allele contrast (C vs. G), co-dominant models (GC vs. GG and CC vs. GG), dominant model (GC+CC vs. GG), and recessive model (CC vs. GG+GC). Heterogeneity between the included studies was calculated using Cochran's Q-test and I2 statistical method (Huedo-Medina et al., 2006). The fixed-effect model was used when the I2-value was less than 40%; otherwise, a random-effects model was used (Mantel and Haenszel, 1959; DerSimonian, 1996). Subgroup analyses were performed according to race (Asian, Caucasian, and African), control design (population base, hospital base), cancer type, subject number, and NOS evaluation. Moreover, meta-regression was performed to interpret the between-group heterogeneity. Cumulative meta-analyses were performed to assess the continuous tendency in the results. Furthermore, sensitivity analyses were conducted to examine the stability of the results by removing each study one by one. Potential publication biases were examined using Egger's linear regression and Begg's funnel plots. Statistical analyses were performed using STATA version 14.0 (Stata Corporation, College Station, TX, USA) (Begg and Mazumdar, 1994; Egger et al., 1997). A two-sided P-value less than 0.05 was considered statistically significant.

Results

Study characteristics

A total of 164 articles were collected. The study selection process was presented in Figure 1. One hundred and fifty-six articles were excluded during the comprehensive screening procedures based on the article titles, abstracts, and full texts. Thus, eight articles including nine independent case-control studies involving 5,224 patients and 5,612 controls were identified and included in this meta-analysis according to our inclusion criteria (Ryan et al., 2012; Zhang, 2012; Kupcinskas et al., 2014a,b; Wang et al., 2014; Zhang et al., 2015; Jiang et al., 2016; Ying et al., 2016). Five case-control studies included 4,430 cases and 4,403 controls in Asian populations(Zhang, 2012; Wang et al., 2014; Zhang et al., 2015; Jiang et al., 2016; Ying et al., 2016), thee case-control studies involved 700 cases and 1,024 controls in Caucasian populations (Ryan et al., 2012; Kupcinskas et al., 2014a,b), and one case-control study involved 94 cases and 185 controls in an African American population (Ryan et al., 2012). For genotyping, four studies used the Taqman method (Ryan et al., 2012; Kupcinskas et al., 2014a,b), three studies used the MassARRAY method (Wang et al., 2014; Jiang et al., 2016; Ying et al., 2016), and the remaining two studies used other methods (SNaPshot, PCR-RFLP) (Zhang, 2012; Zhang et al., 2015). Moreover, five studies focused on colorectal cancer(Ryan et al., 2012; Zhang, 2012; Kupcinskas et al., 2014a; Ying et al., 2016), two studies focused on gastric cancer (Kupcinskas et al., 2014b; Jiang et al., 2016), one study focused on hepatocellular cancer (Wang et al., 2014) and other focused on esophageal cancer (Zhang et al., 2015). The genotype distributions in all control groups were satisfied with the HWE. The characteristics of the included studies are presented in Table 1.

Figure 1

Table 1

Genotype distribution
CaseControl
First authorYearCountry/RegionRacialSource of controlsCaseControlGGGCCCGGGCCCGenotyping methodsP for HWEMAF in controlTypeNOS
Zhang MW2012ChinaAsianPB4604651562218314722890PCR-RFLP0.920.44CRC9
Ryan BM-a2012USACaucasianMixed14524874890871169Taqman0.870.82CRC8
Ryan BM-b2012USAAfricanMixed94185124834289562Taqman0.390.59CRC8
Kupcinskas J-12014LatviaCaucasianHB36335025882501386251Taqman0.110.84GC8
Kupcinskas J-22014LatviaCaucasianHB1924267471381296318Taqman0.160.86CRC7
Wang R2014ChinaAsianPB993992304500189318497177MassARRAY0.480.43HC8
Zhang P2015ChinaAsianPB738882217384137291440151SNaPshot0.480.42ESCC9
Jiang J2016ChinaAsianHB894989278451165296483210MassARRAY0.620.46GC9
Ying HQ2016ChinaAsianHB1,3451,075232690423250512313MassARRAY0.150.53CRC8

Characteristics of included studies on MicroRNA-608 rs4919510 G>C polymorphisms and digestive system cancers risk.

HWE in control.

CRC, Colorectal cancer; GC, Gastric cancer; HC, Hepatocellular cancer; ESCC, Esophageal squamous cell carcinoma; MAF, Minor allele frequency in control group; PB, Population-based; HB, Hospital-based; Mixed, combined PB with HB.

Meta-analysis

Overall, no significant association between the microRNA-608 rs4919510 G>C polymorphism and DSCs risk was observed with the included studies (C vs. G: OR = 1.00, 95% CI = 0.91–1.10, P = 0.98, I2 = 53.4%; GC vs. GG: OR = 1.07, 95% CI = 0.93–1.24, P = 0.35, I2 = 43.0%; CC vs. GG: OR = 1.01, 95% CI = 0.83–1.24, P = 0.90, I2 = 56.4%; GC+CC vs. GG: OR = 1.05, 95% CI = 0.89–1.23, P = 0.58, I2 = 56.4%, (Figure 2); CC vs. GG+GC: OR = 0.99, 95% CI = 0.91–1.09, P = 0.89, I2 = 0%) (Supplementary Figure S1).

Figure 2

Subgroup and meta-regression analysis

The results of race diversity demonstrated some significant protective effects between microRNA-608 rs4919510 G>C polymorphism and DSCs risk in Caucasian populations (C vs. G: OR = 0.82, 95% CI = 0.68–0.99, P = 0.03, I2 = 0%; CC vs. GG: OR = 0.59, 95% CI = 0.36–0.97, P = 0.04, I2 = 0%; GC+CC vs. GG: OR = 0.61, 95% CI = 0.37–0.99, P = 0.05, I2 = 0%, Figure 2). But the other Subgroup analysis based on control design, cancer type, subject number, and NOS evaluation did not present any significant association between microRNA-608 rs4919510 G>C polymorphism and DSCs risk. All these statistical results were presented in Table 2. Moreover, meta-regression were conducted due to the slight heterogeneity between included studies existed, but the statistical test did not find any remarkable factors that contributed to the current heterogeneity (Table 3).

Table 2

N*C vs. GGC vs. GGCC vs. GGGC+CC vs. GGCC vs. GG+GC
OR95% CIPaI2(%)OR95% CIPaI2(%)OR95% CIPaI2(%)OR95% CIPaI2(%)OR95% CIPaI2(%)
Total91.000.91–1.100.9853.41.070.93-1.240.3543.01.010.83–1.240.9056.41.050.89–1.230.5856.40.990.91–1.090.890
RACIAL
 Asian51.040.95–1.150.4061.31.110.94–1.300.1857.51.090.88–1.350.4366.21.100.93–1.300.2567.31.030.93–1.140.6217.6
 Caucasian30.820.68–0.990.0300.650.39–1.090.1000.590.36–0.970.0400.610.37–0.990.0500.840.68–1.040.120
 African11.110.78–1.590.56NA1.180.55–2.520.67NA1.280.58–2.830.54NA1.220.59–2.520.59NA1.120.67–1.890.66NA
DESIGN
 PB31.050.96–1.140.2911.81.060.93–1.210.3801.090.92–1.300.3201.070.94–1.210.3015.31.050.91–1.220.500
 HB40.960.80–1.160.7075.81.020.71–1.470.9173.10.900.57–1.420.6580.00.970.66–1.430.8779.60.970.86–1.100.6337.3
 Mixed20.930.66–1.320.6944.61.020.55–1.910.9400.960.47–1.950.9018.71.000.55–1.811.0000.890.64–1.250.5120.3
TYPE
 CRC51.000.84–1.180.9658.61.110.83–1.500.4847.41.060.75–1.500.7451.0 .01.070.77–1.490.6759.51.020.89–1.160.820
 GC20.900.80–1.010.0800.810.45–1.440.4761.30.790.62–1.000.0537.90.770.44–1.340.3662.60.850.71–1.020.090
 Other21.080.98–1.180.1101.10–0.95–1.280.2001.160.96–1.410.1301.120.97–1.230.1201.100.921.290.310
SUBJECTS
  < 1,00050.900.79–1.010.0800.870.69–1.110.2700.810.61–1.070.1300.850.58–1.070.1700.890.75–1.050.170
 >1,00041.060.96–1.180.2464.01.150.98–1.360.0958.91.140.90–1.440.2869.71.150.96–1.370.1369.01.040.93–1.160.4931.3
NOS EVALUATION
 NOS ≥ 930.990.87–1.12.0.8351.01.040.91–1.190.600.90.950.75–1.230.7748.61.010.93–1.150.8137.60.940.81–1.090.4319.5
 NOS < 961.000.87–1.150.7855.91.070.82–1.410.6352.41.030.76–1.410.8455.61.040.78–1.380.8060.51.020.91–1.150.680

Summary ORs and 95% CI of MicroRNA-608 rs4919510 G>C polymorphisms and digestive system cancers risk.

*

Numbers of comparisons.

a

P-value for OR evaluation.

Table 3

C vs. GGC vs. GGCC vs. GGGC+CC vs. GGCC vs. GG+GC
tPtPtPtPtP
Racial−1.030.337−0.900.400−0.850.425−0.920.388−0.860.417
Design−0.710.5010.080.935−0.460.661−0.200.845−0.990.354
Type0.480.648−0.140.8910.210.8420.050.9600.560.594
Subjects2.010.0841.760.1221.750.1231.860.1051.480.184
NOS evaluation0.210.8410.470.6560.390.7110.330.7480.680.520

Characteristics of included covariate in meta-regression analysis in all five genetic models.

Table 4

CriteriaScore
REPRESENTATIVENESS OF CASES
Consecutive/randomly selected cases with clearly defined sampling frame2
Not consecutive/randomly selected case or without clearly defined sampling frame1
Not described0
SOURCE OF CONTROLS
Population- or Healthy-based2
Hospital-bases1
Not described0
HARDY-WEINBERG EQUILIBRIUM IN CONTROLS
Hardy-Weinberg equilibrium2
Hardy-Weinberg disequilibrium1
Not available0
GENOTYPING EXAMINATION
Genotyping done under “blinded” condition and repeated again2
Genotyping done under “blinded” condition or repeated again1
Unblinded done or not mentioned and unrepeated0
SUBJECTS
Number >1,0001
Number < 1,0000
ASSOCIATION ASSESSMENT
Assess association between genotypes and cancer with appropriate statistics and adjustment for confounders2
Assess association between genotypes and cancer with appropriate statistics and without adjustment for confounders1
Inappropriate statistics used0

Scale for quality evaluation.

Cumulative and sensitivity analyses

Cumulative analysis indicated a consistent tendency accompanied with the continuously published studies (Figure 3 for GC+CC vs. GG model) (Supplementary Figure S2). Sensitive analyses were performed by omitting each study one at a time, and the results did not indicate any significant changes (Figure 4 for GC+CC vs. GG model) (Supplementary Figure S3).

Figure 3

Figure 4

Publication bias

Publication bias was examined with Begg's test and no apparent asymmetry of the funnel plot was found (Figure 5 for GC+CC vs. GG model) (Supplementary Figure S4). These results were confirmed with Egger's test (C vs. G: P = 0.11; GC vs. GG: P = 0.24; CC vs. GG: P = 0.16; GC+CC vs. GG: P = 0.22; CC vs. GG+GC: P = 0.19).

Figure 5

Discussion

MicroRNAs is a family of small, noncoding, evolutionarily conserved endogenous RNAs consisting of 22 nucleotides that can regulate gene expression through a post-transcriptional pathway binding to the 3′-UTR of target mRNAs (Lagos-Quintana et al., 2001). MicroRNAs take part in numerous biologic processes, including cell proliferation, differentiation and apoptosis (Ambros, 2003). Evidences have shown that microRNAs play a composite function for oncogenes and tumor suppressor genes. Some mutations, such as SNPs in microRNAs or their binding site, may alter cancer susceptibility in various tumor types and different populations (Niu et al., 2014, 2015).

MicroRNA-608 is novel small RNA that has been shown to be significantly associated with cancer development, metastasis and recurrence. Wang et al. had found that the expression of miR-608 expression was significantly reduced in hepatocellular carcinoma and presented a negative correlation of the tumor size, differentiation, and so on (Wang K. et al., 2016). Yang et al. reported that the miR-608 could suppress the tumorigenesis of colon cancer cells by targeting NAA10 mRNA for degradation. And the overexpression of miR-608 could decrease NAA10 mRNA and protein levels, and thereby suppressed cell proliferation, migration, and cell-cycle progression (Yang et al., 2016).

Rs4919510 G>C is the most common SNP locus in microRNA-608 and might contribute a different free energy to its binding site of target genes, including INSR,CD4, GHR and RXRB (Landi et al., 2008). According the research of Huang et al., the variant genotypes (CG/GG) and variant G-allele were specifically associated with increased risk of HER2-positive subtype in breast cancer patients. HSF1 is a predicted target of miR-608 and the C-to-G substitution of rs4919510 might weaken the suppression of HSF1 mRNA by miR-608, leading to relatively high expression of HSF1 protein, which could up-regulate the expression of HER2 and facilitated tumorigenesis. The above evidence indicated that the alteration of free energy caused by gene mutation between miR-608 and targeted mRNA would influence the protein expression and change the cancer susceptibility (Huang et al., 2012). Since 2008, a series of case-control studies focused on microRNA-608 rs4919510 G>C polymorphism and DSCs risk were published, but the conclusions were inconsistent. Kupcinskas et al. explored the association between microRNA-608 rs4919510 G>C polymorphism and gastric cancer and found an apparently increased risk of gastric cancer in a Latvian population with GG genotype (GG vs. CC, OR = 2.34, 95% CI, 1.08–5.04) (Kupcinskas et al., 2014b). Conversely, Ying et al. evaluated the microRNA-608 rs4919510 G>C polymorphism in Chinese individuals and concluded that the rs4919510 G-allele and genotype GG were protective factors for stage 0 to II colorectal cancer (CRC) (GG vs. CC, OR = 0.70, 95% CI = 0.55–0.88) (Ying et al., 2016). However, other studies did not find any significant relationship between the microRNA-608 rs4919510 G>C polymorphism and DSCs risk, including these studies by Wang et al. (2014) and Zhang et al. (2015). How can we conduct a more precise assessment of the association between microRNA-608 rs4919510 G>C polymorphism and DSC risk with these inconsistent results? As we know, few studies and small sample size maybe the important reasons due to the confusing results. Therefore, we conducted the meta-analysis with nine published case-control studies to investigate the association between the microRNA-608 rs4919510 G>C polymorphism and DSCs susceptibility. Overall, our meta-analysis did not identify any significant association in all of the genetic models. To our knowledge, the potential real results might be confounded or even conceal by some subgroup factors, such as race diversity, control design, cancer type, subject number, and NOS evaluation. So, we subsequently conducted stratified analyses by these subgroup factors. The results indicate that the microRNA-608 rs4919510 G>C polymorphism play an important protective effect against DSCs development in Caucasian population only, but not in other subgroup analysis. It was supposed that the racial diversity maybe the critical factor that contributed to the different result. However, according to current results in general population and subgroups, the precise biological mechanism of the association between microRNA-608 rs4919510 G>C polymorphism and DSCs risk were remains unclear. To our minds, the possible explanation was that the microRNA-608 rs4919510 G>C polymorphism not participated in cancer susceptibility directly, but play a connected and synergism role during the development of cancer. All above results and speculations were needed to confirm with new molecular and epidemiological studies in various races.

In 2017, Liu et al. (2016) and Wu et al. (2017) published two meta-analysis assessing the association between the microRNA-608 rs4919510 G>C polymorphism and cancer risk. Their researches included 10 and 18 case-control studies respectively, which consisted of thyroid cancer, breast cancer, lung cancer and others. Liu et al. indicated that the presence of CC genotypes might be protective against tumorigenesis, especially in Caucasians. And Wu et al. suggested that CG genotypes would increase the cancer risk in the Chinese populations. However, the results of their meta-analysis focused on the general cancer, and the association between the microRNA-608 rs4919510 G>C polymorphism and DSCs risk was not conducted. To our knowledge, this is the first meta-analysis to explore the association between the microRNA-608 rs4919510 G>C polymorphism and DSCs risk. Compared with previous meta-analyses, the present meta-analysis used a more scientific retrieval strategy and a larger sample size. Additionally, more rigorous methodology, including cumulative and sensitivity analyses and quality evaluation with modified NOS, were used to estimate the genetic effects of the microRNA-608 rs4919510 G>C polymorphism on carcinogenesis. This provided a more accurate evaluation of the association between the microRNA-608 rs4919510 G>C polymorphism and DSCs risk.

However, some limitations of this analysis could not be avoided and should be addressed. First, some heterogeneity was apparent among the included studies, which might distort the results in the current meta-analysis. The meta-regression was conducted but not found any distinct factor that contributes to the current heterogeneity. Fortunately, the heterogeneity was partly alleviated in the subsequent stratified analysis, such as race diversity and cancer type. Second, this meta-analysis only focused on one SNP locus (microRNA-608 rs4919510 G>C), and the results were calculated without gene-gene and gene-environment risk factors, leading to a failure to interpret the potential interaction mechanisms. Third, almost all of the included studies focused on Asian and Caucasian populations, which would restrict the application of our results to other populations.

In summary, our meta-analysis demonstrated that the microRNA-608 rs4919510 G>C polymorphism might play an important role in DSCs susceptibility. Further studies in different races and regions with larger population sizes are needed to confirm our findings.

Statements

Author contributions

X-FL, J-KS, and Y-MN conceived the study. X-FL, J-KS, and J-WC searched the databases and extracted the data. J-KS, Y-QZ, and ML analyzed the data. X-FL, J-KS, and J-WC wrote the draft of the paper. JZ and Y-MN reviewed the manuscript. All the authors approved the final manuscript.

Funding

This study was supported by the Foundations of the Science and Technology Department of Hubei Province (No. 2016CFB567, 2014CFB364), and the Hubei Province health and family planning scientific research project (No. WJ2017F069, WJ2015Q041) and Taihe Hospital (2016BSQD02). The funders had no roles in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2018.00705/full#supplementary-material

References

  • 1

    AmbrosV. (2003). MicroRNA pathways in flies and worms: growth, death, fat, stress, and timing. Cell113, 673676. 10.1016/S0092-8674(03)00428-8

  • 2

    AmbrosV. (2004). The functions of animal microRNAs. Nature431, 350355. 10.1038/nature02871

  • 3

    BartelD. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell.116, 281297. 10.1016/S0092-8674(04)00045-5

  • 4

    BeggC. B.MazumdarM. (1994). Operating characteristics of a rank correlation test for publication bias. Biometrics50, 10881101. 10.2307/2533446

  • 5

    BoschS. L.van RooijenS. J.BokkerinkG. M.BraamH. J.DerikxL. A.PoortmansP.et al. (2017). Acute toxicity and surgical complications after preoperative (chemo)radiation therapy for rectal cancer in patients with inflammatory bowel disease. Radiother. Oncol.123, 147153. 10.1016/j.radonc.2017.02.009

  • 6

    DerSimonianR. (1996). Meta-analysis in the design and monitoring of clinical trials. Stat. Med.15, 1237-1248. 10.1002/(SICI)1097-0258(19960630)15:12&lt;1237::AID-SIM301&gt;3.0.CO;2-N

  • 7

    EggerM.Davey SmithG.SchneiderM.MinderC. (1997). Bias in meta-analysis detected by a simple, graphical test. BMJ315, 629634. 10.1136/bmj.315.7109.629

  • 8

    ErrenT. C.MorfeldP.FosterR. G.ReiterR. J.GrossJ. V.WestermannI. K. (2016). Sleep and cancer: synthesis of experimental data and meta-analyses of cancer incidence among some 1,500,000 study individuals in 13 countries. Chronobiol. Int.33, 325350. 10.3109/07420528.2016.1149486

  • 9

    GaoC.GaneshB. P.ShiZ.ShahR. R.FultzR.MajorA.et al. (2017). Gut microbe-mediated suppression of inflammation-associated colon carcinogenesis by luminal histamine production. Am. J. Pathol.187, 23232336. 10.1016/j.ajpath.2017.06.011

  • 10

    HsuC. Y.LiuP. H.HoS. Y.HuangY. H.LeeY. H.ChiouY. Y.et al. (2017). Impact of tumor burden on prognostic prediction for patients with terminal stage hepatocellular carcinoma: a nomogram study. PLoS ONE12:e0188031. 10.1371/journal.pone.0188031

  • 11

    HuangA. J.YuK. D.LiJ.FanL.ShaoZ. M. (2012). Polymorphism rs4919510:C>G in mature sequence of human microRNA-608 contributes to the risk of HER2-positive breast cancer but not other subtypes. PLoS ONE7:e35252. 10.1371/journal.pone.0035252

  • 12

    Huedo-MedinaT. B.Sanchez-MecaJ.Marin-MartinezF.BotellaJ. (2006). Assessing heterogeneity in meta-analysis: Q statistic or I2 index?Psychol. Methods11, 193206.10.1037/1082-989X.11.2.193

  • 13

    JarzynskiA.ZajacP.ZebrowskiR.BoguszewskaA.Polz-DacewiczM. (2017). Occurrence of BK virus and human papilloma virus in colorectal cancer. Ann. Agric. Environ. Med.24, 440445. 10.26444/aaem/74648

  • 14

    JayasekaraH.EnglishD. R.HaydonA.HodgeA. M.LynchB. M.RostyC.et al. (2018). Associations of alcohol intake, smoking, physical activity and obesity with survival following colorectal cancer diagnosis by stage, anatomic site and tumor molecular subtype. Int. J. Cancer142, 238250. 10.1002/ijc.31049

  • 15

    JiangJ.JiaZ. F.CaoD. H.WuY. H.SunZ. W.CaoX. Y. (2016). Association of the miR-146a rs2910164 polymorphism with gastric cancer susceptibility and prognosis. Future Oncol.12, 22152226. 10.2217/fon-2016-0224

  • 16

    JinjuvadiaR.SalamiA.LenhartA.JinjuvadiaK.LiangpunsakulS.SalgiaR. (2017). Hepatocellular carcinoma: a decade of hospitalizations and financial burden in the united states. Am. J. Med. Sci.354, 362369. 10.1016/j.amjms.2017.05.016

  • 17

    KupcinskasJ.BruzaiteI.JuzenasS.GyvyteU.JonaitisL.KiudelisG.et al. (2014a). Lack of association between miR-27a, miR-146a, miR-196a-2, miR-492 and miR-608 gene polymorphisms and colorectal cancer. Sci. Rep.4:5993. 10.1038/srep05993

  • 18

    KupcinskasJ.WexT.LinkA.LejaM.BruzaiteI.SteponaitieneR.et al. (2014b). Gene polymorphisms of micrornas in helicobacter pylori-induced high risk atrophic gastritis and gastric cancer. PLoS ONE9:e87467. 10.1371/journal.pone.0087467

  • 19

    Lagos-QuintanaM.RauhutR.LendeckelW.TuschlT. (2001). Identification of novel genes coding for small expressed RNAs. Science294, 853858. 10.1126/science.1064921

  • 20

    LandiD.GemignaniF.BaraleR.LandiS. (2008). A catalog of polymorphisms falling in microRNA-binding regions of cancer genes. DNA Cell Biol.27, 3543. 10.1089/dna.2007.0650

  • 21

    LeeR. C.FeinbaumR. L.AmbrosV. (1993). The, C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell75, 843854.

  • 22

    LiuH.ZhouY.LiuQ.XiaoG.WangB.LiW.et al. (2016). Association of miR-608 rs4919510 polymorphism and cancer risk: a meta-analysis based on 13,664 subjects. Oncotarget.8, 3702337031. 10.18632/oncotarget.9509

  • 23

    MantelN.HaenszelW. (1959). Statistical aspects of the analysis of data from retrospective studies of disease. J. Natl. Cancer Inst.22, 719748.

  • 24

    MartiniS.NairV.PatelS. R.EichingerF.NelsonR. G.WeilE. J.et al. (2013). From single nucleotide polymorphism to transcriptional mechanism: a model for FRMD3 in diabetic nephropathy. Diabetes62, 26052612. 10.2337/db12-1416

  • 25

    MoherD.LiberatiA.TetzlaffJ.AltmanD. G. (2009). Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. BMJ339:b2535. 10.1136/bmj.b2535

  • 26

    MotoyamaS.SatoY.SasakiT.WakitaA.KawakitaY.LiuJ.et al. (2017). Efficacy and safety of neoadjuvant chemoradiotherapy following esophagectomy with japanese-style extended 3-field lymphadenectomy for thoracic esophageal cancer. Anticancer Res.37, 58375843. 10.21873/anticanres.12027

  • 27

    NederlofN.SlamanA. E.van HagenP.van der GaastA.SlankamenacK.GisbertzS. S.et al. (2016). Using the comprehensive complication index to assess the impact of neoadjuvant chemoradiotherapy on complication severity after esophagectomy for cancer. Ann. Surg. Oncol.23, 39643971. 10.1245/s10434-016-5291-3

  • 28

    NiuY. M.DuX. Y.LuM. Y.XuQ. L.LuoJ.ShenM. (2015). Significant association between functional microRNA polymorphisms and head and neck cancer susceptibility: a comprehensive meta-analysis. Sci. Rep.5:12972. 10.1038/srep12972

  • 29

    NiuY.YuanH.ShenM.LiH.HuY.ChenN. (2014). Association between cyclooxygenase-2 gene polymorphisms and head and neck squamous cell carcinoma risk. J. Craniofac. Surg.25, 333337. 10.1097/SCS.0000000000000372

  • 30

    RyanB. M.McClaryA. C.ValeriN.RobinsonD.PaoneA.BowmanE. D.et al. (2012). rs4919510 in hsa-mir-608 is associated with outcome but not risk of colorectal cancer. PLoS ONE7:e36306. 10.1371/journal.pone.0036306

  • 31

    SiegelR. L.MillerK. D.JemalA. (2017). Cancer statistics, 2017. CA Cancer J. Clin.67, 730. 10.3322/caac.21387

  • 32

    TessitoreA.CicciarelliG.Del VecchioF.GaggianoA.VerzellaD.FischiettiM.et al. (2014). MicroRNAs in the DNA damage/repair network and cancer. Int. J. Genomics2014:820248. 10.1155/2014/820248

  • 33

    TorreL. A.BrayF.SiegelR. L.FerlayJ.Lortet-TieulentJ.JemalA. (2015). Global cancer statistics, 2012. CA Cancer J. Clin.65, 87108. 10.3322/caac.21262

  • 34

    ValeriN.GaspariniP.FabbriM.BraconiC.VeroneseA.LovatF.et al. (2010). Modulation of mismatch repair and genomic stability by miR-155. Proc. Natl. Acad. Sci. U.S.A.107, 69826987. 10.1073/pnas.1002472107

  • 35

    WangK.LiangQ.WeiL.ZhangW.ZhuP. (2016). MicroRNA-608 acts as a prognostic marker and inhibits the cell proliferation in hepatocellular carcinoma by macrophage migration inhibitory factor. Tumour Biol.37, 38233830. 10.1007/s13277-015-4213-5

  • 36

    WangR.ZhangJ.MaY.ChenL.GuoS.ZhangX.et al. (2014). Association study of miR149 rs2292832 and miR608 rs4919510 and the risk of hepatocellular carcinoma in a largescale population. Mol. Med. Rep.10, 27362744. 10.3892/mmr.2014.2536

  • 37

    WangZ.XueY.WangP.ZhuJ.MaJ. (2016). MiR-608 inhibits the migration and invasion of glioma stem cells by targeting macrophage migration inhibitory factor. Oncol. Rep.35, 27332742. 10.3892/or.2016.4652

  • 38

    WuS.YuanW.ShenY.LuX.LiY.TianT.et al. (2017). The miR-608 rs4919510 polymorphism may modify cancer susceptibility based on type. Tumour Biol.39:1010428317703819. 10.1177/1010428317703819

  • 39

    YangH.LiQ.NiuJ.LiB.JiangD.WanZ.et al. (2016). microRNA-342-5p and miR-608 inhibit colon cancer tumorigenesis by targeting NAA10. Oncotarget7, 27092720. 10.18632/oncotarget.6458

  • 40

    YeY.WangK. K.GuJ.YangH.LinJ.AjaniJ. A.et al. (2008). Genetic variations in microRNA-related genes are novel susceptibility loci for esophageal cancer risk. Cancer Prev. Res. (Phila).1, 460469. 10.1158/1940-6207.CAPR-08-0135

  • 41

    YingH. Q.PengH. X.HeB. S.PanY. Q.WangF.SunH. L.et al. (2016). MiR-608, pre-miR-124-1 and pre-miR26a-1 polymorphisms modify susceptibility and recurrence-free survival in surgically resected CRC individuals. Oncotarget7, 7586575873. 10.18632/oncotarget.12422

  • 42

    ZhangM. W. (2012). Study on the Associations of Lifestyle-Related Factors, Genetic Variants in miRNA Encoding Regions and miRNA Binding Sites With Colorectal Cancer Risk. Zhejiang: Zhejiang University.

  • 43

    ZhangP.WangJ.LuT.WangX.ZhengY.GuoS.et al. (2015). miR-449b rs10061133 and miR-4293 rs12220909 polymorphisms are associated with decreased esophageal squamous cell carcinoma in a chinese population. Tumour Biol.36, 87898795. 10.1007/s13277-015-3422-2

Summary

Keywords

microRNA-608, rs4919510, cancer, polymorphism, susceptibility

Citation

Li X-F, Song J-K, Cai J-W, Zeng Y-Q, Li M, Zhu J and Niu Y-M (2018) No Association Between MicroRNA-608 rs4919510 G>C Polymorphism and Digestive System Cancers Susceptibility: A Meta-Analysis Based on 10,836 Individuals. Front. Physiol. 9:705. doi: 10.3389/fphys.2018.00705

Received

05 March 2018

Accepted

22 May 2018

Published

07 June 2018

Volume

9 - 2018

Edited by

Brian James Morris, University of Sydney, Australia

Reviewed by

Hector A. Cabrera-Fuentes, Justus Liebig Universität Gießen, Germany; Andrey Turchinovich, Deutsches Krebsforschungszentrum (DKFZ), Germany

Updates

Copyright

*Correspondence: Yu-Ming Niu

This article was submitted to Integrative Physiology, a section of the journal Frontiers in Physiology

†These authors have contributed equally to this work.

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics