3FunMap: full-sib family functional mapping of dynamic traits

MOTIVATION: Functional mapping that embeds the developmental mechanisms of complex traits shows great power to study the dynamic pattern of genetic effects triggered by individual quantitative trait loci (QTLs). A full-sib family, produced by crossing two heterozygous parents, is characteristic of uncertainties about cross-type at a locus and linkage phase between different loci. Integrating functional mapping into a full-sib family requires a model selection procedure capable of addressing these uncertainties. 3FunMap, written in VC++ 6.0, provides a flexible and extensible platform to perform full-sib functional mapping of dynamic traits. Functions in the package encompass linkage phase determination, marker map construction and the pattern identification of QTL segregation, dynamic tests of QTL effects, permutation tests and numerical simulation. We demonstrate the features of 3FunMap through real data analysis and computer simulation. AVAILABILITY: http://statgen.psu.edu/software.     

A non-stationary model for functional mapping of complex traits

SUMMARY: Understanding the genetic control of growth is fundamental to agricultural, evolutionary and biomedical genetic research. In this article, we present a statistical model for mapping quantitative trait loci (QTL) that are responsible for genetic differences in growth trajectories during ontogenetic development. This model is derived within the maximum likelihood context, implemented with the expectation-maximization algorithm. We incorporate mathematical aspects of growth processes to model the mean vector and structured antedependence models to approximate time-dependent covariance matrices for longitudinal traits. Our model has been employed to map QTL that affect body mass growth trajectories in both male and female mice of an F2 population derived from the Large and Small mouse strains. The results from this model are compared with those from the autoregressive-based functional mapping approach. Based on results from computer simulation studies, we suggest that these two models are alternative to one another and should be used simultaneously for the same dataset.     

A fast algorithm for functional mapping of complex traits

By integrating the underlying developmental mechanisms for the phenotypic formation of traits into a mapping framework, functional mapping has emerged as an important statistical approach for mapping complex traits. In this note, we explore the feasibility of using the simplex algorithm as an alternative to solve the mixture-based likelihood for functional mapping of complex traits. The results from the simplex algorithm are consistent with those from the traditional EM algorithm, but the simplex algorithm has considerably reduced computational times. Moreover, because of its nonderivative nature and easy implementation with current software, the simplex algorithm enjoys an advantage over the EM algorithm in the dynamic modeling and analysis of complex traits.     

Bayesian functional mapping of dynamic quantitative traits

Without consideration of other linked QTLs responsible for dynamic trait, original functional mapping based on a single QTL model is not optimal for analyzing multiple dynamic trait loci. Despite that composite functional mapping incorporates the effects of genetic background outside the tested QTL in mapping model, the arbitrary choice of background markers also impact on the power of QTL detection. In this study, we proposed Bayesian functional mapping strategy that can simultaneously identify multiple QTL controlling developmental patterns of dynamic traits over the genome. Our proposed method fits the change of each QTL effect with the time by Legendre polynomial and takes the residual covariance structure into account using the first autoregressive equation. Also, Bayesian shrinkage estimation was employed to estimate the model parameters. Especially, we specify the gamma distribution as the prior for the first-order auto-regressive coefficient, which will guarantee the convergence of Bayesian sampling. Simulations showed that the proposed method could accurately estimate the QTL parameters and had a greater statistical power of QTL detection than the composite functional mapping. A real data analysis of leaf age growth in rice is used for the demonstration of our method. It shows that our Bayesian functional mapping can detect more QTLs as compared to composite functional mapping.     

High-resolution genetic mapping of complex traits.

Positional cloning requires high-resolution genetic mapping. To plan a positional cloning project, one needs to know how many informative meioses will be required to narrow the search for a disease gene to an acceptably small region. For a simple Mendelian trait studied with linkage analysis, the answer is straightforward. In this paper, we address the situation of a complex trait studied with affected-relative-pair methods. We derive mathematical formulas for the size of an appropriate confidence region, as a function of the relative risk attributable to the gene. Using these results, we provide graphs showing the number of relative pairs required to narrow the gene hunt to an interval of a given size. For example, we show that localizing a gene to 1 cM requires a median of 200 sib pairs for a locus causing a fivefold increased risk to an offspring and 700 sib pairs for a locus causing a twofold increased risk. We discuss the implications of these results for the positional cloning of genes underlying complex traits.     

