Novel and lethal case of cardiac involvement in DNM1L mitochondrial encephalopathy.

Pathogenic DNM1L mutations cause a mitochondrial disorder with a highly variable clinical phenotype characterized by developmental delay, hypotonia, seizures, microcephaly, poor feeding, ocular abnormalities, and dysarthria. We report the case of an 8-month-old female with autosomal dominant, de novo DNM1L c. 1228G>A (p. E410K) mutation and mitochondrial disorder, septo-optic dysplasia, hypotonia, developmental delay, elevated blood lactate, and severe mitochondrial cardiomyopathy leading to nonischemic congestive heart failure and cardiogenic shock resulting in death. This case suggests that cardiac involvement, previously undescribed, can be a clinically important feature of this syndrome and should be screened for at diagnosis.

Genes known to escape X chromosome inactivation predict co-morbid chronic musculoskeletal pain and posttraumatic stress symptom development in women following trauma exposure.

Co-morbid chronic musculoskeletal pain (CMSP) and posttraumatic stress symptoms (PTSS) are frequent sequelae of motor vehicle collision, are associated with greater disability than either outcome alone, and are more prevalent in women than men. In the current study we assessed for evidence that gene transcripts originating from the X chromosome contribute to sex differences in vulnerability to CMSP and PTSS after motor vehicle collision. Nested samples were drawn from a longitudinal study of African American individuals, and CMSP (0-10 numeric rating scale) and PTSS (impact of events scale, revised) outcomes were assessed 6 months following motor vehicle collision. Blood RNA were sequenced (n = 101) and the relationship between X chromosome mRNA expression levels and co-morbid CMSP and PTSS outcomes was evaluated using logistic regression analyses. A disproportionate number of peritraumatic X chromosome mRNA predicting CMSP and PTSS in women were genes previously found to escape X chromosome inactivation (11/40, z = -2.9, p = .004). Secondary analyses assessing gene ontology relationships between these genes identified an enrichment in genes known to influence neuronal plasticity. Further, the relationship of expression of two critical regulators of X chromosome inactivation, X-inactive specific transcript (XIST) and Yin Yang 1 (YY1), was different in women developing CMSP and PTSS. Together, these data suggest that X chromosome genes that escape inactivation may contribute to sex differences in vulnerability to CMSP and PTSS after motor vehicle collision.

Characterization of large structural genetic mosaicism in human autosomes.

Analyses of genome-wide association study (GWAS) data have revealed that detectable genetic mosaicism involving large (>2 Mb) structural autosomal alterations occurs in a fraction of individuals. We present results for a set of 24,849 genotyped individuals (total GWAS set II [TGSII]) in whom 341 large autosomal abnormalities were observed in 168 (0.68%) individuals. Merging data from the new TGSII set with data from two prior reports (the Gene-Environment Association Studies and the total GWAS set I) generated a large dataset of 127,179 individuals; we then conducted a meta-analysis to investigate the patterns of detectable autosomal mosaicism (n = 1,315 events in 925 [0.73%] individuals). Restricting to events >2 Mb in size, we observed an increase in event frequency as event size decreased. The combined results underscore that the rate of detectable mosaicism increases with age (p value = 5.5 × 10(-31)) and is higher in men (p value = 0.002) but lower in participants of African ancestry (p value = 0.003). In a subset of 47 individuals from whom serial samples were collected up to 6 years apart, complex changes were noted over time and showed an overall increase in the proportion of mosaic cells as age increased. Our large combined sample allowed for a unique ability to characterize detectable genetic mosaicism involving large structural events and strengthens the emerging evidence of non-random erosion of the genome in the aging population.

Integration of multiethnic fine-mapping and genomic annotation to prioritize candidate functional SNPs at prostate cancer susceptibility regions.

