Multiple-trait quantitative trait loci analysis using a large mouse sibship

Quantitative trait loci influencing several phenotypes were assessed using a genetically heterogeneous mouse population. The 145 individuals were produced by a cross between (BALB/cJ x C57BL/6J)F1 females and (C3H/HeJ x DBA/2J)F1 males. The population is genetically equivalent to full siblings derived from heterozygous parents, with known linkage phase. Each individual in the population represents a unique combination of alleles from the inbred grandparents. Quantitative phenotypes for eight T cell measures were obtained at 8 and 18 mo of age. Single-marker locus, repeated measures analysis of variance identified nine marker-phenotype associations with an experimentwise significance level of P < 0.05. Six of the eight quantitative phenotypes could be associated with at least one locus having experiment-wide significance. Composite interval, repeated measures analysis of variance identified 13 chromosomal regions with comparisonwise (nominal) significance associations of P < 0.001. The heterozygous-parent cross provides a reproducible, general method for identification of loci associated with quantitative trait phenotypes or repeated phenotypic measures.     

Multivariate multipoint linkage analysis of quantitative trait loci

We studied the nervous systems and tumors of two patients with anti-Yo-associated paraneoplastic cerebellar degeneration (PCD). In both patients the underlying tumor was an ovarian adenocarcinoma that expressed Yo antigens and contained extensive infiltrates of lymphocytes and plasma cells. The major central nervous system findings were a complete loss of cerebellar Purkinje cells with Bergmann astrogliosis. One patient had inflammatory infiltrates in the medulla and pons, and moderate axonal loss and demyelination involving the spinal cord. No inflammatory infiltrates were identified in the cerebrum, cerebellum or brain-stem of the other patient. Using quantitative Western blot analysis, deposits of anti-Yo IgG could not be demonstrated in the nervous system, possibly as a result of the loss of cells expressing Yo antigens. The detection of the anti-Yo antibody as a common marker of PCD in one patient with inflammatory infiltrates and another without infiltrates suggests that some PCD pathologically classified as "non-inflammatory" may represent a final burn-out stage of a cellular immune-mediated disorder. Our findings indicate that Purkinje cells are the main, but not necessarily the exclusive, targets of this disorder.     

Iteratively reweighted least squares mapping of quantitative trait loci

Mapping quantitative trait loci (QTL) is a typical problem of regression with uncertain independent variables because the genotype of a putative QTL is not observed. Rather, the genotype is inferred from marker information. The method of maximum likelihood (ML) methods is considered to be the optimal solution for this problem because the distribution of the unobserved QTL genotype is fully taken into account. The simple linear regression method (REG) is a first-order approximation to ML and usually performs very well. In this study, an iteratively reweighted least squares method (IRWLS) is proposed. The new method is a second-order approximation to ML because both the expectation and the variance of the unobserved QTL genotype are taken into consideration. The IRWLS is developed in the context of a single large outbred family. The properties of IRWLS are demonstrated and compared with REG and ML via replicated Monte Carlo simulations. The conclusions are: (1) when marker information content is high, the three methods perform equally well, but ML and IRWLS outperform REG when marker information content is low and the variance explained by the QTL is high; (2) when the residual distribution is not normal, ML can fail or have low power to detect small QTLs, but REG and IRWLS are robust to non-normality; and (3) when the residual distribution is normal, the performance of IRWLS is almost identical to ML, but the computational speed of IRWLS is many times faster than that of ML.     

Mapping quantitative trait loci for immune capacity in the pig

Immune capacity traits show considerable genetic variation in outbred populations. To identify quantitative trait loci (QTLs) for immune capacity in the pig, various measures of immune function (total and differential leukocyte counts, neutrophil phagocytosis, mitogen-induced proliferation, IL-2 production, and virus induced IFN-alpha production in whole blood cultures, and Ab responses to two Escherichia coli antigens) were determined in 200 F2 animals from a wild pig-Swedish Yorkshire intercross. The pedigree has been typed for 236 genetic markers covering all autosomes, the X chromosome and the X/Y pseudoautosomal region. Through interval mapping using a least-squares method, four QTLs with significant effects were identified; one for total leukocyte counts, one for mitogen-induced proliferation, one for prevaccination levels of Abs to E. coli Ag K88, and one for Ab response to the O149 Ag. In addition, several putative QTLs were indicated. The results from the present study conclusively show that it is possible to identify QTLs for immune capacity traits in outbred pig populations by genome analysis.     

Mapping quantitative trait loci using multiple families of line crosses

