High-resolution quantitative trait locus mapping for body weight in mice by recombinant progeny testing

A major obstacle to the positional cloning of quantitative trait loci (QTLs) lies in resolving genetic factors whose allelic effects are blurred by environmental and background genetic variation. We investigate a fine-mapping approach that combines the use of an interval-specific congenic strain with progeny testing of recombinants for markers flanking a QTL. We apply the approach to map a murine QTL with an approximately 20% effect on growth rate by progeny testing 39 recombinants in a 12 cM region of the X chromosome. We use a likelihood analysis in an attempt to maximize the information on QTL map location and effect. The major X-linked effect is mapped to an approximately 2 cM region flanked by markers about 5 cM apart, outside which LOD support for the QTL drops extremely steeply by about 80. Nearly unambiguous assignment of the QTL genotypic state is obtained for each recombinant. The resolution of individual recombinants in the region is therefore sufficiently high to facilitate the positional cloning of the locus, although progress has been hampered because the genomic region containing the QTL shows an exceptionally low level of polymorphism in comparison with recent studies.     

The genetics of mating recognition between Drosophila simulans and D. sechellia

During courtship, visual and chemical signals are often exchanged between the sexes. The proper exchange of such signals ensures intraspecific recognition. We have examined the genetic basis of interspecific differences in male mating behaviour and pheromone concentration between Drosophila simulans and D. sechellia by using Drosophila simulans/D. sechellia introgression lines. Our results show a majority of quantitative trait loci (QTLs) explaining variation in both male mating behaviour and pheromone concentration to be located on the third chromosome. One QTL found on the third chromosome explains variation in time needed to start courtship and copulation as well as time spent courting. The position of such QTL (approximately 84A-88B) with effects on courtship and copulation aspects of mating includes the candidate sex determination gene doublesex (84E5-6) and Voila (86E1-2), a gene that affects male courtship in D. melanogaster. One additional third chromosome QTL explained variation in 7-tricosene pheromone concentrations among males. The interval mapping position of this QTL (approximately 68E-76E) did not overlap with the position detected for differences in mating behaviour and the intervals did not include candidate genes previously identified as having an effect on D. melanogaster cuticular hydrocarbon production. We did not detect any directionality of the effect of Drosophila sechellia allele introgressions in male mating recognition.     

Dissection and fine mapping of a major QTL for preharvest sprouting resistance in white wheat Rio Blanco

Preharvest sprouting (PHS) is a major constraint to white wheat production. Previously, we mapped quantitative trait loci (QTL) for PHS resistance in white wheat by using a recombinant inbred line (RIL) population derived from the cross Rio Blanco/NW97S186. One QTL, QPhs.pseru-3A, showed a major effect on PHS resistance, and three simple sequence repeat (SSR) markers were mapped in the QTL region. To determine the flanking markers for the QTL and narrow down the QTL to a smaller chromosome region, we developed a new fine mapping population of 1,874 secondary segregating F₂ plants by selfing an F6 RIL (RIL25) that was heterozygous in the three SSR marker loci. Segregation of PHS resistance in the population fitted monogenic inheritance. An additive effect of the QTL played a major role on PHS resistance, but a dominant effect was also observed. Fifty-six recombinants among the three SSR markers were identified in the population and selfed to produce homozygous recombinants or QTL near-isogenic lines (NIL). PHS evaluation of the recombinants delineated the QTL in the region close to Xbarc57 flanked by Xbarc321 and Xbarc12. To saturate the QTL region, 11 amplified fragment length polymorphism (AFLP) markers were mapped in the QTL region with 7 AFLP co-segregated with Xbarc57 by using the NIL population. Dissection of the QTL as a Mendelian factor and saturation of the QTL region with additional markers created a solid foundation for positional cloning of the major QTL.     

The genetic basis of behavioral isolation between Drosophila mauritiana and D. sechellia

