Evidence of Gene–Environment Interactions between Common Breast Cancer Susceptibility Loci and Established Environmental Risk Factors

Various common genetic susceptibility loci have been identified for breast cancer; however, it is unclear how they combine with lifestyle/environmental risk factors to influence risk. We undertook an international collaborative study to assess gene-environment interaction for risk of breast cancer. Data from 24 studies of the Breast Cancer Association Consortium were pooled. Using up to 34,793 invasive breast cancers and 41,099 controls, we examined whether the relative risks associated with 23 single nucleotide polymorphisms were modified by 10 established environmental risk factors (age at menarche, parity, breastfeeding, body mass index, height, oral contraceptive use, menopausal hormone therapy use, alcohol consumption, cigarette smoking, physical activity) in women of European ancestry. We used logistic regression models stratified by study and adjusted for age and performed likelihood ratio tests to assess gene–environment interactions. All statistical tests were two-sided. We replicated previously reported potential interactions between LSP1-rs3817198 and parity (P_interaction = 2.4×10⁻⁶) and between CASP8-rs17468277 and alcohol consumption (P_interaction = 3.1×10⁻⁴). Overall, the per-allele odds ratio (95% confidence interval) for LSP1-rs3817198 was 1.08 (1.01–1.16) in nulliparous women and ranged from 1.03 (0.96–1.10) in parous women with one birth to 1.26 (1.16–1.37) in women with at least four births. For CASP8-rs17468277, the per-allele OR was 0.91 (0.85–0.98) in those with an alcohol intake of <20 g/day and 1.45 (1.14–1.85) in those who drank ≥20 g/day. Additionally, interaction was found between 1p11.2-rs11249433 and ever being parous (P_interaction = 5.3×10⁻⁵), with a per-allele OR of 1.14 (1.11–1.17) in parous women and 0.98 (0.92–1.05) in nulliparous women. These data provide first strong evidence that the risk of breast cancer associated with some common genetic variants may vary with environmental risk factors.     

Assessing Interactions between Common Genetic Variant on 2q35 and Hormone Receptor Status with Breast Cancer Risk: Evidence Based on 26 Studies

Genome-wide association studies have identified 2q35-rs13387042 as a new breast cancer (BC) susceptibility locus in populations of European descent. Since then, the relationship between 2q35-rs13387042 and breast cancer has been reported in various ethnic groups; however, these studies have yielded inconsistent results. To investigate this inconsistency, we performed a meta-analysis of 26 studies involving a total of 101,529 cases and 167,363 controls for 2q35-rs13387042 polymorphism to evaluate its effect on genetic susceptibility for breast cancer. An overall random effects odds ratio of 1.14 (95% CI: 1.11–1.16, P<10⁻⁵) was found for rs13387042-A variant. Significant results were also observed using dominant (OR = 1.14, 95% CI: 1.12–1.17, P<10⁻⁵), recessive (OR = 1.17, 95% CI: 1.13–1.21, P<10⁻⁵) and co-dominant genetic model (heterozygous: OR = 1.15, 95% CI: 1.12–1.19, P<10⁻⁵; homozygous: OR = 1.20, 95% CI: 1.15–1.24, P<10⁻⁵). There was strong evidence of heterogeneity, which largely disappeared after stratification by ethnicity. Significant associations were found in East Asians, and White populations when stratified by ethnicity; while no significant associations were observed in Africans and other ethnic populations. An association was observed for both ER-positive (OR = 1.17, 95% 1.15–1.19; P<10⁻⁵) and ER-negative disease (OR = 1.08, 95% CI: 1.04–1.13; P<10⁻⁴) and both progesterone receptor (PR)-positive (OR = 1.18, 95% CI: 1.15–1.21; P<10⁻⁵) and PR-negative disease (OR = 1.10, 95% CI: 1.05–1.15; P<10⁻⁴). In conclusion, this meta-analysis demonstrated that the A allele of 2q35-rs13387042 is a risk factor associated with increased breast cancer susceptibility.     

Identification of a Functional Genetic Variant at 16q12.1 for Breast Cancer Risk: Results from the Asia Breast Cancer Consortium

