@article {144, title = {A cross-population atlas of genetic associations for 220 human phenotypes.}, journal = {Nat Genet}, volume = {53}, year = {2021}, month = {2021 10}, pages = {1415-1424}, abstract = {

Current genome-wide association studies do not yet capture sufficient diversity in populations and scope of phenotypes. To expand an atlas of genetic associations in non-European populations, we conducted 220 deep-phenotype genome-wide association studies (diseases, biomarkers and medication usage) in BioBank Japan (n = 179,000), by incorporating past medical history and text-mining of electronic medical records. Meta-analyses with the UK Biobank and FinnGen (n = 628,000) identified ~5,000 new loci, which improved the resolution of the genomic map of human traits. This atlas elucidated the landscape of pleiotropy as represented by the major histocompatibility complex locus, where we conducted HLA fine-mapping. Finally, we performed statistical decomposition of matrices of phenome-wide summary statistics, and identified latent genetic components, which pinpointed responsible variants and biological mechanisms underlying current disease classifications across populations. The decomposed components enabled genetically informed subtyping of similar diseases (for example, allergic diseases). Our study suggests a potential avenue for hypothesis-free re-investigation of human diseases through genetics.

}, keywords = {ABO Blood-Group System, Biological Specimen Banks, Genetic Association Studies, Genetic Loci, Genetic Pleiotropy, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Major Histocompatibility Complex, Meta-Analysis as Topic, Mutation, Phenotype}, issn = {1546-1718}, doi = {10.1038/s41588-021-00931-x}, author = {Sakaue, Saori and Kanai, Masahiro and Tanigawa, Yosuke and Karjalainen, Juha and Kurki, Mitja and Koshiba, Seizo and Narita, Akira and Konuma, Takahiro and Yamamoto, Kenichi and Akiyama, Masato and Ishigaki, Kazuyoshi and Suzuki, Akari and Suzuki, Ken and Obara, Wataru and Yamaji, Ken and Takahashi, Kazuhisa and Asai, Satoshi and Takahashi, Yasuo and Suzuki, Takao and Shinozaki, Nobuaki and Yamaguchi, Hiroki and Minami, Shiro and Murayama, Shigeo and Yoshimori, Kozo and Nagayama, Satoshi and Obata, Daisuke and Higashiyama, Masahiko and Masumoto, Akihide and Koretsune, Yukihiro and Ito, Kaoru and Terao, Chikashi and Yamauchi, Toshimasa and Komuro, Issei and Kadowaki, Takashi and Tamiya, Gen and Yamamoto, Masayuki and Nakamura, Yusuke and Kubo, Michiaki and Murakami, Yoshinori and Yamamoto, Kazuhiko and Kamatani, Yoichiro and Palotie, Aarno and Rivas, Manuel A and Daly, Mark J and Matsuda, Koichi and Okada, Yukinori} } @article {132, title = {Disentangling selection on genetically correlated polygenic traits via whole-genome genealogies.}, journal = {Am J Hum Genet}, volume = {108}, year = {2021}, month = {2021 02 04}, pages = {219-239}, abstract = {

We present a full-likelihood method to infer polygenic adaptation from DNA sequence variation and GWAS summary statistics to quantify recent transient directional selection acting on a complex trait. Through simulations of polygenic trait architecture evolution and GWASs, we show the method substantially improves power over current methods. We examine the robustness of the method under stratification, uncertainty and bias in marginal effects, uncertainty in the causal SNPs, allelic heterogeneity, negative selection, and low GWAS sample size. The method can quantify selection acting on correlated traits, controlling for pleiotropy even among traits with strong genetic correlation (|r|=80\%) while retaining high power to attribute selection to the causal trait. When the causal trait is excluded from analysis, selection is attributed to its closest proxy. We discuss limitations of the method, cautioning against strongly causal interpretations of the results, and the possibility of undetectable gene-by-environment (GxE) interactions. We apply the method to 56 human polygenic traits, revealing signals of directional selection on pigmentation, life history, glycated hemoglobin (HbA1c), and other traits. We also conduct joint testing of 137 pairs of genetically correlated traits, revealing widespread correlated response acting on these traits (2.6-fold enrichment, p = 1.5~{\texttimes} 10). Signs of selection on some traits previously reported as adaptive (e.g., educational attainment and hair color) are largely attributable to correlated response (p = 2.9~{\texttimes} 10 and 1.7~{\texttimes} 10, respectively). Lastly, our joint test shows antagonistic selection has increased type 2 diabetes risk and decrease HbA1c (p = 1.5~{\texttimes} 10).

}, keywords = {Computer Simulation, Diabetes Mellitus, Type 2, Evolution, Molecular, Gene-Environment Interaction, Genetic Heterogeneity, Genetic Pleiotropy, Genome, Human, Genome-Wide Association Study, Glycated Hemoglobin A, Humans, Models, Genetic, Multifactorial Inheritance, Phenotype, Polymorphism, Single Nucleotide, Sample Size, Selection, Genetic}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2020.12.005}, author = {Stern, Aaron J and Speidel, Leo and Zaitlen, Noah A and Nielsen, Rasmus} } @article {131, title = {Graphical analysis for phenome-wide causal discovery in genotyped population-scale biobanks.}, journal = {Nat Commun}, volume = {12}, year = {2021}, month = {2021 01 13}, pages = {350}, abstract = {

Causal inference via Mendelian randomization requires making strong assumptions about horizontal pleiotropy, where genetic instruments are connected to the outcome not only through the exposure. Here, we present causal Graphical Analysis Using Genetics (cGAUGE), a pipeline that overcomes these limitations using instrument filters with provable properties. This is achievable by identifying conditional independencies while examining multiple traits. cGAUGE also uses ExSep (Exposure-based Separation), a novel test for the existence of causal pathways that does not require selecting instruments. In simulated data we illustrate how cGAUGE can reduce the empirical false discovery rate by up to 30\%, while retaining the majority of true discoveries. On 96 complex traits from 337,198 subjects from the UK Biobank, our results cover expected causal links and many new ones that were previously suggested by correlation-based observational studies. Notably, we identify multiple risk factors for cardiovascular disease, including red blood cell distribution width.

}, keywords = {Biological Specimen Banks, Cardiovascular Diseases, Causality, Computer Simulation, Gene Regulatory Networks, Genetic Pleiotropy, Genetic Variation, Genome-Wide Association Study, Genotype, Humans, Mendelian Randomization Analysis, Models, Theoretical, Multifactorial Inheritance, Phenotype, Risk Factors}, issn = {2041-1723}, doi = {10.1038/s41467-020-20516-2}, author = {Amar, David and Sinnott-Armstrong, Nasa and Ashley, Euan A and Rivas, Manuel A} } @article {133, title = {RAFFI: Accurate and fast familial relationship inference in large scale biobank studies using RaPID.}, journal = {PLoS Genet}, volume = {17}, year = {2021}, month = {2021 01}, pages = {e1009315}, abstract = {

Inference of relationships from whole-genome genetic data of a cohort is a crucial prerequisite for genome-wide association studies. Typically, relationships are inferred by computing the kinship coefficients (ϕ) and the genome-wide probability of zero IBD sharing (π0) among all pairs of individuals. Current leading methods are based on pairwise comparisons, which may not scale up to very large cohorts (e.g., sample size >1 million). Here, we propose an efficient relationship inference method, RAFFI. RAFFI leverages the efficient RaPID method to call IBD segments first, then estimate the ϕ and π0 from detected IBD segments. This inference is achieved by a data-driven approach that adjusts the estimation based on phasing quality and genotyping quality. Using simulations, we showed that RAFFI is robust against phasing/genotyping errors, admix events, and varying marker densities, and achieves higher accuracy compared to KING, the current leading method, especially for more distant relatives. When applied to the phased UK Biobank data with ~500K individuals, RAFFI is approximately 18 times faster than KING. We expect RAFFI will offer fast and accurate relatedness inference for even larger cohorts.

}, keywords = {Biological Specimen Banks, Genome, Human, Genome-Wide Association Study, Genotyping Techniques, Haplotypes, Humans, Models, Genetic, Pedigree, Polymorphism, Single Nucleotide}, issn = {1553-7404}, doi = {10.1371/journal.pgen.1009315}, author = {Naseri, Ardalan and Shi, Junjie and Lin, Xihong and Zhang, Shaojie and Zhi, Degui} } @article {134, title = {Regulatory genomic circuitry of human disease loci by integrative epigenomics.}, journal = {Nature}, volume = {590}, year = {2021}, month = {2021 02}, pages = {300-307}, abstract = {

Annotating the molecular basis of human disease remains an unsolved challenge, as 93\% of disease loci are non-coding and gene-regulatory annotations are highly incomplete. Here we present EpiMap, a compendium comprising 10,000~epigenomic maps across 800~samples, which we used to define chromatin states, high-resolution enhancers, enhancer modules, upstream regulators and downstream target genes. We used this resource to annotate 30,000~genetic loci that were associated with 540~traits, predicting trait-relevant tissues, putative causal nucleotide variants in enriched tissue enhancers and candidate tissue-specific target genes for each. We partitioned multifactorial traits into tissue-specific contributing factors with distinct functional enrichments and disease comorbidity patterns, and revealed both single-factor monotropic and multifactor pleiotropic loci. Top-scoring loci frequently had multiple predicted driver variants, converging through multiple enhancers with a common target gene, multiple genes in common tissues, or multiple genes and multiple tissues, indicating extensive pleiotropy. Our results demonstrate the importance of dense, rich, high-resolution epigenomic annotations for the investigation of complex traits.

}, keywords = {Chromatin, Disease, Enhancer Elements, Genetic, Epigenesis, Genetic, Epigenomics, Female, Gene Regulatory Networks, Genetic Loci, Genome-Wide Association Study, Humans, Male, Multifactorial Inheritance, Organ Specificity, Reproducibility of Results}, issn = {1476-4687}, doi = {10.1038/s41586-020-03145-z}, author = {Boix, Carles A and James, Benjamin T and Park, Yongjin P and Meuleman, Wouter and Kellis, Manolis} } @article {87, title = {Allelic Heterogeneity at the CRP Locus Identified by Whole-Genome Sequencing in Multi-ancestry Cohorts.}, journal = {Am J Hum Genet}, volume = {106}, year = {2020}, month = {2020 01 02}, pages = {112-120}, abstract = {

Whole-genome sequencing (WGS) can improve assessment of low-frequency and rare variants, particularly in non-European populations that have been underrepresented in existing genomic studies. The genetic determinants of C-reactive protein (CRP), a biomarker of chronic inflammation, have been extensively studied, with existing genome-wide association studies (GWASs) conducted in >200,000 individuals of European ancestry. In order to discover novel loci associated with CRP levels, we examined a multi-ancestry population (n = 23,279) with WGS (\~{}38{\texttimes} coverage) from the Trans-Omics for Precision Medicine (TOPMed) program. We found evidence for eight distinct associations at the CRP locus, including two variants that have not been identified previously (rs11265259 and rs181704186), both of which are non-coding and more common in individuals of African ancestry (\~{}10\% and \~{}1\% minor allele frequency, respectively, and rare or monomorphic in 1000 Genomes populations of East Asian, South Asian, and European ancestry). We show that the minor (G) allele of rs181704186 is associated with lower CRP levels and decreased transcriptional activity and protein binding in~vitro, providing a plausible molecular mechanism for this African ancestry-specific signal. The individuals homozygous for rs181704186-G have a mean CRP level of 0.23~mg/L, in contrast to individuals heterozygous for rs181704186 with mean CRP of 2.97~mg/L and major allele homozygotes with mean CRP of 4.11~mg/L. This study demonstrates the utility of WGS in multi-ethnic populations to drive discovery of complex trait associations of large effect and to identify functional alleles in noncoding regulatory regions.

}, keywords = {African Continental Ancestry Group, Asian Continental Ancestry Group, C-Reactive Protein, Cohort Studies, European Continental Ancestry Group, Gene Frequency, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Linkage Disequilibrium, Polymorphism, Single Nucleotide, Whole Genome Sequencing}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2019.12.002}, author = {Raffield, Laura M and Iyengar, Apoorva K and Wang, Biqi and Gaynor, Sheila M and Spracklen, Cassandra N and Zhong, Xue and Kowalski, Madeline H and Salimi, Shabnam and Polfus, Linda M and Benjamin, Emelia J and Bis, Joshua C and Bowler, Russell and Cade, Brian E and Choi, Won Jung and Comellas, Alejandro P and Correa, Adolfo and Cruz, Pedro and Doddapaneni, Harsha and Durda, Peter and Gogarten, Stephanie M and Jain, Deepti and Kim, Ryan W and Kral, Brian G and Lange, Leslie A and Larson, Martin G and Laurie, Cecelia and Lee, Jiwon and Lee, Seonwook and Lewis, Joshua P and Metcalf, Ginger A and Mitchell, Braxton D and Momin, Zeineen and Muzny, Donna M and Pankratz, Nathan and Park, Cheol Joo and Rich, Stephen S and Rotter, Jerome I and Ryan, Kathleen and Seo, Daekwan and Tracy, Russell P and Viaud-Martinez, Karine A and Yanek, Lisa R and Zhao, Lue Ping and Lin, Xihong and Li, Bingshan and Li, Yun and Dupuis, Jos{\'e}e and Reiner, Alexander P and Mohlke, Karen L and Auer, Paul L} } @article {91, title = {Assessing Digital Phenotyping to Enhance Genetic Studies of Human Diseases.}, journal = {Am J Hum Genet}, volume = {106}, year = {2020}, month = {2020 05 07}, pages = {611-622}, abstract = {

Population-scale biobanks that combine genetic data and high-dimensional phenotyping for a large number of participants provide an exciting opportunity to perform genome-wide association studies (GWAS) to identify genetic variants associated with diverse quantitative traits and diseases. A major challenge for GWAS in population biobanks is ascertaining disease cases from heterogeneous data sources such as hospital records, digital questionnaire responses, or interviews. In this study, we use genetic parameters, including genetic correlation, to evaluate whether GWAS performed using cases in the UK Biobank ascertained from hospital records, questionnaire responses, and family history of disease implicate similar disease genetics across a range of effect sizes. We find that hospital record and questionnaire GWAS largely identify similar genetic effects for many complex phenotypes and that combining together both phenotyping methods improves power to detect genetic associations. We also show that family history GWAS using cases ascertained on family history of disease agrees with combined hospital record and questionnaire GWAS and that family history GWAS has better power to detect genetic associations for some phenotypes. Overall, this work demonstrates that digital phenotyping and unstructured phenotype data can be combined with structured data such as hospital records to identify cases for GWAS in biobanks and improve the ability of such studies to identify genetic associations.

}, keywords = {Asthma, Databases, Factual, Disease, Female, Genetics, Medical, Genome-Wide Association Study, Genotype, Humans, Male, Neoplasms, Phenotype, United Kingdom}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2020.03.007}, author = {DeBoever, Christopher and Tanigawa, Yosuke and Aguirre, Matthew and McInnes, Greg and Lavertu, Adam and Rivas, Manuel A} } @article {100, title = {On the cross-population generalizability of gene expression prediction models.}, journal = {PLoS Genet}, volume = {16}, year = {2020}, month = {2020 08}, pages = {e1008927}, abstract = {

The genetic control of gene expression is a core component of human physiology. For the past several years, transcriptome-wide association studies have leveraged large datasets of linked genotype and RNA sequencing information to create a powerful gene-based test of association that has been used in dozens of studies. While numerous discoveries have been made, the populations in the training data are overwhelmingly of European descent, and little is known about the generalizability of these models to other populations. Here, we test for cross-population generalizability of gene expression prediction models using a dataset of African American individuals with RNA-Seq data in whole blood. We find that the default models trained in large datasets such as GTEx and DGN fare poorly in African Americans, with a notable reduction in prediction accuracy when compared to European Americans. We replicate these limitations in cross-population generalizability using the five populations in the GEUVADIS dataset. Via realistic simulations of both populations and gene expression, we show that accurate cross-population generalizability of transcriptome prediction only arises when eQTL architecture is substantially shared across populations. In contrast, models with non-identical eQTLs showed patterns similar to real-world data. Therefore, generating RNA-Seq data in diverse populations is a critical step towards multi-ethnic utility of gene expression prediction.

}, keywords = {African Americans, Gene Expression Profiling, Genome-Wide Association Study, Humans, Models, Genetic, Quantitative Trait Loci, Reference Standards, RNA-Seq, Transcriptome}, issn = {1553-7404}, doi = {10.1371/journal.pgen.1008927}, author = {Keys, Kevin L and Mak, Angel C Y and White, Marquitta J and Eckalbar, Walter L and Dahl, Andrew W and Mefford, Joel and Mikhaylova, Anna V and Contreras, Mar{\'\i}a G and Elhawary, Jennifer R and Eng, Celeste and Hu, Donglei and Huntsman, Scott and Oh, Sam S and Salazar, Sandra and LeNoir, Michael A and Ye, Jimmie C and Thornton, Timothy A and Zaitlen, Noah and Burchard, Esteban G and Gignoux, Christopher R} } @article {99, title = {Dynamic incorporation of multiple in silico functional annotations empowers rare variant association analysis of large whole-genome sequencing studies at scale.}, journal = {Nat Genet}, volume = {52}, year = {2020}, month = {2020 09}, pages = {969-983}, abstract = {

Large-scale whole-genome sequencing studies have enabled the analysis of rare variants (RVs) associated with complex phenotypes. Commonly used RV association tests have limited scope to leverage variant functions. We propose STAAR (variant-set test for association using annotation information), a scalable and powerful RV association test method that effectively incorporates both variant categories and multiple complementary annotations using a dynamic weighting scheme. For the latter, we introduce {\textquoteright}annotation principal components{\textquoteright}, multidimensional summaries of in silico variant annotations. STAAR accounts for population structure and relatedness and is scalable for analyzing very large cohort and biobank whole-genome sequencing studies of continuous and dichotomous traits. We applied STAAR to identify RVs associated with four lipid traits in 12,316 discovery and 17,822 replication samples from the Trans-Omics for Precision Medicine Program. We discovered and replicated new RV associations, including disruptive missense RVs of NPC1L1 and an intergenic region near APOC1P1 associated with low-density lipoprotein cholesterol.

}, keywords = {Cholesterol, LDL, Computer Simulation, Genetic Predisposition to Disease, Genetic Variation, Genome, Genome-Wide Association Study, Humans, Models, Genetic, Molecular Sequence Annotation, Phenotype, Whole Genome Sequencing}, issn = {1546-1718}, doi = {10.1038/s41588-020-0676-4}, author = {Li, Xihao and Li, Zilin and Zhou, Hufeng and Gaynor, Sheila M and Liu, Yaowu and Chen, Han and Sun, Ryan and Dey, Rounak and Arnett, Donna K and Aslibekyan, Stella and Ballantyne, Christie M and Bielak, Lawrence F and Blangero, John and Boerwinkle, Eric and Bowden, Donald W and Broome, Jai G and Conomos, Matthew P and Correa, Adolfo and Cupples, L Adrienne and Curran, Joanne E and Freedman, Barry I and Guo, Xiuqing and Hindy, George and Irvin, Marguerite R and Kardia, Sharon L R and Kathiresan, Sekar and Khan, Alyna T and Kooperberg, Charles L and Laurie, Cathy C and Liu, X Shirley and Mahaney, Michael C and Manichaikul, Ani W and Martin, Lisa W and Mathias, Rasika A and McGarvey, Stephen T and Mitchell, Braxton D and Montasser, May E and Moore, Jill E and Morrison, Alanna C and O{\textquoteright}Connell, Jeffrey R and Palmer, Nicholette D and Pampana, Akhil and Peralta, Juan M and Peyser, Patricia A and Psaty, Bruce M and Redline, Susan and Rice, Kenneth M and Rich, Stephen S and Smith, Jennifer A and Tiwari, Hemant K and Tsai, Michael Y and Vasan, Ramachandran S and Wang, Fei Fei and Weeks, Daniel E and Weng, Zhiping and Wilson, James G and Yanek, Lisa R and Neale, Benjamin M and Sunyaev, Shamil R and Abecasis, Gon{\c c}alo R and Rotter, Jerome I and Willer, Cristen J and Peloso, Gina M and Natarajan, Pradeep and Lin, Xihong} } @article {121, title = {A fast and scalable framework for large-scale and ultrahigh-dimensional sparse regression with application to the UK Biobank.}, journal = {PLoS Genet}, volume = {16}, year = {2020}, month = {2020 10}, pages = {e1009141}, abstract = {

The UK Biobank is a very large, prospective population-based cohort study across the United Kingdom. It provides unprecedented opportunities for researchers to investigate the relationship between genotypic information and phenotypes of interest. Multiple regression methods, compared with genome-wide association studies (GWAS), have already been showed to greatly improve the prediction performance for a variety of phenotypes. In the high-dimensional settings, the lasso, since its first proposal in statistics, has been proved to be an effective method for simultaneous variable selection and estimation. However, the large-scale and ultrahigh dimension seen in the UK Biobank pose new challenges for applying the lasso method, as many existing algorithms and their implementations are not scalable to large applications. In this paper, we propose a computational framework called batch screening iterative lasso (BASIL) that can take advantage of any existing lasso solver and easily build a scalable solution for very large data, including those that are larger than the memory size. We introduce snpnet, an R package that implements the proposed algorithm on top of glmnet and optimizes for single nucleotide polymorphism (SNP) datasets. It currently supports l1-penalized linear model, logistic regression, Cox model, and also extends to the elastic net with l1/l2 penalty. We demonstrate results on the UK Biobank dataset, where we achieve competitive predictive performance for all four phenotypes considered (height, body mass index, asthma, high cholesterol) using only a small fraction of the variants compared with other established polygenic risk score methods.

}, keywords = {Algorithms, Asthma, Biological Specimen Banks, Body Height, Body Mass Index, Cholesterol, Cohort Studies, Genetics, Population, Genome-Wide Association Study, Genotype, Humans, Logistic Models, Phenotype, Polymorphism, Single Nucleotide, Proportional Hazards Models, United Kingdom}, issn = {1553-7404}, doi = {10.1371/journal.pgen.1009141}, author = {Qian, Junyang and Tanigawa, Yosuke and Du, Wenfei and Aguirre, Matthew and Chang, Chris and Tibshirani, Robert and Rivas, Manuel A and Hastie, Trevor} } @article {74, title = {Identifying causal variants and genes using functional genomics in specialized cell types and contexts.}, journal = {Hum Genet}, volume = {139}, year = {2020}, month = {2020 Jan}, pages = {95-102}, abstract = {

A central goal in human genetics is the identification of variants and genes that influence the risk of polygenic diseases. In the past decade, genome-wide association studies (GWAS) have identified tens of thousands of genetic loci associated with various diseases. Since the majority of such loci lie within non-coding regions and have many candidate variants in linkage disequilibrium, it has been challenging to accurately identify specific causal variants and genes. To aid in their discovery a variety of statistical and experimental approaches have been developed. These approaches often borrow information from functional genomics assays such as ATAC-seq, ChIP-seq and RNA-seq to annotate functional variants and identify regulatory relationships between variants and genes. While such approaches are powerful, given the diversity of cell types and environments, it is paramount to select disease-relevant contexts for follow-up analyses. In this review, we discuss the latest developments, challenges, and best practices for determining the causal mechanisms of polygenic disease risk variants with functional genomics data from specialized cell types.

}, keywords = {Cell Lineage, Genes, Genetic Predisposition to Disease, Genome, Human, Genome-Wide Association Study, Genomics, Humans, Polymorphism, Single Nucleotide, Quantitative Trait Loci}, issn = {1432-1203}, doi = {10.1007/s00439-019-02044-2}, author = {Liu, Boxiang and Montgomery, Stephen B} } @article {113, title = {Improving the trans-ancestry portability of polygenic risk scores by prioritizing variants in predicted cell-type-specific regulatory elements.}, journal = {Nat Genet}, volume = {52}, year = {2020}, month = {2020 12}, pages = {1346-1354}, abstract = {

Poor trans-ancestry portability of polygenic risk scores is a consequence of Eurocentric genetic studies and limited knowledge of shared causal variants. Leveraging regulatory annotations may improve portability by prioritizing functional over tagging variants. We constructed a resource of 707 cell-type-specific IMPACT regulatory annotations by aggregating 5,345 epigenetic datasets to predict binding patterns of 142 transcription factors across 245 cell types. We then partitioned the common SNP heritability of 111 genome-wide association study summary statistics of European (average n ≈ 189,000) and East Asian (average n ≈ 157,000) origin. IMPACT annotations captured consistent SNP heritability between populations, suggesting prioritization of shared functional variants. Variant prioritization using IMPACT resulted in increased trans-ancestry portability of polygenic risk scores from Europeans to East Asians across all 21 phenotypes analyzed (49.9\% mean relative increase in R). Our study identifies a crucial role for functional annotations such as IMPACT to improve the trans-ancestry portability of genetic data.

}, keywords = {Asian Continental Ancestry Group, Base Sequence, Computational Biology, Enhancer Elements, Genetic, European Continental Ancestry Group, Gene Expression Regulation, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Models, Genetic, Molecular Sequence Annotation, Multifactorial Inheritance, Polymorphism, Single Nucleotide}, issn = {1546-1718}, doi = {10.1038/s41588-020-00740-8}, author = {Amariuta, Tiffany and Ishigaki, Kazuyoshi and Sugishita, Hiroki and Ohta, Tazro and Koido, Masaru and Dey, Kushal K and Matsuda, Koichi and Murakami, Yoshinori and Price, Alkes L and Kawakami, Eiryo and Terao, Chikashi and Raychaudhuri, Soumya} } @article {128, title = {A phenome-wide association study of 26 mendelian genes reveals phenotypic expressivity of common and rare variants within the general population.}, journal = {PLoS Genet}, volume = {16}, year = {2020}, month = {2020 11}, pages = {e1008802}, abstract = {

The clinical evaluation of a genetic syndrome relies upon recognition of a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. However, individuals displaying a single phenotype of a complex syndrome may not meet criteria for clinical diagnosis or genetic testing. Here, we present a phenome-wide association study (PheWAS) approach to systematically explore the phenotypic expressivity of common and rare alleles in genes associated with four well-described syndromic diseases (Alagille (AS), Marfan (MS), DiGeorge (DS), and Noonan (NS) syndromes) in the general population. Using human phenotype ontology (HPO) terms, we systematically mapped 60 phenotypes related to AS, MS, DS and NS in 337,198 unrelated white British from the UK Biobank (UKBB) based on their hospital admission records, self-administrated questionnaires, and physiological measurements. We performed logistic regression adjusting for age, sex, and the first 5 genetic principal components, for each phenotype and each variant in the target genes (JAG1, NOTCH2 FBN1, PTPN1 and RAS-opathy genes, and genes in the 22q11.2 locus) and performed a gene burden test. Overall, we observed multiple phenotype-genotype correlations, such as the association between variation in JAG1, FBN1, PTPN11 and SOS2 with diastolic and systolic blood pressure; and pleiotropy among multiple variants in syndromic genes. For example, rs11066309 in PTPN11 was significantly associated with a lower body mass index, an increased risk of hypothyroidism and a smaller size for gestational age, all in concordance with NS-related phenotypes. Similarly, rs589668 in FBN1 was associated with an increase in body height and blood pressure, and a reduced body fat percentage as observed in Marfan syndrome. Our findings suggest that the spectrum of associations of common and rare variants in genes involved in syndromic diseases can be extended to individual phenotypes within the general population.

}, keywords = {Alagille Syndrome, Alleles, Biological Variation, Population, DiGeorge Syndrome, European Continental Ancestry Group, Female, Gene Frequency, Genetic Association Studies, Genetic Predisposition to Disease, Genetic Testing, Genetic Variation, Genome-Wide Association Study, Humans, Male, Marfan Syndrome, Noonan Syndrome, Phenotype, Polymorphism, Single Nucleotide, United Kingdom}, issn = {1553-7404}, doi = {10.1371/journal.pgen.1008802}, author = {Tcheandjieu, Catherine and Aguirre, Matthew and Gustafsson, Stefan and Saha, Priyanka and Potiny, Praneetha and Haendel, Melissa and Ingelsson, Erik and Rivas, Manuel A and Priest, James R} } @article {89, title = {A positively selected FBN1 missense variant reduces height in Peruvian individuals.}, journal = {Nature}, volume = {582}, year = {2020}, month = {2020 06}, pages = {234-239}, abstract = {

On average, Peruvian individuals are among the shortest in the world. Here we show that Native American ancestry is associated with reduced height in an ethnically diverse group of Peruvian individuals, and identify a population-specific, missense variant in the FBN1 gene (E1297G) that is significantly associated with lower height. Each copy of the minor allele (frequency of 4.7\%) reduces height by 2.2~cm (4.4~cm in homozygous individuals). To our knowledge, this is the largest effect size known for a common height-associated variant. FBN1 encodes the extracellular matrix protein fibrillin~1, which is a major structural component of microfibrils. We observed less densely packed fibrillin-1-rich microfibrils with irregular edges in the skin of individuals who were homozygous for G1297 compared with individuals who were homozygous for E1297. Moreover, we show that the E1297G locus is under positive selection in non-African populations, and that the E1297 variant shows subtle evidence of positive selection specifically within the Peruvian population. This variant is also significantly more frequent in coastal Peruvian populations than in populations from the Andes or the Amazon, which suggests that short stature might be the result of adaptation to factors that are associated with the coastal environment in Peru.

}, keywords = {Body Height, Female, Fibrillin-1, Gene Frequency, Genome-Wide Association Study, Heredity, Humans, Indians, South American, Male, Microfibrils, Mutation, Missense, Peru, Selection, Genetic}, issn = {1476-4687}, doi = {10.1038/s41586-020-2302-0}, author = {Asgari, Samira and Luo, Yang and Akbari, Ali and Belbin, Gillian M and Li, Xinyi and Harris, Daniel N and Selig, Martin and Bartell, Eric and Calderon, Roger and Slowikowski, Kamil and Contreras, Carmen and Yataco, Rosa and Galea, Jerome T and Jimenez, Judith and Coit, Julia M and Farro{\~n}ay, Chandel and Nazarian, Rosalynn M and O{\textquoteright}Connor, Timothy D and Dietz, Harry C and Hirschhorn, Joel N and Guio, Heinner and Lecca, Leonid and Kenny, Eimear E and Freeman, Esther E and Murray, Megan B and Raychaudhuri, Soumya} } @article {90, title = {Rare protein-altering variants in ANGPTL7 lower intraocular pressure and protect against glaucoma.}, journal = {PLoS Genet}, volume = {16}, year = {2020}, month = {2020 05}, pages = {e1008682}, abstract = {

Protein-altering variants that are protective against human disease provide in vivo validation of therapeutic targets. Here we use genotyping data from UK Biobank (n = 337,151 unrelated White British individuals) and FinnGen (n = 176,899) to conduct a search for protein-altering variants conferring lower intraocular pressure (IOP) and protection against glaucoma. Through rare protein-altering variant association analysis, we find a missense variant in ANGPTL7 in UK Biobank (rs28991009, p.Gln175His, MAF = 0.8\%, genotyped in 82,253 individuals with measured IOP and an independent set of 4,238 glaucoma patients and 250,660 controls) that significantly lowers IOP (β = -0.53 and -0.67 mmHg for heterozygotes, -3.40 and -2.37 mmHg for homozygotes, P = 5.96 x 10-9 and 1.07 x 10-13 for corneal compensated and Goldman-correlated IOP, respectively) and is associated with 34\% reduced risk of glaucoma (P = 0.0062). In FinnGen, we identify an ANGPTL7 missense variant at a greater than 50-fold increased frequency in Finland compared with other populations (rs147660927, p.Arg220Cys, MAF Finland = 4.3\%), which was genotyped in 6,537 glaucoma patients and 170,362 controls and is associated with a 29\% lower glaucoma risk (P = 1.9 x 10-12 for all glaucoma types and also protection against its subtypes including exfoliation, primary open-angle, and primary angle-closure). We further find three rarer variants in UK Biobank, including a protein-truncating variant, which confer a strong composite lowering of IOP (P = 0.0012 and 0.24 for Goldman-correlated and corneal compensated IOP, respectively), suggesting the protective mechanism likely resides in the loss of interaction or function. Our results support inhibition or down-regulation of ANGPTL7 as a therapeutic strategy for glaucoma.

}, keywords = {Adult, Aged, Aged, 80 and over, Angiopoietin-like Proteins, Biological Specimen Banks, Case-Control Studies, Cohort Studies, Female, Finland, Gene Frequency, Genetic Predisposition to Disease, Genetics, Population, Genome-Wide Association Study, Glaucoma, Humans, Intraocular Pressure, Loss of Function Mutation, Male, Middle Aged, Mutation, Missense, Polymorphism, Single Nucleotide, United Kingdom}, issn = {1553-7404}, doi = {10.1371/journal.pgen.1008682}, author = {Tanigawa, Yosuke and Wainberg, Michael and Karjalainen, Juha and Kiiskinen, Tuomo and Venkataraman, Guhan and Lemmel{\"a}, Susanna and Turunen, Joni A and Graham, Robert R and Havulinna, Aki S and Perola, Markus and Palotie, Aarno and Daly, Mark J and Rivas, Manuel A} } @article {86, title = {A Robust Method Uncovers Significant Context-Specific Heritability in Diverse Complex Traits.}, journal = {Am J Hum Genet}, volume = {106}, year = {2020}, month = {2020 01 02}, pages = {71-91}, abstract = {

Gene-environment interactions (GxE) can be fundamental in applications ranging from functional genomics to precision medicine and is a conjectured source of substantial heritability. However, unbiased methods to profile GxE genome-wide are nascent and, as we show, cannot accommodate general environment variables, modest sample sizes, heterogeneous noise, and binary traits. To address this gap, we propose a simple, unifying mixed model for gene-environment interaction (GxEMM). In simulations and theory, we show that GxEMM can dramatically improve estimates and eliminate false positives when the assumptions of existing methods fail. We apply GxEMM to a range of human and model organism datasets and find broad evidence of context-specific genetic effects, including GxSex, GxAdversity, and GxDisease interactions across thousands of clinical and molecular phenotypes. Overall, GxEMM is broadly applicable for testing and quantifying polygenic interactions, which can be useful for explaining heritability and invaluable for determining biologically relevant environments.

}, keywords = {Adult, Animals, Computer Simulation, Female, Gene-Environment Interaction, Genetic Markers, Genome-Wide Association Study, Humans, Male, Mental Disorders, Middle Aged, Models, Genetic, Multifactorial Inheritance, Phenomics, Phenotype, Rats}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2019.11.015}, author = {Dahl, Andy and Nguyen, Khiem and Cai, Na and Gandal, Michael J and Flint, Jonathan and Zaitlen, Noah} } @article {110, title = {A Transcriptome-Wide Association Study Identifies Candidate Susceptibility Genes for Pancreatic Cancer Risk.}, journal = {Cancer Res}, volume = {80}, year = {2020}, month = {2020 10 15}, pages = {4346-4354}, abstract = {

Pancreatic cancer is among the most well-characterized cancer types, yet a large proportion of the heritability of pancreatic cancer risk remains unclear. Here, we performed a large transcriptome-wide association study to systematically investigate associations between genetically predicted gene expression in normal pancreas tissue and pancreatic cancer risk. Using data from 305 subjects of mostly European descent in the Genotype-Tissue Expression Project, we built comprehensive genetic models to predict normal pancreas tissue gene expression, modifying the UTMOST (unified test for molecular signatures). These prediction models were applied to the genetic data of 8,275 pancreatic cancer cases and 6,723 controls of European ancestry. Thirteen genes showed an association of genetically predicted expression with pancreatic cancer risk at an FDR <= 0.05, including seven previously reported genes (, and ) and six novel genes not yet reported for pancreatic cancer risk [6q27: OR (95\% confidence interval (CI), 1.54 (1.25-1.89); 13q12.13: OR (95\% CI), 0.78 (0.70-0.88); 14q24.3: OR (95\% CI), 1.35 (1.17-1.56); 17q12: OR (95\% CI), 6.49 (2.96-14.27); 17q21.1: OR (95\% CI), 1.94 (1.45-2.58); and 20p13: OR (95\% CI): 1.41 (1.20-1.66)]. The associations for 10 of these genes (, and ) remained statistically significant even after adjusting for risk SNPs identified in previous genome-wide association study. Collectively, this analysis identified novel candidate susceptibility genes for pancreatic cancer that warrant further investigation. SIGNIFICANCE: A transcriptome-wide association analysis identified seven previously reported and six novel candidate susceptibility genes for pancreatic cancer risk.

}, keywords = {Age Factors, Case-Control Studies, European Continental Ancestry Group, Female, Gene Expression Regulation, Neoplastic, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Male, Models, Genetic, Pancreatic Neoplasms, Polymorphism, Single Nucleotide}, issn = {1538-7445}, doi = {10.1158/0008-5472.CAN-20-1353}, author = {Liu, Duo and Zhou, Dan and Sun, Yanfa and Zhu, Jingjing and Ghoneim, Dalia and Wu, Chong and Yao, Qizhi and Gamazon, Eric R and Cox, Nancy J and Wu, Lang} } @article {88, title = {Whole-Genome and RNA Sequencing Reveal Variation and Transcriptomic Coordination in the Developing Human Prefrontal Cortex.}, journal = {Cell Rep}, volume = {31}, year = {2020}, month = {2020 04 07}, pages = {107489}, abstract = {

Gene expression levels vary across developmental stage, cell type, and region in the brain. Genomic variants also contribute to the variation in expression, and some neuropsychiatric disorder loci may exert their effects through this mechanism. To investigate these relationships, we present BrainVar, a unique resource of paired whole-genome and bulk tissue RNA sequencing from the dorsolateral prefrontal cortex of 176 individuals across prenatal and postnatal development. Here we identify common variants that alter gene expression (expression quantitative trait loci [eQTLs]) constantly across development or predominantly during prenatal or postnatal stages. Both "constant" and "temporal-predominant" eQTLs are enriched for loci associated with neuropsychiatric traits and disorders and colocalize with specific variants. Expression levels of more than 12,000 genes rise or fall in a concerted late-fetal transition, with the transitional genes enriched for cell-type-specific genes and neuropsychiatric risk loci, underscoring the importance of cataloging developmental trajectories in understanding cortical physiology and pathology.

}, keywords = {Base Sequence, Brain, Computational Biology, Databases, Genetic, Genetic Predisposition to Disease, Genetic Variation, Genome-Wide Association Study, Genomics, Humans, Phenotype, Polymorphism, Single Nucleotide, Prefrontal Cortex, Quantitative Trait Loci, Sequence Analysis, RNA, Transcriptome, Whole Exome Sequencing, Whole Genome Sequencing}, issn = {2211-1247}, doi = {10.1016/j.celrep.2020.03.053}, author = {Werling, Donna M and Pochareddy, Sirisha and Choi, Jinmyung and An, Joon-Yong and Sheppard, Brooke and Peng, Minshi and Li, Zhen and Dastmalchi, Claudia and Santpere, Gabriel and Sousa, Andr{\'e} M M and Tebbenkamp, Andrew T N and Kaur, Navjot and Gulden, Forrest O and Breen, Michael S and Liang, Lindsay and Gilson, Michael C and Zhao, Xuefang and Dong, Shan and Klei, Lambertus and Cicek, A Ercument and Buxbaum, Joseph D and Adle-Biassette, Homa and Thomas, Jean-Leon and Aldinger, Kimberly A and O{\textquoteright}Day, Diana R and Glass, Ian A and Zaitlen, Noah A and Talkowski, Michael E and Roeder, Kathryn and State, Matthew W and Devlin, Bernie and Sanders, Stephan J and Sestan, Nenad} } @article {59, title = {Adjusting for Principal Components of Molecular Phenotypes Induces Replicating False Positives.}, journal = {Genetics}, volume = {211}, year = {2019}, month = {2019 04}, pages = {1179-1189}, abstract = {

High-throughput measurements of molecular phenotypes provide an unprecedented opportunity to model cellular processes and their impact on disease. These highly structured datasets are usually strongly confounded, creating false positives and reducing power. This has motivated many approaches based on principal components analysis (PCA) to estimate and correct for confounders, which have become indispensable elements of association tests between molecular phenotypes and both genetic and nongenetic factors. Here, we show that these correction approaches induce a bias, and that it persists for large sample sizes and replicates out-of-sample. We prove this theoretically for PCA by deriving an analytic, deterministic, and intuitive bias approximation. We assess other methods with realistic simulations, which show that perturbing any of several basic parameters can cause false positive rate (FPR) inflation. Our experiments show the bias depends on covariate and confounder sparsity, effect sizes, and their correlation. Surprisingly, when the covariate and confounder have [Formula: see text], standard two-step methods all have [Formula: see text]-fold FPR inflation. Our analysis informs best practices for confounder correction in genomic studies, and suggests many false discoveries have been made and replicated in some differential expression analyses.

}, keywords = {Animals, Genome-Wide Association Study, Humans, Models, Genetic, Phenotype, Principal Component Analysis, Quantitative Trait Loci, Reproducibility of Results}, issn = {1943-2631}, doi = {10.1534/genetics.118.301768}, author = {Dahl, Andy and Guillemot, Vincent and Mefford, Joel and Aschard, Hugues and Zaitlen, Noah} } @article {64, title = {A Bayesian framework that integrates multi-omics data and gene networks predicts risk genes from schizophrenia GWAS data.}, journal = {Nat Neurosci}, volume = {22}, year = {2019}, month = {2019 05}, pages = {691-699}, abstract = {

Genome-wide association studies (GWAS) have identified more than 100 schizophrenia (SCZ)-associated loci, but using these findings to illuminate disease biology remains a challenge. Here we present integrative risk gene selector (iRIGS), a Bayesian framework that integrates multi-omics data and gene networks to infer risk genes in GWAS loci. By applying iRIGS to SCZ GWAS data, we predicted a set of high-confidence risk genes, most of which are not the nearest genes to the GWAS index variants. High-confidence risk genes account for a significantly enriched heritability, as estimated by stratified linkage disequilibrium score regression. Moreover, high-confidence risk genes are predominantly expressed in brain tissues, especially prenatally, and are enriched for targets of approved drugs, suggesting opportunities to reposition existing drugs for SCZ. Thus, iRIGS can leverage accumulating functional genomics and GWAS data to advance our understanding of SCZ etiology and potential therapeutics.

}, keywords = {Animals, Bayes Theorem, Disease Models, Animal, Gene Regulatory Networks, Genetic Predisposition to Disease, Genome-Wide Association Study, Genomics, Humans, Mice, Risk Factors, Schizophrenia}, issn = {1546-1726}, doi = {10.1038/s41593-019-0382-7}, author = {Wang, Quan and Chen, Rui and Cheng, Feixiong and Wei, Qiang and Ji, Ying and Yang, Hai and Zhong, Xue and Tao, Ran and Wen, Zhexing and Sutcliffe, James S and Liu, Chunyu and Cook, Edwin H and Cox, Nancy J and Li, Bingshan} } @article {55, title = {Methods for the Analysis and Interpretation for Rare Variants Associated with Complex Traits.}, journal = {Curr Protoc Hum Genet}, volume = {101}, year = {2019}, month = {2019 04}, pages = {e83}, abstract = {

With the advent of Next Generation Sequencing (NGS) technologies, whole genome and whole exome DNA sequencing has become affordable for routine genetic studies. Coupled with improved genotyping arrays and genotype imputation methodologies, it is increasingly feasible to obtain rare genetic variant information in large datasets. Such datasets allow researchers to gain a more complete understanding of the genetic architecture of complex traits caused by rare variants. State-of-the-art statistical methods for the statistical genetics analysis of sequence-based association, including efficient algorithms for association analysis in biobank-scale datasets, gene-association tests, meta-analysis, fine mapping methods that integrate functional genomic dataset, and phenome-wide association studies (PheWAS), are reviewed here. These methods are expected to be highly useful for next generation statistical genetics analysis in the era of precision medicine. {\textcopyright} 2019 by John Wiley \& Sons, Inc.

}, keywords = {Algorithms, Genetic Predisposition to Disease, Genome, Human, Genome-Wide Association Study, Genotype, High-Throughput Nucleotide Sequencing, Humans, Multifactorial Inheritance, Phenotype, Polymorphism, Single Nucleotide, Whole Exome Sequencing, Whole Genome Sequencing}, issn = {1934-8258}, doi = {10.1002/cphg.83}, author = {Weissenkampen, J Dylan and Jiang, Yu and Eckert, Scott and Jiang, Bibo and Li, Bingshan and Liu, Dajiang J} } @article {51, title = {A multi-task convolutional deep neural network for variant calling in single molecule sequencing.}, journal = {Nat Commun}, volume = {10}, year = {2019}, month = {2019 03 01}, pages = {998}, abstract = {

The accurate identification of DNA sequence variants is an important, but challenging task in genomics. It is particularly difficult for single molecule sequencing, which has a per-nucleotide error rate of ~5-15\%. Meeting this demand, we developed Clairvoyante, a multi-task five-layer convolutional neural network model for predicting variant type (SNP or indel), zygosity, alternative allele and indel length from aligned reads. For the well-characterized NA12878 human sample, Clairvoyante achieves 99.67, 95.78, 90.53\% F1-score on 1KP common variants, and 98.65, 92.57, 87.26\% F1-score for whole-genome analysis, using Illumina, PacBio, and Oxford Nanopore data, respectively. Training on a second human sample shows Clairvoyante is sample agnostic and finds variants in less than 2 h on a standard server. Furthermore, we present 3,135 variants that are missed using Illumina but supported independently by both PacBio and Oxford Nanopore reads. Clairvoyante is available open-source ( https://github.com/aquaskyline/Clairvoyante ), with modules to train, utilize and visualize the model.

}, keywords = {Base Sequence, Computational Biology, DNA Mutational Analysis, Genome, Human, Genome-Wide Association Study, Genomics, Genotype, Genotyping Techniques, Humans, INDEL Mutation, Nanopores, Neural Networks, Computer, Polymorphism, Single Nucleotide, Sequence Analysis, DNA, Software}, issn = {2041-1723}, doi = {10.1038/s41467-019-09025-z}, author = {Luo, Ruibang and Sedlazeck, Fritz J and Lam, Tak-Wah and Schatz, Michael C} } @article {57, title = {Opportunities and challenges for transcriptome-wide association studies.}, journal = {Nat Genet}, volume = {51}, year = {2019}, month = {2019 04}, pages = {592-599}, abstract = {

Transcriptome-wide association studies (TWAS) integrate genome-wide association studies (GWAS) and gene expression datasets to identify gene-trait associations. In this Perspective, we explore properties of TWAS as a potential approach to prioritize causal genes at GWAS loci, by using simulations and case studies of literature-curated candidate causal genes for schizophrenia, low-density-lipoprotein cholesterol and Crohn{\textquoteright}s disease. We explore risk loci where TWAS accurately prioritizes the likely causal gene as well as loci where TWAS prioritizes multiple genes, some likely to be non-causal, owing to sharing of expression quantitative trait loci (eQTL). TWAS is especially prone to spurious prioritization with expression data from non-trait-related tissues or cell types, owing to substantial cross-cell-type variation in expression levels and eQTL strengths. Nonetheless, TWAS prioritizes candidate causal genes more accurately than simple baselines. We suggest best practices for causal-gene prioritization with TWAS and discuss future opportunities for improvement. Our results showcase the strengths and limitations of using eQTL datasets to determine causal genes at GWAS loci.

}, keywords = {Crohn Disease, Genetic Predisposition to Disease, Genetic Variation, Genome-Wide Association Study, Humans, Lipoproteins, LDL, Quantitative Trait Loci, Schizophrenia, Transcriptome}, issn = {1546-1718}, doi = {10.1038/s41588-019-0385-z}, author = {Wainberg, Michael and Sinnott-Armstrong, Nasa and Mancuso, Nicholas and Barbeira, Alvaro N and Knowles, David A and Golan, David and Ermel, Raili and Ruusalepp, Arno and Quertermous, Thomas and Hao, Ke and Bj{\"o}rkegren, Johan L M and Im, Hae Kyung and Pasaniuc, Bogdan and Rivas, Manuel A and Kundaje, Anshul} } @article {48, title = {Quantification of frequency-dependent genetic architectures in 25 UK Biobank traits reveals action of negative selection.}, journal = {Nat Commun}, volume = {10}, year = {2019}, month = {2019 02 15}, pages = {790}, abstract = {

Understanding the role of rare variants is important in elucidating the genetic basis of human disease. Negative selection can cause rare variants to have larger per-allele effect sizes than common variants. Here, we develop a method to estimate the minor allele frequency (MAF) dependence of SNP effect sizes. We use a model in which per-allele effect sizes have variance proportional to [p(1 - p)], where p is the MAF and negative values of α imply larger effect sizes for rare variants. We estimate α for 25 UK Biobank diseases and complex traits. All traits produce negative α estimates, with best-fit mean of -0.38 (s.e. 0.02) across traits. Despite larger rare variant effect sizes, rare variants (MAF < 1\%) explain less than 10\% of total SNP-heritability for most traits analyzed. Using evolutionary modeling and forward simulations, we validate the α model of MAF-dependent trait effects and assess plausible values of relevant evolutionary parameters.

}, keywords = {Algorithms, Alleles, Biological Specimen Banks, Gene Frequency, Genome-Wide Association Study, Genotype, Humans, Models, Genetic, Polymorphism, Single Nucleotide, Quantitative Trait, Heritable, Selection, Genetic, United Kingdom}, issn = {2041-1723}, doi = {10.1038/s41467-019-08424-6}, author = {Schoech, Armin P and Jordan, Daniel M and Loh, Po-Ru and Gazal, Steven and O{\textquoteright}Connor, Luke J and Balick, Daniel J and Palamara, Pier F and Finucane, Hilary K and Sunyaev, Shamil R and Price, Alkes L} } @article {35, title = {A common loss-of-function variant is associated with lower vitamin B concentration in African Americans.}, journal = {Blood}, volume = {131}, year = {2018}, month = {2018 06 21}, pages = {2859-2863}, keywords = {African Americans, Female, Genome-Wide Association Study, Humans, Loss of Function Mutation, Male, Polymorphism, Single Nucleotide, Transcobalamins, Vitamin B 12, Vitamin B Deficiency}, issn = {1528-0020}, doi = {10.1182/blood-2018-03-841023}, author = {Hu, Yao and Raffield, Laura M and Polfus, Linda M and Moscati, Arden and Nadkarni, Girish and Preuss, Michael H and Zhong, Xue and Wei, Qiang and Rich, Stephen S and Li, Yun and Wilson, James G and Correa, Adolfo and Loos, Ruth J F and Li, Bingshan and Auer, Paul L and Reiner, Alex P} } @article {44, title = {Fine-mapping and functional studies highlight potential causal variants for rheumatoid arthritis and type 1 diabetes.}, journal = {Nat Genet}, volume = {50}, year = {2018}, month = {2018 10}, pages = {1366-1374}, abstract = {

To define potentially causal variants for autoimmune disease, we fine-mapped 76 rheumatoid arthritis (11,475 cases, 15,870 controls) and type 1 diabetes loci (9,334 cases, 11,111 controls). After sequencing 799 1-kilobase regulatory (H3K4me3) regions within these loci in 568 individuals, we observed accurate imputation for 89\% of common variants. We defined credible sets of <=5 causal variants at 5 rheumatoid arthritis and 10 type 1 diabetes loci. We identified potentially causal missense variants at DNASE1L3, PTPN22, SH2B3, and TYK2, and noncoding variants at MEG3, CD28-CTLA4, and IL2RA. We also identified potential candidate causal variants at SIRPG and TNFAIP3. Using functional assays, we confirmed allele-specific protein binding and differential enhancer activity for three variants: the CD28-CTLA4 rs117701653 SNP, MEG3 rs34552516 indel, and TNFAIP3 rs35926684 indel.

}, keywords = {Alleles, Arthritis, Rheumatoid, Case-Control Studies, CD28 Antigens, Chromosome Mapping, CTLA-4 Antigen, Diabetes Mellitus, Type 1, Gene Frequency, Genetic Loci, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Jurkat Cells, Mutation, Polymorphism, Single Nucleotide, Quantitative Trait Loci, RNA, Long Noncoding, Tumor Necrosis Factor alpha-Induced Protein 3}, issn = {1546-1718}, doi = {10.1038/s41588-018-0216-7}, author = {Westra, Harm-Jan and Mart{\'\i}nez-Bonet, Marta and Onengut-Gumuscu, Suna and Lee, Annette and Luo, Yang and Teslovich, Nikola and Worthington, Jane and Martin, Javier and Huizinga, Tom and Klareskog, Lars and Rantapaa-Dahlqvist, Solbritt and Chen, Wei-Min and Quinlan, Aaron and Todd, John A and Eyre, Steve and Nigrovic, Peter A and Gregersen, Peter K and Rich, Stephen S and Raychaudhuri, Soumya} } @article {40, title = {Functional architecture of low-frequency variants highlights strength of negative selection across coding and non-coding annotations.}, journal = {Nat Genet}, volume = {50}, year = {2018}, month = {2018 11}, pages = {1600-1607}, abstract = {

Common variant heritability has been widely reported to be concentrated in variants within cell-type-specific non-coding functional annotations, but little is known about low-frequency variant functional architectures. We partitioned the heritability of both low-frequency (0.5\%<= minor allele frequency <5\%) and common (minor allele frequency >=5\%) variants in 40 UK Biobank traits across a broad set of functional annotations. We determined that non-synonymous coding variants explain 17 {\textpm} 1\% of low-frequency variant heritability ([Formula: see text]) versus 2.1 {\textpm} 0.2\% of common variant heritability ([Formula: see text]). Cell-type-specific non-coding annotations that were significantly enriched for [Formula: see text] of corresponding traits were similarly enriched for [Formula: see text] for most traits, but more enriched for brain-related annotations and traits. For example, H3K4me3 marks in brain dorsolateral prefrontal cortex explain 57 {\textpm} 12\% of [Formula: see text] versus 12 {\textpm} 2\% of [Formula: see text] for neuroticism. Forward simulations confirmed that low-frequency variant enrichment depends on the mean selection coefficient of causal variants in the annotation, and can be used to predict effect size variance of causal rare variants (minor allele frequency <0.5\%).

}, keywords = {Alleles, Biological Specimen Banks, European Continental Ancestry Group, Gene Frequency, Genetic Variation, Genetics, Population, Genome-Wide Association Study, Humans, Linkage Disequilibrium, Molecular Sequence Annotation, Open Reading Frames, Polymorphism, Single Nucleotide, Selection, Genetic, United Kingdom}, issn = {1546-1718}, doi = {10.1038/s41588-018-0231-8}, author = {Gazal, Steven and Loh, Po-Ru and Finucane, Hilary K and Ganna, Andrea and Schoech, Armin and Sunyaev, Shamil and Price, Alkes L} } @article {39, title = {Imputation-Aware Tag SNP Selection To Improve Power for Large-Scale, Multi-ethnic Association Studies.}, journal = {G3 (Bethesda)}, volume = {8}, year = {2018}, month = {2018 10 03}, pages = {3255-3267}, abstract = {

The emergence of very large cohorts in genomic research has facilitated a focus on genotype-imputation strategies to power rare variant association. These strategies have benefited from improvements in imputation methods and association tests, however little attention has been paid to ways in which array design can increase rare variant association power. Therefore, we developed a novel framework to select tag SNPs using the reference panel of 26 populations from Phase 3 of the 1000 Genomes Project. We evaluate tag SNP performance mean imputed r at untyped sites using leave-one-out internal validation and standard imputation methods, rather than pairwise linkage disequilibrium. Moving beyond pairwise metrics allows us to account for haplotype diversity across the genome for improve imputation accuracy and demonstrates population-specific biases from pairwise estimates. We also examine array design strategies that contrast multi-ethnic cohorts single populations, and show a boost in performance for the former can be obtained by prioritizing tag SNPs that contribute information across multiple populations simultaneously. Using our framework, we demonstrate increased imputation accuracy for rare variants (frequency < 1\%) by 0.5-3.1\% for an array of one million sites and 0.7-7.1\% for an array of 500,000 sites, depending on the population. Finally, we show how recent explosive growth in non-African populations means tag SNPs capture on average 30\% fewer other variants than in African populations. The unified framework presented here will enable investigators to make informed decisions for the design of new arrays, and help empower the next phase of rare variant association for global health.

}, keywords = {Computational Biology, Databases, Nucleic Acid, Ethnic Groups, Genetic Association Studies, Genetics, Population, Genome-Wide Association Study, Humans, Linkage Disequilibrium, Models, Genetic, Polymorphism, Single Nucleotide, Reproducibility of Results, Selection, Genetic}, issn = {2160-1836}, doi = {10.1534/g3.118.200502}, author = {Wojcik, Genevieve L and Fuchsberger, Christian and Taliun, Daniel and Welch, Ryan and Martin, Alicia R and Shringarpure, Suyash and Carlson, Christopher S and Abecasis, Goncalo and Kang, Hyun Min and Boehnke, Michael and Bustamante, Carlos D and Gignoux, Christopher R and Kenny, Eimear E} } @article {28, title = {Medical relevance of protein-truncating variants across 337,205 individuals in the UK Biobank study.}, journal = {Nat Commun}, volume = {9}, year = {2018}, month = {2018 04 24}, pages = {1612}, abstract = {

Protein-truncating variants can have profound effects on gene function and are critical for clinical genome interpretation and generating therapeutic hypotheses, but their relevance to medical phenotypes has not been systematically assessed. Here, we characterize the effect of 18,228 protein-truncating variants across 135 phenotypes from the UK Biobank and find 27 associations between medical phenotypes and protein-truncating variants in genes outside the major histocompatibility complex. We perform phenome-wide analyses and directly measure the effect in homozygous carriers, commonly referred to as "human knockouts," across medical phenotypes for genes implicated as being protective against disease or associated with at least one phenotype in our study. We find several genes with strong pleiotropic or non-additive effects. Our results illustrate the importance of protein-truncating variants in a variety of diseases.

}, keywords = {Databases, Nucleic Acid, Genome-Wide Association Study, Humans, Phenotype, Proteins, Sequence Deletion, United Kingdom}, issn = {2041-1723}, doi = {10.1038/s41467-018-03910-9}, author = {DeBoever, Christopher and Tanigawa, Yosuke and Lindholm, Malene E and McInnes, Greg and Lavertu, Adam and Ingelsson, Erik and Chang, Chris and Ashley, Euan A and Bustamante, Carlos D and Daly, Mark J and Rivas, Manuel A} } @article {12, title = {Multiple phenotype association tests using summary statistics in genome-wide association studies.}, journal = {Biometrics}, volume = {74}, year = {2018}, month = {2018 Mar}, pages = {165-175}, abstract = {

We study in this article jointly testing the associations of a genetic variant with correlated multiple phenotypes using the summary statistics of individual phenotype analysis from Genome-Wide Association Studies (GWASs). We estimated the between-phenotype correlation matrix using the summary statistics of individual phenotype GWAS analyses, and developed genetic association tests for multiple phenotypes by accounting for between-phenotype correlation without the need to access individual-level data. Since genetic variants often affect multiple phenotypes differently across the genome and the between-phenotype correlation can be arbitrary, we proposed robust and powerful multiple phenotype testing procedures by jointly testing a common mean and a variance component in linear mixed models for summary statistics. We computed the p-values of the proposed tests analytically. This computational advantage makes our methods practically appealing in large-scale GWASs. We performed simulation studies to show that the proposed tests maintained correct type I error rates, and to compare their powers in various settings with the existing methods. We applied the proposed tests to a GWAS Global Lipids Genetics Consortium summary statistics data set and identified additional genetic variants that were missed by the original single-trait analysis.

}, keywords = {Analysis of Variance, Computer Simulation, Genome-Wide Association Study, Humans, Linear Models, Lipids, Models, Genetic, Phenotype}, issn = {1541-0420}, doi = {10.1111/biom.12735}, author = {Liu, Zhonghua and Lin, Xihong} } @article {45, title = {Quantifying the Polygenic Contribution to Cutaneous Squamous Cell Carcinoma Risk.}, journal = {J Invest Dermatol}, volume = {138}, year = {2018}, month = {2018 07}, pages = {1507-1510}, abstract = {

Genetic factors play an important role in cutaneous squamous cell carcinoma risk. Genome-wide association studies have identified 21 single nucleotide polymorphisms associated with cutaneous squamous cell carcinoma risk. Yet no studies have attempted to quantify the contribution of heritability to cutaneous squamous cell carcinoma risk by calculating the population attributable risk using a combination of all discovered genetic variants. Using an additive multi-locus linear logistic model, we determined the cumulative association of these 21 genetic regions to cutaneous squamous cell carcinoma population attributable risk. We computed a multi-locus population attributable risk of 62\%, suggesting that if the effects of all the risk alleles were removed from a population, the cutaneous squamous cell carcinoma risk would drop by 62\%. Using stratified analysis, we also examined the impact of sex on polygenic risk score, and found that men have an increased relative risk throughout the spectrum of the polygenic risk score. Quantifying the impact of genetic predisposition on the proportion of cancer cases can guide future research decisions and public health policy planning.

}, keywords = {Carcinoma, Squamous Cell, Female, Genetic Loci, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Linear Models, Logistic Models, Male, Models, Biological, Polymorphism, Single Nucleotide, Risk Assessment, Risk Factors, Sex Factors, Skin Neoplasms}, issn = {1523-1747}, doi = {10.1016/j.jid.2018.01.031}, author = {Sordillo, Joanne E and Kraft, Peter and Wu, Ann Chen and Asgari, Maryam M} } @article {15, title = {Cohort-specific imputation of gene expression improves prediction of warfarin dose for African Americans.}, journal = {Genome Med}, volume = {9}, year = {2017}, month = {2017 Nov 24}, pages = {98}, abstract = {

BACKGROUND: Genome-wide association studies are useful for discovering genotype-phenotype associations but are limited because they require large cohorts to identify a signal, which can be population-specific. Mapping genetic variation to genes improves power and allows the effects of both protein-coding variation as well as variation in expression to be combined into "gene level" effects.

METHODS: Previous work has shown that warfarin dose can be predicted using information from genetic variation that affects protein-coding regions. Here, we introduce a method that improves dose prediction by integrating tissue-specific gene expression. In particular, we use drug pathways and expression quantitative trait loci knowledge to impute gene expression-on the assumption that differential expression of key pathway genes may impact dose requirement. We focus on 116 genes from the pharmacokinetic and pharmacodynamic pathways of warfarin within training and validation sets comprising both European and African-descent individuals.

RESULTS: We build gene-tissue signatures associated with warfarin dose in a cohort-specific manner and identify a signature of 11 gene-tissue pairs that significantly augments the International Warfarin Pharmacogenetics Consortium dosage-prediction algorithm in both populations.

CONCLUSIONS: Our results demonstrate that imputed expression can improve dose prediction and bridge population-specific compositions. MATLAB code is available at https://github.com/assafgo/warfarin-cohort.

}, keywords = {Anticoagulants, Black or African American, Female, Gene Expression, Genome-Wide Association Study, Humans, Male, Tissue Distribution, Warfarin}, issn = {1756-994X}, doi = {10.1186/s13073-017-0495-0}, author = {Gottlieb, Assaf and Daneshjou, Roxana and DeGorter, Marianne and Bourgeois, Stephane and Svensson, Peter J and Wadelius, Mia and Deloukas, Panos and Montgomery, Stephen B and Altman, Russ B} } @article {21, title = {Covariate selection for association screening in multiphenotype genetic studies.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 Dec}, pages = {1789-1795}, abstract = {

Testing for associations in big data faces the problem of multiple comparisons, wherein true signals are difficult to detect on the background of all associations queried. This difficulty is particularly salient in human genetic association studies, in which phenotypic variation is often driven by numerous variants of small effect. The current strategy to improve power to identify these weak associations consists of applying standard marginal statistical approaches and increasing study sample sizes. Although successful, this approach does not leverage the environmental and genetic factors shared among the multiple phenotypes collected in contemporary cohorts. Here we developed covariates for multiphenotype studies (CMS), an approach that improves power when correlated phenotypes are measured on the same samples. Our analyses of real and simulated data provide direct evidence that correlated phenotypes can be used to achieve increases in power to levels often surpassing the power gained by a twofold increase in sample size.

}, keywords = {Algorithms, Genetic Association Studies, Genetic Variation, Genome-Wide Association Study, Genotype, Humans, Models, Genetic, Multivariate Analysis, Phenotype, Reproducibility of Results, Sample Size}, issn = {1546-1718}, doi = {10.1038/ng.3975}, author = {Aschard, Hugues and Guillemot, Vincent and Vilhjalmsson, Bjarni and Patel, Chirag J and Skurnik, David and Ye, Chun J and Wolpin, Brian and Kraft, Peter and Zaitlen, Noah} } @article {22, title = {Current Challenges and New Opportunities for Gene-Environment Interaction Studies of Complex Diseases.}, journal = {Am J Epidemiol}, volume = {186}, year = {2017}, month = {2017 Oct 01}, pages = {753-761}, abstract = {

Recently, many new approaches, study designs, and statistical and analytical methods have emerged for studying gene-environment interactions (G{\texttimes}Es) in large-scale studies of human populations. There are opportunities in this field, particularly with respect to the incorporation of -omics and next-generation sequencing data and continual improvement in measures of environmental exposures implicated in complex disease outcomes. In a workshop called "Current Challenges and New Opportunities for Gene-Environment Interaction Studies of Complex Diseases," held October 17-18, 2014, by the National Institute of Environmental Health Sciences and the National Cancer Institute in conjunction with the annual American Society of Human Genetics meeting, participants explored new approaches and tools that have been developed in recent years for G{\texttimes}E discovery. This paper highlights current and critical issues and themes in G{\texttimes}E research that need additional consideration, including the improved data analytical methods, environmental exposure assessment, and incorporation of functional data and annotations.

}, keywords = {Disease, Gene-Environment Interaction, Genetic Predisposition to Disease, Genome-Wide Association Study, High-Throughput Nucleotide Sequencing, Humans, Software}, issn = {1476-6256}, doi = {10.1093/aje/kwx227}, author = {McAllister, Kimberly and Mechanic, Leah E and Amos, Christopher and Aschard, Hugues and Blair, Ian A and Chatterjee, Nilanjan and Conti, David and Gauderman, W James and Hsu, Li and Hutter, Carolyn M and Jankowska, Marta M and Kerr, Jacqueline and Kraft, Peter and Montgomery, Stephen B and Mukherjee, Bhramar and Papanicolaou, George J and Patel, Chirag J and Ritchie, Marylyn D and Ritz, Beate R and Thomas, Duncan C and Wei, Peng and Witte, John S} } @article {11, title = {Estimating the selective effects of heterozygous protein-truncating variants from human exome data.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 May}, pages = {806-810}, abstract = {

The evolutionary cost of gene loss is a central question in genetics and has been investigated in model organisms and human cell lines. In humans, tolerance of the loss of one or both functional copies of a gene is related to the gene{\textquoteright}s causal role in disease. However, estimates of the selection and dominance coefficients in humans have been elusive. Here we analyze exome sequence data from 60,706 individuals to make genome-wide estimates of selection against heterozygous loss of gene function. Using this distribution of selection coefficients for heterozygous protein-truncating variants (PTVs), we provide corresponding Bayesian estimates for individual genes. We find that genes under the strongest selection are enriched in embryonic lethal mouse knockouts, Mendelian disease-associated genes, and regulators of transcription. Screening by essentiality, we find a large set of genes under strong selection that are likely to have crucial functions but have not yet been thoroughly characterized.

}, keywords = {Algorithms, Animals, Bayes Theorem, Exome, Gene Frequency, Genetic Predisposition to Disease, Genetic Variation, Genome-Wide Association Study, Genotype, Heterozygote, Humans, Mice, Knockout, Models, Genetic, Mutation, Selection, Genetic, Sequence Analysis, DNA}, issn = {1546-1718}, doi = {10.1038/ng.3831}, author = {Cassa, Christopher A and Weghorn, Donate and Balick, Daniel J and Jordan, Daniel M and Nusinow, David and Samocha, Kaitlin E and O{\textquoteright}Donnell-Luria, Anne and MacArthur, Daniel G and Daly, Mark J and Beier, David R and Sunyaev, Shamil R} } @article {19, title = {A genetic stochastic process model for genome-wide joint analysis of biomarker dynamics and disease susceptibility with longitudinal data.}, journal = {Genet Epidemiol}, volume = {41}, year = {2017}, month = {2017 Nov}, pages = {620-635}, abstract = {

Unraveling the underlying biological mechanisms or pathways behind the effects of genetic variations on complex diseases remains one of the major challenges in the post-GWAS (where GWAS is genome-wide association study) era. To further explore the relationship between genetic variations, biomarkers, and diseases for elucidating underlying pathological mechanism, a huge effort has been placed on examining pleiotropic and gene-environmental interaction effects. We propose a novel genetic stochastic process model (GSPM) that can be applied to GWAS and jointly investigate the genetic effects on longitudinally measured biomarkers and risks of diseases. This model is characterized by more profound biological interpretation and takes into account the dynamics of biomarkers during follow-up when investigating the hazards of a disease. We illustrate the rationale and evaluate the performance of the proposed model through two GWAS. One is to detect single nucleotide polymorphisms (SNPs) having interaction effects on type 2 diabetes (T2D) with body mass index (BMI) and the other is to detect SNPs affecting the optimal BMI level for protecting from T2D. We identified multiple SNPs that showed interaction effects with BMI on T2D, including a novel SNP rs11757677 in the CDKAL1 gene (P = 5.77 {\texttimes} 10 ). We also found a SNP rs1551133 located on 2q14.2 that reversed the effect of BMI on T2D (P = 6.70 {\texttimes} 10 ). In conclusion, the proposed GSPM provides a promising and useful tool in GWAS of longitudinal data for interrogating pleiotropic and interaction effects to gain more insights into the relationship between genes, quantitative biomarkers, and risks of complex diseases.

}, keywords = {Biomarkers, Body Mass Index, Computational Biology, Diabetes Mellitus, Type 2, Female, Gene-Environment Interaction, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Male, Models, Genetic, Models, Statistical, Polymorphism, Single Nucleotide, Risk, Stochastic Processes}, issn = {1098-2272}, doi = {10.1002/gepi.22058}, author = {He, Liang and Zhbannikov, Ilya and Arbeev, Konstantin G and Yashin, Anatoliy I and Kulminski, Alexander M} } @article {26, title = {Improved methods for multi-trait fine mapping of pleiotropic risk loci.}, journal = {Bioinformatics}, volume = {33}, year = {2017}, month = {2017 Jan 15}, pages = {248-255}, abstract = {

MOTIVATION: Genome-wide association studies (GWAS) have identified thousands of regions in the genome that contain genetic variants that increase risk for complex traits and diseases. However, the variants uncovered in GWAS are typically not biologically causal, but rather, correlated to the true causal variant through linkage disequilibrium (LD). To discern the true causal variant(s), a variety of statistical fine-mapping methods have been proposed to prioritize variants for functional validation.

RESULTS: In this work we introduce a new approach, fastPAINTOR, that leverages evidence across correlated traits, as well as functional annotation data, to improve fine-mapping accuracy at pleiotropic risk loci. To improve computational efficiency, we describe an new importance sampling scheme to perform model inference. First, we demonstrate in simulations that by leveraging functional annotation data, fastPAINTOR increases fine-mapping resolution relative to existing methods. Next, we show that jointly modeling pleiotropic risk regions improves fine-mapping resolution compared to standard single trait and pleiotropic fine mapping strategies. We report a reduction in the number of SNPs required for follow-up in order to capture 90\% of the causal variants from 23 SNPs per locus using a single trait to 12 SNPs when fine-mapping two traits simultaneously. Finally, we analyze summary association data from a large-scale GWAS of lipids and show that these improvements are largely sustained in real data.

AVAILABILITY AND IMPLEMENTATION: The fastPAINTOR framework is implemented in the PAINTOR v3.0 package which is publicly available to the research community http://bogdan.bioinformatics.ucla.edu/software/paintor CONTACT: gkichaev@ucla.eduSupplementary information: Supplementary data are available at Bioinformatics online.

}, keywords = {Chromosome Mapping, Genetic Diseases, Inborn, Genetic Loci, Genetic Pleiotropy, Genome-Wide Association Study, Genomics, Humans, Linkage Disequilibrium, Lipid Metabolism, Models, Genetic, Polymorphism, Single Nucleotide, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btw615}, author = {Kichaev, Gleb and Roytman, Megan and Johnson, Ruth and Eskin, Eleazar and Lindstr{\"o}m, Sara and Kraft, Peter and Pasaniuc, Bogdan} } @article {23, title = {Incorporation of Biological Knowledge Into the Study of Gene-Environment Interactions.}, journal = {Am J Epidemiol}, volume = {186}, year = {2017}, month = {2017 Oct 01}, pages = {771-777}, abstract = {

A growing knowledge base of genetic and environmental information has greatly enabled the study of disease risk factors. However, the computational complexity and statistical burden of testing all variants by all environments has required novel study designs and hypothesis-driven approaches. We discuss how incorporating biological knowledge from model organisms, functional genomics, and integrative approaches can empower the discovery of novel gene-environment interactions and discuss specific methodological considerations with each approach. We consider specific examples where the application of these approaches has uncovered effects of gene-environment interactions relevant to drug response and immunity, and we highlight how such improvements enable a greater understanding of the pathogenesis of disease and the realization of precision medicine.

}, keywords = {Animals, Disease, Gene-Environment Interaction, Genome-Wide Association Study, Genomics, Humans, Models, Animal, Sequence Analysis, RNA}, issn = {1476-6256}, doi = {10.1093/aje/kwx229}, author = {Ritchie, Marylyn D and Davis, Joe R and Aschard, Hugues and Battle, Alexis and Conti, David and Du, Mengmeng and Eskin, Eleazar and Fallin, M Daniele and Hsu, Li and Kraft, Peter and Moore, Jason H and Pierce, Brandon L and Bien, Stephanie A and Thomas, Duncan C and Wei, Peng and Montgomery, Stephen B} } @article {24, title = {Lessons Learned From Past Gene-Environment Interaction Successes.}, journal = {Am J Epidemiol}, volume = {186}, year = {2017}, month = {2017 Oct 01}, pages = {778-786}, abstract = {

Genetic and environmental factors are both known to contribute to susceptibility to complex diseases. Therefore, the study of gene-environment interaction (G{\texttimes}E) has been a focus of research for several years. In this article, select examples of G{\texttimes}E from the literature are described to highlight different approaches and underlying principles related to the success of these studies. These examples can be broadly categorized as studies of single metabolism genes, genes in complex metabolism pathways, ranges of exposure levels, functional approaches and model systems, and pharmacogenomics. Some studies illustrated the success of studying exposure metabolism for which candidate genes can be identified. Moreover, some G{\texttimes}E successes depended on the availability of high-quality exposure assessment and longitudinal measures, study populations with a wide range of exposure levels, and the inclusion of ethnically and geographically diverse populations. In several examples, large population sizes were required to detect G{\texttimes}Es. Other examples illustrated the impact of accurately defining scale of the interactions (i.e., additive or multiplicative). Last, model systems and functional approaches provided insights into G{\texttimes}E in several examples. Future studies may benefit from these lessons learned.

}, keywords = {Biomedical Research, Disease, Environmental Exposure, Gene-Environment Interaction, Genetic Predisposition to Disease, Genetic Variation, Genome-Wide Association Study, Humans, Models, Biological}, issn = {1476-6256}, doi = {10.1093/aje/kwx230}, author = {Ritz, Beate R and Chatterjee, Nilanjan and Garcia-Closas, Montserrat and Gauderman, W James and Pierce, Brandon L and Kraft, Peter and Tanner, Caroline M and Mechanic, Leah E and McAllister, Kimberly} } @article {8, title = {Limited statistical evidence for shared genetic effects of eQTLs and autoimmune-disease-associated loci in three major immune-cell types.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 Apr}, pages = {600-605}, abstract = {

Most autoimmune-disease-risk effects identified by genome-wide association studies (GWAS) localize to open chromatin with gene-regulatory activity. GWAS loci are also enriched in expression quantitative trait loci (eQTLs), thus suggesting that most risk variants alter gene expression. However, because causal variants are difficult to identify, and cis-eQTLs occur frequently, it remains challenging to identify specific instances of disease-relevant changes to gene regulation. Here, we used a novel joint likelihood framework with higher resolution than that of previous methods to identify loci where autoimmune-disease risk and an eQTL are driven by a single shared genetic effect. Using eQTLs from three major immune subpopulations, we found shared effects in only \~{}25\% of the loci examined. Thus, we show that a fraction of gene-regulatory changes suggest strong mechanistic hypotheses for disease risk, but we conclude that most risk mechanisms are not likely to involve changes in basal gene expression.

}, keywords = {Autoimmune Diseases, Gene Expression, Gene Regulatory Networks, Genetic Predisposition to Disease, Genome-Wide Association Study, Humans, Immunity, Polymorphism, Single Nucleotide, Quantitative Trait Loci}, issn = {1546-1718}, doi = {10.1038/ng.3795}, author = {Chun, Sung and Casparino, Alexandra and Patsopoulos, Nikolaos A and Croteau-Chonka, Damien C and Raby, Benjamin A and De Jager, Philip L and Sunyaev, Shamil R and Cotsapas, Chris} } @article {9, title = {Population- and individual-specific regulatory variation in Sardinia.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 May}, pages = {700-707}, abstract = {

Genetic studies of complex traits have mainly identified associations with noncoding variants. To further determine the contribution of regulatory variation, we combined whole-genome and transcriptome data for 624 individuals from Sardinia to identify common and rare variants that influence gene expression and splicing. We identified 21,183 expression quantitative trait loci (eQTLs) and 6,768 splicing quantitative trait loci (sQTLs), including 619 new QTLs. We identified high-frequency QTLs and found evidence of selection near genes involved in malarial resistance and increased multiple sclerosis risk, reflecting the epidemiological history of Sardinia. Using family relationships, we identified 809 segregating expression outliers (median z score of 2.97), averaging 13.3 genes per individual. Outlier genes were enriched for proximal rare variants, providing a new approach to study large-effect regulatory variants and their relevance to traits. Our results provide insight into the effects of regulatory variants and their relationship to population history and individual genetic risk.

}, keywords = {Alternative Splicing, Chromosome Mapping, Family Health, Female, Gene Expression Profiling, Genetic Predisposition to Disease, Genetic Variation, Genetics, Population, Genome-Wide Association Study, Genotype, Humans, Italy, Male, Polymorphism, Single Nucleotide, Quantitative Trait Loci, Transcription Initiation Site}, issn = {1546-1718}, doi = {10.1038/ng.3840}, author = {Pala, Mauro and Zappala, Zachary and Marongiu, Mara and Li, Xin and Davis, Joe R and Cusano, Roberto and Crobu, Francesca and Kukurba, Kimberly R and Gloudemans, Michael J and Reinier, Frederic and Berutti, Riccardo and Piras, Maria G and Mulas, Antonella and Zoledziewska, Magdalena and Marongiu, Michele and Sorokin, Elena P and Hess, Gaelen T and Smith, Kevin S and Busonero, Fabio and Maschio, Andrea and Steri, Maristella and Sidore, Carlo and Sanna, Serena and Fiorillo, Edoardo and Bassik, Michael C and Sawcer, Stephen J and Battle, Alexis and Novembre, John and Jones, Chris and Angius, Andrea and Abecasis, Gon{\c c}alo R and Schlessinger, David and Cucca, Francesco and Montgomery, Stephen B} } @article {16, title = {Testing for the indirect effect under the null for genome-wide mediation analyses.}, journal = {Genet Epidemiol}, volume = {41}, year = {2017}, month = {2017 Dec}, pages = {824-833}, abstract = {

Mediation analysis helps researchers assess whether part or all of an exposure{\textquoteright}s effect on an outcome is due to an intermediate variable. The indirect effect can help in designing interventions on the mediator as opposed to the exposure and better understanding the outcome{\textquoteright}s mechanisms. Mediation analysis has seen increased use in genome-wide epidemiological studies to test for an exposure of interest being mediated through a genomic measure such as gene expression or DNA methylation (DNAm). Testing for the indirect effect is challenged by the fact that the null hypothesis is composite. We examined the performance of commonly used mediation testing methods for the indirect effect in genome-wide mediation studies. When there is no association between the exposure and the mediator and no association between the mediator and the outcome, we show that these common tests are overly conservative. This is a case that will arise frequently in genome-wide mediation studies. Caution is hence needed when applying the commonly used mediation tests in genome-wide mediation studies. We evaluated the performance of these methods using simulation studies, and performed an epigenome-wide mediation association study in the Normative Aging Study, analyzing DNAm as a mediator of the effect of pack-years on FEV .

}, keywords = {Basic Helix-Loop-Helix Transcription Factors, DNA Methylation, Epigenomics, Genome-Wide Association Study, Humans, Lung Neoplasms, Models, Genetic, Repressor Proteins}, issn = {1098-2272}, doi = {10.1002/gepi.22084}, author = {Barfield, Richard and Shen, Jincheng and Just, Allan C and Vokonas, Pantel S and Schwartz, Joel and Baccarelli, Andrea A and VanderWeele, Tyler J and Lin, Xihong} } @article {25, title = {VEXOR: an integrative environment for prioritization of functional variants in fine-mapping analysis.}, journal = {Bioinformatics}, volume = {33}, year = {2017}, month = {2017 May 01}, pages = {1389-1391}, abstract = {

MOTIVATION: The identification of the functional variants responsible for observed genome-wide association studies (GWAS) signals is one of the most challenging tasks of the post-GWAS research era. Several tools have been developed to annotate genetic variants by their genomic location and potential functional implications. Each of these tools has its own requirements and internal logic, which forces the user to become acquainted with each interface.

RESULTS: From an awareness of the amount of work needed to analyze a single locus, we have built a flexible, versatile and easy-to-use web interface designed to help in prioritizing variants and predicting their potential functional implications. This interface acts as a single-point of entry linking association results with reference tools and relevant experiments.

AVAILABILITY AND IMPLEMENTATION: VEXOR is an integrative web application implemented through the Shiny framework and available at: http://romix.genome.ulaval.ca/vexor.

CONTACT: arnaud.droit@crchuq.ulaval.ca.

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

}, keywords = {Genome, Human, Genome-Wide Association Study, Genomics, Humans, Polymorphism, Genetic, Sequence Analysis, DNA, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btw826}, author = {Lema{\c c}on, Audrey and Joly Beauparlant, Charles and Soucy, Penny and Allen, Jamie and Easton, Douglas and Kraft, Peter and Simard, Jacques and Droit, Arnaud} }