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

Marker-assisted breeding as next-generation strategy for genetic improvement of productivity and quality: can it be realized in cotton?

The dawdling development in genetic improvement of cotton with conventional breeding program is chiefly due to lack of complete knowledge on and precise manipulation of fiber productivity and quality. Naturally available cotton continues to be a resource for the upcoming breeding program, and contemporary technologies to exploit the available natural variation are outlined in this paper for further improvement of fiber. Particularly emphasis is given to application, obstacles, and perspectives of marker-assisted breeding since it appears to be more promising in manipulating novel genes that are available in the cotton germplasm. Deployment of system quantitative genetics in marker-assisted breeding program would be essential to realize its role in cotton. At the same time, role of genetic engineering and in vitro mutagenesis cannot be ruled out in genetic improvement of cotton.

Agronomic characterization, genetic diversity and association analysis of cotton cultivars using simple sequence repeat molecular markers

Cotton is the most important textile plant in the world and is one of the most important crops for the production of oilseed. Because of its worldwide economic importance, new cultivars are constantly being released in the world and consequently in the Greek market, as Greece is the largest producer in Europe. We used simple sequence repeat (SSR) markers for the identification and the phylogenetic analysis of the most widely cultivated cotton cultivars in Greece. Initially, we used 12 pairs of SSR molecular markers for the analysis of 29 cultivars of Gossypium hirsutum and an interspecific hybrid (G. hirsutum x G. barbadense). Of the 12 pairs of SSR primers, 11 amplified polymorphic products, while one pair did not amplify any product. Globally, 17 polymorphic marker loci were identified. Two to four different alleles were amplified at each genomic locus, with a mean of 2.53 alleles per locus. Among the 30 genotypes that we analyzed, the polymorphism information content ranged from 0 to 0.548, with a mean of 0.293. Three main groups were formed among the 30 genotypes when a phylogenetic analysis was performed using UPGMA. Computational analysis of each molecular marker separately showed an association of SSR markers with agronomic traits such as fiber quality. To our knowledge, this is the first in-depth molecular analysis of cotton cultivars grown in Greece using SSR markers. An analysis of association of SSR markers with fiber quality traits of 29 cotton cultivars is reported for the first time.

Genome structure of cotton revealed by a genome-wide SSR genetic map constructed from a BC1 population between gossypium hirsutum and G. barbadense

Background

Cotton, with a large genome, is an important crop throughout the world. A high-density genetic linkage map is the prerequisite for cotton genetics and breeding. A genetic map based on simple polymerase chain reaction markers will be efficient for marker-assisted breeding in cotton, and markers from transcribed sequences have more chance to target genes related to traits. To construct a genome-wide, functional marker-based genetic linkage map in cotton, we isolated and mapped expressed sequence tag-simple sequence repeats (EST-SSRs) from cotton ESTs derived from the A₁, D₅, (AD)₁, and (AD)₂ genome.

Results

A total of 3177 new EST-SSRs developed in our laboratory and other newly released SSRs were used to enrich our interspecific BC₁ genetic linkage map. A total of 547 loci and 911 loci were obtained from our EST-SSRs and the newly released SSRs, respectively. The 1458 loci together with our previously published data were used to construct an updated genetic linkage map. The final map included 2316 loci on the 26 cotton chromosomes, 4418.9 cM in total length and 1.91 cM in average distance between adjacent markers. To our knowledge, this map is one of the three most dense linkage maps in cotton. Twenty-one segregation distortion regions (SDRs) were found in this map; three segregation distorted chromosomes, Chr02, Chr16, and Chr18, were identified with 99.9% of distorted markers segregating toward the heterozygous allele. Functional analysis of SSR sequences showed that 1633 loci of this map (70.6%) were transcribed loci and 1332 loci (57.5%) were translated loci.

Conclusions

This map lays groundwork for further genetic analyses of important quantitative traits, marker-assisted selection, and genome organization architecture in cotton as well as for comparative genomics between cotton and other species. The segregation distorted chromosomes can be a guide to identify segregation distortion loci in cotton. The annotation of SSR sequences identified frequent and rare gene ontology items on each chromosome, which is helpful to discover functions of cotton chromosomes.

Polyploidization altered gene functions in cotton (Gossypium spp.)

Cotton (Gossypium spp.) is an important crop plant that is widely grown to produce both natural textile fibers and cottonseed oil. Cotton fibers, the economically more important product of the cotton plant, are seed trichomes derived from individual cells of the epidermal layer of the seed coat. It has been known for a long time that large numbers of genes determine the development of cotton fiber, and more recently it has been determined that these genes are distributed across At and Dt subgenomes of tetraploid AD cottons. In the present study, the organization and evolution of the fiber development genes were investigated through the construction of an integrated genetic and physical map of fiber development genes whose functions have been verified and confirmed. A total of 535 cotton fiber development genes, including 103 fiber transcription factors, 259 fiber development genes, and 173 SSR-contained fiber ESTs, were analyzed at the subgenome level. A total of 499 fiber related contigs were selected and assembled. Together these contigs covered about 151 Mb in physical length, or about 6.7% of the tetraploid cotton genome. Among the 499 contigs, 397 were anchored onto individual chromosomes. Results from our studies on the distribution patterns of the fiber development genes and transcription factors between the At and Dt subgenomes showed that more transcription factors were from Dt subgenome than At, whereas more fiber development genes were from At subgenome than Dt. Combining our mapping results with previous reports that more fiber QTLs were mapped in Dt subgenome than At subgenome, the results suggested a new functional hypothesis for tetraploid cotton. After the merging of the two diploid Gossypium genomes, the At subgenome has provided most of the genes for fiber development, because it continues to function similar to its fiber producing diploid A genome ancestor. On the other hand, the Dt subgenome, with its non-fiber producing D genome ancestor, provides more transcription factors that regulate the expression of the fiber genes in the At subgenome. This hypothesis would explain previously published mapping results. At the same time, this integrated map of fiber development genes would provide a framework to clone individual full-length fiber genes, to elucidate the physiological mechanisms of the fiber differentiation, elongation, and maturation, and to systematically study the functional network of these genes that interact during the process of fiber development in the tetraploid cottons.