Understanding how species form is a fundamental question in evolutionary biology. Identifying the genetic bases of barriers that prevent gene flow between species provides insight into how speciation occurs. Here, I analyze a poorly understood reproductive isolating barrier, prezygotic reproductive isolation. I perform a genetic analysis of prezygotic isolation between two closely related species of Drosophila, D. mauritiana and D. sechellia. I first confirm the existence of strong behavioral isolation between D. mauritiana females and D. sechellia males. Next, I examine the genetic basis of behavioral isolation by (1) scanning an existing set of introgression lines for chromosomal regions that have a large effect on isolation; and (2) mapping quantitative trait loci (QTL) that underlie behavioral isolation via backcross analysis. In particular, I map QTL that determine whether a hybrid backcross female and a D. sechellia male will mate. I identify a single significant QTL, on the X chromosome, suggesting that few major-effect loci contribute to behavioral isolation between these species. In further work, I refine the map position of the QTL to a small region of the X chromosome.     

Efficient mapping of plant height quantitative trait loci in a sorghum association population with introgressed dwarfing genes

Of the four major dwarfing genes described in sorghum, only Dw3 has been cloned. We used association mapping to characterize the phenotypic effects of the dw3 mutation and to fine map a second, epistatic dwarfing QTL on sorghum chromosome 9 (Sb-HT9.1). Our panel of 378 sorghum inbreds includes 230 sorghum conversion (SC) lines, which are exotic lines that have been introgressed with dwarfing quantitative trait loci (QTL) from a common parent. The causal mutation in dw3 associates with reduced lower internode length and an elongation of the apex, consistent with its role as an auxin efflux carrier. Lines carrying the dw3 mutation display high haplotype homozygosity over several megabases in the Dw3 region, but most markers linked to Dw3 do not associate significantly with plant height due to allele sharing between Dw3 and dw3 individuals. Using markers with a high mutation rate and the dw3 mutation as an interaction term, significant trait associations were detected across a 7-Mb region around Sb-HT9.1, largely due to higher detection power in the SC lines. Conversely, the likely QTL interval for Sb-HT9.1 was reduced to approximately 100 kb, demonstrating that the unique structure of this association panel provides both power and resolution for a genomewide scan.     

Genetics of egg production in Drosophila sechellia

Drosophila sechellia, an island endemic that specializes on a single host plant, has a lower rate of egg production than its generalist sister species D. melanogaster, D. simulans, and D. mauritiana. Earlier work showed that part of this difference in egg production was due to a reduction in the number of ovarioles in D. sechellia relative to its sister species. Here, I extend this earlier work by genetically analyzing the difference in egg production between D. sechellia and D. simulans. In all, 10 genetic markers were used in several interspecific backcrosses to identify chromosome regions that affected the rate of egg production. While previously mapped factors affecting ovariole number appear to impact the rate of egg production, new, non-ovariole factors were also identified. Overall, the difference in egg production between D. sechellia and D. simulans appears to be a polygenic trait. The relationship between these factors and genes involved the adaptation of D. sechellia to its host plant is not yet clear. The data are consistent with the hypothesis that decline in egg production is, in part, a negative pleiotropic effect of genetic changes required for host specialization in D. sechellia, although finer-scale genetic analysis of both traits is needed to truly test this hypothesis.     

Quantitative trait loci for cuticular hydrocarbons associated with sexual isolation between Drosophila simulans and D. sechellia

The identification of genes with large effects on sexual isolation and speciation is an important link between classic evolutionary genetics and molecular biology. Few genes that affect sexual isolation and speciation have been identified, perhaps because many traits influencing sexual isolation are complex behaviors. Cuticular hydrocarbons (CHs) of species of the Drosophila melanogaster group play a large role in sexual isolation by functioning as contact pheromones influencing mate recognition. Some of the genes that play key roles in determining species-specific CHs have been identified. We have performed separate quantitative trait locus (QTL) analyses of 7-tricosene (7-T) and 7,11-heptacosadiene (7,11-HD), the two major female CHs differing between D. simulans and D. sechellia. We find that approximately 40% of the phenotypic variance in each CH is associated with two to four chromosomal regions. A region on the right arm of chromosome 3 contains QTL that affect both traits, but other QTL are in distinct chromosomal regions. Epistatic interactions were detected between two pairs of QTL for 7,11-HD such that if either were homozygous for the D. simulans allele, the fly was similar to D. simulans in phenotype, with a low level of 7,11-HD. We discuss the location of these regions with regard to candidate genes for CH production, including those for desaturases.     

