Gabriel, Stacey

Most genetic variants associated with complex traits are hypothesized to regulate gene expression. To understand the genetics underlying gene expression variability, we characterized 14,324 RNA-sequencing samples from the Trans-Omics for Precision Medicine program and performed expression and splicing quantitative trait locus (e/sQTL) analyses in six tissues and cell types, including whole blood ( n = 6454) and lung ( n = 1291). We detected tens of thousands of secondary cis-e/sQTLs, showing that secondary cis-e/sQTL discovery remains unsaturated. We fine-mapped UK Biobank–derived genome-wide association study (GWAS) signals from 164 traits and identified e/sQTL colocalizations for 10,611 GWAS signals, including 7096 that colocalize with secondary e/sQTLs. Our results suggest that even larger e/sQTL analyses will uncover additional secondary e/sQTLs, further benefiting GWAS interpretation.

Abstract Rationale Dysanapsis refers to a mismatch between airway tree caliber and lung size arising early in life. Dysanapsis assessed by computed tomography (CT) is evident by early adulthood and associated with chronic obstructive pulmonary disease (COPD) risk later in life. Objectives By examining the genetic factors associated with CT-assessed dysanapsis, we aimed to elucidate its molecular underpinnings and physiological significance across the lifespan. Methods We performed a genome-wide association study of CT-assessed dysanapsis in 11,951 adults, including individuals from two population-based and two COPD-enriched studies. We applied colocalization analysis to integrate genome-wide association study and gene expression data from whole blood and lung. Genetic variants associated with dysanapsis were combined into a genetic risk score that was applied to examine association with lung function in children from a population-based birth cohort (n = 1,278) and adults from the UKBiobank (n = 369,157). Measurements and Main Results CT-assessed dysanapsis was associated with genetic variants from 21 independent signals in 19 gene regions, implicating HHIP (hedgehog interacting protein), DSP, and NPNT as potential molecular targets based on colocalization of their expression. A higher dysanapsis genetic risk score was associated with obstructive spirometry among 5-year-old children and among adults in the fifth, sixth, and seventh decades of life. Conclusions CT-assessed dysanapsis is associated with variation in genes previously implicated in lung development, and dysanapsis genetic risk is associated with obstructive lung function from early life through older adulthood. Dysanapsis may represent an endophenotype link between the genetic variations associated with lung function and COPD.

Abstract Exonic variants present some of the strongest links between genotype and phenotype. However, these variants can have significant inter-individual pathogenicity differences, known as variable penetrance. In this study, we propose a model where genetically controlled mRNA splicing modulates the pathogenicity of exonic variants. By first cataloging exonic inclusion from RNA-sequencing data in GTEx V8, we find that pathogenic alleles are depleted on highly included exons. Using a large-scale phased whole genome sequencing data from the TOPMed consortium, we observe that this effect may be driven by common splice-regulatory genetic variants, and that natural selection acts on haplotype configurations that reduce the transcript inclusion of putatively pathogenic variants, especially when limiting to haploinsufficient genes. Finally, we test if this effect may be relevant for autism risk using families from the Simons Simplex Collection, but find that splicing of pathogenic alleles has a penetrance reducing effect here as well. Overall, our results indicate that common splice-regulatory variants may play a role in reducing the damaging effects of rare exonic variants.