Genetic factors play an important role in the etiology of breast cancer. We carried out a multi-stage genome-wide association (GWA) study in over 28,000 cases and controls recruited from 12 studies conducted in Asian and European American women to identify genetic susceptibility loci for breast cancer. After analyzing 684,457 SNPs in 2,073 cases and 2,084 controls in Chinese women, we evaluated 53 SNPs for fast-track replication in an independent set of 4,425 cases and 1,915 controls of Chinese origin. Four replicated SNPs were further investigated in an independent set of 6,173 cases and 6,340 controls from seven other studies conducted in Asian women. SNP rs4784227 was consistently associated with breast cancer risk across all studies with adjusted odds ratios (95% confidence intervals) of 1.25 (1.20−1.31) per allele (P = 3.2×10⁻²⁵) in the pooled analysis of samples from all Asian samples. This SNP was also associated with breast cancer risk among European Americans (per allele OR  = 1.19, 95% CI  = 1.09−1.31, P = 1.3×10⁻⁴, 2,797 cases and 2,662 controls). SNP rs4784227 is located at 16q12.1, a region identified previously for breast cancer risk among Europeans. The association of this SNP with breast cancer risk remained highly statistically significant in Asians after adjusting for previously-reported SNPs in this region. In vitro experiments using both luciferase reporter and electrophoretic mobility shift assays demonstrated functional significance of this SNP. These results provide strong evidence implicating rs4784227 as a functional causal variant for breast cancer in the locus 16q12.1 and demonstrate the utility of conducting genetic association studies in populations with different genetic architectures.     

Ancestry-Shift Refinement Mapping of the C6orf97-ESR1 Breast Cancer Susceptibility Locus

We used an approach that we term ancestry-shift refinement mapping to investigate an association, originally discovered in a GWAS of a Chinese population, between rs2046210[T] and breast cancer susceptibility. The locus is on 6q25.1 in proximity to the C6orf97 and estrogen receptor α (ESR1) genes. We identified a panel of SNPs that are correlated with rs2046210 in Chinese, but not necessarily so in other ancestral populations, and genotyped them in breast cancer case∶control samples of Asian, European, and African origin, a total of 10,176 cases and 13,286 controls. We found that rs2046210[T] does not confer substantial risk of breast cancer in Europeans and Africans (OR = 1.04, P = 0.099, and OR = 0.98, P = 0.77, respectively). Rather, in those ancestries, an association signal arises from a group of less common SNPs typified by rs9397435. The rs9397435[G] allele was found to confer risk of breast cancer in European (OR = 1.15, P = 1.2×10⁻³), African (OR = 1.35, P = 0.014), and Asian (OR = 1.23, P = 2.9×10⁻⁴) population samples. Combined over all ancestries, the OR was 1.19 (P = 3.9×10⁻⁷), was without significant heterogeneity between ancestries (P_het = 0.36) and the SNP fully accounted for the association signal in each ancestry. Haplotypes bearing rs9397435[G] are well tagged by rs2046210[T] only in Asians. The rs9397435[G] allele showed associations with both estrogen receptor positive and estrogen receptor negative breast cancer. Using early-draft data from the 1,000 Genomes project, we found that the risk allele of a novel SNP (rs77275268), which is closely correlated with rs9397435, disrupts a partially methylated CpG sequence within a known CTCF binding site. These studies demonstrate that shifting the analysis among ancestral populations can provide valuable resolution in association mapping.     

A Genome Wide Meta-Analysis Study for Identification of Common Variation Associated with Breast Cancer Prognosis

Objective

Genome wide association studies (GWAs) of breast cancer mortality have identified few potential associations. The concordance between these studies is unclear. In this study, we used a meta-analysis of two prognostic GWAs and a replication cohort to identify the strongest associations and to evaluate the loci suggested in previous studies. We attempt to identify those SNPs which could impact overall survival irrespective of the age of onset.

Methods

To facilitate the meta-analysis and to refine the association signals, SNPs were imputed using data from the 1000 genomes project. Cox-proportional hazard models were used to estimate hazard ratios (HR) in 536 patients from the POSH cohort (Prospective study of Outcomes in Sporadic versus Hereditary breast cancer) and 805 patients from the HEBCS cohort (Helsinki Breast Cancer Study). These hazard ratios were combined using a Mantel-Haenszel fixed effects meta-analysis and a p-value threshold of 5×10⁻⁸ was used to determine significance. Replication was performed in 1523 additional patients from the POSH study.

Results

