Joseph H. Rothstein

Genome-Wide Association of Familial Late-Onset Alzheimer's Disease Replicates BIN1 and CLU and Nominates CUGBP2 in Interaction with APOE

Late-onset Alzheimers disease (LOAD) is the most common form of dementia in the elderly. The National Institute of Aging-Late Onset Alzheimers Disease Family Study and the National Cell Repository for Alzheimers Disease conducted a joint genome-wide association study (GWAS) of multiplex LOAD families (3,839 affected and unaffected individuals from 992 families plus additional unrelated neurologically evaluated normal subjects) using the 610 IlluminaQuad panel. This cohort represents the largest family-based GWAS of LOAD to date, with analyses limited here to the European-American subjects. SNPs near APOE gave highly significant results (e.g., rs2075650, p = 3.2×10−81), but no other genome-wide significant evidence for association was obtained in the full sample. Analyses that stratified on APOE genotypes identified SNPs on chromosome 10p14 in CUGBP2 with genome-wide significant evidence for association within APOE ε4 homozygotes (e.g., rs201119, p = 1.5×10−8). Association in this gene was replicated in an independent sample consisting of three cohorts. There was evidence of association for recently-reported LOAD risk loci, including BIN1 (rs7561528, p = 0.009 with, and p = 0.03 without, APOE adjustment) and CLU (rs11136000, p = 0.023 with, and p = 0.008 without, APOE adjustment), with weaker support for CR1. However, our results provide strong evidence that association with PICALM (rs3851179, p = 0.69 with, and p = 0.039 without, APOE adjustment) and EXOC3L2 is affected by correlation with APOE, and thus may represent spurious association. Our results indicate that genetic structure coupled with ascertainment bias resulting from the strong APOE association affect genome-wide results and interpretation of some recently reported associations. We show that a locus such as APOE, with large effects and strong association with disease, can lead to samples that require appropriate adjustment for this locus to avoid both false positive and false negative evidence of association. We suggest that similar adjustments may also be needed for many other large multi-site studies.

Hormone-receptor expression and ovarian cancer survival: an Ovarian Tumor Tissue Analysis consortium study

BACKGROUND Few biomarkers of ovarian cancer prognosis have been established, partly because subtype-specific associations might be obscured in studies combining all histopathological subtypes. We examined whether tumour expression of the progesterone receptor (PR) and oestrogen receptor (ER) was associated with subtype-specific survival. METHODS 12 studies participating in the Ovarian Tumor Tissue Analysis consortium contributed tissue microarray sections and clinical data to our study. Participants included in our analysis had been diagnosed with invasive serous, mucinous, endometrioid, or clear-cell carcinomas of the ovary. For a patient to be eligible, tissue microarrays, clinical follow-up data, age at diagnosis, and tumour grade and stage had to be available. Clinical data were obtained from medical records, cancer registries, death certificates, pathology reports, and review of histological slides. PR and ER statuses were assessed by central immunohistochemistry analysis done by masked pathologists. PR and ER staining was defined as negative (<1% tumour cell nuclei), weak (1 to <50%), or strong (≥50%). Associations with disease-specific survival were assessed. FINDINGS 2933 women with invasive epithelial ovarian cancer were included: 1742 with high-grade serous carcinoma, 110 with low-grade serous carcinoma, 207 with mucinous carcinoma, 484 with endometrioid carcinoma, and 390 with clear-cell carcinoma. PR expression was associated with improved disease-specific survival in endometrioid carcinoma (log-rank p<0·0001) and high-grade serous carcinoma (log-rank p=0·0006), and ER expression was associated with improved disease-specific survival in endometrioid carcinoma (log-rank p<0·0001). We recorded no significant associations for mucinous, clear-cell, or low-grade serous carcinoma. Positive hormone-receptor expression (weak or strong staining for PR or ER, or both) was associated with significantly improved disease-specific survival in endometrioid carcinoma compared with negative hormone-receptor expression, independent of study site, age, stage, and grade (hazard ratio 0·33, 95% CI 0·21-0·51; p<0·0001). Strong PR expression was independently associated with improved disease-specific survival in high-grade serous carcinoma (0·71, 0·55-0·91; p=0·0080), but weak PR expression was not (1·02, 0·89-1·18; p=0·74). INTERPRETATION PR and ER are prognostic biomarkers for endometrioid and high-grade serous ovarian cancers. Clinical trials, stratified by subtype and biomarker status, are needed to establish whether hormone-receptor status predicts response to endocrine treatment, and whether it could guide personalised treatment for ovarian cancer. FUNDING Carraresi Foundation and others.