Background: Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies. Methods: The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants). Results: Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P =7.4×10 −14 ) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P =1.6×10 −14 ), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P =1.8×10 −6 with the secondary filter), while SERPINC1 did not (minimum P =4.4×10 −5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1 , became significant ( P =4.4×10 −7 using all missense variants with minor allele frequency <0.0005). Conclusions: Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.

Abstract Background: Although the spectrum of genomic alterations in primary, treatment-naïve breast tumors has been described, the genomic landscape of HER2+ MBC remains underexplored. Furthermore, tumor genomic alterations that arise after progression on anti-HER2 therapy are largely unknown. Methods: We prospectively collected metastatic tumor biopsies from patients (pts) enrolled on TBCRC003 (NCT00470704), a phase II study evaluating the combination of lapatinib (L) and T in pts with HER2+ MBC who had varying degrees of prior T exposure. We performed WES on baseline metastatic biopsies and normal DNA from 57 pts. In 36 pts, we also performed WES on pre-treatment primary tumors. Tumors were analyzed for point mutations, insertions/deletions, and copy number alterations. Results: Total accrual was 116 pts. 87 pts were registered in one of two efficacy cohorts: Cohort 1 included pts w no prior T for MBC. Pts with prior adjuvant T were included if the interval from last T to 1st recurrence &gt; 12 months. Cohort 2 included pts with 1-2 prior lines of T for MBC or recurrence within 12 months of adjuvant T. An additional 29 pts were enrolled in a biomarker cohort (Cohort 3). Per-protocol efficacy analyses for 85 pts deemed evaluable are shown below: Objective Response RateClinical Benefit RateMedian Time to ProgressionCohort 150% (90% CI 33.8-66.2%)57.5% (95% CI 40.9-73.0)7.4 monthsCohort 222.2% (90% CI 11.2-37.1%)42.2% (95% CI 27.7-57.8)5.3 months As we previously reported (Wagle et al, ASCO 2014), across 57 metastatic tumors, significant recurrently mutated genes were TP53 (n=30; 53%) and PIK3CA (n=19; 33%). The frequency of mutant TP53 and PIK3CA was not significantly different from 119 primary, treatment-naïve HER2+ tumors sequenced in the TCGA study (50%, p=0.8 and 27%, p=0.5, respectively). Recurrent copy number alterations were also similar to TCGA data. Comparing the 38 pts who received any prior T with the 19 pts who did not, there was no significant difference in the incidence of mutant TP53 (53% vs 53%, p=1.0) and PIK3CA (37% vs 26%, p=0.6). We identified mutations in the HER2 kinase domain in 4/38 pts who received prior T (11%), as compared to 0/19 T-naïve pts. HER2 kinase domain mutations have been identified in ∼2% of HER2-negative cancers but &lt;1% of primary HER2+ cancers. 3 of the mutations were L755S, which has been shown to be resistant to L and sensitive to irreversible HER2 inhibitors. The 4th mutation was D742N, a novel kinase domain mutation. None of the 4 pts with HER2 kinase domain mutations had an objective response, though 1 pt had stable disease for 29 weeks. An analysis comparing paired archival primary tumors and baseline metastatic biopsies from 36 pts to identify genomic alterations acquired or enriched in the metastatic tumors will be presented. Conclusions: We present an analysis of the genomic landscape of HER2+ MBC, including comparisons between matched primary tumors and metastatic biopsies. Somatic HER2 kinase mutations in pts with HER2+ MBC treated with prior T suggests that these mutations may be involved in resistance to T, and may predict poor response to additional anti-HER2 therapy with combined L and T. Novel therapeutic approaches may be required for these pts. Citation Format: Nikhil Wagle, Nancy U Lin, Andrea L Richardson, Ignaty Leshchiner, Ingrid A Mayer, Andres Forero-Torres, Timothy J Hobday, Elizabeth C Dees, Rita Nanda, Mothaffar F Rimawi, Hao Guo, William T Barry, Ron Bose, Wei Shen, Antonio C Wolff, Stacey B Gabriel, Levi A Garraway, Eric P Winer, Ian E Krop. Whole exome sequencing (WES) of HER2+ metastatic breast cancer (MBC) from patients with or without prior trastuzumab (T): A correlative analysis of TBCRC003 [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr PD3-5.

A Fine Balance Intellectual and neurological disabilities can arise from diverse developmental aberrations. Novarino et al. (p. 394 , published online 6 September; see the Perspective by Beaudet ) have now determined the genetic basis for one such disorder for a small group of patients. Exome sequencing led to identification of mutations in a kinase BCKDK ( Branched Chain Ketoacid Dehydrogenase Kinase ) that regulates metabolism of branched-chain amino acids such as valine, leucine, and isoleucine. Mice with homozygous mutations in the BCKDK gene showed developmental and neurological abnormalities resembling those in certain mouse autism models. Analysis of transport mechanisms responsible for carrying amino acids across the blood-brain barrier revealed competition between the branched-chain amino acids and large neutral amino acids. Nutritional supplementation with extra branched-chain amino acids in the diet of mice carrying homozygous mutations in the BCKDK gene normalized their phenotype.