Fine mapping of qhir1 influencing in vivo haploid induction in maize

Production of haploids by the in vivo haploid induction method has now become routine for generating new inbred lines in maize. In previous studies, a major quantitative trait locus (QTL) (qhir1) located in bin 1.04 was detected, explaining up to 66 % of the genotypic variance for haploid induction rate (HIR). Our objectives were to (1) fine-map qhir1 and (2) identify closely linked markers useful for marker-assisted breeding of new inducers. For this purpose, we screened a mapping population of 14,375 F₂ plants produced from a cross between haploid inducer UH400 and non-inducer line 1680 to identify recombinants. Based on sequence information from the B73 reference genome, markers polymorphic between the two parents were developed to conduct fine mapping with these recombinants. A progeny test mapping strategy was applied to accurately determine the HIR of the 14 recombinants identified. Furthermore, F₃ progeny of recombinant F₂ plants were genotyped and in parallel evaluated for HIR. We corroborated earlier studies in that qhir1 has both a significantly positive effect on HIR but also a strong selective disadvantage, as indicated by significant segregation distortion. Altogether, we were able to narrow down the qhir1 locus to a 243 kb region flanked by markers X291 and X263.     

High resolution mapping and identification of new quantitative trait loci (QTL) affecting susceptibility to Marek's disease

Marek's disease (MD) is a lymphoproliferative disease of chickens that costs the poultry industry approximately $1 billion annually. Genetic resistance to MD is gaining increased attention to augment vaccinal control as disease outbreaks occur more frequently. Previously, analysis of a 272 F2 White Leghorn resource population measured for many MD traits and genotyped for 78 microsatellite markers revealed two and four quantitative trait loci (QTL) with significant and suggestive association, respectively, to one or more MD associated traits. Additional genetic markers have since been scored on the MD resource population to increase QTL resolution and genome coverage. Saturation of four of the QTL regions with 17 markers revealed five new QTL while 32 markers extended the genome coverage by 400 + CM and uncovered three more QTL. QTL analysis by single-point and interval mapping algorithms agreed well when marker saturation was approximately 20 CM or less. Currently 127 genetic markers cover approximately 68% of the genome that contain up to 14 MD QTL associated to one or more MD trait; seven at the significant level and seven at the suggestive level. Individually each QTL accounts for 2-10% of the variation and, in general, resistance was dominant although the resistant allele may come from either parental line. This study suggests that a limited number of genomic regions play a major role in the genetic control of MD resistance. Markers linked to these loci may be useful for selection of MD resistant stock by the poultry industry following verification of the association within their breeding populations.     

High-resolution joint linkage disequilibrium and linkage mapping of quantitative trait loci based on sibship data

This paper proposes variance component models for high resolution joint linkage disequilibrium (LD) and linkage mapping of quantitative trait loci (QTL) based on sibship data; this can include population data if independent individuals are treated as single sibships. One application of these models is late onset complex disease gene mapping, when parental data are not available. The models simultaneously incorporate both LD and linkage information. The LD information is contained in mean coefficients of sibship data. The linkage information is contained in the variance-covariance matrices of trait values for sibships with at least two siblings. We derive formulas for calculating the probability of sharing two trait alleles identical by descent (IBD) for sibpairs in interval mapping of QTL; this is the coefficient of dominant variance of the trait covariance of sibpairs on major QTL. To investigate the performance of the formulas, we calculate the numerical values via the formulas and get satisfactory approximations. We compare the power and sample sizes for both LD and linkage mapping. By simulation and theoretical analysis, we compare the results with those of Fulker and Abecasis "AbAw" approach. It is well known that the resolution of linkage analysis can be low for complex disease gene mapping. LD mapping, on the other hand, can increase mapping precision and is useful in high resolution mapping. Linkage analysis is less sensitive to population subdivisions and admixtures. The level of LD is sensitive to population stratification which may easily lead to spurious association. Performing a joint analysis of LD and linkage mapping can help to overcome the limits of both approaches. Moreover, the advantages of the two complementary strategies can be utilized maximally. In practice, linkage analysis may be performed using pedigree data to identify suggestive linkage between markers and trait loci based on a sparse marker map. In the presence of linkage, joint LD and linkage mapping can be carried out to do fine gene mapping based on a dense genetic map using both pedigree and population data. Population and pedigree data of any type can be combined to perform a joint analysis of high resolution LD and linkage mapping of QTL by generalizing the method.     

