@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 {146, title = {Sub-genic intolerance, ClinVar, and the epilepsies: A whole-exome sequencing study of 29,165 individuals.}, journal = {Am J Hum Genet}, volume = {108}, year = {2021}, month = {2021 06 03}, pages = {965-982}, abstract = {

Both mild and severe epilepsies are influenced by variants in the same genes, yet an explanation for the resulting phenotypic variation is unknown. As part of the ongoing Epi25 Collaboration, we performed a whole-exome sequencing analysis of 13,487 epilepsy-affected individuals and 15,678 control individuals. While prior Epi25 studies focused on gene-based collapsing analyses, we asked how the pattern of variation within genes differs by epilepsy type. Specifically, we compared the genetic architectures of severe developmental and epileptic encephalopathies (DEEs) and two generally less severe epilepsies, genetic generalized epilepsy and non-acquired focal epilepsy (NAFE). Our gene-based rare variant collapsing analysis used geographic ancestry-based clustering that included broader ancestries than previously possible and revealed novel associations. Using the missense intolerance ratio (MTR), we found that variants in DEE-affected individuals are in significantly more intolerant genic sub-regions than those in NAFE-affected individuals. Only previously reported pathogenic variants absent in available genomic datasets showed a significant burden in epilepsy-affected individuals compared with control individuals, and the ultra-rare pathogenic variants associated with DEE were located in more intolerant genic sub-regions than variants associated with non-DEE epilepsies. MTR filtering improved the yield of ultra-rare pathogenic variants in affected individuals compared with control individuals. Finally, analysis of variants in genes without a disease association revealed a significant burden of loss-of-function variants in the genes most intolerant to such variation, indicating additional epilepsy-risk genes yet to be discovered. Taken together, our study suggests that genic and sub-genic intolerance are critical characteristics for interpreting the effects of variation in genes that influence epilepsy.

}, keywords = {Case-Control Studies, Cohort Studies, Epilepsy, Exome, Genetic Markers, Genetic Predisposition to Disease, Genetic Testing, Genetic Variation, Humans, Phenotype, Whole Exome Sequencing}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2021.04.009} } @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 {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 {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 {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 {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 {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 {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 {36, title = {Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases.}, journal = {Nat Genet}, volume = {50}, year = {2018}, month = {2018 May}, pages = {693-698}, abstract = {

Horizontal pleiotropy occurs when the variant has an effect on disease outside of its effect on the exposure in Mendelian randomization (MR). Violation of the {\textquoteright}no horizontal pleiotropy{\textquoteright} assumption can cause severe bias in MR. We developed the Mendelian randomization pleiotropy residual sum and outlier (MR-PRESSO) test to identify horizontal pleiotropic outliers in multi-instrument summary-level MR testing. We showed using simulations that the MR-PRESSO test is best suited when horizontal pleiotropy occurs in <50\% of instruments. Next we applied the MR-PRESSO test, along with several other MR tests, to complex traits and diseases and found that horizontal pleiotropy (i) was detectable in over 48\% of significant causal relationships in MR; (ii) introduced distortions in the causal estimates in MR that ranged on average from -131\% to 201\%; (iii) induced false-positive causal relationships in up to 10\% of relationships; and (iv) could be corrected in some but not all instances.

}, keywords = {Disease, Genetic Pleiotropy, Genetic Predisposition to Disease, Genetic Variation, Humans}, issn = {1546-1718}, doi = {10.1038/s41588-018-0099-7}, author = {Verbanck, Marie and Chen, Chia-Yen and Neale, Benjamin and Do, Ron} } @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 Oct}, 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 {29, title = {Phenotype risk scores identify patients with unrecognized Mendelian disease patterns.}, journal = {Science}, volume = {359}, year = {2018}, month = {2018 Mar 16}, pages = {1233-1239}, abstract = {

Genetic association studies often examine features independently, potentially missing subpopulations with multiple phenotypes that share a single cause. We describe an approach that aggregates phenotypes on the basis of patterns described by Mendelian diseases. We mapped the clinical features of 1204 Mendelian diseases into phenotypes captured from the electronic health record (EHR) and summarized this evidence as phenotype risk scores (PheRSs). In an initial validation, PheRS distinguished cases and controls of five Mendelian diseases. Applying PheRS to 21,701 genotyped individuals uncovered 18 associations between rare variants and phenotypes consistent with Mendelian diseases. In 16 patients, the rare genetic variants were associated with severe outcomes such as organ transplants. PheRS can augment rare-variant interpretation and may identify subsets of patients with distinct genetic causes for common diseases.

}, keywords = {Databases, Genetic, DNA Mutational Analysis, Electronic Health Records, Exome, Genetic Association Studies, Genetic Diseases, Inborn, Genetic Predisposition to Disease, Genetic Variation, Humans, Phenotype, Risk Factors}, issn = {1095-9203}, doi = {10.1126/science.aal4043}, author = {Bastarache, Lisa and Hughey, Jacob J and Hebbring, Scott and Marlo, Joy and Zhao, Wanke and Ho, Wanting T and Van Driest, Sara L and McGregor, Tracy L and Mosley, Jonathan D and Wells, Quinn S and Temple, Michael and Ramirez, Andrea H and Carroll, Robert and Osterman, Travis and Edwards, Todd and Ruderfer, Douglas and Velez Edwards, Digna R and Hamid, Rizwan and Cogan, Joy and Glazer, Andrew and Wei, Wei-Qi and Feng, QiPing and Brilliant, Murray and Zhao, Zhizhuang J and Cox, Nancy J and Roden, Dan M and Denny, Joshua C} } @article {45, title = {Quantifying the Polygenic Contribution to Cutaneous Squamous Cell Carcinoma Risk.}, journal = {J Invest Dermatol}, volume = {138}, year = {2018}, month = {2018 Jul}, 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 {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 {14, title = {Human Demographic History Impacts Genetic Risk Prediction across Diverse Populations.}, journal = {Am J Hum Genet}, volume = {100}, year = {2017}, month = {2017 Apr 06}, pages = {635-649}, abstract = {

The vast majority of genome-wide association studies (GWASs) are performed in Europeans, and their transferability to other populations is dependent on many factors (e.g., linkage disequilibrium, allele frequencies, genetic architecture). As medical genomics studies become increasingly large and diverse, gaining insights into population history and consequently the transferability of disease risk measurement is critical. Here, we disentangle recent population history in the widely used 1000 Genomes Project reference panel, with an emphasis on populations underrepresented in medical studies. To examine the transferability of single-ancestry GWASs, we used published summary statistics to calculate polygenic risk scores for eight well-studied phenotypes. We identify directional inconsistencies in all scores; for example, height is predicted to decrease with genetic distance from Europeans, despite robust anthropological evidence that West Africans are as tall as Europeans on average. To gain deeper quantitative insights into GWAS transferability, we developed a complex trait coalescent-based simulation framework considering effects of polygenicity, causal allele frequency divergence, and heritability. As expected, correlations between true and inferred risk are typically highest in the population from which summary statistics were derived. We demonstrate that scores inferred from European GWASs are biased by genetic drift in other populations even when choosing the same causal variants and that biases in any direction are possible and unpredictable. This work cautions that summarizing findings from large-scale GWASs may have limited portability to other populations using standard approaches and highlights the need for generalized risk prediction methods and the inclusion of more diverse individuals in medical genomics.

}, keywords = {Americas, Genetic Predisposition to Disease, Genetics, Medical, Genetics, Population, Haplotypes, Human Genome Project, Humans, Multifactorial Inheritance, Racial Groups}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2017.03.004}, author = {Martin, Alicia R and Gignoux, Christopher R and Walters, Raymond K and Wojcik, Genevieve L and Neale, Benjamin M and Gravel, Simon and Daly, Mark J and Bustamante, Carlos D and Kenny, Eimear E} } @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} }