Identification of Molecular Markers Associated with Verticillium Wilt Resistance in Alfalfa (Medicago Sativa L.) Using High-Resolution Melting

Verticillium wilt, caused by the soilborne fungus, Verticillium alfalfae, is one of the most serious diseases of alfalfa (Medicago sativa L.) worldwide. To identify loci associated with resistance to Verticillium wilt, a bulk segregant analysis was conducted in susceptible or resistant pools constructed from 13 synthetic alfalfa populations, followed by association mapping in two F1 populations consisted of 352 individuals. Simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers were used for genotyping. Phenotyping was done by manual inoculation of the pathogen to replicated cloned plants of each individual and disease severity was scored using a standard scale. Marker-trait association was analyzed by TASSEL. Seventeen SNP markers significantly associated with Verticillium wilt resistance were identified and they were located on chromosomes 1, 2, 4, 7 and 8. SNP markers identified on chromosomes 2, 4 and 7 co-locate with regions of Verticillium wilt resistance loci reported in M. truncatula. Additional markers identified on chromosomes 1 and 8 located the regions where no Verticillium resistance locus has been reported. This study highlights the value of SNP genotyping by high resolution melting to identify the disease resistance loci in tetraploid alfalfa. With further validation, the markers identified in this study could be used for improving resistance to Verticillium wilt in alfalfa breeding programs.     

Genotyping‐by‐sequencing‐based genome‐wide association studies on Verticillium wilt resistance in autotetraploid alfalfa (Medicago sativa L.)

Verticillium wilt (VW) is a fungal disease that causes severe yield losses in alfalfa. The most effective method to control the disease is through the development and use of resistant varieties. The identification of marker loci linked to VW resistance can facilitate breeding for disease‐resistant alfalfa. In the present investigation, we applied an integrated framework of genome‐wide association with genotyping‐by‐sequencing (GBS) to identify VW resistance loci in a panel of elite alfalfa breeding lines. Phenotyping was performed by manual inoculation of the pathogen to healthy seedlings, and scoring for disease resistance was carried out according to the standard test of the North America Alfalfa Improvement Conference (NAAIC). Marker–trait association by linkage disequilibrium identified 10 single nucleotide polymorphism (SNP) markers significantly associated with VW resistance. Alignment of the SNP marker sequences to the M. truncatula genome revealed multiple quantitative trait loci (QTLs). Three, two, one and five markers were located on chromosomes 5, 6, 7 and 8, respectively. Resistance loci found on chromosomes 7 and 8 in the present study co‐localized with the QTLs reported previously. A pairwise alignment (blastn ) using the flanking sequences of the resistance loci against the M. truncatula genome identified potential candidate genes with putative disease resistance function. With further investigation, these markers may be implemented into breeding programmes using marker‐assisted selection, ultimately leading to improved VW resistance in alfalfa.     

Quantitative trait locus analysis of Verticillium wilt resistance in an introgressed recombinant inbred population of Upland cotton

Verticillium wilt (VW) of Upland cotton (Gossypium hirsutum L.) is caused by the soil-borne fungal pathogen Verticillium dahlia Kleb. The availability of VW-resistant cultivars is vital for control of this economically important disease, but there is a paucity of Upland cotton breeding lines and cultivars with a high level of resistance to VW. In general, G. barbadense L. (source of Pima cotton) is more VW-resistant than Upland cotton. However, the transfer of VW resistance from G. barbadense to Upland cotton is challenging because of hybrid breakdown in the F2 and successive generations of interspecific populations. We conducted two replicated greenhouse studies (tests 1 and 2) to assess the heritability of VW resistance to a defoliating V. dahliae isolate and identify genetic markers associated with VW resistance in an Upland cotton recombinant inbred mapping population that has stable introgression from Pima cotton. Disease ratings at the seedling stage on several different days after the first inoculation (DAI) in test 1, as well as the percentages of infected and defoliated leaves at 2 DAI in test 2, were found to be low to moderately heritable, indicating the importance of a replicated progeny test in selection for VW resistance. With a newly constructed linkage map consisting of 882 simple sequence repeat, single nucleotide polymorphism, and resistance gene analog–amplified fragment length polymorphism marker loci, we identified a total of 21 quantitative trait loci (QTLs) on 11 chromosomes and two linkage groups associated with VW resistance at several different DAIs in greenhouse tests 1 and 2. The markers associated with the VW resistance QTLs will facilitate fine mapping and cloning of VW resistance genes and genomics-assisted breeding for VW-resistant cultivars.     

