iDArTs: increasing the value of genomic resources at no cost

For many years, genetic markers have been the building blocks in assembling genomic knowledge. Improved technology and methods for collecting marker data have increased accuracy, increased throughput, and reduced cost. However, common genotyping technology still produces far fewer markers in plant species than in animals and humans. We propose a new type of genetic marker based on the Diversity Arrays Technology (DArT) genotyping system for organisms lacking a reference genetic sequence. These markers are based directly on microarray probe intensity profiles and hence are called iDArTs. They require no additional genotyping beyond screening with a DArT array. Since standard methods of genetic analysis cannot be used with these continuous markers, we develop novel methods for the common bi-parental experimental designs doubled haploids, recombinant inbred lines, and backcrosses. These enable the augmentation of genetic maps with iDArTs and permit quantitative trait locus mapping with both discrete and continuous markers. We use simulation to demonstrate the power of this approach for marker mapping. In addition, we construct maps and perform linkage analysis for these DArT genotypes using the doubled haploid progeny lines from a cross between the wheat cultivars Chara and Glenlea. These methods allow access to a previously untapped genetic resource by extracting additional information from the raw data. With no additional genotyping cost, we are able to double the number of markers mapped and thereby increase genome coverage.     

Identification of RAPD markers linked to the loci controlling erucic acid level in rapeseed

The recent development of the industrial use of rapeseed oil rich in erucic acid has led to increased interest in the improvement of the high-erucic-acid (50–60%) varieties and to research towards genotypes containing a very high erucic acid content. This trait is controlled by two genes with additive effects. The low-erucic-acid trait was relatively easily introduced through backcrosses into various backgrounds because the zero-erucic-acid homozygotes were clearly identified in the segregating populations. To select for high erucic acid level is more difficult because of the partial overlap of the high-erucic-acid homozygous class and the intermediate one, containing heterozygotes. In order to help conventional breeding, RAPD markers were used to map the two genes involved in determining the erucic acid content in a doubled haploid progeny derived from a low x high erucic acid F1 hybrid. The two genes were successfully localized in two independent linkage group, through a QTL approach. A close association was found between individual plant genotypes and the erucic acid content of the doubled haploid progeny, and it was shown that the two genes do not contribute uniformly to the C22:1 level. The value of molecular gene mapping of such a trait in a conventional breeding programme is discussed.Abbreviations BSA bulked segregant analysis - DH doubled haploid - NIL near-isogenic lines - QTL quantitative trait locus - C22:1 erucic acid - TAG triacyl glycerol - SCAR sequence characterized amplified region     

Determining the sample size for co-dominant molecular marker-assisted linkage detection for a monogenic qualitative trait by controlling the type-I and type-II errors in a segregating F2 population

Tests for linkage are usually performed using the lod score method. A critical question in linkage analyses is the choice of sample size. The appropriate sample size depends on the desired type-I error and power of the test. This paper investigates the exact type-I error and power of the lod score method in a segregating F(2) population with co-dominant markers and a qualitative monogenic dominant-recessive trait. For illustration, a disease-resistance trait is considered, where the susceptible allele is recessive. A procedure is suggested for finding the appropriate sample size. It is shown that recessive plants have about twice the information content of dominant plants, so the former should be preferred for linkage detection. In some cases the exact alpha-values for a given nominal alpha may be rather small due to the discrete nature of the sampling distribution in small samples. We show that a gain in power is possible by using exact methods.     

Direct power comparisons between simple LOD scores and NPL scores for linkage analysis in complex diseases.

