Quantitative genetic dissection of complex traits in a QTL-mapping pedigree

This paper summarizes and modifies quantitative genetic analyses on a pedigree used to map genetic factors (i.e., QTLs) underlying a complex trait. The total genetic variance can be exactly estimated based on the F₂ family derived from two homozygous parents for alternative alleles at all QTLs of interest. The parents, F₁ hybrids, and two backcrosses are combined to each parent, and the total number of QTLs and the number of dominant QTLs are estimated under the assumptions of gene association with the two parents, equal gene effect, no linkage, and no epistasis among QTLs. Further relaxation for each of the assumptions are made in detail. The biometric estimator for the QTL number and action mode averaged over the entire genome could provide some basic and complementary information to QTL mapping designed to detect the effect and location of specific genetic factors.     

Broad and narrow heritabilities of quantitative traits in a founder population

Estimation of the components of variance for a quantitative trait allows one to evaluate both the degree to which genetics influences the trait and the trait's underlying genetic architecture. For particular traits, the estimates also may have implications for discriminating between potential models of selection and for choosing an appropriate model for linkage analysis. Using a recently developed method, we estimate the additive and dominance components of variance--or, equivalently, the narrow and broad sense heritabilities--of several traits in the Hutterites, a founder population with extensive genealogical records. As a result of inbreeding and because Hutterite individuals are typically related through multiple lines of descent, we expect that power to detect dominance variance will be increased relative to that in outbred studies. Furthermore, the communal lifestyle of the Hutterites allows us to evaluate the genetic influences in a relatively homogeneous environment. Four phenotypes had a significant dominance variance, resulting in a relatively high broad heritability. We estimated the narrow and broad heritabilities as being, respectively,.36 and.96 for LDL,.51 and 1.0 for serotonin levels, and.45 and.76 for fat free mass (FFM). There was no significant additive component for systolic blood pressure (SBP), resulting in a narrow heritability of 0 and a broad heritability of.45. There were several traits for which we found no significant dominance component, resulting in equal broad and narrow heritability estimates. These traits and their heritabilities are as follows: HDL,.63; triglycerides,.37; diastolic blood pressure,.21; immunoglobulin E,.63; lipoprotein(a),.77; and body-mass index,.54. The large difference between broad and narrow heritabilities for LDL, serotonin, FFM, and SBP are indicative of strong dominance effects in these phenotypes. To our knowledge, this is the first study to report an estimate of heritability for serotonin and to detect a dominance variance for LDL, FFM, and SBP.     

Bay-0 x Shahdara recombinant inbred line population: a powerful tool for the genetic dissection of complex traits in Arabidopsis

Natural genetic variation in Arabidopsis is considerable, but has not yet been used extensively as a source of variants to identify new genes of interest. From the cross between two genetically distant ecotypes, Bay-0 and Shahdara, we generated a Recombinant Inbred Line (RIL) population dedicated to Quantitative Trait Locus (QTL) mapping. A set of 38 physically anchored microsatellite markers was created to construct a robust genetic map from the 420 F6 lines. These markers, evenly distributed throughout the five chromosomes, revealed a remarkable equilibrium in the segregation of parental alleles in the genome. As a model character, we have analysed the genetic basis of variation in flowering time in two different environments. The simultaneous mapping of both large- and small-effect QTLs responsible for this variation explained 90% of the total genotypic variance. Two of the detected QTLs colocalize very precisely with FRIGIDA and FLOWERING LOCUS C genes; we provide information on the polymorphism of genes confirming this hypothesis. Another QTL maps in a region where no QTL had been found previously for this trait. This confirms the accuracy of QTL detection using the Bay-0 x Shahdara RIL population, which constitutes the largest in size available so far in Arabidopsis. As an alternative to mutant analysis, this population represents a powerful tool which is currently being used to undertake the genetic dissection of complex metabolic pathways.     

Re-creation of the genetic composition of a founder population

Human ethnic groups are frequently comprised of two or more founder populations. One of these founding populations is often available for contemporary sampling. We describe a method for reconstructing the composition of a missing founder population using the highly informative haplotypes comprising the HLA system. An application of the method is demonstrated using bone marrow registry samples of African Americans. We use contemporary samples of African Americans and European Americans to derive haplotypes of the West African founder populations. This approach may also be useful for reconstructing ancestral haplotypes for regions elsewhere in the genome.     

Estimation of genetic variability of the founder population in a conservation scheme using microsatellites

In a conservation programme with genealogical records it is possible to estimate the amount of variability of the founder population from a measure of the similarity among the individuals in the current population based on microsatellite markers. Here we compare three available methods and we shown that the one based on the molecular coancestry coefficient should be preferred.     