Although no SNPs achieved genome wide significance, three SNPs have significant association in the replication cohort and combined p-values less than 5.6×10⁻⁶. These SNPs are; rs421379 which is 556 kb upstream of ARRDC3 (HR = 1.49, 95% confidence interval (CI) = 1.27–1.75, P = 1.1×10⁻⁶), rs12358475 which is between ECHDC3 and PROSER2 (HR = 0.75, CI = 0.67–0.85, P = 1.8×10⁻⁶), and rs1728400 which is between LINC00917 and FOXF1.

Conclusions

In a genome wide meta-analysis of two independent cohorts from UK and Finland, we identified potential associations at three distinct loci. Phenotypic heterogeneity and relatively small sample sizes may explain the lack of genome wide significant findings. However, the replication at three SNPs in the validation cohort shows promise for future studies in larger cohorts. We did not find strong evidence for concordance between the few associations highlighted by previous GWAs of breast cancer survival and this study.     

A Common Variant in the SIAH2 Locus Is Associated with Estrogen Receptor-Positive Breast Cancer in the Chinese Han Population

Background

Genome-wide association studies (GWAS) have identified many loci associated with breast cancer risk. These studies have primarily been conducted in populations of European descent.

Objective

To determine whether previously reported susceptibility loci in other ethnic groups are also risk factors for breast cancer in a Chinese population.

Method

We genotyped 21 previously reported single nucleotide polymorphisms (SNPs) within a female Chinese cohort of 1203 breast cancer cases and 2525 healthy controls using the Sequenom iPlex platform. Fourteen SNPs passed the quality control test. These SNPs were subjected to statistical analysis for the entire cohort and were further analyzed for estrogen receptor (ER) status. The associations of the SNPs with disease susceptibility were assessed using logistic regression, adjusting for age. The Bonferroni correction was used to conservatively account for multiple testing, and the threshold for statistical significance was P<3.57×10⁻³ (0.05/14).

Result

Although none of the SNPs showed an overall association with breast cancer, an analysis of the ER status of the breast cancer patients revealed that the SIAH2 locus (rs6788895; P = 5.73×10⁻⁴, odds ratio [OR] = 0.81) is associated with ER-positive breast cancer.

Conclusion

A common variant in the SIAH2 locus is associated with ER-positive breast cancer in the Chinese Han population. The replication of published GWAS results in other ethnic groups provides important information regarding the genetic etiology of breast cancer.     

The role of osteopontin (OPN/SPP1) haplotypes in the susceptibility to Crohn's disease

Background

Osteopontin represents a multifunctional molecule playing a pivotal role in chronic inflammatory and autoimmune diseases. Its expression is increased in inflammatory bowel disease (IBD). The aim of our study was to analyze the association of osteopontin (OPN/SPP1) gene variants in a large cohort of IBD patients.

Methodology/Principal Findings

Genomic DNA from 2819 Caucasian individuals (n = 841 patients with Crohn''s disease (CD), n = 473 patients with ulcerative colitis (UC), and n = 1505 healthy unrelated controls) was analyzed for nine OPN SNPs (rs2728127, rs2853744, rs11730582, rs11739060, rs28357094, rs4754 = p.Asp80Asp, rs1126616 = p.Ala236Ala, rs1126772 and rs9138). Considering the important role of osteopontin in Th17-mediated diseases, we performed analysis for epistasis with IBD-associated IL23R variants and analyzed serum levels of the Th17 cytokine IL-22. For four OPN SNPs (rs4754, rs1126616, rs1126772 and rs9138), we observed significantly different distributions between male and female CD patients. rs4754 was protective in male CD patients (p = 0.0004, OR = 0.69). None of the other investigated OPN SNPs was associated with CD or UC susceptibility. However, several OPN haplotypes showed significant associations with CD susceptibility. The strongest association was found for a haplotype consisting of the 8 OPN SNPs rs2728127-rs2853744-rs11730582-rs11439060-rs28357094-rs112661-rs1126772-rs9138 (omnibus p-value = 2.07×10⁻⁸). Overall, the mean IL-22 secretion in the combined group of OPN minor allele carriers with CD was significantly lower than that of CD patients with OPN wildtype alleles (p = 3.66×10⁻⁵). There was evidence for weak epistasis between the OPN SNP rs28357094 with the IL23R SNP rs10489629 (p = 4.18×10⁻²) and between OPN SNP rs1126616 and IL23R SNP rs2201841 (p = 4.18×10⁻²) but none of these associations remained significant after Bonferroni correction.

