Breast Cancer Risk Genes - Association Analysis in More than 113,000 Women

BACKGROUND

Genetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking.

METHODS

We used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity.

RESULTS

Protein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants.

CONCLUSIONS

The results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.).

A variant of the HTRA1 gene increases susceptibility to age-related macular degeneration

Age-related macular degeneration (AMD) is the most common cause of irreversible vision loss in the developed world and has a strong genetic predisposition. A locus at human chromosome 10q26 affects the risk of AMD, but the precise gene(s) have not been identified. We genotyped 581 AMD cases and 309 normal controls in a Caucasian cohort in Utah. We demonstrate that a single-nucleotide polymorphism, rs11200638, in the promoter region of HTRA1 is the most likely causal variant for AMD at 10q26 and is estimated to confer a population attributable risk of 49.3%. The HTRA1 gene encodes a secreted serine protease. Preliminary analysis of lymphocytes and retinal pigment epithelium from four AMD patients revealed that the risk allele was associated with elevated expression levels of HTRA1 mRNA and protein. We also found that drusen in the eyes of AMD patients were strongly immunolabeled with HTRA1 antibody. Together, these findings support a key role for HTRA1 in AMD susceptibility and identify a potential new pathway for AMD pathogenesis.

A candidate prostate cancer susceptibility gene at chromosome 17p

It is difficult to identify genes that predispose to prostate cancer due to late age at diagnosis, presence of phenocopies within high-risk pedigrees and genetic complexity. A genome-wide scan of large, high-risk pedigrees from Utah has provided evidence for linkage to a locus on chromosome 17p. We carried out positional cloning and mutation screening within the refined interval, identifying a gene, ELAC2, harboring mutations (including a frameshift and a nonconservative missense change) that segregate with prostate cancer in two pedigrees. In addition, two common missense variants in the gene are associated with the occurrence of prostate cancer. ELAC2 is a member of an uncharacterized gene family predicted to encode a metal-dependent hydrolase domain that is conserved among eukaryotes, archaebacteria and eubacteria. The gene product bears amino acid sequence similarity to two better understood protein families, namely the PSO2 (SNM1) DNA interstrand crosslink repair proteins and the 73-kD subunit of mRNA 3' end cleavage and polyadenylation specificity factor (CPSF73).

HOXB13 is a susceptibility gene for prostate cancer: results from the International Consortium for Prostate Cancer Genetics (ICPCG)

Prostate cancer has a strong familial component but uncovering the molecular basis for inherited susceptibility for this disease has been challenging. Recently, a rare, recurrent mutation (G84E) in HOXB13 was reported to be associated with prostate cancer risk. Confirmation and characterization of this finding is necessary to potentially translate this information to the clinic. To examine this finding in a large international sample of prostate cancer families, we genotyped this mutation and 14 other SNPs in or flanking HOXB13 in 2,443 prostate cancer families recruited by the International Consortium for Prostate Cancer Genetics (ICPCG). At least one mutation carrier was found in 112 prostate cancer families (4.6 %), all of European descent. Within carrier families, the G84E mutation was more common in men with a diagnosis of prostate cancer (194 of 382, 51 %) than those without (42 of 137, 30 %), P = 9.9 × 10(-8) [odds ratio 4.42 (95 % confidence interval 2.56-7.64)]. A family-based association test found G84E to be significantly over-transmitted from parents to affected offspring (P = 6.5 × 10(-6)). Analysis of markers flanking the G84E mutation indicates that it resides in the same haplotype in 95 % of carriers, consistent with a founder effect. Clinical characteristics of cancers in mutation carriers included features of high-risk disease. These findings demonstrate that the HOXB13 G84E mutation is present in ~5 % of prostate cancer families, predominantly of European descent, and confirm its association with prostate cancer risk. While future studies are needed to more fully define the clinical utility of this observation, this allele and others like it could form the basis for early, targeted screening of men at elevated risk for this common, clinically heterogeneous cancer.

Interleukin-1 receptor antagonist gene polymorphism and coronary artery disease

BACKGROUND

