eQTL Mapping Research Articles

Expression quantitative trait loci (eQTLs) in microRNA genes are enriched for schizophrenia and bipolar disorder association signals.

Schizophrenia (SZ) and bipolar disorder (BD) have substantial negative impact on the quality of human life. Both, microRNA (miRNA) expression profiling in SZ and BD postmortem brains [and genome-wide association studies (GWAS)] have implicated miRNAs in disease etiology. Here, we aim to determine whether significant GWAS signals observed in the Psychiatric Genetic Consortium (PGC) are enriched for miRNAs. A two-stage approach was used to determine whether association signals from PGC affect miRNAs: (i) statistical assessment of enrichment using a Simes test and sum of squares test (SST) and (ii) biological evidence that quantitative trait loci (eQTL) mapping to known miRNA genes affect their expression in an independent sample of 78 postmortem brains from the Stanley Medical Research Institute. A total of 2567 independent single nucleotide polymorphisms (SNPs) (R2 > 0.8) were mapped locally, within 1 Mb, to all known miRNAs (miRBase v. 21). We show robust enrichment for SZ- and BD-related SNPs with miRNAs using Simes (SZ: p ≤ 0.0023, BD: p ≤ 0.038), which remained significant after adjusting for background inflation in SZ (empirical p = 0.018) and approached significance in BD (empirical p = 0.07). At a false discovery rate of 10%, we identified a total of 32 eQTLs to influence miRNA expression; 11 of these overlapped with BD. Our approach of integrating PGC findings with eQTL results can be used to generate specific hypotheses regarding the role of miRNAs in SZ and BD.

Calling genotypes from public RNA-sequencing data enables identification of genetic variants that affect gene-expression levels.

BackgroundRNA-sequencing (RNA-seq) is a powerful technique for the identification of genetic variants that affect gene-expression levels, either through expression quantitative trait locus (eQTL) mapping or through allele-specific expression (ASE) analysis. Given increasing numbers of RNA-seq samples in the public domain, we here studied to what extent eQTLs and ASE effects can be identified when using public RNA-seq data while deriving the genotypes from the RNA-sequencing reads themselves.MethodsWe downloaded the raw reads for all available human RNA-seq datasets. Using these reads we performed gene expression quantification. All samples were jointly normalized and subjected to a strict quality control. We also derived genotypes using the RNA-seq reads and used imputation to infer non-coding variants. This allowed us to perform eQTL mapping and ASE analyses jointly on all samples that passed quality control. Our results were validated using samples for which DNA-seq genotypes were available.Results4,978 public human RNA-seq runs, representing many different tissues and cell-types, passed quality control. Even though these data originated from many different laboratories, samples reflecting the same cell type clustered together, suggesting that technical biases due to different sequencing protocols are limited. In a joint analysis on the 1,262 samples with high quality genotypes, we identified cis-eQTLs effects for 8,034 unique genes (at a false discovery rate ≤0.05). eQTL mapping on individual tissues revealed that a limited number of samples already suffice to identify tissue-specific eQTLs for known disease-associated genetic variants. Additionally, we observed strong ASE effects for 34 rare pathogenic variants, corroborating previously observed effects on the corresponding protein levels.ConclusionsBy deriving and imputing genotypes from RNA-seq data, it is possible to identify both eQTLs and ASE effects. Given the exponential growth of the number of publicly available RNA-seq samples, we expect this approach will become especially relevant for studying the effects of tissue-specific and rare pathogenic genetic variants to aid clinical interpretation of exome and genome sequencing.Electronic supplementary materialThe online version of this article (doi:10.1186/s13073-015-0152-4) contains supplementary material, which is available to authorized users.

Enrichment of inflammatory bowel disease and colorectal cancer risk variants in colon expression quantitative trait loci.

