A review of statistical methods for expression quantitative trait loci mapping

With high-throughput technologies now widely available, investigators can easily measure thousands of phenotypes for quantitative trait loci (QTL) mapping. Microarray measurements are particularly amenable to QTL mapping, as evidenced by a number of recent studies demonstrating utility across a broad range of biological endeavors. The early success stories have impelled a rapid increase in both the number and complexity of expression QTL (eQTL) experiments. Consequently, there is a need to consider the statistical principles involved in the design and analysis of these experiments and the methods currently being used. In this article we review these principles and methods and discuss the open questions most likely to yield significant progress toward increasing the amount of meaningful information obtained from eQTL mapping experiments.     

Ghost QTL and hotspots in experimental crosses: novel approach for modeling polygenic effects

Ghost quantitative trait loci (QTL) are the false discoveries in QTL mapping, that arise due to the “accumulation” of the polygenic effects, uniformly distributed over the genome. The locations on the chromosome that are strongly correlated with the total of the polygenic effects depend on a specific sample correlation structure determined by the genotypes at all loci. The problem is particularly severe when the same genotypes are used to study multiple QTL, e.g. using recombinant inbred lines or studying the expression QTL. In this case, the ghost QTL phenomenon can lead to false hotspots, where multiple QTL show apparent linkage to the same locus. We illustrate the problem using the classic backcross design and suggest that it can be solved by the application of the extended mixed effect model, where the random effects are allowed to have a nonzero mean. We provide formulas for estimating the thresholds for the corresponding t-test statistics and use them in the stepwise selection strategy, which allows for a simultaneous detection of several QTL. Extensive simulation studies illustrate that our approach eliminates ghost QTL/false hotspots, while preserving a high power of true QTL detection.     

Identification of quantitative trait loci influencing foliar concentrations of terpenes and formylated phloroglucinol compounds in Eucalyptus nitens

Leaves of eucalypt species contain a variety of plant secondary metabolites, including terpenoids and formylated phloroglucinol compounds (FPCs). Both terpene and FPC concentrations are quantitative traits that can show large variation within a population and have been shown to be heritable. The molecular genetic basis of this variation is currently unknown. Progeny from a field trial of a three-generation mapping pedigree of Eucalyptus nitens were assayed for terpenes and FPCs. Quantitative trait loci (QTL) analyses were conducted using a map constructed from 296 markers to locate regions of the genome influencing foliar concentrations of these plant secondary compounds. A large number of significant QTL for 14 traits were located across nine linkage groups, with significant clustering of QTL on linkage groups 7, 8 and 9. As expected, QTL for biosynthetically related compounds commonly colocated, but QTL for unrelated monterpenes and FPCs also mapped closely together. Colocation of these QTL with mapped candidate genes from the various biosynthetic pathways, and subsequent use of these genes in association mapping, will assist in determining the causes of variation in plant secondary metabolites in eucalypts.     

关联分析及其在植物遗传学研究中的应用

植物的很多重要经济性状均属于复杂性状。基于连锁分析的QTL作图是研究复杂性状的有效手段,但其尚存在一定的局限性。随着现代生物学的发展,一种基于连锁不平衡的新剖分复杂性状方法——关联分析法,开始应用于植物遗传学研究。与QTL作图法相比,应用关联分析法具有不需要构建特殊的群体,可同时对多个等位基因进行分析,定位QTL精度可达到单基因水平等优势。该文介绍了关联分析方法学的基础和特性,简述了其在植物遗传学研究中的进展情况,并对其未来发展和在植物遗传学研究中的应用进行了展望。     

QTL mapping and the genetic basis of adaptation: recent developments

