Transcriptome-Wide Association Studies (TWAS): Methodologies, Applications, and Challenges.

Transcriptome-wide association study (TWAS) methodologies aim to identify genetic effects on phenotypes through the mediation of gene transcription. In TWAS, in silico models of gene expression are trained as functions of genetic variants and then applied to genome-wide association study (GWAS) data. This post-GWAS analysis identifies gene-trait associations with high interpretability, enabling follow-up functional genomics studies and the development of genetics-anchored resources.

Two polymorphic gene loci associated with treprostinil dose in pulmonary arterial hypertension.

Objective: Prostacyclin infusion for pulmonary arterial hypertension (PAH) is an effective therapy with varied dosing requirements and clinical response. The major aim of this study was to determine new biologically-based predictors of prostacyclin treatment response heterogeneity.

Methods: Ninety-eight patients with hemodynamically defined PAH at two academic medical centers volunteered for registry studies.

GRIK5 Genetically Regulated Expression Associated with Eye and Vascular Phenomes: Discovery through Iteration among Biobanks, Electronic Health Records, and Zebrafish.

Although the use of model systems for studying the mechanism of mutations that have a large effect is common, we highlight here the ways that zebrafish-model-system studies of a gene, GRIK5, that contributes to the polygenic liability to develop eye diseases have helped to illuminate a mechanism that implicates vascular biology in eye disease. A gene-expression prediction derived from a reference transcriptome panel applied to BioVU, a large electronic health record (EHR)-linked biobank at Vanderbilt University Medical Center, implicated reduced GRIK5 expression in diverse eye diseases. We tested the function of GRIK5 by depletion of its ortholog in zebrafish, and we observed reduced blood vessel numbers and integrity in the eye and increased vascular permeability.

Using an atlas of gene regulation across 44 human tissues to inform complex disease- and trait-associated variation.

We apply integrative approaches to expression quantitative loci (eQTLs) from 44 tissues from the Genotype-Tissue Expression project and genome-wide association study data. About 60% of known trait-associated loci are in linkage disequilibrium with a cis-eQTL, over half of which were not found in previous large-scale whole blood studies. Applying polygenic analyses to metabolic, cardiovascular, anthropometric, autoimmune, and neurodegenerative traits, we find that eQTLs are significantly enriched for trait associations in relevant pathogenic tissues and explain a substantial proportion of the heritability (40-80%).

Evidence of selection on splicing-associated loci in human populations and relevance to disease loci mapping.

We performed a whole-genome scan of genetic variants in splicing regulatory elements (SREs) and evaluated the extent to which natural selection has shaped extant patterns of variation in SREs. We investigated the degree of differentiation of single nucleotide polymorphisms (SNPs) in SREs among human populations and applied long-range haplotype- and multilocus allelic differentiation-based methods to detect selection signatures. We describe an approach, sampling a large number of loci across the genome from functional classes and using the consensus from multiple tests, for identifying candidates for selection signals.

Transcriptomic variation of pharmacogenes in multiple human tissues and lymphoblastoid cell lines.

Variation in the expression level and activity of genes involved in drug disposition and action ('pharmacogenes') can affect drug response and toxicity, especially when in tissues of pharmacological importance. Previous studies have relied primarily on microarrays to understand gene expression differences, or have focused on a single tissue or small number of samples. The goal of this study was to use RNA-sequencing (RNA-seq) to determine the expression levels and alternative splicing of 389 Pharmacogenomics Research Network pharmacogenes across four tissues (liver, kidney, heart and adipose) and lymphoblastoid cell lines, which are used widely in pharmacogenomics studies.

SCAN database: facilitating integrative analyses of cytosine modification and expression QTL.

Functional annotation of genetic variants including single nucleotide polymorphisms (SNPs) and copy number variations (CNV) promises to greatly improve our understanding of human complex traits. Previous transcriptomic studies involving individuals from different global populations have investigated the genetic architecture of gene expression variation by mapping expression quantitative trait loci (eQTL). Functional interpretation of genome-wide association studies (GWAS) has identified enrichment of eQTL in top signals from GWAS of human complex traits.

Cross-disorder genome-wide analyses suggest a complex genetic relationship between Tourette's syndrome and OCD.

Objective: Obsessive-compulsive disorder (OCD) and Tourette's syndrome are highly heritable neurodevelopmental disorders that are thought to share genetic risk factors. However, the identification of definitive susceptibility genes for these etiologically complex disorders remains elusive. The authors report a combined genome-wide association study (GWAS) of Tourette's syndrome and OCD.

Establishment of CYP2D6 reference samples by multiple validated genotyping platforms.