BackgroundGenome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) associated with diseases of the colon including inflammatory bowel diseases (IBD) and colorectal cancer (CRC). However, the functional role of many of these SNPs is largely unknown and tissue-specific resources are lacking. Expression quantitative trait loci (eQTL) mapping identifies target genes of disease-associated SNPs. This study provides a comprehensive eQTL map of distal colonic samples obtained from 40 healthy African Americans and demonstrates their relevance for GWAS of colonic diseases.Results8.4 million imputed SNPs were tested for their associations with 16,252 expression probes representing 12,363 unique genes. 1,941 significant cis-eQTL, corresponding to 122 independent signals, were identified at a false discovery rate (FDR) of 0.01. Overall, among colon cis-eQTL, there was significant enrichment for GWAS variants for IBD (Crohn’s disease [CD] and ulcerative colitis [UC]) and CRC as well as type 2 diabetes and body mass index. ERAP2, ADCY3, INPP5E, UBA7, SFMBT1, NXPE1 and REXO2 were identified as target genes for IBD-associated variants. The CRC-associated eQTL rs3802842 was associated with the expression of C11orf93 (COLCA2). Enrichment of colon eQTL near transcription start sites and for active histone marks was demonstrated, and eQTL with high population differentiation were identified.ConclusionsThrough the comprehensive study of eQTL in the human colon, this study identified novel target genes for IBD- and CRC-associated genetic variants. Moreover, bioinformatic characterization of colon eQTL provides a tissue-specific tool to improve understanding of biological differences in diseases between different ethnic groups.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1292-z) contains supplementary material, which is available to authorized users.

The genetic architecture of gene expression levels in wild baboons.

Primate evolution has been argued to result, in part, from changes in how genes are regulated. However, we still know little about gene regulation in natural primate populations. We conducted an RNA sequencing (RNA-seq)-based study of baboons from an intensively studied wild population. We performed complementary expression quantitative trait locus (eQTL) mapping and allele-specific expression analyses, discovering substantial evidence for, and surprising power to detect, genetic effects on gene expression levels in the baboons. eQTL were most likely to be identified for lineage-specific, rapidly evolving genes; interestingly, genes with eQTL significantly overlapped between baboons and a comparable human eQTL data set. Our results suggest that genes vary in their tolerance of genetic perturbation, and that this property may be conserved across species. Further, they establish the feasibility of eQTL mapping using RNA-seq data alone, and represent an important step towards understanding the genetic architecture of gene expression in primates.

Fast and robust group-wise eQTL mapping using sparse graphical models.

BackgroundGenome-wide expression quantitative trait loci (eQTL) studies have emerged as a powerful tool to understand the genetic basis of gene expression and complex traits. The traditional eQTL methods focus on testing the associations between individual single-nucleotide polymorphisms (SNPs) and gene expression traits. A major drawback of this approach is that it cannot model the joint effect of a set of SNPs on a set of genes, which may correspond to hidden biological pathways.ResultsWe introduce a new approach to identify novel group-wise associations between sets of SNPs and sets of genes. Such associations are captured by hidden variables connecting SNPs and genes. Our model is a linear-Gaussian model and uses two types of hidden variables. One captures the set associations between SNPs and genes, and the other captures confounders. We develop an efficient optimization procedure which makes this approach suitable for large scale studies. Extensive experimental evaluations on both simulated and real datasets demonstrate that the proposed methods can effectively capture both individual and group-wise signals that cannot be identified by the state-of-the-art eQTL mapping methods.ConclusionsConsidering group-wise associations significantly improves the accuracy of eQTL mapping, and the successful multi-layer regression model opens a new approach to understand how multiple SNPs interact with each other to jointly affect the expression level of a group of genes.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-014-0421-z) contains supplementary material, which is available to authorized users.

A systematic heritability analysis of the human whole blood transcriptome.

