Two-trait-locus linkage analysis: a powerful strategy for mapping complex genetic traits.

Recent advances in molecular biology have provided geneticists with ever-increasing numbers of highly polymorphic genetic markers that have made possible linkage mapping of loci responsible for many human diseases. However, nearly all diseases mapped to date follow clear Mendelian, single-locus segregation patterns. In contrast, many common familial diseases such as diabetes, psoriasis, several forms of cancer, and schizophrenia are familial and appear to have a genetic component but do not exhibit simple Mendelian transmission. More complex models are required to explain the genetics of these important diseases. In this paper, we explore two-trait-locus, two-marker-locus linkage analysis in which two trait loci are mapped simultaneously to separate genetic markers. We compare the utility of this approach to standard one-trait-locus, one-marker-locus linkage analysis with and without allowance for heterogeneity. We also compare the utility of the two-trait-locus, two-marker-locus analysis to two-trait-locus, one-marker-locus linkage analysis. For common diseases, pedigrees are often bilineal, with disease genes entering via two or more unrelated pedigree members. Since such pedigrees often are avoided in linkage studies, we also investigate the relative information content of unilineal and bilineal pedigrees. For the dominant-or-recessive and threshold models that we consider, we find that two-trait-locus, two-marker-locus linkage analysis can provide substantially more linkage information, as measured by expected maximum lod score, than standard one-trait-locus, one-marker-locus methods, even allowing for heterogeneity, while, for a dominant-or-dominant generating model, one-locus models that allow for heterogeneity extract essentially as much information as the two-trait-locus methods. For these three models, we also find that bilineal pedigrees provide sufficient linkage information to warrant their inclusion in such studies. We also discuss strategies for assessing the significance of the two linkages assumed in two-trait-locus, two-marker-locus models.     

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.     

Microsatellite-AFLP for genetic mapping of complex polyploids.

In spite of the economical relevance of polyploid crops, genetic mapping of these species has been relatively overlooked. This is because of intrinsic difficulties such as the uncertainty of the chromosome behavior at meiosis I and the need for very large segregating populations. An important, yet underestimated issue, in mapping polyploids is the choice of the molecular marker system. An ideal molecular marker system for polyploid mapping should maximize the percentage of single dose markers (SDMs) detected and the possibility of recognizing allelic markers. In the present work, the marker index for genetic mapping (MIgm) of M-AFLP is compared with that of AFLP and SAMPL. M-AFLPs have the highest MIgm values (22 vs. 18.5 of SAMPL and 9.83 of AFLP) mostly because of their high power to detect polymorphism. Owing to their prevalent codominant inheritance, it is proposed that M-AFLP can be used for the preliminary identification of hom(e)ologous groups.     

Strategies for genetic mapping of categorical traits

The search for efficient and powerful statistical methods and optimal mapping strategies for categorical traits under various experimental designs continues to be one of the main tasks in genetic mapping studies. Methodologies for genetic mapping of categorical traits can generally be classified into two groups, linear and non-linear models. We develop a method based on a threshold model, termed mixture threshold model to handle ordinal (or binary) data from multiple families. Monte Carlo simulations are done to compare its statistical efficiencies and properties of the proposed non-linear model with a linear model for genetic mapping of categorical traits using multiple families. The mixture threshold model has notably higher statistical power than linear models. There may be an optimal sampling strategy (family size vs number of families) in which genetic mapping reaches its maximal power and minimal estimation errors. A single large-sibship family does not necessarily produce the maximal power for detection of quantitative trait loci (QTL) due to genetic sampling of QTL alleles. The QTL allelic model has a marked impact on efficiency of genetic mapping of categorical traits in terms of statistical power and QTL parameter estimation. Compared with a fixed number of QTL alleles (two or four), the model with an infinite number of QTL alleles and normally distributed allelic effects results in loss of statistical power. The results imply that inbred designs (e.g. F₂ or four-way crosses) with a few QTL alleles segregating or reducing number of QTL alleles (e.g. by selection) in outbred populations are desirable in genetic mapping of categorical traits using data from multiple families. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

FunMap: functional mapping of complex traits

