Single Nucleotide Polymorphisms with Cis-Regulatory Effects on Long Non-Coding Transcripts in Human Primary Monocytes

We applied genome-wide allele-specific expression analysis of monocytes from 188 samples. Monocytes were purified from white blood cells of healthy blood donors to detect cis-acting genetic variation that regulates the expression of long non-coding RNAs. We analysed 8929 regions harboring genes for potential long non-coding RNA that were retrieved from data from the ENCODE project. Of these regions, 60% were annotated as intergenic, which implies that they do not overlap with protein-coding genes. Focusing on the intergenic regions, and using stringent analysis of the allele-specific expression data, we detected robust cis-regulatory SNPs in 258 out of 489 informative intergenic regions included in the analysis. The cis-regulatory SNPs that were significantly associated with allele-specific expression of long non-coding RNAs were enriched to enhancer regions marked for active or bivalent, poised chromatin by histone modifications. Out of the lncRNA regions regulated by cis-acting regulatory SNPs, 20% (n = 52) were co-regulated with the closest protein coding gene. We compared the identified cis-regulatory SNPs with those in the catalog of SNPs identified by genome-wide association studies of human diseases and traits. This comparison identified 32 SNPs in loci from genome-wide association studies that displayed a strong association signal with allele-specific expression of non-coding RNAs in monocytes, with p-values ranging from 6.7×10⁻⁷ to 9.5×10⁻⁸⁹. The identified cis-regulatory SNPs are associated with diseases of the immune system, like multiple sclerosis and rheumatoid arthritis.     

Large-scale evaluation of candidate genes identifies associations between VEGF polymorphisms and bladder cancer risk

Common genetic variation could alter the risk for developing bladder cancer. We conducted a large-scale evaluation of single nucleotide polymorphisms (SNPs) in candidate genes for cancer to identify common variants that influence bladder cancer risk. An Illumina GoldenGate assay was used to genotype 1,433 SNPs within or near 386 genes in 1,086 cases and 1,033 controls in Spain. The most significant finding was in the 5′ UTR of VEGF (rs25648, p for likelihood ratio test, 2 degrees of freedom = 1 × 10⁻⁵). To further investigate the region, we analyzed 29 additional SNPs in VEGF, selected to saturate the promoter and 5′ UTR and to tag common genetic variation in this gene. Three additional SNPs in the promoter region (rs833052, rs1109324, and rs1547651) were associated with increased risk for bladder cancer: odds ratio (95% confidence interval): 2.52 (1.06–5.97), 2.74 (1.26–5.98), and 3.02 (1.36–6.63), respectively; and a polymorphism in intron 2 (rs3024994) was associated with reduced risk: 0.65 (0.46–0.91). Two of the promoter SNPs and the intron 2 SNP showed linkage disequilibrium with rs25648. Haplotype analyses revealed three blocks of linkage disequilibrium with significant associations for two blocks including the promoter and 5′ UTR (global p = 0.02 and 0.009, respectively). These findings are biologically plausible since VEGF is critical in angiogenesis, which is important for tumor growth, its elevated expression in bladder tumors correlates with tumor progression, and specific 5′ UTR haplotypes have been shown to influence promoter activity. Associations between bladder cancer risk and other genes in this report were not robust based on false discovery rate calculations. In conclusion, this large-scale evaluation of candidate cancer genes has identified common genetic variants in the regulatory regions of VEGF that could be associated with bladder cancer risk.     

