@article {141, title = {Association between Smoking History and Tumor Mutation Burden in Advanced Non-Small Cell Lung Cancer.}, journal = {Cancer Res}, volume = {81}, year = {2021}, month = {2021 May 01}, pages = {2566-2573}, abstract = {

Lung carcinogenesis is a complex and stepwise process involving accumulation of genetic mutations in signaling and oncogenic pathways via interactions with environmental factors and host susceptibility. Tobacco exposure is the leading cause of lung cancer, but its relationship to clinically relevant mutations and the composite tumor mutation burden (TMB) has not been fully elucidated. In this study, we investigated the dose-response relationship in a retrospective observational study of 931 patients treated for advanced-stage non-small cell lung cancer (NSCLC) between April 2013 and February 2020 at the Dana Farber Cancer Institute and Brigham and Women{\textquoteright}s Hospital. Doubling smoking pack-years was associated with increased and less frequent and mutations, whereas doubling smoking-free months was associated with more frequent . In advanced lung adenocarcinoma, doubling smoking pack-years was associated with an increase in TMB, whereas doubling smoking-free months was associated with a decrease in TMB, after controlling for age, gender, and stage. There is a significant dose-response association of smoking history with genetic alterations in cancer-related pathways and TMB in advanced lung adenocarcinoma. SIGNIFICANCE: This study clarifies the relationship between smoking history and clinically relevant mutations in non-small cell lung cancer, revealing the potential of smoking history as a surrogate for tumor mutation burden.

}, issn = {1538-7445}, doi = {10.1158/0008-5472.CAN-20-3991}, author = {Wang, Xinan and Ricciuti, Biagio and Nguyen, Tom and Li, Xihao and Rabin, Michael S and Awad, Mark M and Lin, Xihong and Johnson, Bruce E and Christiani, David C} } @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 {130, title = {Elevated Polygenic Burden for Autism Spectrum Disorder Is Associated With the Broad Autism Phenotype in Mothers of Individuals With Autism Spectrum Disorder.}, journal = {Biol Psychiatry}, volume = {89}, year = {2021}, month = {2021 03 01}, pages = {476-485}, abstract = {

BACKGROUND: Autism spectrum disorder (ASD) is a multifactorial neurodevelopmental disorder that encompasses a complex and heterogeneous set of traits. Subclinical traits that mirror the core features of ASD, referred to as the broad autism phenotype (BAP), have been documented repeatedly in unaffected relatives and are believed to reflect underlying genetic liability to ASD. The BAP may help inform the etiology of ASD by allowing the stratification of families into more phenotypically and etiologically homogeneous subgroups. This study explores polygenic scores related to the BAP.

METHODS: Phenotypic and genotypic information were obtained from 2614 trios from the Simons Simplex Collection. Polygenic scores of ASD (ASD-PGSs) were generated across the sample to determine the shared genetic overlap between the BAP and ASD. Maternal and paternal ASD-PGSs were explored in relation to BAP traits and their child{\textquoteright}s ASD symptomatology.

RESULTS: Maternal pragmatic language was related to child{\textquoteright}s social communicative atypicalities. In fathers, rigid personality was related to increased repetitive behaviors in children. Maternal (but not paternal) ASD-PGSs were related to the pragmatic language and rigid BAP domains.

CONCLUSIONS: Associations emerged between parent and child phenotypes, with more associations emerging in mothers than in fathers. ASD-PGS associations emerged with BAP in mothers only, highlighting the potential for a female protective factor, and implicating the polygenic etiology of ASD-related phenotypes in the BAP.

}, keywords = {Autism Spectrum Disorder, Autistic Disorder, Child, Fathers, Female, Humans, Male, Mothers, Phenotype}, issn = {1873-2402}, doi = {10.1016/j.biopsych.2020.08.029}, author = {Nayar, Kritika and Sealock, Julia M and Maltman, Nell and Bush, Lauren and Cook, Edwin H and Davis, Lea K and Losh, Molly} } @article {139, title = {GWAS of three molecular traits highlights core genes and pathways alongside a highly polygenic background.}, journal = {Elife}, volume = {10}, year = {2021}, month = {2021 Feb 15}, abstract = {

Genome-wide association studies (GWAS) have been used to study the genetic basis of a wide variety of complex diseases and other traits. We describe UK Biobank GWAS results for three molecular traits-urate, IGF-1, and testosterone-with better-understood biology than most other complex traits. We find that many of the most significant hits are readily interpretable. We observe huge enrichment of associations near genes involved in the relevant biosynthesis, transport, or signaling pathways. We show how GWAS data illuminate the biology of each trait, including differences in testosterone regulation between females and males. At the same time, even these molecular traits are highly polygenic, with many thousands of variants spread across the genome contributing to trait variance. In summary, for these three molecular traits we identify strong enrichment of signal in putative core gene sets, even while most of the SNP-based heritability is driven by a massively polygenic background.

}, issn = {2050-084X}, doi = {10.7554/eLife.58615}, author = {Sinnott-Armstrong, Nasa and Naqvi, Sahin and Rivas, Manuel and Pritchard, Jonathan K} } @article {140, title = {Native American Ancestry and Air Pollution Interact to Impact Bronchodilator Response in Puerto Rican Children with Asthma.}, journal = {Ethn Dis}, volume = {31}, year = {2021}, month = {2021 Winter}, pages = {77-88}, abstract = {

Objective: Asthma is the most common chronic disease in children. Short-acting bronchodilator medications are the most commonly prescribed asthma treatment worldwide, regardless of disease severity. Puerto Rican children display the highest asthma morbidity and mortality of any US population. Alarmingly, Puerto Rican children with asthma display poor bronchodilator drug response (BDR). Reduced BDR may explain, in part, the increased asthma morbidity and mortality observed in Puerto Rican children with asthma. Gene-environment interactions may explain a portion of the heritability of BDR. We aimed to identify gene-environment interactions associated with BDR in Puerto Rican children with asthma.

Setting: Genetic, environmental, and psycho-social data from the Genes-environments and Admixture in Latino Americans (GALA II) case-control study.

Participants: Our discovery dataset consisted of 658 Puerto Rican children with asthma; our replication dataset consisted of 514 Mexican American children with asthma.

Main Outcome Measures: We assessed the association of pairwise interaction models with BDR using ViSEN (Visualization of Statistical Epistasis Networks).

Results: We identified a non-linear interaction between Native American genetic ancestry and air pollution significantly associated with BDR in Puerto Rican children with asthma. This interaction was robust to adjustment for age and sex but was not significantly associated with BDR in our replication population.

Conclusions: Decreased Native American ancestry coupled with increased air pollution exposure was associated with increased BDR in Puerto Rican children with asthma. Our study acknowledges BDR{\textquoteright}s phenotypic complexity, and emphasizes the importance of integrating social, environmental, and biological data to further our understanding of complex disease.

}, issn = {1945-0826}, doi = {10.18865/ed.31.1.77}, author = {Contreras, Mar{\'\i}a G and Keys, Kevin and Maga{\~n}a, Joaquin and Goddard, Pag{\'e} C and Risse-Adams, Oona and Zeiger, Andrew M and Mak, Angel C Y and Samedy-Bates, Lesly-Anne and Neophytou, Andreas M and Lee, Eunice and Thakur, Neeta and Elhawary, Jennifer R and Hu, Donglei and Huntsman, Scott and Eng, Celeste and Hu, Ting and Burchard, Esteban G and White, Marquitta J} } @article {145, title = {Population sequencing data reveal a compendium of mutational processes in the human germ line.}, journal = {Science}, volume = {373}, year = {2021}, month = {2021 08 27}, pages = {1030-1035}, abstract = {

Biological mechanisms underlying human germline mutations remain largely unknown. We statistically decompose variation in the rate and spectra of mutations along the genome using volume-regularized nonnegative matrix factorization. The analysis of a sequencing dataset (TOPMed) reveals nine processes that explain the variation in mutation properties between loci. We provide a biological interpretation for seven of these processes. We associate one process with bulky DNA lesions that are resolved asymmetrically with respect to transcription and replication. Two processes track direction of replication fork and replication timing, respectively. We identify a mutagenic effect of active demethylation primarily acting in regulatory regions and a mutagenic effect of long interspersed nuclear elements. We localize a mutagenic process specific to oocytes from population sequencing data. This process appears transcriptionally asymmetric.

}, keywords = {Algorithms, CpG Islands, DNA Damage, DNA Demethylation, DNA Mutational Analysis, DNA Replication, Genetic Variation, Genome, Human, Germ Cells, Germ-Line Mutation, Humans, Long Interspersed Nucleotide Elements, Mutagenesis, Oocytes, Transcription, Genetic}, issn = {1095-9203}, doi = {10.1126/science.aba7408}, author = {Seplyarskiy, Vladimir B and Soldatov, Ruslan A and Koch, Evan and McGinty, Ryan J and Goldmann, Jakob M and Hernandez, Ryan D and Barnes, Kathleen and Correa, Adolfo and Burchard, Esteban G and Ellinor, Patrick T and McGarvey, Stephen T and Mitchell, Braxton D and Vasan, Ramachandran S and Redline, Susan and Silverman, Edwin and Weiss, Scott T and Arnett, Donna K and Blangero, John and Boerwinkle, Eric and He, Jiang and Montgomery, Courtney and Rao, D C and Rotter, Jerome I and Taylor, Kent D and Brody, Jennifer A and Chen, Yii-Der Ida and de Las Fuentes, Lisa and Hwu, Chii-Min and Rich, Stephen S and Manichaikul, Ani W and Mychaleckyj, Josyf C and Palmer, Nicholette D and Smith, Jennifer A and Kardia, Sharon L R and Peyser, Patricia A and Bielak, Lawrence F and O{\textquoteright}Connor, Timothy D and Emery, Leslie S and Gilissen, Christian and Wong, Wendy S W and Kharchenko, Peter V and Sunyaev, Shamil} } @article {143, title = {Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.}, journal = {Cell}, volume = {184}, year = {2021}, month = {2021 05 13}, pages = {2633-2648.e19}, abstract = {

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

}, keywords = {Coronary Artery Disease, Diabetes Mellitus, Type 1, Diabetes Mellitus, Type 2, Disease, Gene Expression Profiling, Genetic Variation, Humans, Inflammatory Bowel Diseases, Multifactorial Inheritance, Organ Specificity, Population, Quantitative Trait Loci, RNA, Long Noncoding, Transcriptome}, issn = {1097-4172}, doi = {10.1016/j.cell.2021.03.050}, author = {de Goede, Olivia M and Nachun, Daniel C and Ferraro, Nicole M and Gloudemans, Michael J and Rao, Abhiram S and Smail, Craig and Eulalio, Tiffany Y and Aguet, Francois and Ng, Bernard and Xu, Jishu and Barbeira, Alvaro N and Castel, Stephane E and Kim-Hellmuth, Sarah and Park, YoSon and Scott, Alexandra J and Strober, Benjamin J and Brown, Christopher D and Wen, Xiaoquan and Hall, Ira M and Battle, Alexis and Lappalainen, Tuuli and Im, Hae Kyung and Ardlie, Kristin G and Mostafavi, Sara and Quertermous, Thomas and Kirkegaard, Karla and Montgomery, Stephen B} } @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 {118, title = {Synaptic processes and immune-related pathways implicated in Tourette syndrome.}, journal = {Transl Psychiatry}, volume = {11}, year = {2021}, month = {2021 01 18}, pages = {56}, abstract = {

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

}, issn = {2158-3188}, doi = {10.1038/s41398-020-01082-z}, author = {Tsetsos, Fotis and Yu, Dongmei and Sul, Jae Hoon and Huang, Alden Y and Illmann, Cornelia and Osiecki, Lisa and Darrow, Sabrina M and Hirschtritt, Matthew E and Greenberg, Erica and Muller-Vahl, Kirsten R and Stuhrmann, Manfred and Dion, Yves and Rouleau, Guy A and Aschauer, Harald and Stamenkovic, Mara and Schl{\"o}gelhofer, Monika and Sandor, Paul and Barr, Cathy L and Grados, Marco A and Singer, Harvey S and N{\"o}then, Markus M and Hebebrand, Johannes and Hinney, Anke and King, Robert A and Fernandez, Thomas V and Barta, Csaba and Tarnok, Zsanett and Nagy, Peter and Depienne, Christel and Worbe, Yulia and Hartmann, Andreas and Budman, Cathy L and Rizzo, Renata and Lyon, Gholson J and McMahon, William M and Batterson, James R and Cath, Danielle C and Malaty, Irene A and Okun, Michael S and Berlin, Cheston and Woods, Douglas W and Lee, Paul C and Jankovic, Joseph and Robertson, Mary M and Gilbert, Donald L and Brown, Lawrence W and Coffey, Barbara J and Dietrich, Andrea and Hoekstra, Pieter J and Kuperman, Samuel and Zinner, Samuel H and Wagner, Michael and Knowles, James A and Jeremy Willsey, A and Tischfield, Jay A and Heiman, Gary A and Cox, Nancy J and Freimer, Nelson B and Neale, Benjamin M and Davis, Lea K and Coppola, Giovanni and Mathews, Carol A and Scharf, Jeremiah M and Paschou, Peristera and Barr, Cathy L and Batterson, James R and Berlin, Cheston and Budman, Cathy L and Cath, Danielle C and Coppola, Giovanni and Cox, Nancy J and Darrow, Sabrina and Davis, Lea K and Dion, Yves and Freimer, Nelson B and Grados, Marco A and Greenberg, Erica and Hirschtritt, Matthew E and Huang, Alden Y and Illmann, Cornelia and King, Robert A and Kurlan, Roger and Leckman, James F and Lyon, Gholson J and Malaty, Irene A and Mathews, Carol A and McMahon, William M and Neale, Benjamin M and Okun, Michael S and Osiecki, Lisa and Robertson, Mary M and Rouleau, Guy A and Sandor, Paul and Scharf, Jeremiah M and Singer, Harvey S and Smit, Jan H and Sul, Jae Hoon and Yu, Dongmei and Aschauer, Harald Aschauer Harald and Barta, Csaba and Budman, Cathy L and Cath, Danielle C and Depienne, Christel and Hartmann, Andreas and Hebebrand, Johannes and Konstantinidis, Anastasios and Mathews, Carol A and M{\"u}ller-Vahl, Kirsten and Nagy, Peter and N{\"o}then, Markus M and Paschou, Peristera and Rizzo, Renata and Rouleau, Guy A and Sandor, Paul and Scharf, Jeremiah M and Schl{\"o}gelhofer, Monika and Stamenkovic, Mara and Stuhrmann, Manfred and Tsetsos, Fotis and Tarnok, Zsanett and Wolanczyk, Tomasz and Worbe, Yulia and Brown, Lawrence and Cheon, Keun-Ah and Coffey, Barbara J and Dietrich, Andrea and Fernandez, Thomas V and Garcia-Delgar, Blanca and Gilbert, Donald and Grice, Dorothy E and Hagstr{\o}m, Julie and Hedderly, Tammy and Heiman, Gary A and Heyman, Isobel and Hoekstra, Pieter J and Huyser, Chaim and Kim, Young Key and Kim, Young-Shin and King, Robert A and Koh, Yun-Joo and Kook, Sodahm and Kuperman, Samuel and Leventhal, Bennett L and Madruga-Garrido, Marcos and Mir, Pablo and Morer, Astrid and M{\"u}nchau, Alexander and Plessen, Kerstin J and Roessner, Veit and Shin, Eun-Young and Song, Dong-Ho and Song, Jungeun and Tischfield, Jay A and Willsey, A Jeremy and Zinner, Samuel and Aschauer, Harald and Barr, Cathy L and Barta, Csaba and Batterson, James R and Berlin, Cheston and Brown, Lawrence and Budman, Cathy L and Cath, Danielle C and Coffey, Barbara J and Coppola, Giovanni and Cox, Nancy J and Darrow, Sabrina and Davis, Lea K and Depienne, Christel and Dietrich, Andrea and Dion, Yves and Fernandez, Thomas and Freimer, Nelson B and Gilbert, Donald and Grados, Marco A and Greenberg, Erica and Hartmann, Andreas and Hebebrand, Johannes and Heiman, Gary and Hirschtritt, Matthew E and Hoekstra, Pieter and Huang, Alden Y and Illmann, Cornelia and Jankovic, Joseph and King, Robert A and Kuperman, Samuel and Lee, Paul C and Lyon, Gholson J and Malaty, Irene A and Mathews, Carol A and McMahon, William M and M{\"u}ller-Vahl, Kirsten and Nagy, Peter and Neale, Benjamin M and N{\"o}then, Markus M and Okun, Michael S and Osiecki, Lisa and Paschou, Peristera and Rizzo, Renata and Robertson, Mary M and Rouleau, Guy A and Sandor, Paul and Scharf, Jeremiah M and Schl{\"o}gelhofer, Monika and Singer, Harvey S and Stamenkovic, Mara and Stuhrmann, Manfred and Sul, Jae Hoon and Tarnok, Zsanett and Tischfield, Jay and Tsetsos, Fotis and Willsey, A Jeremy and Woods, Douglas and Worbe, Yulia and Yu, Dongmei and Zinner, Samuel} } @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 {85, title = {A brief history of human disease genetics.}, journal = {Nature}, volume = {577}, year = {2020}, month = {2020 01}, pages = {179-189}, abstract = {

A primary goal of human genetics is to identify DNA sequence variants that influence biomedical traits, particularly those related to the onset and progression of human disease. Over the past 25 years, progress in realizing this objective has been transformed by advances in technology, foundational genomic resources and analytical tools, and by access to vast amounts of genotype and phenotype data. Genetic discoveries have substantially improved~our understanding of the mechanisms responsible for many rare and common diseases and driven development of novel preventative and therapeutic strategies. Medical innovation will increasingly focus on delivering care tailored to individual patterns of genetic predisposition.

}, keywords = {Animals, Genetic Testing, Genetic Variation, Genomics, Genotype, Humans, Phenotype, Rare Diseases}, issn = {1476-4687}, doi = {10.1038/s41586-019-1879-7}, author = {Claussnitzer, Melina and Cho, Judy H and Collins, Rory and Cox, Nancy J and Dermitzakis, Emmanouil T and Hurles, Matthew E and Kathiresan, Sekar and Kenny, Eimear E and Lindgren, Cecilia M and MacArthur, Daniel G and North, Kathryn N and Plon, Sharon E and Rehm, Heidi L and Risch, Neil and Rotimi, Charles N and Shendure, Jay and Soranzo, Nicole and McCarthy, Mark I} } @article {126, title = {A common variant in PNPLA3 is associated with age at diagnosis of NAFLD in patients from a multi-ethnic biobank.}, journal = {J Hepatol}, volume = {72}, year = {2020}, month = {2020 06}, pages = {1070-1081}, abstract = {

BACKGROUND \& AIMS: The Ile138Met variant (rs738409) in the PNPLA3 gene has the largest effect on non-alcoholic fatty liver disease (NAFLD), increasing the risk of progression to severe forms of liver disease. It remains unknown if the variant plays a role in age of NAFLD onset. We aimed to determine if rs738409 impacts on the age of NAFLD diagnosis.

METHODS: We applied a novel natural language processing (NLP) algorithm to a longitudinal electronic health records (EHR) dataset of >27,000 individuals with genetic data from a multi-ethnic biobank, defining NAFLD cases (n~= 1,703) and confirming controls (n~= 8,119). We conducted i) a survival analysis to determine if age at diagnosis differed by rs738409 genotype, ii) a receiver operating characteristics analysis to assess the utility of the rs738409 genotype in discriminating NAFLD cases from controls, and iii) a phenome-wide association study (PheWAS) between rs738409 and 10,095 EHR-derived disease diagnoses.

RESULTS: The PNPLA3 G risk allele was associated with: i) earlier age of NAFLD diagnosis, with the strongest effect in Hispanics (hazard ratio 1.33; 95\% CI 1.15-1.53; p <0.0001) among whom a NAFLD diagnosis was 15\% more likely in risk allele carriers vs. non-carriers; ii) increased NAFLD risk (odds ratio 1.61; 95\% CI 1.349-1.73; p <0.0001), with the strongest effect among Hispanics (odds ratio 1.43; 95\% CI 1.28-1.59; p <0.0001); iii) additional liver diseases in a PheWAS (p <4.95~{\texttimes} 10) where the risk variant also associated with earlier age of diagnosis.

