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 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.

Analysis of high-risk pedigrees identifies 11 candidate variants for Alzheimer's disease.

INTRODUCTION: Analysis of sequence data in high-risk pedigrees is a powerful approach to detect rare predisposition variants. METHODS: Rare, shared candidate predisposition variants were identified from exome sequencing 19 Alzheimer's disease (AD)-affected cousin pairs selected from high-risk pedigrees. Variants were further prioritized by risk association in various external datasets. Candidate variants emerging from these analyses were tested for co-segregation to additional affected relatives of the original sequenced pedigree members. RESULTS: AD-affected high-risk cousin pairs contained 564 shared rare variants. Eleven variants spanning 10 genes were prioritized in external datasets: rs201665195 (ABCA7), and rs28933981 (TTR) were previously implicated in AD pathology; rs141402160 (NOTCH3) and rs140914494 (NOTCH3) were previously reported; rs200290640 (PIDD1) and rs199752248 (PIDD1) were present in more than one cousin pair; rs61729902 (SNAP91), rs140129800 (COX6A2, AC026471), and rs191804178 (MUC16) were not present in a longevity cohort; and rs148294193 (PELI3) and rs147599881 (FCHO1) approached significance from analysis of AD-related phenotypes. Three variants were validated via evidence of co-segregation to additional relatives (PELI3, ABCA7, and SNAP91). DISCUSSION: These analyses support ABCA7 and TTR as AD risk genes, expand on previously reported NOTCH3 variant identification, and prioritize seven additional candidate variants.

Colorectal cancer risk based on extended family history and body mass index.

Family history and body mass index (BMI) are well-known risk factors for colorectal cancer (CRC), however, their joint effects are not well described. Using linked data for genealogy, self-reported height and weight from driver's licenses, and the Utah Surveillance, Epidemiology, and End-Results cancer registry, we found that an increasing number of first-degree relatives (FDR) with CRC is associated with higher standardized incidence ratio (SIR) for overweight/obese probands but not for under/normal weight probands. For probands with two CRC-affected FDRs, the SIR = 1.91 (95% CI [0.52, 4.89]) for under/normal weight probands and SIR = 4.31 (95% CI [2.46, 7.00]) for overweight/obese probands. In the absence of CRC-affected FDRs, any number of CRC-affected SDRs did not significantly increase CRC risk for under/normal weight probands, but for overweight/obese probands with at least three CRC-affected SDRs the SIR = 2.68 (95% CI [1.29, 4.93]). In the absence of CRC-affected FDRs and SDRs, any number of CRC-affected third-degree relatives (TDRs) did not increase risk in under/normal weight probands, but significantly elevated risk for overweight/obese probands with at least two CRC-affected TDRs was observed; SIR = 1.32 (95% CI [1.04, 1.65]). For nonsyndromic CRC, maximum midlife BMI affects risk based on family history and should be taken into account for CRC risk communication when possible.

A role for the MEGF6 gene in predisposition to osteoporosis.

Osteoporosis is a common skeletal disorder characterized by deterioration of bone tissue. The set of genetic factors contributing to osteoporosis is not completely specified. High-risk osteoporosis pedigrees were analyzed to identify genes that may confer susceptibility to disease. Candidate predisposition variants were identified initially by whole exome sequencing of affected-relative pairs, approximately cousins, from 10 pedigrees. Variants were filtered on the basis of population frequency, concordance between pairs of cousins, affecting a gene associated with osteoporosis, and likelihood to have functionally damaging, pathogenic consequences. Subsequently, variants were tested for segregation in 68 additional relatives of the index carriers. A rare variant in MEGF6 (rs755467862) showed strong evidence of segregation with the disease phenotype. Predicted protein folding indicated the variant (Cys200Tyr) may disrupt structure of an EGF-like calcium-binding domain of MEGF6. Functional analyses demonstrated that complete loss of the paralogous genes megf6a and megf6b in zebrafish resulted in significant delay of cartilage and bone formation. Segregation analyses, in silico protein structure modeling, and functional assays support a role for MEGF6 in predisposition to osteoporosis.