Conclusions/Significance

Our study identified OPN haplotypes as modifiers of CD susceptibility, while the combined effects of certain OPN variants may modulate IL-22 secretion.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.     

Genome-Wide Testing of Putative Functional Exonic Variants in Relationship with Breast and Prostate Cancer Risk in a Multiethnic Population

Rare variation in protein coding sequence is poorly captured by GWAS arrays and has been hypothesized to contribute to disease heritability. Using the Illumina HumanExome SNP array, we successfully genotyped 191,032 common and rare non-synonymous, splice site, or nonsense variants in a multiethnic sample of 2,984 breast cancer cases, 4,376 prostate cancer cases, and 7,545 controls. In breast cancer, the strongest associations included either SNPs in or gene burden scores for genes LDLRAD1, SLC19A1, FGFBP3, CASP5, MMAB, SLC16A6, and INS-IGF2. In prostate cancer, one of the most associated SNPs was in the gene GPRC6A (rs2274911, Pro91Ser, OR = 0.88, P = 1.3×10⁻⁵) near to a known risk locus for prostate cancer; other suggestive associations were noted in genes such as F13A1, ANXA4, MANSC1, and GP6. For both breast and prostate cancer, several of the most significant associations involving SNPs or gene burden scores (sum of minor alleles) were noted in genes previously reported to be associated with a cancer-related phenotype. However, only one of the associations (rs145889899 in LDLRAD1, p = 2.5×10⁻⁷ only seen in African Americans) for overall breast or prostate cancer risk was statistically significant after correcting for multiple comparisons. In addition to breast and prostate cancer, other cancer-related traits were examined (body mass index, PSA level, and alcohol drinking) with a number of known and potentially novel associations described. In general, these findings do not support there being many protein coding variants of moderate to high risk for breast and prostate cancer with odds ratios over a range that is probably required for protein coding variation to play a truly outstanding role in risk heritability. Very large sample sizes will be required to better define the role of rare and less penetrant coding variation in prostate and breast cancer disease genetics.     

Identification of Novel Genetic Loci Associated with Thyroid Peroxidase Antibodies and Clinical Thyroid Disease

Autoimmune thyroid diseases (AITD) are common, affecting 2-5% of the general population. Individuals with positive thyroid peroxidase antibodies (TPOAbs) have an increased risk of autoimmune hypothyroidism (Hashimoto''s thyroiditis), as well as autoimmune hyperthyroidism (Graves'' disease). As the possible causative genes of TPOAbs and AITD remain largely unknown, we performed GWAS meta-analyses in 18,297 individuals for TPOAb-positivity (1769 TPOAb-positives and 16,528 TPOAb-negatives) and in 12,353 individuals for TPOAb serum levels, with replication in 8,990 individuals. Significant associations (P<5×10⁻⁸) were detected at TPO-rs11675434, ATXN2-rs653178, and BACH2-rs10944479 for TPOAb-positivity, and at TPO-rs11675434, MAGI3-rs1230666, and KALRN-rs2010099 for TPOAb levels. Individual and combined effects (genetic risk scores) of these variants on (subclinical) hypo- and hyperthyroidism, goiter and thyroid cancer were studied. Individuals with a high genetic risk score had, besides an increased risk of TPOAb-positivity (OR: 2.18, 95% CI 1.68–2.81, P = 8.1×10⁻⁸), a higher risk of increased thyroid-stimulating hormone levels (OR: 1.51, 95% CI 1.26–1.82, P = 2.9×10⁻⁶), as well as a decreased risk of goiter (OR: 0.77, 95% CI 0.66–0.89, P = 6.5×10⁻⁴). The MAGI3 and BACH2 variants were associated with an increased risk of hyperthyroidism, which was replicated in an independent cohort of patients with Graves'' disease (OR: 1.37, 95% CI 1.22–1.54, P = 1.2×10⁻⁷ and OR: 1.25, 95% CI 1.12–1.39, P = 6.2×10⁻⁵). The MAGI3 variant was also associated with an increased risk of hypothyroidism (OR: 1.57, 95% CI 1.18–2.10, P = 1.9×10⁻³). This first GWAS meta-analysis for TPOAbs identified five newly associated loci, three of which were also associated with clinical thyroid disease. With these markers we identified a large subgroup in the general population with a substantially increased risk of TPOAbs. The results provide insight into why individuals with thyroid autoimmunity do or do not eventually develop thyroid disease, and these markers may therefore predict which TPOAb-positives are particularly at risk of developing clinical thyroid dysfunction.     

