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

Aspartyl proteases identified as candidate genes of a fiber length QTL, qFLD05, that regulates fiber length in cotton (Gossypium hirsutum L.)

KEY MESSAGE

GhAP genes were identified as the candidates involved in cotton fiber length under the scope of fine mapping a stable fiber length QTL, qFLD05. Moreover, the transcription factor GhWRKY40 positively regulated GhAP3 to decrease fiber length. Fiber length (FL) is an economically important fiber quality trait. Although several genes controlling cotton fiber development have been identified, our understanding of this process remains limited. In this study, an FL QTL (qFLD05) was fine-mapped to a 216.9-kb interval using a secondary F2:3 population derived from the upland hybrid cultivar Ji1518. This mapped genomic segment included 15 coding genes, four of which were annotated as aspartyl proteases (GhAP1-GhAP4). GhAPs were identified as candidates for qFLD05 as the sequence variations in GhAPs were associated with FL deviations in the mapping population, and functional validation of GhAP3 and GhAP4 indicated a longer FL following decreases in their expression levels through virus-induced gene silencing (VIGS). Subsequently, the potential involvement of GhWRKY40 in the regulatory network was revealed: GhWRKY40 positively regulated GhAP3's expression according to transcriptional profiling, VIGS, yeast one-hybrid assays and dual-luciferase experiments. Furthermore, alterations in the expression of the eight previously reported cotton FL-responsive genes from the above three VIGS lines (GhAP3, GhAP4 and GhWRKY40) implied that MYB5_A12 was involved in the GhWRKY40-GhAP network. In short, we unveiled the unprecedented FL regulation roles of GhAPs in cotton, which was possibly further regulated by GhWRKY40. These findings will reveal the genetic basis of FL development associated with qFLD05 and be beneficial for the marker-assisted selection of long-staple cotton.

Development and utility of SSR markers based on Brassica sp. whole-genome in triangle of U

Introduction

Simple sequence repeats (SSR), also known as microsatellites, are crucial molecular markers in both animals and plants. Despite extensive previous research on SSRs, the development of microsatellite markers in Brassica crops remains limited and inefficient.

Methods

Krait software was used to identify microsatellites by genome-wide and marker development based on three recently sequenced basic species of Brassica crops in the triangle of U (Brassica rapa, B. nigra and B. oleracea), as well as three allotetraploids (B. juncea, B. napus and B. carinata) using public databases. Subsequently, the primers and the characteristics of microsatellites for most of them were accordingly designed on each chromosome of each of the six Brassica species, and their physical locations were identified,and the cross-transferability of primers have been carried out. In addition, a B-genome specific SSR marker was screened out.

Results

A total of 79341, 92089, 125443, 173964, 173604, and 222160 SSR loci have been identified from the whole genome sequences of Brassica crops within the triangle of U crops, B. rapa (AA), B. nigra (BB), B. oleracea (CC), B. napus (AACC), B. juncea (AABB) and B. carinata (BBCC), respectively. Comparing the number distribution of the three allotetraploid SSR loci in the three subgenomes AA, BB and CC, results indicate that the allotetraploid species have significant reduction in the number of SSR loci in the genome compared with their basic diploid counterparts. Moreover, we compared the basic species with their corresponding varieties, and found that the microsatellite characters between the allotetraploids and their corresponding basic species were very similar or almost identical. Subsequently, each of the 40 SSR primers was employed to investigate the polymorphism potential of B. rapa (85.27%), B. nigra (81.33%) and B. oleracea (73.45%), and B. rapa was found to have a higher cross-transfer rate among the basic species in the triangle of U. Meanwhile, a B-genome specific SSR marker, BniSSR23228 possessing the (AAGGA)3 sequence characteristics was obtained, and it located in chromosome B3 with a total length of 97 bp.

Discussion

In this study, results suggest that the pattern of distribution may be highly conserved during the differentiation of basic Brassica species and their allotetraploid counterparts. Our data indicated that the allotetraploidization process resulted in a significant reduction in SSR loci in the three subgenomes AA, BB and CC. The reasons may be partial gene dominated chromosomal homologous recombination and rearrangement during the evolution of basic diploid species into allotetraploids. This study provides a basis for future genomics and genetic research on the relatedness of Brassica species.

Transcriptomic datasets of Verticillium wilt resistant and non-resistant Gossypium barbadense varieties during pathogen inoculation

Cotton is a significant cash crop and the primary source of natural fiber globally. Among the numerous diseases encountered in cotton production, Verticillium wilt is one of the most serious, caused by the pathogen Verticillium dahliae (V. dahliae). Unfortunately, there are no effective targeted methods to combat this disease. Genomic resources for Verticillium wilt resistance primarily exist in Gossypium barbadense (G. barbadense). Regrettably, there have been limited transcriptomic comparisons between V. dahliae-resistant and -susceptible varieties of G. barbadense due to the scarcity of susceptible resources. In this study, we conducted a transcriptome analysis on both V. dahliae-resistant and -susceptible varieties of G. barbadense at the 0, 12, 24 and 48 hours after V. dahliae inoculation. This comparative transcriptome analysis yielded high-quality data and offered new insights into the molecular mechanisms underlying cotton's resistance against this destructive pathogen.

Occurrence of simple sequence repeats in cDNA sequences of safflower (Carthamus tinctorius) reveals the importance of SSR-containing genes for cell biology and dynamic response to environmental cues

Safflower (Carthamus tinctorius) is a diploid crop plant belonging to the family Asteraceae and is well known as one of important oilseed crops due to edible oil containing unsaturated fatty acids. In recent years it is gaining increased attention for food, pharmaceutical and industrial uses, and hence the updating its breeding methods is necessary. Genic simple sequence repeats (SSRs) in addition of being desire molecular markers, are supposed to influence gene function and the respective phenotype. This study aimed to identify SSRs in cDNA sequences and further analysis of the functional features of the SSR-containing genes to elucidate their role in biological and cellular processes. We identified 1,841 SSR regions in 1,667 cDNA sequences. Among all types of repeats, trinucleotide repeats were the most abundant (35.7%), followed by hexanucleotide (29.6%) and dinucleotide repeats (22.0%). Thirty five SSR primer pairs were validated by PCR reaction, detected a high rate of polymorphism (>57%) among safflower accessions, physically mapped on safflower genome and could clearly discriminate the cultivated accessions from wild relatives. The cDNA-derived SSR markers are suitable for evaluation of genetic diversity, linkage and association mapping studies and genome-based breeding programmes. Occurrence of SSR repeats in biologically-important classes of proteins such as kinases, transferases and transcription factors was inferred from functional analyses, which along with variability of their repeat copies, can endow the cell and whole organism the flexibility of facing with continuously changing environment, and indicate a structure-based evolution mechanism of the genome which acts as an up-to-dating tool for the cell and whole origanism, which is realized in GO terms such as involvement of most SSR-containing genes in biological, cellular and metabolic processes, especially in response to stimulus, response to stress, interaction to other organisms and defense responses.

Genome-wide quantitative trait loci reveal the genetic basis of cotton fibre quality and yield-related traits in a Gossypium hirsutum recombinant inbred line population

Cotton is widely cultivated globally because it provides natural fibre for the textile industry and human use. To identify quantitative trait loci (QTLs)/genes associated with fibre quality and yield, a recombinant inbred line (RIL) population was developed in upland cotton. A consensus map covering the whole genome was constructed with three types of markers (8295 markers, 5197.17 centimorgans (cM)). Six fibre yield and quality traits were evaluated in 17 environments, and 983 QTLs were identified, 198 of which were stable and mainly distributed on chromosomes 4, 6, 7, 13, 21 and 25. Thirty-seven QTL clusters were identified, in which 92.8% of paired traits with significant medium or high positive correlations had the same QTL additive effect directions, and all of the paired traits with significant medium or high negative correlations had opposite additive effect directions. In total, 1297 genes were discovered in the QTL clusters, 414 of which were expressed in two RNA-Seq data sets. Many genes were discovered, 23 of which were promising candidates. Six important QTL clusters that included both fibre quality and yield traits were identified with opposite additive effect directions, and those on chromosome 13 (qClu-chr13-2) could increase fibre quality but reduce yield; this result was validated in a natural population using three markers. These data could provide information about the genetic basis of cotton fibre quality and yield and help cotton breeders to improve fibre quality and yield simultaneously.

Transcriptome analysis identified aberrant gene expression in pollen developmental pathways leading to CGMS in cotton (Gossypium hirsutum L.)

