Segregating variation in the transcriptome: cis regulation and additivity of effects

Independent effects of cis- and trans-regulatory variation on gene expression in Drosophila melanogaster

Biochemical interactions between cis-regulatory DNA sequences and trans-regulatory gene products suggest that cis- and trans-acting polymorphisms may interact genetically. Here we present a strategy to test this hypothesis by comparing the relative cis-regulatory activity of two alleles in different genetic backgrounds. Of the eight genes surveyed in this study, five were affected by trans-acting variation that altered total transcript levels, two of which were also affected by differences in cis-regulation. The presence of trans-acting variation had no effect on relative cis-regulatory activity, showing that cis-regulatory polymorphisms can function independently of trans-regulatory variation. The frequency of such independent interactions on a genomic scale is yet to be determined.

Molecular ecology of global change

Global environmental change is altering the selection regime for all biota. The key selective factors are altered mean, variance and seasonality of climatic variables and increase in CO(2) concentration itself. We review recent studies that document rapid evolution to global climate change at the phenotypic and genetic level, as a response to shifts in these factors. Among the traits that have changed are photoperiod responses, stress tolerance and traits associated with enhanced dispersal. The genetic basis of two traits with a critical role under climate change, stress tolerance and photoperiod behaviour, is beginning to be understood for model organisms, providing a starting point for candidate gene approaches in targeted nonmodel species. Most studies that have documented evolutionary change are correlative, while selection experiments that manipulate relevant variables are rare. The latter are particularly valuable for prediction because they provide insight into heritable change to simulated future conditions. An important gap is that experimental selection regimes have mostly been testing one variable at a time, while synergistic interactions are likely under global change. The expanding toolbox available to molecular ecologists holds great promise for identifying the genetic basis of many more traits relevant to fitness under global change. Such knowledge, in turn, will significantly advance predictions on global change effects because presence and polymorphism of critical genes can be directly assessed. Moreover, knowledge of the genetic architecture of trait correlations will provide the necessary framework for understanding limits to phenotypic evolution; in particular as lack of critical gene polymorphism or entire pathways, metabolic costs of tolerance and linkage or pleiotropy causing negative trait correlations. Synergism among stressor impacts on organismal function may be causally related to conflict among transcriptomic syndromes specific to stressor types. Because adaptation to changing environment is always contingent upon the spatial distribution of genetic variation, high-resolution estimates of gene flow and hybridization should be used to inform predictions of evolutionary rates.

Gene expression quantitative trait locus analysis of 16 000 barley genes reveals a complex pattern of genome-wide transcriptional regulation

Transcript abundance from cRNA hybridizations to Affymetrix microarrays can be used for simultaneous marker development and genome-wide gene expression quantitative trait locus (eQTL) analysis of crops. We have previously shown that it is easily possible to use Affymetrix expression arrays to profile individuals from a segregating population to accurately identify robust polymorphic molecular genetic markers. We applied the method to identify more than 2000 genetic polymorphisms (transcript derived markers, TDMs) from an experiment involving two commercial varieties of barley (Hordeum vulgare; Steptoe and Morex) and their doubled-haploid progeny. With this set of TDMs, we constructed a genetic map and used it for the genome-wide eQTL analysis of about 16 000 genes in a relatively large population (n = 139). We identified 23 738 significant eQTLs at a genome-wide significance (P </= 0.05), affecting the expression of 12 987 genes. Over a third of these genes with expression variation have only one identified eQTL while the rest have two to six. A large proportion of the quantitatively controlled transcripts appear to be regulated by both cis and trans effects. More than half of the quantitatively controlled transcripts appear to be primarily regulated by cis eQTLs in this population. We show that although there appear to be eQTL hotspots many of these are in chromosomal regions of low recombination, such as genetic centromeres, and so have a high gene density per centimorgan. Some chromosomal regions have a significant excess of eQTL over the number expected from gene density, and many of these are biased towards eQTL for which the allele from one particular parent increases the expression level.

Extensive additivity of gene expression differentiates subspecies of the house mouse

We have studied different subspecies of the house mouse and their reciprocal F(1) hybrids to estimate the within-locus mode of inheritance for subspecies differences in gene expression in three tissues (brain, liver, and testis) of male mice. This study investigates the mode of inheritance in crosses at a larger taxonomic distance than has been previously systematically investigated. We found the vast majority of transcripts to be additively expressed with only a few transcripts showing dominance or overdominance in expression, except for one direction of one cross, which showed large mis-expression in the testis. We suggest that, as time passes, more genes come to influence expression, and if there is no directional dominance, additivity becomes increasingly more likely, up to a threshold beyond which there is F(1) hybrid breakdown. Some previous studies on different organisms have found a large degree of dominance, commonly at shorter taxonomic differences. We surveyed these findings and show that the most consistent association exists between the amount of additivity detected in a study and the expression analysis method (in particular microarray platform), suggesting that at least some of the differences among studies might be methodological. Most studies agree with ours in that within-locus additivity seems to be general mode of inheritance for transcript expression. Differentially expressed transcripts identified in our screen among subspecies of house mice are candidate genes that could be involved in reproductive isolation between these subspecies.

Cis sequence effects on gene expression

Background