Quantitative trait loci (QTL) mapping has been used in a number of evolutionary studies to study the genetic basis of adaptation by mapping individual QTL that explain the differences between differentiated populations and also estimating their effects and interaction in the mapping population. This analysis can provide clues about the evolutionary history of populations and causes of the population differentiation. QTL mapping analysis methods and associated computer programs provide us tools for such an inference on the genetic basis and architecture of quantitative trait variation in a mapping population. Current methods have the capability to separate and localize multiple QTL and estimate their effects and interaction on a quantitative trait. More recent methods have been targeted to provide a comprehensive inference on the overall genetic architecture of multiple traits in a number of environments. This development is important for evolutionary studies on the genetic basis of multiple trait variation, genotype by environment interaction, host–parasite interaction, and also microarray gene expression QTL analysis.     

An integrative approach for the identification of quantitative trait loci

The genetic dissection of complex traits is one of the most difficult and most important challenges facing science today. We discuss here an integrative approach to quantitative trait loci (QTL) mapping in mice. This approach makes use of the wealth of genetic tools available in mice, as well as the recent advances in genome sequence data already available for a number of inbred mouse strains. We have developed mapping strategies that allow a stepwise narrowing of a QTL mapping interval, prioritizing candidate genes for further analysis with the potential of identifying the most probable candidate gene for the given trait. This approach integrates traditional mapping tools, fine mapping tools, sequence-based analysis, bioinformatics and gene expression.     

Deciphering the genomic architecture of the stickleback brain with a novel multilocus gene‐mapping approach

Quantitative traits important to organismal function and fitness, such as brain size, are presumably controlled by many small‐effect loci. Deciphering the genetic architecture of such traits with traditional quantitative trait locus (QTL) mapping methods is challenging. Here, we investigated the genetic architecture of brain size (and the size of five different brain parts) in nine‐spined sticklebacks (Pungitius pungitius) with the aid of novel multilocus QTL‐mapping approaches based on a de‐biased LASSO method. Apart from having more statistical power to detect QTL and reduced rate of false positives than conventional QTL‐mapping approaches, the developed methods can handle large marker panels and provide estimates of genomic heritability. Single‐locus analyses of an F₂ interpopulation cross with 239 individuals and 15 198, fully informative single nucleotide polymorphisms (SNPs) uncovered 79 QTL associated with variation in stickleback brain size traits. Many of these loci were in strong linkage disequilibrium (LD) with each other, and consequently, a multilocus mapping of individual SNPs, accounting for LD structure in the data, recovered only four significant QTL. However, a multilocus mapping of SNPs grouped by linkage group (LG) identified 14 LGs (1–6 depending on the trait) that influence variation in brain traits. For instance, 17.6% of the variation in relative brain size was explainable by cumulative effects of SNPs distributed over six LGs, whereas 42% of the variation was accounted for by all 21 LGs. Hence, the results suggest that variation in stickleback brain traits is influenced by many small‐effect loci. Apart from suggesting moderately heritable (h² ≈ 0.15–0.42) multifactorial genetic architecture of brain traits, the results highlight the challenges in identifying the loci contributing to variation in quantitative traits. Nevertheless, the results demonstrate that the novel QTL‐mapping approach developed here has distinctive advantages over the traditional QTL‐mapping methods in analyses of dense marker panels.     

家禽数量性状基因座定位的研究进展

近二十年来,各种DNA标记技术及相关生物技术的发展和完善为高密度、覆盖面广的连锁图谱的构建及QTL的定位奠定了基础,本就目前世界上建立的几个较有影响的资源家系、各种DNA标记技术、家禽中定位的QTL及存在问题等方面作一综述。     

QUANTITATIVE GENETIC ANALYSIS OF DIVERGENCE IN MALE SECONDARY SEXUAL TRAITS BETWEEN DROSOPHILA SIMULANS AND DROSOPHILA MAURITIANA

