A model for quantitative trait loci mapping,linkage phase,and segregation pattern estimation for a full-sib progeny

Quantitative trait loci (QTL) mapping is an important approach for the study of the genetic architecture of quantitative traits. For perennial species, inbred lines cannot be obtained due to inbreed depression and a long juvenile period. Instead, linkage mapping can be performed by using a full-sib progeny. This creates a complex scenario because both markers and QTL alleles can have different segregation patterns as well as different linkage phases between them. We present a two-step method for QTL mapping using full-sib progeny based on composite interval mapping (i.e., interval mapping with cofactors), considering an integrated genetic map with markers with different segregation patterns and conditional probabilities obtained by a multipoint approach. The model is based on three orthogonal contrasts to estimate the additive effect (one in each parent) and dominance effect. These estimatives are obtained using the EM algorithm. In the first step, the genome is scanned to detect QTL. After, segregation pattern and linkage phases between QTL and markers are estimated. A simulated example is presented to validate the methodology. In general, the new model is more effective than existing approaches, because it can reveal QTL present in a full-sib progeny that segregates in any pattern present and can also identify dominance effects. Also, the inclusion of cofactors provided more statistical power for QTL mapping.     

Contribution of maternal effect QTL to genetic architecture of early growth in mice

Existing approaches to characterizing quantitative trait loci (QTL) utilize a paradigm explicitly focused on the direct effects of genes, where phenotypic variation among individuals is mapped onto genetic variation of those individuals. For many characters, however, the genotype of the mother via its maternal effect accounts for a considerable portion of the genetically based variation in progeny phenotypes. Thus the focus on direct effect QTL may result in an insufficient or misleading characterization of genetic architecture due to the omission of the potentially important source of genetic variance contributed by maternal effects. We analyze the relative contribution of direct and maternal effect (ME) QTL to early growth in mice using a three-generation intercross of the Small (SM/J) and Large (LG/J) inbred mouse lineages. Using interval mapping and composite interval mapping, direct effect (DE) QTL for early growth (change in body mass during the interval from week 1 to 2) were detected in the F(2) generation of the intercross (n = 510), where no maternal genetic effect variance is present (all individuals are progeny of genetically identical F(1) mothers). ME QTL were detected by treating the phenotypes of cross-fostered F(3) pups as a characteristic of their nurse-dam (n = 168 dams with cross-fostered progeny). Five DE QTL, significant at a chromosome wide level (alpha = 0.05), were detected, with two significant at a genome wide level. FourME QTL significant at the chromosome wide level were detected, with three significant at the genome wide level. A model containing only DE QTL accounted for 11.8% of phenotypic variance, while a model containing only ME QTL accounted for 31.5% of the among litter variance in growth. There was no evidence for pleiotropy of DE and ME loci since there was no overlap between loci detected in these two analyses. Epistasis between all pairs of loci was analyzed for both DEs and MEs. Ten pairs of loci showed significant epistasis for MEs (alpha = 0.05 corrected for multiple comparisons) while four pairs showed significant epistasis for DEs on early growth.     

Genome-wide QTL mapping of saltwater tolerance in sibling species of Anopheles (malaria vector) mosquitoes

Although freshwater (FW) is the ancestral habitat for larval mosquitoes, multiple species independently evolved the ability to survive in saltwater (SW). Here, we use quantitative trait locus (QTL) mapping to investigate the genetic architecture of osmoregulation in Anopheles mosquitoes, vectors of human malaria. We analyzed 1134 backcross progeny from a cross between the obligate FW species An. coluzzii, and its closely related euryhaline sibling species An. merus. Tests of 2387 markers with Bayesian interval mapping and machine learning (random forests) yielded six genomic regions associated with SW tolerance. Overlap in QTL regions from both approaches enhances confidence in QTL identification. Evidence exists for synergistic as well as disruptive epistasis among loci. Intriguingly, one QTL region containing ion transporters spans the 2Rop chromosomal inversion that distinguishes these species. Rather than a simple trait controlled by one or a few loci, our data are most consistent with a complex, polygenic mode of inheritance.     

