@article {117, title = {A deoxyribonuclease 1-like 3 genetic variant associates with asthma exacerbations.}, journal = {J Allergy Clin Immunol}, volume = {147}, year = {2021}, month = {2021 Mar}, pages = {1095-1097.e10}, issn = {1097-6825}, doi = {10.1016/j.jaci.2020.09.027}, author = {Herrera-Luis, Esther and Lorenzo-Diaz, Fabian and Samedy-Bates, Lesly-Anne and Eng, Celeste and Villar, Jes{\'u}s and Rodriguez-Santana, Jose R and Burchard, Esteban G and Pino-Yanes, Maria} } @article {114, title = {Genetics of 35 blood and urine biomarkers in the UK Biobank.}, journal = {Nat Genet}, volume = {53}, year = {2021}, month = {2021 02}, pages = {185-194}, abstract = {

Clinical laboratory tests are a critical component of the continuum of care. We evaluate the genetic basis of 35 blood and urine laboratory measurements in the UK Biobank (n = 363,228 individuals). We identify 1,857 loci associated with at least one trait, containing 3,374 fine-mapped associations and additional sets of large-effect (>0.1 s.d.) protein-altering, human leukocyte antigen (HLA) and copy number variant (CNV) associations. Through Mendelian randomization (MR) analysis, we discover 51 causal relationships, including previously known agonistic effects of urate on gout and cystatin C on stroke. Finally, we develop polygenic risk scores (PRSs) for each biomarker and build {\textquoteright}multi-PRS{\textquoteright} models for diseases using 35 PRSs simultaneously, which improved chronic kidney disease, type 2 diabetes, gout and alcoholic cirrhosis genetic risk stratification in an independent dataset (FinnGen; n = 135,500) relative to single-disease PRSs. Together, our results delineate the genetic basis of biomarkers and their causal influences on diseases and improve genetic risk stratification for common diseases.

}, keywords = {Biological Specimen Banks, Biomarkers, Cardiovascular Diseases, Diabetes Mellitus, Type 2, DNA Copy Number Variations, Genetic Pleiotropy, HLA Antigens, Humans, Linkage Disequilibrium, Liver-Specific Organic Anion Transporter 1, Mendelian Randomization Analysis, Polymorphism, Single Nucleotide, Proteins, Renal Insufficiency, Chronic, Serine Endopeptidases, United Kingdom}, issn = {1546-1718}, doi = {10.1038/s41588-020-00757-z}, author = {Sinnott-Armstrong, Nasa and Tanigawa, Yosuke and Amar, David and Mars, Nina and Benner, Christian and Aguirre, Matthew and Venkataraman, Guhan Ram and Wainberg, Michael and Ollila, Hanna M and Kiiskinen, Tuomo and Havulinna, Aki S and Pirruccello, James P and Qian, Junyang and Shcherbina, Anna and Rodriguez, Fatima and Assimes, Themistocles L and Agarwala, Vineeta and Tibshirani, Robert and Hastie, Trevor and Ripatti, Samuli and Pritchard, Jonathan K and Daly, Mark J and Rivas, Manuel A} } @article {98, title = {Genome-wide association study reveals a novel locus for asthma with severe exacerbations in diverse populations.}, journal = {Pediatr Allergy Immunol}, volume = {32}, year = {2021}, month = {2021 Jan}, pages = {106-115}, abstract = {

BACKGROUND: Severe asthma exacerbations are a major cause of asthma morbidity and increased healthcare costs. Several studies have shown racial and ethnic differences in asthma exacerbation rates. We aimed to identify genetic variants associated with severe exacerbations in two high-risk populations for asthma.

METHODS: A genome-wide association study of asthma in children and youth with severe exacerbations was performed in 1283 exacerbators and 2027 controls without asthma of Latino ancestry. Independent suggestive variants (P <= 5 {\texttimes} 10 ) were selected for replication in 448 African Americans exacerbators and 595 controls. Case-only analyses were performed comparing the exacerbators with additional 898 Latinos and 524 African Americans asthma patients without exacerbations, while adjusting by treatment category as a proxy of asthma severity. We analyzed the functionality of associated variants with in silico methods and by correlating genotypes with methylation levels in whole blood in a subset of 473 Latinos.

RESULTS: We identified two genome-wide significant associations for susceptibility to asthma with severe exacerbations, including a novel locus located at chromosome 2p21 (rs4952375, odds ratio = 1.39, P = 3.8 {\texttimes} 10 ), which was also associated with asthma exacerbations in a case-only analysis (odds ratio = 1.25, P = 1.95 {\texttimes} 10 ). This polymorphism is an expression quantitative trait locus of the long intergenic non-protein coding RNA 1913 (LINC01913) in lung tissues (P = 1.3 {\texttimes} 10 ) and influences methylation levels of the protein kinase domain-containing cytoplasmic (PKDCC) gene in whole-blood cells (P = 9.8 {\texttimes} 10 ).

CONCLUSION: We identified a novel susceptibility locus for severe asthma exacerbations in Hispanic/Latino and African American youths with functional effects in gene expression and methylation status of neighboring genes.

}, issn = {1399-3038}, doi = {10.1111/pai.13337}, author = {Herrera-Luis, Esther and Espuela-Ortiz, Antonio and Lorenzo-Diaz, Fabian and Keys, Kevin L and Mak, Angel C Y and Eng, Celeste and Huntsman, Scott and Villar, Jes{\'u}s and Rodriguez-Santana, Jose R and Burchard, Esteban G and Pino-Yanes, Maria} } @article {142, title = {Identification of rare and common regulatory variants in pluripotent cells using population-scale transcriptomics.}, journal = {Nat Genet}, volume = {53}, year = {2021}, month = {2021 03}, pages = {313-321}, abstract = {

Induced pluripotent stem cells (iPSCs) are an established cellular system to study the impact of genetic variants in derived cell types and developmental contexts. However, in their pluripotent state, the disease impact of genetic variants is less well known. Here, we integrate data from 1,367 human iPSC lines to comprehensively map common and rare regulatory variants in human pluripotent cells. Using this population-scale resource, we report hundreds of new colocalization events for human traits specific to iPSCs, and find increased power to identify rare regulatory variants compared with somatic tissues. Finally, we demonstrate how iPSCs enable the identification of causal genes for rare diseases.

}, keywords = {Bardet-Biedl Syndrome, Calcium Channels, Cell Line, Cerebellar Ataxia, DNA Methylation, Gene Expression, Genetic Variation, Humans, Induced Pluripotent Stem Cells, Polymorphism, Single Nucleotide, Proteins, Quantitative Trait Loci, Rare Diseases, Regulatory Sequences, Nucleic Acid, Sequence Analysis, RNA, Whole Genome Sequencing}, issn = {1546-1718}, doi = {10.1038/s41588-021-00800-7}, author = {Bonder, Marc Jan and Smail, Craig and Gloudemans, Michael J and Fr{\'e}sard, Laure and Jakubosky, David and D{\textquoteright}Antonio, Matteo and Li, Xin and Ferraro, Nicole M and Carcamo-Orive, Ivan and Mirauta, Bogdan and Seaton, Daniel D and Cai, Na and Vakili, Dara and Horta, Danilo and Zhao, Chunli and Zastrow, Diane B and Bonner, Devon E and Wheeler, Matthew T and Kilpinen, Helena and Knowles, Joshua W and Smith, Erin N and Frazer, Kelly A and Montgomery, Stephen B and Stegle, Oliver} } @article {135, title = {Polygenic risk modeling with latent trait-related genetic components.}, journal = {Eur J Hum Genet}, year = {2021}, month = {2021 Feb 08}, abstract = {

Polygenic risk models have led to significant advances in understanding complex diseases and their clinical presentation. While polygenic risk scores (PRS) can effectively predict outcomes, they do not generally account for disease subtypes or pathways which underlie within-trait diversity. Here, we introduce a latent factor model of genetic risk based on components from Decomposition of Genetic Associations (DeGAs), which we call the DeGAs polygenic risk score (dPRS). We compute DeGAs using genetic associations for 977 traits and find that dPRS performs comparably to standard PRS while offering greater interpretability. We show how to decompose an individual{\textquoteright}s genetic risk for a trait across DeGAs components, with examples for body mass index (BMI) and myocardial infarction (heart attack) in 337,151 white British individuals in the UK Biobank, with replication in a further set of 25,486 non-British white individuals. We find that BMI polygenic risk factorizes into components related to fat-free mass, fat mass, and overall health indicators like physical activity. Most individuals with high dPRS for BMI have strong contributions from both a fat-mass component and a fat-free mass component, whereas a few "outlier" individuals have strong contributions from only one of the two components. Overall, our method enables fine-scale interpretation of the drivers of genetic risk for complex traits.

}, issn = {1476-5438}, doi = {10.1038/s41431-021-00813-0}, author = {Aguirre, Matthew and Tanigawa, Yosuke and Venkataraman, Guhan Ram and Tibshirani, Rob and Hastie, Trevor and Rivas, Manuel A} } @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 {112, title = {Whole exome sequencing analyses reveal gene-microbiota interactions in the context of IBD.}, journal = {Gut}, volume = {70}, year = {2021}, month = {2021 Feb}, pages = {285-296}, abstract = {

OBJECTIVE: Both the gut microbiome and host genetics are known to play significant roles in the pathogenesis of IBD. However, the interaction between these two factors and its implications in the aetiology of IBD remain underexplored. Here, we report on the influence of host genetics on the gut microbiome in IBD.

DESIGN: To evaluate the impact of host genetics on the gut microbiota of patients with IBD, we combined whole exome sequencing of the host genome and whole genome shotgun sequencing of 1464 faecal samples from 525 patients with IBD and 939 population-based controls. We followed a four-step analysis: (1) exome-wide microbial quantitative trait loci (mbQTL) analyses, (2) a targeted approach focusing on IBD-associated genomic regions and protein truncating variants (PTVs, minor allele frequency (MAF) >5\%), (3) gene-based burden tests on PTVs with MAF <5\% and exome copy number variations (CNVs) with site frequency <1\%, (4) joint analysis of both cohorts to identify the interactions between disease and host genetics.

RESULTS: We identified 12 mbQTLs, including variants in the IBD-associated genes , , and . For example, the decrease of the pathway acetyl-coenzyme A biosynthesis, which is involved in short chain fatty acids production, was associated with variants in the gene (false discovery rate <0.05). Changes in functional pathways involved in the metabolic potential were also observed in participants carrying rare PTVs or CNVs in , and genes. These genes are known for their function in the immune system. Moreover, interaction analyses confirmed previously known IBD disease-specific mbQTLs in .

CONCLUSION: This study highlights that both common and rare genetic variants affecting the immune system are key factors in shaping the gut microbiota in the context of IBD and pinpoints towards potential mechanisms for disease treatment.

}, issn = {1468-3288}, doi = {10.1136/gutjnl-2019-319706}, author = {Hu, Shixian and Vich Vila, Arnau and Gacesa, Ranko and Collij, Valerie and Stevens, Christine and Fu, Jack M and Wong, Isaac and Talkowski, Michael E and Rivas, Manuel A and Imhann, Floris and Bolte, Laura and van Dullemen, Hendrik and Dijkstra, Gerard and Visschedijk, Marijn C and Festen, Eleonora A and Xavier, Ramnik J and Fu, Jingyuan and Daly, Mark J and Wijmenga, Cisca and Zhernakova, Alexandra and Kurilshikov, Alexander and Weersma, Rinse K} } @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 {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 {104, title = {FasTag: Automatic text classification of unstructured medical narratives.}, journal = {PLoS One}, volume = {15}, year = {2020}, month = {2020}, pages = {e0234647}, abstract = {

Unstructured clinical narratives are continuously being recorded as part of delivery of care in electronic health records, and dedicated tagging staff spend considerable effort manually assigning clinical codes for billing purposes. Despite these efforts, however, label availability and accuracy are both suboptimal. In this retrospective study, we aimed to automate the assignment of top-level International Classification of Diseases version 9 (ICD-9) codes to clinical records from human and veterinary data stores using minimal manual labor and feature curation. Automating top-level annotations could in turn enable rapid cohort identification, especially in a veterinary setting. To this end, we trained long short-term memory (LSTM) recurrent neural networks (RNNs) on 52,722 human and 89,591 veterinary records. We investigated the accuracy of both separate-domain and combined-domain models and probed model portability. We established relevant baseline classification performances by training Decision Trees (DT) and Random Forests (RF). We also investigated whether transforming the data using MetaMap Lite, a clinical natural language processing tool, affected classification performance. We showed that the LSTM-RNNs accurately classify veterinary and human text narratives into top-level categories with an average weighted macro F1 score of 0.74 and 0.68 respectively. In the "neoplasia" category, the model trained on veterinary data had a high validation accuracy in veterinary data and moderate accuracy in human data, with F1 scores of 0.91 and 0.70 respectively. Our LSTM method scored slightly higher than that of the DT and RF models. The use of LSTM-RNN models represents a scalable structure that could prove useful in cohort identification for comparative oncology studies. Digitization of human and veterinary health information will continue to be a reality, particularly in the form of unstructured narratives. Our approach is a step forward for these two domains to learn from and inform one another.

}, keywords = {Animals, Automation, Data Mining, Databases as Topic, Humans, Narrative Medicine, Reproducibility of Results, Software, Species Specificity}, issn = {1932-6203}, doi = {10.1371/journal.pone.0234647}, author = {Venkataraman, Guhan Ram and Pineda, Arturo Lopez and Bear Don{\textquoteright}t Walk Iv, Oliver J and Zehnder, Ashley M and Ayyar, Sandeep and Page, Rodney L and Bustamante, Carlos D and Rivas, Manuel A} } @article {107, title = {Genome-Wide Analysis Reveals Mucociliary Remodeling of the Nasal Airway Epithelium Induced by Urban PM.}, journal = {Am J Respir Cell Mol Biol}, volume = {63}, year = {2020}, month = {2020 08}, pages = {172-184}, abstract = {

Air pollution particulate matter <2.5 μm (PM) exposure is associated with poor respiratory outcomes. Mechanisms underlying PM-induced lung pathobiology are poorly understood but likely involve cellular and molecular changes to the airway epithelium. We extracted and chemically characterized the organic and water-soluble components of air pollution PM samples, then determined the whole transcriptome response of human nasal mucociliary airway epithelial cultures to a dose series of PM extracts. We found that PM organic extract (OE), but not water-soluble extract, elicited a potent, dose-dependent transcriptomic response from the mucociliary epithelium. Exposure to a moderate OE dose modified the expression of 424 genes, including activation of aryl hydrocarbon receptor signaling and an IL-1 inflammatory program. We generated an OE-response gene network defined by eight functional enrichment groups, which exhibited high connectivity through , , and . This OE exposure also robustly activated a mucus secretory expression program (>100 genes), which included transcriptional drivers of mucus metaplasia ( and ). Exposure to a higher OE dose modified the expression of 1,240 genes and further exacerbated expression responses observed at the moderate dose, including the mucus secretory program. Moreover, the higher OE dose significantly increased the / gel-forming mucin expression ratio and strongly downregulated ciliated cell expression programs, including key ciliating cell transcription factors (e.g., and ). Chronic OE stimulation induced mucus metaplasia-like remodeling characterized by increases in MUC5AC secretory cells and MUC5AC mucus secretions. This epithelial remodeling may underlie poor respiratory outcomes associated with high PM exposure.

}, issn = {1535-4989}, doi = {10.1165/rcmb.2019-0454OC}, author = {Montgomery, Michael T and Sajuthi, Satria P and Cho, Seung-Hyun and Everman, Jamie L and Rios, Cydney L and Goldfarbmuren, Katherine C and Jackson, Nathan D and Saef, Benjamin and Cromie, Meghan and Eng, Celeste and Medina, Vivian and Elhawary, Jennifer R and Oh, Sam S and Rodriguez-Santana, Jose and Vladar, Eszter K and Burchard, Esteban G and Seibold, Max A} } @article {81, title = {Identification of cancer driver genes based on nucleotide context.}, journal = {Nat Genet}, volume = {52}, year = {2020}, month = {2020 02}, pages = {208-218}, abstract = {

Cancer genomes contain large numbers of somatic mutations but few of these mutations drive tumor development. Current approaches either identify driver genes on the basis of mutational recurrence or approximate the functional consequences of nonsynonymous mutations by using bioinformatic scores. Passenger mutations are enriched in characteristic nucleotide contexts, whereas driver mutations occur in functional positions, which are not necessarily surrounded by a particular nucleotide context. We observed that mutations in contexts that deviate from the characteristic contexts around passenger mutations provide a signal in favor of driver genes. We therefore developed a method that combines this feature with the signals traditionally used for driver-gene identification. We applied our method to whole-exome sequencing data from 11,873 tumor-normal pairs and identified 460 driver genes that clustered into 21 cancer-related pathways. Our study provides a resource of driver genes across 28 tumor types with additional driver genes identified according to mutations in unusual nucleotide contexts.

}, keywords = {Cluster Analysis, Computational Biology, Humans, Mutation, Neoplasms, Nucleotides, Proteins, Whole Exome Sequencing}, issn = {1546-1718}, doi = {10.1038/s41588-019-0572-y}, author = {Dietlein, Felix and Weghorn, Donate and Taylor-Weiner, Amaro and Richters, Andr{\'e} and Reardon, Brendan and Liu, David and Lander, Eric S and Van Allen, Eliezer M and Sunyaev, Shamil 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 {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 {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 {29, title = {Phenotype risk scores identify patients with unrecognized Mendelian disease patterns.}, journal = {Science}, volume = {359}, year = {2018}, month = {2018 03 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 {21, title = {Covariate selection for association screening in multiphenotype genetic studies.}, journal = {Nat Genet}, volume = {49}, year = {2017}, month = {2017 Dec}, pages = {1789-1795}, abstract = {

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

}, keywords = {Algorithms, Genetic Association Studies, Genetic Variation, Genome-Wide Association Study, Genotype, Humans, Models, Genetic, Multivariate Analysis, Phenotype, Reproducibility of Results, Sample Size}, issn = {1546-1718}, doi = {10.1038/ng.3975}, author = {Aschard, Hugues and Guillemot, Vincent and Vilhjalmsson, Bjarni and Patel, Chirag J and Skurnik, David and Ye, Chun J and Wolpin, Brian and Kraft, Peter and Zaitlen, Noah} } @article {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 {13, title = {Negative selection in humans and fruit flies involves synergistic epistasis.}, journal = {Science}, volume = {356}, year = {2017}, month = {2017 May 05}, pages = {539-542}, abstract = {

Negative selection against deleterious alleles produced by mutation influences within-population variation as the most pervasive form of natural selection. However, it is not known whether deleterious alleles affect fitness independently, so that cumulative fitness loss depends exponentially on the number of deleterious alleles, or synergistically, so that each additional deleterious allele results in a larger decrease in relative fitness. Negative selection with synergistic epistasis should produce negative linkage disequilibrium between deleterious alleles and, therefore, an underdispersed distribution of the number of deleterious alleles in the genome. Indeed, we detected underdispersion of the number of rare loss-of-function alleles in eight independent data sets from human and fly populations. Thus, selection against rare protein-disrupting alleles is characterized by synergistic epistasis, which may explain how human and fly populations persist despite high genomic mutation rates.

}, keywords = {Alleles, Animals, Drosophila melanogaster, Epistasis, Genetic, Genetic Fitness, Genome, Human, Genome, Insect, Humans, Linkage Disequilibrium, Mutation Rate, Mutation, Missense, Selection, Genetic}, issn = {1095-9203}, doi = {10.1126/science.aah5238}, author = {Sohail, Mashaal and Vakhrusheva, Olga A and Sul, Jae Hoon and Pulit, Sara L and Francioli, Laurent C and van den Berg, Leonard H and Veldink, Jan H and de Bakker, Paul I W and Bazykin, Georgii A and Kondrashov, Alexey S and Sunyaev, Shamil R} } @article {16, title = {Testing for the indirect effect under the null for genome-wide mediation analyses.}, journal = {Genet Epidemiol}, volume = {41}, year = {2017}, month = {2017 Dec}, pages = {824-833}, abstract = {

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

}, keywords = {Basic Helix-Loop-Helix Transcription Factors, DNA Methylation, Epigenomics, Genome-Wide Association Study, Humans, Lung Neoplasms, Models, Genetic, Repressor Proteins}, issn = {1098-2272}, doi = {10.1002/gepi.22084}, author = {Barfield, Richard and Shen, Jincheng and Just, Allan C and Vokonas, Pantel S and Schwartz, Joel and Baccarelli, Andrea A and VanderWeele, Tyler J and Lin, Xihong} }