Chromosomal assignment of RFLP linkage groups harboring important QTLs on an intraspecific cotton (Gossypium hirsutum L.) Joinmap

Transcriptome profiling, sequence characterization, and SNP-based chromosomal assignment of the EXPANSIN genes in cotton

The knowledge of biological significance associated with DNA markers is very limited in cotton. SNPs are potential functional marker to tag genes of biological importance. Plant expansins are a group of extracellular proteins that directly modify the mechanical properties of cell walls, enable turgor-driven cell extension, and likely affect length and quality of cotton fibers. Here, we report the expression profiles of EXPANSIN transcripts during fiber elongation and the discovery of SNP markers, assess the SNP characteristics, and localize six EXPANSIN A genes to chromosomes. Transcriptome profiling of cotton fiber oligonucleotide microarrays revealed that seven EXPANSIN transcripts were differentially expressed when there was parallel polar elongation during morphogenesis at early stage of fiber development, suggesting that major and minor isoforms perform discrete functions during polar elongation and lateral expansion. Ancestral and homoeologous relationships of the six EXPANSIN A genes were revealed by phylogenetic grouping and comparison to extant A- and D-genome relatives of contemporary AD-genome cottons. The average rate of SNP per nucleotide was 2.35% (one SNP per 43 bp), with 1.74 and 3.99% occurring in coding and noncoding regions, respectively, in the selected genotypes. An unequal evolutionary rate of the EXPANSIN A genes at the subgenome level of tetraploid cotton was recorded. Chromosomal locations for each of the six EXPANSIN A genes were established by gene-specific SNP markers. Results revealed a strategy for discovering SNP markers in a polyploidy species like cotton. These markers could be useful to associate candidate genes with the complex fiber traits in MAS.

A microsatellite-based, gene-rich linkage map reveals genome structure, function and evolution in Gossypium

The mapping of functional genes plays an important role in studies of genome structure, function, and evolution, as well as allowing gene cloning and marker-assisted selection to improve agriculturally important traits. Simple sequence repeats (SSRs) developed from expressed sequence tags (ESTs), EST-SSR (eSSR), can be employed as putative functional marker loci to easily tag corresponding functional genes. In this paper, 2218 eSSRs, 1554 from G. raimondii-derived and 754 from G. hirsutum-derived ESTs, were developed and used to screen polymorphisms to enhance our backbone genetic map in allotetraploid cotton. Of the 1554 G. raimondii-derived eSSRs, 744 eSSRs were able to successfully amplify polymorphisms between our two mapping parents, TM-1 and Hai7124, presenting a polymorphic rate of 47.9%. However, only a 23.9% (159/754) polymorphic rate was produced from G. hirsutum-derived eSSRs. No relationship was observed between the level of polymorphism, motif type, and tissue origin, but the polymorphism appeared to be correlated with repeat type. After integrating these new eSSRs, our enhanced genetic map consists of 1790 loci in 26 linkage groups and covers 3425.8 cM with an average intermarker distance of 1.91 cM. This microsatellite-based, gene-rich linkage map contains 71.96% functional marker loci, of which 87.11% are eSSR loci. There were 132 duplicated loci bridging 13 homeologous At/Dt chromosome pairs. Two reciprocal translocations after polyploidization between A2 and A3, and between A4 and A5, chromosomes were further confirmed. A functional analysis of 975 ESTs producing 1122 eSSR loci tagged in the map revealed that 60% had clear BLASTX hits (<1e(-10)) to the Uniprot database and that 475 were associated mainly with genes belonging to the three major gene ontology categories of biological process, cellular component, and molecular function; many of the ESTs were associated with two or more category functions. The results presented here will provide new insights for future investigations of functional and evolutionary genomics, especially those associated with cotton fiber improvement.

Accumulation of genome-specific transcripts, transcription factors and phytohormonal regulators during early stages of fiber cell development in allotetraploid cotton