Interpretation of biological mechanisms underlying genetic risk associations for prostate cancer is complicated by the relatively large number of risk variants (n = 100) and the thousands of surrogate SNPs in linkage disequilibrium. Here, we combined three distinct approaches: multiethnic fine-mapping, putative functional annotation (based upon epigenetic data and genome-encoded features), and expression quantitative trait loci (eQTL) analyses, in an attempt to reduce this complexity. We examined 67 risk regions using genotyping and imputation-based fine-mapping in populations of European (cases/controls: 8600/6946), African (cases/controls: 5327/5136), Japanese (cases/controls: 2563/4391) and Latino (cases/controls: 1034/1046) ancestry. Markers at 55 regions passed a region-specific significance threshold (P-value cutoff range: 3.9 × 10(-4)-5.6 × 10(-3)) and in 30 regions we identified markers that were more significantly associated with risk than the previously reported variants in the multiethnic sample. Novel secondary signals (P < 5.0 × 10(-6)) were also detected in two regions (rs13062436/3q21 and rs17181170/3p12). Among 666 variants in the 55 regions with P-values within one order of magnitude of the most-associated marker, 193 variants (29%) in 48 regions overlapped with epigenetic or other putative functional marks. In 11 of the 55 regions, cis-eQTLs were detected with nearby genes. For 12 of the 55 regions (22%), the most significant region-specific, prostate-cancer associated variant represented the strongest candidate functional variant based on our annotations; the number of regions increased to 20 (36%) and 27 (49%) when examining the 2 and 3 most significantly associated variants in each region, respectively. These results have prioritized subsets of candidate variants for downstream functional evaluation.

Four Susceptibility Loci for Gallstone Disease Identified in a Meta-analysis of Genome-Wide Association Studies.

BACKGROUND & AIMS: A genome-wide association study (GWAS) of 280 cases identified the hepatic cholesterol transporter ABCG8 as a locus associated with risk for gallstone disease, but findings have not been reported from any other GWAS of this phenotype. We performed a large-scale, meta-analysis of GWASs of individuals of European ancestry with available prior genotype data, to identify additional genetic risk factors for gallstone disease. METHODS: We obtained per-allele odds ratio (OR) and standard error estimates using age- and sex-adjusted logistic regression models within each of the 10 discovery studies (8720 cases and 55,152 controls). We performed an inverse variance weighted, fixed-effects meta-analysis of study-specific estimates to identify single-nucleotide polymorphisms that were associated independently with gallstone disease. Associations were replicated in 6489 cases and 62,797 controls. RESULTS: We observed independent associations for 2 single-nucleotide polymorphisms at the ABCG8 locus: rs11887534 (OR, 1.69; 95% confidence interval [CI], 1.54-1.86; P = 2.44 × 10(-60)) and rs4245791 (OR, 1.27; P = 1.90 × 10(-34)). We also identified and/or replicated associations for rs9843304 in TM4SF4 (OR, 1.12; 95% CI, 1.08-1.16; P = 6.09 × 10(-11)), rs2547231 in SULT2A1 (encodes a sulfoconjugation enzyme that acts on hydroxysteroids and cholesterol-derived sterol bile acids) (OR, 1.17; 95% CI, 1.12-1.21; P = 2.24 × 10(-10)), rs1260326 in glucokinase regulatory protein (OR, 1.12; 95% CI, 1.07-1.17; P = 2.55 × 10(-10)), and rs6471717 near CYP7A1 (encodes an enzyme that catalyzes conversion of cholesterol to primary bile acids) (OR, 1.11; 95% CI, 1.08-1.15; P = 8.84 × 10(-9)). Among individuals of African American and Hispanic American ancestry, rs11887534 and rs4245791 were associated positively with gallstone disease risk, whereas the association for the rs1260326 variant was inverse. CONCLUSIONS: In this large-scale GWAS of gallstone disease, we identified 4 loci in genes that have putative functions in cholesterol metabolism and transport, and sulfonylation of bile acids or hydroxysteroids.