MicroRNA Related Polymorphisms and Breast Cancer Risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88–0.96), rs1052532 (OR 0.97; 95% CI: 0.95–0.99), rs10719 (OR 0.97; 95% CI: 0.94–0.99), rs4687554 (OR 0.97; 95% CI: 0.95–0.99, and rs3134615 (OR 1.03; 95% CI: 1.01–1.05) located in the 3′ UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.     

Genetic Variants in Caveolin-1 and RhoA/ROCK1 Are Associated with Clear Cell Renal Cell Carcinoma Risk in a Chinese Population

Background

The RhoA/ROCK pathway and Caveolin-1 (Cav-1) participate in the process of tumorigenesis in numerous types of cancer. Up-regulation of RhoA/ROCK and Cav-1 expression is considered to be associated with the development and progression of clear cell renal cell carcinoma (ccRCC). We investigated the association between genetic variations of RhoA/ROCK and Cav-1 and the risk of ccRCC in the Chinese population.

Methods

Between May 2004 and March 2014, a total of 1,248 clear cell renal cell carcinoma cases and 1,440 cancer-free controls were enrolled in this hospital-based case-control study. Nine SNPs in RhoA/ROCK and Cav-1 were genotyped using the TaqMan assay.

Result

We found two SNPs (Cav-1 rs1049334 and ROCK1 rs35996865) were significantly associated with the increasing risk of ccRCC (P = 0.002 and P < 0.001 respectively). The analysis of combined risk alleles revealed that patients with 2–4 risk alleles showed a more remarkable growth of ccRCC risk than the patients with 0–1 risk alleles(OR = 1.66, 95%CI = 1.31–2.11, P < 0.001). Younger subjects (P = 0.001, OR = 1.83, 95%CI = 1.30–2.57), higher weight subjects (P = 0.001, OR = 1.76, 95%CI = 1.25–2.47), female subjects (P = 0.007, OR = 1.75, 95% CI = 1.17–2.62), nonsmokers (P < 0.001, OR = 1.67, 95%CI = 1.26–2.23), drinkers (P = 0.025, OR = 1.75, 95% CI = 1.07–2.85), subjects with hypertension (P = 0.025, OR = 1.75, 95% CI = 1.07–2.85) and diabetes (P = 0.026, OR = 4.31, 95% CI = 1.19–15.62) showed a stronger association between the combined risk alleles and the risk of ccRCC by using the stratification analysis. Furthermore, we observed higher Cav-1 mRNA levels in the presence of the rs1049334 A allele in normal renal tissues.

Conclusion

Our results indicate that the two SNPs (Cav-1 rs1049334 and ROCK1 rs35996865) and genotypes with a combination of 2–4 risk alleles were associated with the risk of ccRCC. The functional SNP rs1049334 may affect the risk of ccRCC by altering the expression of Cav-1 and the relevance between the risk effects and the functional impact of this polymorphism needs further validation.     

Genetic Determinants of Pelvic Organ Prolapse among African American and Hispanic Women in the Women’s Health Initiative

Current evidence suggests a multifactorial etiology to pelvic organ prolapse (POP), including genetic predisposition. We conducted a genome-wide association study of POP in African American (AA) and Hispanic (HP) women from the Women’s Health Initiative Hormone Therapy study. Cases were defined as any POP (grades 1–3) or moderate/severe POP (grades 2–3), while controls had grade 0 POP. We performed race-specific multiple logistic regression analyses between SNPs imputed to 1000 genomes in relation to POP (grade 0 vs 1–3; grade 0 vs 2–3) adjusting for age at diagnosis, body mass index, parity, and genetic ancestry. There were 1274 controls and 1427 cases of any POP and 317 cases of moderate/severe POP. Although none of the analyses reached genome-wide significance (p<5x10^-8), we noted variants in several loci that met p<10⁻⁶. In race-specific analysis of grade 0 vs 2–3, intronic SNPs in the CPE gene (rs28573326, OR:2.14; 95% CI 1.62–2.83; p = 1.0x10^-7) were associated with POP in AAs, and SNPs in the gene AL132709.5 (rs1950626, OR:2.96; 95% CI 1.96–4.48, p = 2.6x10^-7) were associated with POP in HPs. Inverse variance fixed-effect meta-analysis of the race-specific results showed suggestive signals for SNPs in the DPP6 gene (rs11243354, OR:1.36; p = 4.2x10^-7) in the grade 0 vs 1–3 analyses and for SNPs around PGBD5 (rs740494, OR:2.17; p = 8.6x10^-7) and SHC3 (rs2209875, OR:0.60; p = 9.3x10^-7) in the grade 0 vs 2–3 analyses. While we did not identify genome-wide significant findings, we document several SNPs reaching suggestive statistical significance. Further interrogation of POP in larger minority samples is warranted.     

Polymorphisms in the mTOR Gene and Risk of Sporadic Prostate Cancer in an Eastern Chinese Population

Background

The mTOR gene regulates cell growth by controlling mRNA translation, ribosome biogenesis, autophagy, and metabolism. Abnormally increased expression of mTOR was associated with carcinogenesis, and its functional single nucleotide polymorphisms (SNPs) may regulate the expression of mTOR and thus contribute to cancer risk.

Methodology/Principal Findings

In a hospital-based case-control study of 1004 prostate cancer (PCa) cases and 1051 cancer-free controls, we genotyped six potentially functional SNPs of mTOR (rs2536 T>C, rs1883965 G>A, rs1034528 G>C, rs17036508 T>C, rs3806317 A>G, and rs2295080 T>G) and assessed their associations with risk of PCa by using logistic regression analysis.

Conclusions/Significances

In the single-locus analysis, we found a significantly increased risk of PCa associated with mTOR rs2536 CT/CC and rs1034528 CG/CC genotypes [adjusted OR = 1.42 (1.13–1.78), P = 0.003 and 1.29 (1.07–1.55), P = 0.007), respectively], compared with their common homozygous genotypes, whereas mTOR rs2295080 GT/GG genotypes were associated with a decreased risk of PCa [adjusted OR = 0.76 (0.64–0.92), P = 0.003], compared with wild-type TT genotypes. In the combined analysis of the six SNPs, we found that individuals carrying two or more adverse genotypes had an increased risk of PCa [adjusted OR = 1.24 (1.04–1.47), P = 0.016], compared with individuals carrying less than two adverse genotypes. In the multiple dimension reduction analysis, body mass index (BMI) was the best one-factor model with the highest CVC (100%) and the lowest prediction error (42.7%) among all seven factors. The model including an interaction among BMI, rs17036508, and rs2536 was the best three-factor model with the highest CVC (100%) and the lowest prediction error of 41.9%. These findings suggested that mTOR SNPs may contribute to the risk of PCa in Eastern Chinese men, but the effect was weak and needs further validation by larger population-based studies.     