How Genome-Wide SNP-SNP Interactions Relate to Nasopharyngeal Carcinoma Susceptibility

This study is the first to use genome-wide association study (GWAS) data to evaluate the multidimensional genetic architecture underlying nasopharyngeal cancer. Since analysis of data from GWAS confirms a close and consistent association between elevated risk for nasopharyngeal carcinoma (NPC) and major histocompatibility complex class 1 genes, our goal here was to explore lesser effects of gene-gene interactions. We conducted an exhaustive genome-wide analysis of GWAS data of NPC, revealing two-locus interactions occurring between single nucleotide polymorphisms (SNPs), and identified a number of suggestive interaction loci which were missed by traditional GWAS analyses. Although none of the interaction pairs we identified passed the genome-wide Bonferroni-adjusted threshold for significance, using independent GWAS data from the same population (Stage 2), we selected 66 SNP pairs in 39 clusters with P<0.01. We identified that in several chromosome regions, multiple suggestive interactions group to form a block-like signal, effectively reducing the rate of false discovery. The strongest cluster of interactions involved the CREB5 gene and a SNP rs1607979 on chromosome 17q22 (P = 9.86×10⁻¹¹) which also show trans-expression quantitative loci (eQTL) association in Chinese population. We then detected a complicated cis-interaction pattern around the NPC-associated HLA-B locus, which is immediately adjacent to copy-number variations implicated in male susceptibility for NPC. While it remains to be seen exactly how and to what degree SNP-SNP interactions such as these affect susceptibility for nasopharyngeal cancer, future research on these questions holds great promise for increasing our understanding of this disease’s genetic etiology, and possibly also that of other gene-related cancers.     

A Genetic Variant in 12q13, a Possible Risk Factor for Bipolar Disorder,Is Associated with Depressive State,Accounting for Stressful Life Events

Genome-wide association studies (GWASs) have identified a number of susceptibility genes for schizophrenia (SCZ) and bipolar disorder (BD). However, the identification of risk genes for major depressive disorder (MDD) has been unsuccessful because the etiology of MDD is more influenced by environmental factors; thus, gene–environment (G×E) interactions are important, such as interplay with stressful life events (SLEs). We assessed the G×E interactions and main effects of genes targeting depressive symptoms. Using a case–control design, 922 hospital staff members were evaluated for depressive symptoms according to Beck Depressive Inventory (BDI; “depression” and “control” groups were classified by scores of 10 in the BDI test), SLEs, and personality. A total of sixty-three genetic variants were selected on the basis of previous GWASs of MDD, SCZ, and BD as well as candidate-gene (SLC6A4, BDNF, DBH, and FKBP5) studies. Logistic regression analysis revealed a marginally significant interaction (genetic variant × SLE) at rs4523957 (P_uncorrected = 0.0034) with depression and a significant association of single nucleotide polymorphism identified from evidence of BD GWAS (rs7296288, downstream of DHH at 12q13.1) with depression as the main effect (P_uncorrected = 9.4×10⁻⁴, P_corrected = 0.0424). We also found that SLEs had a larger impact on depression (odds ratio∼3), as reported previously. These results suggest that DHH plays a possible role in depression etiology; however, variants from MDD or SCZ GWAS evidence or candidate genes showed no significant associations or minimal effects of interactions with SLEs on depression.     

Host Genetic Factors and Vaccine-Induced Immunity to HBV Infection: Haplotype Analysis

