Revealing the architecture of gene regulation: the promise of eQTL studies

Expression quantitative trait loci (eQTL) mapping studies have become a widely used tool for identifying genetic variants that affect gene regulation. In these studies, expression levels are viewed as quantitative traits, and gene expression phenotypes are mapped to particular genomic loci by combining studies of variation in gene expression patterns with genome-wide genotyping. Results from recent eQTL mapping studies have revealed substantial heritable variation in gene expression within and between populations. In many cases, genetic factors that influence gene expression levels can be mapped to proximal (putatively cis) eQTLs and, less often, to distal (putatively trans) eQTLs. Beyond providing great insight into the biology of gene regulation, a combination of eQTL studies with results from traditional linkage or association studies of human disease may help predict a specific regulatory role for polymorphic sites previously associated with disease.     

THE DETECTION OF PLASTICITY GENES IN HETEROGENEOUS ENVIRONMENTS

The molecular genetic mechanisms for phenotypic plasticity across heterogeneous macro- and microenvironments were examined using the Populus genomic map constructed by DNA-based markers. Three hypotheses have been suggested to explain genetic variation in phenotypic response to varying environments (i.e., reaction norm): Lerner's homeostasis, allelic sensitivity, and gene regulation. The homeostasis hypothesis, which predicts that heterozygotes are less sensitive to the environment than homozygotes, was supported for phenotypic plasticity to unpredictable environments (microenvironmental plasticity) at the whole-genome level, but for phenotypic plasticity to predictable environments (macroenvironmental plasticity) the hypothesis was supported only at functioning quantitative trait loci (QTLs). For all growth traits studied, gene regulation was suggested to play a prevailing role in determining the norms of reaction to environments. Indirect evidence for gene regulation is that there tend to be more QTLs with larger effects on the phenotype in optimal growing conditions than suboptimal growing conditions because the expression of these QTLs identified is mediated by regulatory genes. Direct evidence for gene regulation is the identification of some loci that differ from QTLs for trait values within environments and exert an environmentally dependent control over structural gene expression. In this study, fewer environmentally sensitive QTLs were detected that display unparalleled allelic effects across environments. For stem height, there were more regulatory loci and more structural loci (whose expression is determined by gene regulation) affecting phenotypic plasticity than for basal area. It was found that microenvironmental plasticity was likely controlled by different genetic systems than those for macroenvironmental plasticity.     

Advances in quantitative trait analysis in yeast

Understanding the genetic mechanisms underlying complex traits is one of the next frontiers in biology. The budding yeast Saccharomyces cerevisiae has become an important model for elucidating the mechanisms that govern natural genetic and phenotypic variation. This success is partially due to its intrinsic biological features, such as the short sexual generation time, high meiotic recombination rate, and small genome size. Precise reverse genetics technologies allow the high throughput manipulation of genetic information with exquisite precision, offering the unique opportunity to experimentally measure the phenotypic effect of genetic variants. Population genomic and phenomic studies have revealed widespread variation between diverged populations, characteristic of man-made environments, as well as geographic clusters of wild strains along with naturally occurring recombinant strains (mosaics). Here, we review these recent studies and provide a perspective on how these previously unappreciated levels of variation can help to bridge our understanding of the genotype-phenotype gap, keeping budding yeast at the forefront of genetic studies. Not only are quantitative trait loci (QTL) being mapped with high resolution down to the nucleotide, for the first time QTLs of modest effect and complex interactions between these QTLs and between QTLs and the environment are being determined experimentally at unprecedented levels using next generation techniques of deep sequencing selected pools of individuals as well as multi-generational crosses.     

Rare and common regulatory variation in population-scale sequenced human genomes

Population-scale genome sequencing allows the characterization of functional effects of a broad spectrum of genetic variants underlying human phenotypic variation. Here, we investigate the influence of rare and common genetic variants on gene expression patterns, using variants identified from sequencing data from the 1000 genomes project in an African and European population sample and gene expression data from lymphoblastoid cell lines. We detect comparable numbers of expression quantitative trait loci (eQTLs) when compared to genotypes obtained from HapMap 3, but as many as 80% of the top expression quantitative trait variants (eQTVs) discovered from 1000 genomes data are novel. The properties of the newly discovered variants suggest that mapping common causal regulatory variants is challenging even with full resequencing data; however, we observe significant enrichment of regulatory effects in splice-site and nonsense variants. Using RNA sequencing data, we show that 46.2% of nonsynonymous variants are differentially expressed in at least one individual in our sample, creating widespread potential for interactions between functional protein-coding and regulatory variants. We also use allele-specific expression to identify putative rare causal regulatory variants. Furthermore, we demonstrate that outlier expression values can be due to rare variant effects, and we approximate the number of such effects harboured in an individual by effect size. Our results demonstrate that integration of genomic and RNA sequencing analyses allows for the joint assessment of genome sequence and genome function.     

Genome-wide Epigenetic Data Facilitate Understanding of Disease Susceptibility Association Studies

