Mapping quantitative trait loci underlying triploid endosperm traits

Endosperm, which is derived from two polar nuclei fusing with one sperm, is a triploid tissue in cereals. Endosperm tissue determines the grain quality of cereals. Improving grain quality is one of the important breeding objectives in cereals. However, current statistical methods for mapping quantitative trait loci (QTL) under diploid genetic control have not been effective for dealing with endosperm traits because of the complexity of their triploid inheritance. In this paper, we derive for the first time the conditional probabilities of F(3) endosperm QTL genotypes given different flanking marker genotypes in F(2) plants. Using these probabilities, we develop a multiple linear regression method implemented via the iteratively reweighted least-squares (IRWLS) algorithm and a maximum likelihood method (ML) implemented via the expectation-maximization (EM) algorithm to map QTL underlying endosperm traits. We use the mean value of endosperm traits of F(3) seeds as the dependent variable and the expectations of genotypic indicators for additive and dominance effect of a putative QTL flanked by a pair of markers as independent variables for IRWLS mapping. However, if an endosperm trait is measured quantitatively using a single endosperm sample, the ML mapping method can be used to separate the two dominance effects. Efficiency of the methods is verified through extensive Monte Carlo simulation studies. Results of simulation show that the proposed methods provide accurate estimates of both the QTL effects and locations with very high statistical power. With these methods, we are now ready to map endosperm traits, as we can for regular quantitative trait under diploid control.     

Backward Haplotype Transmission Association (BHTA) algorithm - a fast multiple-marker screening method

The mapping of complex traits is one of the most important and central areas of human genetics today. Recent attention has been focused on genome scans using a large number of marker loci. Because complex traits are typically caused by multiple genes, the common approaches of mapping them by testing markers one after another fail to capture the substantial information of interactions among disease loci. Here we propose a backward haplotype transmission association (BHTA) algorithm to address this problem. The algorithm can administer a screening on any disease model when case-parent trio data are available. It identifies the important subset of an original larger marker set by eliminating the markers of least significance, one at a time, after a complete evaluation of its importance. In contrast with the existing methods, three major advantages emerge from this approach. First, it can be applied flexibly to arbitrary markers, regardless of their locations. Second, it takes into account haplotype information; it is more powerful in detecting the multifactorial traits in the presence of haplotypic association. Finally, the proposed method can potentially prove to be more efficient in future genomewide scans, in terms of greater accuracy of gene detection and substantially reduced number of tests required in scans. We illustrate the performance of the algorithm with several examples, including one real data set with 31 markers for a study on the Gilles de la Tourette syndrome. Detailed theoretical justifications are also included, which explains why the algorithm is likely to select the 'correct' markers.     

Population structure and marker–trait associations for pomological traits in peach and nectarine cultivars

Marker–trait associations based on populations from controlled crosses have been established in peach using markers mapped on the peach consensus map. In this study, we explored the utility of unstructured populations for association mapping to determine useful marker–trait associations in peach/nectarine cultivars. We used 94 peach cultivars representing local Spanish and modern cultivars from international breeding programs that are maintained at the Experimental Station of Aula Dei, Spain. This collection was characterized for pomological traits and was screened with 40 SSR markers that span the peach genome. Population structure analysis using STRUCTURE software identified two subpopulations, the local and modern cultivars, with admixture within both groups. The local Spanish cultivars were somewhat less diverse than modern cultivars. Marker–trait associations were determined in TASSEL with and without modelling coefficient of membership (Q) values as covariates. The results showed significant associations with pomological traits. We chose three markers on LG4 because of their proximity to the endoPG locus (freestone–melting flesh) that strongly affects pomological traits. Two genotypes of BPPCT015 marker showed significant associations with harvest date, flavonoids and sorbitol. Also, two genotypes of CPPCT028 showed associations with harvest date, total phenolics, RAC, and total sugars. Finally, two genotypes of endoPG1 showed associations with flesh firmness and total sugars. The analysis of linkage disequilibrium (LD) revealed a high level of LD up to 20 cM, and decay at farther distances. Therefore, association mapping could be a powerful tool for identifying marker–trait associations and would be useful for marker-assisted selection in peach breeding.     

Detection of genes for ordinal traits in nuclear families and a unified approach for association studies