Predictive Value of Liver Enzymes and Inflammatory Biomarkers for the Severity of Liver Fibrosis Stage in HIV/HCV Co-Infected Patients

Objective

The aim of our study was to assess a possible association between plasma inflammatory biomarkers (CRP, IL-6, soluble CD14) and the extent of fibrosis or cirrhosis using a FibroScan® in HIV/HCV co-infected patients.

Methods

This cross-sectional study assessed 60 HIV/HCV co-infected patients who had paired plasma samples and FibroScan® values available. All included patients were controlled for HIV infection (HIV-1 RNA <50 copies/mL) and had detectable HCV RNA levels. Levels of three biomarkers were measured in all samples using commercial ELISA kits. Multivariate logistic regression models identified factors associated with the METAVIR stages of fibrosis (F0–F2 vs. F3–F4).

Results

In univariate logistic regression analyses, in addition to sCD14 (odds ratio [OR] = 3.23, 95% confidence interval [95%CI] = 1.30–7.97, P = 0.01), aspartate aminotransferase (AST), alanine aminotransferase, platelet counts, and CD4 cell counts were associated with the stage of liver fibrosis and, thus, were introduced into the model. However, only AST (OR = 1.06, 95%CI = 1.02–1.10, P = 0.0009) was independently associated with F3–F4 stage liver fibrosis.

Conclusions

In our study of HIV/HCV co-infected patients, sCD14 plasma level, a biomarker of monocyte activation, was not independently associated with the F3–F4 stage of liver fibrosis. We hypothesize that the higher levels of inflammation markers observed in HIV/HCV co-infected patients, compared to HCV mono-infected patients, prevent this association being observed within this population.     

Highly sensitive and specific microRNA expression profiling using BeadArray technology

We have developed a highly sensitive, specific and reproducible method for microRNA (miRNA) expression profiling, using the BeadArray™ technology. This method incorporates an enzyme-assisted specificity step, a solid-phase primer extension to distinguish between members of miRNA families. In addition, a universal PCR is used to amplify all targets prior to array hybridization. Currently, assay probes are designed to simultaneously analyse 735 well-annotated human miRNAs. Using this method, highly reproducible miRNA expression profiles were generated with 100–200 ng total RNA input. Furthermore, very similar expression profiles were obtained with total RNA and enriched small RNA species (R² ≥ 0.97). The method has a 3.5–4 log (10⁵–10⁹ molecules) dynamic range and is able to detect 1.2- to 1.3-fold-differences between samples. Expression profiles generated by this method are highly comparable to those obtained with RT–PCR (R² = 0.85–0.90) and direct sequencing (R = 0.87–0.89). This method, in conjunction with the 96-sample array matrix should prove useful for high-throughput expression profiling of miRNAs in large numbers of tissue samples.     

Connecting the Dots: Potential of Data Integration to Identify Regulatory SNPs in Late-Onset Alzheimer's Disease GWAS Findings