Gene expression during the early stages of fiber cell development and in allopolyploid crops is poorly understood. Here we report computational and expression analyses of 32 789 high-quality ESTs derived from Gossypium hirsutum L. Texas Marker-1 (TM-1) immature ovules (GH_TMO). The ESTs were assembled into 8540 unique sequences including 4036 tentative consensus sequences (TCs) and 4504 singletons, representing approximately 15% of the unique sequences in the cotton EST collection. Compared with approximately 178 000 existing ESTs derived from elongating fibers and non-fiber tissues, GH_TMO ESTs showed a significant increase in the percentage of genes encoding putative transcription factors such as MYB and WRKY and genes encoding predicted proteins involved in auxin, brassinosteroid (BR), gibberellic acid (GA), abscisic acid (ABA) and ethylene signaling pathways. Cotton homologs related to MIXTA, MYB5, GL2 and eight genes in the auxin, BR, GA and ethylene pathways were induced during fiber cell initiation but repressed in the naked seed mutant (N1N1) that is impaired in fiber formation. The data agree with the known roles of MYB and WRKY transcription factors in Arabidopsis leaf trichome development and the well-documented phytohormonal effects on fiber cell development in immature cotton ovules cultured in vitro. Moreover, the phytohormonal pathway-related genes were induced prior to the activation of MYB-like genes, suggesting an important role of phytohormones in cell fate determination. Significantly, AA sub-genome ESTs of all functional classifications including cell-cycle control and transcription factor activity were selectively enriched in G. hirsutum L., an allotetraploid derived from polyploidization between AA and DD genome species, a result consistent with the production of long lint fibers in AA genome species. These results suggest general roles for genome-specific, phytohormonal and transcriptional gene regulation during the early stages of fiber cell development in cotton allopolyploids.

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.

Identification and mapping of microsatellite markers linked to a root-knot nematode resistance gene (rkn1) in Acala NemX cotton (Gossypium hirsutum L.)

Host-plant resistance is the most economic and effective strategy for root-knot nematode (RKN) Meloidogyne incognita control in cotton (Gossypium hirsutum L.). Molecular markers linked to resistance are important for incorporating resistance genes into elite cultivars. To screen for microsatellite markers (SSR) closely linked to RKN resistance in G. hirsutum cv. Acala NemX, F1, F2, BC1F1, and F2:7 recombinant inbred lines (RILs) from intraspecific crosses and an F2 from an interspecific cross with G. barbadense cv. Pima S-7 were used. Screening of 284 SSR markers, which cover all the known identified chromosomes and most linkage groups of cotton, was performed by bulked segregant analysis, revealing informative SSRs. The informative SSRs were then mapped on the above populations. One co-dominant SSR marker CIR316 was identified tightly linked to a major resistance gene (designated as rkn1), producing amplified DNA fragments of approximately 221 bp (CIR316a) and 210 bp (CIR316c) in Acala NemX and susceptible Acala SJ-2, respectively. The linkage between CIR316a marker and resistance gene rkn1 in Acala NemX had an estimated distance of 2.1-3.3 cM depending on the population used. Additional markers, including BNL1231 with loose linkage to rkn1 (map distance 25.1-27.4 cM), BNL1066, and CIR003 allowed the rkn1 gene to be mapped to cotton linkage group A03. This is the first report in cotton with a closely linked major gene locus determining nematode resistance, and informative SSRs may be used for marker-assisted selection.

Cleaved AFLP (cAFLP), a modified amplified fragment length polymorphism analysis for cotton

In certain plant species including cotton (Gossypium hirsutum L. or Gossypium barbadense L.), the level of amplified fragment length polymorphism (AFLP) is relatively low, limiting its utilization in the development of genome-wide linkage maps. We propose the use of frequent restriction enzymes in combination with AFLP to cleave the AFLP fragments, called cleaved AFLP analysis (cAFLP). Using four Upland cotton genotypes (G. hirsutum) and three Pima cotton (G. barbadense), we demonstrated that cAFLP generated 67% and 132% more polymorphic markers than AFLP in Upland and Pima cotton, respectively. This resulted in 15.5 and 25.5 polymorphic cAFLP markers per AFLP primer combination, as compared to 9.1 and 11.0 polymorphic AFLP. The cAFLP-based genetic similarity (GS) is generally lower than the AFLP-based GS, even though both marker systems are overall congruent. In some cases, cAFLP can better resolve genetic relationships between genotypes, rendering a higher discriminatory power. Given the high-resolution power of capillary-based DNA sequencing system, we further propose that AFLP and cAFLP amplicons from the same primer combination can be pooled as one sample before electrophoresis. The combination produced an average of 18.5 and 31.0 polymorphic markers per primer pair in Upland and Pima cotton, respectively. Using several restriction enzyme combinations before pre-selective amplification in combination with various frequent 4 bp-cutters or 6 bp-cutters after selective amplification, the pooled AFLP and cAFLP will provide unlimited number of polymorphic markers for genome-wide mapping and fingerprinting.

