@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 {74, title = {Identifying causal variants and genes using functional genomics in specialized cell types and contexts.}, journal = {Hum Genet}, volume = {139}, year = {2020}, month = {2020 Jan}, pages = {95-102}, abstract = {

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

}, keywords = {Cell Lineage, Genes, Genetic Predisposition to Disease, Genome, Human, Genome-Wide Association Study, Genomics, Humans, Polymorphism, Single Nucleotide, Quantitative Trait Loci}, issn = {1432-1203}, doi = {10.1007/s00439-019-02044-2}, author = {Liu, Boxiang and Montgomery, Stephen B} } @article {88, title = {Whole-Genome and RNA Sequencing Reveal Variation and Transcriptomic Coordination in the Developing Human Prefrontal Cortex.}, journal = {Cell Rep}, volume = {31}, year = {2020}, month = {2020 04 07}, pages = {107489}, abstract = {

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

}, keywords = {Base Sequence, Brain, Computational Biology, Databases, Genetic, Genetic Predisposition to Disease, Genetic Variation, Genome-Wide Association Study, Genomics, Humans, Phenotype, Polymorphism, Single Nucleotide, Prefrontal Cortex, Quantitative Trait Loci, Sequence Analysis, RNA, Transcriptome, Whole Exome Sequencing, Whole Genome Sequencing}, issn = {2211-1247}, doi = {10.1016/j.celrep.2020.03.053}, author = {Werling, Donna M and Pochareddy, Sirisha and Choi, Jinmyung and An, Joon-Yong and Sheppard, Brooke and Peng, Minshi and Li, Zhen and Dastmalchi, Claudia and Santpere, Gabriel and Sousa, Andr{\'e} M M and Tebbenkamp, Andrew T N and Kaur, Navjot and Gulden, Forrest O and Breen, Michael S and Liang, Lindsay and Gilson, Michael C and Zhao, Xuefang and Dong, Shan and Klei, Lambertus and Cicek, A Ercument and Buxbaum, Joseph D and Adle-Biassette, Homa and Thomas, Jean-Leon and Aldinger, Kimberly A and O{\textquoteright}Day, Diana R and Glass, Ian A and Zaitlen, Noah A and Talkowski, Michael E and Roeder, Kathryn and State, Matthew W and Devlin, Bernie and Sanders, Stephan J and Sestan, Nenad} } @article {64, title = {A Bayesian framework that integrates multi-omics data and gene networks predicts risk genes from schizophrenia GWAS data.}, journal = {Nat Neurosci}, volume = {22}, year = {2019}, month = {2019 05}, pages = {691-699}, abstract = {

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

}, keywords = {Animals, Bayes Theorem, Disease Models, Animal, Gene Regulatory Networks, Genetic Predisposition to Disease, Genome-Wide Association Study, Genomics, Humans, Mice, Risk Factors, Schizophrenia}, issn = {1546-1726}, doi = {10.1038/s41593-019-0382-7}, author = {Wang, Quan and Chen, Rui and Cheng, Feixiong and Wei, Qiang and Ji, Ying and Yang, Hai and Zhong, Xue and Tao, Ran and Wen, Zhexing and Sutcliffe, James S and Liu, Chunyu and Cook, Edwin H and Cox, Nancy J and Li, Bingshan} } @article {51, title = {A multi-task convolutional deep neural network for variant calling in single molecule sequencing.}, journal = {Nat Commun}, volume = {10}, year = {2019}, month = {2019 Mar 01}, pages = {998}, abstract = {

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

}, keywords = {Base Sequence, Computational Biology, DNA Mutational Analysis, Genome, Human, Genome-Wide Association Study, Genomics, Genotype, Genotyping Techniques, Humans, INDEL Mutation, Nanopores, Neural Networks, Computer, Polymorphism, Single Nucleotide, Sequence Analysis, DNA, Software}, issn = {2041-1723}, doi = {10.1038/s41467-019-09025-z}, author = {Luo, Ruibang and Sedlazeck, Fritz J and Lam, Tak-Wah and Schatz, Michael C} } @article {31, title = {Functional genomics of stromal cells in chronic inflammatory diseases.}, journal = {Curr Opin Rheumatol}, volume = {30}, year = {2018}, month = {2018 Jan}, pages = {65-71}, abstract = {

PURPOSE OF REVIEW: Stroma is a broad term referring to the connective tissue matrix in which other cells reside. It is composed of diverse cell types with functions such as extracellular matrix maintenance, blood and lymph vessel development, and effector cell recruitment. The tissue microenvironment is determined by the molecular characteristics and relative abundances of different stromal cells such as fibroblasts, endothelial cells, pericytes, and mesenchymal precursor cells. Stromal cell heterogeneity is explained by embryonic developmental lineage, stages of differentiation to other cell types, and activation states. Interaction between immune and stromal cell types is critical to wound healing, cancer, and a wide range of inflammatory diseases. Here, we review recent studies of inflammatory diseases that use functional genomics and single-cell technologies to identify and characterize stromal cell types associated with pathogenesis.

RECENT FINDINGS: High dimensional strategies using mRNA sequencing, mass cytometry, and fluorescence activated cell-sorting with fresh primary tissue samples are producing detailed views of what is happening in diseased tissue in rheumatoid arthritis, inflammatory bowel disease, and cancer. Fibroblasts positive for CD90 (Thy-1) are enriched in the synovium of rheumatoid arthritis patients. Single-cell RNA-seq studies will lead to more discoveries about the stroma in the near future.

SUMMARY: Stromal cells form the microenvironment of inflamed and diseased tissues. Functional genomics is producing an increasingly detailed view of subsets of stromal cells with pathogenic functions in rheumatic diseases and cancer. Future genomics studies will discover disease mechanisms by perturbing molecular pathways with chemokines and therapies known to affect patient outcomes. Functional genomics studies with large sample sizes of patient tissues will identify patient subsets with different disease phenotypes or treatment responses.

}, keywords = {Arthritis, Rheumatoid, Cell Differentiation, Chemokines, Endothelial Cells, Fibroblasts, Flow Cytometry, Genomics, Humans, Inflammatory Bowel Diseases, Mesenchymal Stem Cells, Pericytes, Rheumatic Diseases, RNA, Messenger, Sequence Analysis, RNA, Single-Cell Analysis, Stromal Cells, Synovial Membrane, Thy-1 Antigens}, issn = {1531-6963}, doi = {10.1097/BOR.0000000000000455}, author = {Slowikowski, Kamil and Wei, Kevin and Brenner, Michael B and Raychaudhuri, Soumya} } @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 {7, title = {The Effects of Migration and Assortative Mating on Admixture Linkage Disequilibrium.}, journal = {Genetics}, volume = {205}, year = {2017}, month = {2017 Jan}, pages = {375-383}, abstract = {

Statistical models in medical and population genetics typically assume that individuals assort randomly in a population. While this simplifies model complexity, it contradicts an increasing body of evidence of nonrandom mating in human populations. Specifically, it has been shown that assortative mating is significantly affected by genomic ancestry. In this work, we examine the effects of ancestry-assortative mating on the linkage disequilibrium between local ancestry tracks of individuals in an admixed population. To accomplish this, we develop an extension to the Wright-Fisher model that allows for ancestry-based assortative mating. We show that ancestry-assortment perturbs the distribution of local ancestry linkage disequilibrium (LAD) and the variance of ancestry in a population as a function of the number of generations since admixture. This assortment effect can induce errors in demographic inference of admixed populations when methods assume random mating. We derive closed form formulae for LAD under an assortative-mating model with and without migration. We observe that LAD depends on the correlation of global ancestry of couples in each generation, the migration rate of each of the ancestral populations, the initial proportions of ancestral populations, and the number of generations since admixture. We also present the first direct evidence of ancestry-assortment in African Americans and examine LAD in simulated and real admixed population data of African Americans. We find that demographic inference under the assumption of random mating significantly underestimates the number of generations since admixture, and that accounting for assortative mating using the patterns of LAD results in estimates that more closely agrees with the historical narrative.

}, keywords = {Alleles, Black or African American, Datasets as Topic, Gene Frequency, Genetics, Population, Genomics, Human Migration, Humans, Linkage Disequilibrium, Models, Genetic, Models, Statistical, Polymorphism, Single Nucleotide}, issn = {1943-2631}, doi = {10.1534/genetics.116.192138}, author = {Zaitlen, Noah and Huntsman, Scott and Hu, Donglei and Spear, Melissa and Eng, Celeste and Oh, Sam S and White, Marquitta J and Mak, Angel and Davis, Adam and Meade, Kelly and Brigino-Buenaventura, Emerita and LeNoir, Michael A and Bibbins-Domingo, Kirsten and Burchard, Esteban G and Halperin, Eran} } @article {30, title = {The impact of rare variation on gene expression across tissues.}, journal = {Nature}, volume = {550}, year = {2017}, month = {2017 Oct 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 {26, title = {Improved methods for multi-trait fine mapping of pleiotropic risk loci.}, journal = {Bioinformatics}, volume = {33}, year = {2017}, month = {2017 Jan 15}, pages = {248-255}, abstract = {

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

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

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

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

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

}, keywords = {Animals, Disease, Gene-Environment Interaction, Genome-Wide Association Study, Genomics, Humans, Models, Animal, Sequence Analysis, RNA}, issn = {1476-6256}, doi = {10.1093/aje/kwx229}, author = {Ritchie, Marylyn D and Davis, Joe R and Aschard, Hugues and Battle, Alexis and Conti, David and Du, Mengmeng and Eskin, Eleazar and Fallin, M Daniele and Hsu, Li and Kraft, Peter and Moore, Jason H and Pierce, Brandon L and Bien, Stephanie A and Thomas, Duncan C and Wei, Peng and Montgomery, Stephen B} } @article {25, title = {VEXOR: an integrative environment for prioritization of functional variants in fine-mapping analysis.}, journal = {Bioinformatics}, volume = {33}, year = {2017}, month = {2017 May 01}, pages = {1389-1391}, abstract = {

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

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

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

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

SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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