SSR markers closely associated with genes for resistance to root-knot nematode on chromosomes 11 and 14 of Upland cotton

Unraveling genomic regions and candidate genes for multiple disease resistance in upland cotton using meta-QTL analysis

Verticillium wilt (VW), Fusarium wilt (FW) and Root-knot nematode (RKN) are the main diseases affecting cotton production. However, many reported quantitative trait loci (QTLs) for cotton resistance have not been used for agricultural practices because of inconsistencies in the cotton genetic background. The integration of existing cotton genetic resources can facilitate the discovery of important genomic regions and candidate genes involved in disease resistance. Here, an improved and comprehensive meta-QTL analysis was conducted on 487 disease resistant QTLs from 31 studies in the last two decades. A consensus linkage map with genetic overall length of 3006.59 cM containing 8650 markers was constructed. A total of 28 Meta-QTLs (MQTLs) were discovered, among which nine MQTLs were identified as related to resistance to multiple diseases. Candidate genes were predicted based on public transcriptome data and enriched in pathways related to disease resistance. This study used a method based on the integration of Meta-QTL, known genes and transcriptomics to reveal major genomic regions and putative candidate genes for resistance to multiple diseases, providing a new basis for marker-assisted selection of high disease resistance in cotton breeding.

Evaluation of Meloidogyne Incognita and Rotylenchulus Reniformis Nematode-resistant Cotton Cultivars with Supplemental Corteva Agriscience Nematicides

Meloidogyne incognita- and Rotylenchulus reniformis-resistant new cotton cultivars have recently become available, giving growers a new option in nematode management. The objectives of this study were: (i) to determine the yield potential of the new cultivars PHY 360 W3FE (M. incognita-resistant) and PHY 332 W3FE (R. reniformis-resistant) in nematode-infested fields and (ii) to evaluate the effects of combining the nematicides Reklemel (fluazaindolizine), Vydate C-LV (oxamyl), and the seed treatment BIOST Nematicide 100 (heat killed Burkholderia rinojenses and its non-living spent fermentation media) with resistant cotton cultivars on nematode population levels and lint yield. Field experiments in 2020 and 2021 indicated M. incognita population levels were 73% lower on PHY 360 W3FE (R) and 80% lower for R. reniformis on the PHY 332 W3FE (R) at 40 days after planting. Nematode eggs per gram of root were further reduced an average of 86% after the addition of Reklemel and Vydate C-LV when averaging both cultivars over the two years. Tests with BIOST Nematicide 100 + Reklemel + Vydate C-LV (0.56 + 2.5 L/ha) in both M. incognita and R. reniformis fields produced higher lint yields. Overall, planting PHY 360 W3FE (R) and PHY 332 W3FE (R) improved yields an average of 364 kg/ha while limiting nematode population increases. The addition of the nematicides further increased yields 152 kg/ha of the nematode-resistant cultivars.

Virulence of Two Isolates of Meloidogyne enterolobii (Guava Root-Knot Nematode) from North Carolina on Cotton Lines Resistant to Southern Root-Knot Nematode (M. incognita) and Reniform Nematode (Rotylenchulus reniformis)

Meloidogyne enterolobii [the guava root-knot nematode (RKN)] is an emerging plant-parasitic nematode that poses a threat to Upland cotton (Gossypium hirsutum) production in the southeastern United States. Like other RKN spp., M. enterolobii has a wide host range and proven ability to overcome resistance sources that have helped protect crops from other Meloidogyne spp., including the southern RKN (Meloidogyne incognita). In this study we evaluated the virulence of two North Carolina M. enterolobii isolates on Upland cotton germplasm lines having resistance quantitative trait loci (QTL) to RKN (M240 RNR, MRk-Rn-1) and/or reniform nematode (Rotylenchulus reniformis) (M713 Ren1, MRk-Rn-1) in comparison to their susceptible recurrent parents (DPL61, SG747). Multiple assays using eggs or J2 as inoculum demonstrated that both isolates reproduced equally well on all germplasm lines, producing reproductive factor (RF) values ≥ 6 on the otherwise nematode-resistant lines. Measurements of seedling growth in control and inoculated containers suggested that existing nematode-resistance QTL may offer a level of tolerance to M. enterolobii infection that should be further explored in greenhouse and field environments. Meloidogyne enterolobii infection of SG747 and MRk-Rn-1 showed nearly identical stages of symptom and nematode development over a time-course of 24 days. These data demonstrate that existing RKN and RN resistance QTL available in elite cotton varieties to producers are most likely insufficient in preventing yield loss due to M. enterolobii and that future research should focus on (i) understanding the M. enterolobii-cotton interaction at the molecular level, and (ii) screening novel germplasm collections to identify resistance loci.

Rotation of Cotton (Gossypium hirsutum) Cultivars and Fallow on Yield and Rotylenchulus reniformis

A three-year rotation of cotton (Gossypium hirsutum) cultivars either resistant (R) or susceptible (S) to Rotylenchulus reniformis and fallow (F) was examined for effect on cotton yield and nematode density. In year 1, 2, and 3, the resistant cultivar (DP 2143NR B3XF) yielded 78, 77, and 113% higher than the susceptible cultivar (DP 2044 B3XF). Fallow in year 1 followed by S in year 2 (F1S2) improved yield in year 2 by 24% compared with S1S2, but not as much as R1S2 (41% yield increase over S1S2). One year of fallow followed by R (F1R2) had lower yield in year 2 (11% reduction) than R1R2. The highest yield after three years of these rotations occurred with R1R2R3, followed by R1S2R3 (17% less yield) and F1F2S3 (35% less yield). Rotylenchulus reniformis density in soil averaged 57, 65, and 70% lower (year 1, 2, 3, respectively) in R1R2R3 compared with S1S2S3. In years 1 and 2, LOG10 transformed nematode density (LREN) was lower in F1, and F1F2, than for all other combinations. In year 3, the lowest LREN were associated with R1R2R3, F1S2F3, and F1F2S3. The highest LREN were associated with F1R2S3, F1S2S3, S1S2S3, R1R2S3, and R1S2S3. The combination of higher yield and lower nematode density will be a strong incentive for producers to use the R. reniformis resistant cultivars continuously.

Effects of Fusarium wilt on Cotton Cultivars with and Without Meloidogyne incognita Resistance in Fields