Meta-analysis of cotton fiber quality QTLs across diverse environments in a Gossypium hirsutum x G. barbadense RIL population

BACKGROUND

Cotton fibers (produced by Gossypium species) are the premier natural fibers for textile production. The two tetraploid species, G. barbadense (Gb) and G. hirsutum (Gh), differ significantly in their fiber properties, the former having much longer, finer and stronger fibers that are highly prized. A better understanding of the genetics and underlying biological causes of these differences will aid further improvement of cotton quality through breeding and biotechnology. We evaluated an inter-specific Gh x Gb recombinant inbred line (RIL) population for fiber characteristics in 11 independent experiments under field and glasshouse conditions. Sites were located on 4 continents and 5 countries and some locations were analyzed over multiple years.

RESULTS

The RIL population displayed a large variability for all major fiber traits. QTL analyses were performed on a per-site basis by composite interval mapping. Among the 651 putative QTLs (LOD > 2), 167 had a LOD exceeding permutation based thresholds. Coincidence in QTL location across data sets was assessed for the fiber trait categories strength, elongation, length, length uniformity, fineness/maturity, and color. A meta-analysis of more than a thousand putative QTLs was conducted with MetaQTL software to integrate QTL data from the RIL and 3 backcross populations (from the same parents) and to compare them with the literature. Although the global level of congruence across experiments and populations was generally moderate, the QTL clustering was possible for 30 trait x chromosome combinations (5 traits in 19 different chromosomes) where an effective co-localization of unidirectional (similar sign of additivity) QTLs from at least 5 different data sets was observed. Most consistent meta-clusters were identified for fiber color on chromosomes c6, c8 and c25, fineness on c15, and fiber length on c3.

CONCLUSIONS

Meta-analysis provided a reliable means of integrating phenotypic and genetic mapping data across multiple populations and environments for complex fiber traits. The consistent chromosomal regions contributing to fiber quality traits constitute good candidates for the further dissection of the genetic and genomic factors underlying important fiber characteristics, and for marker-assisted selection.

Natural diversity in flowering responses of Arabidopsis thaliana caused by variation in a tandem gene array

Tandemly arrayed genes that belong to gene families characterize genomes of many organisms. Gene duplication and subsequent relaxation of selection can lead to the establishment of paralogous cluster members that may evolve along different trajectories. Here, we report on the structural variation in MADS AFFECTING FLOWERING 2 (MAF2) gene, one member of the tandemly duplicated cluster of MADS-box-containing transcription factors in Arabidopsis thaliana. The altered gene structure at the MAF2 locus is present as a moderate-frequency polymorphism in Arabidopsis and leads to the extensive diversity in transcript patterns due to alternative splicing. Rearrangements at the MAF2 locus are associated with an early flowering phenotype in BC(5) lines. The lack of suppression of flowering time in a MAF2-insertion line expressing the MAF2-specific artificial miRNA suggests that these MAF2 variants are behaving as loss-of-function alleles. The variation in gene architecture is also associated with segregation distortion, which may have facilitated the spread and the establishment of the corresponding alleles throughout the Eurasian range of the A. thaliana population.

Near-isogenic cotton germplasm lines that differ in fiber-bundle strength have temporal differences in fiber gene expression patterns as revealed by comparative high-throughput profiling

Gene expression profiles of developing cotton (Gossypium hirsutum L.) fibers from two near-isogenic lines (NILs) that differ in fiber-bundle strength, short-fiber content, and in fewer than two genetic loci were compared using an oligonucleotide microarray. Fiber gene expression was compared at five time points spanning fiber elongation and secondary cell wall (SCW) biosynthesis. Fiber samples were collected from field plots in a randomized, complete block design, with three spatially distinct biological replications for each NIL at each time point. Microarray hybridizations were performed in a loop experimental design that allowed comparisons of fiber gene expression profiles as a function of time between the two NILs. Overall, developmental expression patterns revealed by the microarray experiment agreed with previously reported cotton fiber gene expression patterns for specific genes. Additionally, genes expressed coordinately with the onset of SCW biosynthesis in cotton fiber correlated with gene expression patterns of other SCW-producing plant tissues. Functional classification and enrichment analysis of differentially expressed genes between the two NILs revealed that genes associated with SCW biosynthesis were significantly up-regulated in fibers of the high-fiber quality line at the transition stage of cotton fiber development. For independent corroboration of the microarray results, 15 genes were selected for quantitative reverse transcription PCR analysis of fiber gene expression. These analyses, conducted over multiple field years, confirmed the temporal difference in fiber gene expression between the two NILs. We hypothesize that the loci conferring temporal differences in fiber gene expression between the NILs are important regulatory sequences that offer the potential for more targeted manipulation of cotton fiber quality.