Evaluation of linkage disequilibrium and its effect on non-parametric multipoint linkage analysis using two high density single-nucleotide polymorphism mapping panels

Genotype data from the Illumina Linkage III SNP panel (n = 4,720 SNPs) and the Affymetrix 10 k mapping array (n = 11,120 SNPs) were used to test the effects of linkage disequilibrium (LD) between SNPs in a linkage analysis in the Collaborative Study on the Genetics of Alcoholism pedigree collection (143 pedigrees; 1,614 individuals). The average r2 between adjacent markers across the genetic map was 0.099 +/- 0.003 in the Illumina III panel and 0.17 +/- 0.003 in the Affymetrix 10 k array. In order to determine the effect of LD between marker loci in a nonparametric multipoint linkage analysis, markers in strong LD with another marker (r2 > 0.40) were removed (n = 471 loci in the Illumina panel; n = 1,804 loci in the Affymetrix panel) and the linkage analysis results were compared to the results using the entire marker sets. In all analyses using the ALDX1 phenotype, 8 linkage regions on 5 chromosomes (2, 7, 10, 11, X) were detected (peak markers p < 0.01), and the Illumina panel detected an additional region on chromosome 6. Analysis of the same pedigree set and ALDX1 phenotype using short tandem repeat markers (STRs) resulted in 3 linkage regions on 3 chromosomes (peak markers p < 0.01). These results suggest that in this pedigree set, LD between loci with spacing similar to the SNP panels tested may not significantly affect the overall detection of linkage regions in a genome scan. Moreover, since the data quality and information content are greatly improved in the SNP panels over STR genotyping methods, new linkage regions may be identified due to higher information content and data quality in a dense SNP linkage panel.     

Linkage disequilibrium fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

Association mapping of kernel size and milling quality in wheat (Triticum aestivum L.) cultivars

Association mapping is a method for detection of gene effects based on linkage disequilibrium (LD) that complements QTL analysis in the development of tools for molecular plant breeding. In this study, association mapping was performed on a selected sample of 95 cultivars of soft winter wheat. Population structure was estimated on the basis of 36 unlinked simple-sequence repeat (SSR) markers. The extent of LD was estimated on chromosomes 2D and part of 5A, relative to the LD observed among unlinked markers. Consistent LD on chromosome 2D was <1 cM, whereas in the centromeric region of 5A, LD extended for approximately 5 cM. Association of 62 SSR loci on chromosomes 2D, 5A, and 5B with kernel morphology and milling quality was analyzed through a mixed-effects model, where subpopulation was considered as a random factor and the marker tested was considered as a fixed factor. Permutations were used to adjust the threshold of significance for multiple testing within chromosomes. In agreement with previous QTL analysis, significant markers for kernel size were detected on the three chromosomes tested, and alleles potentially useful for selection were identified. Our results demonstrated that association mapping could complement and enhance previous QTL information for marker-assisted selection.     

陆地棉品种和骨干品系黄萎病抗性鉴定

选育和推广抗病品种是防治陆地棉黄萎病的主要措施,为了早日实现多类型、多区域大面积抗病品种的应用,本研究选取107份遗传背景差异较大的种质,利用河北省农林科学院棉花研究所小安舍试验站黄萎病病圃进行了3年黄萎病抗性重复鉴定。鉴定得到抗病品系8个,占7.5%;耐病品种(系)20个,占18.7%。本研究表明,当前被作为育种亲本的抗病品系还太少,需要深入开展抗病遗传机制,以及与其他经济性状协同改良的关系,为陆地棉抗病育种提供理论指导;达到抗病或接近抗病水平的大部分品种(系)来自于海陆野远缘后代,具有外源基因血统,证明了远缘杂交是陆地棉黄萎病抗性改良的有效手段。     