Cotton (Gossypium hirsutum) cultivar trials were conducted in four fields (6 trials total) with Meloidogyne incognita (Mi)/Fusarium oxysporum f. sp. vasinfectum (Fov) from 2019 to 2021. Cotton cultivars were divided into groups based on company/Mi resistance: S = susceptible to Mi; R-FM, R-DP, and R-PHY = resistance to Mi in FiberMax®, Deltapine®, and Phytogen® cultivars, respectively; ST 4946GLB2 (moderate resistance to Mi and observed field tolerance to Fov); and ST 5600B2XF (resistance to Mi). The S and R-FM groups had the highest transformed Mi densities LOG10(Mi + 1) (LMi = 3.22 and 3.01, respectively), while R-DP and R-PHY had the lowest LMi (2.21 and 1.85, respectively). Plant mortality (%) was higher for R-DP (28.1%) than for all other groups except ST 5600B2XF (24.8%). Mi-susceptible cultivars averaged 23.3% mortality. Relative yield (0-1 scale) was higher for ST 4946GLB2 (0.706) and R-PHY (0.635) than for R-DP (0.530), ST 5600B2XF (0.578), and S (0.491). All groups except R-DP averaged higher relative yield than the susceptible cultivars. ST 4946GLB2 had the lowest mortality (16.5%) and highest relative yield, while R-DP cultivars had the highest mortality and no difference in relative yield from the Mi-susceptible cultivars. The group of R-DP cultivars had excellent Mi resistance but were susceptible to Fov. No cultivars were identified with high resistance to Fov.

A GWAS identified a major QTL for resistance to Fusarium wilt (Fusarium oxysporum f. sp. vasinfectum) race 4 in a MAGIC population of Upland cotton and a meta-analysis of QTLs for Fusarium wilt resistance

A major QTL conferring resistance to Fusarium wilt race 4 in a narrow region of chromosome D02 was identified in a MAGIC population of 550 RILs of Upland cotton. Numerous studies have been conducted to investigate the genetic basis of Fusarium wilt (FW, caused by Fusarium oxysporum f. sp. vasinfectum, FOV) resistance using bi-parental and association mapping populations in cotton. In this study, a multi-parent advanced generation inter-cross (MAGIC) population of 550 recombinant inbred lines (RILs), together with their 11 Upland cotton (Gossypium hirsutum) parents, was used to identify QTLs for FOV race 4 (FOV4) resistance. Among the parents, Acala Ultima, M-240 RNR, and Stoneville 474 were the most resistant, while Deltapine Acala 90, Coker 315, and Stoneville 825 were the most susceptible. Twenty-two MAGIC lines were consistently resistant to FOV4. Through a genome-wide association study (GWAS) based on 473,516 polymorphic SNPs, a major FOV4 resistance QTL within a narrow region on chromosomes D02 was detected, allowing identification of 14 candidate genes. Additionally, a meta-analysis of 133 published FW resistance QTLs showed a D subgenome and individual chromosome bias and no correlation between homeologous chromosome pairs. This study represents the first GWAS study using a largest genetic population and the most comprehensive meta-analysis for FW resistance in cotton. The results illustrated that 550 lines were not enough for high resolution mapping to pinpoint a candidate gene, and experimental errors in phenotyping cotton for FW resistance further compromised the accuracy and precision in QTL localization and identification of candidate genes. This study identified FOV4-resistant parents and MAGIC lines, and the first major QTL for FOV4 resistance in Upland cotton, providing useful information for developing FOV4-resistant cultivars and further genomic studies towards identification of causal genes for FOV4 resistance in cotton.

The effect of two QTLs for resistance to Meloidogyne incognita in cotton on nematode egression from roots

Cotton is widely grown in the southern US and Meloidogyne incognita is its most significant pathogen. The germplasm line M-120 RNR is highly resistant to M. incognita due to two resistance QTLs (quantitative trait loci), qMi-C11 and qMi-C14. Both QTLs reduce total egg production, but the QTLs affect M. incognita development at different life stages. The QTLs do not appear to affect initial penetration of M. incognita but genotypes containing qMi-C11 had fewer nematodes in the roots 8 days after inoculation than near isolines without qMi-C11, which may indicate M. incognita egression from roots. Three greenhouse trials were conducted using cotton isolines to determine whether qMi-C11 and qMi-C14 affect egression of M. incognita juveniles from roots. On each of the five sampling dates (4, 6, 8, 10, and 12 DAI), nematodes that egressed from roots were counted and roots were stained to count nematodes that remained in the roots. The effect of resistance QTLs on M. incognita egression from the roots differed among the trials. Nematode egression was consistently numerically greater, but inconsistently statistically different, from plants with both QTLs than from plants with neither QTL. Plants with only one QTL generally did not differ from plants with both QTLs, and the effects of qMi-C11 and qMi-C14 did not differ in any consistent way. In a separate experiment, plants with neither QTL had more eggs per egg mass than did plants with both QTLs, whereas plants with only one QTL had an intermediate number. Root gall size was measured in two trials and no consistent differences in gall size were observed. We conclude that (1) qMi-C11 and qMi-C14 do not stimulate nematode egression from cotton roots, (2) both qMi-C11 and qMi-C14 reduce M. incognita eggs/egg mass, and (3) neither qMi-C11 nor qMi-C14 affect gall size.

Analysis of Cotton Chromosome 11 and 14 Root-Knot Nematode Resistance Quantitative Trait Loci Effects on Root-Knot Nematode Postinfection Development, Egg Mass Formation, and Fecundity

Cotton (Gossypium hirsutum) resistance to root-knot nematode (RKN) (Meloidogyne incognita) is controlled by quantitative trait loci (QTLs) on chromosomes 11 (CHR11) and 14 (CHR14). The individual contributions of these QTLs to resistance are not completely understood. We developed near isogenic lines susceptible at both loci (null), having CHR11 or CHR14 alone, and having both QTLs (CHR11/CHR14). RKN reproduction, postinfection development, egg mass formation, and adult female fecundity were evaluated. Total RKN reproduction was reduced more in CHR14 versus CHR11 but not as greatly as in CHR11/CHR14. Second-stage juvenile (J2) development to the J3 and J4 (J3+J4) life stages was delayed in CHR11, whereas the J2 transition to J3+J4 in CHR14 followed a similar track as in null plants. Development of J3+J4 nematodes to adult females was inhibited in CHR14 at 21 days after inoculation (DAI). Adult female numbers were decreased in CHR11 and CHR14 at 21 and 28 DAI, with CHR11/CHR14 showing an even greater reduction by 28 DAI. The number of egg masses per gram of root at 21, 28, and 35 DAI formed on CHR11 and CHR14 followed a similar track as numbers of adult females. RKN adult female fecundity (eggs/egg mass) was reduced for CHR11 and CHR14 compared with the null at 21 DAI; however, CHR11 eggs/egg mass was only slightly reduced versus the null by 28 DAI. In contrast, CHR14 eggs/egg mass was like CHR11/CHR14, showing a 4-fold decrease compared with CHR11 and the null.