The role of surface characteristics in eliciting humoral encapsulation of foreign bodies in Plasmodium-refractory and -susceptible strains of Anopheles gambiae

A refractory strain of the mosquito, Anopheles gambiae, melanotically encapsulates and kills many species of malaria parasites, whereas susceptible strains allow the parasites to develop normally. To study the role of surface characteristics in eliciting this immune response, 27 types of chromatography beads that differed in matrix type, charge, functional group, and functional group density were assayed for degree of melanotic encapsulation in refractory and susceptible mosquitoes. Overall, two glucan-based matrices, Sephadex (dextran) and cellulose, stimulated the strongest responses, regardless of functional group. Substituting matrix hydroxyl groups with functional groups on Sephadex and cellulose beads decreased the level of encapsulation. These results demonstrate that glucans induce melanotic encapsulation in An. gambiae. Beads with agarose, polystyrene, and acrylic matrices, and most methacrylate-based beads elicited little or no melanization; however, epoxide-methacrylate beads were encapsulated, demonstrating that glucans are not essential for eliciting a response. Comparisons between the two strains demonstrated that refractory mosquitoes melanized many bead types to a greater degree than did susceptible mosquitoes. On this basis, we propose that an important difference between the two strains is that one of the enzymes involved in the melanization pathway functions at a higher level in the refractory strain. Finally, of all beads tested, only 85% substituted CM-Sephadex beads were virtually unmelanized in susceptible mosquitoes but highly melanized in the refractory strain; thus, a specific surface microenvironment is necessary to demonstrate this effect.     

Mapping of quantitative trait loci controlling adaptive traits in coastal Douglas fir. III. Quantitative trait loci-by-environment interactions

Quantitative trait loci (QTL) were mapped in the woody perennial Douglas fir (Pseudotsuga menziesii var. menziesii [Mirb.] Franco) for complex traits controlling the timing of growth initiation and growth cessation. QTL were estimated under controlled environmental conditions to identify QTL interactions with photoperiod, moisture stress, winter chilling, and spring temperatures. A three-generation mapping population of 460 cloned progeny was used for genetic mapping and phenotypic evaluations. An all-marker interval mapping method was used for scanning the genome for the presence of QTL and single-factor ANOVA was used for estimating QTL-by-environment interactions. A modest number of QTL were detected per trait, with individual QTL explaining up to 9.5% of the phenotypic variation. Two QTL-by-treatment interactions were found for growth initiation, whereas several QTL-by-treatment interactions were detected among growth cessation traits. This is the first report of QTL interactions with specific environmental signals in forest trees and will assist in the identification of candidate genes controlling these important adaptive traits in perennial plants.     

Mapping of shoot fly tolerance loci in sorghum using SSR markers

Sorghum (Sorghum bicolor (L.) Moench) is one of the most important crops in the semiarid regions of the world. One of the important biotic constraints to sorghum production in India is the shoot fly which attacks sorghum at the seedling stage. Identification of the genomic regions containing quantitative trait loci (QTLs) for resistance to shoot fly and the linked markers can facilitate sorghum improvement programmes through marker-assisted selection. A simple sequence repeat (SSR) marker- based skeleton linkage map of two linkage groups of sorghum was constructed in a population of 135 recombinant inbred lines (RIL) derived from a cross between IS18551 (resistant to shoot fly) and 296B (susceptible to shoot fly). A total of 14 SSR markers, seven each on linkage groups A and C were mapped. Using data of different shoot fly resistance component traits, one QTL which is common for glossiness, oviposition and dead hearts was detected following composite interval mapping (CIM) on linkage group A. The phenotypic variation explained by this QTL ranged from 3.8%–6.3%. Besides the QTL detected by CIM, two more QTLs were detected following multi-trait composite interval mapping (MCIM), one each on linkage groups A and C for the combinations of traits which were correlated with each other. Results of the present study are novel as we could find out the QTLs governing more than one trait (pleiotropic QTLs). The identification of pleiotropic QTLs will help in improvement of more than one trait at a time with the help of the same linked markers. For all the QTLs, the resistant parent IS18551 contributed resistant alleles.     