Cytokine gene variations are contributory factors in inflammatory pathology. Allele frequencies of interleukin (IL)-1 cluster genes [IL-1A(-889), IL-1B(-511), IL-1B(+3953), IL-1RN Intron 2 VNTR] and tissue necrosis factor (TNF)-alpha gene [TNFA(-308)] were measured in healthy blood donors (healthy control subjects), patients with angiographically normal coronary arteries (patient control subjects), single-vessel coronary disease (SVD), and those with multivessel coronary disease (MVD).

METHODS AND RESULTS

Five hundred fifty-six patients attending for coronary angiography in Sheffield were studied: 130 patient control subjects, 98 SVD, and 328 MVD. Significant associations were tested in an independent population (London) of 350: 57 SVD, 191 MVD, and 102 control subjects. IL-1RN*2 frequency in Sheffield patient control subjects was the same as in 827 healthy control subjects. IL-1RN*2 was significantly overrepresented in Sheffield SVD patients (34% vs 23% in patient control subjects); IL-1RN*2 homozygotes in the SVD population (chi2 carriage=8.490, 1 df, P=0.0036). This effect was present though not quite significant in the London population (P=0. 0603). A summary trend test of the IL-1RN SVD genotype data for Sheffield and London showed a significant association with *2 (P=0. 0024). No significant effect of genotype at IL-1RN was observed in the Sheffield or London MVD populations. Genotype distribution analysis comparing the SVD and MVD populations at IL-1RN showed a highly significant trend (P=0.0007) with the use of pooled data. No significant associations were seen for the other polymorphisms.

CONCLUSIONS

IL-1RN*2 was significantly associated with SVD. A difference in genetic association between SVD and MVD was also apparent.

Predisposition locus for major depression at chromosome 12q22-12q23.2

Major depression disorder is a common psychiatric disease with a major economic impact on society. In many cases, no effective treatment is available. The etiology of major depression is complex, but it is clear that the disease is, to a large extent, determined genetically, especially among individuals with a familial history of major depression, presumably through the involvement of multiple predisposition genes in addition to an environmental component. As a first step toward identification of chromosomal loci contributing to genetic predisposition to major depression, we have conducted a genomewide scan by using 628 microsatellite markers on 1,890 individuals from 110 Utah pedigrees with a strong family history of major depression. We identified significant linkage to major depression in males at marker D12S1300 (multipoint heterogeneity LOD score 4.6; P=.00003 after adjustment for multiple testing). With additional markers, the linkage evidence became highly significant, with the multipoint heterogeneity LOD score at marker D12S1706 increasing to 6.1 (P=.0000007 after adjustment for multiple testing). This study confirms the presence of one or more genes involved in psychiatric diseases on the q arm of chromosome 12 and provides strong evidence for the existence of a sex-specific predisposition gene to major depression at 12q22-q23.2.

Genome-wide association study of follicular lymphoma identifies a risk locus at 6p21.32

To identify susceptibility loci for non-Hodgkin lymphoma subtypes, we conducted a three-stage genome-wide association study. We identified two variants associated with follicular lymphoma at 6p21.32 (rs10484561, combined P = 1.12 x 10(-29) and rs7755224, combined P = 2.00 x 10(-19); r(2) = 1.0), supporting the idea that major histocompatibility complex genetic variation influences follicular lymphoma susceptibility. We also found confirmatory evidence of a previously reported association between chronic lymphocytic leukemia/small lymphocytic lymphoma and rs735665 (combined P = 4.24 x 10(-9)).

An analysis of linkage disequilibrium in the interleukin-1 gene cluster, using a novel grouping method for multiallelic markers