Genome-wide expression quantitative trait locus (eQTL) mapping may reveal common genetic variants regulating gene expression. In addition to mapping eQTLs, we systematically evaluated the heritability of the whole blood transcriptome in 5,626 participants from the Framingham Heart Study. Of all gene expression measurements, about 40 % exhibit evidence of being heritable [hgeneExp(2) > 0, (p < 0.05)], the average heritability was estimated to be 0.13, and 10 % display hgeneExp(2) > 0.2. To identify the role of eQTLs in promoting phenotype differences and disease susceptibility, we investigated the proportion of cis/trans eQTLs in different heritability categories and discovered that genes with higher heritability are more likely to have cis eQTLs that explain large proportions of variance in the expression of the corresponding genes. Single cis eQTLs explain 0.33-0.53 of variance in transcripts on average, whereas single trans eQTLs only explain 0.02-0.07. The top cis eQTLs tend to explain more variance in the corresponding gene when its hgeneExp(2) is greater. Taking body mass index (BMI) as a case study, we cross-linked cis/trans eQTLs with both GWAS SNPs and differentially expressed genes for BMI. We discovered that BMI GWAS SNPs in 16p11.2 (e.g., rs7359397) are associated with several BMI differentially expressed genes in a cis manner (e.g. SULT1A1, SPNS1, and TUFM). These BMI signature genes explain a much larger proportion of variance in BMI than do the GWAS SNPs. Our results shed light on the impact of eQTLs on the heritability of the human whole blood transcriptome and its relations to phenotype differences.

Proper Use of Allele-Specific Expression Improves Statistical Power for cis -eQTL Mapping with RNA-Seq Data

Proper Use of Allele-Specific Expression Improves Statistical Power for cis -eQTL Mapping with RNA-Seq Data

Quality Control and Analysis of NGS RNA Sequencing Data.

Transcriptome sequencing, where RNA is isolated, converted to library of cDNA fragments, and sequenced using next-generation sequencing technology, has become the method of choice for the genome-wide characterization of mRNA levels. It offers a more accurate quantification of transcript levels than array-based methods, but also has the added benefit of allowing the discovery of novel gene/transcripts, alternative splice junctions, and novel RNAs. In addition, RNA sequencing may be used to investigate differential gene expression, allelic imbalance, eQTL mapping, RNA editing, RNA-protein interactions, and alternative splicing. A number of statistical methods and tools are available for differential expression analysis using RNA sequencing data and these are continually being developed and improved to handle more complex experimental designs. This chapter describes an example workflow for the quality control and analysis of raw RNA sequencing reads for the purposes of differential gene expression analysis, followed by pathway/enrichment analysis of significantly different genes. The methods and tools described are just one example of how this analysis can be conducted, but they can be applied to most standard RNA sequencing studies of differential gene expression. The methods covered are based on Illumina HiSeq single-end 50 bp reads. However, all programs used are capable of working with paired-end data, subsequent to minor adaptations.

Fine-mapping the 2q37 and 17q11.2-q22 loci for novel genes and sequence variants associated with a genetic predisposition to prostate cancer.

The 2q37 and 17q12-q22 loci are linked to an increased prostate cancer (PrCa) risk. No candidate gene has been localized at 2q37 and the HOXB13 variant G84E only partially explains the linkage to 17q21-q22 observed in Finland. We screened these regions by targeted DNA sequencing to search for cancer-associated variants. Altogether, four novel susceptibility alleles were identified. Two ZNF652 (17q21.3) variants, rs116890317 and rs79670217, increased the risk of both sporadic and hereditary PrCa (rs116890317: OR = 3.3-7.8, p = 0.003-3.3 × 10(-5) ; rs79670217: OR = 1.6-1.9, p = 0.002-0.009). The HDAC4 (2q37.2) variant rs73000144 (OR = 14.6, p = 0.018) and the EFCAB13 (17q21.3) variant rs118004742 (OR = 1.8, p = 0.048) were overrepresented in patients with familial PrCa. To map the variants within 2q37 and 17q11.2-q22 that may regulate PrCa-associated genes, we combined DNA sequencing results with transcriptome data obtained by RNA sequencing. This expression quantitative trait locus (eQTL) analysis identified 272 single-nucleotide polymorphisms (SNPs) possibly regulating six genes that were differentially expressed between cases and controls. In a modified approach, prefiltered PrCa-associated SNPs were exploited and interestingly, a novel eQTL targeting ZNF652 was identified. The novel variants identified in this study could be utilized for PrCa risk assessment, and they further validate the suggested role of ZNF652 as a PrCa candidate gene. The regulatory regions discovered by eQTL mapping increase our understanding of the relationship between regulation of gene expression and susceptibility to PrCa and provide a valuable starting point for future functional research.

