Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

A Rare Variant in MDH2 (rs111879470) Is Associated with Predisposition to Recurrent Breast Cancer in an Extended High-Risk Pedigree.

A significant fraction of breast cancer recurs, with lethal outcome, but specific genetic variants responsible have yet to be identified. Five cousin pairs with recurrent breast cancer from pedigrees with a statistical excess of recurrent breast cancer were sequenced to identify rare, shared candidate predisposition variants. The candidates were tested for association with breast cancer risk with UKBiobank data. Additional breast cancer cases were assayed for a subset of candidate variants to test for co-segregation. Three-dimensional protein structure prediction methods were used to investigate how the mutation under consideration is predicted to change structural and electrostatic properties in the mutated protein. One hundred and eighty-one rare candidate predisposition variants were shared in at least one cousin pair from a high-risk pedigree. A rare variant in MDH2 was found to segregate with breast-cancer-affected relatives in one extended pedigree. MDH2 is an estrogen-stimulated gene encoding the protein malate dehydrogenase, which catalyzes the reversible oxidation of malate to oxaloacetate. The molecular simulation results strongly suggest that the mutation changes the NAD+ binding pocket electrostatics of MDH2. This small sequencing study, using a powerful approach based on recurrent breast cancer cases from high-risk pedigrees, identified a set of strong candidate variants for inherited predisposition for breast cancer recurrence, including MDH2, which should be pursued in other resources.

Germline EGFR Mutations and Familial Lung Cancer.

PURPOSE: The genomic underpinnings of inherited lung cancer risk are poorly understood. This prospective study characterized the clinical phenotype of patients and families with germline EGFR pathogenic variants (PVs). METHODS: The Investigating Hereditary Risk from T790M study (ClinicalTrials.gov identifier: NCT01754025) enrolled patients with lung cancer whose tumor profiling harbored possible germline EGFR PVs and their relatives, either in person or remotely, providing germline testing and follow-up. RESULTS: A total of 141 participants were enrolled over a 5-year period, 100 (71%) remotely. Based upon previous genotyping, 116 participants from 59 kindreds were tested for EGFR T790M, demonstrating a pattern of Mendelian inheritance with variable lung cancer penetrance. In confirmed or obligate carriers of a germline EGFR PV from 39 different kindreds, 50/91 (55%) were affected with lung cancer with 34/65 (52%) diagnosed by age 60 years. Somatic testing of lung cancers in carriers revealed that 35 of 37 (95%) had an EGFR driver comutation. Among 36 germline carriers without a cancer diagnosis, 15 had computed tomography (CT) imaging and nine had lung nodules, including a 28-year-old with >10 lung nodules. Given geographic enrichment of germline EGFR T790M in the southeast United States, genome-wide haplotyping of 46 germline carriers was performed and identified a 4.1-Mb haplotype shared by 41 (89%), estimated to originate 223-279 years ago. CONCLUSION: To our knowledge, this is the first prospective description of familial EGFR-mutant lung cancer, identifying a recent founder germline EGFR T790M variant enriched in the Southeast United States. The high prevalence of EGFR-driver lung adenocarcinomas and lung nodules in germline carriers supports effort to identify affected patients and family members for investigation of CT-based screening for these high-risk individuals.

High-Risk Pedigree Study Identifies LRBA (rs62346982) as a Likely Predisposition Variant for Prostate Cancer.