Cytochrome P450 2D6 (cytochrome P450, family 2, subfamily D, polypeptide 6 (CYP2D6)), a highly polymorphic drug-metabolizing enzyme, is involved in the metabolism of one-quarter of the most commonly prescribed medications. Here we have applied multiple genotyping methods and Sanger sequencing to assign precise and reproducible CYP2D6 genotypes, including copy numbers, for 48 HapMap samples. Furthermore, by analyzing a set of 50 human liver microsomes using endoxifen formation from N-desmethyl-tamoxifen as the phenotype of interest, we observed a significant positive correlation between CYP2D6 genotype-assigned activity score and endoxifen formation rate (rs = 0.

Integrating cell-based and clinical genome-wide studies to identify genetic variants contributing to treatment failure in neuroblastoma patients.

High-risk neuroblastoma is an aggressive malignancy, with high rates of treatment failure. We evaluated genetic variants associated with in vitro sensitivity to two derivatives of cyclophosphamide for association with clinical response in a separate replication cohort of neuroblastoma patients (n = 2,709). To determine sensitivity, lymphoblastoid cell lines (LCLs) were exposed to increasing concentrations of 4-hydroperoxycyclophosphamide (4HC; n = 422) and phosphoramide mustard (PM; n = 428).

Partitioning the heritability of Tourette syndrome and obsessive compulsive disorder reveals differences in genetic architecture.

The direct estimation of heritability from genome-wide common variant data as implemented in the program Genome-wide Complex Trait Analysis (GCTA) has provided a means to quantify heritability attributable to all interrogated variants. We have quantified the variance in liability to disease explained by all SNPs for two phenotypically-related neurobehavioral disorders, obsessive-compulsive disorder (OCD) and Tourette Syndrome (TS), using GCTA. Our analysis yielded a heritability point estimate of 0.

Genetic variants associated with warfarin dose in African-American individuals: a genome-wide association study.

Background: VKORC1 and CYP2C9 are important contributors to warfarin dose variability, but explain less variability for individuals of African descent than for those of European or Asian descent. We aimed to identify additional variants contributing to warfarin dose requirements in African Americans.

Methods: We did a genome-wide association study of discovery and replication cohorts.

Trans-population analysis of genetic mechanisms of ethnic disparities in neuroblastoma survival.

Background: Black patients with neuroblastoma have a higher prevalence of high-risk disease and worse outcome than white patients. We sought to investigate the relationship between genetic variation and the disparities in survival observed in neuroblastoma.

Methods: The analytic cohort was composed of 2709 patients.

Genome-wide association study of Tourette's syndrome.

Tourette's syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada.

Genome-wide association study of obsessive-compulsive disorder.

Obsessive-compulsive disorder (OCD) is a common, debilitating neuropsychiatric illness with complex genetic etiology. The International OCD Foundation Genetics Collaborative (IOCDF-GC) is a multi-national collaboration established to discover the genetic variation predisposing to OCD. A set of individuals affected with DSM-IV OCD, a subset of their parents, and unselected controls, were genotyped with several different Illumina SNP microarrays.

Genome-wide association and meta-analysis in populations from Starr County, Texas, and Mexico City identify type 2 diabetes susceptibility loci and enrichment for expression quantitative trait loci in top signals.

Aims/hypothesis: We conducted genome-wide association studies (GWASs) and expression quantitative trait loci (eQTL) analyses to identify and characterise risk loci for type 2 diabetes in Mexican-Americans from Starr County, TX, USA.

Method: Using 1.8 million directly interrogated and imputed genotypes in 837 unrelated type 2 diabetes cases and 436 normoglycaemic controls, we conducted Armitage trend tests.

Spoiling the whole bunch: quality control aimed at preserving the integrity of high-throughput genotyping.

False-positive or false-negative results attributable to undetected genotyping errors and confounding factors present a constant challenge for genome-wide association studies (GWAS) given the low signals associated with complex phenotypes and the noise associated with high-throughput genotyping. In the context of the genetics of kidneys in diabetes (GoKinD) study, we identify a source of error in genotype calling and demonstrate that a standard battery of quality-control (QC) measures is not sufficient to detect and/or correct it. We show that, if genotyping and calling are done by plate (batch), even a few DNA samples of marginally acceptable quality can profoundly alter the allele calls for other samples on the plate.

The neutrophil alloantigen HNA-3a (5b) is located on choline transporter-like protein 2 and appears to be encoded by an R>Q154 amino acid substitution.

The molecular basis of the HNA-3a/b (5b/a) leukocyte antigen system has not yet been defined despite evidence that HNA-3a-specific antibodies are particularly prone to cause severe, often fatal, transfusion-related lung injury. We used genome-wide single nucleotide polymorphism scanning and sequencing of DNA from persons of different HNA-3a/b phenotypes to identify a single single nucleotide polymorphism in exon 7 of the CLT2 gene (SLC44A2) that predicts an amino acid substitution in the first extracellular loop of choline transporter-like protein 2, a member of the choline transporter-like protein family of membrane glycoproteins, and correlates perfectly with HNA-3a/b phenotypes (R154 encodes HNA-3a; Q154 encodes HNA-3b). Mass spectrometric analysis of proteins immunoprecipitated from leukocytes by anti-HNA-3a provided direct evidence that anti-HNA-3a recognizes choline transporter-like protein 2.

SCAN: SNP and copy number annotation.

Motivation: Genome-wide association studies (GWAS) generate relationships between hundreds of thousands of single nucleotide polymorphisms (SNPs) and complex phenotypes. The contribution of the traditionally overlooked copy number variations (CNVs) to complex traits is also being actively studied. To facilitate the interpretation of the data and the designing of follow-up experimental validations, we have developed a database that enables the sensible prioritization of these variants by combining several approaches, involving not only publicly available physical and functional annotations but also multilocus linkage disequilibrium (LD) annotations as well as annotations of expression quantitative trait loci (eQTLs).

Identification of type 2 diabetes genes in Mexican Americans through genome-wide association studies.

Objective: The objective of this study was to identify DNA polymorphisms associated with type 2 diabetes in a Mexican-American population.

Research Design And Methods: We genotyped 116,204 single nucleotide polymorphisms (SNPs) in 281 Mexican Americans with type 2 diabetes and 280 random Mexican Americans from Starr County, Texas, using the Affymetrix GeneChip Human Mapping 100K set. Allelic association exact tests were calculated.