To avoid a loss in statistical power as a result of homozygous individuals being selected as parents of a mapping population, one can use multiple families of line crosses for quantitative trait genetic linkage analysis. Two strategies of combining data are investigated: the fixed-model and the random-model strategies. The fixed-model approach estimates and tests the average effect of gene substitution for each parent, while the random-model approach treats each effect of gene substitution as a random variable and directly estimates and tests the variance of gene substitution. Extensive Monte Carlo simulations verify that the two strategies perform equally well, although the random model is preferable in combining data from a large number of families. Simulations also show that there may be an optimal sampling strategy (number of families vs. number of individuals per family) in which QTL mapping reaches its maximum power and minimum estimation error. Deviation from the optimal strategy reduces the efficiency of the method.     

Extreme discordant sib pairs for mapping quantitative trait loci in humans

Analysis of differences between siblings (sib pair analysis) is a standard method of genetic linkage analysis for mapping quantitative trait loci, such as those contributing to hypertension and obesity, in humans. In traditional designs, pairs are selected at random or with one sib having an extreme trait value. The majority of such pairs provide little power to detect linkage; only pairs that are concordant for high values, low values, or extremely discordant pairs (for example, one in the top 10 percent and the other in the bottom 10 percent of the distribution) provide substantial power. Focus on discordant pairs can reduce the amount of genotyping necessary over conventional designs by 10- to 40-fold.     

Interacting quantitative trait loci control phenotypic variation in murine estradiol-regulated responses

The steroid hormone estradiol (E2) elicits a spectrum of systemic and uterotropic responses in vivo. For example, E2 treatment of ovariectomized adult and sexually immature rodents leads to uterine leukocytic infiltration, cell proliferation, and organ growth. E2-regulated growth is also associated with a variety of normal and pathological phenotypes. Historically, the uterine growth response has been used as the key model to understand the molecular and biochemical mechanisms underlying E2-dependent growth. In this study, genome exclusion mapping identified two quantitative trait loci (QTL) in the mouse, Est2 and Est3 on chromosomes 5 and 11, respectively, that control the phenotypic variation in uterine wet weight. Both QTL are linked to a variety of E2-regulated genes, suggesting that they may represent loci within conserved gene complexes that play fundamental roles in mediating the effects of E2. Interaction and multiple trait analyses using the uterine leukocyte response and wet weight suggest that Est4, a QTL on chromosome 10, may encode an interacting factor that influences the quantitative variation in both responses. Our results show that E2-dependent responses can be genetically controlled and that a genetic basis may underlie the variation observed in many E2-dependent phenotypes.     

Using multivariate genetic modeling to detect pleiotropic quantitative trait loci

The gastric mucosa has been regarded as an active site of humoral immunity since the discovery of Helicobacter pylori. The present study was conducted to determine the in vivo activity of gastric B cells in 53 gastric cancer patients. B-cell activity was measured by protein-A plaque assay, in which IgA-, IgM-, and IgG-plaque-forming cells (PFC) were counted. The number of PFC was associated with the stage of cancer, but the response of lymphocytes in a non-tumorous area (NML) and tumor-infiltrating lymphocytes (TIL) differed. PFC in both sites were decreased compared to n0 cancer in n1 lymph node metastasis-positive cancer, while only NML showed raised PFC in n2 + (P < 0.05, vs TIL). Cancer cells penetrating the submucosa caused the PFC of TIL (but not of NML) to decrease. Invasion of the intratumor capillary (V) or lymphatic (Ly) vessels also caused PFC to change, showing differences of Ig class; there was a decrease of PFC in V2 (IgG- and IgM-PFC) and in Ly2 (all Ig-PFC). IgA-PFC in Ly1 differed in TIL (decrease of PFC) and NML (increase). PFC also differed in TIL and NML in cancer cells, as follows: TIL < NML in tubular and poorly differentiated adenocarcinoma and TIL > NML in papillary and signet ring cell adenocarcinoma. Changes in lymph node (LNL) and blood lymphocytes were similar to those in gastric PFC whose IgA value was 10 times as much as that of LNL. The 5-year survival rate was significantly better in patients with lower rather than higher PFC such as 89% vs 68%. Gastric B cells thus appear to be active and to reflect gastric mucosal immunity.     

Advances in statistical methods to map quantitative trait loci in outbred populations