In population- and family-based association studies, it is useful to have some knowledge of the patterns of linkage disequilibrium that exist between markers in candidate regions. When such studies are carried out with multiallelic markers, it is often convenient to group the alleles into a biallelic system, for analysis. In this study, we specifically examined the interleukin-1 (IL-1) gene cluster on chromosome 2, a region containing candidates for many inflammatory and autoimmune disorders. Data were collected on eight markers, four of which were multiallelic. Using these data, we investigated the effect of three allele-grouping strategies, including a novel method, on the detection of linkage disequilibrium. The novel approach, termed the "delta method," measures the deviation from the expected haplotype frequencies under linkage equilibrium, for each allelic combination. This information is then used to group the alleles, in an attempt to avoid the grouping together of alleles at one locus that are in opposite disequilibrium with the same allele at the second locus. The estimate haplotype frequencies (EH) program was used to estimate haplotype frequencies and the disequilibrium measure. In our data it was found that the delta method compared well with the other two strategies. Using this method, we found that there was a reasonable correlation between disequilibrium and physical distance in the region (r=-.540, P=.001, one-tailed). We also identified a common, eight-locus haplotype of the IL-1 gene cluster.

A combined genomewide linkage scan of 1,233 families for prostate cancer-susceptibility genes conducted by the international consortium for prostate cancer genetics

Evidence of the existence of major prostate cancer (PC)-susceptibility genes has been provided by multiple segregation analyses. Although genomewide screens have been performed in over a dozen independent studies, few chromosomal regions have been consistently identified as regions of interest. One of the major difficulties is genetic heterogeneity, possibly due to multiple, incompletely penetrant PC-susceptibility genes. In this study, we explored two approaches to overcome this difficulty, in an analysis of a large number of families with PC in the International Consortium for Prostate Cancer Genetics (ICPCG). One approach was to combine linkage data from a total of 1,233 families to increase the statistical power for detecting linkage. Using parametric (dominant and recessive) and nonparametric analyses, we identified five regions with "suggestive" linkage (LOD score >1.86): 5q12, 8p21, 15q11, 17q21, and 22q12. The second approach was to focus on subsets of families that are more likely to segregate highly penetrant mutations, including families with large numbers of affected individuals or early age at diagnosis. Stronger evidence of linkage in several regions was identified, including a "significant" linkage at 22q12, with a LOD score of 3.57, and five suggestive linkages (1q25, 8q13, 13q14, 16p13, and 17q21) in 269 families with at least five affected members. In addition, four additional suggestive linkages (3p24, 5q35, 11q22, and Xq12) were found in 606 families with mean age at diagnosis of < or = 65 years. Although it is difficult to determine the true statistical significance of these findings, a conservative interpretation of these results would be that if major PC-susceptibility genes do exist, they are most likely located in the regions generating suggestive or significant linkage signals in this large study.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

An interleukin-1 genotype is associated with fatal outcome of meningococcal disease

To determine whether known variants of the interleukin-1 (IL-1) and tumor necrosis factor (TNF) gene families are associated with severe manifestations of meningococcal disease, 276 white patients 4-70 years of age (median, 17 years) were genotyped. All patients had microbiologically proven Neisseria meningitidis infection; 39 died and 237 survived. A significant association (P<.001) was found between fatal outcome and genotype at IL1B (nucleotide position -511). Homozygous individuals, both for the common (1/1) and the rare (2/2) alleles, had increased odds ratios (ORs) for death, compared with heterozygous individuals (1/2): ORs (95% confidence intervals [CIs]) were 3.39 (1.39-8.29) and 7.35 (2.51-21.45), respectively. The mortality rates according to genotype at IL1B (-511) were 18.0% (1/1), 6.1% (1/2), and 32.3% (2/2), compared with 14.2% overall. The composite genotype, consisting of heterozygosity of IL1B (-511) together with homozygosity of the common allele of the IL-1 receptor antagonist gene (IL1RN) at +2018, was significantly associated with survival (P=.018; OR, 7.78 [95% CI, 1. 05-59.05]). There was no association between TNF genotype and fatal outcome. These data suggest that IL-1 genotype influences the severity of meningococcal disease.

Novel genetic association between the corneodesmosin (MHC S) gene and susceptibility to psoriasis