Evidence for multiple loci from a genome scan of autism kindreds

We performed a genome-wide linkage scan using highly polymorphic microsatellite markers. To minimize genetic heterogeneity, we focused on sibpairs meeting the strict diagnosis of autism. In our primary analyses, we observed a strong linkage signal (P=0.0006, 133.16 cM) on chromosome 7q at a location coincident with other linkage studies. When a more relaxed diagnostic criteria was used, linkage evidence at this location was weaker (P=0.01). The sample was stratified into families with only male affected subjects (MO) and families with at least one female affected subject (FC). The strongest signal unique to the MO group was on chromosome 11 (P=0.0009, 83.82 cM), and for the FC group on chromosome 4 (P=0.002, 111.41 cM). We also divided the sample into regression positive and regression negative families. The regression-positive group showed modest linkage signals on chromosomes 10 (P=0.003, 0 cM) and 14 (P=0.005, 104.2 cM). More significant peaks were seen in the regression negative group on chromosomes 3 (P=0.0002, 140.06 cM) and 4 (P=0.0005, 111.41 cM). Finally, we used language acquisition data as a quantitative trait in our linkage analysis and observed a chromosome 9 signal (149.01 cM) of P=0.00006 and an empirical P-value of 0.0008 at the same location. Our work provides strong conformation for an autism locus on 7q and suggestive evidence for several other chromosomal locations. Diagnostic specificity and detailed analysis of the autism phenotype is critical for identifying autism loci.

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.

Multipoint linkage analysis with many multiallelic or dense diallelic markers: Markov chain-Monte Carlo provides practical approaches for genome scans on general pedigrees.

Computations for genome scans need to adapt to the increasing use of dense diallelic markers as well as of full-chromosome multipoint linkage analysis with either diallelic or multiallelic markers. Whereas suitable exact-computation tools are available for use with small pedigrees, equivalent exact computation for larger pedigrees remains infeasible. Markov chain-Monte Carlo (MCMC)-based methods currently provide the only computationally practical option. To date, no systematic comparison of the performance of MCMC-based programs is available, nor have these programs been systematically evaluated for use with dense diallelic markers. Using simulated data, we evaluate the performance of two MCMC-based linkage-analysis programs--lm_markers from the MORGAN package and SimWalk2--under a variety of analysis conditions. Pedigrees consisted of 14, 52, or 98 individuals in 3, 5, or 6 generations, respectively, with increasing amounts of missing data in larger pedigrees. One hundred replicates of markers and trait data were simulated on a 100-cM chromosome, with up to 10 multiallelic and up to 200 diallelic markers used simultaneously for computation of multipoint LOD scores. Exact computation was available for comparison in most situations, and comparison with a perfectly informative marker or interprogram comparison was available in the remaining situations. Our results confirm the accuracy of both programs in multipoint analysis with multiallelic markers on pedigrees of varied sizes and missing-data patterns, but there are some computational differences. In contrast, for large numbers of dense diallelic markers, only the lm_markers program was able to provide accurate results within a computationally practical time. Thus, programs in the MORGAN package are the first available to provide a computationally practical option for accurate linkage analyses in genome scans with both large numbers of diallelic markers and large pedigrees.

Genomewide scan for real-word reading subphenotypes of dyslexia: Novel chromosome 13 locus and genetic complexity