Several methods have been proposed for linkage analysis of complex traits with unknown mode of inheritance. These methods include the LOD score maximized over disease models (MMLS) and the "nonparametric" linkage (NPL) statistic. In previous work, we evaluated the increase of type I error when maximizing over two or more genetic models, and we compared the power of MMLS to detect linkage, in a number of complex modes of inheritance, with analysis assuming the true model. In the present study, we compare MMLS and NPL directly. We simulated 100 data sets with 20 families each, using 26 generating models: (1) 4 intermediate models (penetrance of heterozygote between that of the two homozygotes); (2) 6 two-locus additive models; and (3) 16 two-locus heterogeneity models (admixture alpha = 1.0,.7,.5, and.3; alpha = 1.0 replicates simple Mendelian models). For LOD scores, we assumed dominant and recessive inheritance with 50% penetrance. We took the higher of the two maximum LOD scores and subtracted 0.3 to correct for multiple tests (MMLS-C). We compared expected maximum LOD scores and power, using MMLS-C and NPL as well as the true model. Since NPL uses only the affected family members, we also performed an affecteds-only analysis using MMLS-C. The MMLS-C was both uniformly more powerful than NPL for most cases we examined, except when linkage information was low, and close to the results for the true model under locus heterogeneity. We still found better power for the MMLS-C compared with NPL in affecteds-only analysis. The results show that use of two simple modes of inheritance at a fixed penetrance can have more power than NPL when the trait mode of inheritance is complex and when there is heterogeneity in the data set.     

Should maize doubled haploids be induced among F<Subscript>1</Subscript> or F<Subscript>2</Subscript> plants?

Maize (Zea mays L.) doubled haploid lines are typically produced from F₁ plants. Studies have suggested that the low frequency of recombinants in doubled haploids may reduce the response to selection. My objective was to determine if, for sustaining long-term response, doubled haploids should be induced in F₁ or F₂ plants during maize inbred development. In simulation experiments, I examined the response to multiple cycles of testcross selection among doubled haploid lines derived from F₁ plants (denoted by DH), doubled haploid lines derived from F₂ plants (DH_F2), and recombinant inbred (RI) lines derived by single-seed descent. For a trait controlled by 100 or more quantitative trait loci (QTL), the cumulative responses to selection were up to 4–6% larger among DH_F2 lines than among DH lines. The cumulative responses were up to 5–8% larger among RI lines than among DH lines. The QTL become unlinked as the number of QTL in a finite genome decreases, and the responses among RI, DH, and DH_F2 lines were equal or nearly equal when only 20 QTL controlled the trait. Metabolic-flux epistasis reduced the differences in the response among RI, DH, and DH_F2 lines. Overall, the results indicated that doubled haploids should be induced from F₂ plants rather than from F₁ plants. If year-round nurseries are used and new F₁ crosses for inbred development are initially created on a speculative basis, the development of doubled haploids from F₂ rather than F₁ plants should not cause a delay in inbred development.     

Region-Based Association Analysis of Human Quantitative Traits in Related Individuals

Regional-based association analysis instead of individual testing of each SNP was introduced in genome-wide association studies to increase the power of gene mapping, especially for rare genetic variants. For regional association tests, the kernel machine-based regression approach was recently proposed as a more powerful alternative to collapsing-based methods. However, the vast majority of existing algorithms and software for the kernel machine-based regression are applicable only to unrelated samples. In this paper, we present a new method for the kernel machine-based regression association analysis of quantitative traits in samples of related individuals. The method is based on the GRAMMAR+ transformation of phenotypes of related individuals, followed by use of existing kernel machine-based regression software for unrelated samples. We compared the performance of kernel-based association analysis on the material of the Genetic Analysis Workshop 17 family sample and real human data by using our transformation, the original untransformed trait, and environmental residuals. We demonstrated that only the GRAMMAR+ transformation produced type I errors close to the nominal value and that this method had the highest empirical power. The new method can be applied to analysis of related samples by using existing software for kernel-based association analysis developed for unrelated samples.     

Using molecular markers to map multiple quantitative trait loci: models for backcross,recombinant inbred,and doubled haploid progeny