A splice mutation in the PHKG1 gene causes high glycogen content and low meat quality in pig skeletal muscle.

Glycolytic potential (GP) in skeletal muscle is economically important in the pig industry because of its effect on pork processing yield. We have previously mapped a major quantitative trait loci (QTL) for GP on chromosome 3 in a White Duroc × Erhualian F2 intercross. We herein performed a systems genetic analysis to identify the causal variant underlying the phenotype QTL (pQTL). We first conducted genome-wide association analyses in the F2 intercross and an F19 Sutai pig population. The QTL was then refined to an 180-kb interval based on the 2-LOD drop method. We then performed expression QTL (eQTL) mapping using muscle transcriptome data from 497 F2 animals. Within the QTL interval, only one gene (PHKG1) has a cis-eQTL that was colocolizated with pQTL peaked at the same SNP. The PHKG1 gene encodes a catalytic subunit of the phosphorylase kinase (PhK), which functions in the cascade activation of glycogen breakdown. Deep sequencing of PHKG1 revealed a point mutation (C>A) in a splice acceptor site of intron 9, resulting in a 32-bp deletion in the open reading frame and generating a premature stop codon. The aberrant transcript induces nonsense-mediated decay, leading to lower protein level and weaker enzymatic activity in affected animals. The mutation causes an increase of 43% in GP and a decrease of>20% in water-holding capacity of pork. These effects were consistent across the F2 and Sutai populations, as well as Duroc × (Landrace × Yorkshire) hybrid pigs. The unfavorable allele exists predominantly in Duroc-derived pigs. The findings provide new insights into understanding risk factors affecting glucose metabolism, and would greatly contribute to the genetic improvement of meat quality in Duroc related pigs.

Mapping eQTL networks with mixed graphical Markov models.

Expression quantitative trait loci (eQTL) mapping constitutes a challenging problem due to, among other reasons, the high-dimensional multivariate nature of gene-expression traits. Next to the expression heterogeneity produced by confounding factors and other sources of unwanted variation, indirect effects spread throughout genes as a result of genetic, molecular, and environmental perturbations. From a multivariate perspective one would like to adjust for the effect of all of these factors to end up with a network of direct associations connecting the path from genotype to phenotype. In this article we approach this challenge with mixed graphical Markov models, higher-order conditional independences, and q-order correlation graphs. These models show that additive genetic effects propagate through the network as function of gene-gene correlations. Our estimation of the eQTL network underlying a well-studied yeast data set leads to a sparse structure with more direct genetic and regulatory associations that enable a straightforward comparison of the genetic control of gene expression across chromosomes. Interestingly, it also reveals that eQTLs explain most of the expression variability of network hub genes.

Novel distal eQTL analysis demonstrates effect of population genetic architecture on detecting and interpreting associations.

Mapping expression quantitative trait loci (eQTL) has identified genetic variants associated with transcription rates and has provided insight into genotype–phenotype associations obtained from genome-wide association studies (GWAS). Traditional eQTL mapping methods present significant challenges for the multiple-testing burden, resulting in a limited ability to detect eQTL that reside distal to the affected gene. To overcome this, we developed a novel eQTL testing approach, “network-based, large-scale identification of distal eQTL” (NetLIFT), which performs eQTL testing based on the pairwise conditional dependencies between genes’ expression levels. When applied to existing data from yeast segregants, NetLIFT replicated most previously identified distal eQTL and identified 46% more genes with distal effects compared to local effects. In liver data from mouse lines derived through the Collaborative Cross project, NetLIFT detected 5744 genes with local eQTL while 3322 genes had distal eQTL. This analysis revealed founder-of-origin effects for a subset of local eQTL that may contribute to previously described phenotypic differences in metabolic traits. In human lymphoblastoid cell lines, NetLIFT was able to detect 1274 transcripts with distal eQTL that had not been reported in previous studies, while 2483 transcripts with local eQTL were identified. In all species, we found no enrichment for transcription factors facilitating eQTL associations; instead, we found that most trans-acting factors were annotated for metabolic function, suggesting that genetic variation may indirectly regulate multigene pathways by targeting key components of feedback processes within regulatory networks. Furthermore, the unique genetic history of each population appears to influence the detection of genes with local and distal eQTL.

