Citations

Mixture analysis is an emerging statistical tool in epidemiological research that seeks to estimate the health effects associated with mixtures of several exposures. This approach acknowledges that individuals experience many simultaneous exposures and it can estimate the relative importance of components in the mixture. Health effects due to mixtures may vary over space driven by to political, demographic, environmental, or other differences. In such cases, estimating a global mixture effect without accounting for spatial variation would induce bias in effect estimates and potentially lower statistical power. To date, no methods have been developed to estimate spatially varying chemical mixture effects. We developed a Bayesian spatially varying mixture model that estimates spatially varying mixture effects and the importance weights of components in the mixture, while adjusting for covariates. We demonstrate the efficacy of the model through a simulation study that varies the number of mixtures (one and two) and spatial pattern (global, one‐dimensional, radial) and magnitude of mixture effects, showing that the model is able to accurately reproduce the spatial pattern of mixture effects across a diverse set of scenarios. Finally, we apply our model to a multi‐center case‐control study of non‐Hodgkin lymphoma (NHL) in Detroit, Iowa, Los Angeles, and Seattle. We identify significant spatially varying positive and inverse associations with NHL for two mixtures of pesticides in Iowa and do not find strong spatial effects at the other three centers. In conclusion, the Bayesian spatially varying mixture model represents a novel method for modeling spatial variation in mixture effects.

The exposome is an ideal in public health research that posits that individuals experience risk for adverse health outcomes from a wide variety of sources over their lifecourse. There have been increases in data collection in the various components of the exposome, but novel statistical methods are needed that capture multiple dimensions of risk at once. We introduce a Bayesian index low‐rank kriging (LRK) multiple membership model (MMM) to simultaneously estimate the health effects of one or more groups of exposures, the relative importance of exposure components, and cumulative spatial risk over time using residential histories. The model employs an MMM to consider all residential locations for subjects weighted by duration and LRK to increase computational efficiency. We demonstrate the performance of the Bayesian index LRK‐MMM through a simulation study, showing that the model accurately and consistently estimates the health effects of one or several group indices and has high power to identify a region of elevated spatial risk due to unmeasured environmental exposures. Finally, we apply our model to data from a multicenter case‐control study of non‐Hodgkin lymphoma (NHL), finding a significant positive association between one index of pesticides and risk for NHL in Iowa. Additionally, we find an area of significantly elevated spatial risk for NHL in Los Angeles. In conclusion, our Bayesian index LRK‐MMM represents a step forward toward bringing the ideals of the exposome into practice for environmental risk analyzes.

Abstract Background and Aims Although many studies revealed transcriptomic subtypes of HCC, concordance of the subtypes are not fully examined. We aim to examine a consensus of transcriptomic subtypes and correlate them with clinical outcomes. Approach and Results By integrating 16 previously established genomic signatures for HCC subtypes, we identified five clinically and molecularly distinct consensus subtypes. STM (STeM) is characterized by high stem cell features, vascular invasion, and poor prognosis. CIN (Chromosomal INstability) has moderate stem cell features, but high genomic instability and low immune activity. IMH (IMmune High) is characterized by high immune activity. BCM (Beta‐Catenin with high Male predominance) is characterized by prominent β‐catenin activation, low miRNA expression, hypomethylation, and high sensitivity to sorafenib. DLP (Differentiated and Low Proliferation) is differentiated with high hepatocyte nuclear factor 4A activity. We also developed and validated a robust predictor of consensus subtype with 100 genes and demonstrated that five subtypes were well conserved in patient‐derived xenograft models and cell lines. By analyzing serum proteomic data from the same patients, we further identified potential serum biomarkers that can stratify patients into subtypes. Conclusions Five HCC subtypes are correlated with genomic phenotypes and clinical outcomes and highly conserved in preclinical models, providing a framework for selecting the most appropriate models for preclinical studies.