Psoriasis is an inflammatory skin disease of unknown origin, but with a clear genetic component. The strongest genetic association has been found with the major histocompatibility complex (MHC) region, and specifically between susceptibility to familial early onset psoriasis and human leukocyte antigen (HLA)-Cw6. The basis of this association of the HLA-C locus with disease pathogenesis is, however, not clear, and it is possible that other genes, or a combination of genes, in the HLA region are of functional importance. The MHC S gene is expressed specifically in keratinocyte differentiation and, being located 160 kb telomeric of HLA-C, is a plausible candidate gene. We analysed the allelic distribution of two polymorphisms in the MHC S gene (at +619 and +1243) in a case-control association study. We could confirm a significant association between psoriasis and HLA-Cw6 [odds ratio (OR) = 7.75]. No association was found between disease (or any subtypes) and the MHC S gene polymorphism at position +619, despite its close proximity to HLA-C and the strong linkage disequilibrium between the loci. However, a significant trend with the rarer allele at MHC S (+1243) and psoriasis was detected in the overall data set (OR = 2. 66; P = 2 [times] 10(-)9). This effect was most pronounced in the type 1a (early onset) psoriatics (OR = 3.43). Furthermore, homozygosity for the associated allele at MHC S (+1243) increased the risk of disease over that for carriage of HLA-Cw6 alone (OR = 9. 38), suggesting that allele 2 of MHC S (+1243) provides an additional risk in psoriasis susceptibility. The strong association found here, coupled with the biological involvement of the MHC S gene product corneodesmosin in skin physiology, implicates this locus (or a haplotype across HLA-C and MHC S ) in the impaired desquamation characteristic of psoriasis.

Meta-analysis of genome-wide association studies discovers multiple loci for chronic lymphocytic leukemia

Chronic lymphocytic leukemia (CLL) is a common lymphoid malignancy with strong heritability. To further understand the genetic susceptibility for CLL and identify common loci associated with risk, we conducted a meta-analysis of four genome-wide association studies (GWAS) composed of 3,100 cases and 7,667 controls with follow-up replication in 1,958 cases and 5,530 controls. Here we report three new loci at 3p24.1 (rs9880772, EOMES, P=2.55 × 10(-11)), 6p25.2 (rs73718779, SERPINB6, P=1.97 × 10(-8)) and 3q28 (rs9815073, LPP, P=3.62 × 10(-8)), as well as a new independent SNP at the known 2q13 locus (rs9308731, BCL2L11, P=1.00 × 10(-11)) in the combined analysis. We find suggestive evidence (P<5 × 10(-7)) for two additional new loci at 4q24 (rs10028805, BANK1, P=7.19 × 10(-8)) and 3p22.2 (rs1274963, CSRNP1, P=2.12 × 10(-7)). Pathway analyses of new and known CLL loci consistently show a strong role for apoptosis, providing further evidence for the importance of this biological pathway in CLL susceptibility.

Genomewide transmission/disequilibrium testing--consideration of the genotypic relative risks at disease loci

Genomewide association studies are set to become the tool of the future for detection of small-effect genes in complex diseases. It will therefore be necessary to calculate sufficient sample sizes with which to perform them. In this paper I illustrate how to calculate the required number of families for general genotypic relative risks (GRRs). I show the superior sensitivity of the genomewide association study over the standard genomewide affected-sib-pair linkage analysis, for a range of different underlying GRR patterns. I also illustrate the extent of change in the sample sizes that is necessary for a genomewide association analysis depending on the pattern of the GRRs at the disease locus. In many cases, the comparative numbers of families required under different genetic mechanisms vary by several orders of magnitude. These sometimes dramatic differences have important implications for the determination of the size of the collection of samples prior to analysis and for the types of effects that are likely--and unlikely--to be detected by such an analysis.

Assessment of polygenic architecture and risk prediction based on common variants across fourteen cancers

Genome-wide association studies (GWAS) have led to the identification of hundreds of susceptibility loci across cancers, but the impact of further studies remains uncertain. Here we analyse summary-level data from GWAS of European ancestry across fourteen cancer sites to estimate the number of common susceptibility variants (polygenicity) and underlying effect-size distribution. All cancers show a high degree of polygenicity, involving at a minimum of thousands of loci. We project that sample sizes required to explain 80% of GWAS heritability vary from 60,000 cases for testicular to over 1,000,000 cases for lung cancer. The maximum relative risk achievable for subjects at the 99th risk percentile of underlying polygenic risk scores (PRS), compared to average risk, ranges from 12 for testicular to 2.5 for ovarian cancer. We show that PRS have potential for risk stratification for cancers of breast, colon and prostate, but less so for others because of modest heritability and lower incidence.