The relationship between commercial cotton cultivars with varying Meloidogyne incognita resistance genes and yield

Small plot cotton cultivar trials (12 trials) were conducted from 2016 to 2019 in fields infested with Meloidogyne incognita. Entries in these trials included commercial cultivars with partial and high resistance to M. incognita, as well as cultivars with no known resistance. Different resistant groups were created based on different cotton seed companies and their descriptions of the M. incognita resistant cultivars. Groups were none (susceptible); partial resistance found in Stoneville or Fibermax cultivars (PR-FM/ST); partial resistance found in PhytoGen cultivars (PR-PHY); resistance (unknown gene(s)) in Deltapine cultivars (NR-DP); and highly resistant cultivars homozygous for RK1 and RK2 resistant genes in PhytoGen cultivars (HR-PHY). The highest lint yields using a mixed model analysis were found in the PR-FM/ST (1,396 kg lint/ha), HR-PHY (1,327 kg lint/ha), and PR-PHY (1,314 kg lint/ha) groups. Yield for NR-DP (1,234 kg lint/ha) was not different (p > 0.05) than yield for susceptible cultivars (1,243 kg lint/ha). If the older resistant cultivars from Deltapine and PhytoGen (those with only Roundup Ready® herbicide technology) were removed from the analysis, then HR-PHY yields increased by 133 kg of lint/ha to 1,460 kg lint/ha and NR-DP yields remained approximately unchanged (1,227 kg lint/ha). Newer HR-PHY had much improved yield over the first HR-PHY cultivars. Newer HR-PHY averaged 17% higher yield than the susceptible group. LOG10 (M. incognita eggs/500 cm3 soil + 1) were highest for the susceptible cultivars (3.2), followed by PR-FM/ST (2.6), NR-DP (2.4), PR-PHY (2.1), and lowest with HR-PHY (1.4). The newer HR-PHY cultivars (those with ENLIST® herbicide technology) combine excellent yields (17% higher than susceptible cultivars) with high (96%) suppression of M. incognita.

Small plot cotton cultivar trials (12 trials) were conducted from 2016 to 2019 in fields infested with Meloidogyne incognita. Entries in these trials included commercial cultivars with partial and high resistance to M. incognita, as well as cultivars with no known resistance. Different resistant groups were created based on different cotton seed companies and their descriptions of the M. incognita resistant cultivars. Groups were none (susceptible); partial resistance found in Stoneville or Fibermax cultivars (PR-FM/ST); partial resistance found in PhytoGen cultivars (PR-PHY); resistance (unknown gene(s)) in Deltapine cultivars (NR-DP); and highly resistant cultivars homozygous for RK1 and RK2 resistant genes in PhytoGen cultivars (HR-PHY). The highest lint yields using a mixed model analysis were found in the PR-FM/ST (1,396 kg lint/ha), HR-PHY (1,327 kg lint/ha), and PR-PHY (1,314 kg lint/ha) groups. Yield for NR-DP (1,234 kg lint/ha) was not different (p > 0.05) than yield for susceptible cultivars (1,243 kg lint/ha). If the older resistant cultivars from Deltapine and PhytoGen (those with only Roundup Ready® herbicide technology) were removed from the analysis, then HR-PHY yields increased by 133 kg of lint/ha to 1,460 kg lint/ha and NR-DP yields remained approximately unchanged (1,227 kg lint/ha). Newer HR-PHY had much improved yield over the first HR-PHY cultivars. Newer HR-PHY averaged 17% higher yield than the susceptible group. LOG10 (M. incognita eggs/500 cm3 soil + 1) were highest for the susceptible cultivars (3.2), followed by PR-FM/ST (2.6), NR-DP (2.4), PR-PHY (2.1), and lowest with HR-PHY (1.4). The newer HR-PHY cultivars (those with ENLIST® herbicide technology) combine excellent yields (17% higher than susceptible cultivars) with high (96%) suppression of M. incognita.

Sequence Composition of Bacterial Chromosome Clones in a Transgressive Root-Knot Nematode Resistance Chromosome Region in Tetraploid Cotton

Plants evolve innate immunity including resistance genes to defend against pest and pathogen attack. Our previous studies in cotton (Gossypium spp.) revealed that one telomeric segment on chromosome (Chr) 11 in G. hirsutum cv. Acala NemX (rkn1 locus) contributed to transgressive resistance to the plant parasitic nematode Meloidogyne incognita, but the highly homologous segment on homoeologous Chr 21 had no resistance contribution. To better understand the resistance mechanism, a bacterial chromosome (BAC) library of Acala N901 (Acala NemX resistance source) was used to select, sequence, and analyze BAC clones associated with SSR markers in the complex rkn1 resistance region. Sequence alignment with the susceptible G. hirsutum cv. TM-1 genome indicated that 23 BACs mapped to TM-1-Chr11 and 18 BACs mapped to TM-1-Chr 21. Genetic and physical mapping confirmed less BAC sequence (53-84%) mapped with the TM-1 genome in the rkn1 region on Chr 11 than to the homologous region (>89%) on Chr 21. A 3.1-cM genetic distance between the rkn1 flanking markers CIR316 and CIR069 was mapped in a Pima S-7 × Acala NemX RIL population with a physical distance ∼1 Mbp in TM-1. NCBI Blast and Gene annotation indicated that both Chr 11 and Chr 21 harbor resistance gene-rich cluster regions, but more multiple homologous copies of Resistance (R) proteins and of adjacent transposable elements (TE) are present within Chr 11 than within Chr 21. (CC)-NB-LRR type R proteins were found in the rkn1 region close to CIR316, and (TIR)-NB-LRR type R proteins were identified in another resistance rich region 10 cM from CIR 316 (∼3.1 Mbp in the TM-1 genome). The identified unique insertion/deletion in NB-ARC domain, different copies of LRR domain, multiple copies or duplication of R proteins, adjacent protein kinases, or TE in the rkn1 region on Chr 11 might be major factors contributing to complex recombination and transgressive resistance.