Background In participants with concomitant chronic coronary disease and advanced chronic kidney disease (CKD), the effect of treatment strategies on the timing of dialysis initiation is not well characterized. Methods and Results In ISCHEMIA‐CKD (International Study of Comparative Health Effectiveness With Medical and Invasive Approaches–Chronic Kidney Disease), 777 participants with advanced CKD and moderate or severe ischemia were randomized to either an initial invasive or conservative management strategy. Herein, we compare the proportion of randomized participants with non–dialysis‐requiring CKD at baseline (n=362) who initiated dialysis and compare the time to dialysis initiation between invasive versus conservative management arms. Using multivariable Cox regression analysis, we also sought to identify the effect of invasive versus conservative chronic coronary disease management strategies on dialysis initiation. At a median follow‐up of 23 months (25th–75th interquartile range, 14–32 months), dialysis was initiated in 18.9% of participants (36/190) in the invasive strategy and 16.9% of participants (29/172) in the conservative strategy (P= 0.22). The median time to dialysis initiation was 6.0 months (interquartile range, 3.0–16.0 months) in the invasive group and 18.2 months (interquartile range, 12.2–25.0 months) in the conservative group (P= 0.004), with no difference in procedural acute kidney injury rates between the groups (7.8% versus 5.4%; P =0.26). Baseline clinical factors associated with earlier dialysis initiation were lower baseline estimated glomerular filtration rate (hazard ratio [HR] associated with 5‐unit decrease, 2.08 [95% CI, 1.72–2.56]; P <0.001), diabetes (HR, 2.30 [95% CI, 1.28–4.13]; P =0.005), hypertension (HR, 7.97 [95% CI, 1.09–58.21]; P =0.041), and Hispanic ethnicity (HR, 2.34 [95% CI, 1.22–4.47]; P =0.010). Conclusions In participants with non–dialysis‐requiring CKD in ISCHEMIA‐CKD, randomization to an invasive chronic coronary disease management strategy (relative to a conservative chronic coronary disease management strategy) is associated with an accelerated time to initiation of maintenance dialysis for kidney failure. Registration URL: https://www.clinicaltrials.gov ; Unique identifier: NCT01985360.

Abstract The Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) in collaboration with the Colorectal Cancer Family Registry (CCFR) aims to identify genetic variants and environmental risk factors that impact colorectal cancer (CRC). Over 30 studies from North America, Europe, and Australia participate in GECCO. These studies have collected clinical, epidemiological, and survival data, as well as blood and tumor biospecimens, for over 80,000 CRC cases and controls. The current study aims to conduct targeted deep sequencing of tumors and matching normal DNA to identify recurrent and novel somatic and germline variants in 4,200 CRC cases. To achieve this goal, an AmpliSeq targeted sequencing panel of 1.12 Mbp was constructed to sequence the coding regions of 190 significantly mutated genes identified from whole exome sequencing datasets generated by the Nurses’ Health Study and Health Professional’s Follow-up Study, and The Cancer Genome Atlas. The panel also covers coding regions of 15 genes with germline high penetrance mutations in CRC, 54 regions to detect CRC-related copy number alterations (CNAs), and microsatellite and homopolymer repeat regions to identify defective DNA mismatch repair. Primers were also included to detect Fusobacterium nucleatum DNA in tumor biopsies, as F. nucleatum is thought to promote CRC carcinogenesis. Sequencing of the DNA libraries on Illumina HiSeq 2500 produced a mean coverage of greater than 500X for tumor DNA and 100X for normal DNA, with &gt;85% of the bases covered at the target at 50x. So far, targeted sequencing of &gt;1,500 DNA samples from CRC tumors and normal tissues has identified recurrent and novel somatic mutations, germline genetic variants, and hypermutation status of the tumors due to defective DNA mismatch repair or pathogenic mutations in the POLE gene. Targeted sequencing has also allowed quantification of the F. nucleatum DNA in tumor biopsies; the results were validated by a multiplex QPCR assay. At the AACR annual meeting, we will present targeted sequencing results generated from the first two GECCO-participating studies (n=1,300 cases). These data will be valuable for future association testing of somatic mutations, CNAs, hypermutation status, and F. nucleatum with germline genetic variants, lifestyle, and environmental risk factors and survival. This large study will allow development of better strategies for diagnosis, treatment, and prevention of CRC. Citation Format: Syed H. Zaidi, Wei Sun, Jeroen Huyghe, Catherine S. Grasso, Quang Trinh, Charles Connolly, Amy French, Jasmine Mu, Marios Giannakis, Eve Shinbrot, Ivan Borozan, Michael J. Quist, Hermann Brenner, Daniel Buchanan, Peter Campbell, Andrew Chan, Jenny Chang-Claude, Vincent Ferretti, Charles Fuchs, Andrea Gsur, Marc Gunter, Tabitha Harrison, Michael Hoffmeister, Wen-Yi Huang, Paul Krzyzanowski, Stephen Lee, Mathieu Lemire, Jessica Miller, Danielle Pasternack, Cherie Teney, Elaine Mardis, Polly Newcomb, Lincoln Stein, Lee Timms, David Wheeler, Christina Yung, Niha Zubair, Levi Garraway, Shuji Ogino, Li Hsu, Steven Gallinger, Stephen Thibodeau, Thomas Hudson, Ulrike Peters. Targeted deep sequencing of colorectal tumor tissues to study associations of tumor subtypes with germline genetic, lifestyle, and environmental risk factors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1286. doi:10.1158/1538-7445.AM2017-1286