Principal component analysis for selection of optimal SNP-sets that capture intragenic genetic variation

Candidate gene association studies often utilize one single nucleotide polymorphism (SNP) for analysis, with an initial report typically not being replicated by subsequent studies. The failure to replicate may result from incomplete or poor identification of disease-related variants or haplotypes, possibly due to naive SNP selection. A method for identification of linkage disequilibrium (LD) groups and selection of SNPs that capture sufficient intra-genic genetic diversity is described. We assume all SNPs with minor allele frequency above a pre-determined frequency have been identified. Principal component analysis (PCA) is applied to evaluate multivariate SNP correlations to infer groups of SNPs in LD (LD-groups) and to establish an optimal set of group-tagging SNPs (gtSNPs) that provide the most comprehensive coverage of intra-genic diversity while minimizing the resources necessary to perform an informative association analysis. This PCA method differs from haplotype block (HB) and haplotype-tagging SNP (htSNP) methods, in that an LD-group of SNPs need not be a contiguous DNA fragment. Results of the PCA method compared well with existing htSNP methods while also providing advantages over those methods, including an indication of the optimal number of SNPs needed. Further, evaluation of the method over multiple replicates of simulated data indicated PCA to be a robust method for SNP selection. Our findings suggest that PCA may be a powerful tool for establishing an optimal SNP set that maximizes the amount of genetic variation captured for a candidate gene using a minimal number of SNPs.

Genome-wide association analysis implicates dysregulation of immunity genes in chronic lymphocytic leukaemia

Several chronic lymphocytic leukaemia (CLL) susceptibility loci have been reported; however, much of the heritable risk remains unidentified. Here we perform a meta-analysis of six genome-wide association studies, imputed using a merged reference panel of 1,000 Genomes and UK10K data, totalling 6,200 cases and 17,598 controls after replication. We identify nine risk loci at 1p36.11 (rs34676223, P=5.04 × 10-13), 1q42.13 (rs41271473, P=1.06 × 10-10), 4q24 (rs71597109, P=1.37 × 10-10), 4q35.1 (rs57214277, P=3.69 × 10-8), 6p21.31 (rs3800461, P=1.97 × 10-8), 11q23.2 (rs61904987, P=2.64 × 10-11), 18q21.1 (rs1036935, P=3.27 × 10-8), 19p13.3 (rs7254272, P=4.67 × 10-8) and 22q13.33 (rs140522, P=2.70 × 10-9). These new and established risk loci map to areas of active chromatin and show an over-representation of transcription factor binding for the key determinants of B-cell development and immune response.

Meta association of colorectal cancer confirms risk alleles at 8q24 and 18q21

BACKGROUND

Genome-wide association studies of colorectal cancer (CRC) have identified genetic variants that reproducibly associate with CRC. Associations of 12 single nucleotide polymorphisms at 8q24, 9p24, and 18q21 (SMAD7) and CRC were investigated in a three-center collaborative study including two U.K. case-control cohorts (Sheffield and Leeds) and a U.S. case-control study of CRC cases from high-risk Utah pedigrees.

METHODS

Our combined resource included 1,092 CRC case subjects and 1,060 age- and sex-matched controls. Meta statistics and Monte Carlo significance testing using Genie software provided a valid combined analysis of our mixed independent and related case-control resource. We also evaluated whether these associations differed by sex, age at diagnosis, family history, or tumor site.

RESULTS