Male sterility (induced or natural) is a potential tool for commercial hybrid seed production in different crops. Despite numerous endeavors to understand the physiological, hereditary, and molecular cascade of events governing CMS in cotton, the exact biological process controlling sterility and fertility reconstruction remains obscure. During current study, RNA-Seq using Ion Torrent S5 platform is carried out to identify 'molecular portraits' in floral buds among the Cytoplasmic Genic Male Sterility (CGMS) line, its near-isogenic maintainer, and restorer lines. A total of 300, 438 and 455 genes were differentially expressed in CGMS, Maintainer, and Restorer lines respectively. The functional analysis using AgriGo revealed suppression in the pathways involved in biogenesis and metabolism of secondary metabolites which play an important role in pollen and anther maturation. Enrichment analysis showed dearth related to pollen and anther's development in sterile line, including anomalous expression of genes and transcription factors that have a role in the development of the reproductive organ, abnormal cytoskeleton formation, defects in cell wall formation. The current study found aberrant expression of DYT1, AMS and cytochrome P450 genes involved in tapetum formation, pollen development, pollen exine and anther cuticle formation associated to male sterility as well as fertility restoration of CGMS. In the current study, more numbers of DEGs were found on Chromosome D05 and A05 as compared to other chromosomes. Expression pattern analysis of fourteen randomly selected genes using qRT-PCR showed high concurrence with gene expression profile of RNA-Seq analysis accompanied by a strong correlation of 0.82. The present study provides an important support for future studies in identifying interaction between cyto-nuclear molecular portraits, to accelerate functional genomics and molecular breeding related to cytoplasmic male sterility studies in cotton.

Identification of Loci and Candidate Genes Responsible for Fiber Length in Upland Cotton (Gossypium hirsutum L.) via Association Mapping and Linkage Analyses

Fiber length (FL) is an important fiber quality trait in cotton. Although many fiber quality quantitative trait loci (QTL) responsible for FL have been identified, most cannot be applied to breeding programs, mainly due to unstable environments or large confidence intervals. In this study, we combined a genome-wide association study (GWAS) and linkage mapping to identify and validate high-quality QTLs responsible for FL. For the GWAS, we developed 93,250 high-quality single-nucleotide polymorphism (SNP) markers based on 355 accessions, and the FL was measured in eight different environments. For the linkage mapping, we constructed an F ₂ population from two extreme accessions. The high-density linkage maps spanned 3,848.29 cM, with an average marker interval of 1.41 cM. In total, 14 and 13 QTLs were identified in the association and linkage mapping analyses, respectively. Most importantly, a major QTL on chromosome D03 identified in both populations explained more than 10% of the phenotypic variation (PV). Furthermore, we found that a sucrose synthesis-related gene (Gh_D03G1338) was associated with FL in this QTL region. The RNA-seq data showed that Gh_D03G1338 was highly expressed during the fiber development stage, and the qRT-PCR analysis showed significant expression differences between the long fiber and short fiber varieties. These results suggest that Gh_D03G1338 may determine cotton fiber elongation by regulating the synthesis of sucrose. Favorable QTLs and candidate genes should be useful for increasing fiber quality in cotton breeding.

A high-density genetic map of extra-long staple cotton (Gossypium barbadense) constructed using genotyping-by-sequencing based single nucleotide polymorphic markers and identification of fiber traits-related QTL in a recombinant inbred line population

Background

Gossypium barbadense (Sea Island, Egyptian or Pima cotton) cotton has high fiber quality, however, few studies have investigated the genetic basis of its traits using molecular markers. Genome complexity reduction approaches such as genotyping-by-sequencing have been utilized to develop abundant markers for the construction of high-density genetic maps to locate quantitative trait loci (QTLs).

Results

The Chinese G. barbadense cultivar 5917 and American Pima S-7 were used to develop a recombinant inbred line (RIL) population with 143 lines. The 143 RILs together with their parents were tested in three replicated field tests for lint yield traits (boll weight and lint percentage) and fiber quality traits (fiber length, fiber elongation, fiber strength, fiber uniformity and micronaire) and then genotyped using GBS to develop single-nucleotide polymorphism (SNP) markers. A high-density genetic map with 26 linkage groups (LGs) was constructed using 3557 GBS SNPs spanning a total genetic distance of 3076.23 cM at an average density of 1.09 cM between adjacent markers. A total of 42 QTLs were identified, including 24 QTLs on 12 LGs for fiber quality and 18 QTLs on 7 LGs for lint yield traits, with LG1 (9 QTLs), LG10 (7 QTLs) and LG14 (6 QTLs) carrying more QTLs. Common QTLs for the same traits and overlapping QTLs for different traits were detected. Each individual QTLs explained 0.97 to 20.7% of the phenotypic variation.

Conclusions

This study represents one of the first genetic mapping studies on the fiber quality and lint yield traits in a RIL population of G. barbadense using GBS-SNPs. The results provide important information for the subsequent fine mapping of QTLs and the prediction of candidate genes towards map-based cloning and marker-assisted selection in cotton.

Genetic diversity, linkage disequilibrium, and association mapping analyses of Gossypium barbadense L. germplasm

Limited polymorphism and narrow genetic base, due to genetic bottleneck through historic domestication, highlight a need for comprehensive characterization and utilization of existing genetic diversity in cotton germplasm collections. In this study, 288 worldwide Gossypium barbadense L. cotton germplasm accessions were evaluated in two diverse environments (Uzbekistan and USA). These accessions were assessed for genetic diversity, population structure, linkage disequilibrium (LD), and LD-based association mapping (AM) of fiber quality traits using 108 genome-wide simple sequence repeat (SSR) markers. Analyses revealed structured population characteristics and a high level of intra-variability (67.2%) and moderate interpopulation differentiation (32.8%). Eight percent and 4.3% of markers revealed LD in the genome of the G. barbadense at critical values of r2 ≥ 0.1 and r2 ≥ 0.2, respectively. The LD decay was on average 24.8 cM at the threshold of r2 ≥ 0.05. LD retained on average distance of 3.36 cM at the threshold of r2 ≥ 0.1. Based on the phenotypic evaluations in the two diverse environments, 100 marker loci revealed a strong association with major fiber quality traits using mixed linear model (MLM) based association mapping approach. Fourteen marker loci were found to be consistent with previously identified quantitative trait loci (QTLs), and 86 were found to be new unreported marker loci. Our results provide insights into the breeding history and genetic relationship of G. barbadense germplasm and should be helpful for the improvement of cotton cultivars using molecular breeding and omics-based technologies.

Association mapping analysis of fiber yield and quality traits in Upland cotton (Gossypium hirsutum L.)

Fiber yield and quality are the most important traits for Upland cotton (Gossypium hirsutum L.). Identifying high yield and good fiber quality genes are the prime concern of researchers in cotton breeding. Association mapping offers an alternative and powerful method for detecting those complex agronomic traits. In this study, 198 simple sequence repeats (SSRs) were used to screen markers associated with fiber yield and quality traits with 302 elite Upland cotton accessions that were evaluated in 12 locations representing the Yellow River and Yangtze River cotton growing regions of China. Three subpopulations were found after the estimation of population structure. The pair-wise kinship values varied from 0 to 0.867. Only 1.59% of the total marker locus pairs showed significant linkage disequilibrium (LD, p < 0.001). The genome-wide LD decayed within the genetic distance of ~30 to 32 cM at r ² = 0.1, and decreased to ~1 to 2 cM at r ² = 0.2, indicating the potential for association mapping. Analysis based on a mixed linear model detected 57 significant (p < 0.01) marker-trait associations, including seven associations for fiber length, ten for fiber micronaire, nine for fiber strength, eight for fiber elongation, five for fiber uniformity index, five for fiber uniformity ratio, six for boll weight and seven for lint percent, for a total of 35 SSR markers, of which 11 markers were associated with more than one trait. Among marker-trait associations, 24 associations coincided with the previously reported quantitative trait loci (QTLs), the remainder were newly identified QTLs/genes. The QTLs identified in this study will potentially facilitate improvement of fiber yield and quality in the future cotton molecular breeding programs.

Fine mapping and candidate gene analysis of qFL-chr1, a fiber length QTL in cotton

A fiber length QTL, qFL-chr1, was fine mapped to a 0.9 cM interval of cotton chromosome 1. Two positional candidate genes showed positive correlation between gene expression level and fiber length. Prior analysis of a backcross-self mapping population derived from a cross between Gossypium hirsutum L. and G. barbadense L. revealed a QTL on chromosome 1 associated with increased fiber length (qFL-chr1), which was confirmed in three independent populations of near-isogenic introgression lines (NIILs). Here, a single NIIL, R01-40-08, was used to develop a large population segregating for the target region. Twenty-two PCR-based polymorphic markers used to genotype 1672 BC₄F₂ plants identified 432 recombinants containing breakpoints in the target region. Substitution mapping using 141 informative recombinants narrowed the position of qFL-chr1 to a 1.0-cM interval between SSR markers MUSS084 and CIR018. To exclude possible effects of non-target introgressions on fiber length, different heterozygous BC₄F₃ plants introgressed between SSR markers NAU3384 and CGR5144 were selected to develop sub-NILs. The qFL-chr1 was further mapped at 0.9-cM interval between MUSS422 and CIR018 by comparisons of sub-NIL phenotype, and increased fiber length by ~1 mm. The 2.38-Mb region between MUSS422 and CIR018 in G. barbadense contained 19 annotated genes. Expression levels of two of these genes, GOBAR07705 (encoding 1-aminocyclopropane-1-carboxylate synthase) and GOBAR25992 (encoding amino acid permease), were positively correlated with fiber length in a small F₂ population, supporting these genes as candidates for qFL-chr1.

A Genetic Map Between Gossypium hirsutum and the Brazilian Endemic G. mustelinum and Its Application to QTL Mapping