From the editor’s desk

Improving cotton fibre quality is a major breeding goal for Upland cotton in China. To investigate the genetic mechanisms of fibre quality, a diverse panel of 403 Upland cotton accessions was grown, and the fibre quality traits were measured in six different environments. Genotyping was performed with genomewide simple sequence repeats. A total of 201 markers were polymorphic and generated 394 allele loci, and 403 accessions were arranged into two subgroups using Structure software. Of the marker loci, 18.94% showed significant linkage disequilibrium (\(P< 0.05\)). A mixed linear model in association mapping showed that 51 associations were significant between 39 polymorphic loci and five fibre quality traits, according to best linear unbiased prediction, and in at least three of six environments. Of the 39 associated marker loci, 12 were coincident with previous studies. There were 41 typical accessions identified as containing favourable allele loci related to fibre quality traits. The identified favourable QTL alleles and typical accessions for fibre quality are excellent genetic resources for future cotton breeding in China.     

Genome-wide association mapping of blood cell traits in mice

Genetic variations in blood cell parameters can impact clinical traits. We report here the mapping of blood cell traits in a panel of 100 inbred strains of mice of the Hybrid Mouse Diversity Panel (HMDP) using genome-wide association (GWA). We replicated a locus previously identified in using linkage analysis in several genetic crosses for mean corpuscular volume (MCV) and a number of other red blood cell traits on distal chromosome 7. Our peak for SNP association to MCV occurred in a linkage disequilibrium (LD) block spanning from 109.38 to 111.75 Mb that includes Hbb-b1, the likely causal gene. Altogether, we identified five loci controlling red blood cell traits (on chromosomes 1, 7, 11, 12, and 16), and four of these correspond to loci for red blood cell traits reported in a recent human GWA study. For white blood cells, including granulocytes, monocytes, and lymphocytes, a total of six significant loci were identified on chromosomes 1, 6, 8, 11, 12, and 15. An average of ten candidate genes were found at each locus and those were prioritized by examining functional variants in the HMDP such as missense and expression variants. These results provide intermediate phenotypes and candidate loci for genetic studies of atherosclerosis and cancer as well as inflammatory and immune disorders in mice.     

Genetic analysis of Verticillium wilt resistance in a backcross inbred line population and a meta-analysis of quantitative trait loci for disease resistance in cotton

Background

Verticillium wilt (VW) and Fusarium wilt (FW), caused by the soil-borne fungi Verticillium dahliae and Fusarium oxysporum f. sp. vasinfectum, respectively, are two most destructive diseases in cotton production worldwide. Root-knot nematodes (Meloidogyne incognita, RKN) and reniform nematodes (Rotylenchulus reniformis, RN) cause the highest yield loss in the U.S. Planting disease resistant cultivars is the most cost effective control method. Numerous studies have reported mapping of quantitative trait loci (QTLs) for disease resistance in cotton; however, very few reliable QTLs were identified for use in genomic research and breeding.

Results

This study first performed a 4-year replicated test of a backcross inbred line (BIL) population for VW resistance, and 10 resistance QTLs were mapped based on a 2895 cM linkage map with 392 SSR markers. The 10 VW QTLs were then placed to a consensus linkage map with other 182 VW QTLs, 75 RKN QTLs, 27 FW QTLs, and 7 RN QTLs reported from 32 publications. A meta-analysis of QTLs identified 28 QTL clusters including 13, 8 and 3 QTL hotspots for resistance to VW, RKN and FW, respectively. The number of QTLs and QTL clusters on chromosomes especially in the A-subgenome was significantly correlated with the number of nucleotide-binding site (NBS) genes, and the distribution of QTLs between homeologous A- and D- subgenome chromosomes was also significantly correlated.

Conclusions

