Search

Quantitative trait loci (QTL) hotspots (genomic locations affecting many traits) are a common feature in genetical genomics studies and are biologically interesting since they may harbor critical regulators. Therefore, statistical procedures to assess the significance of hotspots are of key importance. One approach, randomly allocating observed QTL across the genomic locations separately by trait, implicitly assumes all traits are uncorrelated. Recently, an empirical test for QTL hotspots was proposed on the basis of the number of traits that exceed a predetermined LOD value, such as the standard permutation LOD threshold. The permutation null distribution of the maximum number of traits across all genomic locations preserves the correlation structure among the phenotypes, avoiding the detection of spurious hotspots due to nongenetic correlation induced by uncontrolled environmental factors and unmeasured variables. However, by considering only the number of traits above a threshold, without accounting for the magnitude of the LOD scores, relevant information is lost. In particular, biologically interesting hotspots composed of a moderate to small number of traits with strong LOD scores may be neglected as nonsignificant. In this article we propose a quantile-based permutation approach that simultaneously accounts for the number and the LOD scores of traits within the hotspots. By considering a sliding scale of mapping thresholds, our method can assess the statistical significance of both small and large hotspots. Although the proposed approach can be applied to any type of heritable high-volume "omic" data set, we restrict our attention to expression (e)QTL analysis. We assess and compare the performances of these three methods in simulations and we illustrate how our approach can effectively assess the significance of moderate and small hotspots with strong LOD scores in a yeast expression data set.

We extensively phenotyped 6000 Arabidopsis plants with experimentally perturbed DNA methylomes as well as a diverse panel of natural accessions in a common garden. We found that alterations in DNA methylation not only caused heritable phenotypic diversity but also produced heritability patterns closely resembling those of the natural accessions. Our findings indicate that epigenetically induced and naturally occurring variation in complex traits share part of their polygenic architecture and may offer complementary adaptation routes in ecological settings.

Alternative splicing is considered a major mechanism for creating multicellular diversity from a limited repertoire of genes. Here, we performed the first study of genetic variation controlling alternative splicing patterns by comprehensively identifying quantitative trait loci affecting the differential expression of transcript isoforms in a large recombinant inbred population of Caenorhabditis elegans, using a new generation of whole-genome very-high-density oligonucleotide microarrays. Using 60 experimental lines, we were able to detect 435 genes with substantial heritable variation, of which 36% were regulated at a distance (in trans). Nonetheless, we find only a very small number of examples of heritable variation in alternative splicing (22 transcripts), and most of these genes colocalize with the associated genomic loci. Our findings suggest that the regulatory mechanism of alternative splicing in C. elegans is robust toward genetic variation at the genome-wide scale, which is in striking contrast to earlier observations in humans.

Combining global gene-expression profiling and genetic analysis of natural allelic variation (genetical genomics) has great potential in dissecting the genetic pathways underlying complex phenotypes. Efficient use of microarrays is paramount in experimental design as the cost of conducting this type of study is high. For those organisms where recombinant inbred lines are available for mapping, the "distant pair design" maximizes the number of informative contrasts over all marker loci. Here, we describe an extension of this design, named the "optimal pair design," for use with F2 crosses between outbred lines. The performance of this design is investigated by simulation and compared to several other two-color microarray designs. We show that, for a given number of microarrays, the optimal pair design outperforms all other designs considered for detection of expression quantitative trait loci (eQTL) with additive effects by linkage analysis. We also discuss the suitability of this design for outbred crosses in organisms with large genomes and for detection of dominance.