Systematic annotation of celiac disease loci refines pathological pathways and suggests a genetic explanation for increased interferon-gamma levels.

Although genome-wide association studies and fine mapping have identified 39 non-HLA loci associated with celiac disease (CD), it is difficult to pinpoint the functional variants and susceptibility genes in these loci. We applied integrative approaches to annotate and prioritize functional single nucleotide polymorphisms (SNPs), genes and pathways affected in CD. CD-associated SNPs were intersected with regulatory elements categorized by the ENCODE project to prioritize functional variants, while results from cis-expression quantitative trait loci (eQTL) mapping in 1469 blood samples were combined with co-expression analyses to prioritize causative genes. To identify the key cell types involved in CD, we performed pathway analysis on RNA-sequencing data from different immune cell populations and on publicly available expression data on non-immune tissues. We discovered that CD SNPs are significantly enriched in B-cell-specific enhancer regions, suggesting that, besides T-cell processes, B-cell responses play a major role in CD. By combining eQTL and co-expression analyses, we prioritized 43 susceptibility genes in 36 loci. Pathway and tissue-specific expression analyses on these genes suggested enrichment of CD genes in the Th1, Th2 and Th17 pathways, but also predicted a role for four genes in the intestinal barrier function. We also discovered an intricate transcriptional connectivity between CD susceptibility genes and interferon-γ, a key effector in CD, despite the absence of CD-associated SNPs in the IFNG locus. Using systems biology, we prioritized the CD-associated functional SNPs and genes. By highlighting a role for B cells in CD, which classically has been described as a T-cell-driven disease, we offer new insights into the mechanisms and pathways underlying CD.

Genetical and comparative genomics of Brassica under altered Ca supply identifies Arabidopsis Ca-transporter orthologs.

Although Ca transport in plants is highly complex, the overexpression of vacuolar Ca(2+) transporters in crops is a promising new technology to improve dietary Ca supplies through biofortification. Here, we sought to identify novel targets for increasing plant Ca accumulation using genetical and comparative genomics. Expression quantitative trait locus (eQTL) mapping to 1895 cis- and 8015 trans-loci were identified in shoots of an inbred mapping population of Brassica rapa (IMB211 × R500); 23 cis- and 948 trans-eQTLs responded specifically to altered Ca supply. eQTLs were screened for functional significance using a large database of shoot Ca concentration phenotypes of Arabidopsis thaliana. From 31 Arabidopsis gene identifiers tagged to robust shoot Ca concentration phenotypes, 21 mapped to 27 B. rapa eQTLs, including orthologs of the Ca(2+) transporters At-CAX1 and At-ACA8. Two of three independent missense mutants of BraA.cax1a, isolated previously by targeting induced local lesions in genomes, have allele-specific shoot Ca concentration phenotypes compared with their segregating wild types. BraA.CAX1a is a promising target for altering the Ca composition of Brassica, consistent with prior knowledge from Arabidopsis. We conclude that multiple-environment eQTL analysis of complex crop genomes combined with comparative genomics is a powerful technique for novel gene identification/prioritization.

Identification and Validation of Genetic Variants that Influence Transcription Factor and Cell Signaling Protein Levels

Many genetic variants associated with human disease have been found to be associated with alterations in mRNA expression. Although it is commonly assumed that mRNA expression changes will lead to consequent changes in protein levels, methodological challenges have limited our ability to test the degree to which this assumption holds true. Here, we further developed the micro-western array approach and globally examined relationships between human genetic variation and cellular protein levels. We collected more than 250,000 protein level measurements comprising 441 transcription factor and signaling protein isoforms across 68 Yoruba (YRI) HapMap lymphoblastoid cell lines (LCLs) and identified 12 cis and 160 trans protein level QTLs (pQTLs) at a false discovery rate (FDR) of 20%. Whereas up to two thirds of cis mRNA expression QTLs (eQTLs) were also pQTLs, many pQTLs were not associated with mRNA expression. Notably, we replicated and functionally validated a trans pQTL relationship between the KARS lysyl-tRNA synthetase locus and levels of the DIDO1 protein. This study demonstrates proof of concept in applying an antibody-based microarray approach to iteratively measure the levels of human proteins and relate these levels to human genome variation and other genomic data sets. Our results suggest that protein-based mechanisms might functionally buffer genetic alterations that influence mRNA expression levels and that pQTLs might contribute phenotypic diversity to a human population independently of influences on mRNA expression.