Resistance Quantitative Trait Loci qMi-C11 and qMi-C14 in Cotton Have Different Effects on the Development of Meloidogyne incognita, the Southern Root-Knot Nematode

Quantitative trait loci (QTLs) qMi-C11 and qMi-C14 impart a high level of resistance to Meloidogyne incognita in cotton. Breeders had previously backcrossed both QTLs into the susceptible Coker 201 to create the highly resistant M-120 RNR, and we crossed Coker 201 and M-120 RNR to create near-isogenic lines with either qMi-C11 or qMi-C14. Previous work suggests different modes of action for qMi-C11 and qMi-C14. To document individual and combined effects of the QTLs on nematode development and reproduction, Coker 201 (neither QTL), M-120 RNR (both QTLs), CH11 near isoline (qMi-C11), and CH14 near isoline (qMi-C14) were inoculated with M. incognita. At 4, 8, 12, 16, 20, 25, and 30 days after inoculation (DAI), roots were stained to observe nematode developmental stages (second-stage juvenile [J2], swollen second-stage juvenile [SJ2], third-stage juvenile [J3], fourth-stage juvenile [J4], and female), and the number of galls was counted. At 20, 25, 30, and 40 DAI, M. incognita eggs were harvested and counted. At 30 DAI, 80% of the nematodes on Coker 201 were female compared with 50, 40, and 33% females on CH14, CH11, and M-120 RNR, respectively, and greater proportions of nematodes remained in J2 in M-120 RNR (41%), CH11 (58%), and CH14 (27%) than in Coker 201 (9%). More nematodes progressed to J3 or J4 on Coker 201 and CH14 than on CH11 or M-120 RNR. Coker 201 and CH14 had more galls than M-120 RNR. Coker 201 had more eggs than the other genotypes at 30 DAI. Nematode development beyond J2 or SJ2 was significantly reduced by qMi-C11, and development beyond J3 or J4 was significantly reduced by qMi-C14. This study confirms that qMi-C11 and qMi-C14 act at different times and have different effects on the development of M. incognita, and therefore, they have different modes of action.

A novel variant of Gh_D02G0276 is required for root-knot nematode resistance on chromosome 14 (D02) in Upland cotton

MAGIC population sequencing and virus-induced gene silencing identify Gh_D02G0276 as a novel root-knot nematode resistance gene on chromosome 14 in Upland cotton. The southern root-knot nematode [RKN; Meloidogyne incognita (Kofoid & White)] remains the primary yield-limiting biotic stress to Upland cotton (Gossypium hirsutum L.) throughout the southeastern USA. While useful genetic markers have been developed for two major RKN resistance loci on chromosomes 11 (A11) and 14 (D02), these markers are not completely effective because the causative genes have not been identified. Here, we sequenced 550 recombinant inbred lines (RILs) from a multi-parent advanced generation intercross (MAGIC) population to identify five RILs that had informative recombinations near the D02-RKN resistance locus. The RKN resistance phenotypes of these five RILs narrowed the D02-RKN locus to a 30-kb region with four candidate genes. We conducted virus-induced gene silencing (VIGS) on each of these genes and found that Gh_D02G0276 was required for suppression of RKN egg production conferred by the Chr. D02 resistance gene. The resistant lines all possessed an allele of Gh_D02G0276 that showed non-synonymous mutations and was prematurely truncated. Furthermore, a Gh_D02G0276-specific marker for the resistance allele variant was able to identify RKN-resistant germplasm from a collection of 367 cotton accessions. The Gh_D02G0276 peptide shares similarity with domesticated hAT-like transposases with additional novel N- and C-terminal domains that resemble the target of known RKN effector molecules and a prokaryotic motif, respectively. The truncation in the resistance allele results in a loss of a plant nuclear gene-specific C-terminal motif, potentially rendering this domain antigenic due to its high homology with bacterial proteins. The conclusive identification of this RKN resistance gene opens new avenues for understanding plant resistance mechanisms to RKN as well as opportunities to develop more efficient marker-assisted selection in cotton breeding programs.

Quantitative Trait Loci Mapping of Multiple Independent Loci for Resistance to Fusarium oxysporum f. sp. vasinfectum Races 1 and 4 in an Interspecific Cotton Population

Fusarium wilt, caused by the soilborne fungal pathogen Fusarium oxysporum f. sp. vasinfectum, is a vascular disease of cotton (Gossypium spp.). F. oxysporum f. sp. vasinfectum race 1 (FOV1) causes major plant injury and yield loss in G. hirsutum cultivars with coinfection with root-knot nematode (Meloidogyne incognita), while F. oxysporum f. sp. vasinfectum race 4 (FOV4) causes plant damage without nematode coinfection in G. hirsutum and in G. barbadense cultivars. Quantitative trait loci (QTL) analysis of the interspecific cross G. barbadense Pima S-7 × G. hirsutum Acala NemX revealed separate multiple loci determining resistance to FOV1 and FOV4, confirming that race specificity occurs in F. oxysporum f. sp. vasinfectum. Based on the area under the disease progress stairs, six major QTLs on chromosomes (Chrs) 1, 2, 12, 15 (2), and 21 contributing 7 to 15% to FOV1 resistance and two major QTLs on Chrs 14 and 17 contributing 12 to 33% to FOV4 resistance were identified. Minor-effect QTLs contributing to resistance to both FOV1 and FOV4 were also identified. These results define and establish a pathosystem of race-specific resistance under polygenic control. This research also validates the importance of previously reported markers and chromosome regions and adds new information for the location of F. oxysporum f. sp. vasinfectum resistance genes. Some F₈ recombinant inbred lines have resistance to both FOV1 and FOV4 and also to root-knot nematode, providing multiple resistance sources for breeding.

MicroRNA-target gene responses to root knot nematode (Meloidogyne incognita) infection in cotton (Gossypium hirsutum L.)