Expression QTL-based analyses reveal candidate causal genes and loci across five tumor types.

The majority of trait-associated loci discovered through genome-wide association studies are located outside of known protein coding regions. Consequently, it is difficult to ascertain the mechanism underlying these variants and to pinpoint the causal alleles. Expression quantitative trait loci (eQTLs) provide an organizing principle to address both of these issues. eQTLs are genetic loci that correlate with RNA transcript levels. Large-scale data sets such as the Cancer Genome Atlas (TCGA) provide an ideal opportunity to systematically evaluate eQTLs as they have generated multiple data types on hundreds of samples. We evaluated the determinants of gene expression (germline variants and somatic copy number and methylation) and performed cis-eQTL analyses for mRNA expression and miRNA expression in five tumor types (breast, colon, kidney, lung and prostate). We next tested 149 known cancer risk loci for eQTL effects, and observed that 42 (28.2%) were significantly associated with at least one transcript. Lastly, we described a fine-mapping strategy for these 42 eQTL target-gene associations based on an integrated strategy that combines the eQTL level of significance and the regulatory potential as measured by DNaseI hypersensitivity. For each of the risk loci, our analyses suggested 1 to 81 candidate causal variants that may be prioritized for downstream functional analysis. In summary, our study provided a comprehensive landscape of the genetic determinants of gene expression in different tumor types and ranked the genes and loci for further functional assessment of known cancer risk loci.

Genome-wide association study identifies multiple loci associated with bladder cancer risk.

Candidate gene and genome-wide association studies (GWAS) have identified 11 independent susceptibility loci associated with bladder cancer risk. To discover additional risk variants, we conducted a new GWAS of 2422 bladder cancer cases and 5751 controls, followed by a meta-analysis with two independently published bladder cancer GWAS, resulting in a combined analysis of 6911 cases and 11 814 controls of European descent. TaqMan genotyping of 13 promising single nucleotide polymorphisms with P < 1 × 10(-5) was pursued in a follow-up set of 801 cases and 1307 controls. Two new loci achieved genome-wide statistical significance: rs10936599 on 3q26.2 (P = 4.53 × 10(-9)) and rs907611 on 11p15.5 (P = 4.11 × 10(-8)). Two notable loci were also identified that approached genome-wide statistical significance: rs6104690 on 20p12.2 (P = 7.13 × 10(-7)) and rs4510656 on 6p22.3 (P = 6.98 × 10(-7)); these require further studies for confirmation. In conclusion, our study has identified new susceptibility alleles for bladder cancer risk that require fine-mapping and laboratory investigation, which could further understanding into the biological underpinnings of bladder carcinogenesis.

Functional variants at the 11q13 risk locus for breast cancer regulate cyclin D1 expression through long-range enhancers.

Analysis of 4,405 variants in 89,050 European subjects from 41 case-control studies identified three independent association signals for estrogen-receptor-positive tumors at 11q13. The strongest signal maps to a transcriptional enhancer element in which the G allele of the best candidate causative variant rs554219 increases risk of breast cancer, reduces both binding of ELK4 transcription factor and luciferase activity in reporter assays, and may be associated with low cyclin D1 protein levels in tumors. Another candidate variant, rs78540526, lies in the same enhancer element. Risk association signal 2, rs75915166, creates a GATA3 binding site within a silencer element. Chromatin conformation studies demonstrate that these enhancer and silencer elements interact with each other and with their likely target gene, CCND1.

Identification of a common variant with potential pleiotropic effect on risk of inflammatory bowel disease and colorectal cancer.