The sibling species Drosophila simulans and D. mauritiana differ significantly in a number of male secondary sexual traits, providing an ideal system for genetic analysis of interspecific morphological divergence. In the experiment reported here, F₁ hybrids from a cross of two inbred lines were backcrossed in both directions and about 200 flies from each backcross were scored for several traits (bristle numbers and cuticle areas), as well as 18 markers distributed throughout the genome. Each trait was analyzed by composite interval mapping to identify quantitative trait loci (QTL) and estimate their effects. For each trait, from one to eight loci were detected, with more divergent traits showing evidence for greater numbers of QTL. Estimates of additive effects varied widely, with a range of 0.4 to 4.1 environmental standard deviation units and an average of 2.2 units. There was substantial evidence for nonadditive effects, since the magnitude of estimates often differed significantly between the two backcrosses. The sign of the estimated effect differed among QTL for bristle traits, but not for cuticle area traits, suggesting that these two types of trait may have undergone different types of selection. Finally, several similarities were found between different traits in the estimated positions of QTL, suggesting that pleiotropy and/or linkage of QTL may have been important in the evolution of these traits.     

Quantile-based permutation thresholds for quantitative trait Loci hotspots

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.     

关联分析及其在植物遗传学研究中的应用

植物的很多重要经济性状均属于复杂性状。基于连锁分析的QTL作图是研究复杂性状的有效手段, 但其尚存在一定的局限性。随着现代生物学的发展, 一种基于连锁不平衡的新剖分复杂性状方法--关联分析法, 开始应用于植物遗传学研究。与QTL作图法相比, 应用关联分析法具有不需要构建特殊的群体, 可同时对多个等位基因进行分析, 定位QTL精度可达到单基因水平等优势。该文介绍了关联分析方法学的基础和特性, 简述了其在植物遗传学研究中的进展情况, 并对其未来发展和在植物遗传学研究中的应用进行了展望。     

Regulatory hotspots are associated with plant gene expression under varying soil phosphorus supply in Brassica rapa

Gene expression is a quantitative trait that can be mapped genetically in structured populations to identify expression quantitative trait loci (eQTL). Genes and regulatory networks underlying complex traits can subsequently be inferred. Using a recently released genome sequence, we have defined cis- and trans-eQTL and their environmental response to low phosphorus (P) availability within a complex plant genome and found hotspots of trans-eQTL within the genome. Interval mapping, using P supply as a covariate, revealed 18,876 eQTL. trans-eQTL hotspots occurred on chromosomes A06 and A01 within Brassica rapa; these were enriched with P metabolism-related Gene Ontology terms (A06) as well as chloroplast- and photosynthesis-related terms (A01). We have also attributed heritability components to measures of gene expression across environments, allowing the identification of novel gene expression markers and gene expression changes associated with low P availability. Informative gene expression markers were used to map eQTL and P use efficiency-related QTL. Genes responsive to P supply had large environmental and heritable variance components. Regulatory loci and genes associated with P use efficiency identified through eQTL analysis are potential targets for further characterization and may have potential for crop improvement.     

Quantitative trait loci mapping of phenotypic plasticity and genotype-environment interactions in plant and insect performance

Community genetic studies generally ignore the plasticity of the functional traits through which the effect is passed from individuals to the associated community. However, the ability of organisms to be phenotypically plastic allows them to rapidly adapt to changing environments and plasticity is commonly observed across all taxa. Owing to the fitness benefits of phenotypic plasticity, evolutionary biologists are interested in its genetic basis, which could explain how phenotypic plasticity is involved in the evolution of species interactions. Two current ideas exist: (i) phenotypic plasticity is caused by environmentally sensitive loci associated with a phenotype; (ii) phenotypic plasticity is caused by regulatory genes that simply influence the plasticity of a phenotype. Here, we designed a quantitative trait loci (QTL) mapping experiment to locate QTL on the barley genome associated with barley performance when the environment varies in the presence of aphids, and the composition of the rhizosphere. We simultaneously mapped aphid performance across variable rhizosphere environments. We mapped main effects, QTL × environment interaction (QTL×E), and phenotypic plasticity (measured as the difference in mean trait values) for barley and aphid performance onto the barley genome using an interval mapping procedure. We found that QTL associated with phenotypic plasticity were co-located with main effect QTL and QTL×E. We also located phenotypic plasticity QTL that were located separately from main effect QTL. These results support both of the current ideas of how phenotypic plasticity is genetically based and provide an initial insight into the functional genetic basis of how phenotypically plastic traits may still be important sources of community genetic effects.