Genetic mapping of new cotton fiber loci using EST-derived microsatellites in an interspecific recombinant inbred line cotton population

There is an immediate need for a high-density genetic map of cotton anchored with fiber genes to facilitate marker-assisted selection (MAS) for improved fiber traits. With this goal in mind, genetic mapping with a new set of microsatellite markers [comprising both simple (SSR) and complex (CSR) sequence repeat markers] was performed on 183 recombinant inbred lines (RILs) developed from the progeny of the interspecific cross Gossypium hirsutum L. cv. TM1 x Gossypium barbadense L. Pima 3-79. Microsatellite markers were developed using 1557 ESTs-containing SSRs (> or = 10 bp) and 5794 EST-containing CSRs (> or = 12 bp) obtained from approximately 14,000 consensus sequences derived from fiber ESTs generated from the cultivated diploid species Gossypium arboreum L. cv AKA8401. From a total of 1232 EST-derived SSR (MUSS) and CSR (MUCS) primer-pairs, 1019 (83%) successfully amplified PCR products from a survey panel of six Gossypium species; 202 (19.8%) were polymorphic between the G. hirsutum L. and G. barbadense L. parents of the interspecific mapping population. Among these polymorphic markers, only 86 (42.6%) showed significant sequence homology to annotated genes with known function. The chromosomal locations of 36 microsatellites were associated with 14 chromosomes and/or 13 chromosome arms of the cotton genome by hypoaneuploid deficiency analysis, enabling us to assign genetic linkage groups (LG) to specific chromosomes. The resulting genetic map consists of 193 loci, including 121 new fiber loci not previously mapped. These fiber loci were mapped to 19 chromosomes and 11 LG spanning 1277 cM, providing approximately 27% genome coverage. Preliminary quantitative trait loci analysis suggested that chromosomes 2, 3, 15, and 18 may harbor genes for traits related to fiber quality. These new PCR-based microsatellite markers derived from cotton fiber ESTs will facilitate the development of a high-resolution integrated genetic map of cotton for structural and functional study of fiber genes and MAS of genes that enhance fiber quality.

Genetic mapping of new cotton fiber loci using EST-derived microsatellites in an interspecific recombinant inbred line cotton population

There is an immediate need for a high-density genetic map of cotton anchored with fiber genes to facilitate marker-assisted selection (MAS) for improved fiber traits. With this goal in mind, genetic mapping with a new set of microsatellite markers [comprising both simple (SSR) and complex (CSR) sequence repeat markers] was performed on 183 recombinant inbred lines (RILs) developed from the progeny of the interspecific cross Gossypium hirsutum L. cv. TM1 x Gossypium barbadense L. Pima 3-79. Microsatellite markers were developed using 1557 ESTs-containing SSRs (> or = 10 bp) and 5794 EST-containing CSRs (> or = 12 bp) obtained from approximately 14,000 consensus sequences derived from fiber ESTs generated from the cultivated diploid species Gossypium arboreum L. cv AKA8401. From a total of 1232 EST-derived SSR (MUSS) and CSR (MUCS) primer-pairs, 1019 (83%) successfully amplified PCR products from a survey panel of six Gossypium species; 202 (19.8%) were polymorphic between the G. hirsutum L. and G. barbadense L. parents of the interspecific mapping population. Among these polymorphic markers, only 86 (42.6%) showed significant sequence homology to annotated genes with known function. The chromosomal locations of 36 microsatellites were associated with 14 chromosomes and/or 13 chromosome arms of the cotton genome by hypoaneuploid deficiency analysis, enabling us to assign genetic linkage groups (LG) to specific chromosomes. The resulting genetic map consists of 193 loci, including 121 new fiber loci not previously mapped. These fiber loci were mapped to 19 chromosomes and 11 LG spanning 1277 cM, providing approximately 27% genome coverage. Preliminary quantitative trait loci analysis suggested that chromosomes 2, 3, 15, and 18 may harbor genes for traits related to fiber quality. These new PCR-based microsatellite markers derived from cotton fiber ESTs will facilitate the development of a high-resolution integrated genetic map of cotton for structural and functional study of fiber genes and MAS of genes that enhance fiber quality.