Sampling nucleotide diversity in cotton

BACKGROUND

Cultivated cotton is an annual fiber crop derived mainly from two perennial species, Gossypium hirsutum L. or upland cotton, and G. barbadense L., extra long-staple fiber Pima or Egyptian cotton. These two cultivated species are among five allotetraploid species presumably derived monophyletically between G. arboreum and G. raimondii. Genomic-based approaches have been hindered by the limited variation within species. Yet, population-based methods are being used for genome-wide introgression of novel alleles from G. mustelinum and G. tomentosum into G. hirsutum using combinations of backcrossing, selfing, and inter-mating. Recombinant inbred line populations between genetics standards TM-1, (G. hirsutum) x 3-79 (G. barbadense) have been developed to allow high-density genetic mapping of traits.

RESULTS

This paper describes a strategy to efficiently characterize genomic variation (SNPs and indels) within and among cotton species. Over 1000 SNPs from 270 loci and 279 indels from 92 loci segregating in G. hirsutum and G. barbadense were genotyped across a standard panel of 24 lines, 16 of which are elite cotton breeding lines and 8 mapping parents of populations from six cotton species. Over 200 loci were genetically mapped in a core mapping population derived from TM-1 and 3-79 and in G. hirsutum breeding germplasm.

CONCLUSION

In this research, SNP and indel diversity is characterized for 270 single-copy polymorphic loci in cotton. A strategy for SNP discovery is defined to pre-screen loci for copy number and polymorphism. Our data indicate that the A and D genomes in both diploid and tetraploid cotton remain distinct from each such that paralogs can be distinguished. This research provides mapped DNA markers for intra-specific crosses and introgression of exotic germplasm in cotton.

Using three overlapped RILs to dissect genetically clustered QTL for fiber strength on Chro.D8 in Upland cotton

Fiber strength is an important trait among cotton fiber qualities due to ongoing changes in spinning technology. Major quantitative trait loci (QTL) for fiber quality enable molecular marker-assisted selection (MAS) to effectively improve fiber quality of cotton cultivars. We previously identified a major QTL for fiber strength derived from 7235 in Upland cotton. In the present study, in order to fine-map fiber strength QTL, we chose three recombinant inbred lines (RIL), 7TR-133, 7TR-132, and 7TR-214, developed from a cross between 7235 and TM-1 for backcrossing to TM-1 to develop three large mapping populations. Phenotypic data for fiber strength traits were collected in Nanjing (JES/NAU) and Xinjiang (BES/XJ) in 2006 and 2007. Three simple sequence repeat (SSR) genetic linkage maps on Chro.24(D8) were constructed using these three backcrossed populations. The SSR genetic maps were constructed using 907 individuals in (7TR-133 x TM-1)F(2) (Pop A), 670 in (7TR-132 x TM-1)F(2) (Pop B), and 940 in (7TR-214 x TM-1)F(2) (Pop C). The average distance between SSR loci was 0.62, 1.7, and 0.56 cM for the three maps. MapQTL 5 software detected five-clustered QTL (2.5 < LOD < 29.8) on Chro.D8 for fiber strength following analysis of three RIL backcrossed F(2)/F(2:3) progenies at JES/NAU and BES/XJ over 2 years. Five QTL for fiber strength exhibited a total phenotypic variance (PV) of 28.8-59.6%.

A new interspecific, Gossypium hirsutum x G. barbadense, RIL population: towards a unified consensus linkage map of tetraploid cotton

We report the development of a new interspecific cotton recombinant inbred line (RIL) population of 140 lines deriving from an interspecific cross between Gossypium hirsutum (Gh) and G. barbadense (Gb), using the same two parents that have served for the construction of a BC(1) map and for the marker-assisted backcross selection program underway at CIRAD. Two marker systems, microsatellites and AFLPs, were used. An important feature of the RIL population was its marked segregation distortion with a genome-wide bias to Gh alleles (parental genome ratio is 71/29). The RIL map displays an excellent colinearity with the BC(1) map, although it is severely contracted in terms of map size. Existence of 255 loci in common (between 6 and 14 per chromosome) allowed the integration of the two data sets. A consensus BC(1)-RIL map based upon 215 individuals (75 BC1 + 140 RIL) was built. It consisted of 1,745 loci, spanned 3,637 cM, intermediate between the sizes of the two component maps, and constituted a solid framework to cross align cotton maps using common markers. The new RIL population will be further exploited for fiber property QTL mapping and eQTL mapping.

Development and analysis of EST-SSRs for flax (Linum usitatissimum L.)