SUMMARY: FunMap is a Web-based user interface designed to map quantitative trait loci (QTL) affecting function-valued traits or infinite-dimensional traits in well-structured pedigrees or natural populations. User input includes three files: longitudinal trait data, marker genotypes and/or a linkage map. This software allows for a systematic genome-wide scan and significance test of QTL throughout the map. The dynamic change of QTL effects during the time course of growth is automatically drawn, from which specific biological hypotheses regarding the genetic control mechanisms of growth and development can be tested. AVAILABILITY: http://web.biostat.ufl.edu/~cma/genetics/software.html     

Prospects for admixture mapping of complex traits

Admixture mapping extends to human populations the principles that underlie linkage analysis of an experimental cross. For detecting genes that contribute to ethnic variation in disease risk, admixture mapping has greater statistical power than family-linkage studies. In comparison with association studies, admixture mapping requires far fewer markers to search the genome and is less affected by allelic heterogeneity. Statistical-analysis programs for admixture mapping are now available, and a genomewide panel of markers for admixture mapping in populations formed by West African-European admixture has been assembled. Some of the remaining technical challenges include the ability to ensure that the statistical methods are robust and to develop marker panels for other admixed populations. Where admixed populations and panels of markers informative for ancestry are available, admixture mapping can be applied to localize genes that contribute to ethnic variation in any measurable trait.     

Genotype matrix mapping: searching for quantitative trait loci interactions in genetic variation in complex traits.

In order to reveal quantitative trait loci (QTL) interactions and the relationship between various interactions in complex traits, we have developed a new QTL mapping approach, named genotype matrix mapping (GMM), which searches for QTL interactions in genetic variation. The central approach in GMM is the following. (1) Each tested marker is given a virtual matrix, named a genotype matrix (GM), containing intersecting lines and rows equal to the total allele number for that marker in the population analyzed. (2) QTL interactions are then estimated and compared through virtual networks among the GMs. To evaluate the contribution of marker combinations to a quantitative phenotype, the GMM method divides the samples into two non-overlapping subclasses, S(0) and S(1); the former contains the samples that have a specific genotype pattern to be evaluated, and the latter contains samples that do not. Based on this division, the F-measure is calculated as an index of significance. With the GMM method, we extracted significant marker combinations consisting of one to three interacting markers. The results indicated there were multiple QTL interactions affecting the phenotype (flowering date). GMM will be a valuable approach to identify QTL interactions in genetic variation of a complex trait within a variety of organisms.     

Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits

SUMMARY: A website for performing power calculations for the design of linkage and association genetic mapping studies of complex traits. AVAILABILITY: The package is made available athttp://statgen.iop.kcl.ac.uk/gpc/.     

High-resolution genetic mapping using the Mouse Diversity outbred population

The JAX Diversity Outbred population is a new mouse resource derived from partially inbred Collaborative Cross strains and maintained by randomized outcrossing. As such, it segregates the same allelic variants as the Collaborative Cross but embeds these in a distinct population architecture in which each animal has a high degree of heterozygosity and carries a unique combination of alleles. Phenotypic diversity is striking and often divergent from phenotypes seen in the founder strains of the Collaborative Cross. Allele frequencies and recombination density in early generations of Diversity Outbred mice are consistent with expectations based on simulations of the mating design. We describe analytical methods for genetic mapping using this resource and demonstrate the power and high mapping resolution achieved with this population by mapping a serum cholesterol trait to a 2-Mb region on chromosome 3 containing only 11 genes. Analysis of the estimated allele effects in conjunction with complete genome sequence data of the founder strains reduced the pool of candidate polymorphisms to seven SNPs, five of which are located in an intergenic region upstream of the Foxo1 gene.     

A dissection model for mapping complex traits

Many quantitative traits are composites of other traits that contribute differentially to genetic variation. Quantitative trait locus (QTL) mapping of these composite traits can benefit by incorporating the mechanistic process of how their formation is mediated by the underlying components. We propose a dissection model by which to map these interconnected components traits under a joint likelihood setting. The model can test how a composite trait is determined by pleiotropic QTLs for its component traits or jointly by different sets of QTLs each responsible for a different component. The model can visualize the pattern of time‐varying genetic effects for individual components and their impacts on composite traits. The dissection model was used to map two composite traits, stemwood volume growth decomposed into its stem height, stem diameter and stem form components for Euramerican poplar adult trees, and total lateral root length constituted by its average lateral root length and lateral root number components for Euphrates poplar seedlings. We found the pattern of how QTLs for different components contribute to phenotypic variation in composite traits. The detailed understanding of the genetic machineries of composite traits will not only help in the design of molecular breeding in plants and animals, but also shed light on the evolutionary processes of quantitative traits under natural selection.     