Complex traits such as susceptibility to diseases are determined in part by variants at multiple genetic loci. Genome-wide association studies can identify these loci, but most phenotype-associated variants lie distal to protein-coding regions and are likely involved in regulating gene expression. Understanding how these genetic variants affect complex traits depends on the ability to predict and test the function of the genomic elements harboring them. Community efforts such as the ENCODE Project provide a wealth of data about epigenetic features associated with gene regulation. These data enable the prediction of testable functions for many phenotype-associated variants.     

Evolution of pleiotropy: epistatic interaction pattern supports a mechanistic model underlying variation in genotype-phenotype map

The genotype-phenotype (GP) map consists of developmental and physiological mechanisms mapping genetic onto phenotypic variation. It determines the distribution of heritable phenotypic variance on which selection can act. Comparative studies of morphology as well as of gene regulatory networks show that the GP map itself evolves, yet little is known about the actual evolutionary mechanisms involved. The study of such mechanisms requires exploring the variation in GP maps at the population level, which presently is easier to quantify by statistical genetic methods rather than by regulatory network structures. We focus on the evolution of pleiotropy, a major structural aspect of the GP map. Pleiotropic genes affect multiple traits and underlie genetic covariance between traits, often causing evolutionary constraints. Previous quantitative genetic studies have demonstrated population-level variation in pleiotropy in the form of loci, at which genotypes differ in the genetic covariation between traits. This variation can potentially fuel evolution of the GP map under selection and/or drift. Here, we propose a developmental mechanism underlying population genetic variation in covariance and test its predictions. Specifically, the mechanism predicts that the loci identified as responsible for genetic variation in pleiotropy are involved in trait-specific epistatic interactions. We test this prediction for loci affecting allometric relationships between traits in an advanced intercross between inbred mouse strains. The results consistently support the prediction. We further find a high degree of sign epistasis in these interactions, which we interpret as an indication of adaptive gene complexes within the diverged parental lines.     

Determining causality and consequence of expression quantitative trait loci

Expression quantitative trait loci (eQTLs) are currently the most abundant and systematically-surveyed class of functional consequence for genetic variation. Recent genetic studies of gene expression have identified thousands of eQTLs in diverse tissue types for the majority of human genes. Application of this large eQTL catalog provides an important resource for understanding the molecular basis of common genetic diseases. However, only now has both the availability of individuals with full genomes and corresponding advances in functional genomics provided the opportunity to dissect eQTLs to identify causal regulatory variants. Resolving the properties of such causal regulatory variants is improving understanding of the molecular mechanisms that influence traits and guiding the development of new genome-scale approaches to variant interpretation. In this review, we provide an overview of current computational and experimental methods for identifying causal regulatory variants and predicting their phenotypic consequences.     

To clone or not to clone plant QTLs: present and future challenges

Recent technical advancements and refinement of analytical methods have enabled the loci (quantitative trait loci, QTLs) responsible for the genetic control of quantitative traits to be dissected molecularly. To date, most plant QTLs have been cloned using a positional cloning approach following identification in experimental crosses. In some cases, an association between sequence variation at a candidate gene and a phenotype has been established by analysing existing genetic accessions. These strategies can be refined using appropriate genetic materials and the latest developments in genomics platforms. We foresee that although QTL analysis and cloning addressing naturally occurring genetic variation should shed light on mechanisms of plant adaptation, a greater emphasis on approaches relying on mutagenesis and candidate gene validation is likely to accelerate the pace of discovering the genes underlying QTLs.     

Identifying the molecular basis of QTLs: eQTLs add a new dimension

Natural genetic variation within plant species is at the core of plant science ranging from agriculture to evolution. Whereas much progress has been made in mapping quantitative trait loci (QTLs) controlling this natural variation, the elucidation of the underlying molecular mechanisms has remained a bottleneck. Recent systems biology tools have significantly shortened the time required to proceed from a mapped locus to testing of candidate genes. These tools enable research on natural variation to move from simple reductionistic studies focused on individual genes to integrative studies connecting molecular variation at multiple loci with physiological consequences. This review focuses on recent examples that demonstrate how expression QTL data can be used for gene discovery and exploited to untangle complex regulatory networks.     

Aberrant Gene Expression in Humans

Gene expression as an intermediate molecular phenotype has been a focus of research interest. In particular, studies of expression quantitative trait loci (eQTL) have offered promise for understanding gene regulation through the discovery of genetic variants that explain variation in gene expression levels. Existing eQTL methods are designed for assessing the effects of common variants, but not rare variants. Here, we address the problem by establishing a novel analytical framework for evaluating the effects of rare or private variants on gene expression. Our method starts from the identification of outlier individuals that show markedly different gene expression from the majority of a population, and then reveals the contributions of private SNPs to the aberrant gene expression in these outliers. Using population-scale mRNA sequencing data, we identify outlier individuals using a multivariate approach. We find that outlier individuals are more readily detected with respect to gene sets that include genes involved in cellular regulation and signal transduction, and less likely to be detected with respect to the gene sets with genes involved in metabolic pathways and other fundamental molecular functions. Analysis of polymorphic data suggests that private SNPs of outlier individuals are enriched in the enhancer and promoter regions of corresponding aberrantly-expressed genes, suggesting a specific regulatory role of private SNPs, while the commonly-occurring regulatory genetic variants (i.e., eQTL SNPs) show little evidence of involvement. Additional data suggest that non-genetic factors may also underlie aberrant gene expression. Taken together, our findings advance a novel viewpoint relevant to situations wherein common eQTLs fail to predict gene expression when heritable, rare inter-individual variation exists. The analytical framework we describe, taking into consideration the reality of differential phenotypic robustness, may be valuable for investigating complex traits and conditions.     