Among the seven tetraploid cotton species, little is known about transmission genetics and genome organization in Gossypium mustelinum, the species most distant from the source of most cultivated cotton, G. hirsutum. In this research, an F₂ population was developed from an interspecific cross between G. hirsutum and G. mustelinum (HM). A genetic linkage map was constructed mainly using simple sequence repeat (SSRs) and restriction fragment length polymorphism (RFLP) DNA markers. The arrangements of most genetic loci along the HM chromosomes were identical to those of other tetraploid cotton species. However, both major and minor structural rearrangements were also observed, for which we propose a parsimony-based model for structural divergence of tetraploid cottons from common ancestors. Sequences of mapped markers were used for alignment with the 26 scaffolds of the G. hirsutum draft genome, and showed high consistency. Quantitative trait locus (QTL) mapping of fiber elongation in advanced backcross populations derived from the same parents demonstrated the value of the HM map. The HM map will serve as a valuable resource for QTL mapping and introgression of G. mustelinum alleles into G. hirsutum, and help clarify evolutionary relationships between the tetraploid cotton genomes.

Structure, evolution, and comparative genomics of tetraploid cotton based on a high-density genetic linkage map

A high-density linkage map was constructed using 1,885 newly obtained loci and 3,747 previously published loci, which included 5,152 loci with 4696.03 cM in total length and 0.91 cM in mean distance. Homology analysis in the cotton genome further confirmed the 13 expected homologous chromosome pairs and revealed an obvious inversion on Chr10 or Chr20 and repeated inversions on Chr07 or Chr16. In addition, two reciprocal translocations between Chr02 and Chr03 and between Chr04 and Chr05 were confirmed. Comparative genomics between the tetraploid cotton and the diploid cottons showed that no major structural changes exist between D_T and D chromosomes but rather between A_T and A chromosomes. Blast analysis between the tetraploid cotton genome and the mixed genome of two diploid cottons showed that most AD chromosomes, regardless of whether it is from the A_T or D_T genome, preferentially matched with the corresponding homologous chromosome in the diploid A genome, and then the corresponding homologous chromosome in the diploid D genome, indicating that the diploid D genome underwent converted evolution by the diploid A genome to form the D_T genome during polyploidization. In addition, the results reflected that a series of chromosomal translocations occurred among Chr01/Chr15, Chr02/Chr14, Chr03/Chr17, Chr04/Chr22, and Chr05/Chr19.

Genome Wide SSR High Density Genetic Map Construction from an Interspecific Cross of Gossypium hirsutum × Gossypium tomentosum

A high density genetic map was constructed using F2 population derived from an interspecific cross of G. hirsutum × G. tomentosum. The map consisted of 3093 marker loci distributed across all the 26 chromosomes and covered 4365.3 cM of cotton genome with an average inter-marker distance of 1.48 cM. The maximum length of chromosome was 218.38 cM and the minimum was 122.09 cM with an average length of 167.90 cM. A sub-genome covers more genetic distance (2189.01 cM) with an average inter loci distance of 1.53 cM than D sub-genome which covers a length of 2176.29 cM with an average distance of 1.43 cM. There were 716 distorted loci in the map accounting for 23.14% and most distorted loci were distributed on D sub-genome (25.06%), which were more than on A sub-genome (21.23%). In our map 49 segregation hotspots (SDR) were distributed across the genome with more on D sub-genome as compared to A genome. Two post-polyploidization reciprocal translocations of "A2/A3 and A4/A5" were suggested by seven pairs of duplicate loci. The map constructed through these studies is one of the three densest genetic maps in cotton however; this is the first dense genome wide SSR interspecific genetic map between G. hirsutum and G. tomentosum.

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.

High Resolution Consensus Mapping of Quantitative Trait Loci for Fiber Strength, Length and Micronaire on Chromosome 25 of the Upland Cotton (Gossypium hirsutum L.)

Cotton (Gossypium hirsutum L.) is an important agricultural crop that provides renewable natural fiber resources for the global textile industry. Technological developments in the textile industry and improvements in human living standards have increased the requirement for supplies and better quality cotton. Upland cotton 0-153 is an elite cultivar harboring strong fiber strength genes. To conduct quantitative trait locus (QTL) mapping for fiber quality in 0-153, we developed a population of 196 recombinant inbred lines (RILs) from a cross between 0-153 and sGK9708. The fiber quality traits in 11 environments were measured and a genetic linkage map of chromosome 25 comprising 210 loci was constructed using this RIL population, mainly using simple sequence repeat markers and single nucleotide polymorphism markers. QTLs were identified across diverse environments using the composite interval mapping method. A total of 37 QTLs for fiber quality traits were identified on chromosome 25, of which 17 were stably expressed in at least in two environments. A stable fiber strength QTL, qFS-chr25-4, which was detected in seven environments and was located in the marker interval between CRI-SNP120491 and BNL2572, could explain 6.53%-11.83% of the observed phenotypic variations. Meta-analysis also confirmed the above QTLs with previous reports. Application of these QTLs could contribute to improving fiber quality and provide information for marker-assisted selection.

The Negative Correlation between Fiber Color and Quality Traits Revealed by QTL Analysis

Naturally existing colored cotton was far from perfection due to having genetic factors for lower yield, poor fiber quality and monotonous color. These factors posed a challenge to colored cotton breeding and innovation. To identify novel quantitative trait loci (QTL) for fiber color along with understanding of correlation between fiber color and quality in colored cotton, a RIL and two F2 populations were generated from crosses among Zong128 (Brown fiber cotton) and two white fiber cotton lines which were then analyzed in four environments. Two stable and major QTLs (qLC-7-1, qFC-7-1) for fiber lint and fuzz color were detected accounting for 16.01%-59.85% of the phenotypic variation across multiple generations and environments. Meanwhile, some minor QTLs were also identified on chromosomes 5, 14, 21 and 24 providing low phenotypic variation (<5%) from only F2 populations, not from the RILs population. Especially, a multiple-effect locus for fiber color and quality has been detected between flanking markers NAU1043 and NAU3654 on chromosome 7 (A genome) over multiple environments. Of which, qLC-7-1, qFC-7-1 were responsible for positive effects and improved fiber color in offsprings. Meanwhile, the QTLs (qFL-7-1, qFU-7-1, qFF-7-1, qFE-7-1, and qFS-7-1) for fiber quality had negative effects and explained 2.19%-8.78% of the phenotypic variation. This multiple-effect locus for fiber color and quality may reveal the negative correlation between the two types of above traits, so paving the way towards cotton genetic improvement.

Genome-wide mining, characterization, and development of microsatellite markers in gossypium species

Although much research has been conducted to characterize microsatellites and develop markers, the distribution of microsatellites remains ambiguous and the use of microsatellite markers in genomic studies and marker-assisted selection is limited. To identify microsatellites for cotton research, we mined 100,290, 83,160, and 56,937 microsatellites with frequencies of 41.2, 49.1, and 74.8 microsatellites per Mb in the recently sequenced Gossypium species: G. hirsutum, G. arboreum, and G. raimondii, respectively. The distributions of microsatellites in their genomes were non-random and were positively and negatively correlated with genes and transposable elements, respectively. Of the 77,996 developed microsatellite markers, 65,498 were physically anchored to the 26 chromosomes of G. hirsutum with an average marker density of 34 markers per Mb. We confirmed 67,880 (87%) universal and 7,705 (9.9%) new genic microsatellite markers. The polymorphism was estimated in above three species by in silico PCR and validated with 505 markers in G. hirsutum. We further predicted 8,825 polymorphic microsatellite markers within G. hirsutum acc. TM-1 and G. barbadense cv. Hai7124. In our study, genome-wide mining and characterization of microsatellites, and marker development were very useful for the saturation of the allotetraploid genetic linkage map, genome evolution studies and comparative genome mapping.

Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes

Background

SNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes.

Results

We construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis.

Conclusions

We report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Electronic supplementary material

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

Constructing a high-density linkage map for Gossypium hirsutum × Gossypium barbadense and identifying QTLs for lint percentage

To introgress the good fiber quality and yield from Gossypium barbadense into a commercial Upland cotton variety, a high-density simple sequence repeat (SSR) genetic linkage map was developed from a BC1 F1 population of Gossypium hirsutum × Gossypium barbadense. The map comprised 2,292 loci and covered 5115.16 centiMorgan (cM) of the cotton AD genome, with an average marker interval of 2.23 cM. Of the marker order for 1,577 common loci on this new map, 90.36% agrees well with the marker order on the D genome sequence genetic map. Compared with five published high-density SSR genetic maps, 53.14% of marker loci were newly discovered in this map. Twenty-six quantitative trait loci (QTLs) for lint percentage (LP) were identified on nine chromosomes. Nine stable or common QTLs could be used for marker-assisted selection. Fifty percent of the QTLs were from G. barbadense and increased LP by 1.07%-2.41%. These results indicated that the map could be used for screening chromosome substitution segments from G. barbadense in the Upland cotton background, identifying QTLs or genes from G. barbadense, and further developing the gene pyramiding effect for improving fiber yield and quality.

Exploring triacylglycerol biosynthetic pathway in developing seeds of Chia (Salvia hispanica L.): a transcriptomic approach