There is growing interest in genomewide association analysis using single-nucleotide polymorphisms (SNPs), because traditional linkage studies are not as powerful in identifying genes for common, complex diseases. Tests for linkage disequilibrium have been developed for binary and quantitative traits. However, since many human conditions and diseases are measured in an ordinal scale, methods need to be developed to investigate the association of genes and ordinal traits. Thus, in the current report we propose and derive a score test statistic that identifies genes that are associated with ordinal traits when gametic disequilibrium between a marker and trait loci exists. Through simulation, the performance of this new test is examined for both ordinal traits and quantitative traits. The proposed statistic not only accommodates and is more powerful for ordinal traits, but also has similar power to that of existing tests when the trait is quantitative. Therefore, our proposed statistic has the potential to serve as a unified approach to identifying genes that are associated with any trait, regardless of how the trait is measured. We further demonstrated the advantage of our test by revealing a significant association (P = 0.00067) between alcohol dependence and a SNP in the growth-associated protein 43.     

Genetic diversity among and within Iranian and non-Iranian barely (Hordeum vulgare L.) genotypes using SSR and storage proteins markers

The objectives of this study were evaluation of genetic diversity and marker–trait association of 64 barley (Hordeum vulgare L.) genotypes using hordeins and simple sequence repeats (SSRs) markers under optimal moisture and drought stress conditions. Moreover, to evaluate the response of barley genotypes to drought stress, five drought tolerance indices were calculated. SSRs and hordeins generated clear patterns with high polymorphism. SSRs and hordeins analysis provided us with useful information on the level of polymorphism and diversity in barley. Marker–trait associations were studied for 22 agronomic traits using 122 SSR markers (obtained from 14 primer pairs) and 51 hordeins bands in 64 barley genotypes under both normal and stress conditions. Phenotypic traits strongly associated with SSRs were also strongly associated with hordeins. Generally, we believed that at least some of these markers would be informative and validated and can be used in marker-assisted selection (MAS) under drought stress.     

Genetic design and statistical power of nested association mapping in maize

We investigated the genetic and statistical properties of the nested association mapping (NAM) design currently being implemented in maize (26 diverse founders and 5000 distinct immortal genotypes) to dissect the genetic basis of complex quantitative traits. The NAM design simultaneously exploits the advantages of both linkage analysis and association mapping. We demonstrated the power of NAM for high-power cost-effective genome scans through computer simulations based on empirical marker data and simulated traits with different complexities. With common-parent-specific (CPS) markers genotyped for the founders and the progenies, the inheritance of chromosome segments nested within two adjacent CPS markers was inferred through linkage. Genotyping the founders with additional high-density markers enabled the projection of genetic information, capturing linkage disequilibrium information, from founders to progenies. With 5000 genotypes, 30-79% of the simulated quantitative trait loci (QTL) were precisely identified. By integrating genetic design, natural diversity, and genomics technologies, this new complex trait dissection strategy should greatly facilitate endeavors to link molecular variation with phenotypic variation for various complex traits.     

An Effective Method to Identify Heritable Components from Multivariate Phenotypes

Multivariate phenotypes may be characterized collectively by a variety of low level traits, such as in the diagnosis of a disease that relies on multiple disease indicators. Such multivariate phenotypes are often used in genetic association studies. If highly heritable components of a multivariate phenotype can be identified, it can maximize the likelihood of finding genetic associations. Existing methods for phenotype refinement perform unsupervised cluster analysis on low-level traits and hence do not assess heritability. Existing heritable component analytics either cannot utilize general pedigrees or have to estimate the entire covariance matrix of low-level traits from limited samples, which leads to inaccurate estimates and is often computationally prohibitive. It is also difficult for these methods to exclude fixed effects from other covariates such as age, sex and race, in order to identify truly heritable components. We propose to search for a combination of low-level traits and directly maximize the heritability of this combined trait. A quadratic optimization problem is thus derived where the objective function is formulated by decomposing the traditional maximum likelihood method for estimating the heritability of a quantitative trait. The proposed approach can generate linearly-combined traits of high heritability that has been corrected for the fixed effects of covariates. The effectiveness of the proposed approach is demonstrated in simulations and by a case study of cocaine dependence. Our approach was computationally efficient and derived traits of higher heritability than those by other methods. Additional association analysis with the derived cocaine-use trait identified genetic markers that were replicated in an independent sample, further confirming the utility and advantage of the proposed approach.     

Genetic distance predicts trait differentiation at the subpopulation but not the individual level in eelgrass,Zostera marina