Late-onset Alzheimer''s disease (LOAD) is a multifactorial disorder with over twenty loci associated with disease risk. Given the number of genome-wide significant variants that fall outside of coding regions, it is possible that some of these variants alter some function of gene expression rather than tagging coding variants that alter protein structure and/or function. RegulomeDB is a database that annotates regulatory functions of genetic variants. In this study, we utilized RegulomeDB to investigate potential regulatory functions of lead single nucleotide polymorphisms (SNPs) identified in five genome-wide association studies (GWAS) of risk and age-at onset (AAO) of LOAD, as well as SNPs in LD (r²≥0.80) with the lead GWAS SNPs. Of a total 614 SNPs examined, 394 returned RegulomeDB scores of 1–6. Of those 394 variants, 34 showed strong evidence of regulatory function (RegulomeDB score <3), and only 3 of them were genome-wide significant SNPs (ZCWPW1/rs1476679, CLU/rs1532278 and ABCA7/rs3764650). This study further supports the assumption that some of the non-coding GWAS SNPs are true associations rather than tagged associations and demonstrates the application of RegulomeDB to GWAS data.     

Polymorphisms in a Putative Enhancer at the 10q21.2 Breast Cancer Risk Locus Regulate NRBF2 Expression

Genome-wide association studies have identified SNPs near ZNF365 at 10q21.2 that are associated with both breast cancer risk and mammographic density. To identify the most likely causal SNPs, we fine mapped the association signal by genotyping 428 SNPs across the region in 89,050 European and 12,893 Asian case and control subjects from the Breast Cancer Association Consortium. We identified four independent sets of correlated, highly trait-associated variants (iCHAVs), three of which were located within ZNF365. The most strongly risk-associated SNP, rs10995201 in iCHAV1, showed clear evidence of association with both estrogen receptor (ER)-positive (OR = 0.85 [0.82–0.88]) and ER-negative (OR = 0.87 [0.82–0.91]) disease, and was also the SNP most strongly associated with percent mammographic density. iCHAV2 (lead SNP, chr10: 64,258,684:D) and iCHAV3 (lead SNP, rs7922449) were also associated with ER-positive (OR = 0.93 [0.91–0.95] and OR = 1.06 [1.03–1.09]) and ER-negative (OR = 0.95 [0.91–0.98] and OR = 1.08 [1.04–1.13]) disease. There was weaker evidence for iCHAV4, located 5′ of ADO, associated only with ER-positive breast cancer (OR = 0.93 [0.90–0.96]). We found 12, 17, 18, and 2 candidate causal SNPs for breast cancer in iCHAVs 1–4, respectively. Chromosome conformation capture analysis showed that iCHAV2 interacts with the ZNF365 and NRBF2 (more than 600 kb away) promoters in normal and cancerous breast epithelial cells. Luciferase assays did not identify SNPs that affect transactivation of ZNF365, but identified a protective haplotype in iCHAV2, associated with silencing of the NRBF2 promoter, implicating this gene in the etiology of breast cancer.     

Diversity and Community Composition of Methanogenic Archaea in the Rumen of Scottish Upland Sheep Assessed by Different Methods

Ruminal archaeomes of two mature sheep grazing in the Scottish uplands were analysed by different sequencing and analysis methods in order to compare the apparent archaeal communities. All methods revealed that the majority of methanogens belonged to the Methanobacteriales order containing the Methanobrevibacter, Methanosphaera and Methanobacteria genera. Sanger sequenced 1.3 kb 16S rRNA gene amplicons identified the main species of Methanobrevibacter present to be a SGMT Clade member Mbb. millerae (≥91% of OTUs); Methanosphaera comprised the remainder of the OTUs. The primers did not amplify ruminal Thermoplasmatales-related 16S rRNA genes. Illumina sequenced V6–V8 16S rRNA gene amplicons identified similar Methanobrevibacter spp. and Methanosphaera clades and also identified the Thermoplasmatales-related order as 13% of total archaea. Unusually, both methods concluded that Mbb. ruminantium and relatives from the same clade (RO) were almost absent. Sequences mapping to rumen 16S rRNA and mcrA gene references were extracted from Illumina metagenome data. Mapping of the metagenome data to16S rRNA gene references produced taxonomic identification to Order level including 2–3% Thermoplasmatales, but was unable to discriminate to species level. Mapping of the metagenome data to mcrA gene references resolved 69% to unclassified Methanobacteriales. Only 30% of sequences were assigned to species level clades: of the sequences assigned to Methanobrevibacter, most mapped to SGMT (16%) and RO (10%) clades. The Sanger 16S amplicon and Illumina metagenome mcrA analyses showed similar species richness (Chao1 Index 19–35), while Illumina metagenome and amplicon 16S rRNA analysis gave lower richness estimates (10–18). The values of the Shannon Index were low in all methods, indicating low richness and uneven species distribution. Thus, although much information may be extracted from the other methods, Illumina amplicon sequencing of the V6–V8 16S rRNA gene would be the method of choice for studying rumen archaeal communities.     