There is evidence for contribution of inherited factors to prostate cancer, and more specifically to lethal prostate cancer, but few responsible genes/variants have been identified. We examined genetic sequence data for 51 affected cousin pairs who each died from prostate cancer and who were members of high-risk prostate cancer pedigrees in order to identify rare variants shared by the cousins as candidate predisposition variants. Candidate variants were tested for association with prostate cancer risk in UK Biobank data. Candidate variants were also assayed in 1195 additional sampled Utah prostate cancer cases. We used 3D protein structure prediction methods to analyze structural changes and provide insights into mechanisms of pathogenicity. Almost 4000 rare (<0.005) variants were identified as shared in the 51 affected cousin pairs. One candidate variant was also significantly associated with prostate cancer risk among the 840 variants with data in UK Biobank, in the gene LRBA (p = 3.2 × 10-5; OR = 2.09). The rare risk variant in LRBA was observed to segregate in five pedigrees. The overall predicted structures of the mutant protein do not show any significant overall changes upon mutation, but the mutated structure loses a helical structure for the two residues after the mutation. This unique analysis of closely related individuals with lethal prostate cancer, who were members of high-risk prostate cancer pedigrees, has identified a strong set of candidate predisposition variants which should be pursued in independent studies. Validation data for a subset of the candidates identified are presented, with strong evidence for a rare variant in LRBA.

A rare FGF5 candidate variant (rs112475347) for predisposition to nonsquamous, nonsmall-cell lung cancer.

A unique approach with rare resources was used to identify candidate variants predisposing to familial nonsquamous nonsmall-cell lung cancers (NSNSCLC). We analyzed sequence data from NSNSCLC-affected cousin pairs belonging to high-risk lung cancer pedigrees identified in a genealogy of Utah linked to statewide cancer records to identify rare, shared candidate predisposition variants. Variants were tested for association with lung cancer risk in UK Biobank. Evidence for linkage with lung cancer was also reviewed in families from the Genetic Epidemiology of Lung Cancer Consortium. Protein prediction modeling compared the mutation with reference. We sequenced NSNSCLC-affected cousin pairs from eight high-risk lung cancer pedigrees and identified 66 rare candidate variants shared in the cousin pairs. One variant in the FGF5 gene also showed significant association with lung cancer in UKBiobank. This variant was observed in 3/163 additional sampled Utah lung cancer cases, 2 of whom were related in another independent pedigree. Modeling of the predicted protein predicted a second binding site for SO4 that may indicate binding differences. This unique study identified multiple candidate predisposition variants for NSNSCLC, including a rare variant in FGF5 that was significantly associated with lung cancer risk and that segregated with lung cancer in the two pedigrees in which it was observed. FGF5 is an oncogenic factor in several human cancers, and the mutation found here (W81C) changes the binding ability of heparan sulfate to FGF5, which might lead to its deregulation. These results support FGF5 as a potential NSNSCLC predisposition gene and present additional candidate predisposition variants.

Comorbidity and Cancer Disease Rates among Those at High-Risk for Alzheimer's Disease: A Population Database Analysis.

(1) Importance: Alzheimer's disease (AD) is complex and only partially understood. Analyzing the relationship between other more treatable or preventable diseases and AD may help in the prevention and the eventual development of treatments for AD. Risk estimation in a high-risk population, rather than a population already affected with AD, may reduce some bias in risk estimates. (2) Objective: To examine the rates of various comorbidities and cancers in individuals at high-risk for AD, but without a clinical diagnosis, relative to individuals from the same population with normal AD risk. (3) Design, Setting, and Participants: We conducted a study using data from the Utah Population Database (UPDB). The UPDB contains linked data from the Utah Cancer Registry, Utah death certificates, the Intermountain Health patient population, and the University of Utah Health patient population. Subjects were selected based on the availability of ancestral data, linked health information, and self-reported biometrics. (4) Results: In total, 75,877 participants who were estimated to be at high risk for AD based on family history, but who did not have an active AD diagnosis, were analyzed. A lower incidence of diabetes (RR = 0.95, 95% CI [0.92,0.97], p < 0.001), hypertension (RR = 0.97, 95% CI [0.95,0.99], p < 0.001), and heart disease (RR = 0.95, 95% CI [0.93,0.98], p < 0.001) was found. There was no difference in rates of cerebrovascular disease or other forms of dementia. Of the 15 types of cancer analyzed: breast (RR = 1.23, 95% CI [1.16, 1.30], p < 0.001); colorectal (RR = 1.30, 95% CI [1.21, 1.39], p < 0.001); kidney (RR = 1.49, 95% CI (1.29, 1.72), p < 0.001); lung (RR = 1.25, 95% CI [1.13, 1.37], p < 0.001); non-Hodgkin's Lymphoma (RR = 1.29, 95% CI [1.15, 1.44], p < 0.001); pancreas (RR = 1.34, 95% CI [1.16, 1.55], p < 0.001); stomach (RR = 1.59, 95% CI [1.36, 1.86], p < 0.001); and bladder (RR = 1.40, 95% CI [1.25, 1.56], p < 0.001), cancers were observed in significant excess among individuals at high-risk for AD after correction for multiple testing. (5) Conclusions and Relevance: Since age is the greatest risk factor for the development of AD, individuals who reach more advanced ages are at increased risk of developing AD. Consistent with this, people with fewer comorbidities earlier in life are more likely to reach an age where AD becomes a larger risk. Our findings show that individuals at high risk for AD have a decreased incidence of various other diseases. This is further supported by our finding that our high-risk group was also found to have an increased incidence of various cancers, which also increase in risk with age. There is the possibility that a more meaningful or etiological relationship exists among these various comorbidities. Further research into the etiological relationship between AD and these comorbidities may elucidate these possible interactions.