Hepatitis B virus (HBV) infection remains a significant health burden world-wide, although vaccines help decrease this problem. We previously identified associations of single nucleotide polymorphisms in several candidate genes with vaccine-induced peak antibody level (anti-HBs), which is predictive of long-term vaccine efficacy and protection against infection and persistent carriage; here we report on a haplotype-based analysis. A total of 688 SNPs from 117 genes were examined for a two, three and four sliding window haplotype analysis in a Gambian cohort. Analysis was performed on 197 unrelated individuals, 454 individuals from 174 families, and the combined sample (N = 651). Global and individual haplotype association tests were carried out (adjusted for covariates), employing peak anti-HBs level as outcome. Five genes (CD44, CD58, CDC42, IL19 and IL1R1) had at least one significant haplotype in the unrelated or family analysis as well as the combined analysis. Previous single locus results were confirmed for CD44 (combined global p = 9.1×10⁻⁵ for rs353644-rs353630-rs7937602) and CD58 (combined global p = 0.008 for rs1414275-rs11588376-rs1016140). Haplotypes in CDC42, IL19 and IL1R1 also associated with peak anti-HBs level. We have identified strong haplotype effects on HBV vaccine-induced antibody level in five genes, three of which, CDC42, IL19 and IL1R1, did not show evidence of association in a single SNP analyses and corroborated the majority of these effects in two datasets. The haplotype analysis identified associations with HBV vaccine-induced immunity in several new genes.     

Genetic Susceptibility to Refractive Error: Association of Vasoactive Intestinal Peptide Receptor 2 (VIPR2) with High Myopia in Chinese

Myopia is the most common ocular disease worldwide. We investigated the association of high myopia with the common single nucleotide polymorphisms (SNPs) of five candidate genes – early growth response 1 (EGR1), v-fos FBJ murine osteosarcoma viral oncogene homolog (FOS), jun oncogene (JUN), vasoactive intestinal peptide (VIP), and vasoactive intestinal peptide receptor 2 (VIPR2). We recruited 1200 unrelated Chinese subjects with 600 cases (spherical equivalent ≤−8.00 diopters) and 600 controls (spherical equivalent within ±1.00 diopter). A discovery sample set was formed from 300 cases and 300 controls, and a replication sample set from the remaining samples. Tag SNPs were genotyped for the discovery sample set, and the most significant haplotypes and their constituent SNPs were followed up with the replication sample set. The allele and haplotype frequencies in cases and controls were compared by logistic regression adjusted for sex and age to give P _a values, and multiple comparisons were corrected by permutation test to give P _aemp values. Odd ratios (OR) were calculated accordingly. In the discovery phase, EGR1, JUN and VIP did not show any significant association while FOS and VIPR2 demonstrated significant haplotype association with high myopia. In the replication phase, the haplotype association for VIPR2 was successfully replicated, but not FOS. In analysis combining both sample sets, the most significant association signals of VIPR2 were the single marker rs2071625 (P _a = 0.0008, P _aemp = 0.0046 and OR = 0.75) and the 4-SNP haplotype window rs2071623-rs2071625-rs2730220-rs885863 (omnibus test, P _a = 9.10e-10 and P _aemp = 0.0001) with one protective haplotype (GGGG: P _aemp = 0.0002 and OR = 0.52) and one high-risk haplotype (GAGA: P _aemp = 0.0027 and OR = 4.68). This 4-SNP haplotype window was the most significant in all sample sets examined. This is the first study to suggest a role of VIPR2 in the genetic susceptibility to high myopia. EGR1, JUN, FOS and VIP are unlikely to be important in predisposing humans to high myopia.     

A Genome-Wide Association Study Reveals Variants in ARL15 that Influence Adiponectin Levels

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P≤5×10⁻⁸). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P≤0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2×10⁻¹⁹ for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9×10⁻⁸, n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5×10⁻⁶, n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2×10⁻³, n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

Genetic Predisposition to In Situ and Invasive Lobular Carcinoma of the Breast

Invasive lobular breast cancer (ILC) accounts for 10–15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09–1.18), P = 6.0×10⁻¹⁰; P-het for ILC vs IDC ER+ tumors = 1.8×10⁻⁴). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P<0.05. Two SNPs showed significantly stronger associations for ILC than LCIS (rs2981579/10q26/FGFR2, P-het = 0.04 and rs889312/5q11/MAP3K1, P-het = 0.03); and two showed stronger associations for LCIS than ILC (rs6678914/1q32/LGR6, P-het = 0.001 and rs1752911/6q14, P-het = 0.04). In addition, seven of the 75 known loci showed significant differences between ER+ tumors with IDC and ILC histology, three of these showing stronger associations for ILC (rs11249433/1p11, rs2981579/10q26/FGFR2 and rs10995190/10q21/ZNF365) and four associated only with IDC (5p12/rs10941679; rs2588809/14q24/RAD51L1, rs6472903/8q21 and rs1550623/2q31/CDCA7). In conclusion, we have identified one novel lobular breast cancer specific predisposition polymorphism at 7q34, and shown for the first time that common breast cancer polymorphisms predispose to LCIS. We have shown that many of the ER+ breast cancer predisposition loci also predispose to ILC, although there is some heterogeneity between ER+ lobular and ER+ IDC tumors. These data provide evidence for overlapping, but distinct etiological pathways within ER+ breast cancer between morphological subtypes.     