Epigenome-Wide DNA Methylation in Hearing Ability: New Mechanisms for an Old Problem

Epigenetic regulation of gene expression has been shown to change over time and may be associated with environmental exposures in common complex traits. Age-related hearing impairment is a complex disorder, known to be heritable, with heritability estimates of 57–70%. Epigenetic regulation might explain the observed difference in age of onset and magnitude of hearing impairment with age. Epigenetic epidemiology studies using unrelated samples can be limited in their ability to detect small effects, and recent epigenetic findings in twins underscore the power of this well matched study design. We investigated the association between venous blood DNA methylation epigenome-wide and hearing ability. Pure-tone audiometry (PTA) and Illumina HumanMethylation array data were obtained from female twin volunteers enrolled in the TwinsUK register. Two study groups were explored: first, an epigenome-wide association scan (EWAS) was performed in a discovery sample (n = 115 subjects, age range: 47–83 years, Illumina 27 k array), then replication of the top ten associated probes from the discovery EWAS was attempted in a second unrelated sample (n = 203, age range: 41–86 years, Illumina 450 k array). Finally, a set of monozygotic (MZ) twin pairs (n = 21 pairs) within the discovery sample (Illumina 27 k array) was investigated in more detail in an MZ discordance analysis. Hearing ability was strongly associated with DNA methylation levels in the promoter regions of several genes, including TCF25 (cg01161216, p = 6.6×10⁻⁶), FGFR1 (cg15791248, p = 5.7×10⁻⁵) and POLE (cg18877514, p = 6.3×10⁻⁵). Replication of these results in a second sample confirmed the presence of differential methylation at TCF25 (p(replication) = 6×10⁻⁵) and POLE (p(replication) = 0.016). In the MZ discordance analysis, twins'' intrapair difference in hearing ability correlated with DNA methylation differences at ACP6 (cg01377755, r = −0.75, p = 1.2×10⁻⁴) and MEF2D (cg08156349, r = −0.75, p = 1.4×10⁻⁴). Examination of gene expression in skin, suggests an influence of differential methylation on expression, which may account for the variation in hearing ability with age.     

Fine Mapping for Weaver Syndrome in Brown Swiss Cattle and the Identification of 41 Concordant Mutations across NRCAM,PNPLA8 and CTTNBP2

Bovine Progressive Degenerative Myeloencephalopathy (Weaver Syndrome) is a recessive neurological disease that has been observed in the Brown Swiss cattle breed since the 1970’s in North America and Europe. Bilateral hind leg weakness and ataxia appear in afflicted animals at 6 to 18 months of age, and slowly progresses to total loss of hind limb control by 3 to 4 years of age. While Weaver has previously been mapped to Bos taurus autosome (BTA) 4∶46–56 Mb and a diagnostic test based on the 6 microsatellite (MS) markers is commercially available, neither the causative gene nor mutation has been identified; therefore misdiagnosis can occur due to recombination between the diagnostic MS markers and the causative mutation. Analysis of 34,980 BTA 4 SNPs genotypes derived from the Illumina BovineHD assay for 20 Brown Swiss Weaver carriers and 49 homozygous normal bulls refined the Weaver locus to 48–53 Mb. Genotyping of 153 SNPs, identified from whole genome sequencing of 10 normal and 10 carrier animals, across a validation set of 841 animals resulted in the identification of 41 diagnostic SNPs that were concordant with the disease. Except for one intergenic SNP all are associated with genes expressed in nervous tissues: 37 distal to NRCAM, one non-synonymous (serine to asparagine) in PNPLA8, one synonymous and one non-synonymous (lysine to glutamic acid) in CTTNBP2. Haplotype and imputation analyses of 7,458 Brown Swiss animals with Illumina BovineSNP50 data and the 41 diagnostic SNPs resulted in the identification of only one haplotype concordant with the Weaver phenotype. Use of this haplotype and the diagnostic SNPs more accurately identifies Weaver carriers in both Brown Swiss purebred and influenced herds.     