Multivariate eQTL mapping uncovers functional variation on the X-chromosome associated with complex disease traits

Very few studies have investigated the associations between genetic polymorphisms and gene expression on the X-chromosome. This is a major bottleneck when conducting functional follow-up studies of trait-associated variants, as those identified in genome-wide association studies (GWAS). We used a multivariate approach to test the association between individual single nucleotide polymorphisms (SNPs) and exon expression levels measured in 356 Epstein–Barr virus-transformed lymphoblastoid cell lines (LCLs) from the Geuvadis RNA sequencing project to identify SNPs associated with variation in gene expression on the X-chromosome, which we refer to as eSNPs. At an FDR of 5 %, we discovered 548 independent [linkage disequilibrium (LD) r ² < 0.1] eSNPs on the X-chromosome. Of these, 35 were in LD (r ² > 0.8) with previously published disease- or trait-associated variants identified through GWAS. One of the strongest eSNPs identified was rs35975601, which was associated with F8A1 expression (p value = 3 × 10^?20) and was in LD with a type 1 diabetes risk variant. Additionally, we identified a number of genes for which eSNPs were in LD with multiple diseases or traits, including DNASE1L1 which was mapped to bilirubin levels, type 1 diabetes and schizophrenia. Our results also indicate that multivariate exon-level analysis provides a more powerful approach than univariate gene-level analysis, particularly when SNPs influence the expression of different exons with different magnitude and/or direction of effect. The associations identified in our study may provide new insights into the molecular process by which gene expression may contribute to trait variation or disease risk in humans.     

Approaches to mapping genetically correlated complex traits

Our Markov chain Monte Carlo (MCMC) methods were used in linkage analyses of the Framingham Heart Study data using all available pedigrees. Our goal was to detect and map loci associated with covariate-adjusted traits log triglyceride (lnTG) and high-density lipoprotein cholesterol (HDL) using multipoint LOD score analysis, Bayesian oligogenic linkage analysis and identity-by-descent (IBD) scoring methods. Each method used all marker data for all markers on a chromosome. Bayesian linkage analysis detected a linkage signal on chromosome 7 for lnTG and HDL, corroborating previously published results. However, these results were not replicated in a classical linkage analysis of the data or by using IBD scoring methods.We conclude that Bayesian linkage analysis provides a powerful paradigm for mapping trait loci but interpretation of the Bayesian linkage signals is subjective. In the absence of a LOD score method accommodating genetically complex traits and linkage heterogeneity, validation of these signals remains elusive.     

Use of population isolates for mapping complex traits

Geneticists have repeatedly turned to population isolates for mapping and cloning Mendelian disease genes. Population isolates possess many advantages in this regard. Foremost among these is the tendency for affected individuals to share ancestral haplotypes derived from a handful of founders. These haplotype signatures have guided scientists in the fine mapping of scores of rare disease genes. The past successes with Mendelian disorders using population isolates have prompted unprecedented interest among medical researchers in both the public and private sectors. Despite the obvious genetic and environmental complications, geneticists have targeted several population isolates for mapping genes for complex diseases.     

Epistasis,complex traits,and mapping genes

Using a three-locus model wherein two loci regulate a third, candidate locus, I examine physiological epistasis from the gene's eye view of the regulated locus. I show that, depending upon genetic background at the regulatory loci, an allele at the candidate locus can be dominant, additive, recessive, neutral, over-dominant, or under-dominant in its effects on fitness. This kind of variation in allelic effect caused by variation in genetic background from population to population, from time to time in the same population, or sample to sample makes finding and mapping the genes underlying a complex phenotype difficult. The rate of evolution of such genes can also be slowed, especially in genetically subdivided metapopulations with migration. Nevertheless, understanding how variation in genetic background causes variation in allelic effects permits the genetic architecture of such complex traits to be dissected into the interacting component genes. While some backgrounds diminish allelic effects and make finding and mapping genes difficult, other backgrounds enhance allelic effects and facilitate gene mapping.     