Statistical methods to map quantitative trait loci (QTL) in outbred populations are reviewed, extensions and applications to human and plant genetic data are indicated, and areas for further research are identified. Simple and computationally inexpensive methods include (multiple) linear regression of phenotype on marker genotypes and regression of squared phenotypic differences among relative pairs on estimated proportions of identity-by-descent at a locus. These methods are less suited for genetic parameter estimation in outbred populations but allow the determination of test statistic distributions via simulation or data permutation; however, further inferences including confidence intervals of QTL location require the use of Monte Carlo or bootstrap sampling techniques. A method which is intermediate in computational requirements is residual maximum likelihood (REML) with a covariance matrix of random QTL effects conditional on information from multiple linked markers. Testing for the number of QTLs on a chromosome is difficult in a classical framework. The computationally most demanding methods are maximum likelihood and Bayesian analysis, which take account of the distribution of multilocus marker-QTL genotypes on a pedigree and permit investigators to fit different models of variation at the QTL. The Bayesian analysis includes the number of QTLs on a chromosome as an unknown.     

Detection of quantitative trait loci in outbred populations with incomplete marker data

AcrA protein is a component of the multi-drug efflux complex AcrAB-TolC of Escherichia coli. Judged by the hypersusceptibility phenotype of acrA mutants, the AcrAB-TolC system pumps out an extraordinarily wide variety of antibiotics, chemotherapeutic agents, detergents and dyes. This complex traverses both the inner and outer membranes of E. coli and catalyzes efflux of the drugs directly into the medium. The coordinated operation of the inner membrane transporter AcrB and outer membrane channel TolC is thought to be mediated by AcrA. The latter is a lipoprotein located in the periplasmic space. We show here that a lipid-deficient derivative of AcrA is functionally active as demonstrated by the complementation of the hypersusceptibility phenotype of the acrA mutant. Purified non-lipidated and intact forms of AcrA were able to restore, with similar efficiency, the activity of AcrA-dependent efflux of erythromycin in Ca2+-sucrose-treated E. coli cells. Using analytical ultracentrifugation and dynamic light scattering techniques we determined hydrodynamic properties of the non-lipidated AcrA and found that AcrA exists in solution as a highly asymmetric monomeric molecule with an axial ratio of 8. This elongated shape of AcrA is compatible with the hypothesis that this protein spans the periplasmic space coordinating the concerted operation of inner and outer membrane components of the complex.     

A two-stage half-sib design for mapping quantitative trait loci in food animals

Daughter and granddaughter half-sib designs for mapping quantitative trait loci were modified to increase experimental power. This new design includes a two-stage procedure, in contrast to conventional one-step half-sib designs. In stage 1, a few progeny of each sire are genotyped for marker loci. Based on the analyses of stage 1 data, some sires are chosen to continue genotyping more progeny for stage 2. When multiple chromosomes are under investigation, chromosomes and sires for stage 2 are selected based on the analysis of stage 1 data. Sire selection results in increased frequency of heterozygous genotypes of interest in stage 2 if the markers are linked to those genes. Chromosome selection can increase the proportion of chromosomes with segregating quantitative trait loci in stage 2 if not all of the chromosomes evaluated in stage 1 have segregating quantitative trait loci. Numerical results indicated that two-stage half-sib designs are generally more powerful than conventional designs when 1) the noncentrality parameter is moderate or larger, 2) larger quantitative trait loci are mapped using tightly linked markers in larger families, and 3) variation is large in numbers and sizes of segregating quantitative trait loci among the chromosomes evaluated in stage 1.     

A general Monte Carlo method for mapping multiple quantitative trait loci

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.     

Confirmation of quantitative trait loci for ethanol sensitivity in long-sleep and short-sleep mice

PURPOSE: To determine the efficacy of the combination of cisplatin, fluorouracil, and high-dose l-leucovorin (PFL) as organ-preserving induction therapy followed by radiotherapy in untreated patients with advanced squamous cell carcinoma of the head and neck. PATIENTS AND METHODS: This was a phase II study of PFL in 47 patients with resectable stage III (n = 20) and IV (n = 27) M0 squamous cell carcinoma of the head and neck, including larynx (n = 20), hypopharynx (n = 14), and oropharynx (n = 13). The PFL regimen consisted of cisplatin 25 mg/m2 on days 1 through 5, fluorouracil 800 mg/m2 CI on days 2 through 6, and l-leucovorin 250 mg/m2 on days 1 through 6, all by continuous intravenous infusion every 21 to 28 days for three courses. The primary study endpoint was initial response to and local disease control rate with PFL as induction chemotherapy, with an aim to confirm the previously reported complete response rate of 60% to 70%. RESULTS: Of 47 patients enrolled, 46 were evaluable for response to PFL, 14 (30%) achieved a complete response, and 25 (54%) achieved a partial response, for an overall response rate of 84%. Of 39 patients evaluable for response after radiation therapy, 27 (69%) achieved a complete response and 11 (28%) a partial response. Local disease control was achieved in 37 of 46 (80%). Grade 3 or 4 toxic effects occurred frequently, with neutropenia in 27 (59%) of 46 evaluable patients, thrombocytopenia in 30%, mucositis in 41%, diarrhea in 13%, and nausea/ vomiting in 13%, but there were no treatment-related deaths. With a median follow-up of 35 months there have been nine recurrences (four local/regional and five distant) and 17 deaths (12 in patients with disease progression and five not directly related to the primary tumor). Second primary tumors have developed in six patients. At 3 years 62% of the patients remain alive with no disease progression, and the 3-year survival estimate with preserved organ function is 66%. CONCLUSION: PFL induction chemotherapy produced only a modest complete response rate, possibly due to suboptimal dose intensity, and was associated with substantial, although not life-threatening, toxicity. Newer regimens and treatment modalities are still needed in the management of advanced squamous cell carcinoma of the head and neck.     

