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.     

Mapping quantitative trait loci for carcass and meat quality traits in a wild boar x Large White intercross

An intercross between wild boar and a domestic Large White pig population was used to map quantitative trait loci (QTL) for body proportions, weight of internal organs, carcass composition, and meat quality. The results concerning growth traits and fat deposition traits have been reported elsewhere. In the present study, all 200 F2 animals, their parents, and their grandparents were genotyped for 236 markers. The marker genotypes were used to calculate the additive and dominance coefficients at fixed positions in the genome of each F2 animal, and the trait values were regressed onto these coefficients in intervals of 1 cM. In addition, the effect of proportion of wild boar alleles was tested for each chromosome. Significant QTL effects were found for percentage lean meat and percentage lean meat plus bone in various cuts, proportion of bone in relation to lean meat in ham, muscle area, and carcass length. The significant QTL were located on chromosomes 2, 3, 4, and 8. Each QTL explained 9 to 16% of the residual variance of the traits. Gene action for most QTL was largely additive. For meat quality traits, there were no QTL that reached the significance threshold. However, the average proportion of wild boar alleles across the genome had highly significant effects on reflectance and drip loss. The results show that there are several chromosome regions with a considerable effect on carcass traits in pigs.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Estimation of effects of quantitative trait loci in large complex pedigrees

A method was derived to estimate effects of quantitative trait loci (QTL) using incomplete genotype information in large outbreeding populations with complex pedigrees. The method accounts for background genes by estimating polygenic effects. The basic equations used are very similar to the usual linear mixed model equations for polygenic models, and segregation analysis was used to estimate the probabilities of the QTL genotypes for each animal. Method R was used to estimate the polygenic heritability simultaneously with the QTL effects. Also, initial allele frequencies were estimated. The method was tested in a simulated data set of 10,000 animals evenly distributed over 10 generations, where 0, 400 or 10,000 animals were genotyped for a candidate gene. In the absence of selection, the bias of the QTL estimates was < 2%. Selection biased the estimate of the Aa genotype slightly, when zero animals were genotyped. Estimates of the polygenic heritability were 0.251 and 0.257, in absence and presence of selection, respectively, while the simulated value was 0.25. Although not tested in this study, marker information could be accommodated by adjusting the transmission probabilities of the genotypes from parent to offspring according to the marker information. This renders a QTL mapping study in large multi-generation pedigrees possible.     

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.     

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.     

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.     

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.     

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.     

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.     

Quantitative trait loci for murine growth

Body size is an archetypal quantitative trait with variation due to the segregation of many gene loci, each of relatively minor effect, and the environment. We examine the effects of quantitative trait loci (QTLs) on age-specific body weights and growth in the F2 intercross of the LG/J and SM/J strains of inbred mice. Weekly weights (1-10 wk) and 75 microsatellite genotypes were obtained for 535 mice. Interval mapping was used to locate and measure the genotypic effects of QTLs on body weight and growth. QTL effects were detected on 16 of the 19 autosomes with several chromosomes carrying more than one QTL. The number of QTLs for age-specific weights varied from seven at 1 week to 17 at 10 wk. The QTLs were each of relatively minor, subequal effect. QTLs affecting early and late growth were generally distinct, mapping to different chromosomal locations indicating separate genetic and physiological systems for early and later murine growth.     

Common quantitative trait loci for alcohol-related behaviors and CNS neurotensin measures: voluntary ethanol consumption

The alcohol breath test (ABT) is evaluated for variability in response to changes in physiological parameters. The ABT was originally developed in the 1950s, at a time when understanding of pulmonary physiology was quite limited. Over the past decade, physiological studies have shown that alcohol is exchanged entirely within the conducting airways via diffusion from the bronchial circulation. This is in sharp contrast to the old idea that alcohol exchanges in the alveoli in a manner similar to the lower solubility respiratory gases (O2 and CO2). The airway alcohol exchange process is diffusion (airway tissue) and perfusion (bronchial circulation) limited. The dynamics of airway alcohol exchange results in a positively sloped exhaled alveolar plateau that contributes to considerable breathing pattern-dependent variation in measured breath alcohol concentration measurements.