MicroRNAs (miRNAs) are a large class of small regulatory RNA molecules, however no study has been performed to elucidate the role of miRNAs in cotton (Gossypium hirsutum) response to the root knot nematode (RKN, Meloidogyne incognita) infection. We selected 28 miRNAs and 8 miRNA target genes to investigate the miRNA-target gene response to M. incognita infection. Our results show that RKN infection significantly affected the expression of several miRNAs and their targeted genes. After 10 days of RKN infection, expression fold changes on miRNA expressions ranged from down-regulated by 33% to upregulated by 406%; meanwhile the expression levels of miRNA target genes were 45.8% to 231%. Three miRNA-target pairs, miR159-MYB, miR319-TCP4 and miR167-ARF8, showed inverse expression patterns between gene targets and their corresponded miRNAs, suggesting miRNA-mediated gene regulation in cotton roots in response to RKN infection.

Quantitative trait loci mapping of Meloidogyne incognita and M. hapla resistance in a recombinant inbred line population of soybean

A recombinant inbred line population of soybean (Glycine max) was utilised to identify the quantitative trait loci (QTLs) determining the response to infection by two root-knot nematode species, Meloidogyne incognita and M. hapla, in glasshouse assays. QTL analysis detected seven major and four minor QTLs on seven soybean chromosomes ((Chrs) 1, 7, 8, 10, 14, 18, 20) explaining 6-41% phenotypic variance (PVE) for M. incognita root response and nematode reproduction. Three of the major QTLs, on Chrs 7, 10 and 18, were confirmed in previous reports and two major QTLs on Chrs 14 and 20 were detected for the first time. The QTL analysis with M. hapla provides the first report of a major QTL region mapped on Chr 7, explaining 70-82% PVE in M. hapla root response and nematode reproduction. These novel identified QTLs with flanking markers will be helpful in marker-assisted breeding for nematode resistance in soybean.

QTL Analysis of Transgressive Nematode Resistance in Tetraploid Cotton Reveals Complex Interactions in Chromosome 11 Regions

Transgressive segregation in cotton (Gossypium spp.) provides an important approach to enhance resistance to the major pest root-knot nematode (RKN) Meloidogyne incognita. Our previous studies reported transgressive RKN resistance in an intraspecific Gossypium hirsutum resistant NemX × susceptible SJ-2 recombinant inbred line (RIL) population and early generations of interspecific cross Gossypium barbadense (susceptible Pima S-7) × G. hirsutum (NemX). However, the underlying functional mechanisms for this phenomenon are not known. In this study, the region of RKN resistance gene rkn1 on chromosome (Chr) 11 and its homoeologous Chr 21 was fine mapped with G. raimondii D₅ genome reference sequence. Transgressive resistance was found in the later generation of a new RIL population F_2:7 (Pima S-7 × NemX) and one interspecific F₂ (susceptible Pima S-7 × susceptible SJ-2). QTL analysis revealed similar contributions to root-galling and egg-production resistance phenotypes associated with SSR marker CIR316 linked to resistance gene rkn1 in NemX on Chr 11 in all seven populations analyzed. In testcross NemX × F₁ (Pima S-7 × SJ-2) marker allele CIR069-271 from Pima S-7 linked to CIR316 contributed 63% of resistance to galling phenotype in the presence of rkn1. Similarly, in RIL population F_2:8 (NemX × SJ-2), SJ-2 markers closely linked to CIR316 contributed up to 82% of resistance to root-galling. These results were confirmed in BC₁F₁ SJ-2 × F₁ (NemX × SJ-2), F₂ (NemX × SJ-2), and F₂ (Pima S-7 × SJ-2) populations in which up to 44, 36, and 15% contribution in resistance to galling was found, respectively. Transgressive segregation for resistance was universal in all intra- and inter-specific populations, although stronger transgressive resistance occurred in later than in early generations in the intraspecific cross compared with the interspecific cross. Transgressive effects on progeny from susceptible parents are possibly provided in the rkn1 resistance region of chromosome 11 by tandemly arrayed allele (TAA) or gene (TAG) interactions contributing to transgressive resistance. Complex TAA and TAG recombination and interactions in the rkn1 resistance region provide three genes and a model to study disease and transgressive resistance in polyploid plants, and novel genotypes for plant breeding.

QTL mapping for flowering-time and photoperiod insensitivity of cotton Gossypium darwinii Watt

Most wild and semi-wild species of the genus Gossypium are exhibit photoperiod-sensitive flowering. The wild germplasm cotton is a valuable source of genes for genetic improvement of modern cotton cultivars. A bi-parental cotton population segregating for photoperiodic flowering was developed by crossing a photoperiod insensitive irradiation mutant line with its pre-mutagenesis photoperiodic wild-type G. darwinii Watt genotype. Individuals from the F2 and F3 generations were grown with their parental lines and F1 hybrid progeny in the long day and short night summer condition (natural day-length) of Uzbekistan to evaluate photoperiod sensitivity, i.e., flowering-time during the seasons 2008-2009. Through genotyping the individuals of this bi-parental population segregating for flowering-time, linkage maps were constructed using 212 simple-sequence repeat (SSR) and three cleaved amplified polymorphic sequence (CAPS) markers. Six QTLs directly associated with flowering-time and photoperiodic flowering were discovered in the F2 population, whereas eight QTLs were identified in the F3 population. Two QTLs controlling photoperiodic flowering and duration of flowering were common in both populations. In silico annotations of the flanking DNA sequences of mapped SSRs from sequenced cotton (G. hirsutum L.) genome database has identified several potential 'candidate' genes that are known to be associated with regulation of flowering characteristics of plants. The outcome of this research will expand our understanding of the genetic and molecular mechanisms of photoperiodic flowering. Identified markers should be useful for marker-assisted selection in cotton breeding to improve early flowering characteristics.

Coupling of MIC-3 overexpression with the chromosomes 11 and 14 root-knot nematode (RKN) (Meloidogyne incognita) resistance QTLs provides insights into the regulation of the RKN resistance response in Upland cotton (Gossypium hirsutum)