Summary To maximize parameter estimation efficiency and statistical power and to estimate epistasis, the parameters of multiple quantitative trait loci (QTLs) must be simultaneously estimated. If multiple QTL affect a trait, then estimates of means of QTL genotypes from individual locus models are statistically biased. In this paper, I describe methods for estimating means of QTL genotypes and recombination frequencies between marker and quantitative trait loci using multilocus backcross, doubled haploid, recombinant inbred, and testcross progeny models. Expected values of marker genotype means were defined using no double or multiple crossover frequencies and flanking markers for linked and unlinked quantitative trait loci. The expected values for a particular model comprise a system of nonlinear equations that can be solved using an interative algorithm, e.g., the Gauss-Newton algorithm. The solutions are maximum likelihood estimates when the errors are normally distributed. A linear model for estimating the parameters of unlinked quantitative trait loci was found by transforming the nonlinear model. Recombination frequency estimators were defined using this linear model. Certain means of linked QTLs are less efficiently estimated than means of unlinked QTLs.     

QTL analysis: unreliability and bias in estimation procedures

Several statistical methods which employ multiple marker data are currently available for the analysis of quantitative trait loci (QTL) in experimental populations. Although comparable estimates of QTL location and effects have been obtained by these methods, using simulated and real data sets, their accuracy and reliability have not been extensively investigated. The present study specifically examines the merit of using F₂ and doubled haploid populations for locating QTL and estimating their effects. Factors which may affect accuracy and reliability of QTL mapping, such as the number and position of the markers available, the accuracy of the marker locations and the size of the experimental population used, are considered. These aspects are evaluated for QTL of differing heritabilities and locations along the chromosome.A population of 300 F₂ individuals and 150 doubled haploid lines gave estimates of QTL position and effect which were comparable, albeit extremely unreliable. Even for a QTL of high heritability (10%), the confidence interval was 35 cM. There was little increase in reliability to be obtained from using 300, rather than 200, F₂ individuals and 100 doubled haploid lines gave similar results to 150. QTL estimates were not significantly improved either by using the expected, rather than the observed, marker positions or by using a dense map of markers rather than a sparse map. A QTL which was asymmetrically located in the linkage group resulted in inaccurate estimates of QTL position which were seriously biassed at low heritability of the QTL. In a population of 300 F₂ individuals the bias increased from 4 cM to 20 cM, for a QTL with 10% and 2% heritability respectively.     

On the assessment of statistical significance in disease-gene discovery.

One of the major challenges facing genome-scan studies to discover disease genes is the assessment of the genomewide significance. The assessment becomes particularly challenging if the scan involves a large number of markers collected from a relatively small number of meioses. Typically, this assessment has two objectives: to assess genomewide significance under the null hypothesis of no linkage and to evaluate true-positive and false-positive prediction error rates under alternative hypotheses. The distinction between these goals allows one to formulate the problem in the well-established paradigm of statistical hypothesis testing. Within this paradigm, we evaluate the traditional criterion of LOD score 3.0 and a recent suggestion of LOD score 3.6, using the Monte Carlo simulation method. The Monte Carlo experiments show that the type I error varies with the chromosome length, with the number of markers, and also with sample sizes. For a typical setup with 50 informative meioses on 50 markers uniformly distributed on a chromosome of average length (i.e., 150 cM), the use of LOD score 3.0 entails an estimated chromosomewide type I error rate of.00574, leading to a genomewide significance level >.05. In contrast, the corresponding type I error for LOD score 3.6 is.00191, giving a genomewide significance level of slightly <.05. However, with a larger sample size and a shorter chromosome, a LOD score between 3.0 and 3.6 may be preferred, on the basis of proximity to the targeted type I error. In terms of reliability, these two LOD-score criteria appear not to have appreciable differences. These simulation experiments also identified factors that influence power and reliability, shedding light on the design of genome-scan studies.     

Verification of marker–trait associations in biparental winter barley (<Emphasis Type="Italic">Hordeum vulgare</Emphasis> L.) DH populations