Identification and genomic analysis of pedigrees with exceptional longevity identifies candidate rare variants.

BACKGROUND: Longevity as a phenotype entails living longer than average and typically includes living without chronic age-related diseases. Recently, several common genetic components to longevity have been identified. This study aims to identify additional genetic variants associated with longevity using unique and powerful analyses of pedigrees with a statistical excess of healthy elderly individuals identified in the Utah Population Database (UPDB). METHODS: From an existing biorepository of Utah pedigrees, six independent cousin pairs were selected from four extended pedigrees that exhibited an excess of healthy elderly individuals; whole exome sequencing (WES) was performed on two elderly individuals from each pedigree who were either first cousins or first cousins once removed. Rare (<.01 population frequency) variants shared by at least one elderly cousin pair in a region likely to be identical by descent were identified as candidates. Ingenuity Variant Analysis was used to prioritize putative causal variants based on quality control, frequency, and gain or loss of function. The variant frequency was compared in healthy cohorts and in an Alzheimer's disease cohort. Remaining variants were filtered based on their presence in genes reported to have an effect on the aging process, aging of cells, or the longevity process. Validation of these candidate variants included tests of segregation on other elderly relatives. RESULTS: Fifteen rare candidate genetic variants spanning 17 genes shared within cousins were identified as having passed prioritization criteria. Of those variants, six were present in genes that are known or predicted to affect the aging process: rs78408340 (PAM), rs112892337 (ZFAT), rs61737629 (ESPL1), rs141903485 (CEBPE), rs144369314 (UTP4), and rs61753103 (NUP88 and RABEP1). ESPL1 rs61737629 and CEBPE rs141903485 show additional evidence of segregation with longevity in expanded pedigree analyses (p-values = .001 and .0001, respectively). DISCUSSION: This unique pedigree analysis efficiently identified several novel rare candidate variants that may affect the aging process and added support to seven genes that likely contribute to longevity. Further analyses showed evidence for segregation for two rare variants, ESPL1 rs61737629 and CEBPE rs141903485, in the original longevity pedigrees in which they were initially observed. These candidate genes and variants warrant further investigation.

A population-based study of testicular cancer risk among children and young adults from Norway and Utah, USA.

Similar family-based cancer and genealogy data from Norway and Utah allowed comparisons of the incidence of testicular cancer (TC), and exploration of the role of Scandinavian ancestry and family history of TC in TC risk. Our study utilizes data from the Utah Population Database and Norwegian Population Registers. All males born during 1951-2015 were followed for TC until the age of 29 years. A total of 1,974,287 and 832,836 males were born in Norway and Utah, respectively, of whom 2,686 individuals were diagnosed with TC in Norway and 531 in Utah. The incidence per year of TC in Norway (10.6) was twice that observed in Utah (5.1) for males born in the last period (1980-1984). The incidence rates of TC in Utah did not differ according to the presence or absence of Scandinavian ancestry (p = 0.669). Having a brother diagnosed with TC was a strong risk factor for TC among children born in Norway and Utah, with HR = 9.87 (95% CI 5.68-17.16) and 6.02 (95% CI 4.80-7.55), respectively; with even higher HR observed among the subset of children in Utah with Scandinavian ancestry (HR = 12.30, 95% CI 6.78-22.31). A clear difference in TC incidence among individuals born in Norway and descendants of Scandinavian people born in Utah was observed. These differences in TC rates point to the possibility of environmental influence. Family history of TC is a strong risk factor for developing TC in both populations.

Risk of pelvic organ prolapse treatment based on extended family history.