Large-scale East-Asian eQTL mapping reveals novel candidate genes for LD mapping and the genomic landscape of transcriptional effects of sequence variants.

Profiles of sequence variants that influence gene transcription are very important for understanding mechanisms that affect phenotypic variation and disease susceptibility. Using genotypes at 1.4 million SNPs and a comprehensive transcriptional profile of 15,454 coding genes and 6,113 lincRNA genes obtained from peripheral blood cells of 298 Japanese individuals, we mapped expression quantitative trait loci (eQTLs). We identified 3,804 cis-eQTLs (within 500 kb from target genes) and 165 trans-eQTLs (>500 kb away or on different chromosomes). Cis-eQTLs were often located in transcribed or adjacent regions of genes; among these regions, 5′ untranslated regions and 5′ flanking regions had the largest effects. Epigenetic evidence for regulatory potential accumulated in public databases explained the magnitude of the effects of our eQTLs. Cis-eQTLs were often located near the respective target genes, if not within genes. Large effect sizes were observed with eQTLs near target genes, and effect sizes were obviously attenuated as the eQTL distance from the gene increased. Using a very stringent significance threshold, we identified 165 large-effect trans-eQTLs. We used our eQTL map to assess 8,069 disease-associated SNPs identified in 1,436 genome-wide association studies (GWAS). We identified genes that might be truly causative, but GWAS might have failed to identify for 148 out of the GWAS-identified SNPs; for example, TUFM (P = 3.3E-48) was identified for inflammatory bowel disease (early onset); ZFP90 (P = 4.4E-34) for ulcerative colitis; and IDUA (P = 2.2E-11) for Parkinson's disease. We identified four genes (P<2.0E-14) that might be related to three diseases and two hematological traits; each expression is regulated by trans-eQTLs on a different chromosome than the gene.

A genome-wide survey of CD4+ lymphocyte regulatory genetic variants identifies novel asthma genes

Genome-wide association studies have yet to identify the majority of genetic variants involved in asthma. We hypothesized that expression quantitative trait locus (eQTL) mapping can identify novel asthma genes by enabling prioritization of putative functional variants for association testing. We evaluated 6706 cis-acting expression-associated variants (eSNPs) identified through a genome-wide eQTL survey of CD4(+) lymphocytes for association with asthma. eSNPs were tested for association with asthma in 359 asthmatic patients and 846 control subjects from the Childhood Asthma Management Program, with verification by using family-based testing. Significant associations were tested for replication in 579 parent-child trios with asthma from Costa Rica. Further functional validation was performed by using formaldehyde-assisted isolation of regulatory elements (FAIRE) quantitative PCR and chromatin immunoprecipitation PCR in lung-derived epithelial cell lines (Beas-2B and A549) and Jurkat cells, a leukemia cell line derived from T lymphocytes. Cis-acting eSNPs demonstrated associations with asthma in both cohorts. We confirmed the previously reported association of ORMDL3/GSDMB variants with asthma (combined P = 2.9 × 10(-8)). Reproducible associations were also observed for eSNPs in 3 additional genes: fatty acid desaturase 2 (FADS2; P = .002), N-acetyl-α-D-galactosaminidase (NAGA; P = .0002), and Factor XIII, A1 (F13A1; P = .0001). Subsequently, we demonstrated that FADS2 mRNA is increased in CD4(+) lymphocytes in asthmatic patients and that the associated eSNPs reside within DNA segments with histone modifications that denote open chromatin status and confer enhancer activity. Our results demonstrate the utility of eQTL mapping in the identification of novel asthma genes and provide evidence for the importance of FADS2, NAGA, and F13A1 in the pathogenesis of asthma.