Distinct germline genetic susceptibility profiles identified for common non-Hodgkin lymphoma subtypes.

Lymphoma risk is elevated for relatives with common non-Hodgkin lymphoma (NHL) subtypes, suggesting shared genetic susceptibility across subtypes. To evaluate the extent of mutual heritability among NHL subtypes and discover novel loci shared among subtypes, we analyzed data from eight genome-wide association studies within the InterLymph Consortium, including 10,629 cases and 9505 controls. We utilized Association analysis based on SubSETs (ASSET) to discover loci for subsets of NHL subtypes and evaluated shared heritability across the genome using Genome-wide Complex Trait Analysis (GCTA) and polygenic risk scores. We discovered 17 genome-wide significant loci (P < 5 × 10-8) for subsets of NHL subtypes, including a novel locus at 10q23.33 (HHEX) (P = 3.27 × 10-9). Most subset associations were driven primarily by only one subtype. Genome-wide genetic correlations between pairs of subtypes varied broadly from 0.20 to 0.86, suggesting substantial heterogeneity in the extent of shared heritability among subtypes. Polygenic risk score analyses of established loci for different lymphoid malignancies identified strong associations with some NHL subtypes (P < 5 × 10-8), but weak or null associations with others. Although our analyses suggest partially shared heritability and biological pathways, they reveal substantial heterogeneity among NHL subtypes with each having its own distinct germline genetic architecture.

Is There an Inherited Contribution to Risk for Sporadic Unilateral Vestibular Schwannoma? Evidence of Familial Clustering.

OBJECT: Unlike the autosomal dominant inheritance of neurofibromatosis 2, there are no known inherited risk factors for sporadic, unilateral vestibular schwannoma (VS), which comprise most VS cases. The authors tested a hypothesis positing a genetic contribution to predisposition to these lesions by analyzing familial clustering of cases. METHODS: Familial clustering of individuals with unilateral VS was analyzed in two independent genealogical resources with linked diagnosis data: the Veterans Health Administration Genealogy Resource and the Utah Population Database. Tests for excess relatedness, estimation of relative risks (RRs) in close and distant relatives, and identification of pedigrees with a significant excess of unilateral VS among descendants were performed. RESULTS: The average pairwise relatedness of the Veterans Health Administration Genealogy Resource VS cases significantly exceeded the expected relatedness ( p = 0.016), even when close relationships were ignored ( p = 0.002). RR for third- and fifth-degree relatives developing VS were significantly elevated (RR, 60.83; p = 0.0005; 95% confidence interval [CI], 7.37-219.73) and (RR, 11.88; p = 0.013; 95% CI, 1.44-42.90), respectively. No VS-affected first-, second-, or fourth-degree relatives were observed. In the Utah Population Database population, no first- or second-degree relatives with VS were observed. RR for fifth-degree relatives developing VS was significantly elevated (RR, 2.23; p = 0.009; 95% CI, 1.15-3.90). CONCLUSION: These results provide strong evidence for an inherited predisposition to sporadic, unilateral VS. This study exhibits the value of genealogical resources with linked medical data for examining hypotheses regarding inherited predisposition. The high-risk unilateral VS pedigrees identified in two independent resources provide a powerful means of pursuing predisposition gene identification.