Sequence and transcriptional variability within and between individuals are typically studied independently. The joint analysis of sequence and gene expression variation (genetical genomics) provides insight into the role of linked sequence variation in the regulation of gene expression. We investigated the role of sequence variation in cis on gene expression (cis sequence effects) in a group of genes commonly studied in cancer research in lymphoblastoid cell lines. We estimated the proportion of genes exhibiting cis sequence effects and the proportion of gene expression variation explained by cis sequence effects using three different analytical approaches, and compared our results to the literature.

Results

We generated gene expression profiling data at N = 697 candidate genes from N = 30 lymphoblastoid cell lines for this study and used available candidate gene resequencing data at N = 552 candidate genes to identify N = 30 candidate genes with sufficient variance in both datasets for the investigation of cis sequence effects. We used two additive models and the haplotype phylogeny scanning approach of Templeton (Tree Scanning) to evaluate association between individual SNPs, all SNPs at a gene, and diplotypes, with log-transformed gene expression. SNPs and diplotypes at eight candidate genes exhibited statistically significant (p < 0.05) association with gene expression. Using the literature as a "gold standard" to compare 14 genes with data from both this study and the literature, we observed 80% and 85% concordance for genes exhibiting and not exhibiting significant cis sequence effects in our study, respectively.

Conclusion

Based on analysis of our results and the extant literature, one in four genes exhibits significant cis sequence effects, and for these genes, about 30% of gene expression variation is accounted for by cis sequence variation. Despite diverse experimental approaches, the presence or absence of significant cis sequence effects is largely supported by previously published studies.

Allele-specific expression patterns reveal biases and embryo-specific parent-of-origin effects in hybrid maize

We employed allele-specific expression (ASE) analyses to document biased allelic expression in maize (Zea mays). A set of 316 quantitative ASE assays were used to profile the relative allelic expression in seedling tissue derived from five maize hybrids. The different hybrids included in this study exhibit a range of heterosis levels; however, we did not observe differences in the frequencies of allelic bias. Allelic biases in gene expression were consistently observed for approximately 50% of the genes assayed in hybrid seedlings. The relative proportion of genes that exhibit cis- or trans-acting regulatory variation was very similar among the different genotypes. The cis-acting regulatory variation was more prevalent and resulted in greater expression differences than trans-acting regulatory variation for these genes. The ASE assays were further used to compare the relative expression of the B73 and Mo17 alleles in three tissue types (seedling, immature ear, and embryo) derived from reciprocal hybrids. These comparisons provided evidence for tissue-specific cis-acting variation and for a slight maternal expression bias in approximately 20% of genes in embryo tissue. Collectively, these data provide evidence for prevalent cis-acting regulatory variation that contributes to biased allelic expression between genotypes and between tissues.

Expression evolution in yeast genes of single-input modules is mainly due to changes in trans-acting factors

Both cis- and trans-regulatory mutations contribute to gene expression divergence within and between species. To estimate their relative contributions, we examined two yeast strains, BY (a laboratory strain) and RM (a wild strain), for their gene-expression divergence by microarray. Using these data and published ChIP-chip data, we obtained a set of single-regulator-regulated genes that showed expression divergence between BY and RM. We randomly selected 50 of these genes for further study. We developed a step-by-step approach to assess the relative contributions of cis- and trans-variations to expression divergence by using pyrosequencing to quantify the mRNA levels of the BY and RM alleles in the same culture (co-culture) and in hybrid diploids. Forty genes showed expression divergence between the two strains in co-culture, and pyrosequencing of the BY/RM hybrid diploids showed that 45% (18/40) can be attributed to differences in trans-acting factors alone, 17.5% (7/40) mainly to trans-variations, 20% (8/40) to both cis- and trans-acting factors, 7.5% (3/40) mainly to cis-variations, and 10% (4/40) to cis-acting factors alone. In addition, we replaced the BY promoter by the RM promoter in each of 10 BY genes that were found from our microarray data to have expression divergence between BY and RM, and in each case our quantitative PCR analysis revealed a cis effect of the promoter replacement on gene expression. In summary, our study suggests that trans-acting factors play the major role in expression evolution between yeast strains, but the role of cis variation is also important.

Global eQTL mapping reveals the complex genetic architecture of transcript-level variation in Arabidopsis

The genetic architecture of transcript-level variation is largely unknown. The genetic determinants of transcript-level variation were characterized in a recombinant inbred line (RIL) population (n = 211) of Arabidopsis thaliana using whole-genome microarray analysis and expression quantitative trait loci (eQTL) mapping of transcript levels as expression traits (e-traits). Genetic control of transcription was highly complex: one-third of the quantitatively controlled transcripts/e-traits were regulated by cis-eQTL, and many trans-eQTL mapped to hotspots that regulated hundreds to thousands of e-traits. Several thousand eQTL of large phenotypic effect were detected, but almost all (93%) of the 36,871 eQTL were associated with small phenotypic effects (R(2) < 0.3). Many transcripts/e-traits were controlled by multiple eQTL with opposite allelic effects and exhibited higher heritability in the RILs than their parents, suggesting nonadditive genetic variation. To our knowledge, this is the first large-scale global eQTL study in a relatively large plant mapping population. It reveals that the genetic control of transcript level is highly variable and multifaceted and that this complexity may be a general characteristic of eukaryotes.