Survey of the Heritability and Sparse Architecture of Gene Expression Traits across Human Tissues.

Understanding the genetic architecture of gene expression traits is key to elucidating the underlying mechanisms of complex traits. Here, for the first time, we perform a systematic survey of the heritability and the distribution of effect sizes across all representative tissues in the human body. We find that local h2 can be relatively well characterized with 59% of expressed genes showing significant h2 (FDR < 0.1) in the DGN whole blood cohort. However, current sample sizes (n ≤ 922) do not allow us to compute distal h2. Bayesian Sparse Linear Mixed Model (BSLMM) analysis provides strong evidence that the genetic contribution to local expression traits is dominated by a handful of genetic variants rather than by the collective contribution of a large number of variants each of modest size. In other words, the local architecture of gene expression traits is sparse rather than polygenic across all 40 tissues (from DGN and GTEx) examined. This result is confirmed by the sparsity of optimal performing gene expression predictors via elastic net modeling. To further explore the tissue context specificity, we decompose the expression traits into cross-tissue and tissue-specific components using a novel Orthogonal Tissue Decomposition (OTD) approach. Through a series of simulations we show that the cross-tissue and tissue-specific components are identifiable via OTD. Heritability and sparsity estimates of these derived expression phenotypes show similar characteristics to the original traits. Consistent properties relative to prior GTEx multi-tissue analysis results suggest that these traits reflect the expected biology. Finally, we apply this knowledge to develop prediction models of gene expression traits for all tissues. The prediction models, heritability, and prediction performance R2 for original and decomposed expression phenotypes are made publicly available (https://github.com/hakyimlab/PrediXcan).

Poly-omic prediction of complex traits: OmicKriging.

High-confidence prediction of complex traits such as disease risk or drug response is an ultimate goal of personalized medicine. Although genome-wide association studies have discovered thousands of well-replicated polymorphisms associated with a broad spectrum of complex traits, the combined predictive power of these associations for any given trait is generally too low to be of clinical relevance. We propose a novel systems approach to complex trait prediction, which leverages and integrates similarity in genetic, transcriptomic, or other omics-level data. We translate the omic similarity into phenotypic similarity using a method called Kriging, commonly used in geostatistics and machine learning. Our method called OmicKriging emphasizes the use of a wide variety of systems-level data, such as those increasingly made available by comprehensive surveys of the genome, transcriptome, and epigenome, for complex trait prediction. Furthermore, our OmicKriging framework allows easy integration of prior information on the function of subsets of omics-level data from heterogeneous sources without the sometimes heavy computational burden of Bayesian approaches. Using seven disease datasets from the Wellcome Trust Case Control Consortium (WTCCC), we show that OmicKriging allows simple integration of sparse and highly polygenic components yielding comparable performance at a fraction of the computing time of a recently published Bayesian sparse linear mixed model method. Using a cellular growth phenotype, we show that integrating mRNA and microRNA expression data substantially increases performance over either dataset alone. Using clinical statin response, we show improved prediction over existing methods. We provide an R package to implement OmicKriging (http://www.scandb.org/newinterface/tools/OmicKriging.html).

Association of the GGCX (CAA)16/17 repeat polymorphism with higher warfarin dose requirements in African Americans.

OBJECTIVE: Little is known about genetic contributors to higher than usual warfarin dose requirements, particularly for African Americans. This study tested the hypothesis that the γ-glutamyl carboxylase (GGCX) genotype contributes to warfarin dose requirements greater than 7.5 mg/day in an African American population. METHODS: A total of 338 African Americans on a stable dose of warfarin were enrolled. The GGCX rs10654848 (CAA)n, rs12714145 (G>A), and rs699664 (p.R325Q); VKORC1 c.-1639G>A and rs61162043; and CYP2C9*2, *3, *5, *8, *11, and rs7089580 genotypes were tested for their association with dose requirements greater than 7.5 mg/day alone and in the context of other variables known to influence dose variability. RESULTS: The GGCX rs10654848 (CAA)16 or 17 repeat occurred at a frequency of 2.6% in African Americans and was overrepresented among patients requiring greater than 7.5 mg/day versus those who required lower doses (12 vs. 3%, P=0.003; odds ratio 4.0, 95% confidence interval, 1.5-10.5). The GGCX rs10654848 genotype remained associated with high dose requirements on regression analysis including age, body size, and VKORC1 genotype. On linear regression, the GGCX rs10654848 genotype explained 2% of the overall variability in warfarin dose in African Americans. An examination of the GGCX rs10654848 genotype in warfarin-treated Caucasians revealed a (CAA)16 repeat frequency of only 0.27% (P=0.008 compared with African Americans). CONCLUSION: These data support the GGCX rs10654848 genotype as a predictor of higher than usual warfarin doses in African Americans, who have a 10-fold higher frequency of the (CAA)16/17 repeat compared with Caucasians.

The composition of the microbiota modulates allograft rejection.

Transplantation is the only cure for end-stage organ failure, but without immunosuppression, T cells rapidly reject allografts. While genetic disparities between donor and recipient are major determinants of the kinetics of transplant rejection, little is known about the contribution of environmental factors. Because colonized organs have worse transplant outcome than sterile organs, we tested the influence of host and donor microbiota on skin transplant rejection. Compared with untreated conventional mice, pretreatment of donors and recipients with broad-spectrum antibiotics (Abx) or use of germ-free (GF) donors and recipients resulted in prolonged survival of minor antigen-mismatched skin grafts. Increased graft survival correlated with reduced type I IFN signaling in antigen-presenting cells (APCs) and decreased priming of alloreactive T cells. Colonization of GF mice with fecal material from untreated conventional mice, but not from Abx-pretreated mice, enhanced the ability of APCs to prime alloreactive T cells and accelerated graft rejection, suggesting that alloimmunity is modulated by the composition of microbiota rather than the quantity of bacteria. Abx pretreatment of conventional mice also delayed rejection of major antigen-mismatched skin and MHC class II-mismatched cardiac allografts. This study demonstrates that Abx pretreatment prolongs graft survival, suggesting that targeting microbial constituents is a potential therapeutic strategy for enhancing graft acceptance.

Novel single nucleotide polymorphism in CYP2C9 is associated with changes in warfarin clearance and CYP2C9 expression levels in African Americans.

Warfarin is a widely used anticoagulant whose active S-enantiomer is primarily metabolized by the CYP2C9 enzyme. The CYP2C9*2 and CYP2C9*3 alleles are associated with lower warfarin dose requirement and decreased enzyme activity. In contrast, we previously identified a novel single-nucleotide polymorphism (SNP) (rs7089580A > T) in CYP2C9 that is associated with higher warfarin dose requirement in African Americans (AAs). In this study, we examine the effect of rs7089580 on warfarin pharmacokinetics and CYP2C9 expression in 63 AA patients and 32 AA liver tissues, respectively. We found oral clearance of S-warfarin to be higher among carriers of the minor rs7089580 allele (T) compared with wild-type homozygotes (3.73 ± 1.46 vs 2.95 ± 1.39 mL/min; P = 0.04). CYP2C9 messenger RNA expression in liver tissue was also higher among A/T and T/T genotypes compared with A/A (P < 0.02). Our findings indicate that rs7089580 is associated with higher S-warfarin clearance and CYP2C9 expression and may help explain the higher dose requirement of warfarin in AAs. Furthermore, rs7089580 is in complete linkage disequilibrium with the promoter SNP rs12251841 in AAs, which may provide a biologically plausible explanation for the observed effect on CYP2C9 expression levels. Given the many clinically relevant substrates of CYP2C9, identifying polymorphisms that affect expression levels and metabolism across ethnicities is essential for individualization of doses with a narrow therapeutic index.

Commensal Bifidobacterium promotes antitumor immunity and facilitates anti-PD-L1 efficacy.

T cell infiltration of solid tumors is associated with favorable patient outcomes, yet the mechanisms underlying variable immune responses between individuals are not well understood. One possible modulator could be the intestinal microbiota. We compared melanoma growth in mice harboring distinct commensal microbiota and observed differences in spontaneous antitumor immunity, which were eliminated upon cohousing or after fecal transfer. Sequencing of the 16S ribosomal RNA identified Bifidobacterium as associated with the antitumor effects. Oral administration of Bifidobacterium alone improved tumor control to the same degree as programmed cell death protein 1 ligand 1 (PD-L1)-specific antibody therapy (checkpoint blockade), and combination treatment nearly abolished tumor outgrowth. Augmented dendritic cell function leading to enhanced CD8(+) T cell priming and accumulation in the tumor microenvironment mediated the effect. Our data suggest that manipulating the microbiota may modulate cancer immunotherapy.

A gene-based association method for mapping traits using reference transcriptome data.

Genome-wide association studies (GWAS) have identified thousands of variants robustly associated with complex traits. However, the biological mechanisms underlying these associations are, in general, not well understood. We propose a gene-based association method called PrediXcan that directly tests the molecular mechanisms through which genetic variation affects phenotype. The approach estimates the component of gene expression determined by an individual's genetic profile and correlates 'imputed' gene expression with the phenotype under investigation to identify genes involved in the etiology of the phenotype. Genetically regulated gene expression is estimated using whole-genome tissue-dependent prediction models trained with reference transcriptome data sets. PrediXcan enjoys the benefits of gene-based approaches such as reduced multiple-testing burden and a principled approach to the design of follow-up experiments. Our results demonstrate that PrediXcan can detect known and new genes associated with disease traits and provide insights into the mechanism of these associations.