The Effect of Sex Hormone Deficiency on the Incidence of Rotator Cuff Repair: Analysis of a Large Insurance Database.

BACKGROUND: The purpose of the present study was to analyze the association between sex hormone deficiency and rotator cuff repair (RCR) with use of data from a large United States insurance database. METHODS: A retrospective analysis of insured subjects from the Truven Health MarketScan database was conducted, collecting data for RCR cases as well as controls matched for age, sex, and years in the database. Multivariable logistic regression models adjusted for matching variables were utilized to compare RCR status with estrogen deficiency status and testosterone deficiency status. These associations were confirmed with use of data from the Veterans Genealogy Project database, with which the relative risk of RCR was estimated for patients with and without sex hormone deficiency. RESULTS: The odds of RCR for female patients with estrogen deficiency were 48% higher (odds ratio, 1.48; 95% confidence interval, 1.44 to 1.51; p < 0.001) than for those without estrogen deficiency. The odds of RCR for males with testosterone deficiency were 89% higher (odds ratio, 1.89; 95% confidence interval, 1.82 to 1.96; p < 0.001) than for those without testosterone deficiency. Within the Veterans Genealogy Project database, the relative risk of estrogen deficiency among RCR patients was 2.58 (95% confidence interval, 2.15 to 3.06; p < 0.001) and the relative risk of testosterone deficiency was 3.05 (95% confidence interval, 2.67 to 3.47; p < 0.001). CONCLUSIONS: Sex hormone deficiency was significantly associated with RCR. Future prospective studies will be necessary to understand the pathophysiology of rotator cuff disease as it relates to sex hormones. LEVEL OF EVIDENCE: Prognostic Level IV. See Instructions for Authors for a complete description of levels of evidence.

A whole-exome case-control association study to characterize the contribution of rare coding variation to pancreatic cancer risk.

Pancreatic cancer is a deadly disease that accounts for approximately 5% of cancer deaths worldwide, with a dismal 5-year survival rate of 10%. Known genetic risk factors explain only a modest proportion of the heritable risk of pancreatic cancer. We conducted a whole-exome case-control sequencing study in 1,591 pancreatic cancer cases and 2,134 cancer-free controls of European ancestry. In our gene-based analysis, ATM ranked first, with a genome-wide significant p value of 1 × 10-8. The odds ratio for protein-truncating variants in ATM was 24, which is substantially higher than prior estimates, although ours includes a broad 95% confidence interval (4.0-1000). SIK3 was the second highest ranking gene (p = 3.84 × 10-6, false discovery rate or FDR = 0.032). We observed nominally significant association signals in several genes of a priori interest, including BRCA2 (p = 4.3 × 10-4), STK11 (p = 0.003), PALB2 (p = 0.019), and TP53 (p = 0.037), and reported risk estimates for known pathogenic variants and variants of uncertain significance (VUS) in these genes. The rare variants in established susceptibility genes explain approximately 24% of log familial relative risk, which is comparable to the contribution from established common susceptibility variants (17%). In conclusion, this study provides new insights into the genetic susceptibility of pancreatic cancer, refining rare variant risk estimates in known pancreatic cancer susceptibility genes and identifying SIK3 as a novel candidate susceptibility gene. This study highlights the prominent importance of ATM truncating variants and the underappreciated role of VUS in pancreatic cancer etiology.

Early-onset colorectal cancer risk extends to second- and third-degree relatives.