Cost-effective designs for linkage disequilibrium mapping of complex traits

The current development of densely spaced collections of single nucleotide polymorphisms (SNPs) will lead to genomewide association studies for a wide range of diseases in many different populations. Determinations of the appropriate number of SNPs to genotype involve a balancing of power and cost. Several variables are important in these determinations. We show that there are different combinations of sample size and marker density that can be expected to achieve the same power. Within certain bounds, investigators can choose between designs with more subjects and fewer markers or those with more markers and fewer subjects. Which designs are more cost-effective depends on the cost of phenotyping versus the cost of genotyping. We show that, under the assumption of a set cost for genotyping, one can calculate a "threshold cost" for phenotyping; when phenotyping costs per subject are less than this threshold, designs with more subjects will be more cost-effective than designs with more markers. This framework for determining a cost-effective study will aid in the planning of studies, especially if there are choices to be made with respect to phenotyping methods or study populations.     

TASSEL: software for association mapping of complex traits in diverse samples

Association analyses that exploit the natural diversity of a genome to map at very high resolutions are becoming increasingly important. In most studies, however, researchers must contend with the confounding effects of both population and family structure. TASSEL (Trait Analysis by aSSociation, Evolution and Linkage) implements general linear model and mixed linear model approaches for controlling population and family structure. For result interpretation, the program allows for linkage disequilibrium statistics to be calculated and visualized graphically. Database browsing and data importation is facilitated by integrated middleware. Other features include analyzing insertions/deletions, calculating diversity statistics, integration of phenotypic and genotypic data, imputing missing data and calculating principal components.     

FLOSS: flexible ordered subset analysis for linkage mapping of complex traits

The FLOSS software package is a flexible framework for ordered subset analysis. FLOSS is specifically designed for use with the Merlin linkage analysis package, but FLOSS can be used with any linkage analysis software package that reports NPL Z-scores for each locus and family. When FLOSS is used with the Merlin linkage analysis package, one can use either non-parametric Z-scores or Kong and Cox linear allele sharing model LOD scores. Monte Carlo P-values are calculated using a permutation test with an efficient Besag-Clifford sequential stopping rule. FLOSS also has a flexible tool for assigning family covariate scores from Merlin input files. FLOSS includes user documentation and is written in Java for easy portability. The FLOSS source code is documented and designed to be extensible.     

Next-generation mapping of complex traits with phenotype-based selection and introgression

Finding the genes underlying complex traits is difficult. We show that new sequencing technology combined with traditional genetic techniques can efficiently identify genetic regions underlying a complex and quantitative behavioral trait. As a proof of concept we used phenotype-based introgression to backcross loci that control innate food preference in Drosophila simulans into the genomic background of D. sechellia, which expresses the opposite preference. We successfully mapped D. simulans introgression regions in a small mapping population (30 flies) with whole-genome resequencing using light coverage (～1×). We found six loci contributing to D. simulans food preference, one of which overlaps a previously discovered allele. This approach is applicable to many systems, does not rely on laborious marker development or genotyping, does not require existing high quality reference genomes, and needs only small mapping populations. Because introgression is used, researchers can scale mapping population size, replication, and number of backcross generations to their needs. Finally, in contrast to more widely used mapping techniques like F(2) bulk-segregant analysis, our method produces near-isogenic lines that can be kept and reused indefinitely.     

Bayesian mapping of quantitative trait loci for complex binary traits

A complex binary trait is a character that has a dichotomous expression but with a polygenic genetic background. Mapping quantitative trait loci (QTL) for such traits is difficult because of the discrete nature and the reduced variation in the phenotypic distribution. Bayesian statistics are proved to be a powerful tool for solving complicated genetic problems, such as multiple QTL with nonadditive effects, and have been successfully applied to QTL mapping for continuous traits. In this study, we show that Bayesian statistics are particularly useful for mapping QTL for complex binary traits. We model the binary trait under the classical threshold model of quantitative genetics. The Bayesian mapping statistics are developed on the basis of the idea of data augmentation. This treatment allows an easy way to generate the value of a hypothetical underlying variable (called the liability) and a threshold, which in turn allow the use of existing Bayesian statistics. The reversible jump Markov chain Monte Carlo algorithm is used to simulate the posterior samples of all unknowns, including the number of QTL, the locations and effects of identified QTL, genotypes of each individual at both the QTL and markers, and eventually the liability of each individual. The Bayesian mapping ends with an estimation of the joint posterior distribution of the number of QTL and the locations and effects of the identified QTL. Utilities of the method are demonstrated using a simulated outbred full-sib family. A computer program written in FORTRAN language is freely available on request.     