CONCLUSION: Given the role of the rs738409 in NAFLD diagnosis age, our results suggest that stratifying risk within populations known to have an enhanced risk of liver disease, such as Hispanic carriers of the rs738409 variant, would be effective in earlier identification of those who would benefit most from early NAFLD prevention and treatment strategies.

LAY SUMMARY: Despite clear associations between the PNPLA3 rs738409 variant and elevated risk of progression from non-alcoholic fatty liver disease (NAFLD) to more severe forms of liver disease, it remains unknown if PNPLA3 rs738409 plays a role in the age of NAFLD onset. Herein, we found that this risk variant is associated with an earlier age of NAFLD and other liver disease diagnoses; an observation most pronounced in Hispanic Americans. We conclude that PNPLA3 rs738409 could be used to better understand liver disease risk within vulnerable populations and identify patients that may benefit from early prevention strategies.

}, issn = {1600-0641}, doi = {10.1016/j.jhep.2020.01.029}, author = {Walker, Ryan W and Belbin, Gillian M and Sorokin, Elena P and Van Vleck, Tielman and Wojcik, Genevieve L and Moscati, Arden and Gignoux, Christopher R and Cho, Judy and Abul-Husn, Noura S and Nadkarni, Girish and Kenny, Eimear E and Loos, Ruth J F} } @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 {122, title = {Human Demographic History Impacts Genetic Risk Prediction across Diverse Populations.}, journal = {Am J Hum Genet}, volume = {107}, year = {2020}, month = {2020 10 01}, pages = {788-789}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2020.08.020}, 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 {105, title = {Lung Function in African American Children with Asthma Is Associated with Novel Regulatory Variants of the KIT Ligand and Gene-By-Air-Pollution Interaction.}, journal = {Genetics}, volume = {215}, year = {2020}, month = {2020 07}, pages = {869-886}, abstract = {

Baseline lung function, quantified as forced expiratory volume in the first second of exhalation (FEV), is a standard diagnostic criterion used by clinicians to identify and classify lung diseases. Using whole-genome sequencing data from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine project, we identified a novel genetic association with FEV on chromosome 12 in 867 African American children with asthma ( = 1.26 {\texttimes} 10, β = 0.302). Conditional analysis within 1 Mb of the tag signal (rs73429450) yielded one major and two other weaker independent signals within this peak. We explored statistical and functional evidence for all variants in linkage disequilibrium with the three independent signals and yielded nine variants as the most likely candidates responsible for the association with FEV Hi-C data and expression QTL analysis demonstrated that these variants physically interacted with (KIT ligand, also known as ), and their minor alleles were associated with increased expression of the gene in nasal epithelial cells. Gene-by-air-pollution interaction analysis found that the candidate variant rs58475486 interacted with past-year ambient sulfur dioxide exposure ( = 0.003, β = 0.32). This study identified a novel protective genetic association with FEV, possibly mediated through , in African American children with asthma. This is the first study that has identified a genetic association between lung function and , which has established a role in orchestrating allergic inflammation in asthma.

}, issn = {1943-2631}, doi = {10.1534/genetics.120.303231}, author = {Mak, Angel C Y and Sajuthi, Satria and Joo, Jaehyun and Xiao, Shujie and Sleiman, Patrick M and White, Marquitta J and Lee, Eunice Y and Saef, Benjamin and Hu, Donglei and Gui, Hongsheng and Keys, Kevin L and Lurmann, Fred and Jain, Deepti and Abecasis, Goncalo and Kang, Hyun Min and Nickerson, Deborah A and Germer, Soren and Zody, Michael C and Winterkorn, Lara and Reeves, Catherine and Huntsman, Scott and Eng, Celeste and Salazar, Sandra and Oh, Sam S and Gilliland, Frank D and Chen, Zhanghua and Kumar, Rajesh and Mart{\'\i}nez, Fernando D and Wu, Ann Chen and Ziv, Elad and Hakonarson, Hakon and Himes, Blanca E and Williams, L Keoki and Seibold, Max A and Burchard, Esteban G} } @article {103, title = {Multi-Ethnic Genome-Wide Association Study of Decomposed Cardioelectric Phenotypes Illustrates Strategies to Identify and Characterize Evidence of Shared Genetic Effects for Complex Traits.}, journal = {Circ Genom Precis Med}, volume = {13}, year = {2020}, month = {2020 08}, pages = {e002680}, abstract = {

BACKGROUND: We examined how expanding electrocardiographic trait genome-wide association studies to include ancestrally diverse populations, prioritize more precise phenotypic measures, and evaluate evidence for shared genetic effects enabled the detection and characterization of loci.

METHODS: We decomposed 10 seconds, 12-lead electrocardiograms from 34 668 multi-ethnic participants (15\% Black; 30\% Hispanic/Latino) into 6 contiguous, physiologically distinct (P wave, PR segment, QRS interval, ST segment, T wave, and TP segment) and 2 composite, conventional (PR interval and QT interval) interval scale traits and conducted multivariable-adjusted, trait-specific univariate genome-wide association studies using 1000-G imputed single-nucleotide polymorphisms. Evidence of shared genetic effects was evaluated by aggregating meta-analyzed univariate results across the 6 continuous electrocardiographic traits using the combined phenotype adaptive sum of powered scores test.

RESULTS: We identified 6 novels (, and ) and 87 known loci (adaptive sum of powered score test <5{\texttimes}10). Lead single-nucleotide polymorphism rs3211938 at was common in Blacks (minor allele frequency=10\%), near monomorphic in European Americans, and had effects on the QT interval and TP segment that ranked among the largest reported to date for common variants. The other 5 novel loci were observed when evaluating the contiguous but not the composite electrocardiographic traits. Combined phenotype testing did not identify novel electrocardiographic loci unapparent using traditional univariate approaches, although this approach did assist with the characterization of known loci.

CONCLUSIONS: Despite including one-third as many participants as published electrocardiographic trait genome-wide association studies, our study identified 6 novel loci, emphasizing the importance of ancestral diversity and phenotype resolution in this era of ever-growing genome-wide association studies.

}, issn = {2574-8300}, doi = {10.1161/CIRCGEN.119.002680}, author = {Baldassari, Antoine R and Sitlani, Colleen M and Highland, Heather M and Arking, Dan E and Buyske, Steve and Darbar, Dawood and Gondalia, Rahul and Graff, Misa and Guo, Xiuqing and Heckbert, Susan R and Hindorff, Lucia A and Hodonsky, Chani J and Ida Chen, Yii-Der and Kaplan, Robert C and Peters, Ulrike and Post, Wendy and Reiner, Alex P and Rotter, Jerome I and Shohet, Ralph V and Seyerle, Amanda A and Sotoodehnia, Nona and Tao, Ran and Taylor, Kent D and Wojcik, Genevieve L and Yao, Jie and Kenny, Eimear E and Lin, Henry J and Soliman, Elsayed Z and Whitsel, Eric A and North, Kari E and Kooperberg, Charles and Avery, Christy L} } @article {93, title = {Neurophysiological and Genetic Findings in Patients With Juvenile Myoclonic Epilepsy.}, journal = {Front Integr Neurosci}, volume = {14}, year = {2020}, month = {2020}, pages = {45}, abstract = {

Objective: Transcranial magnetic stimulation (TMS), a non-invasive procedure, stimulates the cortex evaluating the central motor pathways. The response is called motor evoked potential (MEP). Polyphasia results when the response crosses the baseline more than twice (zero crossing). Recent research shows MEP polyphasia in patients with generalized genetic epilepsy (GGE) and their first-degree relatives compared with controls. Juvenile Myoclonic Epilepsy (JME), a GGE type, is not well studied regarding polyphasia. In our study, we assessed polyphasia appearance probability with TMS in JME patients, their healthy first-degree relatives and controls. Two genetic approaches were applied to uncover genetic association with polyphasia.

Methods: 20 JME patients, 23 first-degree relatives and 30 controls underwent TMS, obtaining 10-15 MEPs per participant. We evaluated MEP mean number of phases, proportion of MEP trials displaying polyphasia for each subject and variability between groups. Participants underwent whole exome sequencing (WES) via trio-based analysis and two-case scenario. Extensive bioinformatics analysis was applied.

Results: We identified increased polyphasia in patients (85\%) and relatives (70\%) compared to controls (47\%) and significantly higher mean number of zero crossings (i.e., occurrence of phases) (patients 1.49, relatives 1.46, controls 1.22; < 0.05). Trio-based analysis revealed a candidate polymorphism, p.Glu270del,in , in JME patients and their relatives presenting polyphasia. Sanger sequencing analysis in remaining participants showed no significant association. In two-case scenario, a machine learning approach was applied in variants identified from odds ratio analysis and risk prediction scores were obtained for polyphasia. The results revealed 61 variants of which none was associated with polyphasia. Risk prediction scores indeed showed lower probability for non-polyphasic subjects on having polyphasia and higher probability for polyphasic subjects on having polyphasia.

Conclusion: Polyphasia was present in JME patients and relatives in contrast to controls. Although no known clinical symptoms are linked to polyphasia this neurophysiological phenomenon is likely due to common cerebral electrophysiological abnormality. We did not discover direct association between genetic variants obtained and polyphasia. It is likely these genetic traits alone cannot provoke polyphasia, however, this predisposition combined with disturbed brain-electrical activity and tendency to generate seizures may increase the risk of developing polyphasia, mainly in patients and relatives.

}, issn = {1662-5145}, doi = {10.3389/fnint.2020.00045}, author = {Stefani, Stefani and Kousiappa, Ioanna and Nicolaou, Nicoletta and Papathanasiou, Eleftherios S and Oulas, Anastasis and Fanis, Pavlos and Neocleous, Vassos and Phylactou, Leonidas A and Spyrou, George M and Papacostas, Savvas S} } @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 {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 {94, title = {Transcriptomic signatures across human tissues identify functional rare genetic variation.}, journal = {Science}, volume = {369}, year = {2020}, month = {2020 09 11}, abstract = {

Rare genetic variants are abundant across the human genome, and identifying their function and phenotypic impact is a major challenge. Measuring aberrant gene expression has aided in identifying functional, large-effect rare variants (RVs). Here, we expanded detection of genetically driven transcriptome abnormalities by analyzing gene expression, allele-specific expression, and alternative splicing from multitissue RNA-sequencing data, and demonstrate that each signal informs unique classes of RVs. We developed Watershed, a probabilistic model that integrates multiple genomic and transcriptomic signals to predict variant function, validated these predictions in additional cohorts and through experimental assays, and used them to assess RVs in the UK Biobank, the Million Veterans Program, and the Jackson Heart Study. Our results link thousands of RVs to diverse molecular effects and provide evidence to associate RVs affecting the transcriptome with human traits.

}, keywords = {Genetic Variation, Genome, Human, Humans, Multifactorial Inheritance, Organ Specificity, Transcriptome}, issn = {1095-9203}, doi = {10.1126/science.aaz5900}, author = {Ferraro, Nicole M and Strober, Benjamin J and Einson, Jonah and Abell, Nathan S and Aguet, Francois and Barbeira, Alvaro N and Brandt, Margot and Bucan, Maja and Castel, Stephane E and Davis, Joe R and Greenwald, Emily and Hess, Gaelen T and Hilliard, Austin T and Kember, Rachel L and Kotis, Bence and Park, YoSon and Peloso, Gina and Ramdas, Shweta and Scott, Alexandra J and Smail, Craig and Tsang, Emily K and Zekavat, Seyedeh M and Ziosi, Marcello and Ardlie, Kristin G and Assimes, Themistocles L and Bassik, Michael C and Brown, Christopher D and Correa, Adolfo and Hall, Ira and Im, Hae Kyung and Li, Xin and Natarajan, Pradeep and Lappalainen, Tuuli and Mohammadi, Pejman and Montgomery, Stephen B and Battle, Alexis} } @article {109, title = {Type 2 and interferon inflammation regulate SARS-CoV-2 entry factor expression in the airway epithelium.}, journal = {Nat Commun}, volume = {11}, year = {2020}, month = {2020 10 12}, pages = {5139}, abstract = {

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, an emerging virus that utilizes host proteins ACE2 and TMPRSS2 as entry factors. Understanding the factors affecting the pattern and levels of expression of these genes is important for deeper understanding of SARS-CoV-2 tropism and pathogenesis. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci for both ACE2 and TMPRSS2, that vary in frequency across world populations. We find TMPRSS2 is part of a mucus secretory network, highly upregulated by type 2 (T2) inflammation through the action of interleukin-13, and that the interferon response to respiratory viruses highly upregulates ACE2 expression. IL-13 and virus infection mediated effects on ACE2 expression were also observed at the protein level in the airway epithelium. Finally, we define airway responses to common coronavirus infections in children, finding that these infections generate host responses similar to other viral species, including upregulation of IL6 and ACE2. Our results reveal possible mechanisms influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

}, keywords = {Angiotensin-Converting Enzyme 2, Betacoronavirus, Child, Coronavirus Infections, COVID-19, Epithelial Cells, Gene Expression Profiling, Gene Expression Regulation, Genetic Variation, Host-Pathogen Interactions, Humans, Inflammation, Interferons, Interleukin-13, Middle Aged, Nasal Mucosa, Pandemics, Peptidyl-Dipeptidase A, Pneumonia, Viral, SARS-CoV-2, Serine Endopeptidases, Virus Internalization}, issn = {2041-1723}, doi = {10.1038/s41467-020-18781-2}, author = {Sajuthi, Satria P and DeFord, Peter and Li, Yingchun and Jackson, Nathan D and Montgomery, Michael T and Everman, Jamie L and Rios, Cydney L and Pruesse, Elmar and Nolin, James D and Plender, Elizabeth G and Wechsler, Michael E and Mak, Angel C Y and Eng, Celeste and Salazar, Sandra and Medina, Vivian and Wohlford, Eric M and Huntsman, Scott and Nickerson, Deborah A and Germer, Soren and Zody, Michael C and Abecasis, Goncalo and Kang, Hyun Min and Rice, Kenneth M and Kumar, Rajesh and Oh, Sam and Rodriguez-Santana, Jose and Burchard, Esteban G and Seibold, Max A} } @article {115, title = {Type 2 and interferon inflammation strongly regulate SARS-CoV-2 related gene expression in the airway epithelium.}, journal = {bioRxiv}, year = {2020}, month = {2020 Apr 10}, abstract = {

Coronavirus disease 2019 (COVID-19) outcomes vary from asymptomatic infection to death. This disparity may reflect different airway levels of the SARS-CoV-2 receptor, ACE2, and the spike protein activator, TMPRSS2. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci (eQTL) for both and , that vary in frequency across world populations. Importantly, we find is part of a mucus secretory network, highly upregulated by T2 inflammation through the action of interleukin-13, and that interferon response to respiratory viruses highly upregulates expression. Finally, we define airway responses to coronavirus infections in children, finding that these infections upregulate while also stimulating a more pronounced cytotoxic immune response relative to other respiratory viruses. Our results reveal mechanisms likely influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

}, doi = {10.1101/2020.04.09.034454}, author = {Sajuthi, Satria P and DeFord, Peter and Jackson, Nathan D and Montgomery, Michael T and Everman, Jamie L and Rios, Cydney L and Pruesse, Elmar and Nolin, James D and Plender, Elizabeth G and Wechsler, Michael E and Mak, Angel Cy and Eng, Celeste and Salazar, Sandra and Medina, Vivian and Wohlford, Eric M and Huntsman, Scott and Nickerson, Deborah A and Germer, Soren and Zody, Michael C and Abecasis, Goncalo and Kang, Hyun Min and Rice, Kenneth M and Kumar, Rajesh and Oh, Sam and Rodriguez-Santana, Jose and Burchard, Esteban G and Seibold, Max A} } @article {66, title = {Components of genetic associations across 2,138 phenotypes in the UK Biobank highlight adipocyte biology.}, journal = {Nat Commun}, volume = {10}, year = {2019}, month = {2019 Sep 06}, pages = {4064}, abstract = {

Population-based biobanks with genomic and dense phenotype data provide opportunities for generating effective therapeutic hypotheses and understanding the genomic role in disease predisposition. To characterize latent components of genetic associations, we apply truncated singular value decomposition (DeGAs) to matrices of summary statistics derived from genome-wide association analyses across 2,138 phenotypes measured in 337,199 White British individuals in the UK Biobank study. We systematically identify key components of genetic associations and the contributions of variants, genes, and phenotypes to each component. As an illustration of the utility of the approach to inform downstream experiments, we report putative loss of function variants, rs114285050 (GPR151) and rs150090666 (PDE3B), that substantially contribute to obesity-related traits and experimentally demonstrate the role of these genes in adipocyte biology. Our approach to dissect components of genetic associations across the human phenome will accelerate biomedical hypothesis generation by providing insights on previously unexplored latent structures.

}, issn = {2041-1723}, doi = {10.1038/s41467-019-11953-9}, author = {Tanigawa, Yosuke and Li, Jiehan and Justesen, Johanne M and Horn, Heiko and Aguirre, Matthew and DeBoever, Christopher and Chang, Chris and Narasimhan, Balasubramanian and Lage, Kasper and Hastie, Trevor and Park, Chong Y and Bejerano, Gill and Ingelsson, Erik and Rivas, Manuel A} } @article {50, title = {Efficient Variant Set Mixed Model Association Tests for Continuous and Binary Traits in Large-Scale Whole-Genome Sequencing Studies.}, journal = {Am J Hum Genet}, volume = {104}, year = {2019}, month = {2019 02 07}, pages = {260-274}, abstract = {

With advances in whole-genome sequencing (WGS) technology, more advanced statistical methods for testing genetic association with rare variants are being developed. Methods in which variants are grouped for analysis are also known as variant-set, gene-based, and aggregate unit tests. The burden test and sequence kernel association test (SKAT) are two widely used variant-set tests, which were originally developed for samples of unrelated individuals and later have been extended to family data with known pedigree structures. However, computationally efficient and powerful variant-set tests are needed to make analyses tractable in large-scale WGS studies with complex study samples. In this paper, we propose the variant-set mixed model association tests (SMMAT) for continuous and binary traits using the generalized linear mixed model framework. These tests can be applied to large-scale WGS studies involving samples with population structure and relatedness, such as in the National Heart, Lung, and Blood Institute{\textquoteright}s Trans-Omics for Precision Medicine (TOPMed) program. SMMATs share the same null model for different variant sets, and a virtue of this null model, which includes covariates only, is that it needs to be fit only once for all tests in each genome-wide analysis. Simulation studies show that all the proposed SMMATs correctly control type I error rates for both continuous and binary traits in the presence of population structure and relatedness. We also illustrate our tests in a real data example of analysis of plasma fibrinogen levels in the TOPMed program (n = 23,763), using the Analysis Commons, a cloud-based computing platform.

}, keywords = {Chromosomes, Human, Pair 4, Cloud Computing, Female, Fibrinogen, Genetic Association Studies, Genetics, Population, Humans, Male, Models, Genetic, National Heart, Lung, and Blood Institute (U.S.), Precision Medicine, Research Design, Time Factors, United States, Whole Genome Sequencing}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2018.12.012}, author = {Chen, Han and Huffman, Jennifer E and Brody, Jennifer A and Wang, Chaolong and Lee, Seunggeun and Li, Zilin and Gogarten, Stephanie M and Sofer, Tamar and Bielak, Lawrence F and Bis, Joshua C and Blangero, John and Bowler, Russell P and Cade, Brian E and Cho, Michael H and Correa, Adolfo and Curran, Joanne E and de Vries, Paul S and Glahn, David C and Guo, Xiuqing and Johnson, Andrew D and Kardia, Sharon and Kooperberg, Charles and Lewis, Joshua P and Liu, Xiaoming and Mathias, Rasika A and Mitchell, Braxton D and O{\textquoteright}Connell, Jeffrey R and Peyser, Patricia A and Post, Wendy S and Reiner, Alex P and Rich, Stephen S and Rotter, Jerome I and Silverman, Edwin K and Smith, Jennifer A and Vasan, Ramachandran S and Wilson, James G and Yanek, Lisa R and Redline, Susan and Smith, Nicholas L and Boerwinkle, Eric and Borecki, Ingrid B and Cupples, L Adrienne and Laurie, Cathy C and Morrison, Alanna C and Rice, Kenneth M and Lin, Xihong} } @article {42, title = {Error-prone bypass of DNA lesions during lagging-strand replication is a common source of germline and cancer mutations.}, journal = {Nat Genet}, volume = {51}, year = {2019}, month = {2019 01}, pages = {36-41}, abstract = {

Studies in experimental systems have identified a multitude of mutational mechanisms including DNA replication infidelity and DNA damage followed by inefficient repair or replicative bypass. However, the relative contributions of these mechanisms to human germline mutation remain unknown. Here, we show that error-prone damage bypass on the lagging strand plays a major role in human mutagenesis. Transcription-coupled DNA repair removes lesions on the transcribed strand; lesions on the non-transcribed strand are preferentially converted into mutations. In human polymorphism we detect a striking similarity between mutation types predominant on the non-transcribed strand and on the strand lagging during replication. Moreover, damage-induced mutations in cancers accumulate asymmetrically with respect to the direction of replication, suggesting that DNA lesions are resolved asymmetrically. We experimentally demonstrate that replication delay greatly attenuates the mutagenic effect of ultraviolet irradiation, confirming that replication converts DNA damage into mutations. We estimate that at least 10\% of human mutations arise due to DNA damage.

}, keywords = {Cells, Cultured, DNA, DNA Damage, DNA Repair, DNA Replication, Germ-Line Mutation, Humans, Mutagenesis, Neoplasms, Polymorphism, Single Nucleotide, Transcription, Genetic}, issn = {1546-1718}, doi = {10.1038/s41588-018-0285-7}, author = {Seplyarskiy, Vladimir B and Akkuratov, Evgeny E and Akkuratova, Natalia and Andrianova, Maria A and Nikolaev, Sergey I and Bazykin, Georgii A and Adameyko, Igor and Sunyaev, Shamil R} } @article {53, title = {Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies.}, journal = {Elife}, volume = {8}, year = {2019}, month = {2019 Mar 21}, abstract = {

Genetic predictions of height differ among human populations and these differences have been interpreted as evidence of polygenic adaptation. These differences were first detected using SNPs genome-wide significantly associated with height, and shown to grow stronger when large numbers of sub-significant SNPs were included, leading to excitement about the prospect of analyzing large fractions of the genome to detect polygenic adaptation for multiple traits. Previous studies of height have been based on SNP effect size measurements in the GIANT Consortium meta-analysis. Here we repeat the analyses in the UK Biobank, a much more homogeneously designed study. We show that polygenic adaptation signals based on large numbers of SNPs below genome-wide significance are extremely sensitive to biases due to uncorrected population stratification. More generally, our results imply that typical constructions of polygenic scores are sensitive to population stratification and that population-level differences should be interpreted with caution.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor{\textquoteright}s assessment is that all the issues have been addressed (see decision letter).

}, issn = {2050-084X}, doi = {10.7554/eLife.39702}, author = {Sohail, Mashaal and Maier, Robert M and Ganna, Andrea and Bloemendal, Alex and Martin, Alicia R and Turchin, Michael C and Chiang, Charleston Wk and Hirschhorn, Joel and Daly, Mark J and Patterson, Nick and Neale, Benjamin and Mathieson, Iain and Reich, David and Sunyaev, Shamil R} } @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 {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 05}, 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 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 {43, title = {Functional equivalence of genome sequencing analysis pipelines enables harmonized variant calling across human genetics projects.}, journal = {Nat Commun}, volume = {9}, year = {2018}, month = {2018 10 02}, pages = {4038}, abstract = {

Hundreds of thousands of human whole genome sequencing (WGS) datasets will be generated over the next few years. These data are more valuable in aggregate: joint analysis of genomes from many sources increases sample size and statistical power. A central challenge for joint analysis is that different WGS data processing pipelines cause substantial differences in variant calling in combined datasets, necessitating computationally expensive reprocessing. This approach is no longer tenable given the scale of current studies and data volumes. Here, we define WGS data processing standards that allow different groups to produce functionally equivalent (FE) results, yet still innovate on data processing pipelines. We present initial FE pipelines developed at five genome centers and show that they yield similar variant calling results and produce significantly less variability than sequencing replicates. This work alleviates a key technical bottleneck for genome aggregation and helps lay the foundation for community-wide human genetics studies.

}, keywords = {Genome, Human, Human Genetics, Humans, Whole Genome Sequencing}, issn = {2041-1723}, doi = {10.1038/s41467-018-06159-4}, author = {Regier, Allison A and Farjoun, Yossi and Larson, David E and Krasheninina, Olga and Kang, Hyun Min and Howrigan, Daniel P and Chen, Bo-Juen and Kher, Manisha and Banks, Eric and Ames, Darren C and English, Adam C and Li, Heng and Xing, Jinchuan and Zhang, Yeting and Matise, Tara and Abecasis, Gon{\c c}alo R and Salerno, Will and Zody, Michael C and Neale, Benjamin M and Hall, Ira M} } @article {61, title = {Genetic diversity in populations across Latin America: implications for population and medical genetic studies.}, journal = {Curr Opin Genet Dev}, volume = {53}, year = {2018}, month = {2018 12}, pages = {98-104}, abstract = {

Hispanic/Latino (H/L) populations, although linked by culture and aspects of shared history, reflect the complexity of history and migration influencing the Americas. The original settlement by indigenous Americans, followed by postcolonial admixture from multiple continents, has yielded localized genetic patterns. In addition, numerous H/L populations appear to have signatures of pre-colonization and post-colonization bottlenecks, indicating that tens of millions of H/Ls may harbor signatures of founder effects today. Based on both population and medical genetic findings we highlight the extreme differentiation across the Americas, providing evidence for why H/Ls should not be considered a single population in modern human genetics. We highlight the need for additional sampling of understudied H/L groups, and ramifications of these findings for genomic medicine in one-tenth of the world{\textquoteright}s population.

}, keywords = {African Continental Ancestry Group, European Continental Ancestry Group, Genetic Variation, Genetics, Medical, Genetics, Population, Genome, Human, Genomics, Humans, Latin America}, issn = {1879-0380}, doi = {10.1016/j.gde.2018.07.006}, author = {Belbin, Gillian M and Nieves-Col{\'o}n, Maria A and Kenny, Eimear E and Moreno-Estrada, Andres and Gignoux, Christopher R} } @article {34, title = {A synthetic-diploid benchmark for accurate variant-calling evaluation.}, journal = {Nat Methods}, volume = {15}, year = {2018}, month = {2018 08}, pages = {595-597}, abstract = {

Existing benchmark datasets for use in evaluating variant-calling accuracy are constructed from a consensus of known short-variant callers, and they are thus biased toward easy regions that are accessible by these algorithms. We derived a new benchmark dataset from the de novo PacBio assemblies of two fully homozygous human cell lines, which provides a relatively more accurate and less biased estimate of small-variant-calling error rates in a realistic context.

}, keywords = {Algorithms, Benchmarking, Cell Line, Tumor, Databases, Genetic, Diploidy, Female, Genetic Variation, Genome, Human, Homozygote, Humans, Hydatidiform Mole, Pregnancy, Synthetic Biology, Uterine Neoplasms, Whole Genome Sequencing}, issn = {1548-7105}, doi = {10.1038/s41592-018-0054-7}, author = {Li, Heng and Bloom, Jonathan M and Farjoun, Yossi and Fleharty, Mark and Gauthier, Laura and Neale, Benjamin and MacArthur, Daniel} } @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 {18, title = {Genetic effects on gene expression across human tissues.}, journal = {Nature}, volume = {550}, year = {2017}, month = {2017 Oct 11}, pages = {204-213}, abstract = {

Characterization of the molecular function of the human genome and its variation across individuals is essential for identifying the cellular mechanisms that underlie human genetic traits and diseases. The Genotype-Tissue Expression (GTEx) project aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body, many of which are not easily accessible. Here we describe genetic effects on gene expression levels across 44 human tissues. We find that local genetic variation affects gene expression levels for the majority of genes, and we further identify inter-chromosomal genetic effects for 93 genes and 112 loci. On the basis of the identified genetic effects, we characterize patterns of tissue specificity, compare local and distal effects, and evaluate the functional properties of the genetic effects. We also demonstrate that multi-tissue, multi-individual data can be used to identify genes and pathways affected by human disease-associated variation, enabling a mechanistic interpretation of gene regulation and the genetic basis of disease.

}, keywords = {Alleles, Chromosomes, Human, Disease, Female, Gene Expression Profiling, Gene Expression Regulation, Genetic Variation, Genome, Human, Genotype, Humans, Male, Organ Specificity, Quantitative Trait Loci}, issn = {1476-4687}, doi = {10.1038/nature24277}, author = {Battle, Alexis and Brown, Christopher D and Engelhardt, Barbara E and Montgomery, Stephen B} } @article {32, title = {Genetic identification of a common collagen disease in puerto ricans via identity-by-descent mapping in a health system.}, journal = {Elife}, volume = {6}, year = {2017}, month = {2017 09 12}, abstract = {

Achieving confidence in the causality of a disease locus is a complex task that often requires supporting data from both statistical genetics and clinical genomics. Here we describe a combined approach to identify and characterize a genetic disorder that leverages distantly related patients in a health system and population-scale mapping. We utilize genomic data to uncover components of distant pedigrees, in the absence of recorded pedigree information, in the multi-ethnic Bio biobank in New York City. By linking to medical records, we discover a locus associated with both elevated genetic relatedness and extreme short stature. We link the gene, , with a little-known genetic disease, previously thought to be rare and recessive. We demonstrate that disease manifests in both heterozygotes and homozygotes, indicating a common collagen disorder impacting up to 2\% of individuals of Puerto Rican ancestry, leading to a better understanding of the continuum of complex and Mendelian disease.

}, keywords = {Adolescent, Adult, Aged, Child, Collagen Diseases, Female, Fibrillar Collagens, Genotype, Heterozygote, Hispanic Americans, Homozygote, Humans, Male, Middle Aged, Molecular Epidemiology, Multigene Family, Musculoskeletal Diseases, New York City, Pedigree, Whole Genome Sequencing, Young Adult}, issn = {2050-084X}, doi = {10.7554/eLife.25060}, author = {Belbin, Gillian Morven and Odgis, Jacqueline and Sorokin, Elena P and Yee, Muh-Ching and Kohli, Sumita and Glicksberg, Benjamin S and Gignoux, Christopher R and Wojcik, Genevieve L and Van Vleck, Tielman and Jeff, Janina M and Linderman, Michael and Schurmann, Claudia and Ruderfer, Douglas and Cai, Xiaoqiang and Merkelson, Amanda and Justice, Anne E and Young, Kristin L and Graff, Misa and North, Kari E and Peters, Ulrike and James, Regina and Hindorff, Lucia and Kornreich, Ruth and Edelmann, Lisa and Gottesman, Omri and Stahl, Eli Ea and Cho, Judy H and Loos, Ruth Jf and Bottinger, Erwin P and Nadkarni, Girish N and Abul-Husn, Noura S and Kenny, Eimear E} } @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 {30, title = {The impact of rare variation on gene expression across tissues.}, journal = {Nature}, volume = {550}, year = {2017}, month = {2017 10 11}, pages = {239-243}, abstract = {

Rare genetic variants are abundant in humans and are expected to contribute to individual disease risk. While genetic association studies have successfully identified common genetic variants associated with susceptibility, these studies are not practical for identifying rare variants. Efforts to distinguish pathogenic variants from benign rare variants have leveraged the genetic code to identify deleterious protein-coding alleles, but no analogous code exists for non-coding variants. Therefore, ascertaining which rare variants have phenotypic effects remains a major challenge. Rare non-coding variants have been associated with extreme gene expression in studies using single tissues, but their effects across tissues are unknown. Here we identify gene expression outliers, or individuals showing extreme expression levels for a particular gene, across 44 human tissues by using combined analyses of whole genomes and multi-tissue RNA-sequencing data from the Genotype-Tissue Expression (GTEx) project v6p release. We find that 58\% of underexpression and 28\% of overexpression outliers have nearby conserved rare variants compared to 8\% of non-outliers. Additionally, we developed RIVER (RNA-informed variant effect on regulation), a Bayesian statistical model that incorporates expression data to predict a regulatory effect for rare variants with higher accuracy than models using genomic annotations alone. Overall, we demonstrate that rare variants contribute to large gene expression changes across tissues and provide an integrative method for interpretation of rare variants in individual genomes.

}, keywords = {Bayes Theorem, Female, Gene Expression Profiling, Genetic Variation, Genome, Human, Genomics, Genotype, Humans, Male, Models, Genetic, Organ Specificity, Sequence Analysis, RNA}, issn = {1476-4687}, doi = {10.1038/nature24267}, author = {Li, Xin and Kim, Yungil and Tsang, Emily K and Davis, Joe R and Damani, Farhan N and Chiang, Colby and Hess, Gaelen T and Zappala, Zachary and Strober, Benjamin J and Scott, Alexandra J and Li, Amy and Ganna, Andrea and Bassik, Michael C and Merker, Jason D and Hall, Ira M and Battle, Alexis and Montgomery, Stephen B} } @article {38, title = {Linkage disequilibrium-dependent architecture of human complex traits shows action of negative selection.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 Oct}, pages = {1421-1427}, abstract = {

Recent work has hinted at the linkage disequilibrium (LD)-dependent architecture of human complex traits, where SNPs with low levels of LD (LLD) have larger per-SNP heritability. Here we analyzed summary statistics from 56 complex traits (average N = 101,401) by extending stratified LD score regression to continuous annotations. We determined that SNPs with low LLD have significantly larger per-SNP heritability and that roughly half of this effect can be explained by functional annotations negatively correlated with LLD, such as DNase I hypersensitivity sites (DHSs). The remaining signal is largely driven by our finding that more recent common variants tend to have lower LLD and to explain more heritability (P = 2.38 {\texttimes} 10); the youngest 20\% of common SNPs explain 3.9 times more heritability than the oldest 20\%, consistent with the action of negative selection. We also inferred jointly significant effects of other LD-related annotations and confirmed via forward simulations that they jointly predict deleterious effects.

}, keywords = {Alleles, Chi-Square Distribution, Datasets as Topic, Genetic Fitness, Genetic Variation, Humans, Linkage Disequilibrium, Models, Genetic, Molecular Sequence Annotation, Multifactorial Inheritance, Polymorphism, Single Nucleotide, Selection, Genetic}, issn = {1546-1718}, doi = {10.1038/ng.3954}, author = {Gazal, Steven and Finucane, Hilary K and Furlotte, Nicholas A and Loh, Po-Ru and Palamara, Pier Francesco and Liu, Xuanyao and Schoech, Armin and Bulik-Sullivan, Brendan and Neale, Benjamin M and Gusev, Alexander and Price, Alkes L} } @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} }