In previous genome-wide association studies, marker–trait associations for grain yield and additional traits of agronomic importance were identified in the German winter barley (Hordeum vulgare L.) breeding gene pool. In the present study, seven doubled haploid populations segregating for the relevant alleles at the associated loci were used to get information whether these marker–trait associations can be verified in biparental populations and reliably used in applied barley breeding. The doubled haploid populations were phenotyped in field trials at two to five locations each in 1 year and genotyped by 40 trait-associated single nucleotide polymorphisms using an Illumina VeraCode GoldenGate assay. Large phenotypic variation was observed for all traits within at least one doubled haploid population. For 19 out of 58 marker–trait associations tested, the phenotypic means of both marker classes were significantly (p ≤ 0.005) different, thus confirming the association of the respective marker and the quantitative trait locus detected. For example, doubled haploid lines derived from a cross of ‘Malta’ × ‘Goldmine’ carrying different marker alleles differed by 0.41 t/ha in mean grain yield. The 19 (out of 58) marker–trait associations verified correspond to 10 (out of 27) genomic regions. Markers that were verified to be associated with a quantitative trait locus can be implemented directly in winter barley breeding for the selection of parental lines and marker-assisted pedigree selection.     

Monte Carlo simulations on marker grouping and ordering

Four global algorithms, maximum likelihood (ML), sum of adjacent LOD score (SALOD), sum of adjacent recombinant fractions (SARF) and product of adjacent recombinant fraction (PARF), and one approximation algorithm, seriation (SER), were used to compare the marker ordering efficiencies for correctly given linkage groups based on doubled haploid (DH) populations. The Monte Carlo simulation results indicated the marker ordering powers for the five methods were almost identical. High correlation coefficients were greater than 0.99 between grouping power and ordering power, indicating that all these methods for marker ordering were reliable. Therefore, the main problem for linkage analysis was how to improve the grouping power. Since the SER approach provided the advantage of speed without losing ordering power, this approach was used for detailed simulations. For more generality, multiple linkage groups were employed, and population size, linkage cutoff criterion, marker spacing pattern (even or uneven), and marker spacing distance (close or loose) were considered for obtaining acceptable grouping powers. Simulation results indicated that the grouping power was related to population size, marker spacing distance, and cutoff criterion. Generally, a large population size provided higher grouping power than small population size, and closely linked markers provided higher grouping power than loosely linked markers. The cutoff criterion range for achieving acceptable grouping power and ordering power differed for varying cases; however, combining all situations in this study, a cutoff criterion ranging from 50 cM to 60 cM was recommended for achieving acceptable grouping power and ordering power for different cases.     

Recent advances in human quantitative-trait-locus mapping: comparison of methods for selected sibling pairs

During the past few years, there has been a great deal of new work on methods for mapping quantitative-trait loci by use of sibling pairs and sibships. There are several new methods based on linear regression, as well as several more that are based on score statistics. In theory, most of the new methods should be relatively robust to violations of distributional assumptions and to selected sampling, but, in practice, there has been little evaluation of how the methods perform on selected samples. We survey most of the new regression-based statistics and score statistics and propose a few minor variations on the score statistics. We use simulation to evaluate the type I error and the power of all of the statistics, considering (a) population samples of sibling pairs and (b) sibling pairs ascertained on the basis of at least one sibling with a trait value in the top 10% of the distribution. Most of the statistics have correct type I error for selected samples. The statistics proposed by Xu et al. and by Sham and Purcell are generally the most powerful, along with one of our score statistic variants. Even among the methods that are most powerful for "nice" data, some are more robust than others to non-Gaussian trait models and/or misspecified trait parameters.     

An improved method for quantitative trait loci detection and identification of within-line segregation in F2 intercross designs

We present a new flexible, simple, and powerful genome-scan method (flexible intercross analysis, FIA) for detecting quantitative trait loci (QTL) in experimental line crosses. The method is based on a pure random-effects model that simultaneously models between- and within-line QTL variation for single as well as epistatic QTL. It utilizes the score statistic and thereby facilitates computationally efficient significance testing based on empirical significance thresholds obtained by means of permutations. The properties of the method are explored using simulations and analyses of experimental data. The simulations showed that the power of FIA was as good as, or better than, Haley-Knott regression and that FIA was rather insensitive to the level of allelic fixation in the founders, especially for pedigrees with few founders. A chromosome scan was conducted for a meat quality trait in an F(2) intercross in pigs where a mutation in the halothane (Ryanodine receptor, RYR1) gene with a large effect on meat quality was known to segregate in one founder line. FIA obtained significant support for the halothane-associated QTL and identified the base generation allele with the mutated allele. A genome scan was also performed in a previously analyzed chicken F(2) intercross. In the chicken intercross analysis, four previously detected QTL were confirmed at a 5% genomewide significance level, and FIA gave strong evidence (P < 0.01) for two of these QTL to be segregating within the founder lines. FIA was also extended to account for epistasis and using simulations we show that the method provides good estimates of epistatic QTL variance even for segregating QTL. Extensions of FIA and its applications on other intercross populations including backcrosses, advanced intercross lines, and heterogeneous stocks are also discussed.     