Linkage disequilibrium and the mapping of complex human traits.

The potential value of haplotypes defined by several single nucleotide polymorphisms has attracted recent interest. With sufficient linkage disequilibrium (LD), haplotypes could be used in association studies to map common alleles that might influence the susceptibility to common diseases, as well as for reconstructing the evolution of the genome. It has been proposed that a globally useful resource need only be based on high frequency variants, identified from a few modest samples. Rapid progress has been made in quantifying the pattern of human LD and haplotypes defined by such common variants within and among populations. However, the quality and utility of the proposed LD-based resource could be seriously compromised if important sampling and analytical factors are overlooked in its design. The LD map should be based on adequately justified criteria defined by sound population genetic principles.     

QTLNetwork: mapping and visualizing genetic architecture of complex traits in experimental populations

SUMMARY: QTLNetwork is a software package for mapping and visualizing the genetic architecture underlying complex traits for experimental populations derived from a cross between two inbred lines. It can simultaneously map quantitative trait loci (QTL) with individual effects, epistasis and QTL-environment interaction. Currently, it is able to handle data from F(2), backcross, recombinant inbred lines and double-haploid populations, as well as populations from specific mating designs (immortalized F(2) and BC(n)F(n) populations). The Windows version of QTLNetwork was developed with a graphical user interface. Alternatively, the command-line versions have the facility to be run in other prevalent operating systems, such as Linux, Unix and MacOS. AVAILABILITY: http://ibi.zju.edu.cn/software/qtlnetwork.     

Molecular dissection of complex traits in autopolyploids: mapping QTLs affecting sugar yield and related traits in sugarcane

Mapping quantitative trait loci (QTLs) for sugar yield and related traits will provide essential information for sugarcane improvement through marker-assisted selection. Two sugarcane segregating populations derived from interspecific crosses between Saccharum offinarum and Saccharum spontaneum with 264 and 239 individuals, respectively, were evaluated in three replications each for field performance from 1994 to 1996 at Weslaco, Texas. These two populations were analyzed for a total of 735 DNA marker loci to seek QTLs for sugar yield, pol, stalk weight, stalk number, fiber content and ash content. Among the 102 significant associations found between these six traits and DNA markers, 61 could be located on sugarcane linkage maps, while the other 41 were associated with unlinked DNA markers. Fifty of the 61 mapped QTLs were clustered in 12 genomic regions of seven sugarcane homologous groups. Many cases in which QTLs from different genotypes mapped to corresponding locations suggested that at least some of the QTLs on the same cluster might be different allelic forms of the same genes. With a few exceptions that explained part of the transgressive segregation observed for particular traits, the allele effects of most QTLs were consistent with the parental phenotype from which the allele was derived. Plants with a high sugar yield possessed a large number of positive QTLs for sugar yield components and a minimal number of negative QTLs. This indicates the potential effectiveness of marker-assisted selection for sugar yield in sugarcane.     

A nonlinear mixed-effect mixture model for functional mapping of dynamic traits

Functional mapping has emerged as a next-generation statistical tool for mapping quantitative trait loci (QTL) that affect complex dynamic traits. In this article, we incorporated the idea of nonlinear mixed-effect (NLME) models into the mixture-based framework of functional mapping, aimed to generalize the spectrum of applications for functional mapping. NLME-based functional mapping, implemented with the linearization algorithm based on the first-order Taylor expansion, can provide reasonable estimates of QTL genotypic-specific curve parameters (fixed effect) and the between-individual variation of these parameters (random effect). Results from simulation studies suggest that the NLME-based model is more general than traditional functional mapping. The new model can be useful for the identification of the ontogenetic patterns of QTL genetic effects during time course.