Although genome-wide association studies (GWAS) have separately identified many genetic susceptibility loci for ulcerative colitis (UC), Crohn's disease (CD) and colorectal cancer (CRC), there has been no large-scale examination for pleiotropy, or shared genetic susceptibility, for these conditions. We used logistic regression modeling to examine the associations of 181 UC and CD susceptibility variants previously identified by GWAS with risk of CRC using data from the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. We also examined associations of significant variants with clinical and molecular characteristics in a subset of the studies. Among 11794 CRC cases and 14190 controls, rs11676348, the susceptibility single nucleotide polymorphism (SNP) for UC, was significantly associated with reduced risk of CRC (P = 7E-05). The multivariate-adjusted odds ratio of CRC with each copy of the T allele was 0.93 (95% CI 0.89-0.96). The association of the SNP with risk of CRC differed according to mucinous histological features (P heterogeneity = 0.008). In addition, the (T) allele was associated with lower risk of tumors with Crohn's-like reaction but not tumors without such immune infiltrate (P heterogeneity = 0.02) and microsatellite instability-high (MSI-high) but not microsatellite stable or MSI-low tumors (P heterogeneity = 0.03). The minor allele (T) in SNP rs11676348, located downstream from CXCR2 that has been implicated in CRC progression, is associated with a lower risk of CRC, particularly tumors with a mucinous component, Crohn's-like reaction and MSI-high. Our findings offer the promise of risk stratification of inflammatory bowel disease patients for complications such as CRC.

Common germline polymorphisms associated with breast cancer-specific survival.

INTRODUCTION: Previous studies have identified common germline variants nominally associated with breast cancer survival. These associations have not been widely replicated in further studies. The purpose of this study was to evaluate the association of previously reported SNPs with breast cancer-specific survival using data from a pooled analysis of eight breast cancer survival genome-wide association studies (GWAS) from the Breast Cancer Association Consortium. METHODS: A literature review was conducted of all previously published associations between common germline variants and three survival outcomes: breast cancer-specific survival, overall survival and disease-free survival. All associations that reached the nominal significance level of P value <0.05 were included. Single nucleotide polymorphisms that had been previously reported as nominally associated with at least one survival outcome were evaluated in the pooled analysis of over 37,000 breast cancer cases for association with breast cancer-specific survival. Previous associations were evaluated using a one-sided test based on the reported direction of effect. RESULTS: Fifty-six variants from 45 previous publications were evaluated in the meta-analysis. Fifty-four of these were evaluated in the full set of 37,954 breast cancer cases with 2,900 events and the two additional variants were evaluated in a reduced sample size of 30,000 samples in order to ensure independence from the previously published studies. Five variants reached nominal significance (P <0.05) in the pooled GWAS data compared to 2.8 expected under the null hypothesis. Seven additional variants were associated (P <0.05) with ER-positive disease. CONCLUSIONS: Although no variants reached genome-wide significance (P <5 x 10(-8)), these results suggest that there is some evidence of association between candidate common germline variants and breast cancer prognosis. Larger studies from multinational collaborations are necessary to increase the power to detect associations, between common variants and prognosis, at more stringent significance levels.

A genome-wide association study of early menopause and the combined impact of identified variants.

Early menopause (EM) affects up to 10% of the female population, reducing reproductive lifespan considerably. Currently, it constitutes the leading cause of infertility in the western world, affecting mainly those women who postpone their first pregnancy beyond the age of 30 years. The genetic aetiology of EM is largely unknown in the majority of cases. We have undertaken a meta-analysis of genome-wide association studies (GWASs) in 3493 EM cases and 13 598 controls from 10 independent studies. No novel genetic variants were discovered, but the 17 variants previously associated with normal age at natural menopause as a quantitative trait (QT) were also associated with EM and primary ovarian insufficiency (POI). Thus, EM has a genetic aetiology which overlaps variation in normal age at menopause and is at least partly explained by the additive effects of the same polygenic variants. The combined effect of the common variants captured by the single nucleotide polymorphism arrays was estimated to account for ∼30% of the variance in EM. The association between the combined 17 variants and the risk of EM was greater than the best validated non-genetic risk factor, smoking.

Autologous breast reconstruction: preoperative magnetic resonance angiography for perforator flap vessel mapping.