Trend Test for Overdispersed Proportions

The Cochran-Armitage test has commonly been used for a trend test in binomial proportions. The quasi-likelihood method provides a simple approach to model extra-binomial proportions. Two versions of the score and Wald tests using different parameterizations for the extra-binomial variance were investigated: one in terms of intercluster correlation, and another in terms of variance. The Monte Carlo simulation was used to evaluate the performance of the each version of the score test and the Wald test, and the Cochran-Armitage test. The simulation shows that the Cochran-Armitage test has the proper size only for the binomial sample data, and the test is no longer valid when applied to the extra-binomial data. The Wald test is more likely to exceed the nominal level than the score test under either intercluster correlation model or variance model. Both score tests performed very well even with the binomial data; the tests control the type I error and in the meantime maintain the power of detecting the dose effects. Based on the design considered in this paper, the two scores test are comparable. The score test based on the intercluster correlations model seems better controlling the Type I error but appears less powerful than that based on the variance model. An example from a developmental toxicity experiment is given.     

A nonparametric bootstrap method for testing close linkage vs. pleiotropy of coincident quantitative trait loci.

A novel method using the nonparametric bootstrap is proposed for testing whether a quantitative trait locus (QTL) at one chromosomal position could explain effects on two separate traits. If the single-QTL hypothesis is accepted, pleiotropy could explain the effect on two traits. If it is rejected, then the effects on two traits are due to linked QTLs. The method can be used in conjunction with several QTL mapping methods as long as they provide a straightforward estimate of the number of QTLs detectable from the data set. A selection step was introduced in the bootstrap procedure to reduce the conservativeness of the test of close linkage vs. pleiotropy, so that the erroneous rejection of the null hypothesis of pleiotropy only happens at a frequency equal to the nominal type I error risk specified by the user. The approach was assessed using computer simulations and proved to be relatively unbiased and robust over the range of genetic situations tested. An example of its application on a real data set from a saline stress experiment performed on a recombinant population of wheat (Triticum aestivum L. ) doubled haploid lines is also provided.     

MOD-score analysis with simple pedigrees: an overview of likelihood-based linkage methods

A MOD-score analysis, in which the parametric LOD score is maximized with respect to the trait-model parameters, can be a powerful method for the mapping of complex traits. With affected sib pairs, it has been shown before that MOD scores asymptotically follow a mixture of chi(2) distributions with 2, 1 and 0 degrees of freedom under the null hypothesis of no linkage. In that context, a MOD-score analysis yields some (albeit limited) information regarding the trait-model parameters, and there is a chance for an increased power compared to a simple LOD-score analysis. Here, it is shown that with unilineal affected relative pairs, MOD scores asymptotically follow a mixture of chi(2) distributions with 1 and 0 degrees of freedom under the null hypothesis, that is, the same distribution as followed by simple LOD scores. No information regarding the trait model can be obtained in this setting, and no power is gained when compared to a LOD-score analysis. An outlook to larger pedigrees is given. The number of degrees of freedom underlying the null distribution of MOD scores, that depends on the type of pedigrees studied, corresponds to the number of explored dimensions related to power and to the number of parameters that can jointly be estimated.     

Identification of genes for complex disease using longitudinal phenotypes