Targeted identification of markers linked to malaria and filarioid nematode parasite resistance genes in the mosquito Aedes aegypti.

Quantitative trait loci (QTL) have been identified for competence of the mosquito Aedes aegypti to transmit the avian malaria parasite Plasmodium gallinaceum and the human filarial parasite Brugia malayi. Efforts towards the map-based cloning of the associated genes are limited by the availability of genetic markers for fine-scale mapping of the QTL positions. Two F2 mosquito populations were subjected to bulked segregant analysis to identify random amplified polymorphic DNA (RAPD)-PCR fragments linked with the major QTL determining susceptibility to both parasites. Individual mosquitoes for the bulks were selected on the basis of their genotypes at restriction fragment length polymorphism (RFLP) loci tightly linked with the QTL. Pool-positive RAPD fragments were cloned and evaluated as RFLP markers. Of the 62 RAPD/RFLP fragments examined, 10 represented low-copy number sequences. Five of these clones were linked with the major QTL for P. gallinaceum susceptibility (pgs1), of which one clone mapped within the flanking markers that define the QTL interval. The remaining five clones were linked with the major QTL for B. malayi susceptibility (fsb1), and again one clone mapped within the flanking markers that define the QTL interval. In addition, nine RAPD/RFLP fragments were isolated that seem to be of non-mosquito origin.     

Genetic basis of adaptation: flowering time in Arabidopsis thaliana

 We have mapped QTLs (quantitative trait loci) for an adaptive trait, flowering time, in a selfing annual, Arabidopsis thaliana. To obtain a mapping population we made a cross between an early-summer, annual strain, Li-5, and an individual from a late over-wintering natural population, Naantali. From the backcross to Li-5 298 progeny were grown, of which 93 of the most extreme individuals were genotyped. The data were analysed with both interval mapping and composite interval mapping methods to reveal one major and six minor QTLs, with at least one QTL on each of the five chromosomes. The QTL on chromosome 4 was a major one with an effect of 17.3 days on flowering time and explaining 53.4% of the total variance. The others had effects of at most 6.5 days, and they accounted for only small portions of the variance. Epistasis was indicated between one pair of the QTLs. The result of finding one major QTL and little epistasis agrees with previous studies on flowering time in Arabidopsis thaliana and other species. That several QTLs were found was expected considering the large number of possible candidate loci. In the light of the suggested genetic models of gene action at the candidate loci, epistasis was to be expected. The data showed that major QTLs for adaptive traits can be detected in non-domesticated species. Received: 15 January 1997/Accepted: 21 February 1997     

Quantitative trait locus mapping of pyrethroid resistance in Colorado potato beetle, Leptinotarsa decemlineata (Say) (Coleoptera: Chrysomelidae)

Quantitative trait locus (QTL) mapping is a valuable new tool for locating genomic regions that underlie variation in important traits such as insecticide resistance. Because QTL mapping complements a candidate gene strategy for understanding the genetic architecture of important traits, it may also facilitate the identification of genes causing important variation. After mapping the QTL locations, markers closely linked to QTL can be used for genetic analysis of population structure and to measure the spread and increase of resistance-causing QTL alleles. In this study, QTL influencing resistance to the pyrethroid insecticide esfenvalerate were mapped in the Colorado potato beetle Leptinotarsa decemlineata (Say) (CPB). Three QTL contributing to esfenvalerate resistance were identified from a mapping population of 79 individuals analyzed at 90 marker loci. One QTL had a large effect and two QTL had smaller effects. The major QTL occurs on the X chromosome, overlapping the position of a candidate gene (Leptinotarsa decemlineata Voltage sensitive sodium channel [LdVssc1]) previously implicated in pyrethroid resistance. Resistance-increasing alleles at the two minor-effect QTL originated with the susceptible parent, suggesting that alleles of small effect may be segregating in susceptible populations. Comparison of the New York population from which the susceptible parent originated with a more-susceptible population from North Carolina suggests that the minor-effect loci identified here may explain some of the variation in tolerance observed among susceptible populations. DNA sequencing of a portion of LdVssc1 shows that the resistance-conferring allele from the resistant parent does not contain the kdr mutation previously found in CPB and typically observed in other insects that are resistant to pyrethroid insecticides because of changes in this gene.     