BACKGROUND: Selection of a vascular pedicle for autologous breast reconstruction is time consuming and depends on visual evaluation during the surgery. Preoperative imaging of donor site for mapping the perforator artery anatomy greatly improves the efficiency of perforator selection and significantly reduces the operative time. In this article, we present our experience with magnetic resonance angiography (MRA) for perforator vessel mapping including MRA technique and interpretation. METHODS: We have performed over 400 MRA examinations from August 2008 to August 2013 at our institution for preoperative imaging of donor site for mapping the perforator vessel anatomy. Using our optimized imaging protocol with blood pool magnetic resonance imaging contrast agents, multiple donor sites can be imaged in a single MRA examination. Following imaging using the postprocessing and reporting tool, we estimated incidence of commonly used perforators for autologous breast reconstruction. RESULTS: In our practice, anterior abdominal wall tissue is the most commonly used donor site for perforator flap breast reconstruction and deep inferior epigastric artery perforators are the most commonly used vascular pedicle. A thigh flap, based on the profunda femoral artery perforator has become the second most used flap at our institution. In addition, MRA imaging also showed evidence of metastatic disease in 4% of our patient subset. CONCLUSION: Our MRA technique allows the surgeons to confidently assess multiple donor sites for the best perforator and flap design. In conclusion, a well-performed MRA with specific postprocessing provides an accurate method for mapping perforator vessel, at the same time avoiding ionizing radiation.

Genome-wide interaction study of smoking and bladder cancer risk.

Bladder cancer is a complex disease with known environmental and genetic risk factors. We performed a genome-wide interaction study (GWAS) of smoking and bladder cancer risk based on primary scan data from 3002 cases and 4411 controls from the National Cancer Institute Bladder Cancer GWAS. Alternative methods were used to evaluate both additive and multiplicative interactions between individual single nucleotide polymorphisms (SNPs) and smoking exposure. SNPs with interaction P values < 5 × 10(-) (5) were evaluated further in an independent dataset of 2422 bladder cancer cases and 5751 controls. We identified 10 SNPs that showed association in a consistent manner with the initial dataset and in the combined dataset, providing evidence of interaction with tobacco use. Further, two of these novel SNPs showed strong evidence of association with bladder cancer in tobacco use subgroups that approached genome-wide significance. Specifically, rs1711973 (FOXF2) on 6p25.3 was a susceptibility SNP for never smokers [combined odds ratio (OR) = 1.34, 95% confidence interval (CI) = 1.20-1.50, P value = 5.18 × 10(-) (7)]; and rs12216499 (RSPH3-TAGAP-EZR) on 6q25.3 was a susceptibility SNP for ever smokers (combined OR = 0.75, 95% CI = 0.67-0.84, P value = 6.35 × 10(-) (7)). In our analysis of smoking and bladder cancer, the tests for multiplicative interaction seemed to more commonly identify susceptibility loci with associations in never smokers, whereas the additive interaction analysis identified more loci with associations among smokers-including the known smoking and NAT2 acetylation interaction. Our findings provide additional evidence of gene-environment interactions for tobacco and bladder cancer.

Exploring the genetic architecture of circulating 25-hydroxyvitamin D.