Graph-regularized dual Lasso for robust eQTL mapping.

Motivation: As a promising tool for dissecting the genetic basis of complex traits, expression quantitative trait loci (eQTL) mapping has attracted increasing research interest. An important issue in eQTL mapping is how to effectively integrate networks representing interactions among genetic markers and genes. Recently, several Lasso-based methods have been proposed to leverage such network information. Despite their success, existing methods have three common limitations: (i) a preprocessing step is usually needed to cluster the networks; (ii) the incompleteness of the networks and the noise in them are not considered; (iii) other available information, such as location of genetic markers and pathway information are not integrated.Results: To address the limitations of the existing methods, we propose Graph-regularized Dual Lasso (GDL), a robust approach for eQTL mapping. GDL integrates the correlation structures among genetic markers and traits simultaneously. It also takes into account the incompleteness of the networks and is robust to the noise. GDL utilizes graph-based regularizers to model the prior networks and does not require an explicit clustering step. Moreover, it enables further refinement of the partial and noisy networks. We further generalize GDL to incorporate the location of genetic makers and gene-pathway information. We perform extensive experimental evaluations using both simulated and real datasets. Experimental results demonstrate that the proposed methods can effectively integrate various available priori knowledge and significantly outperform the state-of-the-art eQTL mapping methods.Availability: Software for both C++ version and Matlab version is available at http://www.cs.unc.edu/∼weicheng/.Contact: weiwang@cs.ucla.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Genome-Wide Identification of Expression Quantitative Trait Loci (eQTLs) in Human Heart

In recent years genome-wide association studies (GWAS) have uncovered numerous chromosomal loci associated with various electrocardiographic traits and cardiac arrhythmia predisposition. A considerable fraction of these loci lie within inter-genic regions. The underlying trait-associated variants likely reside in regulatory regions and exert their effect by modulating gene expression. Hence, the key to unraveling the molecular mechanisms underlying these cardiac traits is to interrogate variants for association with differential transcript abundance by expression quantitative trait locus (eQTL) analysis. In this study we conducted an eQTL analysis of human heart. For a total of 129 left ventricular samples that were collected from non-diseased human donor hearts, genome-wide transcript abundance and genotyping was determined using microarrays. Each of the 18,402 transcripts and 897,683 SNP genotypes that remained after pre-processing and stringent quality control were tested for eQTL effects. We identified 771 eQTLs, regulating 429 unique transcripts. Overlaying these eQTLs with cardiac GWAS loci identified novel candidates for studies aimed at elucidating the functional and transcriptional impact of these loci. Thus, this work provides for the first time a comprehensive eQTL map of human heart: a powerful and unique resource that enables systems genetics approaches for the study of cardiac traits.

Genetic dissection of the Drosophila melanogaster female head transcriptome reveals widespread allelic heterogeneity.

Modern genetic mapping is plagued by the “missing heritability” problem, which refers to the discordance between the estimated heritabilities of quantitative traits and the variance accounted for by mapped causative variants. One major potential explanation for the missing heritability is allelic heterogeneity, in which there are multiple causative variants at each causative gene with only a fraction having been identified. The majority of genome-wide association studies (GWAS) implicitly assume that a single SNP can explain all the variance for a causative locus. However, if allelic heterogeneity is prevalent, a substantial amount of genetic variance will remain unexplained. In this paper, we take a haplotype-based mapping approach and quantify the number of alleles segregating at each locus using a large set of 7922 eQTL contributing to regulatory variation in the Drosophila melanogaster female head. Not only does this study provide a comprehensive eQTL map for a major community genetic resource, the Drosophila Synthetic Population Resource, but it also provides a direct test of the allelic heterogeneity hypothesis. We find that 95% of cis-eQTLs and 78% of trans-eQTLs are due to multiple alleles, demonstrating that allelic heterogeneity is widespread in Drosophila eQTL. Allelic heterogeneity likely contributes significantly to the missing heritability problem common in GWAS studies.