Abstract The completion of The Cancer Genome Atlas (TCGA) project for colorectal cancer (CRC) has enabled a new, focused phase of sequencing tumor samples for which genome-wide genetic, epidemiological, clinical and lifestyle data have been collected. By combining somatic mutational profiles with these aforementioned data, we may identify and report effective prevention and treatment approaches for a broader population of individuals. The advent of targeted deep sequencing approaches using DNA obtained from formalin fixed paraffin embedded tissues has made possible the genetic characterization of the large numbers of patients needed to make such an effort relevant. As a first step, we describe a custom gene panel designed from large-scale studies for targeted deep sequencing, and its application to over 4,200 CR tumors, collected by the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). This study is designed to identify recurrent somatic mutations and copy number alterations (CNAs) for association testing with germline genetic and lifestyle and environmental risk factors for CRC, and thereby identify relevant approaches to impact cancer prevention. The targeted CRC panel includes 205 genes. Genes were primarily selected as significantly mutated genes identified from the Nurses’ Health Study and Health Professional's Follow-up Study (n = 700), and TCGA (n = 525). We also included 15 genes associated with high penetrance germline mutations and augmented the list to include genes in pathways of somatically altered genes, identified by literature review, from public databases and known to be associated with loss of heterozygosity. For these 205 genes, amplification primers were designed to include all coding regions of transcripts that are listed in the UCSC Genome Browser. For regions with CNAs, the TCGA dataset was analyzed to include regions with greater than or equal to 2 copy focal gains or losses that were found in more than 4 or 3 tumors, respectively. Candidate genes in regions with significant CNAs from the TCGA CRC publication (Nature 2012) were also included. For CNAs, 6 to 12 amplicons were designed for each of the 56 selected regions (32 gains and 24 losses). Additional target regions include: 1) 25 microsatellite loci previously used to identify defective DNA mismatch repair and 212 homoploymer repeats. These include microsatellite loci recommended by the NCI Consensus Panel for identifying MSI; 2) amelogenin (for gender); and 3) Fusobacterium to detect a putative CRC-associated pathogen in tumor biopsies. At the AACR annual meeting, we expect to present results for deep sequencing (∼1,000X) of the first 1,000 CR tumors, including any preliminarily identified pathways and subtypes that may provide the basis for association testing with germline genetic and lifestyle and environmental risk factors needed for inferring better approaches to prevention and treatment of CRC. Citation Format: Syed H. Zaidi, Catie Grasso, Jasmine Mu, Eve Shinbrot, Marios Giannakis, Charles Connolly, Ivan Borozan, Hermann Brenner, Peter Campbell, Andrew Chan, Jenny Chang-Claude, Mengmeng Du, Vincent Ferretti, Amy French, Charles Fuchs, Steven Gallinger, Levi Garraway, Andrea Gsur, Marc Gunter, Tabitha Harrison, Michael Hoffmeister, Li Hsu, Wen-Yi Huang, Jeroen Huyghe, Mathieu Lemire, Elaine Mardis, John McPherson, Polly Newcomb, Lincoln Stein, Wei Sun, Lee Timms, Quang Trinh, David Wheeler, Christina Yung, Niha Zubair, Shuji Ogino, Stephen Thibodeau, Ulrike Peters, Thomas Hudson. Linking the molecular profile of colorectal tumors to germline genetic and environmental risk factors. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 5221.