Ten VW resistance QTL identified in a 4-year replicated study have added useful information to the understanding of the genetic basis of VW resistance in cotton. Twenty-eight disease resistance QTL clusters and 24 hotspots identified from a total of 306 QTLs and linked SSR markers provide important information for marker-assisted selection and high resolution mapping of resistance QTLs and genes. The non-overlapping of most resistance QTL hotspots for different diseases indicates that their resistances are controlled by different genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1682-2) contains supplementary material, which is available to authorized users.     

Identification of Associated SSR Markers for Yield Component and Fiber Quality Traits Based on Frame Map and Upland Cotton Collections

Detecting QTLs (quantitative trait loci) that enhance cotton yield and fiber quality traits and accelerate breeding has been the focus of many cotton breeders. In the present study, 359 SSR (simple sequence repeat) markers were used for the association mapping of 241 Upland cotton collections. A total of 333 markers, representing 733 polymorphic loci, were detected. The average linkage disequilibrium (LD) decay distances were 8.58 cM (r² > 0.1) and 5.76 cM (r² > 0.2). 241 collections were arranged into two subgroups using STRUCTURE software. Mixed linear modeling (MLM) methods (with population structure (Q) and relative kinship matrix (K)) were applied to analyze four phenotypic datasets obtained from four environments (two different locations and two years). Forty-six markers associated with the number of bolls per plant (NB), boll weight (BW), lint percentage (LP), fiber length (FL), fiber strength (FS) and fiber micornaire value (FM) were repeatedly detected in at least two environments. Of 46 associated markers, 32 were identified as new association markers, and 14 had been previously reported in the literature. Nine association markers were near QTLs (at a distance of less than 1–2 LD decay on the reference map) that had been previously described. These results provide new useful markers for marker-assisted selection in breeding programs and new insights for understanding the genetic basis of Upland cotton yields and fiber quality traits at the whole-genome level.     

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.     

Significant evidence for linkage disequilibrium over a 5-cM region among Afrikaners

We explore the extent of deviations from Hardy-Weinberg equilibrium (HWE) at a marker locus and linkage disequilibrium (LD) between pairs of marker loci in the Afrikaner population of South Africa. DNA samples were used for genotyping of 23 loci on six chromosomes. The samples were collected from 91 healthy unrelated Afrikaner adults. Exact tests were used to determine evidence for deviations from HWE at a single marker locus or LD between pairs of marker loci. At the 0.05 level of significance, evidence was found for deviation from HWE at only one of the 23 loci. At the same level of significance, LD was found among 8 of the 34 intrachromosomal pairs of loci. On chromosome 21, there was evidence for LD (P = 0.02) between a pair of loci with a genetic distance of 5.51 cM. On chromosome 2, there was evidence for LD between a pair of loci with a genetic distance of 5.28 cM (P = 0.002) and a pair of loci with a genetic distance of 3.68 cM (P = 0.0004). Detailed analysis of LD for one locus pair indicated that only a few of all alleles participated in the LD and that strong LD was most often positive. Our findings indicate that Afrikaans-speaking Afrikaners represent one of those special populations deemed particularly suitable for disequilibrium mapping.     

Wide-cross whole-genome radiation hybrid mapping of the cotton (Gossypium barbadense L.) genome

Whole-genome radiation hybrid mapping has been applied extensively to human and certain animal species, but little to plants. We recently demonstrated an alternative mapping approach in cotton (Gossypium hirsutum L.), based on segmentation by 5-krad γ-irradiation and derivation of wide-cross whole-genome radiation hybrids (WWRHs). However, limitations observed at the 5-krad level suggested that higher doses might be advantageous. Here, we describe the development of an improved second-generation WWRH panel after higher dose irradiation and compare the resulting map to the 5-krad map. The genome of G. hirsutum (n=26) was used to rescue the radiation-segmented genome of G. barbadense (n=26) introduced via 8- and 12-krad γ-irradiated pollen. Viable seedlings were not recovered after 12-krad irradiation, but 8-krad irradiation permitted plant recovery and construction of a 92-member WWRH mapping panel. Assessment of 31 SSR marker loci from four chromosomes revealed that the 8-krad panel has a marker retention frequency of ca. 76%, which is approximately equivalent to the rate of loss in a low-dose animal radiation hybrid panel. Retention frequencies of loci did not depart significantly from independence when compared between the A and D subgenomes, or according to positions along individual chromosomes. WWRH maps of chromosomes 10 and 17 were generated by the maximum likelihood RHMAP program and the general retention model. The resulting maps bolster evidence that WWRH mapping complements traditional linkage mapping and works in cotton, and that the 8-krad panel complements the 5-krad panel by offering higher rates of chromosome breakages, lower marker retention frequency, and more retention patterns. Electronic Supplementary Material Supplementary material is available for this article at     