BACKGROUND: Family history of pelvic organ prolapse among first-degree relatives is an established risk factor for pelvic organ prolapse; however, consideration of the constellation of family history that extends to distant relationships allows for more accurate determination of risk and may improve pelvic organ prolapse risk prediction estimates. OBJECTIVE: The purpose of this study was to assess risk for pelvic organ prolapse treatment based on varying family histories of pelvic organ prolapse and included number and types of affected relatives, ages of relatives at pelvic organ prolapse treatment, and whether the family history is of maternal or paternal origin. STUDY DESIGN: This was a retrospective, population-based study that involved the Utah Population Database, which is a population resource that includes extensive genealogy information linked to medical records. The study population included 453,522 total women: 4628 women with a diagnosis of treated (surgical or pessary) pelvic organ prolapse and their 15,530 first-degree relatives; 33,782 second-degree relatives, and 66,469 third-degree relatives. We estimated relative risk of treated pelvic organ prolapse based on specific family history constellations. RESULTS: Relative risk estimates increased with a family history of increasing numbers of treated first-degree relatives with pelvic organ prolapse (first-degree relatives, ≥1 [relative risk, 2.36; 95% confidence interval, 2.15-2.58], first-degree relatives, ≥2 [relative risk, 3.79; 95% confidence interval, 2.65-5.24], and first-degree relatives, ≥3 [relative risk, 6.26; 95% confidence interval, 1.29-18.30]). Having a family history of ≥3 affected third-degree relatives (eg, first cousins) and no affected first- or second-degree relatives was similar in risk to having 1 affected first-degree relative. Relative risk estimates decreased with increasing age of treatment for first-degree family members. Risks in individuals with a positive maternal family history for pelvic organ prolapse were consistently higher than risks in individuals with equivalent paternal family history, but paternal inheritance still played a role. Approximately 4% of the total studied female population was found to have a >2-fold risk of being treated for pelvic organ prolapse and is considered high-risk based on their family history. CONCLUSION: We provide estimates for treated pelvic organ prolapse based on an extensive family history of pelvic organ prolapse using a large population-based sample. Risk for treated pelvic organ prolapse increased with increasing numbers of affected close and distant female relatives, earlier age of pelvic organ prolapse treatment in relatives, and maternal inheritance. These risk estimates may be useful for genetic studies and investigation of risk reduction strategies in those at highest risk for pelvic organ prolapse.

Childhood central nervous system tumors and leukemia: Incidence and familial risk. A comparative population-based study in Utah and Norway.

BACKGROUND: In this study, we aimed to evaluate incidence rates and family risk of the most common childhood cancers, tumors in the central nervous system (CNS), and leukemia among individuals from Norway and individuals with Scandinavian ancestry living in Utah. METHODS: We used the Utah Population Database and the Norwegian National Population Register linked to Cancer registries to identify cancers in children born between 1966 and 2015 and their first-degree relatives. We calculated incidence rates and hazards ratios. RESULTS: The overall incidence of CNS tumors increased with consecutive birth cohorts similarly in Utah and Norway (both P < 0.001). Incidence rates of leukemia were more stable and similar in both Utah and in Norway with 4.6/100 000 person-years among children (<15 years) born in the last cohort. A family history of CNS tumors was significantly associated with risk of childhood CNS tumors in Utah HR = 3.05 (95% CI 1.80-5.16) and Norway HR = 2.87 (95% CI 2.20-3.74). In Norway, children with a first-degree relative diagnosed with leukemia had high risk of leukemia (HR = 2.39, 95% CI 1.61-3.55). CONCLUSION: Despite geographical distance and assumed large lifestyle differences, two genetically linked pediatric populations show similar incidences of CNS tumors and leukemia in the period 1966-2015. CNS tumors and leukemia aggregated in families in both countries.

Predictors of Response Outcomes for Research Recruitment Through a Central Cancer Registry: Evidence From 17 Recruitment Efforts for Population-Based Studies.