A set of 146,611 expressed sequence tags (ESTs) were generated from 10 flax cDNA libraries. After assembly, a total of 11,166 contigs and 11,896 singletons were mined for the presence of putative simple sequence repeats (SSRs) and yielded 806 (3.5%) non-redundant sequences which contained 851 putative SSRs. This is equivalent to one EST-SSR per 16.5 kb of sequence. Trinucleotide motifs were the most abundant (76.9%), followed by dinucleotides (13.9%). Tetra-, penta- and hexanucleotide motifs represented <10% of the SSRs identified. A total of 83 SSR motifs were identified. Motif (TTC/GAA)n was the most abundant (10.2%) followed by (CTT/AAG)n (8.7%), (TCT/AGA)n (8.6%), (CT/AG)n (6.7%) and (TC/GA)n (5.3%). A total of 662 primer pairs were designed, of which 610 primer pairs yielded amplicons in a set of 23 flax accessions. Polymorphism between the accessions was found for 248 primer pairs which detected a total of 275 EST-SSR loci. Two to seven alleles were detected per marker. The polymorphism information content value for these markers ranged from 0.08 to 0.82 and averaged 0.35. The 635 alleles detected by the 275 polymorphic EST-SSRs were used to study the genetic relationship of 23 flax accessions. Four major clusters and two singletons were observed. Sub-clusters within the main clusters correlated with the pedigree relationships amongst accessions. The EST-SSRs developed herein represent the first large-scale development of SSR markers in flax. They have potential to be used for the development of genetic and physical maps, quantitative trait loci mapping, genetic diversity studies, association mapping and fingerprinting cultivars for example.

QTL mapping in A-genome diploid Asiatic cotton and their congruence analysis with AD-genome tetraploid cotton in genus Gossypium

Asiatic cotton (Gossypium arboreum L.) is an Old World cultivated cotton species. The sinense race was planted extensively in China. Due to the advances in spinning technology during the last century, the species was replaced by the New World allotetraploid cotton G. hirsutum L. Gossypium arboreum is still grown in India and Pakistan and also used as an elite in current cotton breeding programs. In addition, G. arboreum serves as a model for genomic research in Gossypium. In the present study, we generated an A-genome diploid cotton intraspecific genetic map including 264 SSR loci with three morphological markers mapped to 13 linkage groups. The map spans 2,508.71 cM with an average distance of 9.4 cM between adjacent loci. A population containing 176 F(2:3) families was used to perform quantitative trait loci (QTL) mapping for 17 phenotypes using Multiple QTL Model (MQM) of MapQTL ver 5.0. Overall, 108 QTLs were detected on 13 chromosomes. Thirty-one QTLs for yield and its components were detected in the F2 population. Forty-one QTLs for yield and its components were detected in the F(2:3) families with a total of 43 QTLs for fiber qualities. Two QTLs for seed cotton weight/plant and lint index and three QTLs for seed index were consistently detected both in F2 and F(2:3). Most QTLs for fiber qualities and yields were located at the same interval or neighboring intervals. These results indicated that the negative correlation between fiber qualities and yield traits may result from either pleiotropic effect of one gene or linkage effects of multiple closely linked genes.

SSRs and INDELs mined from the sunflower EST database: abundance, polymorphisms, and cross-taxa utility

Simple sequence repeats (SSRs) are abundant and frequently highly polymorphic in transcribed sequences and widely targeted for marker development in eukaryotes. Sunflower (Helianthus annuus) transcript assemblies were built and mined to identify SSRs and insertions-deletions (INDELs) for marker development, comparative mapping, and other genomics applications in sunflower. We describe the spectrum and frequency of SSRs identified in the sunflower EST database, a catalog of 16,643 EST-SSRs, a collection of 484 EST-SSR and 43 EST-INDEL markers developed from common sunflower ESTs, polymorphisms of the markers among the parents of several intraspecific and interspecific mapping populations, and the transferability of the markers to closely and distantly related species in the Compositae. Of 17,904 unigenes in the transcript assembly, 1,956 (10.9%) harbored one or more SSRs with repeat counts of n > or = 5. EST-SSR markers were 1.6-fold more polymorphic among exotic than elite genotypes and 0.7-fold less polymorphic than non-genic SSR markers. Of 466 EST-SSR or INDEL markers screened for cross-species amplification and polymorphisms, 413 (88.6%) amplified alleles from one or more wild species (H. argophyllus, H. tuberosus, H. anomalus, H. paradoxus, and H. deserticola), whereas 69 (14.8%) amplified alleles from safflower (Carthamus tinctorius) and 67 (14.4%) amplified alleles from lettuce (Lactuca sativa); hence, only a fraction were transferable to distantly related genera in the Compositae, whereas most were transferable to wild relatives of H. annuus. Several thousand additional SSRs were identified in the EST database and supply a wealth of templates for EST-SSR marker development in sunflower.

Assessing genetic diversity in Gossypium arboreum L. cultivars using genomic and EST-derived microsatellites

The cultivated diploid, Gossypium arboreum L., (A genome) is an invaluable genetic resource for improving modern tetraploid cotton (G. hirsutum L. and G. barbadense L.) cultivars. The objective of this research is to select a set of informative and robust microsatellites for studying genetic relationships among accessions of geographically diverse G. arboreum cultivars. From more than 1,500 previously developed simple sequence repeat (SSR) markers, 115 genomic (BNL) and EST-derived (MUCS and MUSS) markers were used to evaluate the allelic diversity of a core panel of G. arboreum accessions. These SSR data enabled advanced genome analyses. A set of 25 SSRs were selected based both upon their high level of informativeness (PIC > or = 0.50) and the production of clear PCR bands on agarose gels. Subsequently, 96 accessions representing a wide spectrum of diversity of G. arboreum cultivars were analyzed with these markers. The 25 SSR loci revealed 75 allelic variants (polymorphisms) ranging from 2 to 4 alleles per locus. The Neighborjoining (NJ) method, based on genetic dissimilarities, revealed that cultivars from geographically adjacent countries tend to cluster together. Outcomes of this research should be useful in decreasing redundancy of effort and in constructing a core collection of G. arboreum, important for efficient use of this genetic resource in cotton breeding.