Genetic analysis of MIC-3 transgene with RKN resistance QTLs provides insight into the resistance regulatory mechanism and provides a framework for testing additional hypotheses. Resistance to root-knot nematode (RKN) (Meloidogyne incognita) in Upland cotton (Gossypium hirsutum) is mediated by two major quantitative trait loci (QTL) located on chromosomes 11 and 14. The MIC-3 (Meloidogyne Induced Cotton3) protein accumulates specifically within the immature galls of RKN-resistant plants that possess these QTLs. Recently, we showed that MIC-3 overexpression in an RKN-susceptible cotton genotype suppressed RKN egg production but not RKN-induced root galling. In this study, the MIC-3 overexpression construct T-DNA in the single-copy transgenic line '14-7-1' was converted into a codominant molecular marker that allowed the marker assisted selection of F2:3 cotton lines, derived from a cross between 14-7-1 and M-240 RNR, having all possible combinations of the chromosomes 11 and 14 QTLs with and without the MIC-3 overexpression construct. Root-knot nematode reproduction (eggs g(-1) root) and severity of RKN-induced root galling were assessed in these lines. We discovered that the addition of MIC-3 overexpression suppressed RKN reproduction in lines lacking both resistance QTLs and in lines having only the chromosome 14 QTL, suggesting an additive effect of the MIC-3 construct with this QTL. In contrast, MIC-3 overexpression did not improve resistance in lines having the single chromosome 11 QTL or in lines having both resistance QTLs, suggesting an epistatic interaction between the chromosome 11 QTL and the MIC-3 construct. Overexpression of MIC-3 did not affect the severity of RKN-induced root galling regardless of QTL genotype. These data provide new insights into the relative order of action of the chromosomes 11 and 14 QTLs and their potential roles in regulating MIC-3 expression as part of the RKN resistance response.

Fine mapping and identification of candidate genes for a QTL affecting Meloidogyne incognita reproduction in Upland cotton

Background

The southern root-knot nematode (Meloidogyne incognita; RKN) is one of the most important economic pests of Upland cotton (Gossypium hirsutum L.). Host plant resistance, the ability of a plant to suppress nematode reproduction, is the most economical, practical, and environmentally sound method to provide protection against this subterranean pest. The resistant line Auburn 623RNR and a number of elite breeding lines derived from it remain the most important source of root-knot nematode (RKN) resistance. Prior genetic analysis has identified two epistatically interacting RKN resistance QTLs, qMi-C11 and qMi-C14, affecting gall formation and RKN reproduction, respectively.

Results

We developed a genetic population segregating only for the qMi-C14 locus and evaluated the genetic effects of this QTL on RKN resistance in the absence of the qMi-C11 locus. The qMi-C14 locus had a LOD score of 12 and accounted for 24.5 % of total phenotypic variation for egg production. In addition to not being significantly associated with gall formation, this locus had a lower main effect on RKN reproduction than found in our previous study, which lends further support to evidence of epistasis with qMi-C11 in imparting RKN resistance in the Auburn 623RNR source. The locus qMi-C14 was fine-mapped with the addition of 16 newly developed markers. By using the reference genome sequence of G. raimondii, we identified 20 candidate genes encoding disease resistance protein homologs in the newly defined 2.3 Mb region flanked by two SSR markers. Resequencing of an RKN resistant and susceptible G. hirsutum germplasm revealed non-synonymous mutations in only four of the coding regions of candidate genes, and these four genes are consequently of high interest.

Conclusions

Our mapping results validated the effects of the qMi-C14 resistance locus, delimiting the QTL to a smaller region, and identified tightly linked SSR markers to improve the efficiency of marker-assisted selection. The candidate genes identified warrant functional studies that will help in identifying and characterizing the actual qMi-C14 defense gene(s) against root-knot nematodes.

Electronic supplementary material

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

Reproduction of Meloidogyne incognita Race 3 on Flue-cured Tobacco Homozygous for Rk1 and/or Rk2 Resistance Genes

Most commercial tobacco cultivars possess the Rk1 resistance gene to races 1 and 3 of Meloidogyne incognita and race 1 of Meloidogyne arenaria, which has caused a shift in population prevalence in Virginia tobacco fields toward other species and races. A number of cultivars now also possess the Rk2 gene for root-knot resistance. Experiments were conducted in 2013 to 2014 to examine whether possessing both Rk1 and Rk2 increases resistance to a variant of M. incognita race 3 compared to either gene alone. Greenhouse trials were arranged in a completely randomized design with Coker 371-Gold (C371G; susceptible), NC 95 and SC 72 (Rk1Rk1), T-15-1-1 (Rk2Rk2), and STNCB-2-28 and NOD 8 (Rk1Rk1 and Rk2Rk2). Each plant was inoculated with 5,000 root-knot nematode eggs; data were collected 60 d postinoculation. Percent galling and numbers of egg masses and eggs were counted, the latter being used to calculate the reproductive index on each host. Despite variability, entries with both Rk1 and Rk2 conferred greater resistance to a variant of M. incognita race 3 than plants with Rk1 or Rk2 alone. Entries with Rk1 alone were successful in reducing root galling and nematode reproduction compared to the susceptible control. Entry T-15-1-1 did not reduce galling compared to the susceptible control but often suppressed reproduction.

Sequence composition of BAC clones and SSR markers mapped to Upland cotton chromosomes 11 and 21 targeting resistance to soil-borne pathogens

Genetic and physical framework mapping in cotton (Gossypium spp.) were used to discover putative gene sequences involved in resistance to common soil-borne pathogens. Chromosome (Chr) 11 and its homoeologous Chr 21 of Upland cotton (G. hirsutum) are foci for discovery of resistance (R) or pathogen-induced R (PR) genes underlying QTLs involved in response to root-knot nematode (Meloidogyne incognita), reniform nematode (Rotylenchulus reniformis), Fusarium wilt (Fusarium oxysporum f.sp. vasinfectum), Verticillium wilt (Verticillium dahliae), and black root rot (Thielaviopsis basicola). Simple sequence repeat (SSR) markers and bacterial artificial chromosome (BAC) clones from a BAC library developed from the Upland cotton Acala Maxxa were mapped on Chr 11 and Chr 21. DNA sequence through Gene Ontology (GO) of 99 of 256 Chr 11 and 109 of 239 Chr 21 previously mapped SSRs revealed response elements to internal and external stimulus, stress, signaling process, and cell death. The reconciliation between genetic and physical mapping of gene annotations from new DNA sequences of 20 BAC clones revealed 467 (Chr 11) and 285 (Chr 21) G. hirsutum putative coding sequences, plus 146 (Chr 11) and 98 (Chr 21) predicted genes. GO functional profiling of Unigenes uncovered genes involved in different metabolic functions and stress response elements (SRE). Our results revealed that Chrs 11 and 21 harbor resistance gene rich genomic regions. Sequence comparisons with the ancestral diploid D5 (G. raimondii), A2 (G. arboreum) and domesticated tetraploid TM-1 AD1 (G. hirsutum) genomes revealed abundance of transposable elements and confirmed the richness of resistance gene motifs in these chromosomes. The sequence information of SSR markers and BAC clones and the genetic mapping of BAC clones provide enhanced genetic and physical frameworks of resistance gene-rich regions of the cotton genome, thereby aiding discovery of R and PR genes and breeding for resistance to cotton diseases.

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.