Using the simulated data set from Genetic Analysis Workshop 13, we explored the advantages of using longitudinal data in genetic analyses. The weighted average of the longitudinal data for each of seven quantitative phenotypes were computed and analyzed. Genome screen results were then compared for these longitudinal phenotypes and the results obtained using two cross-sectional designs: data collected near a single age (45 years) and data collected at a single time point. Significant linkage was obtained for nine regions (LOD scores ranging from 5.5 to 34.6) for six of the phenotypes. Using cross-sectional data, LOD scores were slightly lower for the same chromosomal regions, with two regions becoming nonsignificant and one additional region being identified. The magnitude of the LOD score was highly correlated with the heritability of each phenotype as well as the proportion of phenotypic variance due to that locus. There were no false-positive linkage results using the longitudinal data and three false-positive findings using the cross-sectional data. The three false positive results appear to be due to the kurtosis in the trait distribution, even after removing extreme outliers. Our analyses demonstrated that the use of simple longitudinal phenotypes was a powerful means to detect genes of major to moderate effect on trait variability. In only one instance was the power and heritability of the trait increased by using data from one examination. Power to detect linkage can be improved by identifying the most heritable phenotype, ensuring normality of the trait distribution and maximizing the information utilized through novel longitudinal designs for genetic analysis.     

An Evaluation of Beta-Binomial and Dirichlet-Trinomial Models for the Analysis of Reproductive and Developmental Effects

The common endpoints for the evaluation of reproductive and developmental toxic effects are the number of dead/resorbed fetuses, the number of malformed fetuses, and the number of normal fetuses for each litter. The joint distribution of the three endpoints could be modelled by a Dirichlettrinomial distribution or by a product of two-beta-binomial distributions. A simulation experiment is used to investigate the biases of the maximum likelihood estimate (MLE) for the probability of adverse effects under the Dirichlet-trinomial model and the beta-binomial model. Also, the type I errors and powers of the likelihood ratio test for comparing the difference between treatment and control are evaluated for the two underlying models. In estimation, the two MLE's are comparable, the bias estimates are small. In testing, the likelihood ratio test is generally more powerful under the Dirichlet-trinomial model than the beta-binomial model. The type I error rate is greater than the nominal level using the Dirichlet-trinomial model in some cases, when the data are generated from the two-beta-binomial model, and it is less than the nominal level using the beta-binomial model in other cases, when the data are generated from the Dirichlet-trinomial model.     

Influence of anther pretreatment and culture medium composition on the production of barley doubled haploids from model and low responding cultivars

Different pretreatments were given to anthers of barley before culturing, and their effects assessed on the frequency of embryos and green doubled haploid plants produced. Mannitol pretreatment was better than cold pretreatment for some low responding cultivars. Optimal concentration of mannitol for pretreatment depended on cultivar. Low responding genotypes needed a higher concentration of mannitol than responsive ones. The addition of Ficoll to liquid medium increased the number of embryos and green plants. The influence of the growth regulators 2,4-D and TIBA was assayed using ten cultivars of barley grown in Spain. The anti-auxin TIBA gave good embryo production with some of the low responding cultivars. Two row-type cultivars always produced higher number of embryos and green plantlets than six row-type. The application of these modifications to 10 F₁ hybrids with potential agronomic value, allowed the production of almost 1000 doubled haploid plants from only 3500 anthers. Up to two doubled haploid plants per flower were produced from the cross Monlon × Sonja. This revised version was published online in June 2006 with corrections to the Cover Date.     

The application of the entropy-based statistic for genomic association study of QTL

An entropy-based statistic TPE has been proposed for genomic association study for disease-susceptibility locus.The statistic TPE may be directly adopted and/or extended to quantitative-trait locus (QTL)mapping for quantitative traits.In this article,the statistic TPE was extended and applied to quantitative trait for association analysis of QTL by means of selective genotyping.The statistical properties (the type I error rate and the power) were examined under a range of parameters and population-sampling strategies (e.g.,various genetic models,various heritabilities,and various sample-selection threshold values) by simulation studies.The results indicated that the statistic Tee is robust and powerful for genomic association study of QTL.A simulation study based on the haplotype frequencies of 10 single nucleotide polymorphisms (SNPs) of angiotensin-I converting enzyme genes was conducted to evaluate the performance of the statistic TPE for genetic association study.