Characteristics and analysis of simple sequence repeats in the cotton genome based on a linkage map constructed from a BC1 population between Gossypium hirsutum and G. barbadense

In the past decade, several molecular maps of cotton have been constructed using diverse DNA molecular markers and mapping populations. In this study, an interspecific linkage map of allotetraploid cotton was developed using a BC1 population ((Gossypium hirsutum x G. barbadense) x G. hirsutum). This map was genome-wide and was based entirely on simple sequence repeat (SSR) markers. Forty-four linkage groups were assigned to 26 chromosomes, with 917 loci spanning 5452.2 cM of the genome. The average distance between loci was 5.9 cM, providing uniform coverage of the A subgenome and D subgenome. Characteristics of this map were analyzed in detail, including the distributions of genomic SSRs, expressed sequence tag (EST)-SSRs, and distorted markers. Furthermore, the relationships between motif characteristics (size, type, length) and the level of polymorphism in EST-SSRs were also surveyed. The results showed that tetranucleotide and dinucleotide repeats had similar levels of polymorphism, and ACAT, AC, and ACT repeats had the highest polymorphism rates. Loci with lengths of 27 bp, 33 bp, and 24 bp were more likely to be polymorphic. This work will provide information to assist in designing future EST-SSRs.

Gene-based microsatellite development for mapping and phylogeny studies in eggplant

Background

Eggplant (Solanum melongena L.) is a member of the Solanaceae family. In spite of its widespread cultivation and nutritional and economic importance, its genome has not as yet been extensively investigated. Few analyses have been carried out to determine the genetic diversity of eggplant at the DNA level, and linkage relationships have not been well characterised. As for the other Solanaceae crop species (potato, tomato and pepper), the level of intra-specific polymorphism appears to be rather limited, and so it is important that an effort is made to develop more informative DNA markers to make progress in understanding the genetics of eggplant and to advance its breeding. The aim of the present work was to develop a set of functional microsatellite (SSR) markers, via an in silico analysis of publicly available DNA sequence.

Results

From >3,300 genic DNA sequences, 50 SSR-containing candidates suitable for primer design were recovered. Of these, 39 were functional, and were then applied to a panel of 44 accessions, of which 38 were cultivated eggplant varieties, and six were from related Solanum species. The usefulness of the SSR assays for diversity analysis and taxonomic discrimination was demonstrated by constructing a phylogeny based on SSR polymorphisms, and by the demonstration that most were also functional when tested with template from tomato, pepper and potato. As a results of BLASTN analyses, several eggplant SSRs were found to have homologous counterparts in the phylogenetically related species, which carry microsatellite motifs in the same position.

Conclusion

The set of eggplant EST-SSR markers was informative for phylogenetic analysis and genetic mapping. Since EST-SSRs lie within expressed sequence, they have the potential to serve as perfect markers for genes determining variation in phenotype. Their high level of transferability to other Solanaceae species can be used to provide anchoring points for the integration of genetic maps across species.

SSRs and INDELs mined from the sunflower EST database: abundance, polymorphisms, and cross-taxa utility

Simple sequence repeats (SSRs) are abundant and frequently highly polymorphic in transcribed sequences and widely targeted for marker development in eukaryotes. Sunflower (Helianthus annuus) transcript assemblies were built and mined to identify SSRs and insertions-deletions (INDELs) for marker development, comparative mapping, and other genomics applications in sunflower. We describe the spectrum and frequency of SSRs identified in the sunflower EST database, a catalog of 16,643 EST-SSRs, a collection of 484 EST-SSR and 43 EST-INDEL markers developed from common sunflower ESTs, polymorphisms of the markers among the parents of several intraspecific and interspecific mapping populations, and the transferability of the markers to closely and distantly related species in the Compositae. Of 17,904 unigenes in the transcript assembly, 1,956 (10.9%) harbored one or more SSRs with repeat counts of n > or = 5. EST-SSR markers were 1.6-fold more polymorphic among exotic than elite genotypes and 0.7-fold less polymorphic than non-genic SSR markers. Of 466 EST-SSR or INDEL markers screened for cross-species amplification and polymorphisms, 413 (88.6%) amplified alleles from one or more wild species (H. argophyllus, H. tuberosus, H. anomalus, H. paradoxus, and H. deserticola), whereas 69 (14.8%) amplified alleles from safflower (Carthamus tinctorius) and 67 (14.4%) amplified alleles from lettuce (Lactuca sativa); hence, only a fraction were transferable to distantly related genera in the Compositae, whereas most were transferable to wild relatives of H. annuus. Several thousand additional SSRs were identified in the EST database and supply a wealth of templates for EST-SSR marker development in sunflower.

Gene-rich islands for fiber development in the cotton genome