Extent and genome-wide distribution of linkage disequilibrium in commercial maize germplasm

Association mapping is based on linkage disequilibrium (LD) resulting from historical recombinations and helps understanding the genetic basis of complex traits. Many factors affect LD and, therefore, it must be determined empirically in the germplasm under investigation to examine the prospects of successful genome-wide association mapping. The objectives of our study were to (1) examine the extent of LD with simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers in 1,537 commercial maize inbred lines belonging to four heterotic pools, (2) compare the LD patterns determined by these two marker types, (3) evaluate the number of SNP markers needed to perform genome-wide association analyses, and (4) investigate temporal trends of LD. Mean values of the squared correlation coefficient ( $ \bar{R} $ ) were almost identical for unlinked, linked, and adjacent SSR marker pairs. In contrast, $ \bar{R} $ values were lowest for the unlinked SNP loci and highest for the SNPs within amplicons. LD decay varied across the different heterotic pools and the individual chromosomes. The SSR markers employed in the present study are not adequate for association analysis, because of insufficient marker density for the germplasm evaluated. Based on the decay of LD in the various heterotic pools, we would need between 4,000 and 65,000 SNP markers to detect with a reasonable power associations with rather large quantitative trait loci (QTL). A much higher marker density is required to identify QTL with smaller effects. However, not only the total number of markers but also their distribution among and along the chromosomes are primordial for undertaking powerful association analyses.     

Cotton genome mapping with new microsatellites from Acala ‘Maxxa’ BAC-ends

Fine mapping and positional cloning will eventually improve with the anchoring of additional markers derived from genomic clones such as BACs. From 2,603 new BAC-end genomic sequences from Gossypium hirsutum Acala ‘Maxxa’, 1,316 PCR primer pairs (designated as MUSB) were designed to flank microsatellite or simple sequence repeat motif sequences. Most (1164 or 88%) MUSB primer pairs successfully amplified DNA from three species of cotton with an average of three amplicons per marker and 365 markers (21%) were polymorphic between G. hirsutum and G. barbadense. An interspecific RIL population developed from the above two entries was used to map 433 marker loci and 46 linkage groups with a genetic distance of 2,126.3 cM covering approximately 45% of the cotton genome and an average distance between two loci of 4.9 cM. Based on genome-specific chromosomes identified in G. hirsutum tetraploid (A and D), 56.9% of the coverage was located on the A subgenome while 39.7% was assigned to the D subgenome in the genetic map, suggesting that the A subgenome may be more polymorphic and recombinationally active than originally thought. The linkage groups were assigned to 23 of the 26 chromosomes. This is the first genetic map in which the linkage groups A01 and A02/D03 have been assigned to specific chromosomes. In addition the MUSB-derived markers from BAC-end sequences markers allows fine genetic and QTL mapping of important traits and for the first time provides reconciliation of the genetic and physical maps. Limited QTL analyses suggested that loci on chromosomes 2, 3, 12, 15 and 18 may affect variation in fiber quality traits. The original BAC clones containing the newly mapped MUSB that tag the QTLs provide critical DNA regions for the discovery of gene sequences involved in biological processes such as fiber development and pest resistance in cotton. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Recombinational Landscape and Population Genomics of Caenorhabditis elegans