Abstract Background: Recent genome-wide association studies (GWAS) have identified 28 susceptibility loci for Lung cancer (LC), 32 loci for smoking behavior, 10 loci for chronic obstructive pulmonary disease (COPD), 63 loci for pulmonary function and related phenotypes, totaling 107 GWAS susceptibility loci (as of November 2014). Given the common variants that have been found to be associated with LC in GWAS, targeted-region exome sequencing provides a cost-effective approach for further investigation of high-priority regions of the genome and has the great potential to identify rare causal variants in GWAS loci. Patients and Methods: Using an extreme phenotype approach, we selected 48 sporadic LC patients reporting heavy smoking histories with the presence of severe COPD in whom environmental factors are considered overwhelming, and 54 familial LC cases from families with at least three first-degree relatives with LC who are likely enriched for genomic signal, to search for the disease-causing rare mutations within the target 107 GWAS loci using exome sequencing data. Results: By focusing on exome profiles of these target 107 loci, we identified two disease-causing rare mutations on 6p21.32 TNXB p. Arg504His and 10q25.1 CCDC147 p.Arg696Cys. The homozygous Arg504His mutation presented in four familial and one sporadic LC cases. Notably, the minor allele frequency (MAF) of this variant in Caucasians from 1000 Genomes is very rare (&lt;0.001). The heterozygous CCDC147 Arg696Cys variant were identified in two familial cases and one sporadic case. The MAF for this variant in the Caucasians from 1000 Genomes is 0.004. Both deleterious changes lead to truncation of the putative proteins, which disrupts the predicted domain and alters the cellular localization. We also observed several other suggestive rare mutations on G7C, TNS1, DBH1, and CHRNA5. Conclusion: Our target exome sequencing results demonstrate novel highly disruptive LC risk-conferring TNXB and CCDC147 mutations. Since the rare variants have appreciable effects on disease risk prediction, families and patients may warrant screening for mutations in TNXB and CCDC147 (if validated) to assess their potential LC risk. Citation Format: Yanhong Liu, Farrah Kheradmand, Michael Scheurer, Caleb Davis, David Wheeler, Edwin Silverman, Shete Sanjay, Spiridon Tsavachidis, Georgina Armstrong, Elena Kupert, Marshall Anderson, Yafang Li, Claudio Pikielny, Joan E. Bailey-Wilson, Ming You, Colette Gaba, Mariza DeAndrade, Diptasri Mandal, Claire Simpson, Susan Pinney, Christopher Amos, Margaret Spitz. Target exome sequencing for disease-causing rare mutations in familial and sporadic lung cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4600. doi:10.1158/1538-7445.AM2015-4600

Abstract Background: Hepatocellular carcinomas (HCCs) are heterogenous primary liver cancers with pleomorphic histologic features. We present a comprehensive genetic, genomic, transcriptomic and proteomic analysis of 203 HCCs in The Cancer Genome Atlas Program (TCGA) on behalf of the TCGA Liver analysis working group to provide a deeper understanding of the molecular features of HCCs, to classify them in a clinically-relevant manner, and to provide a public resource that identifies potential targets for emerging therapies.Results: Mutation analyses confirmed the presence of TERT promoter mutations in 40% of HCCs. The analysis also revealed the first reported germline mutation in the TERT promoter in a 30 year old male, possibly reflecting an inherited predisposition to HCC. TERT promoter mutations were enriched in hepatitis C virus infected subjects. Analyses of hepatitis B virus (HBV)-induced HCCs demonstrated recurrent integrations of HBV in the TERT locus. Thus, HBV integration may substitute for TERT promoter mutations as an alternate mechanism of telomerase activation. Other tumor suppressor genes inactivated by mutation included the cell cycle genes TP53 (32% of cases), RB1 (8%), and PTEN (3%), the Wnt/beta catenin pathway gene AXIN1 (4%), and the chromatin remodeling genes ARID1A (9%), BAP1 (6%), ARID1B (4%), ARID2 (4%), and PBRM1 (3%). The oncogene CTNNB1 was mutated in 27% of HCCs. Mutations were also found in the MTOR pathway genes, MTOR (4%), TSC2 (4%) and TSC1 (3%). The patterns of mutational signatures observed included signatures characteristic of the known carcinogens aristolochic acid and aflatoxin B1. Clustering of mRNA, miRNA, or protein expression, DNA hypermethylation and DNA copy number profiles identified five, four, two, five and three clusters, respectively. Analysis of mRNA expression data identified five groups. In particular, Group 1 was associated with Asian race, higher histologic grade, and micro- and macrovascular invasion. Integrated clustering showed a grouping of HCCs into three superclusters: iCluster 1 had a low frequency of CDKN2A silencing, lack of CTNNB1 mutation, a normal-like methylation pattern, and overexpression of proliferation marker genes; iClusters 2 and 3 had a high frequency of CDKN2A silencing, CTNNB1 mutation, and overexpression of oxidative response genes; iCluster 3 tumors had a higher degree of chromosomal instability, high frequency of 17p loss, and a highly altered methylation profile. Copy number aberrations included gains in the MYC, VEGFA, CCND1, TERT, and MET loci and losses in the TP53, RB1, ERFFI1, CDKN2A, ARID1A and PTEN loci. Discussion: Based on integrated analysis of genome, transcriptome, methylome and proteome we showed that HCCs naturally separate into three to five distinct groups with associated molecular and clinical characteristics. Importantly, we identified a novel germline mutation in the TERT promoter that may be associated with an inherited risk for HCC. Note: This abstract was not presented at the meeting. Citation Format: Lewis R. Roberts, David A. Wheeler. Comprehensive integrative characterization of hepatocellular carcinoma: The TCGA HCC project. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 3745. doi:10.1158/1538-7445.AM2015-3745