The primary circulating form of vitamin D is 25-hydroxy vitamin D (25(OH)D), a modifiable trait linked with a growing number of chronic diseases. In addition to environmental determinants of 25(OH)D, including dietary sources and skin ultraviolet B (UVB) exposure, twin- and family-based studies suggest that genetics contribute substantially to vitamin D variability with heritability estimates ranging from 43% to 80%. Genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) located in four gene regions associated with 25(OH)D. These SNPs collectively explain only a fraction of the heritability in 25(OH)D estimated by twin- and family-based studies. Using 25(OH)D concentrations and GWAS data on 5,575 subjects drawn from five cohorts, we hypothesized that genome-wide data, in the form of (1) a polygenic score comprised of hundreds or thousands of SNPs that do not individually reach GWAS significance, or (2) a linear mixed model for genome-wide complex trait analysis, would explain variance in measured circulating 25(OH)D beyond that explained by known genome-wide significant 25(OH)D-associated SNPs. GWAS identified SNPs explained 5.2% of the variation in circulating 25(OH)D in these samples and there was little evidence additional markers significantly improved predictive ability. On average, a polygenic score comprised of GWAS-identified SNPs explained a larger proportion of variation in circulating 25(OH)D than scores comprised of thousands of SNPs that were on average, nonsignificant. Employing a linear mixed model for genome-wide complex trait analysis explained little additional variability (range 0-22%). The absence of a significant polygenic effect in this relatively large sample suggests an oligogenetic architecture for 25(OH)D.

A meta-analysis of genome-wide association studies of breast cancer identifies two novel susceptibility loci at 6q14 and 20q11.

Genome-wide association studies (GWAS) of breast cancer defined by hormone receptor status have revealed loci contributing to susceptibility of estrogen receptor (ER)-negative subtypes. To identify additional genetic variants for ER-negative breast cancer, we conducted the largest meta-analysis of ER-negative disease to date, comprising 4754 ER-negative cases and 31 663 controls from three GWAS: NCI Breast and Prostate Cancer Cohort Consortium (BPC3) (2188 ER-negative cases; 25 519 controls of European ancestry), Triple Negative Breast Cancer Consortium (TNBCC) (1562 triple negative cases; 3399 controls of European ancestry) and African American Breast Cancer Consortium (AABC) (1004 ER-negative cases; 2745 controls). We performed in silico replication of 86 SNPs at P ≤ 1 × 10(-5) in an additional 11 209 breast cancer cases (946 with ER-negative disease) and 16 057 controls of Japanese, Latino and European ancestry. We identified two novel loci for breast cancer at 20q11 and 6q14. SNP rs2284378 at 20q11 was associated with ER-negative breast cancer (combined two-stage OR = 1.16; P = 1.1 × 10(-8)) but showed a weaker association with overall breast cancer (OR = 1.08, P = 1.3 × 10(-6)) based on 17 869 cases and 43 745 controls and no association with ER-positive disease (OR = 1.01, P = 0.67) based on 9965 cases and 22 902 controls. Similarly, rs17530068 at 6q14 was associated with breast cancer (OR = 1.12; P = 1.1 × 10(-9)), and with both ER-positive (OR = 1.09; P = 1.5 × 10(-5)) and ER-negative (OR = 1.16, P = 2.5 × 10(-7)) disease. We also confirmed three known loci associated with ER-negative (19p13) and both ER-negative and ER-positive breast cancer (6q25 and 12p11). Our results highlight the value of large-scale collaborative studies to identify novel breast cancer risk loci.

Genome-wide association study of endometrial cancer in E2C2.

Endometrial cancer (EC), a neoplasm of the uterine epithelial lining, is the most common gynecological malignancy in developed countries and the fourth most common cancer among US women. Women with a family history of EC have an increased risk for the disease, suggesting that inherited genetic factors play a role. We conducted a two-stage genome-wide association study of Type I EC. Stage 1 included 5,472 women (2,695 cases and 2,777 controls) of European ancestry from seven studies. We selected independent single-nucleotide polymorphisms (SNPs) that displayed the most significant associations with EC in Stage 1 for replication among 17,948 women (4,382 cases and 13,566 controls) in a multiethnic population (African America, Asian, Latina, Hawaiian and European ancestry), from nine studies. Although no novel variants reached genome-wide significance, we replicated previously identified associations with genetic markers near the HNF1B locus. Our findings suggest that larger studies with specific tumor classification are necessary to identify novel genetic polymorphisms associated with EC susceptibility.