Chia (Salvia hispanica L.), a member of the mint family (Lamiaceae), is a rediscovered crop with great importance in health and nutrition and is also the highest known terrestrial plant source of heart-healthy omega-3 fatty acid, alpha linolenic acid (ALA). At present, there is no public genomic information or database available for this crop, hindering research on its genetic improvement through genomics-assisted breeding programs. The first comprehensive analysis of the global transcriptome profile of developing Salvia hispanica L. seeds, with special reference to lipid biosynthesis is presented in this study. RNA from five different stages of seed development was extracted and sequenced separately using the Illumina GAIIx platform. De novo assembly of processed reads in the pooled transcriptome using Trinity yielded 76,014 transcripts. The total transcript length was 66,944,462 bases (66.9 Mb), with an average length of approximately 880 bases. In the molecular functions category of Gene Ontology (GO) terms, ATP binding and nucleotide binding were found to be the most abundant and in the biological processes category, the metabolic process and the regulation of transcription-DNA-dependent and oxidation-reduction process were abundant. From the EuKaryotic Orthologous Groups of proteins (KOG) classification, the major category was “Metabolism” (31.97%), of which the most prominent class was ‘carbohydrate metabolism and transport’ (5.81% of total KOG classifications) followed by ‘secondary metabolite biosynthesis transport and catabolism’ (5.34%) and ‘lipid metabolism’ (4.57%). A majority of the candidate genes involved in lipid biosynthesis and oil accumulation were identified. Furthermore, 5596 simple sequence repeats (SSRs) were identified. The transcriptome data was further validated through confirmative PCR and qRT-PCR for select lipid genes. Our study provides insight into the complex transcriptome and will contribute to further genome-wide research and understanding of chia. The identified novel UniGenes will facilitate gene discovery and creation of genomic resource for this crop.

Distribution and evolution of cotton fiber development genes in the fibreless Gossypium raimondii genome

Cotton fiber represents the largest single cell in plants and they serve as models to study cell development. This study investigated the distribution and evolution of fiber Unigenes anchored to recombination hotspots between tetraploid cotton (Gossypium hirsutum) At and Dt subgenomes, and within a parental diploid cotton (Gossypium raimondii) D genome. Comparative analysis of At vs D and Dt vs D showed that 1) the D genome provides many fiber genes after its merger with another parental diploid cotton (Gossypium arboreum) A genome although the D genome itself does not produce any spinnable fiber; 2) similarity of fiber genes is higher between At vs D than between Dt vs D genomic hotspots. This is the first report that fiber genes have higher similarity between At and D than between Dt and D. The finding provides new insights into cotton genomic regions that would facilitate genetic improvement of natural fiber properties.

A reference genetic linkage map of apomictic Hieracium species based on expressed markers derived from developing ovule transcripts

BACKGROUND AND AIMS

Apomixis in plants generates clonal progeny with a maternal genotype through asexual seed formation. Hieracium subgenus Pilosella (Asteraceae) contains polyploid, highly heterozygous apomictic and sexual species. Within apomictic Hieracium, dominant genetic loci independently regulate the qualitative developmental components of apomixis. In H. praealtum, LOSS OF APOMEIOSIS (LOA) enables formation of embryo sacs without meiosis and LOSS OF PARTHENOGENESIS (LOP) enables fertilization-independent seed formation. A locus required for fertilization-independent endosperm formation (AutE) has been identified in H. piloselloides. Additional quantitative loci appear to influence the penetrance of the qualitative loci, although the controlling genes remain unknown. This study aimed to develop the first genetic linkage maps for sexual and apomictic Hieracium species using simple sequence repeat (SSR) markers derived from expressed transcripts within the developing ovaries.

METHODS

RNA from microdissected Hieracium ovule cell types and ovaries was sequenced and SSRs were identified. Two different F1 mapping populations were created to overcome difficulties associated with genome complexity and asexual reproduction. SSR markers were analysed within each mapping population to generate draft linkage maps for apomictic and sexual Hieracium species.

KEY RESULTS

A collection of 14 684 Hieracium expressed SSR markers were developed and linkage maps were constructed for Hieracium species using a subset of the SSR markers. Both the LOA and LOP loci were successfully assigned to linkage groups; however, AutE could not be mapped using the current populations. Comparisons with lettuce (Lactuca sativa) revealed partial macrosynteny between the two Asteraceae species.

CONCLUSIONS

A collection of SSR markers and draft linkage maps were developed for two apomictic and one sexual Hieracium species. These maps will support cloning of controlling genes at LOA and LOP loci in Hieracium and should also assist with identification of quantitative loci that affect the expressivity of apomixis. Future work will focus on mapping AutE using alternative populations.

Development of chromosome-specific markers with high polymorphism for allotetraploid cotton based on genome-wide characterization of simple sequence repeats in diploid cottons (Gossypium arboreum L. and Gossypium raimondii Ulbrich)

Background

Tetraploid cotton contains two sets of homologous chromosomes, the At- and Dt-subgenomes. Consequently, many markers in cotton were mapped to multiple positions during linkage genetic map construction, posing a challenge to anchoring linkage groups and mapping economically-important genes to particular chromosomes. Chromosome-specific markers could solve this problem. Recently, the genomes of two diploid species were sequenced whose progenitors were putative contributors of the At- and Dt-subgenomes to tetraploid cotton. These sequences provide a powerful tool for developing chromosome-specific markers given the high level of synteny among tetraploid and diploid cotton genomes. In this study, simple sequence repeats (SSRs) on each chromosome in the two diploid genomes were characterized. Chromosome-specific SSRs were developed by comparative analysis and proved to distinguish chromosomes.

Results

A total of 200,744 and 142,409 SSRs were detected on the 13 chromosomes of Gossypium arboreum L. and Gossypium raimondii Ulbrich, respectively. Chromosome-specific SSRs were obtained by comparing SSR flanking sequences from each chromosome with those from the other 25 chromosomes. The average was 7,996 per chromosome. To confirm their chromosome specificity, these SSRs were used to distinguish two homologous chromosomes in tetraploid cotton through linkage group construction. The chromosome-specific SSRs and previously-reported chromosome markers were grouped together, and no marker mapped to another homologous chromosome, proving that the chromosome-specific SSRs were unique and could distinguish homologous chromosomes in tetraploid cotton. Because longer dinucleotide AT-rich repeats were the most polymorphic in previous reports, the SSRs on each chromosome were sorted by motif type and repeat length for convenient selection. The primer sequences of all chromosome-specific SSRs were also made publicly available.

Conclusion

Chromosome-specific SSRs are efficient tools for chromosome identification by anchoring linkage groups to particular chromosomes during genetic mapping and are especially useful in mapping of qualitative-trait genes or quantitative trait loci with just a few markers. The SSRs reported here will facilitate a number of genetic and genomic studies in cotton, including construction of high-density genetic maps, positional gene cloning, fingerprinting, and genetic diversity and comparative evolutionary analyses among Gossypium species.

Electronic supplementary material

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

Molecular characterization of the Gossypium Diversity Reference Set of the US National Cotton Germplasm Collection

A core marker set containing markers developed to be informative within a single commercial cotton species can elucidate diversity structure within a multi-species subset of the Gossypium germplasm collection. An understanding of the genetic diversity of cotton (Gossypium spp.) as represented in the US National Cotton Germplasm Collection is essential to develop strategies for collecting, conserving, and utilizing these germplasm resources. The US collection is one of the largest world collections and includes not only accessions with improved yield and fiber quality within cultivated species, but also accessions possessing sources of abiotic and biotic stress resistance often found in wild species. We evaluated the genetic diversity of a subset of 272 diploid and 1,984 tetraploid accessions in the collection (designated the Gossypium Diversity Reference Set) using a core set of 105 microsatellite markers. Utility of the core set of markers in differentiating intra-genome variation was much greater in commercial tetraploid genomes (99.7 % polymorphic bands) than in wild diploid genomes (72.7 % polymorphic bands), and may have been influenced by pre-selection of markers for effectiveness in the commercial species. Principal coordinate analyses revealed that the marker set differentiated interspecific variation among tetraploid species, but was only capable of partially differentiating among species and genomes of the wild diploids. Putative species-specific marker bands in G. hirsutum (73) and G. barbadense (81) were identified that could be used for qualitative identification of misclassifications, redundancies, and introgression within commercial tetraploid species. The results of this broad-scale molecular characterization are essential to the management and conservation of the collection and provide insight and guidance in the use of the collection by the cotton research community in their cotton improvement efforts.

Mapping genomic loci for cotton plant architecture, yield components, and fiber properties in an interspecific (Gossypium hirsutum L. × G. barbadense L.) RIL population