Dyslexia is a common learning disability exhibited as a delay in acquiring reading skills despite adequate intelligence and instruction. Reading single real words (real‐word reading, RWR) is especially impaired in many dyslexics. We performed a genome scan, using variance components (VC) linkage analysis and Bayesian Markov chain Monte Carlo (MCMC) joint segregation and linkage analysis, for three quantitative measures of RWR in 108 multigenerational families, with follow up of the strongest signals with parametric LOD score analyses. We used single‐word reading efficiency (SWE) to assess speed and accuracy of RWR, and word identification (WID) to assess accuracy alone. Adjusting SWE for WID provided a third measure of RWR efficiency. All three methods of analysis identified a strong linkage signal for SWE on chromosome 13q. Based on multipoint analysis with 13 markers we obtained a MCMC intensity ratio (IR) of 53.2 (chromosome‐wide P < 0.004), a VC LOD score of 2.29, and a parametric LOD score of 2.94, based on a quantitative‐trait model from MCMC segregation analysis (SA). A weaker signal for SWE on chromosome 2q occurred in the same location as a significant linkage peak seen previously in a scan for phonological decoding. MCMC oligogenic SA identified three models of transmission for WID, which could be assigned to two distinct linkage peaks on chromosomes 12 and 15. Taken together, these results indicate a locus for efficiency and accuracy of RWR on chromosome 13, and a complex model for inheritance of RWR accuracy with loci on chromosomes 12 and 15.

Identification of novel susceptibility loci for Guam neurodegenerative disease: Challenges of genome scans in genetic isolates

Amyotrophic lateral sclerosis/parkinsonism-dementia complex (ALS/PDC) is a fatal neurodegenerative disease found in the Chamorro people of Guam and other Pacific Island populations. The etiology is unknown, although both genetic and environmental factors appear important. To identify loci for ALS/PDC, we conducted both genome-wide linkage and association analyses, using approximately 400 microsatellite markers, in the largest sample assembled to date, comprising a nearly complete sample of all living and previously sampled deceased cases. A single, large, complex pedigree was ascertained from a village on Guam, with smaller families and a case-control sample ascertained from the rest of Guam by population-based neurological screening and archival review. We found significant evidence for two regions with novel ALS/PDC loci on chromosome 12 and supportive evidence for the involvement of the MAPT region on chromosome 17. D12S1617 on 12p gave the strongest evidence of linkage (maximum LOD score, Z(max) = 4.03) in our initial scan, with additional support in the complete case-control sample in the form of evidence of allelic association at this marker and another nearby marker. D12S79 on 12q also provided significant evidence of linkage (Z(max) = 3.14) with support from flanking markers. Our results suggest that ALS/PDC may be influenced by as many as three loci, while illustrating challenges that are intrinsic in genetic analyses of isolated populations, as well as analytical strategies that are useful in this context. Elucidation of the genetic basis of ALS/PDC should improve our understanding of related neurodegenerative disorders including Alzheimer disease, Parkinson disease, frontotemporal dementia and ALS.

Breast Cancer Risk and Mammographic Density Assessed with Semiautomated and Fully Automated Methods and BI-RADS

Purpose To compare three metrics of breast density on full-field digital mammographic (FFDM) images as predictors of future breast cancer risk. Materials and Methods This institutional review board-approved study included 125 women with invasive breast cancer and 274 age- and race-matched control subjects who underwent screening FFDM during 2004-2013 and provided informed consent. The percentage of density and dense area were assessed semiautomatically with software (Cumulus 4.0; University of Toronto, Toronto, Canada), and volumetric percentage of density and dense volume were assessed automatically with software (Volpara; Volpara Solutions, Wellington, New Zealand). Clinical Breast Imaging Reporting and Data System (BI-RADS) classifications of breast density were extracted from mammography reports. Odds ratios and 95% confidence intervals (CIs) were estimated by using conditional logistic regression stratified according to age and race and adjusted for body mass index, parity, and menopausal status, and the area under the receiver operating characteristic curve (AUC) was computed. Results The adjusted odds ratios and 95% CIs for each standard deviation increment of the percentage of density, dense area, volumetric percentage of density, and dense volume were 1.61 (95% CI: 1.19, 2.19), 1.49 (95% CI: 1.15, 1.92), 1.54 (95% CI: 1.12, 2.10), and 1.41 (95% CI: 1.11, 1.80), respectively. Odds ratios for women with extremely dense breasts compared with those with scattered areas of fibroglandular density were 2.06 (95% CI: 0.85, 4.97) and 2.05 (95% CI: 0.90, 4.64) for BI-RADS and Volpara density classifications, respectively. Clinical BI-RADS was more accurate (AUC, 0.68; 95% CI: 0.63, 0.74) than Volpara (AUC, 0.64; 95% CI: 0.58, 0.70) and continuous measures of percentage of density (AUC, 0.66; 95% CI: 0.60, 0.72), dense area (AUC, 0.66; 95% CI: 0.60, 0.72), volumetric percentage of density (AUC, 0.64; 95% CI: 0.58, 0.70), and density volume (AUC, 0.65; 95% CI: 0.59, 0.71), although the AUC differences were not statistically significant. Conclusion Mammographic density on FFDM images was positively associated with breast cancer risk by using the computer assisted methods and BI-RADS. BI-RADS classification was as accurate as computer-assisted methods for discrimination of patients from control subjects.