A new approach to the problem of multiple comparisons in the genetic dissection of complex traits.

Saturated genetic marker maps are being used to map individual genes affecting quantitative traits. Controlling the "experimentwise" type-I error severely lowers power to detect segregating loci. For preliminary genome scans, we propose controlling the "false discovery rate," that is, the expected proportion of true null hypotheses within the class of rejected null hypotheses. Examples are given based on a granddaughter design analysis of dairy cattle and simulated backcross populations. By controlling the false discovery rate, power to detect true effects is not dependent on the number of tests performed. If no detectable genes are segregating, controlling the false discovery rate is equivalent to controlling the experimentwise error rate. If quantitative loci are segregating in the population, statistical power is increased as compared to control of the experimentwise type-I error. The difference between the two criteria increases with the increase in the number of false null hypotheses. The false discovery rate can be controlled at the same level whether the complete genome or only part of it has been analyzed. Additional levels of contrasts, such as multiple traits or pedigrees, can be handled without the necessity of a proportional decrease in the critical test probability.     

Evaluation of LEXF/FXLE rat recombinant inbred strains for genetic dissection of complex traits.

Recombinant inbred (RI) strains are formed from an outcross between two well-characterized inbred stains followed by at least 20 generations of inbreeding. RI strains can be utilized for the analysis of many complex phenotypic traits. The LEXF/FXLE RI strain set consists of 34 RI strains derived by reciprocal crossing of LE/Stm and F344/Stm. Here we report on genetic dissections of complex traits using this RI set and their parental strains. We have developed strain distribution patterns for 232 informative simple sequence length polymorphism markers. The framework map covers the rat genome except for chromosome Y. Seventy-six phenotype parameters, which included physiological and behavioral traits, were examined for these RI lines. Quantitative trait locus (QTL) analysis of these parameters revealed 27 significant and 91 suggestive QTLs, illustrating the potential of this RI resource for the detection of underlying gene functions for various phenotypes. Although this RI set was originally developed to study susceptibility to chemical-induced tumors, it has been shown to be equally powerful for a wide spectrum of traits. The LEXF/FXLE RI strains have been deposited at the National Bio Resource Project for the Rat in Japan and are maintained under specific pathogen-free conditions. They are available at http://www.anim.med.kyoto-u.ac.jp/nbr.     

The future of path analysis, segregation analysis, and combined models for genetic dissection of complex traits.

Appropriate combined models are discussed for the analysis of complex traits. It is argued that combined models may be necessary for optimally extracting the information from family studies. It is further argued that, especially as we face genes with much smaller effects, our ability to find these genes will depend on how precisely and accurately we are able to model the interrelationships. We need these newer models and methods for optimally extracting the information from family data, and we also need to reorient ourselves as to how we interpret the very information extracted. It is projected that path and segregation analysis, as seen in terms of combined models, will be useful in the new millennium.     

Production of congenic mouse strains carrying NOD-derived diabetogenic genetic intervals: An approach for the genetic dissection of complex traits

Insulin-dependent (Type 1) diabetes (IDD) in the NOD mouse is inherited as a complex polygenic trait making the identification of susceptibility genes difficult. Currently none of the non-MHC IDD susceptibility genes in NOD have been identified. In this paper we describe the congenic mouse approach that we are using for the dissection of complex traits, such as IDD. We produced a series of six congenic strains carrying NOD-derived diabetogenic genomic intervals, which were previously identified by linkage analysis, on a resistant background. These congenic strains were produced for the purpose of characterizing the function of each of these genes, alone and in combinations, in IDD pathogenesis and to allow fine mapping of the NOD IDD susceptibility genes. Histological examination of pancreata from 6 to 8-month-old congenic mice reveals that intervals on Chromosomes (Chrs) 1 and 17, but not 3, 6, and 11, contain NOD-derived genes that can increase the trafficking of mononuclear cells into the pancreas. Insulitis was observed only very rarely, even in older congenic mice, indicating that multiple genes are required for this phenotype. These results demonstrate the utility of this congenic approach for the study of complex genetic traits. Received: 1 September 1995 / Accepted: 20 December 1995     

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.     

The split of the Arara population: comparison of genetic drift and founder effect