Multiple-trait mapping of quantitative trait loci after selective genotyping using logistic regression

Experiments to map QTL usually measure several traits, and not uncommonly genotype only those animals that are extreme for some trait(s). Analysis of selectively genotyped, multiple-trait data presents special problems, and most simple methods lead to biased estimates of the QTL effects. The use of logistic regression to estimate QTL effects is described, where the genotype is treated as the dependent variable and the phenotype as the independent variable. In this way selection on phenotype does not bias the results. If normally distributed errors are assumed, the logistic-regression analysis is almost equivalent to a maximum-likelihood analysis, but can be carried out with standard statistical packages. Analysis of a simulated half-sib experiment shows that logistic regression can estimate the effect and position of a QTL without bias and confirms the increased power achieved by multiple-trait analysis.     

Statistical methods for mapping quantitative trait loci from a dense set of markers

Lander and Botstein introduced statistical methods for searching an entire genome for quantitative trait loci (QTL) in experimental organisms, with emphasis on a backcross design and QTL having only additive effects. We extend their results to intercross and other designs, and we compare the power of the resulting test as a function of the magnitude of the additive and dominance effects, the sample size and intermarker distances. We also compare three methods for constructing confidence regions for a QTL: likelihood regions, Bayesian credible sets, and support regions. We show that with an appropriate evaluation of the coverage probability a support region is approximately a confidence region, and we provide a theroretical explanation of the empirical observation that the size of the support region is proportional to the sample size, not the square root of the sample size, as one might expect from standard statistical theory.     

Selection with recurrent backcrossing to develop congenic lines for quantitative trait loci analysis

SEWALL WRIGHT suggested that genes of large effect on a quantitative trait could be isolated by recurrent backcrossing with selection on the trait. Loci [quantitative trait loci (QTL)] at which the recurrent and nonrecurrent lines have genes of different large effect on the trait would remain segregating, while other loci would become fixed for the gene carried by the recurrent parent. If the recurrent line is inbred and the backcrossing and selection is conducted in a series of replicate lines, in each of which only one backcross parent is selected for each generation, the lines will become congenic to the recurrent parent except for the QTL of large effect and closely linked regions of the genome, and these regions can be identified using a dense set of markers that differ between the parental lines. Such lines would be particularly valuable for subsequent fine-scale mapping and gene cloning; but by chance, even QTL of large effect will be lost from some lines. The probability that QTL of specified effect remain segregating is computed as a function of its effect on the trait, the intensity of selection, and the number of generations of backcrossing. Analytical formulas are given for one or two loci, and simulation is used for more. It is shown that the method could have substantial discriminating ability and thus potential practical value.     

Identification of quantitative trait loci influencing wood specific gravity in an outbred pedigree of loblolly pine

We report the identification of quantitative trait loci (QTL) influencing wood specific gravity (WSG) in an outbred pedigree of loblolly pine (Pinus taeda L.). QTL mapping in an outcrossing species is complicated by the presence of multiple alleles (> 2) at QTL and marker loci. Multiple alleles at QTL allow the examination of interaction among alleles at QTL (deviation from additive gene action). Restriction fragment length polymorphism (RFLP) marker genotypes and wood specific gravity phenotypes were determined for 177 progeny. Two RFLP linkage maps were constructed, representing maternal and paternal parent gamete segregations as inferred from diploid progeny RFLP genotypes. RFLP loci segregating for multiple alleles were vital for aligning the two maps. Each RFLP locus was assayed for cosegregation with WSG QTL using analysis of variance (ANOVA). Five regions of the genome contained one or more RFLP loci showing differences in mean WSG at or below the P = 0.05 level for progeny as grouped by RFLP genotype. One region contained a marker locus (S6a) whose QTL-associated effects were highly significant (P > 0.0002). Marker S6a segregated for multiple alleles, a prerequisite for determining the number of alleles segregating at the linked QTL and analyzing the interactions among QTL alleles. The QTL associated with marker S6a appeared to be segregating for multiple alleles which interacted with each other and with environments. No evidence for digenic epistasis was found among the five QTL.     