Whole-genome profiling of gene expression in a segregating population has the potential to identify the regulatory consequences of natural allelic variation. Costs of such studies are high and require that resources--microarrays and population--are used as efficiently as possible. We show that current studies can be improved significantly by a new design for two-color microarrays. Our "distant pair design" profiles twice as many individuals as there are arrays, cohybridizes individuals with dissimilar genomes, gives more weight to known regulatory loci if wished, and therewith maximizes the power for decomposing expression variation into regulatory factors. It can also exploit a large population (larger than twice the number of available microarrays) as a useful resource to select the most dissimilar pairs of individuals from. Our approach identifies more regulatory factors than alternative strategies do in computer simulations for realistic genome sizes, and similar promising results are obtained in an application on Arabidopsis thaliana. Our results will aid the design and analysis of future studies on gene expression and will help to shed more light on gene regulatory networks.

Short-oligonucleotide arrays typically contain multiple probes per gene. In genetical genomics applications a statistical model for the individual probe signals can help in separating "true" differential mRNA expression from "ghost" effects caused by polymorphisms, misdesigned probes, and batch effects. It can also help in detecting alternative splicing, start, or termination.

Epistasis is a common and important phenomenon, as indicated by results from a number of recent experiments. Unfortunately, the discovery of epistatic quantitative trait loci (QTL) is difficult since one must search for multiple QTL simultaneously in two or more dimensions. Such a multidimensional search necessitates many statistical tests, and a high statistical threshold must be adopted to avoid false positives. Furthermore, the large number of (interaction) parameters in comparison with the number of observations results in a serious danger of overfitting and overinterpretation of the data. In this article we present a new statistical framework for mapping epistasis in inbred line crosses. It is based on reducing the high dimensionality of the problem in two ways. First, epistatic QTL are mapped in a one-dimensional genome scan for high interactions between QTL and the genetic background. Second, the dimension of the search is bounded by penalized likelihood methods. We use simulated backcross data to illustrate the new approach.

The chromatin loop model predicts that genes within the same chromatin domain exhibit coordinated regulation. We here present the first direct experimental support for this model in plants. Two reporter genes, the E. coli beta-glucuronidase gene and the firefly luciferase gene, driven by different promoters, were placed between copies of the chicken lysozyme A element, a member of the matrix-associated region (MAR) group of chromatin boundary elements, and introduced in tobacco (Nicotiana tabacum). In plants carrying A elements, quantitative enzyme activities and mRNA levels of both genes show high correlations compared to control plants. The A element thus creates an artificial chromatin domain that yields coordinated expression. Surprisingly, enzyme activities correlated poorly with their respective mRNA levels. We hypothesize that this indicates the occurrence of "error pipelines" in data generation: systematic errors of a given analytical method will point in the same direction and cancel out in correlation analysis, resulting in better correlations. In combining different methods of analysis, however, such errors do not cancel out and as a result relevant correlations can be masked. Such error pipelines will have to be taken into account when different types of (e.g., whole-genome) data sets are combined in quantitative analyses.

A mixture model approach is presented for the mapping of one or more quantitative trait loci (QTLs) in complex populations. In order to exploit the full power of complete linkage maps the simultaneous likelihood of phenotype and a multilocus (all markers and putative QTLs) genotype is computed. Maximum likelihood estimation in our mixture models is implemented via an Expectation-Maximization algorithm: exact, stochastic or Monte Carlo EM by using a simple and flexible Gibbs sampler. Parameters include allele frequencies of markers and QTLs, discrete or normal effects of biallelic or multiallelic QTLs, and homogeneous or heterogeneous residual variances. As an illustration a dairy cattle data set consisting of twenty half-sib families has been reanalyzed. We discuss the potential which our and other approaches have for realistic multiple-QTL analyses in complex populations.