Multi-environment mapping and meta-analysis of 100-seed weight in soybean

100-Seed weight (100-SW) of soybean is an important but complicated quantitative trait to yield. This study was focus on the quantitative trait loci (QTLs) of soybean 100-SW from 2006 to 2010, using recombination inbred lines population that was derived from a cross between Charleston and Dongnong 594. A total of 23 QTLs for 100-SW were detected in the linkage group C2, D1a, F, G and O. Nine QTLs were identified by composite interval mapping including one QTL with the minimum confidence interval (CI) of 1.3?cM, while 14 QTLs by multiple interval mapping. Furthermore, 94 reported QTLs of 100-SW were integrated with our QTL mapping results using BioMercator. As a result, 15 consensus QTLs and their corresponding markers were identified. The minimum CI was reduced to 1.52?cM by the combination of meta-analysis. These findings may merit fine-mapping of these QTL in soybean.     

Identification of a major quantitative trait locus conditioning resistance to greenbug biotype E in sorghum PI 550610 using simple sequence repeat markers

Greenbug, Schizaphis graminum (Rondani), represents the most important pest insect of sorghum, Sorghum bicolor (L.) Moench, in the Great Plains of the United States. Biotype E is the most widespread and dominant type not only in sorghum and wheat, Triticum aestivum L., fields, but also on many noncultivated grass species. This study was designed to determine sorghum accession PI 550610 resistance to greenbug biotype E, to map the resistance quantitative trait loci (QTLs) by using an established simple sequence repeat (SSR) linkage map and to identify SSR markers closely linked to the major resistance QTLs. In greenhouse screening tests, seedlings of PI 550610 showed strong resistance to the greenbug at a level similar to resistant accession PI550607. For QTL mapping, one F2 population containing 277 progeny and one population containing 233 F2:3 families derived from Westland A line x PI 550610 were used to genotype 132 polymorphic SSR markers and to phenotype seedling resistance to greenbug feeding. Phenotypic evaluation of sorghum seedling damage at 7, 12, 17, and 21 d postinfestation in the F2:3 families revealed that resistance variation was normally distributed. Single marker analysis indicated 16 SSRs spread over five chromosomes were significant for greenbug resistance. Composite interval and multiple interval mapping procedures indicated that a major QTL resided in the interval of 6.8 cM between SSR markers Xtxp358 and Xtxp289 on SBI-09. The results will be valuable in the development of new greenbug biotype E resistant sorghum cultivars and for the further characterization of major genes by map-based cloning.     

Quantitative trait loci for the circadian clock in Neurospora crassa

Neurospora crassa has been a model organism for the study of circadian clocks for the past four decades. Among natural accessions of Neurospora crassa, there is significant variation in clock phenotypes. In an attempt to investigate natural allelic variants contributing to quantitative variation, we used a quantitative trait loci mapping approach to analyze three independent mapping populations whose progenitors were collected from geographically isolated locations. Two circadian clock phenotypes, free-running period and entrained phase, were evaluated in the 188 F(1) progeny of each mapping population. To identify the clock QTL, we applied two QTL mapping analyses: composite interval mapping (CIM) and Bayesian multiple QTL analysis (BMQ). When controlling false positive rates < or =0.05, BMQ appears to be the more sensitive of the two approaches. BMQ confirmed most of the QTL from CIM (18 QTL) and identified 23 additional QTL. While 13 QTL colocalize with previously identified clock genes, we identified 30 QTL that were not linked with any previously characterized clock genes. These are candidate regions where clock genes may be located and are expected to lead to new insights in clock regulation.     

Quantitative trait loci that control vector competence for dengue-2 virus in the mosquito Aedes aegypti