At 8q24, we observed two independent significant associations at single nucleotide polymorphisms located in two different risk regions of 8q24: rs6983267 in region 3 [P(trend) = 0.01; per allele odds ratio (OR), 1.17; 95% confidence intervals (95% CI), 1.03-1.32] and rs10090154 in region 5 (P(trend) = 0.05; per allele OR, 1.24; 95% CI, 1.01-1.51). At 18q21, associations were observed in distal colon tumors but not in proximal or rectal cancers: rs4939827 (P(trend) = 0.007; per allele OR, 0.77; 95% CI, 0.64-0.93; case-case p(diff) = 0.03) and rs12953717 (P(trend) = 0.01; per allele OR, 1.27; 95% CI, 1.06-1.52). We were unable to detect any associations at 9p24 with CRC.

CONCLUSIONS

Our investigation confirms that variants across multiple risk regions of 8q24 are associated with CRC, and that associations at 18q21 differ by tumor site.

Genome-wide linkage analyses of extended Utah pedigrees identifies loci that influence recurrent, early-onset major depression and anxiety disorders

Major depressive disorder (MDD) is a common, clinically heterogeneous disorder often found comorbid with other disorders. We studied recurrent, early-onset MDD (MDD-RE) and anxiety disorders in combination to define powerful phenotypes for genetic study. We used 87 large, extended Utah pedigrees to investigate linkage to 3 phenotypes: "MDD-RE;" "MDD-RE or anxiety;" and "MDD-RE and anxiety;" where in the latter definition the disorders must appear comorbid within an individual. Pedigrees ranged in size from 2 to 6 generations and contained 3 to 42 individuals affected with MDD or anxiety (718 total). In primary analyses, we identified three regions with at least suggestive genome-wide evidence for linkage on chromosomes 3centr, 7p, and 18q. Both 7p and 18q are replication findings for related phenotypes. The best linkage evidence was for a novel locus at 3p12.3-q12.3 (LOD = 3.88, "MDD-RE or anxiety") and 18q21.33-q22.2 (LOD = 3.75, "MDD-RE and anxiety"), a well-established susceptibility locus for bipolar disorder. In our secondary sex-specific analyses, we identified two further regions of interest on chromosomes 4q and 15q. Using linked pedigrees, we localized 3centr and 18q to 9.8 and 12.2 cM, respectively, with potential for further localization with the addition of markers in specific pedigrees. Our success in replication and novel locus identification illustrates the utility of large extended pedigrees for common disorders, such as MDD. Further, it supports the hypothesis that MDD and anxiety disorders have over-lapping genetic etiologies and suggests that comorbid diagnoses may be useful in defining more genetically homogeneous forms of MDD for linkage mapping.

Classification tree analysis: a statistical tool to investigate risk factor interactions with an example for colon cancer (United States)

OBJECTIVE

Classification tree analysis is a potentially powerful tool for investigating multilevel interactions. Within the context of colon cancer etiology it may help identify disease pathways and evaluate important interactions of risk factors.

METHODS

We apply classification tree analysis as a statistical method to investigate interactions of risk factors for colon cancer. We use data collected from a population-based case-control study of newly diagnosed cases of colon cancer (N = 4403 cases and controls).

RESULTS

Our results indicate that, as expected, there are many factors that influence colon cancer risk, and that they interact on many levels. We find that the most important factor is the utilization of aspirin and/or non-steroidal anti-inflammatory drugs (NSAID), with those taking this medication having lower risk. Family history appears as a level two modifying factor when NSAID are not used, whereas Western diet is the second factor when NSAID are taken. The final tree has six levels, contains several modifying factors and correctly classifies case or control status for 60.8% (95% CI 59.4-62.2) of all individuals.

CONCLUSIONS

Our results suggest that risk factors work together to determine disease risk. By accounting for interactions between risk factors we become better able to dissect disease pathways and determine those risk factors that increase susceptibility to disease. Our results highlight the importance of designing studies so that interactions can be addressed.

Shared genomic segment analysis. Mapping disease predisposition genes in extended pedigrees using SNP genotype assays