Six transgenic tobacco lines, each homozygous for the beta-glucuronidase (GUS) gene at a different locus, and wild type were selfed and intercrossed to evaluate GUS activity in all possible hemizygous, homozygous and dihybrid combinations of GUS alleles. The transgenic lines are characterized by their GUS activity (two low, three intermediate, one high), T-DNA complexity (four single-copy, two more complex single-locus) and the presence of the chicken lysozyme matrix-associated region (MAR) around the full T-DNA (two lines). Gene action and interaction was analyzed by weighted linear regression with parameters for additivity, dominance and epistasis. The analysis showed that each of the four single-copy lines acted fully additively. In contrast, the two complex single-locus lines showed classical single-locus overdominance and were epistatic dominant over all other GUS alleles. The latter is manifested in severe suppression of GUS activity in dihybrid lines, irrespective of the presence of MAR elements around the GUS gene. Such elements apparently do not protect against epistatic dominance. The quantitative data suggested that the epistatic dominance and overdominance are based on the same molecular mechanism. Our approach of a genetic analysis of quantitative variation in well-characterized transgenic lines provides a powerful tool to gain insight into complex plant traits.

In this paper we address the mapping of multiple quantitative trait loci (QTLs) in line crosses for which the genetic data are highly incomplete. Such complicated situations occur, for instance, when dominant markers are used or when unequally informative markers are used in experiments with outbred populations. We describe a general and flexible Monte Carlo expectation-maximization (Monte Carlo EM) algorithm for fitting multiple-QTL models to such data. Implementation of this algorithm is straightforward in standard statistical software, but computation may take much time. The method may be generalized to cope with more complex models for animal and human pedigrees. A practical example is presented, where a three-QTL model is adopted in an outbreeding situation with dominant markers. The example is concerned with the linkage between randomly amplified polymorphic DNA (RAPD) markers and QTLs for partial resistance to Fusarium oxysporum in lily.

Although the interval mapping method is widely used for mapping quantitative trait loci (QTLs), it is not very well suited for mapping multiple QTLs. Here, we present the results of a computer simulation to study the application of exact and approximate models for multiple QTLs. In particular, we focus on an automatic two-stage procedure in which in the first stage "important" markers are selected in multiple regression on markers. In the second stage a QTL is moved along the chromosomes by using the preselected markers as cofactors, except for the markers flanking the interval under study. A refined procedure for cases with large numbers of marker cofactors is described. Our approach will be called MQM mapping, where MQM is an acronym for "multiple-QTL models" as well as for "marker-QTL-marker." Our simulation work demonstrates the great advantage of MQM mapping compared to interval mapping in reducing the chance of a type I error (i.e., a QTL is indicated at a location where actually no QTL is present) and in reducing the chance of a type II error (i.e., a QTL is not detected).

A very general method is described for multiple linear regression of a quantitative phenotype on genotype [putative quantitative trait loci (QTLs) and markers] in segregating generations obtained from line crosses. The method exploits two features, (a) the use of additional parental and F1 data, which fixes the joint QTL effects and the environmental error, and (b) the use of markers as cofactors, which reduces the genetic background noise. As a result, a significant increase of QTL detection power is achieved in comparison with conventional QTL mapping. The core of the method is the completion of any missing genotypic (QTL and marker) observations, which is embedded in a general and simple expectation maximization (EM) algorithm to obtain maximum likelihood estimates of the model parameters. The method is described in detail for the analysis of an F2 generation. Because of the generality of the approach, it is easily applicable to other generations, such as backcross progenies and recombinant inbred lines. An example is presented in which multiple QTLs for plant height in tomato are mapped in an F2 progeny, using additional data from the parents and their F1 progeny.

The interval mapping method is widely used for the mapping of quantitative trait loci (QTLs) in segregating generations derived from crosses between inbred lines. The efficiency of detecting and the accuracy of mapping multiple QTLs by using genetic markers are much increased by employing multiple QTL models instead of the single QTL models (and no QTL models) used in interval mapping. However, the computational work involved with multiple QTL models is considerable when the number of QTLs is large. In this paper it is proposed to combine multiple linear regression methods with conventional interval mapping. This is achieved by fitting one QTL at a time in a given interval and simultaneously using (part of) the markers as cofactors to eliminate the effects of additional QTLs. It is shown that the proposed method combines the easy computation of the single QTL interval mapping method with much of the efficiency and accuracy of multiple QTL models.