Least squares interval mapping of quantitative trait loci under the infinitesimal genetic model in outbred populations

Genetic marker and phenotypic data for a quantitative trait were simulated on 20 paternal half-sib families with 100 progeny to investigate properties of within-family-regression interval mapping of a postulated single quantitative trait locus (QTL) in a marker interval under the infinitesimal genetic model, which has been the basis of the application of quantitative genetics to genetic improvement programs, and to investigate use of the infinitesimal model as null hypothesis in testing for presence of a major QTL. Genetic effects on the marked chromosome were generated based on a major gene model, which simulated a central biallelic QTL, or based on 101 biallelic QTL of equal effect, which approximated the infinitesimal model. The marked chromosome contained 0, 3.3%, 13.3%, or 33.3% of genetic variance and heritability was 0.25 or 0.70. Under the polygenic model with 3.3% of genetic variance on the marked chromosome, which corresponds to the infinitesimal model for the bovine, significant QTL effects were found for individual families. Correlations between estimates of QTL effects and true chromosome substitution effects were 0.29 and 0.47 for heritabilities of 0.25 and 0.70 but up to 0.85 with 33.3% of polygenic variance on the marked chromosome. These results illustrate the potential of marker-assisted selection even under the infinitesimal genetic model. Power of tests for presence of QTL was substantially reduced when the polygenic model with 3.3% of genetic variance on the chromosome was used as a null hypothesis. The ability to determine whether genetic variance on a chromosome was contributed by a single QTL of major effect or a large number of QTL with minor effects, corresponding to the infinitesimal model, was limited.     

Polygenic mutation in Drosophila melanogaster: genetic interactions between selection lines and candidate quantitative trait loci

We have investigated genetic interactions between spontaneous mutations affecting abdominal and sternopleural bristle number that have accumulated in 12 long-term selection lines derived from an inbred strain, and mutations at 14 candidate bristle number quantitative trait loci. The quantitative test for complementation was to cross the selection lines to an inbred wild-type strain (the control cross) and to a derivative of the control strain into which the mutant allele at the candidate locus to be tested was substituted (the tester strain). Genetic interactions between spontaneous mutations affecting bristle number and the candidate locus mutations were common, and in several cases the interaction effects were different in males and females. Analyses of variance of the (tester- control) differences among and within groups of replicate lines selected in the same direction for the same trait showed significant group effects for several candidate loci. Genetically, the interactions could be caused by allelism of, and/ or epistasis between, spontaneous mutations in the selection lines and the candidate locus mutations. It is possible that much of the response to selection was from new mutations at candidate bristle number quantitative trait loci, and that for some of these loci, mutation rates were high.     

Mapping quantitative trait loci with dominant and missing markers in various crosses from two inbred lines

Bacterial resistance to beta-lactam antibiotics, a clinically worrying and recurrent problem, is often due to the production of beta-lactamases, enzymes that efficiently hydrolyze the amide bond of the beta-lactam nucleus. Imipenem and other carbapenems escape the activity of most active site serine beta-lactamases and have therefore become very popular drugs for antibacterial chemotherapy in the hospital environment. Their usefulness is, however, threatened by the appearance of new beta-lactamases that efficiently hydrolyze them. This study is focused on the structure and properties of two recently described class A carbapenemases, produced by Serratia marcescens and Enterobacter cloacae strains and leads to a better understanding of the specificity of beta-lactamases. In turn, this will contribute to the design of better antibacterial drugs. Three-dimensional models of the two class A carbapenemases were constructed by homology modeling. They suggested the presence, near the active site of the enzymes, of a disulfide bridge (C69-C238) whose existence was experimentally confirmed. Kinetic parameters were measured with the purified Sme-1 carbapenemase, and an attempt was made to explain its specific substrate profile by analyzing the structures of minimized Henri-Michaelis complexes and comparing them to those obtained for the "classical" TEM-1 beta-lactamase. The peculiar substrate profile of the carbapenemases appears to be strongly correlated with the presence of the disulfide bridge between C69 and C238.