Association mapping of agronomic traits on chromosome 2A of wheat

Association mapping is a method to test the association between molecular markers and quantitative trait loci (QTL) based on linkage disequilibrium (LD). In this study, the collection of 108 wheat germplasm accessions form China were evaluated for their plant heights, spike length, spikelets per spike, grains per spike, thousand kernel weight and spikelets density in 3 years at three locations. And they were genotyped with 85 SSR markers and 40 EST-SSR markers. The population structure was inferred on the basis of unlinked 48 SSR markers and 40 EST-SSR markers. The extent of LD on chromosome 2A was 2.3 cM. Association of 37 SSR loci on chromosomes 2A with six agronomic traits was analysed with a mixed linear model. A total of 14 SSR loci were significantly associated with agronomic traits. Some of the associated markers were located in the QTL region detected in previous linkage mapping analysis. Our results demonstrated that association mapping can enhance QTL information and achieves higher resolution with short LD extent.     

Precise mapping Fhb5, a major QTL conditioning resistance to Fusarium infection in bread wheat (Triticum aestivum L.)

Qfhi.nau-5A is a major quantitative trait locus (QTL) against Fusarium graminearum infection in the resistant wheat germplasm Wangshuibai. Genetic analysis using BC(3)F(2) and BC(4)F(2) populations, derived from selfing two near-isogenic lines (NIL) heterozygous at Qfhi.nau-5A that were developed, respectively, with Mianyang 99-323 and PH691 as the recurrent parent, showed that Qfhi.nau-5A inherited like a single dominant gene. This QTL was thus designated as Fhb5. To fine map it, these two backcross populations and a recombinant inbred line (RIL) population derived from Nanda2419?×?Wangshuibai were screened for recombinants occurring between its two flanking markers Xbarc56 and Xbarc100. Nineteen NIL recombinants were identified from the two backcross populations and nine from the RIL population. In the RIL recombinant selection process, selection against Fhb4 present in the RIL population was incorporated. Genotyping these recombinant lines with ten markers mapping to the Xbarc56-Xbarc100 interval revealed four types of Mianyang 99-323-derived NIL recombinants, three types of PH691-derived NIL recombinants, and four types of RIL recombinants. In different field trials, the percentage of infected spikes of these lines displayed a distinct two-peak distribution. The more resistant class had over 55% less infection than the susceptible class. Common to these resistant genotypes, the 0.3-cM interval flanked by Xgwm304 and Xgwm415 or one of these two loci was derived from Wangshuibai, while none of the susceptible recombinants had Wangshuibai chromatin in this interval. This interval harboring Fhb5 was mapped to the pericentromeric C-5AS3-0.75 bin through deletion bin mapping. The precise localization of Fhb5 will facilitate its utilization in marker-assisted wheat breeding programs.     

Multipoint genetic mapping of quantitative trait loci with dominant markers in outbred populations

We present a multipoint algorithm for mapping quantitative trait loci (QTLs) using dominant markers. The algorithm is designed for outbred populations and is particularly suited for large families. The algorithm works with either codominant or dominant markers, either of which may be interspersed within the same linkage map. Concurrently, the algorithm also partitions dominance variance at the QTL. Computer simulations show that with large families, QTL mapping with dominant markers can be almost as powerful as with bi-allelic, codominant markers. Yet despite this, other situations show a large standard deviation in the estimate of the QTL position, thus making QTL mapping with dominant markers in outbred populations a useful detection tool, albeit limited in its resolution. This revised version was published online in July 2006 with corrections to the Cover Date.     

Simultaneous detection and fine mapping of quantitative trait loci in mice using heterogeneous stocks