A case series of thyroid fine needle aspiration biopsies diagnosed as follicular neoplasm with oncocytic features.

Oncocytic lesions of the thyroid are a heterogeneous group encompassing nonneoplastic and neoplastic entities ranging from benign to malignant and have traditionally been classified as separate entities in thyroid pathology. To illustrate the diversity of these thyroid lesions, we describe three cases of fine needle aspiration biopsies (FNAB) diagnosed as Bethesda Category IV: Follicular neoplasm, oncocytic type, under the 2017 Bethesda System for Reporting Thyroid Cytopathology (TBSRTC), with ThyroSeq v3 molecular testing and subsequent surgical excision.

Utility and performance of cell blocks in cerebrospinal fluid cytology: Experience at two teaching hospitals.

BACKGROUND: Cytology cell blocks (CBs) are not routinely made for cerebrospinal fluid (CSF) specimens. The goal of this study was to identify when CSF CB preparation improves diagnostic performance. MATERIALS AND METHODS: Under institutional review board approval, a retrospective review of CSF cytology cases was conducted at a tertiary university-based hospital and an affiliated county hospital. Patient history, CSF volume, final diagnosis, use of stains, and whether the CB was contributory was determined from the cytopathology report. CSF nucleated cell count data was obtained from the medical record. RESULTS: A total of 69 CSF specimens with CBs from January 2006 to March 2023 were identified from 61 patients. The median CSF volume was 8 mL (interquartile range, 4-13 mL; range, 1-800 mL), with immunohistochemical stains performed on 29 (42%) cases. Per cytology report, CB was contributory in 23 cases (33%), not contributory in 34 cases (49%), and not discussed in 12 cases (17%). The median volume was 8 mL for cases in which CB was contributory, not contributory, or not discussed. There was no difference in average nucleated cell counts between cases in which CB was contributory versus not contributory (73.9 vs. 40.0, p = .175). CONCLUSIONS: CBs for CSF samples were contributory in a subset (33%) of cases. The authors were unable to identify any specific pre-analytic factors, including specimen volume and average nucleated cell counts, for cases in which CB was contributory. Further evaluation is needed to identify if there are scenarios in which CSF CBs should be routinely prepared.

Genome-wide association study identifies novel alleles associated with risk of cutaneous basal cell carcinoma and squamous cell carcinoma.

We conducted a genome-wide association study on cutaneous basal cell carcinoma (BCC) among 2045 cases and 6013 controls of European ancestry, with follow-up replication in 1426 cases and 4845 controls. A non-synonymous SNP in the MC1R gene (rs1805007 encoding Arg151Cys substitution), a previously well-documented pigmentation gene, showed the strongest association with BCC risk in the discovery set (rs1805007[T]: OR (95% CI) for combined discovery set and replication set [1.55 (1.45-1.66); P= 4.3 × 10(-17)]. We identified that an SNP rs12210050 at 6p25 near the EXOC2 gene was associated with an increased risk of BCC [rs12210050[T]: combined OR (95% CI), 1.24 (1.17-1.31); P= 9.9 × 10(-10)]. In the locus on 13q32 near the UBAC2 gene encoding ubiquitin-associated domain-containing protein 2, we also identified a variant conferring susceptibility to BCC [rs7335046 [G]; combined OR (95% CI), 1.26 (1.18-1.34); P= 2.9 × 10(-8)]. We further evaluated the associations of these two novel SNPs (rs12210050 and rs7335046) with squamous cell carcinoma (SCC) risk as well as melanoma risk. We found that both variants, rs12210050[T] [OR (95% CI), 1.35 (1.16-1.57); P= 7.6 × 10(-5)] and rs7335046 [G] [OR (95% CI), 1.21 (1.02-1.44); P= 0.03], were associated with an increased risk of SCC. These two variants were not associated with melanoma risk. We conclude that 6p25 and 13q32 are novel loci conferring susceptibility to non-melanoma skin cancer.