Cotton QTLdb: a cotton QTL database for QTL analysis, visualization, and comparison between Gossypium hirsutum and G. hirsutum × G. barbadense populations

KEY MESSAGE

A specialized database currently containing more than 2200 QTL is established, which allows graphic presentation, visualization and submission of QTL. In cotton quantitative trait loci (QTL), studies are focused on intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations. These two populations are commercially important for the textile industry and are evaluated for fiber quality, yield, seed quality, resistance, physiological, and morphological trait QTL. With meta-analysis data based on the vast amount of QTL studies in cotton it will be beneficial to organize the data into a functional database for the cotton community. Here we provide a tool for cotton researchers to visualize previously identified QTL and submit their own QTL to the Cotton QTLdb database. The database provides the user with the option of selecting various QTL trait types from either the G. hirsutum or G. hirsutum × G. barbadense populations. Based on the user's QTL trait selection, graphical representations of chromosomes of the population selected are displayed in publication ready images. The database also provides users with trait information on QTL, LOD scores, and explained phenotypic variances for all QTL selected. The CottonQTLdb database provides cotton geneticist and breeders with statistical data on cotton QTL previously identified and provides a visualization tool to view QTL positions on chromosomes. Currently the database (Release 1) contains 2274 QTLs, and succeeding QTL studies will be updated regularly by the curators and members of the cotton community that contribute their data to keep the database current. The database is accessible from http://www.cottonqtldb.org.

Cotton QTLdb: a cotton QTL database for QTL analysis, visualization, and comparison between Gossypium hirsutum and G. hirsutum × G. barbadense populations

KEY MESSAGE

A specialized database currently containing more than 2200 QTL is established, which allows graphic presentation, visualization and submission of QTL. In cotton quantitative trait loci (QTL), studies are focused on intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations. These two populations are commercially important for the textile industry and are evaluated for fiber quality, yield, seed quality, resistance, physiological, and morphological trait QTL. With meta-analysis data based on the vast amount of QTL studies in cotton it will be beneficial to organize the data into a functional database for the cotton community. Here we provide a tool for cotton researchers to visualize previously identified QTL and submit their own QTL to the Cotton QTLdb database. The database provides the user with the option of selecting various QTL trait types from either the G. hirsutum or G. hirsutum × G. barbadense populations. Based on the user's QTL trait selection, graphical representations of chromosomes of the population selected are displayed in publication ready images. The database also provides users with trait information on QTL, LOD scores, and explained phenotypic variances for all QTL selected. The CottonQTLdb database provides cotton geneticist and breeders with statistical data on cotton QTL previously identified and provides a visualization tool to view QTL positions on chromosomes. Currently the database (Release 1) contains 2274 QTLs, and succeeding QTL studies will be updated regularly by the curators and members of the cotton community that contribute their data to keep the database current. The database is accessible from http://www.cottonqtldb.org.

Overexpression of MIC-3 indicates a direct role for the MIC gene family in mediating Upland cotton (Gossypium hirsutum) resistance to root-knot nematode (Meloidogyne incognita)

Transgene-based analysis of the MIC-3 gene provides the first report of a cotton gene having a direct role in mediating cotton resistance to root-knot nematode. Major quantitative trait loci have been mapped to Upland cotton (Gossypium hirsutum L.) chromosomes 11 and 14 that govern the highly resistant phenotype in response to infection by root-knot nematode (RKN; Meloidogyne incognita); however, nearly nothing is known regarding the underlying molecular determinants of this RKN-resistant phenotype. Multiple lines of circumstantial evidence have strongly suggested that the MIC (Meloidogyne Induced Cotton) gene family plays an integral role in mediating cotton resistance to RKN. In this report, we demonstrate that overexpression of MIC-3 in the RKN-susceptible genetic background Coker 312 reduces RKN egg production by ca. 60-75 % compared to non-transgenic controls and transgene-null sibling lines. MIC-3 transcript and protein overexpression were confirmed in root tissues of multiple independent transgenic lines with each line showing a similar level of increased resistance to RKN. In contrast to RKN fecundity, transgenic lines showed RKN-induced root galling similar to the susceptible controls. In addition, we determined that this effect of MIC-3 overexpression was specific to RKN as no effect was observed on reniform nematode (Rotylenchulus reniformis) reproduction. Transgenic lines did not show obvious alterations in growth, morphology, flowering, or fiber quality traits. Gene expression analyses showed that MIC-3 transcript levels in uninfected transgenic roots exceeded levels observed in RKN-infected roots of naturally resistant plants and that overexpression did not alter the regulation of native MIC genes in the genome. These results are the first report describing a direct role for a specific gene family in mediating cotton resistance to a plant-parasitic nematode.

A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations

KEY MESSAGE

Based on 1075 and 1059 QTL from intraspecific Upland and interspecific Upland × Pima populations, respectively, the identification of QTL clusters and hotspots provides a useful resource for cotton breeding. Mapping of quantitative trait loci (QTL) is a pre-requisite of marker-assisted selection for crop yield and quality. Recent meta-analysis of QTL in tetraploid cotton (Gossypium spp.) has identified regions of the genome with high concentrations of QTL for various traits called clusters and specific trait QTL called hotspots or meta-QTL (mQTL). However, the meta-analysis included all population types of Gossypium mixing both intraspecific G. hirsutum and interspecific G. hirsutum × G. barbadense populations. This study used 1,075 QTL from 58 publications on intraspecific G. hirsutum and 1,059 QTL from 30 publications on G. hirsutum × G. barbadense populations to perform a comprehensive comparative analysis of QTL clusters and hotspots between the two populations for yield, fiber and seed quality, and biotic and abiotic stress tolerance. QTL hotspots were further analyzed for mQTL within the hotspots using Biomercator V3 software. The ratio of QTL between the two population types was proportional yet differences in hotspot type and placement were observed between the two population types. However, on some chromosomes QTL clusters and hotspots were similar between the two populations. This shows that there are some universal QTL regions in the cultivated tetraploid cotton which remain consistent and some regions which differ between population types. This study for the first time elucidates the similarities and differences in QTL clusters and hotspots between intraspecific and interspecific populations, providing an important resource to cotton breeding programs in marker-assisted selection .