When recruiting research participants through central cancer registries, high response fractions help ensure population-based representation. We conducted multivariable mixed-effects logistic regression to identify case and study characteristics associated with making contact with and obtaining cooperation of Utah cancer cases using data from 17 unique recruitment efforts undertaken by the Utah Cancer Registry (2007-2016) on behalf of the following studies: A Population-Based Childhood Cancer Survivors Cohort Study in Utah, Comparative Effectiveness Analysis of Surgery and Radiation for Prostate Cancer (CEASAR Study), Costs and Benefits of Follow-up Care for Adolescent and Young Adult Cancers, Study of Exome Sequencing for Head and Neck Cancer Susceptibility Genes, Genetic Epidemiology of Chronic Lymphocytic Leukemia, Impact of Remote Familial Colorectal Cancer Risk Assessment and Counseling (Family CARE Project), Massively Parallel Sequencing for Familial Colon Cancer Genes, Medullary Thyroid Carcinoma (MTC) Surveillance Study, Osteosarcoma Surveillance Study, Prostate Cancer Outcomes Study, Risk Education and Assessment for Cancer Heredity Project (REACH Project), Study of Shared Genomic Segment Analysis and Tumor Subtyping in High-Risk Breast-Cancer Gene Pedigrees, Study of Shared Genomic Segment Analysis for Localizing Multiple Myeloma Genes. Characteristics associated with lower odds of contact included Hispanic ethnicity (odds ratio (OR) = 0.34, 95% confidence interval (CI): 0.27, 0.41), nonwhite race (OR = 0.46, 95% CI: 0.35, 0.60), and younger age at contact. Years since diagnosis was inversely associated with making contact. Nonwhite race and age ≥60 years had lower odds of cooperation. Study features with lower odds of cooperation included longitudinal design (OR = 0.50, 95% CI: 0.41, 0.61) and study brochures (OR = 0.70, 95% CI: 0.54, 0.90). Increased odds of cooperation were associated with including a questionnaire (OR = 3.19, 95% CI: 1.54, 6.59), postage stamps (OR = 1.60, 95% CI: 1.21, 2.12), and incentives (OR = 1.62, 95% CI: 1.02, 2.57). Among cases not responding after the first contact, odds of eventual response were lower when >10 days elapsed before subsequent contact (OR = 0.71, 95% CI: 0.59, 0.85). Obtaining high response is challenging, but study features identified in this analysis support better results when recruiting through central cancer registries.

Breast cancer histologic subtypes show excess familial clustering.

BACKGROUND: The inherited predisposition to developing specific histologic subtypes of invasive breast carcinoma has been incompletely investigated. By using a large, population-based database, the authors sought to investigate familial clustering of breast cancer by histologic subtype. METHODS: By using the Utah Population Database, which links genealogy records to the National Cancer Institute's statewide Surveillance, Epidemiology, and End Results cancer registry, the authors identified patients with breast cancer by histology and tested for evidence of shared genetic predisposition to histologic specific subtypes by examining pairwise relatedness and estimating the relative risk (RR) among first-degree, second-degree, and third-degree relatives. RESULTS: The authors identified 23,629 individuals in the Utah Population Database who had at least 3 generations of genealogy and at least 1 primary breast cancer, 2883 (12.2%) of which were specific histologic subtypes other than invasive ductal carcinoma (including inflammatory [n = 178], lobular [n = 1688], and mucinous [n = 542]). Statistically significant excess distant relatedness was identified for the mucinous subtype (P = .011) as well as for inflammatory breast cancers (P = .024). The RR for breast cancer of any histology in second-degree relatives was significantly increased for patients with inflammatory (RR, 1.32; 95% CI, 1.02-1.68; P = .03), lobular (RR, 1.36; 95% CI, 1.25-1.47; P < .001), and mucinous (RR, 1.27; 95% CI, 1.12-1.44; P = .00021) subtypes. CONCLUSIONS: These findings provide evidence for significant familial clustering within histological subtypes for lobular, mucinous, and inflammatory breast carcinomas. Further research is required to identify the underlying genetic variants responsible for the increased risk. Studies of high-risk pedigrees segregating a specific histologic subtype could be a powerful design for predisposition gene identification.