Genome-Wide Association Study Reveals Genetic Architecture of Eating Behavior in Pigs and Its Implications for Humans Obesity by Comparative Mapping

This study was aimed at identifying genomic regions controlling feeding behavior in Danish Duroc boars and its potential implications for eating behavior in humans. Data regarding individual daily feed intake (DFI), total daily time spent in feeder (TPD), number of daily visits to feeder (NVD), average duration of each visit (TPV), mean feed intake per visit (FPV) and mean feed intake rate (FR) were available for 1130 boars. All boars were genotyped using the Illumina Porcine SNP60 BeadChip. The association analyses were performed using the GenABEL package in the R program. Sixteen SNPs were found to have moderate genome-wide significance (p<5E-05) and 76 SNPs had suggestive (p<5E-04) association with feeding behavior traits. MSI2 gene on chromosome (SSC) 14 was very strongly associated with NVD. Thirty-six SNPs were located in genome regions where QTLs have previously been reported for behavior and/or feed intake traits in pigs. The regions: 64–65 Mb on SSC 1, 124–130 Mb on SSC 8, 63–68 Mb on SSC 11, 32–39 Mb and 59–60 Mb on SSC 12 harbored several signifcant SNPs. Synapse genes (GABRR2, PPP1R9B, SYT1, GABRR1, CADPS2, DLGAP2 and GOPC), dephosphorylation genes (PPM1E, DAPP1, PTPN18, PTPRZ1, PTPN4, MTMR4 and RNGTT) and positive regulation of peptide secretion genes (GHRH, NNAT and TCF7L2) were highly significantly associated with feeding behavior traits. This is the first GWAS to identify genetic variants and biological mechanisms for eating behavior in pigs and these results are important for genetic improvement of pig feed efficiency. We have also conducted pig-human comparative gene mapping to reveal key genomic regions and/or genes on the human genome that may influence eating behavior in human beings and consequently affect the development of obesity and metabolic syndrome. This is the first translational genomics study of its kind to report potential candidate genes for eating behavior in humans.     

Replication and Predictive Value of SNPs Associated with Melanoma and Pigmentation Traits in a Southern European Case-Control Study

Background

Genetic association studies have revealed numerous polymorphisms conferring susceptibility to melanoma. We aimed to replicate previously discovered melanoma-associated single-nucleotide polymorphisms (SNPs) in a Greek case-control population, and examine their predictive value.

Methods

Based on a field synopsis of genetic variants of melanoma (MelGene), we genotyped 284 patients and 284 controls at 34 melanoma-associated SNPs of which 19 derived from GWAS. We tested each one of the 33 SNPs passing quality control for association with melanoma both with and without accounting for the presence of well-established phenotypic risk factors. We compared the risk allele frequencies between the Greek population and the HapMap CEU sample. Finally, we evaluated the predictive ability of the replicated SNPs.

Results

Risk allele frequencies were significantly lower compared to the HapMap CEU for eight SNPs (rs16891982 – SLC45A2, rs12203592 – IRF4, rs258322 – CDK10, rs1805007 – MC1R, rs1805008 - MC1R, rs910873 - PIGU, rs17305573- PIGU, and rs1885120 - MTAP) and higher for one SNP (rs6001027 – PLA2G6) indicating a different profile of genetic susceptibility in the studied population. Previously identified effect estimates modestly correlated with those found in our population (r = 0.72, P<0.0001). The strongest associations were observed for rs401681-T in CLPTM1L (odds ratio [OR] 1.60, 95% CI 1.22–2.10; P = 0.001), rs16891982-C in SCL45A2 (OR 0.51, 95% CI 0.34–0.76; P = 0.001), and rs1805007-T in MC1R (OR 4.38, 95% CI 2.03–9.43; P = 2×10⁻⁵). Nominally statistically significant associations were seen also for another 5 variants (rs258322-T in CDK10, rs1805005-T in MC1R, rs1885120-C in MYH7B, rs2218220-T in MTAP and rs4911442-G in the ASIP region). The addition of all SNPs with nominal significance to a clinical non-genetic model did not substantially improve melanoma risk prediction (AUC for clinical model 83.3% versus 83.9%, p = 0.66).

Conclusion

Overall, our study has validated genetic variants that are likely to contribute to melanoma susceptibility in the Greek population.