Functional mapping - how to map and study the genetic architecture of dynamic complex traits

The development of any organism is a complex dynamic process that is controlled by a network of genes as well as by environmental factors. Traditional mapping approaches for analysing phenotypic data measured at a single time point are too simple to reveal the genetic control of developmental processes. A general statistical mapping framework, called functional mapping, has been proposed to characterize, in a single step, the quantitative trait loci (QTLs) or nucleotides (QTNs) that underlie a complex dynamic trait. Functional mapping estimates mathematical parameters that describe the developmental mechanisms of trait formation and expression for each QTL or QTN. The approach provides a useful quantitative and testable framework for assessing the interplay between gene actions or interactions and developmental changes.     

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.     

High-resolution genetic mapping of mammalian motor activity levels in mice

The generation of motor activity levels is under tight neural control to execute essential behaviors, such as movement toward food or for social interaction. To identify novel neurobiological mechanisms underlying motor activity levels, we studied a panel of chromosome substitution (CS) strains derived from mice with high (C57BL/6J strain) or low motor activity levels (A/J strain) using automated home cage behavioral registration. In this study, we genetically mapped the expression of baseline motor activity levels (horizontal distance moved) to mouse chromosome 1. Further genetic mapping of this trait revealed an 8.3-Mb quantitative trait locus (QTL) interval. This locus is distinct from the QTL interval for open-field anxiety-related motor behavior on this chromosome. By data mining, an existing phenotypic and genotypic data set of 2445 genetically heterogeneous mice ( http://gscan.well.ox.ac.uk/ ), we confirmed linkage to the peak marker at 79 970 253 bp and refined the QTL to a 312-kb interval containing a single gene ( A830043J08Rik ). Sequence analysis showed a nucleotide deletion in the 3' untranslated region of the Riken gene. Genome-wide microarray gene expression profiling in brains of discordant F₂ individuals from CS strain 1 showed a significant upregulation of Epha4 in low-active F₂ individuals. Inclusion of a genetic marker for Epha4 confirmed that this gene is located outside of the QTL interval. Both Epha4 and A830043J08Rik are expressed in brain motor circuits, and similar to Epha4 mutants, we found linkage between reduced motor neurons number and A/J chromosome 1. Our findings provide a novel QTL and a potential downstream target underlying motor circuitry development and the expression of physical activity levels.     

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.     

Combining information within and between pedigrees for mapping complex traits.

This paper is concerned with efficient strategies for gene mapping using pedigrees containing small numbers of affecteds and identity-by-descent data from closely spaced markers throughout the genome. Particular attention is paid to additive traits involving phenocopies and/or locus heterogeneity. For a sample of pedigrees containing a particular configuration of affecteds, e.g., pairs of siblings together with a first cousin, we use a likelihood analysis to find 1-df statistics that are very efficient over a broad range of penetrances and allele frequencies. We identify configurations of affecteds that are particularly powerful for detecting linkage, and we show how pedigrees containing different numbers and configurations of affecteds can be efficiently combined in an overall test statistic.     

High-resolution mapping and genetic characterization of the Lazy-2 gravitropic mutant of tomato.

Mutation of the Lazy-2 (Lz-2) gene in tomato (Lycopersicon esculentum mill.) produces a phytochrome-dependent reversal of shoot gravitropism, providing a unique genetic resource for investigating how signals from light modulate gravitropism. We mapped the Lz-2 gene using RFLPs and a PCR-based technique to assess the feasibility of positional cloning. Analysis of a 1338 plant backcross population between L. esculentum and L. pennellii placed Lz-2 within a 1.2 cM interval on chromosome 5, 0.4 cM from TG504-CT201A interval. The inabililty to resolve these markers indicates that Lz-2 resides in a centromeric region in which recombination is highly suppressed. Lazy-2 is tightly linked to but does not encode the gene for ACC4, an enzyme involved in ethylene biosynthesis. We also observed that Lz-2 is partially dominant under certain conditions and stages of development.     

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.