The impact of cis-acting polymorphisms on the human phenotype

Cis-acting polymorphisms that affect gene expression are now known to be frequent, although the extent and mechanisms by which such variation affects the human phenotype are, as yet, only poorly understood. Key signatures of cis-acting variation are differences in gene expression that are tightly associated with regulatory SNPs or expression Quantitative Trait Loci (eQTL) and an imbalance of allelic expression (AEI) in heterozygous samples. Such cis-acting sequence differences appear often to have been under selection within and between populations and are also thought to be important in speciation. Here we describe the example of lactase persistence. In medical research, variants that affect regulation in cis have been implicated in both monogenic and polygenic disorders, and in the metabolism of drugs. In this review we suggest that by further understanding common regulatory variations and how they interact with other genetic and environmental variables it will be possible to gain insight into important mechanisms behind complex disease, with the potential to lead to new methods of diagnosis and treatments.     

Quantitative trait nucleotide analysis using Bayesian model selection

Although much attention has been given to statistical genetic methods for the initial localization and fine mapping of quantitative trait loci (QTLs), little methodological work has been done to date on the problem of statistically identifying the most likely functional polymorphisms using sequence data. In this paper we provide a general statistical genetic framework, called Bayesian quantitative trait nucleotide (BQTN) analysis, for assessing the likely functional status of genetic variants. The approach requires the initial enumeration of all genetic variants in a set of resequenced individuals. These polymorphisms are then typed in a large number of individuals (potentially in families), and marker variation is related to quantitative phenotypic variation using Bayesian model selection and averaging. For each sequence variant a posterior probability of effect is obtained and can be used to prioritize additional molecular functional experiments. An example of this quantitative nucleotide analysis is provided using the GAW12 simulated data. The results show that the BQTN method may be useful for choosing the most likely functional variants within a gene (or set of genes). We also include instructions on how to use our computer program, SOLAR, for association analysis and BQTN analysis.     

Modifier effects between regulatory and protein-coding variation

Genome-wide associations have shown a lot of promise in dissecting the genetics of complex traits in humans with single variants, yet a large fraction of the genetic effects is still unaccounted for. Analyzing genetic interactions between variants (epistasis) is one of the potential ways forward. We investigated the abundance and functional impact of a specific type of epistasis, namely the interaction between regulatory and protein-coding variants. Using genotype and gene expression data from the 210 unrelated individuals of the original four HapMap populations, we have explored the combined effects of regulatory and protein-coding single nucleotide polymorphisms (SNPs). We predict that about 18% (1,502 out of 8,233 nsSNPs) of protein-coding variants are differentially expressed among individuals and demonstrate that regulatory variants can modify the functional effect of a coding variant in cis. Furthermore, we show that such interactions in cis can affect the expression of downstream targets of the gene containing the protein-coding SNP. In this way, a cis interaction between regulatory and protein-coding variants has a trans impact on gene expression. Given the abundance of both types of variants in human populations, we propose that joint consideration of regulatory and protein-coding variants may reveal additional genetic effects underlying complex traits and disease and may shed light on causes of differential penetrance of known disease variants.     

Epistatic selection between coding and regulatory variation in human evolution and disease

Interaction (nonadditive effects) between genetic variants has been highlighted as an important mechanism underlying phenotypic variation, but the discovery of genetic interactions in humans has proved difficult. In this study, we show that the spectrum of variation in the human genome has been shaped by modifier effects of cis-regulatory variation on the functional impact of putatively deleterious protein-coding variants. We analyzed 1000 Genomes population-scale resequencing data from Europe (CEU [Utah residents with Northern and Western European ancestry from the CEPH collection]) and Africa (YRI [Yoruba in Ibadan, Nigeria]) together with gene expression data from arrays and RNA sequencing for the same samples. We observed an underrepresentation of derived putatively functional coding variation on the more highly expressed regulatory haplotype, which suggests stronger purifying selection against deleterious coding variants that have increased penetrance because of their regulatory background. Furthermore, the frequency spectrum and impact size distribution of common regulatory polymorphisms (eQTLs) appear to be shaped in order to minimize the selective disadvantage of having deleterious coding mutations on the more highly expressed haplotype. Interestingly, eQTLs explaining common disease GWAS signals showed an enrichment of putative epistatic effects, suggesting that some disease associations might arise from interactions increasing the penetrance of rare coding variants. In conclusion, our results indicate that regulatory and coding variants often modify the functional impact of each other. This specific type of genetic interaction is detectable from sequencing data in a genome-wide manner, and characterizing these joint effects might help us understand functional mechanisms behind genetic associations to human phenotypes-including both Mendelian and common disease.