Kooperberg, Charles

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 Pancreatic cancer disproportionately affects Black individuals in the United States, but they have limited representation in genetic studies of pancreatic ductal adenocarcinoma (PDAC). To address this gap, we performed admixture mapping and genome-wide association analysis (GWAS) in genetically inferred African ancestry individuals (1,030 cases and 889 controls). Admixture mapping identified three regions with a significantly higher proportion of African ancestry in cases compared to controls (5q33.3, 10p1, 22q12.3). GWAS identified a genome-wide significant association at 5p15.33 ( CLPTM1L , rs383009:T>C, T Allele Frequency=0.51, OR:1.45, P value=1.24×10 -8 ), a locus previously associated with PDAC. Known loci at 5p15.33, 7q32.3, 8q24.21 and 7q25.1 also replicated (P value <0.01). Multi-ancestral fine-mapping identified two potential causal SNPs (rs3830069 and rs2735940) at 5p15.33. Collectively these findings identified novel PDAC risk loci and expanded our understanding of this deadly cancer in underrepresented populations, emphasizing the multifactorial nature of PDAC risk including inherited genetic and non-genetic factors. Statement of Significance To understand how genetic variation contributes to PDAC risk in Black people in North American, we studied individuals of genetically-inferred African ancestry. We identified novel risk loci and differences in the contribution of known loci. This demonstrates that ancestry-informed genetic analyses improve our understanding of PDAC risk and enhances discovery.

Abstract Background The associations between different types of diabetes, characterized by distinct pathophysiology and genetic architecture, and pancreatic ductal adenocarcinoma (PDAC) risk are not understood. Methods We investigated associations of genetic susceptibility to type 2 diabetes (T2D), 8 T2D mechanistic clusters, type 1 diabetes (T1D), and maturity-onset diabetes of the young (MODY) with PDAC risk. We used genome-wide association study (GWAS) summary-level statistics for T2D (242 283 cases, 1 569 734 controls), T1D (18 942 cases, 501 638 controls), and PDAC (10 244 cases and 360 535 controls) in individuals of European ancestry. Results Two-sample Mendelian randomization (MR) using the Robust Adjusted Profile Score (MR-RAPS) method indicated that genetically predicted T2D was associated with PDAC risk (OR = 1.10; 95% CI = 1.05 to 1.15), particularly the T2D obesity (OR = 1.28; 95% CI = 1.15 to 1.42) and lipodystrophy (OR = 1.25; 95% CI = 1.03 to 1.51) clusters. No association was observed for T1D with PDAC risk (OR = 1.01; 95% CI = 0.99 to 1.02). Pathway/gene-set analysis using the summary-based Adaptive Rank Truncated Product (sARTP) method revealed a significant association between the MODY gene-sets and PDAC risk (P = 1.5 × 10−8), which remained after excluding 20 known PDAC GWAS loci (P = 7.6 × 10−4). HNF1A, FOXA3, and HNF4A were the top contributing genes after excluding the previously identified GWAS loci regions. Conclusions Our results from this genetic association study support that T2D, particularly the obesity and lipodystrophy mechanistic clusters, and MODY genomic susceptibility regions play a role in the etiology of PDAC.