Comparison of marker types and map assumptions using Markov chain Monte Carlo-based linkage analysis of COGA data

We performed multipoint linkage analysis of the electrophysiological trait ECB21 on chromosome 4 in the full pedigrees provided by the Collaborative Study on the Genetics of Alcoholism (COGA). Three Markov chain Monte Carlo (MCMC)-based approaches were applied to the provided and re-estimated genetic maps and to five different marker panels consisting of microsatellite (STRP) and/or SNP markers at various densities. We found evidence of linkage near the GABRB1 STRP using all methods, maps, and marker panels. Difficulties encountered with SNP panels included convergence problems and demanding computations.

Dose-Response Association of CD8+ Tumor-Infiltrating Lymphocytes and Survival Time in High-Grade Serous Ovarian Cancer.

Importance Cytotoxic CD8+ tumor-infiltrating lymphocytes (TILs) participate in immune control of epithelial ovarian cancer; however, little is known about prognostic patterns of CD8+ TILs by histotype and in relation to other clinical factors. Objective To define the prognostic role of CD8+ TILs in epithelial ovarian cancer. Design, Setting, and Participants This was a multicenter observational, prospective survival cohort study of the Ovarian Tumor Tissue Analysis Consortium. More than 5500 patients, including 3196 with high-grade serous ovarian carcinomas (HGSOCs), were followed prospectively for over 24 650 person-years. Exposures Following immunohistochemical analysis, CD8+ TILs were identified within the epithelial components of tumor islets. Patients were grouped based on the estimated number of CD8+ TILs per high-powered field: negative (none), low (1-2), moderate (3-19), and high (≥20). CD8+ TILs in a subset of patients were also assessed in a quantitative, uncategorized manner, and the functional form of associations with survival was assessed using penalized B-splines. Main Outcomes and Measures Overall survival time. Results The final sample included 5577 women; mean age at diagnosis was 58.4 years (median, 58.2 years). Among the 5 major invasive histotypes, HGSOCs showed the most infiltration. CD8+ TILs in HGSOCs were significantly associated with longer overall survival; median survival was 2.8 years for patients with no CD8+ TILs and 3.0 years, 3.8 years, and 5.1 years for patients with low, moderate, or high levels of CD8+ TILs, respectively (P value for trend = 4.2 × 10−16). A survival benefit was also observed among women with endometrioid and mucinous carcinomas, but not for those with the other histotypes. Among HGSOCs, CD8+ TILs were favorable regardless of extent of residual disease following cytoreduction, known standard treatment, and germline BRCA1 pathogenic mutation, but were not prognostic for BRCA2 mutation carriers. Evaluation of uncategorized CD8+ TIL counts showed a near-log-linear functional form. Conclusions and Relevance This study demonstrates the histotype-specific nature of immune infiltration and provides definitive evidence for a dose-response relationship between CD8+ TILs and HGSOC survival. That the extent of infiltration is prognostic, not merely its presence or absence, suggests that understanding factors that drive infiltration will be the key to unraveling outcome heterogeneity in this cancer.