Quantitative trait loci (QTL) affecting the ability of the mosquito Aedes aegypti to become infected with dengue-2 virus were mapped in an F(1) intercross. Dengue-susceptible A. aegypti aegypti were crossed with dengue refractory A. aegypti formosus. F(2) offspring were analyzed for midgut infection and escape barriers. In P(1) and F(1) parents and in 207 F(2) individuals, regions of 14 cDNA loci were analyzed with single-strand conformation polymorphism analysis to identify and orient linkage groups with respect to chromosomes I-III. Genotypes were also scored at 57 RAPD-SSCP loci, 5 (TAG)(n) microsatellite loci, and 6 sequence-tagged RAPD loci. Dengue infection phenotypes were scored in 86 F(2) females. Two QTL for a midgut infection barrier were detected with standard and composite interval mapping on chromosomes II and III that accounted for approximately 30% of the phenotypic variance (sigma(2)(p)) in dengue infection and these accounted for 44 and 56%, respectively, of the overall genetic variance (sigma(2)(g)). QTL of minor effect were detected on chromosomes I and III, but these were not detected with composite interval mapping. Evidence for a QTL for midgut escape barrier was detected with standard interval mapping but not with composite interval mapping on chromosome III.     

Quantitative Trait Loci Mapping of Leaf Morphological Traits and Chlorophyll Content in Cultivated Tetraploid Cotton

Genetic mapping provides a powerful tool for quantitative trait loci （QTL） analysis at the molecular level. A simple sequence repeat （SSR） genetic map containing 590 markers and a BCI population from two cultivated tetraploid cotton （Gossypium hirsutum L.） cultivars, namely TM-1 and Hai 7124 （G. barbadense L.）, were used to map and analyze QTL using the composite interval mapping （CIM） method. Thirty one QTLs, 10 for lobe length, 13 for lobe width, six for lobe angle, and two for leaf chlorophyll content, were detected on 15 chromosomes or linkage groups at logarithm of odds （LOD）≥2.0, of which 15 were found for leaf morphology at LOD≥3.0. The genetic effects of the QTL were estimated. These results are fundamental for marker-assisted selection （MAS） of these traits in tetraploid cotton breeding.     

Interval Mapping of Multiple Quantitative Trait Loci

The interval mapping method is widely used for the mapping of quantitative trait loci (QTLs) in segregating generations derived from crosses between inbred lines. The efficiency of detecting and the accuracy of mapping multiple QTLs by using genetic markers are much increased by employing multiple QTL models instead of the single QTL models (and no QTL models) used in interval mapping. However, the computational work involved with multiple QTL models is considerable when the number of QTLs is large. In this paper it is proposed to combine multiple linear regression methods with conventional interval mapping. This is achieved by fitting one QTL at a time in a given interval and simultaneously using (part of) the markers as cofactors to eliminate the effects of additional QTLs. It is shown that the proposed method combines the easy computation of the single QTL interval mapping method with much of the efficiency and accuracy of multiple QTL models.     

Quantitative genetics of vector competence for La Crosse virus and body size in Ochlerotatus hendersoni and Ochlerotatus triseriatus interspecific hybrids

La Crosse virus is a leading cause of pediatric encephalitis in the United States. The mosquito Ochlerotatus triseriatus is an efficient vector for La Crosse virus, whereas the closely related O. hendersoni transmits only at very low rates. Quantitative trait loci (QTL) affecting the ability to orally transmit this virus and adult body size were identified in 164 F(2) female individuals from interspecific crosses of O. hendersoni females and O. triseriatus males using a combination of composite interval mapping (CIM), interval mapping (IM) for binary traits, and single-marker mapping. For oral transmission (OT), no genome locations exceeded the 95% experimentwise threshold for declaring a QTL using IM, but single-marker analysis identified four independent regions significantly associated with OT that we considered as tentative QTL. With two QTL, an increase in OT was associated with alleles from the refractory vector, O. hendersoni, and likely reflect epistatic interactions between genes that were uncovered by our interspecific crosses. For body size, two QTL were identified using CIM and a third tentative QTL was identified using single-marker analysis. The genome regions associated with body size also contain three QTL controlling OT, suggesting that these regions contain either single genes with pleiotropic effects or multiple linked genes independently determining each trait.     

Quantitative trait loci in <Emphasis Type="Italic">Anopheles gambiae</Emphasis> controlling the encapsulation response against <Emphasis Type="Italic">Plasmodium cynomolgi</Emphasis> Ceylon