Family history of cancer and risk of paediatric and young adult's testicular cancer: A Norwegian cohort study.

BACKGROUND: The aim of this study was to examine the association of a family history of cancer with the risk of testicular cancer in young adults. METHODS: This is a prospective cohort study including 1,974,287 males born 1951-2015, of whom 2686 were diagnosed with TC before the age of 30. RESULTS: A history of TC in male relatives was significantly associated with a diagnosis of TC among children and young adults, including brothers (6.3-fold), sons (4.7-fold), fathers (4.4-fold), paternal uncles (2.0-fold) and maternal uncles (1.9-fold). Individuals with a father diagnosed with a carcinoma or sarcoma showed an elevated risk (1.1-fold and 1.8-fold, respectively). A family history of mesothelioma was positively associated with a risk of TC [(father (2.8-fold), mother (4.6-fold) and maternal uncles and aunt (4.4-fold)]. Elevated risks were also observed when siblings were diagnosed with malignant melanoma (1.4-fold). The risk of TC was also increased when fathers (11.1-fold), paternal (4.9-fold) and maternal uncles and aunts (4.6-fold) were diagnosed with malignant neuroepithelial-tumours. CONCLUSION: We found an increased risk of TC among children and young adults with a family history of TC, carcinoma, mesothelioma, sarcoma, malignant melanoma and malignant neuroepithelial tumours. Hereditary cancer syndromes might underlie some of the associations reported in this study.

REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants.

The vast majority of coding variants are rare, and assessment of the contribution of rare variants to complex traits is hampered by low statistical power and limited functional data. Improved methods for predicting the pathogenicity of rare coding variants are needed to facilitate the discovery of disease variants from exome sequencing studies. We developed REVEL (rare exome variant ensemble learner), an ensemble method for predicting the pathogenicity of missense variants on the basis of individual tools: MutPred, FATHMM, VEST, PolyPhen, SIFT, PROVEAN, MutationAssessor, MutationTaster, LRT, GERP, SiPhy, phyloP, and phastCons. REVEL was trained with recently discovered pathogenic and rare neutral missense variants, excluding those previously used to train its constituent tools. When applied to two independent test sets, REVEL had the best overall performance (p < 10-12) as compared to any individual tool and seven ensemble methods: MetaSVM, MetaLR, KGGSeq, Condel, CADD, DANN, and Eigen. Importantly, REVEL also had the best performance for distinguishing pathogenic from rare neutral variants with allele frequencies <0.5%. The area under the receiver operating characteristic curve (AUC) for REVEL was 0.046-0.182 higher in an independent test set of 935 recent SwissVar disease variants and 123,935 putatively neutral exome sequencing variants and 0.027-0.143 higher in an independent test set of 1,953 pathogenic and 2,406 benign variants recently reported in ClinVar than the AUCs for other ensemble methods. We provide pre-computed REVEL scores for all possible human missense variants to facilitate the identification of pathogenic variants in the sea of rare variants discovered as sequencing studies expand in scale.

Population-based relative risks for specific family history constellations of breast cancer.

PURPOSE: Using a large resource linking genealogy with decades of cancer data, a non-traditional approach was used to estimate individualized risk for breast cancer (BC) based on specific family history extending to first cousins, providing a clearer picture of the contribution of various aspects of both close and distant combinations of affected relatives. METHODS: RRs for BC were estimated in 640,366 females for a representative set of breast cancer family history constellations that included number of first- (FDR), second-(SDR), and third-degree relatives (TDR), maternal and paternal relatives, and age at earliest diagnosis in a relative. RESULTS: RRs for first-degree relatives of BC cases ranged from 1.61 (= 1 FDR affected, CI 1.56, 1.67) to 5.00 (≥ 4 FDRs affected, CI 3.35, 7.18). RRs for second-degree relatives of probands with 0 affected FDRs ranged from 1.04 (= 1 SDR affected, CI 1.00, 1.08) to 1.71 (≥ 4 SDRs affected, CI 1.26, 2.27) and for second-degree relatives of probands with exactly 1 FDR from 1.54 (0 SDRs affected, CI 1.47, 1.61) to 4.78 (≥ 5 SDRs; CI 2.47, 8.35). RRs for third-degree relatives with no closer relatives affected were significantly elevated over population risk for probands with ≥ 5 affected TDRs RR = 1.32, CI 1.11, 1.57). CONCLUSIONS: The majority of females in the Utah resource had a positive family history of BC in FDRs to TDRs. Presence of any number of affected FDRs or SDRs significantly increased risk for BC over population risk; and more than four TDRs, even with no affected FDRs or SDRs, significantly increased risk over population risk. Risk prediction derived from the specific and extended family history constellation of affected relatives allows identification of females at increased risk even when they do not have a conventionally defined high-risk family; these risks could be a powerful, efficient tool to individualize cancer screening and prevention.