Abstract Background: Pancreatic cancer is a leading cause of cancer-related death globally. Risk factors for pancreatic cancer include common genetic variants and potentially heavy alcohol consumption. We assessed if genetic variants modify the association between heavy alcohol consumption and pancreatic cancer risk. Methods: We conducted a genome-wide interaction analysis of single-nucleotide polymorphisms (SNP) by heavy alcohol consumption (more than three drinks per day) for pancreatic cancer in European ancestry populations from genome-wide association studies. Our analysis included 3,707 cases and 4,167 controls from case–control studies and 1,098 cases and 1,162 controls from cohort studies. Fixed-effect meta-analyses were conducted. Results: A potential novel region of association on 10p11.22, lead SNP rs7898449 (interaction P value (Pinteraction) = 5.1 × 10−8 in the meta-analysis; Pinteraction = 2.1 × 10−9 in the case–control studies; Pinteraction = 0.91 in the cohort studies), was identified. An SNP correlated with this lead SNP is an expression quantitative trait locus for the neuropilin 1 gene. Of the 17 genomic regions with genome-wide significant evidence of association with pancreatic cancer in prior studies, we observed suggestive evidence that heavy alcohol consumption modified the association for one SNP near LINC00673, rs11655237 on 17q25.1 (Pinteraction = 0.004). Conclusions: We identified a novel genomic region that may be associated with pancreatic cancer risk in conjunction with heavy alcohol consumption located near an expression quantitative trait locus for neuropilin 1, a protein that plays an important role in the development and progression of pancreatic cancer. Impact: This work can provide insights into the etiology of pancreatic cancer, particularly in heavy drinkers.

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross‐sectional and prospective association analyses of blood‐derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole‐genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P <0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1‐SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04–1.12; P <0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (β=0.091; P =0.11) or in the reverse direction (β=−0.012; P =0.076). Additional bidirectional Mendelian randomization analyses revealed that low‐density lipoprotein cholesterol had a causal effect on mtDNA CN (β=−0.084; P <0.001), but the reverse direction was not significant ( P =0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low‐density lipoprotein cholesterol level ( P =0.52), whereas there was a strong direct causal effect of higher low‐density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (β=−0.092; P <0.001). Conclusions Our findings indicate that high low‐density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

BACKGROUND: Premature menopause is a risk factor for accelerated cardiovascular aging, but underlying mechanisms remain incompletely understood. This study investigated the role of leukocyte telomere length (LTL), a marker of cellular aging and genomic instability, in the association of premature menopause with cardiovascular disease. METHODS: Participants from the UK Biobank and Women’s Health Initiative with complete reproductive history and LTL measurements were included. Primary analyses tested the association between age at menopause and LTL using multivariable-adjusted linear regression. Secondary analyses stratified women by history of gynecologic surgery. Mendelian randomization was used to infer causal relationships between LTL and age at natural menopause. Multivariable-adjusted Cox regression and mediation analyses tested the joint associations of premature menopause and LTL with incident coronary artery disease. RESULTS: This study included 130 254 postmenopausal women (UK Biobank: n=122 224; Women’s Health Initiative: n=8030), of whom 4809 (3.7%) had experienced menopause before age 40. Earlier menopause was associated with shorter LTL (meta-analyzed ß=−0.02 SD/5 years of earlier menopause [95% CI, −0.02 to −0.01]; P =7.2×10 -12 ). This association was stronger and significant in both cohorts for women with natural/spontaneous menopause (meta-analyzed ß=−0.04 SD/5 years of earlier menopause [95% CI, −0.04 to −0.03]; P <2.2×10 -16 ) and was independent of hormone therapy use. Mendelian randomization supported a causal association of shorter genetically predicted LTL with earlier age at natural menopause. LTL and age at menopause were independently associated with incident coronary artery disease, and mediation analyses indicated small but significant mediation effects of LTL in the association of menopausal age with coronary artery disease. CONCLUSIONS: Earlier age at menopause is associated with shorter LTL, especially among women with natural menopause. Accelerated telomere shortening may contribute to the heightened cardiovascular risk associated with premature menopause.

Abstract Background: There are conflicting data on whether nonalcoholic fatty liver disease (NAFLD) is associated with susceptibility to pancreatic cancer. Using Mendelian randomization (MR), we investigated the relationship between genetic predisposition to NAFLD and risk for pancreatic cancer. Methods: Data from genome-wide association studies (GWAS) within the Pancreatic Cancer Cohort Consortium (PanScan; cases n = 5,090, controls n = 8,733) and the Pancreatic Cancer Case Control Consortium (PanC4; cases n = 4,163, controls n = 3,792) were analyzed. We used data on 68 genetic variants with four different MR methods [inverse variance weighting (IVW), MR-Egger, simple median, and penalized weighted median] separately to predict genetic heritability of NAFLD. We then assessed the relationship between each of the four MR methods and pancreatic cancer risk, using logistic regression to calculate ORs and 95% confidence intervals (CI), adjusting for PC risk factors, including obesity and diabetes. Results: No association was found between genetically predicted NAFLD and pancreatic cancer risk in the PanScan or PanC4 samples [e.g., PanScan, IVW OR, 1.04; 95% confidence interval (CI), 0.88–1.22; MR-Egger OR, 0.89; 95% CI, 0.65–1.21; PanC4, IVW OR, 1.07; 95% CI, 0.90–1.27; MR-Egger OR, 0.93; 95% CI, 0.67–1.28]. None of the four MR methods indicated an association between genetically predicted NAFLD and pancreatic cancer risk in either sample. Conclusions: Genetic predisposition to NAFLD is not associated with pancreatic cancer risk. Impact: Given the close relationship between NAFLD and metabolic conditions, it is plausible that any association between NAFLD and pancreatic cancer might reflect host metabolic perturbations (e.g., obesity, diabetes, or metabolic syndrome) and does not necessarily reflect a causal relationship between NAFLD and pancreatic cancer.

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.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

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.