Background  

Anopheles gambiae females are the world's most successful vectors of human malaria. However, a fraction of these mosquitoes is refractory to Plasmodium development. L3-5, a laboratory selected refractory strain, encapsulates transforming ookinetes/early oocysts of a wide variety of Plasmodium species. Previous studies on these mosquitoes showed that one major (Pen1) and two minor (Pen2, Pen3) autosomal dominant quantitative trait loci (QTLs) control the melanotic encapsulation response against P. cynomolgi B, a simian malaria originating in Malaysia.     

Mapping of QTLs governing agronomic and yield traits in chickpea

Chickpea is one of the most important leguminous cool season food crops, cultivated prevalently in South Asia and Middle East. The main objective of this study was to identify quantitative trait loci (QTLs) associated with seven agronomic and yield traits in two recombinant inbred line populations of chickpea derived from the crosses JG62 × Vijay (JV population) and Vijay × ICC4958 (VI population) from at least three environments. Single locus QTL analysis involved composite interval mapping (CIM) for individual traits and multiple-trait composite interval mapping (MCIM) for correlated traits to detect pleiotropic QTLs. Two-locus analysis was conducted to identify the main effect QTLs (M-QTLs), epistatic QTLs (E-QTLs) and QTL × environment interactions. Through CIM analysis, a total of 106 significant QTLs (41 in JV and 65 in VI populations) were identified for the seven traits, of which one QTL each for plant height and days to maturity was common in both the populations. Six pleiotropic QTLs that were consistent over the environments were also identified. LG2 in JV and LG1a in VI contained at least one QTL for each trait. Hence, concentrating on these LGs in molecular breeding programs is most likely to bring simultaneous improvement in these traits.     

A new QTL for resistance to Fusarium ear rot in maize

Understanding the inheritance of resistance to Fusarium ear rot is a basic prerequisite for an efficient resistance breeding in maize. In this study, 250 recombinant inbred lines (RILs) along with their resistant (BT-1) and susceptible (N6) parents were planted in Zhengzhou with three replications in 2007 and 2008. Each line was artificially inoculated using the nail-punch method. Significant genotypic variation in response to Fusarium ear rot was detected in both years. Based on a genetic map containing 207 polymorphic simple sequence repeat (SSR) markers with average genetic distances of 8.83?cM, the ear rot resistance quantitative trait loci (QTL) were analyzed by composite interval mapping with a mixed model (MCIM) across the environments. In total, four QTL were detected on chromosomes 3, 4, 5, and 6. The resistance allele at each of these four QTL was contributed by resistant parent BT-1, and accounted for 2.5-10.2% of the phenotypic variation. However, no significant epistasis interaction effect was detected after a two-dimensional genome scan. Among the four QTL, one QTL with the largest effect on chromosome 4 (bin 4.06) can be suggested to be a new locus for resistance to Fusarium ear rot, which broadens the genetic base for resistance to the disease and can be used for further genetic improvement in maize-breeding programs.     

Multiple Trait Analysis of Genetic Mapping for Quantitative Trait Loci

We present in this paper models and statistical methods for performing multiple trait analysis on mapping quantitative trait loci (QTL) based on the composite interval mapping method. By taking into account the correlated structure of multiple traits, this joint analysis has several advantages, compared with separate analyses, for mapping QTL, including the expected improvement on the statistical power of the test for QTL and on the precision of parameter estimation. Also this joint analysis provides formal procedures to test a number of biologically interesting hypotheses concerning the nature of genetic correlations between different traits. Among the testing procedures considered are those for joint mapping, pleiotropy, QTL by environment interaction, and pleiotropy vs. close linkage. The test of pleiotropy (one pleiotropic QTL at a genome position) vs. close linkage (multiple nearby nonpleiotropic QTL) can have important implications for our understanding of the nature of genetic correlations between different traits in certain regions of a genome and also for practical applications in animal and plant breeding because one of the major goals in breeding is to break unfavorable linkage. Results of extensive simulation studies are presented to illustrate various properties of the analyses.