Ecological studies often assume that genetically similar individuals will be more similar in phenotypic traits, such that genetic diversity can serve as a proxy for trait diversity. Here, we explicitly test the relationship between genetic relatedness and trait distance using 40 eelgrass (Zostera marina) genotypes from five sites within Bodega Harbor, CA. We measured traits related to nutrient uptake, morphology, biomass and growth, photosynthesis, and chemical deterrents for all genotypes. We used these trait measurements to calculate a multivariate pairwise trait distance for all possible genotype combinations. We then estimated pairwise relatedness from 11 microsatellite markers. We found significant trait variation among genotypes for nearly every measured trait; however, there was no evidence of a significant correlation between pairwise genetic relatedness and multivariate trait distance among individuals. However, at the subpopulation level (sites within a harbor), genetic (F_ST) and trait differentiation were positively correlated. Our work suggests that pairwise relatedness estimated from neutral marker loci is a poor proxy for trait differentiation between individual genotypes. It remains to be seen whether genomewide measures of genetic differentiation or easily measured “master” traits (like body size) might provide good predictions of overall trait differentiation.     

Dissecting the correlation structure of a bivariate phenotype: Common genes or shared environment?

High correlations between two quantitative traits may be either due to common genetic factors or common environmental factors or a combination of both. In this study, we develop statistical methods to extract the genetic contribution to the total correlation between the components of a bivariate phenotype. Using data on bivariate phenotypes and marker genotypes for sib-pairs, we propose a test for linkage between a common QTL and a marker locus based on the conditional cross-sib trait correlations (trait 1 of sib 1—trait 2 of sib 2 and conversely) given the identity-by-descent (i.b.d.) sharing at the marker locus. We use Monte-Carlo simulations to evaluate the performance of the proposed test under different trait parameters and quantitative trait distributions. An application of the method is illustrated using data on two alcohol-related phenotypes from a project on the collaborative study on the genetics of alcoholism.     

Complex inheritance and localizing disease genes.

Most methods for localizing genes underlying complex traits work under the implicit or explicit assumption of a single disease gene with the possible exception of heterogeneity, that is, different disease genes in different families. We discuss current single-locus and multi-locus methods. Novel approaches are proposed that take into account all marker loci over the genome. A simple example is given for an unconventional statistic, i.e. the mean of allele sharing over all markers on a chromosome.     

A test of transmission/disequilibrium for quantitative traits in pedigree data, by multiple regression.

The transmission/disequilibrium (TD) test (TDT), proposed, by Spielman et al., for binary traits is a powerful method for detection of linkage between a marker locus and a disease locus, in the presence of allelic association. As a test for linkage disequilibrium, the TDT makes the assumption that any allelic association present is due to linkage. Allison proposed a series of TD-type tests for quantitative traits and calculated their power, assuming that the marker locus is the disease locus. All these tests assume that the observations are independent, and therefore they are applicable, as a test for linkage, only for nuclear-family data. In this report, we propose a regression-based TD-type test for linkage between a marker locus and a quantitative trait locus, using information on the parent-to-offspring transmission status of the associated allele at the marker locus. This method does not require independence of observations, thus allowing for analysis of pedigree data as well, and allows adjustment for covariates. We investigate the statistical power and validity of the test by simulating markers at various recombination fractions from the disease locus.     

A recoding scheme for X-linked and pseudoautosomal loci to be used with computer programs for autosomal LOD-score analysis

We present a recoding scheme that allows for a parametric multipoint X-chromosomal linkage analysis of dichotomous traits in the context of a computer program for autosomes that can use trait models with imprinting. Furthermore, with this scheme, it is possible to perform a joint multipoint analysis of X-linked and pseudoautosomal loci. It is required that (1) the marker genotypes of all female nonfounders are available and that (2) there are no male nonfounders who have daughters in the pedigree. The second requirement does not apply if the trait locus is pseudoautosomal. The X-linked marker loci are recorded by adding a dummy allele to the males' hemizygous genotypes. For modelling an X-linked trait locus, five different liability classes are defined, in conjunction with a paternal imprinting model for male nonfounders. The formulation aims at the mapping of a diallelic trait locus relative to an arbitrary number of codominant markers with known genetic distances, in cases where a program for a genuine X-chromosomal analysis is not available.     

Marker assisted selection in crop plants

Genetic mapping of major genes and quantitative traits loci (QTLs) for many important agricultural traits is increasing the integration of biotechnology with the conventional breeding process. Exploitation of the information derived from the map position of traits with agronomical importance and of the linked molecular markers, can be achieved through marker assisted selection (MAS) of the traits during the breeding process. However, empirical applications of this procedure have shown that the success of MAS depends upon several factors, including the genetic base of the trait, the degree of the association between the molecular marker and the target gene, the number of individuals that can be analyzed and the genetic background in which the target gene has to be transferred. MAS for simply inherited traits is gaining increasing importance in breeding programs, allowing an acceleration of the breeding process. Traits related to disease resistance to pathogens and to the quality of some crop products are offering some important examples of a possible routinary application of MAS. For more complex traits, like yield and abiotic stress tolerance, a number of constraints have determined severe limitations on an efficient utilization of MAS in plant breeding, even if there are a few successful applications in improving quantitative traits. Recent advances in genotyping technologies together with comparative and functional genomic approaches are providing useful tools for the selection of genotypes with superior agronomical performancies.     