Abstract Background There are several suspected environmental risk factors for non-Hodgkin lymphoma (NHL). The associations between NHL and environmental chemical exposures have typically been evaluated for individual chemicals (i.e., one-by-one). Objectives We determined the association between a mixture of 27 correlated chemicals measured in house dust and NHL risk. Methods We conducted a population-based case–control study of NHL in four National Cancer Institute–Surveillance, Epidemiology, and End Results centers—Detroit, Michigan; Iowa; Los Angeles County, California; and Seattle, Washington—from 1998 to 2000. We used weighted quantile sum (WQS) regression to model the association of a mixture of chemicals and risk of NHL. The WQS index was a sum of weighted quartiles for 5 polychlorinated biphenyls (PCBs), 7 polycyclic aromatic hydrocarbons (PAHs), and 15 pesticides. We estimated chemical mixture weights and effects for study sites combined and for each site individually, and also for histologic subtypes of NHL. Results The WQS index was statistically significantly associated with NHL overall [odds ratio (OR) = 1.30; 95% CI: 1.08, 1.56; p = 0.006; for one quartile increase] and in the study sites of Detroit (OR = 1.71; 95% CI: 1.02, 2.92; p = 0.045), Los Angeles (OR = 1.44; 95% CI: 1.00, 2.08; p = 0.049), and Iowa (OR = 1.76; 95% CI: 1.23, 2.53; p = 0.002). The index was marginally statistically significant in Seattle (OR = 1.39; 95% CI: 0.97, 1.99; p = 0.071). The most highly weighted chemicals for predicting risk overall were PCB congener 180 and propoxur. Highly weighted chemicals varied by study site; PCBs were more highly weighted in Detroit, and pesticides were more highly weighted in Iowa. Conclusions An index of chemical mixtures was significantly associated with NHL. Our results show the importance of evaluating chemical mixtures when studying cancer risk. Citation Czarnota J, Gennings C, Colt JS, De Roos AJ, Cerhan JR, Severson RK, Hartge P, Ward MH, Wheeler DC. 2015. Analysis of environmental chemical mixtures and non-Hodgkin lymphoma risk in the NCI-SEER NHL Study. Environ Health Perspect 123:965–970; http://dx.doi.org/10.1289/ehp.1408630

ABSTRACT Purpose We conducted a pilot randomized clinical trial of office‐based active vision therapy for the treatment of childhood amblyopia to determine the feasibility of conducting a full‐scale randomized clinical trial. Methods A training and certification program and manual of procedures were developed to certify therapists to administer a standardized vision therapy program in ophthalmology and optometry offices consisting of weekly visits for 16 weeks. Nineteen children, aged 7 to less than 13 years, with amblyopia (20/40–20/100) were randomly assigned to receive either 2 hours of daily patching with active vision therapy or 2 hours of daily patching with placebo vision therapy. Results Therapists in diverse practice settings were successfully trained and certified to perform standardized vision therapy in strict adherence with protocol. Subjects completed 85% of required weekly in‐office vision therapy visits. Eligibility criteria based on age, visual acuity, and stereoacuity, designed to identify children able to complete a standardized vision therapy program and judged likely to benefit from this treatment, led to a high proportion of screened subjects being judged ineligible, resulting in insufficient recruitment. There were difficulties in retrieving adherence data for the computerized home therapy procedures. Conclusions This study demonstrated that a 16‐week treatment trial of vision therapy was feasible with respect to maintaining protocol adherence; however, recruitment under the proposed eligibility criteria, necessitated by the standardized approach to vision therapy, was not successful. A randomized clinical trial of in‐office vision therapy for the treatment of amblyopia would require broadening of the eligibility criteria and improved methods to gather objective data regarding the home therapy. A more flexible approach that customizes vision therapy based on subject age, visual acuity, and stereopsis might be required to allow enrollment of a broader group of subjects.