Recombination rate and linkage disequilibrium, the latter a function of population genomic processes, are the critical parameters for mapping by linkage and association, and their patterns in Caenorhabditis elegans are poorly understood. We performed high-density SNP genotyping on a large panel of recombinant inbred advanced intercross lines (RIAILs) of C. elegans to characterize the landscape of recombination and, on a panel of wild strains, to characterize population genomic patterns. We confirmed that C. elegans autosomes exhibit discrete domains of nearly constant recombination rate, and we show, for the first time, that the pattern holds for the X chromosome as well. The terminal domains of each chromosome, spanning about 7% of the genome, exhibit effectively no recombination. The RIAILs exhibit a 5.3-fold expansion of the genetic map. With median marker spacing of 61 kb, they are a powerful resource for mapping quantitative trait loci in C. elegans. Among 125 wild isolates, we identified only 41 distinct haplotypes. The patterns of genotypic similarity suggest that some presumed wild strains are laboratory contaminants. The Hawaiian strain, CB4856, exhibits genetic isolation from the remainder of the global population, whose members exhibit ample evidence of intercrossing and recombining. The population effective recombination rate, estimated from the pattern of linkage disequilibrium, is correlated with the estimated meiotic recombination rate, but its magnitude implies that the effective rate of outcrossing is extremely low, corroborating reports of selection against recombinant genotypes. Despite the low population, effective recombination rate and extensive linkage disequilibrium among chromosomes, which are techniques that account for background levels of genomic similarity, permit association mapping in wild C. elegans strains.     

Population structure and linkage disequilibrium unravelled in tetraploid potato

Association mapping is considered to be an important alternative strategy for the identification of quantitative trait loci (QTL) as compared to traditional QTL mapping. A necessary prerequisite for association analysis to succeed is detailed information regarding hidden population structure and the extent of linkage disequilibrium. A collection of 430 tetraploid potato cultivars, comprising two association panels, has been analysed with 41 AFLP^® and 53 SSR primer combinations yielding 3364 AFLP fragments and 653 microsatellite alleles, respectively. Polymorphism information content values and detected number of alleles for the SSRs studied illustrate that commercial potato germplasm seems to be equally diverse as Latin American landrace material. Genome-wide linkage disequilibrium (LD)—reported for the first time for tetraploid potato—was observed up to approximately 5 cM using r ² higher than 0.1 as a criterion for significant LD. Within-group LD, however, stretched on average twice as far when compared to overall LD. A Bayesian approach, a distance-based hierarchical clustering approach as well as principal coordinate analysis were adopted to enquire into population structure. Groups differing in year of market release and market segment (starch, processing industry and fresh consumption) were repeatedly detected. The observation of LD up to 5 cM is promising because the required marker density is not likely to disable the possibilities for association mapping research in tetraploid potato. Population structure appeared to be weak, but strong enough to demand careful modelling of genetic relationships in subsequent marker-trait association analyses. There seems to be a good chance that linkage-based marker-trait associations can be identified at moderate marker densities.     

Mapping and quantitative trait loci analysis of verticillium wilt resistance genes in cotton

Verticillium wilt is one of the most serious constraints to cotton production in almost all of the cotton-growing countries. In this study, ＂XinLuZaol＂ （XLZl）, a susceptible cultivar Gossypium hirsutum L. and ＂Hai7124＂ （H7124）, a resistant line G. barbadense, and their F2：3 families were used to map and study the disease index induced by verticillium wilt. A total of 430 SSR loci were mapped into 41 linkage groups; the map spanned 3 745.9 cM and the average distance between adjacent loci was 8.71 cM. Four and five quantitative trait loci （QTLs） were detected based on the disease index investigated on July 22 and August 24 in 2004, respectively. These nine QTLs explained 10.63-28.83% of the phenotypic variance, six of them were located on the D sub-genome. Two QTLs located in the same marker intervals may partly explain the significant correlation of the two traits. QTLs explaining large phenotypic variation were identified in this study, which may be quite useful in cotton anti-disease breeding.