INTRODUCTION: Family history is a risk factor for colorectal cancer (CRC), however whether family history also contributes to non-syndromic early-onset CRC is unknown. METHODS: We estimated risk to first-, second-, and third-degree relatives of early-onset CRC cases in the Utah Pedigree Database. RESULTS: We observed elevated risks beyond RR = 2.0 for early-onset CRC among first- and second-degree relatives of early-onset CRC cases, with RRs of 6.0 and 3.1, respectively. DISCUSSION: Relatives of early-onset CRC cases are at higher risk of both early-onset CRC and CRC at any age, and the location is not necessarily similar to the affected relative.

A Rare Variant in ERF (rs144812092) Predisposes to Prostate and Bladder Cancers in an Extended Pedigree.

Pairs of related bladder cancer cases who belong to pedigrees with an excess of bladder cancer were sequenced to identify rare, shared variants as candidate predisposition variants. Candidate variants were tested for association with bladder cancer risk. A validated variant was assayed for segregation to other related cancer cases, and the predicted protein structure of this variant was analyzed. This study of affected bladder cancer relative pairs from high-risk pedigrees identified 152 bladder cancer predisposition candidate variants. One variant in ERF (ETS Repressing Factor) was significantly associated with bladder cancer risk in an independent population, was observed to segregate with bladder and prostate cancer in relatives, and showed evidence for altering the function of the associated protein. This finding of a rare variant in ERF that is strongly associated with bladder and prostate cancer risk in an extended pedigree both validates ERF as a cancer predisposition gene and shows the continuing value of analyzing affected members of high-risk pedigrees to identify and validate rare cancer predisposition variants.

Differential methylation of G-protein coupled receptor signaling genes in gastrointestinal neuroendocrine tumors.

Neuroendocrine tumors (NETs) of the small intestine undergo large chromosomal and methylation changes. The objective of this study was to identify methylation differences in NETs and consider how the differentially methylated genes may impact patient survival. Genome-wide methylation and chromosomal copy number variation (CNV) of NETs from the small intestine and appendix were measured. Tumors were divided into three molecular subtypes according to CNV results: chromosome 18 loss (18LOH), Multiple CNV, and No CNV. Comparison of 18LOH tumors with MultiCNV and NoCNV tumors identified 901 differentially methylated genes. Genes from the G-protein coupled receptor (GPCR) pathways are statistically overrepresented in the differentially methylated genes. One of the highlighted genes from the GPCR pathway is somatostatin (SST), a clinical target for NETs. Patient survival based on low versus high methylation in all samples identified four significant genes (p < 0.05) OR2S2, SMILR, RNU6-653P, and AC010543.1. Within the 18LOH molecular subtype tumors, survival differences were identified in high versus low methylation of 24 genes. The most significant is TRHR (p < 0.01), a GPCR with multiple FDA-approved drugs. By separating NETs into different molecular subtypes based on chromosomal changes, we find that multiple GPCRs and their ligands appear to be regulated through methylation and correlated with survival. These results suggest opportunities for better treatment strategies for NETs based on molecular features.

An intronic variant in the CELF4 gene is associated with risk for colorectal cancer.

BACKGROUND: Germline predisposition variants associated with colorectal cancer (CRC) have been identified but all are not yet identified. We sought to identify the responsible predisposition germline variant in an extended high-risk CRC pedigree that exhibited evidence of linkage to the 18q12.2 region (TLOD = +2.81). METHODS: DNA from two distantly related carriers of the hypothesized predisposition haplotype on 18q12.2 was sequenced to identify candidate variants. The candidate rare variants shared by the related sequenced subjects were screened in 3,094 CRC cases and 5x population-matched controls from UKBiobank to test for association. Further segregation of the variant was tested via Taqman assay in other sampled individuals in the pedigree. RESULTS: Analysis of whole genome sequence data for the two related hypothesized predisposition haplotype carriers, restricted to the shared haplotype boundaries, identified multiple (n = 6) rare candidate non-coding variants that were tested for association with CRC risk in UKBiobank. A rare intronic variant ofCELF4 gene, rs568643870, was significantly associated with CRC (p = 0.004, OR = 5.0), and segregated with CRC in other members of the linked pedigree. CONCLUSION: Evidence of segregation in a high-risk pedigree, case-control association in an external dataset, and identification of additional CRC-affected carriers in the linked pedigree support a role for a rareCELF4 intronic variant in CRC risk.