AA9int: SNP interaction pattern search using non-hierarchical additive model set.

Motivation: The use of single nucleotide polymorphism (SNP) interactions to predict complex diseases is getting more attention during the past decade, but related statistical methods are still immature. We previously proposed the SNP Interaction Pattern Identifier (SIPI) approach to evaluate 45 SNP interaction patterns/patterns. SIPI is statistically powerful but suffers from a large computation burden. For large-scale studies, it is necessary to use a powerful and computation-efficient method. The objective of this study is to develop an evidence-based mini-version of SIPI as the screening tool or solitary use and to evaluate the impact of inheritance mode and model structure on detecting SNP-SNP interactions. Results: We tested two candidate approaches: the 'Five-Full' and 'AA9int' method. The Five-Full approach is composed of the five full interaction models considering three inheritance modes (additive, dominant and recessive). The AA9int approach is composed of nine interaction models by considering non-hierarchical model structure and the additive mode. Our simulation results show that AA9int has similar statistical power compared to SIPI and is superior to the Five-Full approach, and the impact of the non-hierarchical model structure is greater than that of the inheritance mode in detecting SNP-SNP interactions. In summary, it is recommended that AA9int is a powerful tool to be used either alone or as the screening stage of a two-stage approach (AA9int+SIPI) for detecting SNP-SNP interactions in large-scale studies. Availability and implementation: The 'AA9int' and 'parAA9int' functions (standard and parallel computing version) are added in the SIPI R package, which is freely available at https://linhuiyi.github.io/LinHY_Software/. Supplementary information: Supplementary data are available at Bioinformatics online.

The contribution of the rs55705857 G allele to familial cancer risk as estimated in the Utah population database.

BACKGROUND: IDH1/2 mutated glioma has been associated with a germline risk variant, the rs55705857 G allele. The Utah Population Database (UPDB), a computerized genealogy of people in Utah, is a unique resource to evaluate cancer risk in related individuals. METHODS: One hundred and two individuals with IDH1/2 mutant or 1p/19q co-deleted glioma were genotyped and linked to the UPDB. DNA came from blood (21), tumor tissue (43), or both (38). We determined congruence between somatic and germline samples and estimated the relative risk for developing cancer to first and second-degree relatives of G and A allele carriers at rs55705857. RESULTS: Somatic (glioma) DNA had 85.7% sensitivity (CI 57.2-98.2%) and 95.8% specificity (CI 78.9-99.89%) for germline rs55705857 G allele. Forty-one patients were linked to pedigrees in the UPDB with at least three generations of data. First-degree relatives of rs55705857 G allele carriers were at significantly increased risk for developing cancer (RR = 1.72, p = 0.045, CI 1.02-2.94), and specifically for oligodendroglioma (RR = 57.61, p = 0.017, CI 2.96-320.98) or prostate cancer (RR = 4.10, p = 0.008, CI 1.62-9.58); relatives of individuals without the G allele were not at increased risk. Second-degree relatives of G allele carriers also had significantly increased risk for developing cancer (RR = 1.50, p = 0.007, CI 1.15-2.01). CONCLUSIONS: Tumor DNA may approximate genotype at the rs55705857 locus. We confirmed this locus confers an increased risk of all cancers and especially of oligodendroglioma. No increased cancer or brain tumor risk is seen in family members of individuals without the high-risk G allele.