We examine the utility of high density genotype assays for predisposition gene localization using extended pedigrees. Results for the distribution of the number and length of genomic segments shared identical by descent among relatives previously derived in the context of genomic mismatch scanning are reviewed in the context of dense single nucleotide polymorphism maps. We use long runs of loci at which cases share a common allele identically by state to localize hypothesized predisposition genes. The distribution of such runs under the hypothesis of no genetic effect is evaluated by simulation. Methods are illustrated by analysis of an extended prostate cancer pedigree previously reported to show significant linkage to chromosome 1p23. Our analysis establishes that runs of simple single locus statistics can be powerful, tractable and robust for finding DNA shared between relatives, and that extended pedigrees offer powerful designs for gene detection based on these statistics.

Germline Lysine-Specific Demethylase 1 (LSD1/KDM1A) Mutations Confer Susceptibility to Multiple Myeloma

Given the frequent and largely incurable occurrence of multiple myeloma, identification of germline genetic mutations that predispose cells to multiple myeloma may provide insight into disease etiology and the developmental mechanisms of its cell of origin, the plasma cell (PC). Here, we identified familial and early-onset multiple myeloma kindreds with truncating mutations in lysine-specific demethylase 1 (LSD1/KDM1A), an epigenetic transcriptional repressor that primarily demethylates histone H3 on lysine 4 and regulates hematopoietic stem cell self-renewal. In addition, we found higher rates of germline truncating and predicted deleterious missense KDM1A mutations in patients with multiple myeloma unselected for family history compared with controls. Both monoclonal gammopathy of undetermined significance (MGUS) and multiple myeloma cells have significantly lower KDM1A transcript levels compared with normal PCs. Transcriptome analysis of multiple myeloma cells from KDM1A mutation carriers shows enrichment of pathways and MYC target genes previously associated with myeloma pathogenesis. In mice, antigen challenge followed by pharmacologic inhibition of KDM1A promoted PC expansion, enhanced secondary immune response, elicited appearance of serum paraprotein, and mediated upregulation of MYC transcriptional targets. These changes are consistent with the development of MGUS. Collectively, our findings show that KDM1A is the first autosomal-dominant multiple myeloma germline predisposition gene providing new insights into its mechanistic roles as a tumor suppressor during post-germinal center B-cell differentiation.Significance: KDM1A is the first germline autosomal dominant predisposition gene identified in multiple myeloma and provides new insights into multiple myeloma etiology and the mechanistic role of KDM1A as a tumor suppressor during post-germinal center B-cell differentiation. Cancer Res; 78(10); 2747-59. ©2018 AACR.

Confirmation of chromosome 7q11 locus for predisposition to intracranial aneurysm

A significant linkage of intracranial aneurysm (IA) has recently been reported to chromosomal region 7q11 (MLS=3.22) in a genomic search of 85 Japanese nuclear families with at least two affected siblings (104 sib pairs). This region contains the elastin gene (ELN, OMIM 130160), which is a functional candidate gene for IA. We have replicated this finding through linkage analyses in 13 extended pedigrees from Utah, comprising 39 IA cases. We genotyped three markers flanking ELN and performed two-point and multipoint parametric analyses, employing simple dominant and recessive models. Analyses utilizing a recessive affecteds-only model yielded significant confirmation of linkage to the region (best evidence, multipoint TLOD=2.34, at D7S2421, corrected P=0.001). This study is the first to confirm the linkage of the 7q11 locus for IA.

Pathology of Tumors Associated With Pathogenic Germline Variants in 9 Breast Cancer Susceptibility Genes

IMPORTANCE

Rare germline genetic variants in several genes are associated with increased breast cancer (BC) risk, but their precise contributions to different disease subtypes are unclear. This information is relevant to guidelines for gene panel testing and risk prediction.

OBJECTIVE

To characterize tumors associated with BC susceptibility genes in large-scale population- or hospital-based studies.

DESIGN, SETTING, AND PARTICIPANTS

The multicenter, international case-control analysis of the BRIDGES study included 42 680 patients and 46 387 control participants, comprising women aged 18 to 79 years who were sampled independently of family history from 38 studies. Studies were conducted between 1991 and 2016. Sequencing and analysis took place between 2016 and 2021.

EXPOSURES

Protein-truncating variants and likely pathogenic missense variants in ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53.

MAIN OUTCOMES AND MEASURES