KLK3 SNP-SNP interactions for prediction of prostate cancer aggressiveness.

Risk classification for prostate cancer (PCa) aggressiveness and underlying mechanisms remain inadequate. Interactions between single nucleotide polymorphisms (SNPs) may provide a solution to fill these gaps. To identify SNP-SNP interactions in the four pathways (the angiogenesis-, mitochondria-, miRNA-, and androgen metabolism-related pathways) associated with PCa aggressiveness, we tested 8587 SNPs for 20,729 cases from the PCa consortium. We identified 3 KLK3 SNPs, and 1083 (P < 3.5 × 10-9) and 3145 (P < 1 × 10-5) SNP-SNP interaction pairs significantly associated with PCa aggressiveness. These SNP pairs associated with PCa aggressiveness were more significant than each of their constituent SNP individual effects. The majority (98.6%) of the 3145 pairs involved KLK3. The 3 most common gene-gene interactions were KLK3-COL4A1:COL4A2, KLK3-CDH13, and KLK3-TGFBR3. Predictions from the SNP interaction-based polygenic risk score based on 24 SNP pairs are promising. The prevalence of PCa aggressiveness was 49.8%, 21.9%, and 7.0% for the PCa cases from our cohort with the top 1%, middle 50%, and bottom 1% risk profiles. Potential biological functions of the identified KLK3 SNP-SNP interactions were supported by gene expression and protein-protein interaction results. Our findings suggest KLK3 SNP interactions may play an important role in PCa aggressiveness.

Rare Germline Variants in ATM Predispose to Prostate Cancer: A PRACTICAL Consortium Study.

BACKGROUND: Germline ATM mutations are suggested to contribute to predisposition to prostate cancer (PrCa). Previous studies have had inadequate power to estimate variant effect sizes. OBJECTIVE: To precisely estimate the contribution of germline ATM mutations to PrCa risk. DESIGN, SETTING, AND PARTICIPANTS: We analysed next-generation sequencing data from 13 PRACTICAL study groups comprising 5560 cases and 3353 controls of European ancestry. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Variant Call Format files were harmonised, annotated for rare ATM variants, and classified as tier 1 (likely pathogenic) or tier 2 (potentially deleterious). Associations with overall PrCa risk and clinical subtypes were estimated. RESULTS AND LIMITATIONS: PrCa risk was higher in carriers of a tier 1 germline ATM variant, with an overall odds ratio (OR) of 4.4 (95% confidence interval [CI]: 2.0-9.5). There was also evidence that PrCa cases with younger age at diagnosis (<65 yr) had elevated tier 1 variant frequencies (pdifference = 0.04). Tier 2 variants were also associated with PrCa risk, with an OR of 1.4 (95% CI: 1.1-1.7). CONCLUSIONS: Carriers of pathogenic ATM variants have an elevated risk of developing PrCa and are at an increased risk for earlier-onset disease presentation. These results provide information for counselling of men and their families. PATIENT SUMMARY: In this study, we estimated that men who inherit a likely pathogenic mutation in the ATM gene had an approximately a fourfold risk of developing prostate cancer. In addition, they are likely to develop the disease earlier.

Two-stage Study of Familial Prostate Cancer by Whole-exome Sequencing and Custom Capture Identifies 10 Novel Genes Associated with the Risk of Prostate Cancer.