A quantitative trait locus (QTL) mapping was conducted to better understand the genetic control of plant architecture (PA), yield components (YC), and fiber properties (FP) in the two cultivated tetraploid species of cotton (Gossypium hirsutum L. and G. barbadense L.). One hundred and fifty-nine genomic regions were identified on a saturated genetic map of more than 2,500 SSR and SNP markers, constructed with an interspecific recombinant inbred line (RIL) population derived from the genetic standards of the respective cotton species (G. hirsutum acc. TM-1 × G. barbadense acc. 3-79). Using the single nonparametric and MQM QTL model mapping procedures, we detected 428 putative loci in the 159 genomic regions that confer 24 cotton traits in three diverse production environments [College Station F&B Road (FB), TX; Brazos Bottom (BB), TX; and Shafter (SH), CA]. These putative QTL loci included 25 loci for PA, 60 for YC, and 343 for FP, of which 3, 12, and 60, respectively, were strongly associated with the traits (LOD score ≥ 3.0). Approximately 17.7 % of the PA putative QTL, 32.9 % of the YC QTL, and 48.3 % of the FP QTL had trait associations under multiple environments. The At subgenome (chromosomes 1-13) contributed 72.7 % of loci for PA, 46.2 % for YC, and 50.4 % for FP while the Dt subgenome (chromosomes 14-26) contributed 27.3 % of loci for PA, 53.8 % for YC, and 49.6 % for FP. The data obtained from this study augment prior evidence of QTL clusters or gene islands for specific traits or biological functions existing in several non-homoeologous cotton chromosomes. DNA markers identified in the 159 genomic regions will facilitate further dissection of genetic factors underlying these important traits and marker-assisted selection in cotton.

Genetic diversity and population structure among pea (Pisum sativum L.) cultivars as revealed by simple sequence repeat and novel genic markers

Field pea (Pisum sativum L.) is an important cool season legume crop widely grown around the world. This research provides a basis for selection of pea germplasm across geographical regions in current and future breeding and genetic mapping efforts for pea improvement. Eleven novel genic markers were developed from pea expressed sequence tag (EST) sequences having significant similarity with gene calls from Medicago truncatula spanning at least one intron. In this study, 96 cultivars widely grown or used in breeding programs in the USA and Canada were analyzed for genetic diversity using 31 microsatellite or simple sequence repeat (SSR) and 11 novel EST-derived genic markers. The polymorphic information content varied from 0.01-0.56 among SSR markers and 0.04-0.43 among genic markers. The results showed that SSR and EST-derived genic markers displayed one or more highly reproducible, multi-allelic, and easy to score loci ranging from 200 to 700 bp in size. Genetic diversity was assessed through unweighted neighbor-joining method, and 96 varieties were grouped into three main clusters based on the dissimilarity matrix. Four subpopulations were determined through STRUCTURE analysis with no significant geographic separation of the subpopulations. The findings of the present study can be used to select diverse genotypes to be used as parents of crosses aimed for breeding improved pea cultivars.

Genetic mapping of fiber color genes on two brown cotton cultivars in Xinjiang

In the present study, genetic linkage analysis was carried out to map the fiber color loci Lc 1 and Lc 2 on two brown cotton cultivars with SSR and EST-SSR markers in the reference map by F2 segregation populations. The Lc 1 locus carried by Xincaimian6 (Gossypium hirsutum L.) was flanked by the marker NAU2862 and NAU1043 on the long arm of Chromosome 07, with genetic distance 7.8 cM and 3.8 cM, respectively. The Lc 2 carried by Xincaimian 5 (Gossypium hirsutum L.) was flanked by the marker NAU5433 and NAU2968 on the short arm of Chromosome 06, with genetic distance 4.4 cM and 7.4 cM respectively. Moreover, the marker NAU3735 and marker NAU5434 co-segregated with the Lc 1 and the Lc 2 locus, respectively. The results of marker association studies with these two loci provides the basic information for the final isolation of these important genes in colored cotton, and these linkage markers also could facilitate application of marker assisted selection in the future.

Large scale in-silico identification and characterization of simple sequence repeats (SSRs) from de novo assembled transcriptome of Catharanthus roseus (L.) G. Don

Transcriptomic data of C. roseus offering ample sequence resources for providing better insights into gene diversity: large resource of genic SSR markers to accelerate genomic studies and breeding in Catharanthus . Next-generation sequencing is an efficient system for generating high-throughput complete transcripts/genes and developing molecular markers. We present here the transcriptome sequencing of a 26-day-old Catharanthus roseus seedling tissue using Illumina GAIIX platform that resulted in a total of 3.37 Gb of nucleotide sequence data comprising 29,964,104 reads which were de novo assembled into 26,581 unigenes. Based on similarity searches 58 % of the unigenes were annotated of which 13,580 unique transcripts were assigned 5016 gene ontology terms. Further, 7,687 of the unigenes were found to have Cluster of Orthologous Group classifications, and 4,006 were assigned to 289 Kyoto Encyclopedia of Genes and Genome pathways. Also, 5,221 (19.64 %) of transcripts were distributed to 81 known transcription factor (TF) families. In-silico analysis of the transcriptome resulted in identification of 11,004 SSRs in 26.62 % transcripts from which 2,520 SSR markers were designed which exhibited a non-random pattern of distribution. The most abundant was the trinucleotide repeats (AAG/CTT) followed by the dinucleotide repeats (AG/CT). Location specific analysis of SSRs revealed that SSRs were preferentially associated with the 5'-UTRs with a predicted role in regulation of gene expression. A PCR validation of a set of 48 primers revealed 97.9 % successful amplification, and 76.6 % of them showed polymorphism across different Catharanthus species as well as accessions of C. roseus. In summary, this study will provide an insight into understanding the seedling development and resources for novel gene discovery and SSR development for utilization in marker-assisted selective breeding in C. roseus.

Characterization of Withania somnifera leaf transcriptome and expression analysis of pathogenesis-related genes during salicylic acid signaling

Withania somnifera (L.) Dunal is a valued medicinal plant with pharmaceutical applications. The present study was undertaken to analyze the salicylic acid induced leaf transcriptome of W. somnifera. A total of 45.6 million reads were generated and the de novo assembly yielded 73,523 transcript contig with average transcript contig length of 1620 bp. A total of 71,062 transcripts were annotated and 53,424 of them were assigned GO terms. Mapping of transcript contigs to biological pathways revealed presence of 182 pathways. Seventeen genes representing 12 pathogenesis-related (PR) families were mined from the transcriptome data and their pattern of expression post 17 and 36 hours of salicylic acid treatment was documented. The analysis revealed significant up-regulation of all families of PR genes by 36 hours post treatment except WsPR10. The relative fold expression of transcripts ranged from 1 fold to 6,532 fold. The two families of peroxidases including the lignin-forming anionic peroxidase (WsL-PRX) and suberization-associated anionic peroxidase (WsS-PRX) recorded maximum expression of 377 fold and 6532 fold respectively, while the expression of WsPR10 was down-regulated by 14 fold. Additionally, the most stable reference gene for normalization of qRT-PCR data was also identified. The effect of SA on the accumulation of major secondary metabolites of W. somnifera including withanoside V, withaferin A and withanolide A was also analyzed and an increase in content of all the three metabolites were detected. This is the first report on expression patterns of PR genes during salicylic acid signaling in W. somnifera.

Inheritance and QTL mapping of Fusarium wilt race 4 resistance in cotton

Diseases such as Fusarium wilt [Fusarium oxysporum f.sp. vasinfectum (FOV) Atk. Sny & Hans] represent expanding threats to cotton production. Integrating disease resistance into high-yielding, high-fiber quality cotton (Gossypium spp.) cultivars is one of the most important objectives in cotton breeding programs worldwide. In this study, we conducted a comprehensive analysis of gene action in cotton governing FOV race 4 resistance by combining conventional inheritance and quantitative trait loci (QTL) mapping with molecular markers. A set of diverse cotton populations was generated from crosses encompassing multiple genetic backgrounds. FOV race 4 resistance was investigated using seven parents and their derived populations: three intraspecific (G. hirsutum × G. hirsutum L. and G. barbadense × G. barbadense L.) F1 and F2; five interspecific (G. hirsutum × G. barbadense) F1 and F2; and one RIL. Parents and populations were evaluated for disease severity index (DSI) of leaves, and vascular stem and root staining (VRS) in four greenhouse and two field experiments. Initially, a single resistance gene (Fov4) model was observed in F2 populations based on inheritance of phenotypes. This single Fov4 gene had a major dominant gene action and conferred resistance to FOV race 4 in Pima-S6. The Fov4 gene appears to be located near a genome region on chromosome 14 marked with a QTL Fov4-C14 1 , which made the biggest contribution to the FOV race 4 resistance of the generated F2 progeny. Additional genetic and QTL analyses also identified a set of 11 SSR markers that indicated the involvement of more than one gene and gene interactions across six linkage groups/chromosomes (3, 6, 8, 14, 17, and 25) in the inheritance of FOV race 4 resistance. QTLs detected with minor effects in these populations explained 5-19 % of the DSI or VRS variation. Identified SSR markers for the resistance QTLs with major and minor effects will facilitate for the first time marker-assisted selection for the introgression of FOV race 4 resistance into elite cultivars during the breeding process.

Transcriptome analysis reveals ginsenosides biosynthetic genes, microRNAs and simple sequence repeats in Panax ginseng C. A. Meyer

BACKGROUND

Panax ginseng C. A. Meyer is one of the most widely used medicinal plants. Complete genome information for this species remains unavailable due to its large genome size. At present, analysis of expressed sequence tags is still the most powerful tool for large-scale gene discovery. The global expressed sequence tags from P. ginseng tissues, especially those isolated from stems, leaves and flowers, are still limited, hindering in-depth study of P. ginseng.

RESULTS