Cotton fiber is an economically important seed trichome and the world's leading natural fiber used in the manufacture of textiles. As a step toward elucidating the genomic organization and distribution of gene networks responsible for cotton fiber development, we investigated the distribution of fiber genes in the cotton genome. Results revealed the presence of gene-rich islands for fiber genes with a biased distribution in the tetraploid cotton (Gossypium hirsutum L.) genome that was also linked to discrete fiber developmental stages based on expression profiles. There were 3 fiber gene-rich islands associated with fiber initiation on chromosome 5, 3 islands for the early to middle elongation stage on chromosome 10, 3 islands for the middle to late elongation stage on chromosome 14, and 1 island on chromosome 15 for secondary cell wall deposition, for a total of 10 fiber gene-rich islands. Clustering of functionally related gene clusters in the cotton genome displaying similar transcriptional regulation indicates an organizational hierarchy with significant implications for the genetic enhancement of particular fiber quality traits. The relationship between gene-island distribution and functional expression profiling suggests for the first time the existence of functional coupling gene clusters in the cotton genome.

A preliminary analysis of genome structure and composition in Gossypium hirsutum

Background

Upland cotton has the highest yield, and accounts for > 95% of world cotton production. Decoding upland cotton genomes will undoubtedly provide the ultimate reference and resource for structural, functional, and evolutionary studies of the species. Here, we employed GeneTrek and BAC tagging information approaches to predict the general composition and structure of the allotetraploid cotton genome.

Results

142 BAC sequences from Gossypium hirsutum cv. Maxxa were downloaded and confirmed. These BAC sequence analysis revealed that the tetraploid cotton genome contains over 70,000 candidate genes with duplicated gene copies in homoeologous A- and D-subgenome regions. Gene distribution is uneven, with gene-rich and gene-free regions of the genome. Twenty-one percent of the 142 BACs lacked genes. BAC gene density ranged from 0 to 33.2 per 100 kb, whereas most gene islands contained only one gene with an average of 1.5 genes per island. Retro-elements were found to be a major component, first an enriched LTR/gypsy and second LTR/copia. Most LTR retrotransposons were truncated and in nested structures. In addition, 166 polymorphic loci amplified with SSRs developed from 70 BAC clones were tagged on our backbone genetic map. Seventy-five percent (125/166) of the polymorphic loci were tagged on the D-subgenome. By comprehensively analyzing the molecular size of amplified products among tetraploid G. hirsutum cv. Maxxa, acc. TM-1, and G. barbadense cv. Hai7124, and diploid G. herbaceum var. africanum and G. raimondii, 37 BACs, 12 from the A- and 25 from the D-subgenome, were further anchored to their corresponding subgenome chromosomes. After a large amount of genes sequence comparison from different subgenome BACs, the result showed that introns might have no contribution to different subgenome size in Gossypium.

Conclusion

This study provides us with the first glimpse of cotton genome complexity and serves as a foundation for tetraploid cotton whole genomesequencing in the future.

Cotton (Gossypium spp.) R2R3-MYB transcription factors SNP identification, phylogenomic characterization, chromosome localization, and linkage mapping

R2R3-MYB transcription factors of plants are involved in the regulation of trichome length and density. Several of them are differentially expressed during initiation and elongation of cotton fibers. We report sequence phylogenomic characterization of the six MYB genes, their chromosomal localization, and linkage mapping via SNP marker in AD-genome cotton (2n = 52). Phylogenetic grouping and comparison to At- and Dt-genome putative ancestral diploid species of allotetraploid cotton facilitated differentiation between genome-specific polymorphisms (GSPs) and marker-suitable locus-specific polymorphisms (LSPs). The SNP frequency averaged one per 77 bases overall, and one per 106 and 30 bases in coding and non-coding regions, respectively. SNP-based multivariate relationships conformed to independent evolution of the six MYB homoeologous loci in the four tetraploid species. Nucleotide diversity analysis indicated that the six MYB loci evolved more quickly in the Dt- than At-genome. The greater variation in the Dt-D genome comparisons than that in At-A genome comparisons showed no significant bias among synonymous substitution, non-synonymous substitution, and nucleotide change in non-coding regions. SNPs were concordantly mapped by deletion analysis and linkage mapping, which confirmed their value as candidate gene markers and indicated the reliability of the SNP discovery strategy in tetraploid cotton species. We consider that these SNPs may be useful for genetic dissection of economically important fiber and yield traits because of the role of these genes in fiber development.

Simple sequence repeat genetic linkage maps of A-genome diploid cotton (Gossypium arboreum)

This study introduces the construction of the first intraspecific genetic linkage map of the A-genome diploid cotton with newly developed simple sequence repeat (SSR) markers using 189 F(2) plants derived from the cross of two Asiatic cotton cultivars (Gossypium arboreum L.) Jianglingzhongmian x Zhejiangxiaoshanlüshu. Polymorphisms between the two parents were detected using 6 092 pairs of SSR primers. Two-hundred and sixty-eight pairs of SSR primers with better polymorphisms were picked out to analyze the F(2) population. In total, 320 polymorphic bands were generated and used to construct a linkage map with JoinMap3.0. Two-hundred and sixty-seven loci, including three phenotypic traits were mapped at a logarithms of odds ratio (LOD) > or = 3.0 on 13 linkage groups. The total length of the map was 2 508.71 cM, and the average distance between adjacent markers was 9.40 cM. Chromosome assignments were according to the association of linkages with our backbone tetraploid specific map using the 89 similar SSR loci. Comparisons among the 13 suites of orthologous linkage groups revealed that the A-genome chromosomes are largely collinear with the A(t) and D(t) sub-genome chromosomes. Chromosomes associated with inversions suggested that allopolyploidization was accompanied by homologous chromosomal rearrangement. The inter-chromosomal duplicated loci supply molecular evidence that the A-genome diploid Asiatic cotton is paleopolyploid.