BACKGROUND: Family history of prostate cancer (PCa) is a well-known risk factor, and both common and rare genetic variants are associated with the disease. OBJECTIVE: To detect new genetic variants associated with PCa, capitalizing on the role of family history and more aggressive PCa. DESIGN, SETTING, AND PARTICIPANTS: A two-stage design was used. In stage one, whole-exome sequencing was used to identify potential risk alleles among affected men with a strong family history of disease or with more aggressive disease (491 cases and 429 controls). Aggressive disease was based on a sum of scores for Gleason score, node status, metastasis, tumor stage, prostate-specific antigen at diagnosis, systemic recurrence, and time to PCa death. Genes identified in stage one were screened in stage two using a custom-capture design in an independent set of 2917 cases and 1899 controls. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Frequencies of genetic variants (singly or jointly in a gene) were compared between cases and controls. RESULTS AND LIMITATIONS: Eleven genes previously reported to be associated with PCa were detected (ATM, BRCA2, HOXB13, FAM111A, EMSY, HNF1B, KLK3, MSMB, PCAT1, PRSS3, and TERT), as well as an additional 10 novel genes (PABPC1, QK1, FAM114A1, MUC6, MYCBP2, RAPGEF4, RNASEH2B, ULK4, XPO7, and THAP3). Of these 10 novel genes, all but PABPC1 and ULK4 were primarily associated with the risk of aggressive PCa. CONCLUSIONS: Our approach demonstrates the advantage of gene sequencing in the search for genetic variants associated with PCa and the benefits of sampling patients with a strong family history of disease or an aggressive form of disease. PATIENT SUMMARY: Multiple genes are associated with prostate cancer (PCa) among men with a strong family history of this disease or among men with an aggressive form of PCa.

Genome-wide analysis of high-risk primary brain cancer pedigrees identifies PDXDC1 as a candidate brain cancer predisposition gene.

BACKGROUND: There is evidence for an inherited contribution to primary brain cancer. Linkage analysis of high-risk brain cancer pedigrees has identified candidate regions of interest in which brain cancer predisposition genes are likely to reside. METHODS: Genome-wide linkage analysis was performed in a unique set of 11 informative, extended, high-risk primary brain cancer pedigrees identified in a population genealogy database, which include from 2 to 6 sampled, related primary brain cancer cases. Access to formalin-fixed paraffin embedded tissue samples archived in a biorepository allowed analysis of extended pedigrees. RESULTS: Individual high-risk pedigrees were singly informative for linkage at multiple regions. Suggestive evidence for linkage was observed on chromosomes 2, 3, 14, and 16. The chromosome 16 region in particular contains a promising candidate gene, pyridoxal-dependent decarboxylase domain-containing 1 (PDXDC1), with prior evidence for involvement with glioblastoma from other previously reported experimental settings, and contains the lead single nucleotide polymorphism (rs3198697) from the linkage analysis of the chromosome 16 region. CONCLUSIONS: Pedigrees with a statistical excess of primary brain cancers have been identified in a unique genealogy resource representing the homogeneous Utah population. Genome-wide linkage analysis of these pedigrees has identified a potential candidate predisposition gene, as well as multiple candidate regions that could harbor predisposition loci, and for which further analysis is suggested.

FANCM c5791C>T stopgain mutation (rs144567652) is a familial colorectal cancer risk factor.

PURPOSE: While familial aggregation of colorectal cancer (CRC) is recognized, the majority of the germline predisposition factors remain unidentified, and many high-risk CRC pedigrees remain unexplained by known risk variants. Fanconi Anemia genes have been recognized to be associated with cancer risk. Notably, FANCM (OMIM 609644) variants have been reported to confer risk for CRC and breast cancer. METHODS: Exome sequencing of CRC-affected cousins in a set of 47 independent extended high-risk CRC pedigrees identified a candidate set of rare, shared variants. Variants were tested for association with risk in 744 Utah CRC cases and 1525 controls, and for segregation with CRC in affected relatives. RESULTS: A FANCM stopgain variant was observed in two CRC-affected cousin pairs, each from an independent Utah high-risk pedigree, and yielded a nonsignificant, but elevated OR = 2.05 in a set of Utah cases and controls. Segregation of the variant to other related CRC-affected cases was observed in the two extended pedigrees. CONCLUSION: A rare stopgain variant in FANCM (rs144567652) that is recognized as a breast cancer predisposition variant, and that has previously been proposed, but not confirmed, as a CRC predisposition variant, is validated here as a risk factor for familial CRC.