Two 454 pyrosequencing runs generated a total of 2,423,076 reads from P. ginseng roots, stems, leaves and flowers. The high-quality reads from each of the tissues were independently assembled into separate and shared contigs. In the separately assembled database, 45,849, 6,172, 4,041 and 3,273 unigenes were only found in the roots, stems, leaves and flowers database, respectively. In the jointly assembled database, 178,145 unigenes were observed, including 86,609 contigs and 91,536 singletons. Among the 178,145 unigenes, 105,522 were identified for the first time, of which 65.6% were identified in the stem, leaf or flower cDNA libraries of P. ginseng. After annotation, we discovered 223 unigenes involved in ginsenoside backbone biosynthesis. Additionally, a total of 326 potential cytochrome P450 and 129 potential UDP-glycosyltransferase sequences were predicted based on the annotation results, some of which may encode enzymes responsible for ginsenoside backbone modification. A BLAST search of the obtained high-quality reads identified 14 potential microRNAs in P. ginseng, which were estimated to target 100 protein-coding genes, including transcription factors, transporters and DNA binding proteins, among others. In addition, a total of 13,044 simple sequence repeats were identified from the 178,145 unigenes.

CONCLUSIONS

This study provides global expressed sequence tags for P. ginseng, which will contribute significantly to further genome-wide research and analyses in this species. The novel unigenes identified here enlarge the available P. ginseng gene pool and will facilitate gene discovery. In addition, the identification of microRNAs and the prediction of targets from this study will provide information on gene transcriptional regulation in P. ginseng. Finally, the analysis of simple sequence repeats will provide genetic makers for molecular breeding and genetic applications in this species.

C, Percy R, Stewart JM. Genetic diversity and population structure of cotton (Gossypium spp.) of the New World assessed by SSR markers

A global analysis of cotton (Gossypium spp.) genetic diversity is the first step to understanding its geographical distribution, dissemination, genetic relatedness, and population structure. To assess the genetic diversity and population structure in Gossypium species, 111 cotton accessions representing five allotetraploids (AD1–AD5 genomes), 23 Asiatic diploids of the Old World (A1 and A2 genomes), and 82 diploids of the New World subgenus Houzingenia (D1–D11 genomes) species were assessed using simple sequence repeats (SSR) markers with wide genome coverage. The mean genetic distance (GD) between the two most important New World tetraploid cottons (Upland (Gossypium hirsutum L.) and Pima (Gossypium barbadense L.)) was 0.39. Among the three shrub type sections (Houzingenia, Integrifolia, and Caducibracteolata) and three arborescent sections (Erioxylum, Selera, and Austroamericana), the GD ranged between 0.19 and 0.41. Phylogenetic analyses clustered all species into distinct phylogenetic groups, which were consistent with genomic origin, evolutionary history, and geographic distribution or ecotypes of these accessions, suggesting the existence of clear structured strata. With all of the genomes, the highest statistical analysis of Structure test through measurements of ad hoc (ΔK) occurred at K = 2, with group Q1 with the Old World diploid A genomes and with group Q2 with all the New World diploids of the D genome. AD genome accessions shared nearly equal alleles from both Q1 and Q2 groups. With all of the diploids of the New World D genomes, the highest value of ΔK occurred at K = 5. These results are consistent with the fundamental knowledge of tetraploid AD-genome formation and the rapid radiation of the American diploid cotton linage that took place somewhere in southwestern Mexico, followed by a differentiation–speciation during angiosperm evolution. In addition, SSR markers provide an alternative solution for distinguishing phylogenetic relationships between accessions of different ecotypes and for elucidating population structure of cottons of the New World.

Characterization of expressed sequence tags from developing fibers of Gossypium barbadense and evaluation of insertion-deletion variation in tetraploid cultivated cotton species

Background

Cotton is the leading fiber crop worldwide. Gossypium barbadense is an important species of cotton because of its extra-long staple fibers with superior luster and silkiness. However, a systematic analysis and utilization of cDNA sequences from G. barbadense fiber development remains understudied.

Results

A total of 21,079 high quality sequences were generated from two non-normalized cDNA libraries prepared by using a mixture of G. barbadense Hai7124 fibers and ovules. After assembly processing, a set of 8,653 unigenes were obtained. Of those, 7,786 were matched to known proteins and 7,316 were assigned to functional categories. The molecular functions of these unigenes were mostly related to binding and catalytic activity, and carbohydrate, amino acid, and energy metabolisms were major contributors among the subsets of metabolism. Sequences comparison between G. barbadense and G. hirsutum revealed that 8,245 unigenes from G. barbadense were detected the similarity with those released publicly in G. hirsutum, however, the remaining 408 sequences had no hits against G. hirsutum unigenes database. Furthermore, 13,275 putative ESTs InDels loci involved in the orthologous and/or homoeologous differences between/within G. barbadense and G. hirsutum were discovered by in silico analyses, and 2,160 InDel markers were developed by ESTs with more than five insertions or deletions. By gel electrophoresis combined with sequencing verification, 71.11% candidate InDel loci were reconfirmed orthologous and/or homoeologous loci polymorphisms using G. hirsutum acc TM-1 and G. barbadense cv Hai7124. Blastx result showed among 2,160 InDel loci, 81 with significant function similarity with known genes associated with secondary wall synthesis process, indicating the important roles in fiber quality in tetraploid cultivated cotton species.

Conclusion

Sequence comparisons and InDel markers development will lay the groundwork for promoting the identification of genes related to superior agronomic traits, genetic differentiation and comparative genomic studies between G. hirsutum and G. barbadense.

Genetic mapping and characteristics of genes specifically or preferentially expressed during fiber development in cotton

Cotton fiber is an ideal model to study cell elongation and cell wall construction in plants. During fiber development, some genes and proteins have been reported to be specifically or preferentially expressed. Mapping of them will reveal the genomic distribution of these genes, and will facilitate selection in cotton breeding. Based on previous reports, we designed 331 gene primers and 164 protein primers, and used single-strand conformation polymorphism (SSCP) to map and integrate them into our interspecific BC₁ linkage map. This resulted in the mapping of 57 loci representing 51 genes or proteins on 22 chromosomes. For those three markers which were tightly linked with quantitative trait loci (QTLs), the QTL functions obtained in this study and gene functions reported in previous reports were consistent. Reverse transcription-polymerase chain reaction (RT-PCR) analysis of 52 polymorphic functional primers showed that 21 gene primers and 17 protein primers had differential expression between Emian22 (Gossypium hirsutum) and 3–79 (G. barbadense). Both RT-PCR and quantitative real-time PCR (qRT-PCR) analyses of the three markers tightly linked with QTLs were consistent with QTL analysis and field experiments. Gene Ontology (GO) categorization revealed that almost all 51 mapped genes belonged to multiple categories that contribute to fiber development, indicating that fiber development is a complex process regulated by various genes. These 51 genes were all specifically or preferentially expressed during fiber cell elongation and secondary wall biosynthesis. Therefore, these functional gene-related markers would be beneficial for the genetic improvement of cotton fiber length and strength.

Mapping quantitative trait loci for lint yield and fiber quality across environments in a Gossypium hirsutum × Gossypium barbadense backcross inbred line population

Identification of stable quantitative trait loci (QTLs) across different environments and mapping populations is a prerequisite for marker-assisted selection (MAS) for cotton yield and fiber quality. To construct a genetic linkage map and to identify QTLs for fiber quality and yield traits, a backcross inbred line (BIL) population of 146 lines was developed from a cross between Upland cotton (Gossypium hirsutum) and Egyptian cotton (Gossypium barbadense) through two generations of backcrossing using Upland cotton as the recurrent parent followed by four generations of self pollination. The BIL population together with its two parents was tested in five environments representing three major cotton production regions in China. The genetic map spanned a total genetic distance of 2,895 cM and contained 392 polymorphic SSR loci with an average genetic distance of 7.4 cM per marker. A total of 67 QTLs including 28 for fiber quality and 39 for yield and its components were detected on 23 chromosomes, each of which explained 6.65-25.27% of the phenotypic variation. Twenty-nine QTLs were located on the At subgenome originated from a cultivated diploid cotton, while 38 were on the Dt subgenome from an ancestor that does not produce spinnable fibers. Of the eight common QTLs (12%) detected in more than two environments, two were for fiber quality traits including one for fiber strength and one for uniformity, and six for yield and its components including three for lint yield, one for seedcotton yield, one for lint percentage and one for boll weight. QTL clusters for the same traits or different traits were also identified. This research represents one of the first reports using a permanent advanced backcross inbred population of an interspecific hybrid population to identify QTLs for fiber quality and yield traits in cotton across diverse environments. It provides useful information for transferring desirable genes from G. barbadense to G. hirsutum using MAS.

Mapping of quantitative trait loci for oil content in cottonseed kernel