An integrated genetic and physical map of homoeologous chromosomes 12 and 26 in Upland cotton (G. hirsutum L.)

Background

Upland cotton (G. hirsutum L.) is the leading fiber crop worldwide. Genetic improvement of fiber quality and yield is facilitated by a variety of genomics tools. An integrated genetic and physical map is needed to better characterize quantitative trait loci and to allow for the positional cloning of valuable genes. However, developing integrated genomic tools for complex allotetraploid genomes, like that of cotton, is highly experimental. In this report, we describe an effective approach for developing an integrated physical framework that allows for the distinguishing between subgenomes in cotton.

Results

A physical map has been developed with 220 and 115 BAC contigs for homeologous chromosomes 12 and 26, respectively, covering 73.49 Mb and 34.23 Mb in physical length. Approximately one half of the 220 contigs were anchored to the At subgenome only, while 48 of the 115 contigs were allocated to the Dt subgenome only. Between the two chromosomes, 67 contigs were shared with an estimated overall physical similarity between the two chromosomal homeologs at 40.0 %. A total of 401 fiber unigenes plus 214 non-fiber unigenes were located to chromosome 12 while 207 fiber unigenes plus 183 non-fiber unigenes were allocated to chromosome 26. Anchoring was done through an overgo hybridization approach and all anchored ESTs were functionally annotated via blast analysis.

Conclusion

This integrated genomic map describes the first pair of homoeologous chromosomes of an allotetraploid genome in which BAC contigs were identified and partially separated through the use of chromosome-specific probes and locus-specific genetic markers. The approach used in this study should prove useful in the construction of genome-wide physical maps for polyploid plant genomes including Upland cotton. The identification of Gene-rich islands in the integrated map provides a platform for positional cloning of important genes and the targeted sequencing of specific genomic regions.

A transgressive segregation factor (RKN2) in Gossypium barbadense for nematode resistance clusters with gene rkn1 in G. hirsutum

Host plant resistance is an important strategy for managing root-knot nematode (Meloidogyne incognita) in cotton (Gossypium L.). Here we report evidence for enhanced resistance in interspecific crosses resulting from transgressive segregation of clustered gene loci. Recently, a major gene, rkn1, on chromosome 11 for resistance to M. incognita in cv. Acala NemX was identified using an intraspecific G. hirsutum cross with susceptible cv. Acala SJ-2. Using interspecific crosses of Acala NemXxsusceptible G. barbadense cv. Pima S-7, F1, F2, F2:3, backcross, and testcross Acala NemXxF1 (Pima S-7xSJ-2), parental entries and populations were inoculated in greenhouse tests with M. incognita. Genetic analyses based on nematode-induced root galling and nematode egg production on roots, and molecular marker analysis of the segregating interspecific populations revealed that gene rkn1 interacted with a gene (designated as RKN2) in susceptible Pima S-7 to produce a highly resistant phenotype. RKN2 did not confer resistance in Pima S-7, but when combined with rkn1 (genotype Aa or aa), high levels of resistance were produced in the F1 and segregating F2, F3, and BC1F1 populations. One SSR marker MUCS088 was identified tightly linked to RKN2 within 4.4 cM in a NemXxF1 (Pima S-7xSJ-2) testcross population. Using mapped SSR markers and interspecific segregating populations, MUCS088 linked to the transgressive gene from the susceptible parent and was located in the vicinity of rkn1 on chromosome 11. Diverse genome analyses among A and D genome diploid and tetraploid cottons revealed that marker MUCS088 (165 and 167 bp) is derived from G. arboreum, A2 diploid genome. These results demonstrated that a highly susceptible parent contributed to nematode resistance via transgressive segregation. Derived highly resistant lines can be used as improved resistance sources in cotton breeding, and MUCS088 can be used to monitor RKN2 introgression in diverse populations. The close genomic location of the transgressive resistance determinants provides an important model system for studying transgressive segregation and epistasis in plants.

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.

Exploitation of pepper EST-SSRs and an SSR-based linkage map

As genome and cDNA sequencing projects progress, a tremendous amount of sequence information is becoming publicly available. These sequence resources can be exploited for gene discovery and marker development. Simple sequence repeat (SSR) markers are among the most useful because of their great variability, abundance, and ease of analysis. By in silico analysis of 10,232 non-redundant expressed sequence tags (ESTs) in pepper as a source of SSR markers, 1,201 SSRs were found, corresponding to one SSR in every 3.8 kb of the ESTs. Eighteen percent of the SSR-ESTs were dinucleotide repeats, 66.0% were trinucleotide, 7.7% tetranucleotide, and 8.2% pentanucleotide; AAG (14%) and AG (12.4%) motifs were the most abundant repeat types. Based on the flanking sequences of these 1,201 SSRs, 812 primer pairs that satisfied melting temperature conditions and PCR product sizes were designed. 513 SSRs (63.1%) were successfully amplified and 150 of them (29.2%) showed polymorphism between Capsicum annuum 'TF68' and C. chinense 'Habanero'. Dinucleotide SSRs and EST-SSR markers containing AC-motifs were the most polymorphic. Polymorphism increased with repeat length and repeat number. The polymorphic EST-SSRs were mapped onto the previously generated pepper linkage map, using 107 F(2) individuals from an interspecific cross of TF68 x Habanero. One-hundred and thirtynine EST-SSRs were located on the linkage map in addition to 41 previous SSRs and 63 RFLP markers, forming 14 linkage groups (LGs) and spanning 2,201.5 cM. The EST-SSR markers were distributed over all the LGs. This SSR-based map will be useful as a reference map in Capsicum and should facilitate the use of molecular markers in pepper breeding.