Re-evaluation of the inheritance for root-knot nematode resistance in the Upland cotton germplasm line M-120 RNR revealed two epistatic QTLs conferring resistance

We report a second major QTL for root-knot nematode resistance in the highly resistant Upland cotton line M-120RNR and show epistasis between two resistant QTLs with different mechanisms conferring resistance. In an earlier study, we identified a major QTL on Chromosome 11 associated with resistance to root-knot nematode in the M-120 RNR Upland cotton line (Gossypium hirsutum L.) of the Auburn 623 RNR source. Herein, we re-evaluated the genetics of the resistance to root-knot nematode in the M-120 RNR × Pima S-6 population by linkage mapping using recently published SSR markers. The QTL analysis detected two regions significantly associated with the resistance phenotype. In addition to the QTL previously identified on Chromosome 11 (qMi-C11), a major QTL was identified on Chromosome 14 (qMi-C14). The resistance locus on qMi-C11 originated from the Clevewilt parent, while the qMi-C14 locus originated from the other resistant parent, Mexico Wild Jack Jones. The qMi-C14 locus had logarithms of odds score of 17 and accounted for 45 % of the total phenotype variation in egg production. It was also associated with galling index, but the percent variation explained was only 6 %, suggesting that the qMi-C11 locus had a much stronger effect on root gall suppression than egg production, while the qMi-C14 locus had a stronger effect on egg production than galling. The results also suggest that the transgressive segregation observed in the development of Auburn 623 RNR was due to the pyramiding of at least two main effect QTLs as well as an additive-by-additive epistatic effects between the two resistant loci. The SSRs markers tightly linked to the qMi-C11 and qMi-C14 loci will greatly facilitate the improvement of RKN resistance in cotton via marker-assisted breeding.

A comprehensive meta QTL analysis for fiber quality, yield, yield related and morphological traits, drought tolerance, and disease resistance in tetraploid cotton

BACKGROUND

The study of quantitative trait loci (QTL) in cotton (Gossypium spp.) is focused on traits of agricultural significance. Previous studies have identified a plethora of QTL attributed to fiber quality, disease and pest resistance, branch number, seed quality and yield and yield related traits, drought tolerance, and morphological traits. However, results among these studies differed due to the use of different genetic populations, markers and marker densities, and testing environments. Since two previous meta-QTL analyses were performed on fiber traits, a number of papers on QTL mapping of fiber quality, yield traits, morphological traits, and disease resistance have been published. To obtain a better insight into the genome-wide distribution of QTL and to identify consistent QTL for marker assisted breeding in cotton, an updated comparative QTL analysis is needed.

RESULTS

In this study, a total of 1,223 QTL from 42 different QTL studies in Gossypium were surveyed and mapped using Biomercator V3 based on the Gossypium consensus map from the Cotton Marker Database. A meta-analysis was first performed using manual inference and confirmed by Biomercator V3 to identify possible QTL clusters and hotspots. QTL clusters are composed of QTL of various traits which are concentrated in a specific region on a chromosome, whereas hotspots are composed of only one trait type. QTL were not evenly distributed along the cotton genome and were concentrated in specific regions on each chromosome. QTL hotspots for fiber quality traits were found in the same regions as the clusters, indicating that clusters may also form hotspots.

CONCLUSIONS

Putative QTL clusters were identified via meta-analysis and will be useful for breeding programs and future studies involving Gossypium QTL. The presence of QTL clusters and hotspots indicates consensus regions across cultivated tetraploid Gossypium species, environments, and populations which contain large numbers of QTL, and in some cases multiple QTL associated with the same trait termed a hotspot. This study combines two previous meta-analysis studies and adds all other currently available QTL studies, making it the most comprehensive meta-analysis study in cotton to date.

Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin-streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N(50) of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F(2) population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome ('A' or 'D') each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.

QTL analysis for transgressive resistance to root-knot nematode in interspecific cotton (Gossypium spp.) progeny derived from susceptible parents

The southern root-knot nematode (RKN, Meloidogyne incognita) is a major soil-inhabiting plant parasite that causes significant yield losses in cotton (Gossypium spp.). Progeny from crosses between cotton genotypes susceptible to RKN produced segregants in subsequent populations which were highly resistant to this parasite. A recombinant inbred line (RIL) population of 138 lines developed from a cross between Upland cotton TM-1 (G. hirsutum L.) and Pima 3–79 (G. barbadense L.), both susceptible to RKN, was used to identify quantitative trait loci (QTLs) determining responses to RKN in greenhouse infection assays with simple sequence repeat (SSR) markers. Compared to both parents, 53.6% and 52.1% of RILs showed less (P<0.05) root-galling index (GI) and had lower (P<0.05) nematode egg production (eggs per gram root, EGR). Highly resistant lines (transgressive segregants) were identified in this RIL population for GI and/or EGR in two greenhouse experiments. QTLs were identified using the single-marker analysis nonparametric mapping Kruskal-Wallis test. Four major QTLs located on chromosomes 3, 4, 11, and 17 were identified to account for 8.0 to 12.3% of the phenotypic variance (R²) in root-galling. Two major QTLs accounting for 9.7% and 10.6% of EGR variance were identified on chromosomes 14 and 23 (P<0.005), respectively. In addition, 19 putative QTLs (P<0.05) accounted for 4.5–7.7% of phenotypic variance (R²) in GI, and 15 QTLs accounted for 4.2–7.3% of phenotypic variance in EGR. In lines with alleles positive for resistance contributed by both parents in combinations of two to four QTLs, dramatic reductions of >50% in both GI and EGR were observed. The transgressive segregants with epistatic effects derived from susceptible parents indicate that high levels of nematode resistance in cotton may be attained by pyramiding positive alleles using a QTL mapping approach.

Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin-streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N(50) of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F(2) population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome ('A' or 'D') each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.