Oil content in cottonseed is a major quality trait which when improved through breeding could enhance the competitiveness of cottonseed oil among other vegetable oils. Cottonseed oil content is a quantitative trait controlled by genes in the tetraploid embryo and tetraploid maternal plant genomes, and the knowledge of quantitative trait loci (QTLs) and the genetic effects related to oil content in both genomes could facilitate the improvement in its quality and quantity. However, till date, QTL mapping and genetic analysis related to this trait in cotton have only been conducted in the tetraploid embryo genome. In the current experiment, an IF(2) population of cottonseed kernels from the random crossing of 188 intraspecific recombinant inbred lines which were derived from the hybrid of two parents, HS46 and MARCABUCAG8US-1-88, were used to simultaneously locate QTLs for oil content in the embryo and maternal plant genomes. The four QTLs found to be associated with oil content in cottonseed were: qOC-18-1 on chromosome 18; qOC-LG-11 on linkage group 11; qOC-18-2 on chromosome 18; and qOC-22 on chromosome 22. At a high selection threshold of 0.05, there was strong evidence linking the QTLs above the oil content in cottonseed. Embryo additive and dominant effects from the tetraploid embryo genome, as well as maternal additive effects from the tetraploid maternal plant genome were found to be significant contributors to genetic variation in cottonseed oil content.

Molecular diversity, genomic constitution, and QTL mapping of fiber quality by mapped SSRs in introgression lines derived from Gossypium hirsutum × G. darwinii Watt

Because the genetic basis of current upland cotton cultivars is narrow, exploring new germplasm resources and discovering novel alleles relevant to important agronomic traits have become two of the most important themes in the field of cotton research. In this study, G. darwinii Watt, a wild cotton species, was crossed with four upland cotton cultivars with desirable traits. A total of 105 introgression lines (ILs) were successfully obtained. By using 310 mapped SSRs evenly distributed across the interspecific linkage map of G. hirsutum × G. barbadense, these 105 ILs and their corresponding parents were analyzed. A total of 278 polymorphic loci were detected among the 105 ILs, and the average length of introgression segments accumulated to 333.5 cM, accounting for 6.7 % of the whole genome. These lines included many variations. However, high similarity coefficients existed between lines, even between those derived from different parents. Finally, all the ILs and their upland cotton parents were used for association mapping of fiber quality in three environments. A total of 40 SSRs were found to be associated with five fiber quality indexes (P < 0.05) with some being detected in multiple environments and traits. The contribution rate for trait variation was 6.31 % on average, ranging from 2.00 to 14.79 %. This study develops novel ILs for cotton genetics and breeding and provides the basis for future research into fine mapping of genes related to fiber quality, analyses of candidate gene association, and molecular marker-assisted breeding.

Genomes for jeans: cotton genomics for engineering superior fiber

Twenty years ago, scientists predicted that better understanding of fiber development would lead to novel ways to engineer superior cotton fiber. Advances in genetic resources, DNA markers, DNA sequence information, and gene expression data have indeed provided new insights into fiber initiation, elongation and maturation. Many exciting applications of this knowledge offer the potential to select better cotton genotypes more effectively in mainstream breeding programs or engineer genotypes with improved agronomic and/or quality traits. Here, we discuss recent progress in understanding genes involved in fiber development, and their regulation and manipulation to engineer improved fibers. Better understanding of quantitative trait loci/gene interactions that influence fiber quality and yield may help to tailor superior cotton genotypes to diverse environments.

Development and validation of genic-SSR markers in sesame by RNA-seq

Background

Sesame (Sesamum indicum L.) is one of the most important oil crops; however, a lack of useful molecular markers hinders current genetic research. We performed transcriptome sequencing of samples from different sesame growth and developmental stages, and mining of genic-SSR markers to identify valuable markers for sesame molecular genetics research.

Results

In this study, 75 bp and 100 bp paired-end RNA-seq was used to sequence 24 cDNA libraries, and 42,566 uni-transcripts were assembled from more than 260 million filtered reads. The total length of uni-transcript sequences was 47.99 Mb, and 7,324 SSRs (SSRs ≥15 bp) and 4,440 SSRs (SSRs ≥18 bp) were identified. On average, there was one genic-SSR per 6.55 kb (SSRs ≥15 bp) or 10.81 kb (SSRs ≥18 bp). Among perfect SSRs (≥18 bp), di-nucleotide motifs (48.01%) were the most abundant, followed by tri- (20.96%), hexa- (25.37%), penta- (2.97%), tetra- (2.12%), and mono-nucleotides (0.57%). The top four motif repeats were (AG/CT)n [1,268 (34.51%)], (CA/TG)n [281 (7.65%)], (AT/AT)n [215 (5.85%)], and (GAA/TTC)n [131 (3.57%)]. A total of 2,164 SSR primer pairs were identified in the 4,440 SSR-containing sequences (≥18 bp), and 300 SSR primer pairs were randomly chosen for validation. These SSR markers were amplified and validated in 25 sesame accessions (24 cultivated accessions, one wild species). 276 (92.0%) primer pairs yielded PCR amplification products in 24 cultivars. Thirty two primer pairs (11.59%) exhibited polymorphisms. Moreover, 203 primer pairs (67.67%) yielded PCR amplicons in the wild accession and 167 (60.51%) were polymorphic between species. A UPGMA dendrogram based on genetic similarity coefficients showed that the correlation between genotype and geographical source was low and that the genetic basis of sesame in China is narrow, as previously reported. The 32 polymorphic primer pairs were validated using an F₂ mapping population; 18 primer pairs exhibited polymorphisms between the parents, and 14 genic-SSRs could be integrated into 9 main linkage groups.

Conclusions

2,164 genic-SSR markers have been developed in sesame using transcriptome sequencing. 276 of 300 validated primer pairs successfully yielded PCR amplicons in 24 cultivated sesame accessions. These markers increase current SSR marker resources and will greatly benefit genetic diversity, qualitative and quantitative trait mapping and marker-assisted selection studies in sesame.

Increasing cotton genome coverage with polymorphic SSRs as revealed by SSCP

Simple sequence repeat (SSR) markers are widely used in plant genetics and breeding. However, there are many SSR markers that do not reveal polymorphism in cotton. Traditional SSR genotyping methods only provide information on product sizes. This leaves many marker polymorphism undetected, thus, lowering the utility of SSRs. In the present study, monomorphic SSRs between two mapping parents, 'Emian22' and 3-79, were subjected to single-strand conformation polymorphism (SSCP) analysis to reveal polymorphism. Of the 4194 monomorphic SSR primer pairs, 158 pairs (3.77%) showed polymorphism and revealed 174 polymorphic loci. Sequence analysis showed that the differences in PCR products between the mapping parents were solely due to base transition or transversion, which was in agreement with SSCP principles. SSCP also revealed SSRs with motifs of AT/TA and GAA/CTT were more polymorphic in dinucleotides and trinucleotides, respectively. Genetic mapping integrated 160 loci into our interspecific BC(1) linkage map, 5 of which associated with QTLs related to cotton fiber quality. The technique discussed in the present study enables us to detect polymorphism of monomorphic SSRs, and increase the utilization efficiency of the existing SSR primers.

EST-derived genic molecular markers: development and utilization for generating an advanced transcript map of chickpea

Well-saturated linkage maps especially those based on expressed sequence tag (EST)-derived genic molecular markers (GMMs) are a pre-requisite for molecular breeding. This is especially true in important legumes such as chickpea where few simple sequence repeats (SSR) and even fewer GMM-based maps have been developed. Therefore, in this study, 2,496 ESTs were generated from chickpea seeds and utilized for the development of 487 novel EST-derived functional markers which included 125 EST-SSRs, 151 intron targeted primers (ITPs), 109 expressed sequence tag polymorphisms (ESTPs), and 102 single nucleotide polymorphisms (SNPs). Whereas ESTSSRs, ITPs, and ESTPs were developed by in silico analysis of the developed EST sequences, SNPs were identified by allele resequencing and their genotyping was performedusing the Illumina GoldenGate Assay. Parental polymorphism was analyzed between C. arietinum ICC4958 and C. reticulatum PI489777, parents of the reference chickpea mapping population, using a total of 872 markers: 487 new gene-based markers developed in this study along with 385 previously published markers, of which 318 (36.5%) were found to be polymorphic and were used for genotyping. The genotypic data were integrated with the previously published data of 108 markers and an advanced linkage map was generated that contained 406 loci distributed on eight linkage groups that spanned 1,497.7 cM. The average marker density was 3.68 cM and the average number of markers per LG was 50.8. Among the mapped markers, 303 new genomic locations were defined that included 177 gene-based and 126 gSSRs (genomic SSRs) thereby producing the most advanced gene-rich map of chickpea solely based on co-dominant markers.

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 of New Candidate Gene and EST-Based Molecular Markers for Gossypium Species

New source of molecular markers accelerate the efforts in improving cotton fiber traits and aid in developing high-density integrated genetic maps. We developed new markers based on candidate genes and G. arboreum EST sequences that were used for polymorphism detection followed by genetic and physical mapping. Nineteen gene-based markers were surveyed for polymorphism detection in 26 Gossypium species. Cluster analysis generated a phylogenetic tree with four major sub-clusters for 23 species while three species branched out individually. CAP method enhanced the rate of polymorphism of candidate gene-based markers between G. hirsutum and G. barbadense. Two hundred A-genome based SSR markers were designed after datamining of G. arboreum EST sequences (Mississippi Gossypium arboreum  EST-SSR: MGAES). Over 70% of MGAES markers successfully produced amplicons while 65 of them demonstrated polymorphism between the parents of G. hirsutum and G. barbadense RIL population and formed 14 linkage groups. Chromosomal localization of both candidate gene-based and MGAES markers was assisted by euploid and hypoaneuploid CS-B analysis. Gene-based and MGAES markers were highly informative as they were designed from candidate genes and fiber transcriptome with a potential to be integrated into the existing cotton genetic and physical maps.