QTL mapping for resistance to root-knot nematodes in the M-120 RNR Upland cotton line (Gossypium hirsutum L.) of the Auburn 623 RNR source

Root-knot nematodes Meloidogyne incognita (Kofoid and White) can cause severe yield loss in cotton (Gossypium hirsutum L.). The objectives of this study were to determine the inheritance and genomic location of genes conferring root-knot nematode resistance in M-120 RNR, a highly resistant G. hirsutum line with the Auburn 623 RNR source of resistance. Utilizing two interspecific F(2) populations developed from the same M-120 RNR by Gossypium barbadense (cv. Pima S-6) cross, genome-wide scanning with RFLP markers revealed a marker on Chromosome 7 and two on Chromosome 11 showing significant association with the resistant phenotype. The association was confirmed using SSR markers with the detection of a minor and a major dominant QTL on Chromosome 7 and 11, respectively. Combined across the two populations, the major QTL on Chromosome 11 Mi-C11 had a LOD score of 19.21 (9.69 and 9.61 for Pop1 and Pop2, respectively) and accounted for 63.7% (52.6 and 65.56% for Pop1 and Pop2, respectively) of the total phenotypic variation. The minor QTL locus on Chromosome 7 Mi ( 1 ) -C07 had a LOD score of 3.48 and accounted for 7.7% of the total phenotypic variation in the combined dataset but was detected in only one population. The allele from the M-120 RNR parent contributed to increased resistance in the Mi-C11 locus, but surprisingly, the Pima S-6 allele contributed to increased resistance in the Mi-C07 locus. The M-120 RNR allele in the Mi-C11 locus, derived from the Auburn 623 RNR, is likely to have originated from the Clevewilt 6 cultivar. Results from this study indicated that the SSR marker CIR316 may replace the laborious greenhouse screening in breeding programs to identify genotypes resistant to M. incognita.

SSR mining in coffee tree EST databases: potential use of EST–SSRs as markers for the Coffea genus

Expressed sequence tags (ESTs) from Coffea canephora leaves and fruits were used to search for types and frequencies of simple sequence repeats (EST-SSRs) with a motif length of 1-6 bp. From a non-redundant (NR) EST set of 5,534 potential unigenes, 6.8% SSR-containing sequences were identified, with an average density of one SSR every 7.73 kb of EST sequences. Trinucleotide repeats were found to be the most abundant (34.34%), followed by di- (25.75%) and hexa-nucleotide (22.04%) motifs. The development of unique genic SSR markers was optimized by a computational approach which allowed us to eliminate redundancy in the original EST set and also to test the specificity of each pair of designed primers. Twenty-five EST-SSRs were developed and used to evaluate cross-species transferability in the Coffea genus. The orthology was supported by the amplicon sequence similarity and the amplification patterns. The >94% identity of flanking sequences revealed high sequence conservation across the Coffea genus. A high level of polymorphic loci was obtained regardless of the species considered (from 75% for C. liberica to 86% for C. canephora). Moreover, the polymorphism revealed by EST-SSR was similar to that exposed by genomic SSR. It is concluded that Coffea ESTs are a valuable resource for microsatellite mining. EST-SSR markers developed from C. canephora sequences can be easily transferred to other Coffea species for which very little molecular information is available. They constitute a set of conserved orthologous markers, which would be ideal for assessing genetic diversity in coffee trees as well as for cross-referencing transcribed sequences in comparative genomics studies.

CMD: a Cotton Microsatellite Database resource for Gossypium genomics

BACKGROUND

The Cotton Microsatellite Database (CMD) http://www.cottonssr.org is a curated and integrated web-based relational database providing centralized access to publicly available cotton microsatellites, an invaluable resource for basic and applied research in cotton breeding.

DESCRIPTION

At present CMD contains publication, sequence, primer, mapping and homology data for nine major cotton microsatellite projects, collectively representing 5,484 microsatellites. In addition, CMD displays data for three of the microsatellite projects that have been screened against a panel of core germplasm. The standardized panel consists of 12 diverse genotypes including genetic standards, mapping parents, BAC donors, subgenome representatives, unique breeding lines, exotic introgression sources, and contemporary Upland cottons with significant acreage. A suite of online microsatellite data mining tools are accessible at CMD. These include an SSR server which identifies microsatellites, primers, open reading frames, and GC-content of uploaded sequences; BLAST and FASTA servers providing sequence similarity searches against the existing cotton SSR sequences and primers, a CAP3 server to assemble EST sequences into longer transcripts prior to mining for SSRs, and CMap, a viewer for comparing cotton SSR maps.

CONCLUSION

The collection of publicly available cotton SSR markers in a centralized, readily accessible and curated web-enabled database provides a more efficient utilization of microsatellite resources and will help accelerate basic and applied research in molecular breeding and genetic mapping in Gossypium spp.

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.