We describe a method to simultaneously detect and fine map quantitative trait loci (QTL) that is especially suited to the mapping of modifier loci in mouse mutant models. The method exploits the high level of historical recombination present in a heterogeneous stock (HS), an outbred population of mice derived from known founder strains. The experimental design is an F(2) cross between the HS and a genetically distinct line, such as one carrying a knockout or transgene. QTL detection is performed by a standard genome scan with approximately 100 markers and fine mapping by typing the same animals using densely spaced markers over those candidate regions detected by the scan. The analysis uses an extension of the dynamic-programming technique employed previously to fine map QTL in HS mice. We show by simulation that a QTL accounting for 5% of the total variance can be detected and fine mapped with >50% probability to within 3 cM by genotyping approximately 1500 animals.     

Enhanced efficiency of quantitative trait loci mapping analysis based on multivariate complexes of quantitative traits

An approach to increase the efficiency of mapping quantitative trait loci (QTL) was proposed earlier by the authors on the basis of bivariate analysis of correlated traits. The power of QTL detection using the log-likelihood ratio (LOD scores) grows proportionally to the broad sense heritability. We found that this relationship holds also for correlated traits, so that an increased bivariate heritability implicates a higher LOD score, higher detection power, and better mapping resolution. However, the increased number of parameters to be estimated complicates the application of this approach when a large number of traits are considered simultaneously. Here we present a multivariate generalization of our previous two-trait QTL analysis. The proposed multivariate analogue of QTL contribution to the broad-sense heritability based on interval-specific calculation of eigenvalues and eigenvectors of the residual covariance matrix allows prediction of the expected QTL detection power and mapping resolution for any subset of the initial multivariate trait complex. Permutation technique allows chromosome-wise testing of significance for the whole trait complex and the significance of the contribution of individual traits owing to: (a) their correlation with other traits, (b) dependence on the chromosome in question, and (c) both a and b. An example of application of the proposed method on a real data set of 11 traits from an experiment performed on an F(2)/F(3) mapping population of tetraploid wheat (Triticum durum x T. dicoccoides) is provided.     

Overview of LASSO-related penalized regression methods for quantitative trait mapping and genomic selection

Quantitative trait loci (QTL)/association mapping aims at finding genomic loci associated with the phenotypes, whereas genomic selection focuses on breeding value prediction based on genomic data. Variable selection is a key to both of these tasks as it allows to (1) detect clear mapping signals of QTL activity, and (2) predict the genome-enhanced breeding values accurately. In this paper, we provide an overview of a statistical method called least absolute shrinkage and selection operator (LASSO) and two of its generalizations named elastic net and adaptive LASSO in the contexts of QTL mapping and genomic breeding value prediction in plants (or animals). We also briefly summarize the Bayesian interpretation of LASSO, and the inspired hierarchical Bayesian models. We illustrate the implementation and examine the performance of methods using three public data sets: (1) North American barley data with 127 individuals and 145 markers, (2) a simulated QTLMAS XII data with 5,865 individuals and 6,000 markers for both QTL mapping and genomic selection, and (3) a wheat data with 599 individuals and 1,279 markers only for genomic selection.     

Genetics of Reproductive Isolation in the Drosophila Simulans Clade: Complex Epistasis Underlying Hybrid Male Sterility

We have analyzed the sterility associated with introgressions of the distal one-fourth of the X chromosome from either Drosophila mauritiana or Drosophila sechellia into the genome of Drosophila simulans using a series of visible and DNA markers. Because in Drosophila hybrids, male sterility is usually complete and is often tightly linked with each of several markers used in crosses, a simple genetic basis has generally been assumed. In our low resolution mapping experiment, we were not able to reject the null hypothesis that a single gene, introgressed from either D. mauritiana or D. sechellia, is the cause of male sterility. High resolution mapping, however, reveals a much more complex picture. At least three distinct factors from D. mauritiana, or two from D. sechellia, were identified that need to be jointly present to confer full sterility. Each individual factor by itself is relatively ineffective in causing sterility, or even a partial spermatogenic defect. Moreover, there appear to be more sterility factors on comparable introgressions from D. mauritiana than from D. sechellia. On the basis of these observations, we propose a model which suggests that multilocus weak allele interactions are a very common cause of reproductive incompatibility between closely related species. We also present theoretical argument and empirical evidence against extrapolating the results of within-species analysis to interpret the genetic basis of species differences. The implications of this model on the theories of evolution of species differences and the attempt to understand the mechanisms of hybrid sterility/inviability at the molecular level are discussed.     