The intrinsic-like BC subtypes as defined by estrogen receptor, progesterone receptor, and ERBB2 (formerly known as HER2) status, and tumor grade; morphology; size; stage; lymph node involvement; subtype-specific odds ratios (ORs) for carrying protein-truncating variants and pathogenic missense variants in the 9 BC susceptibility genes.

RESULTS

The mean (SD) ages at interview (control participants) and diagnosis (cases) were 55.1 (11.9) and 55.8 (10.6) years, respectively; all participants were of European or East Asian ethnicity. There was substantial heterogeneity in the distribution of intrinsic subtypes by gene. RAD51C, RAD51D, and BARD1 variants were associated mainly with triple-negative disease (OR, 6.19 [95% CI, 3.17-12.12]; OR, 6.19 [95% CI, 2.99-12.79]; and OR, 10.05 [95% CI, 5.27-19.19], respectively). CHEK2 variants were associated with all subtypes (with ORs ranging from 2.21-3.17) except for triple-negative disease. For ATM variants, the association was strongest for the hormone receptor (HR)+ERBB2- high-grade subtype (OR, 4.99; 95% CI, 3.68-6.76). BRCA1 was associated with increased risk of all subtypes, but the ORs varied widely, being highest for triple-negative disease (OR, 55.32; 95% CI, 40.51-75.55). BRCA2 and PALB2 variants were also associated with triple-negative disease. TP53 variants were most strongly associated with HR+ERBB2+ and HR-ERBB2+ subtypes. Tumors occurring in pathogenic variant carriers were of higher grade. For most genes and subtypes, a decline in ORs was observed with increasing age. Together, the 9 genes were associated with 27.3% of all triple-negative tumors in women 40 years or younger.

CONCLUSIONS AND RELEVANCE

The results of this case-control study suggest that variants in the 9 BC risk genes differ substantially in their associated pathology but are generally associated with triple-negative and/or high-grade disease. Knowing the age and tumor subtype distributions associated with individual BC genes can potentially aid guidelines for gene panel testing, risk prediction, and variant classification and guide targeted screening strategies.

Association of the interleukin 1 receptor antagonist gene with ulcerative colitis in Northern European Caucasians

BACKGROUND AND AIMS

An association between the allele 2 of the interleukin 1 receptor antagonist gene variable number tandem repeats polymorphism in intron 2 and ulcerative colitis was first reported in 1994. Subsequent studies in Caucasian Northern European patients have not confirmed this, although trends towards an association were observed. The lack of statistical significance could reflect inadequate power. In this study the association was reassessed in a large independent set of well characterised Caucasian patients and a meta-analysis of reported patient series was performed.

PATIENTS AND METHODS

A total of 320 patients with endoscopically and histologically confirmed ulcerative colitis (124 pancolitis, 196 left sided and distal disease) and 827 ethnically matched controls were genotyped at polymorphic sites in the interleukin 1 receptor antagonist gene. Carriage rates were compared using chi(2) statistics. A meta-analysis of this and seven previous studies in North European Caucasian patients was performed using the Mantel-Haenszel chi(2) test.

RESULTS

Patients had a significantly increased carriage rate of allele 2 compared with controls (52% v 45%; odds ratio 1.3 (95% confidence interval (CI) 1.01-1.7); p=0.04). The allele 2 carriage rate was highest in extensive colitis (carriage rate 56%; odds ratio 1.5 (95% CI 1.1-2.3) p=0.02) and in individuals who had undergone colectomy (carriage rate 55%; odds ratio 1.5 (95% CI 0.95-2.4); p=0.08). Meta-analysis of all eight studies showed a significant association between carriage of allele 2 and ulcerative colitis (odds ratio 1.23 (95% CI 1.04-1.45); p=0.01).

CONCLUSIONS

The association of the interleukin 1 receptor antagonist gene polymorphism with ulcerative colitis is confirmed. The association is minor and confers only a small risk to an individual but will contribute a high attributable risk in a population due to the high allelic frequency. Accurate phenotypic characterisation defines more homogeneous subsets of patients, such as those with extensive disease, in whom the association is greater.