The total genetic diversity of the Amerindian population is as high as that observed for other continental human populations because a large contribution from variation among tribes makes up for the low variation within tribes. This is attributed mainly to genetic drift acting on small isolated populations. However, a small founder population with a low genetic diversity is another factor that may contribute to the low intratribal diversity. Small founder populations seem to be a frequent event in the formation of new tribes among the Amerindians, but this event is usually not well recorded. In this paper, we analyze the genetic diversity of the Arara of Laranjal village and the Arara of Iriri village, with respect to seven tandem repeat autosomic segments (D1S80, ApoB, D4S43, vW1, vW2, F13A1 and D12S67), two Y-chromosome-specific polymorphisms (DYS19 and DYS199), and mitochondrial DNA (mtDNA) markers (restriction fragment length polymorphisms and sequencing of a segment of the D loop region). The occurrence of a single Y chromosome and mtDNA haplotype, and only 1-4 alleles of the autosomic loci investigated, corroborates historic and demographic records that the Arara of Iriri were founded by a single couple of siblings who came from the Arara of Laranjal, the largest group. Notwithstanding this fact, the genetic distance and the molecular variance between the two Arara villages were greater than those observed between them and other Amazonian tribes, suggesting that the microevolutionary process among Brazilian Amerindians may be misinterpreted if historic demographic data are not considered.     

Toward understanding genetic mechanisms of complex traits in rice

Rice is the primary carbohydrate staple cereal feeding the world population. Many genes, known as quantitative trait loci (QTLs), con-trol most of the agronomically important traits in rice. The identification of QTLs controlling agricultural traits is vital to increase yield and meet the needs of the increasing human population, but the progress met with challenges due to complex QTL inheritance. To date,many QTLs have been detected in rice, including those responsible for yield and grain quality; salt, drought and submergence tolerance;disease and insect resistance; and nutrient utilization efficiency. Map-based cloning techniques have enabled scientists to successfully fine map and clone approximately seventeen QTLs for several traits. Additional in-depth functional analyses and characterizations of these genes will provide valuable assistance in rice molecular breeding.     

Impact of founder population, drift and selection on the genetic diversity of a recently translocated tree population

Recently established, temperate tree populations combine a high level of differentiation for adaptive traits, suggesting rapid genetic evolution, with a high level of genetic diversity within population, suggesting a limited impact of genetic drift and purifying selection. To study experimentally the evolutionary forces in a recently established population, we assessed the spatial and temporal patterns of genetic diversity within a disjunct population of Cedrus atlantica established 140 years ago in south-eastern France from a North African source. The population is expanding through natural regeneration. Three generations were sampled, including founder trees. We analysed 12 isozyme loci, three of which were previously found in tight association with selected genes, and quantitative traits. No bottleneck effect was detected in the founder generation, but a simple test of allelic association revealed an initial disequilibrium which disappeared in the following generations. The impact of genetic drift during secondary evolution was limited, as suggested by the weak temporal differentiation. The genetic load was not reduced after 3 generations, and the quantitative variation for adaptive traits did not change either. Thus, initial genetic changes first proceed from a rapid re-organisation of the diversity through mating and recombination, whereas genetic erosion through drift and selection is delayed due to temporal and spatial stochasticity. Two life-history traits of trees contribute to slowing down the processes of genetic erosion: perenniality and large spatial scale. Thus, one would expect recently established tree populations to have a higher diversity than older ones, which seems in accordance with experimental surveys.     

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.     

Mapping the genetic architecture of complex traits in experimental populations

SUMMARY: Understanding how interactions among set of genes affect diverse phenotypes is having a greater impact on biomedical research, agriculture and evolutionary biology. Mapping and characterizing the isolated effects of single quantitative trait locus (QTL) is a first step, but we also need to assemble networks of QTLs and define non-additive interactions (epistasis) together with a host of potential environmental modulators. In this article, we present a full-QTL model with which to explore the genetic architecture of complex trait in multiple environments. Our model includes the effects of multiple QTLs, epistasis, QTL-by-environment interactions and epistasis-by-environment interactions. A new mapping strategy, including marker interval selection, detection of marker interval interactions and genome scans, is used to evaluate putative locations of multiple QTLs and their interactions. All the mapping procedures are performed in the framework of mixed linear model that are flexible to model environmental factors regardless of fix or random effects being assumed. An F-statistic based on Henderson method III is used for hypothesis tests. This method is less computationally greedy than corresponding likelihood ratio test. In each of the mapping procedures, permutation testing is exploited to control for genome-wide false positive rate, and model selection is used to reduce ghost peaks in F-statistic profile. Parameters of the full-QTL model are estimated using a Bayesian method via Gibbs sampling. Monte Carlo simulations help define the reliability and efficiency of the method. Two real-world phenotypes (BXD mouse olfactory bulb weight data and rice yield data) are used as exemplars to demonstrate our methods. AVAILABILITY: A software package is freely available at http://ibi.zju.edu.cn/software/qtlnetwork