A high-density simple sequence repeat and single nucleotide polymorphism genetic map of the tetraploid cotton genome

Genetic linkage maps play fundamental roles in understanding genome structure, explaining genome formation events during evolution, and discovering the genetic bases of important traits. A high-density cotton (Gossypium spp.) genetic map was developed using representative sets of simple sequence repeat (SSR) and the first public set of single nucleotide polymorphism (SNP) markers to genotype 186 recombinant inbred lines (RILs) derived from an interspecific cross between Gossypium hirsutum L. (TM-1) and G. barbadense L. (3-79). The genetic map comprised 2072 loci (1825 SSRs and 247 SNPs) and covered 3380 centiMorgan (cM) of the cotton genome (AD) with an average marker interval of 1.63 cM. The allotetraploid cotton genome produced equivalent recombination frequencies in its two subgenomes (At and Dt). Of the 2072 loci, 1138 (54.9%) were mapped to 13 At-subgenome chromosomes, covering 1726.8 cM (51.1%), and 934 (45.1%) mapped to 13 Dt-subgenome chromosomes, covering 1653.1 cM (48.9%). The genetically smallest homeologous chromosome pair was Chr. 04 (A04) and 22 (D04), and the largest was Chr. 05 (A05) and 19 (D05). Duplicate loci between and within homeologous chromosomes were identified that facilitate investigations of chromosome translocations. The map augments evidence of reciprocal rearrangement between ancestral forms of Chr. 02 and 03 versus segmental homeologs 14 and 17 as centromeric regions show homeologous between Chr. 02 (A02) and 17 (D02), as well as between Chr. 03 (A03) and 14 (D03). This research represents an important foundation for studies on polyploid cottons, including germplasm characterization, gene discovery, and genome sequence assembly.

Genome-wide functional analysis of the cotton transcriptome by creating an integrated EST database

A total of 28,432 unique contigs (25,371 in consensus contigs and 3,061 as singletons) were assembled from all 268,786 cotton ESTs currently available. Several in silico approaches [comparative genomics, Blast, Gene Ontology (GO) analysis, and pathway enrichment by Kyoto Encyclopedia of Genes and Genomes (KEGG)] were employed to investigate global functions of the cotton transcriptome. Cotton EST contigs were clustered into 5,461 groups with a maximum cluster size of 196 members. A total of 27,956 indel mutants and 149,616 single nucleotide polymorphisms (SNPs) were identified from consensus contigs. Interestingly, many contigs with significantly high frequencies of indels or SNPs encode transcription factors and protein kinases. In a comparison with six model plant species, cotton ESTs show the highest overall similarity to grape. A total of 87 cotton miRNAs were identified; 59 of these have not been reported previously from experimental or bioinformatics investigations. We also predicted 3,260 genes as miRNAs targets, which are associated with multiple biological functions, including stress response, metabolism, hormone signal transduction and fiber development. We identified 151 and 4,214 EST-simple sequence repeats (SSRs) from contigs and raw ESTs respectively. To make these data widely available, and to facilitate access to EST-related genetic information, we integrated our results into a comprehensive, fully downloadable web-based cotton EST database (www.leonxie.com).

Differential gene expression and associated QTL mapping for cotton yield based on a cDNA-AFLP transcriptome map in an immortalized F2

cDNA-AFLP techniques have found new applications in recent years. Currently, the methodology is used to establish differential gene expression and construct linkage maps. In the present study, a transcriptome map based on cDNA-AFLP techniques was constructed using an immortalized F(2) (IF(2)) population of 171 lines. The lines were derived from intercrosses between 180 recombinant inbred lines (RILs) of the cotton hybrid Xiangzamian 2 (Gossypium. hirsutum L.). A total of 302 transcriptome-derived fragments (TDFs) were mapped onto 26 linkage groups that covered 2,477.06 cM in length with an average distance of 8.23 cM between two markers. Seventy-one QTL for yield and yield component traits were detected by CIM procedures based on four environments, with 13 QTL identified in at least two environments. Some TDFs co-located with yield QTL were subsequently sequenced and analyzed by online homology searches. Potential candidate genes for yield and yield component traits were found to encode proteins involved in DNA replication, transcription, translation, and biosynthesis regulation. Furthermore, genes regulating metabolic processes signal transduction, transport, and structural components of organelles were identified. Correlation analysis between expression patterns of TDFs and trait performance detected six TDFs positively correlated to both yield and yield heterosis: six TDFs positively correlated to yield, and seven TDFs to yield heterosis. These TDFs have potential for cloning the functional genes responsible for each corresponding trait and have future value in marker-assisted selection.

Microsatellite isolation and marker development in carrot - genomic distribution, linkage mapping, genetic diversity analysis and marker transferability across Apiaceae

BACKGROUND

The Apiaceae family includes several vegetable and spice crop species among which carrot is the most economically important member, with ~21 million tons produced yearly worldwide. Despite its importance, molecular resources in this species are relatively underdeveloped. The availability of informative, polymorphic, and robust PCR-based markers, such as microsatellites (or SSRs), will facilitate genetics and breeding of carrot and other Apiaceae, including integration of linkage maps, tagging of phenotypic traits and assisting positional gene cloning. Thus, with the purpose of isolating carrot microsatellites, two different strategies were used; a hybridization-based library enrichment for SSRs, and bioinformatic mining of SSRs in BAC-end sequence and EST sequence databases. This work reports on the development of 300 carrot SSR markers and their characterization at various levels.

RESULTS

Evaluation of microsatellites isolated from both DNA sources in subsets of 7 carrot F2 mapping populations revealed that SSRs from the hybridization-based method were longer, had more repeat units and were more polymorphic than SSRs isolated by sequence search. Overall, 196 SSRs (65.1%) were polymorphic in at least one mapping population, and the percentage of polymophic SSRs across F2 populations ranged from 17.8 to 24.7. Polymorphic markers in one family were evaluated in the entire F2, allowing the genetic mapping of 55 SSRs (38 codominant) onto the carrot reference map. The SSR loci were distributed throughout all 9 carrot linkage groups (LGs), with 2 to 9 SSRs/LG. In addition, SSR evaluations in carrot-related taxa indicated that a significant fraction of the carrot SSRs transfer successfully across Apiaceae, with heterologous amplification success rate decreasing with the target-species evolutionary distance from carrot. SSR diversity evaluated in a collection of 65 D. carota accessions revealed a high level of polymorphism for these selected loci, with an average of 19 alleles/locus and 0.84 expected heterozygosity.

CONCLUSIONS

The addition of 55 SSRs to the carrot map, together with marker characterizations in six other mapping populations, will facilitate future comparative mapping studies and integration of carrot maps. The markers developed herein will be a valuable resource for assisting breeding, genetic, diversity, and genomic studies of carrot and other Apiaceae.

Mapping Fusarium wilt race 1 resistance genes in cotton by inheritance, QTL and sequencing composition

Knowledge of the inheritance of disease resistance and genomic regions housing resistance (R) genes is essential to prevent expanding pathogen threats such as Fusarium wilt [Fusarium oxysporum f.sp. vasinfectum (FOV) Atk. Sny & Hans] in cotton (Gossypium spp.). We conducted a comprehensive study combining conventional inheritance, genetic and quantitative trait loci (QTL) mapping, QTL marker-sequence composition, and genome sequencing to examine the distribution, structure and organization of disease R genes to race 1 of FOV in the cotton genome. Molecular markers were applied to F(2) and recombinant inbred line (RIL) interspecific mapping populations from the crosses Pima-S7 (G. barbadense L.) × 'Acala NemX' (G. hirsutum L.) and Upland TM-1 (G. hirsutum) × Pima 3-79 (G. barbadense), respectively. Three greenhouse tests and one field test were used to obtain sequential estimates of severity index (DSI) of leaves, and vascular stem and root staining (VRS). A single resistance gene model was observed for the F(2) population based on inheritance of phenotypes. However, additional inheritance analyses and QTL mapping indicated gene interactions and inheritance from nine cotton chromosomes, with major QTLs detected on five chromosomes [Fov1-C06, Fov1-C08, (Fov1-C11 ( 1 ) and Fov1-C11 ( 2)) , Fov1-C16 and Fov1-C19 loci], explaining 8-31% of the DSI or VRS variation. The Fov1-C16 QTL locus identified in the F(2) and in the RIL populations had a significant role in conferring FOV race 1 resistance in different cotton backgrounds. Identified molecular markers may have important potential for breeding effective FOV race 1 resistance into elite cultivars by marker-assisted selection. Reconciliation between genetic and physical mapping of gene annotations from marker-DNA and new DNA sequences of BAC clones tagged with the resistance-associated QTLs revealed defenses genes induced upon pathogen infection and gene regions rich in disease-response elements, respectively. These offer candidate gene targets for Fusarium wilt resistance response in cotton and other host plants.