Associations of Polymorphisms in Histidine Decarboxylase,Histamine N-Methyltransferase and Histamine Receptor H3 Genes with Breast Cancer

We previously found that genetic polymorphisms in gene coding for histamine H₄ receptors were related to the risk and malignant degree of breast cancer. The roles of polymorphisms in other histamine-related genes, such as histidine decarboxylase (HDC), histamine N-methyltransferase (HNMT) and histamine H₃ receptor (HRH3), remain unexplored. The aim of this study is to analyze the clinical associations of polymorphisms in HDC, HNMT and HRH3 with breast cancer. Two hundred and one unrelated Chinese Han breast cancer patients and 205 ethnicity-matched health controls were recruited for case-control investigation. Genomic DNA from the participants was extracted and 5 single nucleotide polymorphisms (SNPs) in HDC, HNMT and HRH3 were genotyped. We found that polymorphisms of HNMT and HRH3 were irrelevant with breast cancer in the present study. However, the T allele of rs7164386 in HDC significantly decreased the risk of breast cancer (adjusted odds ratios [ORs], 0.387; 95% confidence intervals [CIs], 0.208–0.720; P = 0.003). Furthermore, for HDC haplotypes, the CG haplotype of rs7164386-rs7182203 was more frequent among breast cancer patients (adjusted OR, 1.828; 95% CI, 1.218–2.744; P = 0.004) while the TG haplotype was more frequent among health controls (adjusted OR, 0.351; 95% CI, 0.182–0.678; P = 0.002). These findings indicated that polymorphisms of HDC gene were significantly associated with breast cancer in Chinese Han population and may be novel diagnostic or therapeutic targets for breast cancer. Further studies with larger participants worldwide are still needed for conclusion validation.     

Multi-Variant Pathway Association Analysis Reveals the Importance of Genetic Determinants of Estrogen Metabolism in Breast and Endometrial Cancer Susceptibility

Despite the central role of estrogen exposure in breast and endometrial cancer development and numerous studies of genes in the estrogen metabolic pathway, polymorphisms within the pathway have not been consistently associated with these cancers. We posit that this is due to the complexity of multiple weak genetic effects within the metabolic pathway that can only be effectively detected through multi-variant analysis. We conducted a comprehensive association analysis of the estrogen metabolic pathway by interrogating 239 tagSNPs within 35 genes of the pathway in three tumor samples. The discovery sample consisted of 1,596 breast cancer cases, 719 endometrial cancer cases, and 1,730 controls from Sweden; and the validation sample included 2,245 breast cancer cases and 1,287 controls from Finland. We performed admixture maximum likelihood (AML)–based global tests to evaluate the cumulative effect from multiple SNPs within the whole metabolic pathway and three sub-pathways for androgen synthesis, androgen-to-estrogen conversion, and estrogen removal. In the discovery sample, although no single polymorphism was significant after correction for multiple testing, the pathway-based AML global test suggested association with both breast (p _global = 0.034) and endometrial (p _global = 0.052) cancers. Further testing revealed the association to be focused on polymorphisms within the androgen-to-estrogen conversion sub-pathway, for both breast (p _global = 0.008) and endometrial cancer (p _global = 0.014). The sub-pathway association was validated in the Finnish sample of breast cancer (p _global = 0.015). Further tumor subtype analysis demonstrated that the association of the androgen-to-estrogen conversion sub-pathway was confined to postmenopausal women with sporadic estrogen receptor positive tumors (p _global = 0.0003). Gene-based AML analysis suggested CYP19A1 and UGT2B4 to be the major players within the sub-pathway. Our study indicates that the composite genetic determinants related to the androgen–estrogen conversion are important for the induction of two hormone-associated cancers, particularly for the hormone-driven breast tumour subtypes.