Regression association analysis of fruit traits with molecular markers in cherries

The use of molecular markers supports the study of genetic marker–trait association of biological and agronomic interest in diverse genetic material. In this research, association between simple sequence repeat (SSR) and random amplified polymorphic DNA (RAPD) markers with fruit traits were investigated in two collections of cherries by applying multiple regression analysis (MRA). Thirty-eight SSR alleles and 135 RAPD fragments were found associated with 14 of affecting fruit traits. Some of SSR and RAPD markers were associated with more than one fruit trait in MRA. Such an association may arise due to pleiotropic effect of the linked quantitative trait locus on different traits. For example, some SSR and RAPD markers were associated with all four traits including fruit cracking, fruit firmness, total soluble solid (TSS) and fruit shape. Also, some markers had correlations with all four characters of TSS, anthocyanin, fruit skin color and fruit flesh color, indicating a significant correlation among these traits. Therefore, it is possible to use these markers along with morphological traits in cherry breeding programs for identification of suitable parents to produce mapping populations and hybrid cultivars. Also, these results could be useful in marker-assisted breeding programs when no other genetic information is available.     

Some practical considerations for using RFLP markers to aid in selection during inbreeding of maize

Summary If molecular markers are to be routinely used in maize (Zea mays L.) breeding for selection of quantitative trait loci (QTL), then consistent marker-trait associations across breeding populations are needed, as are efficient methods for weighting information from different markers. Given 15 restriction fragment length polymorphism (RFLP) markers associated with grain yield in testcrosses of 220 [BS11(FR)C7 x FRMol7] F₂ individuals to FRB73, separate weighting schemes were attempted in order to maximize the frequency of favorable marker genotypes associated with increased grain yield in selected F₂ individuals and F₂:S₄ Unes. The following principles were apparent: (1) Differential weighting among markers, in addition to weighting individual marker genotypes on the basis of associated mean effects, should be emphasized when using markers to select in breeding populations. This is due to limited population sizes that can readily be handled. (2) Relatively few markers may need to be used to screen segregating populations (e.g., F₂) of limited size for loci affecting complex traits, such as combining ability for grain yield, assuming prior knowledge of marker-QTL associations. Markers given greatest weight (largest estimates of associated effects) will determine most selections. (3) When marker-based selection is among individuals at higher levels of inbreeding (e.g., S₄) within selected families, more markers need to be used in screening because those associated with relatively small effects have an increased chance of affecting selection.These results suggest a qualitative approach for utilizing RFLP markers to aid in selection of complex traits in commercial hybrid maize breeding programs. Commercial research programs produce thousands of crosses each year aimed at inbred line development. Discovery of molecular markers with consistent QTL associations across breeding populations and close QTL linkages would allow for rapid screening of new F₂ populations at a few key markers. Early elimination of individuals with undesirable genotypes would reduce the extent of hybrid performance testing necessary during later stages of inbreeding.     

Mapping quantitative trait loci using binned genotypes

Precise mapping of quantitative trait loci(QTLs)is critical for assessing genetic effects and identifying candidate genes for quantitative traits.Interval and composite interval mappings have been the methods of choice for several decades,which have provided tools for identifying genomic regions harboring causal genes for quantitative traits.Historically,the concept was developed on the basis of sparse marker maps where genotypes of loci within intervals could not be observed.Currently,genomes of many organisms have been saturated with markers due to the new sequencing technologies.Genotyping by sequencing usually generates hundreds of thousands of single nucleotide polymorphisms(SNPs),which often include the causal polymorphisms.The concept of interval no longer exists,prompting the necessity of a norm change in QTL mapping technology to make use of the high-volume genomic data.Here we developed a statistical method and a software package to map QTLs by binning markers into haplotype blocks,called bins.The new method detects associations of bins with quantitative traits.It borrows the mixed model methodology with a polygenic control from genome-wide association studies(GWAS)and can handle all kinds of experimental populations under the linear mixed model(LMM)framework.We tested the method using both simulated data and data from populations of rice.The results showed that this method has higher power than the current methods.An R package named binQTL is available from GitHub.