gsSKAT: Rapid gene set analysis and multiple testing correction for rare-variant association studies using weighted linear kernels.

Next-generation sequencing technologies have afforded unprecedented characterization of low-frequency and rare genetic variation. Due to low power for single-variant testing, aggregative methods are commonly used to combine observed rare variation within a single gene. Causal variation may also aggregate across multiple genes within relevant biomolecular pathways. Kernel-machine regression and adaptive testing methods for aggregative rare-variant association testing have been demonstrated to be powerful approaches for pathway-level analysis, although these methods tend to be computationally intensive at high-variant dimensionality and require access to complete data. An additional analytical issue in scans of large pathway definition sets is multiple testing correction. Gene set definitions may exhibit substantial genic overlap, and the impact of the resultant correlation in test statistics on Type I error rate control for large agnostic gene set scans has not been fully explored. Herein, we first outline a statistical strategy for aggregative rare-variant analysis using component gene-level linear kernel score test summary statistics as well as derive simple estimators of the effective number of tests for family-wise error rate control. We then conduct extensive simulation studies to characterize the behavior of our approach relative to direct application of kernel and adaptive methods under a variety of conditions. We also apply our method to two case-control studies, respectively, evaluating rare variation in hereditary prostate cancer and schizophrenia. Finally, we provide open-source R code for public use to facilitate easy application of our methods to existing rare-variant analysis results.

Prostate cancer risk regions at 8q24 and 17q24 are differentially associated with somatic TMPRSS2:ERG fusion status.

Molecular and epidemiological differences have been described between TMPRSS2:ERG fusion-positive and fusion-negative prostate cancer (PrCa). Assuming two molecularly distinct subtypes, we have examined 27 common PrCa risk variants, previously identified in genome-wide association studies, for subtype specific associations in a total of 1221 TMPRSS2:ERG phenotyped PrCa cases. In meta-analyses of a discovery set of 552 cases with TMPRSS2:ERG data and 7650 unaffected men from five centers we have found support for the hypothesis that several common risk variants are associated with one particular subtype rather than with PrCa in general. Risk variants were analyzed in case-case comparisons (296 TMPRSS2:ERG fusion-positive versus 256 fusion-negative cases) and an independent set of 669 cases with TMPRSS2:ERG data was established to replicate the top five candidates. Significant differences (P < 0.00185) between the two subtypes were observed for rs16901979 (8q24) and rs1859962 (17q24), which were enriched in TMPRSS2:ERG fusion-negative (OR = 0.53, P = 0.0007) and TMPRSS2:ERG fusion-positive PrCa (OR = 1.30, P = 0.0016), respectively. Expression quantitative trait locus analysis was performed to investigate mechanistic links between risk variants, fusion status and target gene mRNA levels. For rs1859962 at 17q24, genotype dependent expression was observed for the candidate target gene SOX9 in TMPRSS2:ERG fusion-positive PrCa, which was not evident in TMPRSS2:ERG negative tumors. The present study established evidence for the first two common PrCa risk variants differentially associated with TMPRSS2:ERG fusion status. TMPRSS2:ERG phenotyping of larger studies is required to determine comprehensive sets of variants with subtype-specific roles in PrCa.

Family history of cancer and the risk of childhood solid tumours: a Norwegian nationwide register-based cohort study.

This corrects the article DOI: 10.1038/bjc.2017.85.

Gene and pathway level analyses of germline DNA-repair gene variants and prostate cancer susceptibility using the iCOGS-genotyping array.

This corrects the article DOI: 10.1038/bjc.2016.50.