Linkage disequilibrium mapping of quantitative trait loci under truncation selection

As a dense map of single nucleotide polymorphism (SNP) markers are available, population-based linkage disequilibrium (LD) mapping or association study is becoming one of the major tools for identifying quantitative trait loci (QTL) and for fine gene mapping. However, in many cases, LD between the marker and trait locus is not very strong. Approaches that maximize the potential of detecting LD will be essential for the success of LD mapping of QTL. In this paper, we propose two strategies for increasing the probability of detecting LD: (1) phenotypic selection and (2) haplotype LD mapping. To provide the foundations for LD mapping of QTL under selection, we develop analytic tools for assessing the impact of phenotypic selection on allele and haplotype frequencies, and LD under three trait models: single trait locus, two unlinked trait loci, and two linked trait loci with or without epistasis. In addition to a traditional chi(2) test, which compares the difference in allele or haplotype frequencies in the selected sample and population sample, we present multiple regression methods for LD mapping of QTL, and investigate which methods are effective in employing phenotypic selection for QTL mapping. We also develop a statistical framework for investigating and comparing the power of the single marker and multilocus haplotype test for LD mapping of QTL. Finally, the proposed methods are applied to mapping QTL influencing variation in systolic blood pressure in an isolated Chinese population.     

Pleiotropic quantitative trait loci contribute to population divergence in traits associated with life-history variation in Mimulus guttatus

Evolutionary biologists seek to understand the genetic basis for multivariate phenotypic divergence. We constructed an F2 mapping population (N = 539) between two distinct populations of Mimulus guttatus. We measured 20 floral, vegetative, and life-history characters on parents and F1 and F2 hybrids in a common garden experiment. We employed multitrait composite interval mapping to determine the number, effect, and degree of pleiotropy in quantitative trait loci (QTL) affecting divergence in floral, vegetative, and life-history characters. We detected 16 QTL affecting floral traits; 7 affecting vegetative traits; and 5 affecting selected floral, vegetative, and life-history traits. Floral and vegetative traits are clearly polygenic. We detected a few major QTL, with all remaining QTL of small effect. Most detected QTL are pleiotropic, implying that the evolutionary shift between these annual and perennial populations is constrained. We also compared the genetic architecture controlling floral trait divergence both within (our intraspecific study) and between species, on the basis of a previously published analysis of M. guttatus and M. nasutus. Eleven of our 16 floral QTL map to approximately the same location in the interspecific map based on shared, collinear markers, implying that there may be a shared genetic basis for floral divergence within and among species of Mimulus.     

The role of pedigree information in combined linkage disequilibrium and linkage mapping of quantitative trait loci in a general complex pedigree

Combined linkage disequilibrium and linkage (LDL) mapping can exploit historical as well as recent and observed recombinations in a recorded pedigree. We investigated the role of pedigree information in LDL mapping and the performance of LDL mapping in general complex pedigrees. We compared using complete and incomplete genotypic data, spanning 5 or 10 generations of known pedigree, and we used bi- or multiallelic markers that were positioned at 1- or 5-cM intervals. Analyses carried out with or without pedigree information were compared. Results were compared with linkage mapping in some of the data sets. Linkage mapping or LDL mapping with sparse marker spacing ( approximately 5 cM) gave a poorer mapping resolution without considering pedigree information compared to that with considering pedigree information. The difference was bigger in a pedigree of more generations. However, LDL mapping with closely linked markers ( approximately 1 cM) gave a much higher mapping resolution regardless of using pedigree information. This study shows that when marker spacing is dense and there is considerable linkage disequilibrium generated from historical recombinations between flanking markers and QTL, the loss of power due to ignoring pedigree information is negligible and mapping resolution is very high.