Microsatellite isolation and characterization in sunflower (Helianthus annuus L.)

Modifying of calcareous soil with some acidifying materials and its effect on Helianthus annuus (L.) growth

Calcareous soils contain a high concentration of calcium carbonate (CaCO₃), which influences soil properties related to plant growth. Humic acid (HA) and ammonium molybdate (AM) were added as treatments for calcareous soils at concentrations of 0.1, 0.5 and 1 g/l respectively. The pots were divided into three groups. The first set of groups were irrigated with AM, while the second set of groups were irrigated with HA. As a control, the third group was irrigated using only tap water. Many soil properties and plant characteristics were measured during the experiment. The results showed that most of the studied treatments aided to increase organic carbon of calcareous soil and improved sunflower height, leaf area and shoot and root biomass. All investigated treatments significantly enhanced carbohydrates content in the sunflower shoots, except the treatment with 0.1 g/l AM, while only the with AM (under all studied concentrations) significantly enhanced carbohydrates content in roots higher than untreated. Proteins content in the shoots and roots of sunflower significantly increased when treated only with 1 g/l HA higher than control. The amino acid content of sunflower roots enhanced when treated with 0.1 and 1 m/l HA and 0.5 g/l AM Evidently, acidifying materials enhanced the calcareous soil and increased productivity.

Sunflower Hybrid Breeding: From Markers to Genomic Selection

In sunflower, molecular markers for simple traits as, e.g., fertility restoration, high oleic acid content, herbicide tolerance or resistances to Plasmopara halstedii, Puccinia helianthi, or Orobanche cumana have been successfully used in marker-assisted breeding programs for years. However, agronomically important complex quantitative traits like yield, heterosis, drought tolerance, oil content or selection for disease resistance, e.g., against Sclerotinia sclerotiorum have been challenging and will require genome-wide approaches. Plant genetic resources for sunflower are being collected and conserved worldwide that represent valuable resources to study complex traits. Sunflower association panels provide the basis for genome-wide association studies, overcoming disadvantages of biparental populations. Advances in technologies and the availability of the sunflower genome sequence made novel approaches on the whole genome level possible. Genotype-by-sequencing, and whole genome sequencing based on next generation sequencing technologies facilitated the production of large amounts of SNP markers for high density maps as well as SNP arrays and allowed genome-wide association studies and genomic selection in sunflower. Genome wide or candidate gene based association studies have been performed for traits like branching, flowering time, resistance to Sclerotinia head and stalk rot. First steps in genomic selection with regard to hybrid performance and hybrid oil content have shown that genomic selection can successfully address complex quantitative traits in sunflower and will help to speed up sunflower breeding programs in the future. To make sunflower more competitive toward other oil crops higher levels of resistance against pathogens and better yield performance are required. In addition, optimizing plant architecture toward a more complex growth type for higher plant densities has the potential to considerably increase yields per hectare. Integrative approaches combining omic technologies (genomics, transcriptomics, proteomics, metabolomics and phenomics) using bioinformatic tools will facilitate the identification of target genes and markers for complex traits and will give a better insight into the mechanisms behind the traits.

Main and epistatic QTL analyses for Sclerotinia Head Rot resistance in sunflower

Sclerotinia Head Rot (SHR), a disease caused by Sclerotinia sclerotiorum, is one of the most limiting factors in sunflower production. In this study, we identified genomic loci associated with resistance to SHR to support the development of assisted breeding strategies. We genotyped 114 Recombinant Inbred Lines (RILs) along with their parental lines (PAC2 -partially resistant-and RHA266 -susceptible-) by using a 384 single nucleotide polymorphism (SNP) Illumina Oligo Pool Assay to saturate a sunflower genetic map. Subsequently, we tested these lines for SHR resistance using assisted inoculations with S. sclerotiorum ascospores. We also conducted a randomized complete-block assays with three replicates to visually score disease incidence (DI), disease severity (DS), disease intensity (DInt) and incubation period (IP) through four field trials (2010-2014). We finally assessed main effect quantitative trait loci (M-QTLs) and epistatic QTLs (E-QTLs) by composite interval mapping (CIM) and mixed-model-based composite interval mapping (MCIM), respectively. As a result of this study, the improved map incorporates 61 new SNPs over candidate genes. We detected a broad range of narrow sense heritability (h2) values (1.86-59.9%) as well as 36 M-QTLs and 13 E-QTLs along 14 linkage groups (LGs). On LG1, LG10, and LG15, we repeatedly detected QTLs across field trials; which emphasizes their putative effectiveness against SHR. In all selected variables, most of the identified QTLs showed high determination coefficients, associated with moderate to high heritability values. Using markers shared with previous Sclerotinia resistance studies, we compared the QTL locations in LG1, LG2, LG8, LG10, LG11, LG15 and LG16. This study constitutes the largest report of QTLs for SHR resistance in sunflower. Further studies focusing on the regions in LG1, LG10, and LG15 harboring the detected QTLs are necessary to identify causal alleles and contribute to unraveling the complex genetic basis governing the resistance.

Association mapping unveils favorable alleles for grain iron and zinc concentrations in lentil (Lens culinaris subsp. culinaris)

Lentil is a major cool-season grain legume grown in South Asia, West Asia, and North Africa. Populations in developing countries of these regions have micronutrient deficiencies; therefore, breeding programs should focus more on improving the micronutrient content of food. In the present study, a set of 96 diverse germplasm lines were evaluated at three different locations in India to examine the variation in iron (Fe) and zinc (Zn) concentration and identify simple sequence repeat (SSR) markers that associate with the genetic variation. The genetic variation among genotypes of the association mapping (AM) panel was characterized using a genetic distance-based and a general model-based clustering method. The model-based analysis identified six subpopulations, which satisfactorily explained the genetic structure of the AM panel. AM analysis identified three SSRs (PBALC 13, PBALC 206, and GLLC 563) associated with grain Fe concentration explaining 9% to 11% of phenotypic variation and four SSRs (PBALC 353, SSR 317-1, PLC 62, and PBALC 217) were associated with grain Zn concentration explaining 14%, to 21% of phenotypic variation. These identified SSRs exhibited consistent performance across locations. These candidate SSRs can be used in marker-assisted genetic improvement for developing Fe and Zn fortified lentil varieties. Favorable alleles and promising genotypes identified in this study can be utilized for lentil biofortification.

Population structure and genetic diversity characterization of a sunflower association mapping population using SSR and SNP markers

BACKGROUND

Argentina has a long tradition of sunflower breeding, and its germplasm is a valuable genetic resource worldwide. However, knowledge of the genetic constitution and variability levels of the Argentinean germplasm is still scarce, rendering the global map of cultivated sunflower diversity incomplete. In this study, 42 microsatellite loci and 384 single nucleotide polymorphisms (SNPs) were used to characterize the first association mapping population used for quantitative trait loci mapping in sunflower, along with a selection of allied open-pollinated and composite populations from the germplasm bank of the National Institute of Agricultural Technology of Argentina. The ability of different kinds of markers to assess genetic diversity and population structure was also evaluated.

RESULTS

The analysis of polymorphism in the set of sunflower accessions studied here showed that both the microsatellites and SNP markers were informative for germplasm characterization, although to different extents. In general, the estimates of genetic variability were moderate. The average genetic diversity, as quantified by the expected heterozygosity, was 0.52 for SSR loci and 0.29 for SNPs. Within SSR markers, those derived from non-coding regions were able to capture higher levels of diversity than EST-SSR. A significant correlation was found between SSR and SNP- based genetic distances among accessions. Bayesian and multivariate methods were used to infer population structure. Evidence for the existence of three different genetic groups was found consistently across data sets (i.e., SSR, SNP and SSR + SNP), with the maintainer/restorer status being the most prevalent characteristic associated with group delimitation.

CONCLUSION

The present study constitutes the first report comparing the performance of SSR and SNP markers for population genetics analysis in cultivated sunflower. We show that the SSR and SNP panels examined here, either used separately or in conjunction, allowed consistent estimations of genetic diversity and population structure in sunflower breeding materials. The generated knowledge about the levels of diversity and population structure of sunflower germplasm is an important contribution to this crop breeding and conservation.

Exploring genetic variability within lentil (Lens culinaris Medik.) and across related legumes using a newly developed set of microsatellite markers

Lentil (Lens culinaris Medik.) is an economically important grain legume, yet the genetic and genomic resources remain largely uncharacterized and unexploited in this crop. Microsatellites have become markers of choice for crop improvement applications. Hence, simple sequence repeat (SSR) markers were developed for lentil through the construction of genomic library enriched for GA/CT motifs. As a result 122 functional SSR primer pairs were developed from 151 microsatellite loci and validated in L. culinaris cv. Precoz. Thirty three SSR markers were utilized for the analysis of genetic relationships between cultivated and wild species of Lens and related legumes. A total of 123 alleles were amplified at 33 loci ranging from 2-5 alleles with an average of 3.73 alleles per locus. Polymorphic information content (PIC) for all the loci ranged from 0.13 to 0.99 with an average of 0.66 per locus. Varied levels of cross genera transferability were obtained ranging from 69.70 % across Pisum sativum to 12.12 % across Vigna radiata. The UPGMA based dendrogram was able to establish the uniqueness of each genotype and grouped them into two major clusters clearly resolving the genetic relationships within lentil and related species. The new set of SSR markers reported here were efficient and highly polymorphic and would add to the existing repertoire of lentil SSR markers to be utilized in molecular breeding. Moreover, the improved knowledge about intra- and inter-specific genetic relationships would facilitate germplasm utilization for lentil improvement.

Toward a molecular cytogenetic map for cultivated sunflower (Helianthus annuus L.) by landed BAC/BIBAC clones

Conventional karyotypes and various genetic linkage maps have been established in sunflower (Helianthus annuus L., 2n = 34). However, the relationship between linkage groups and individual chromosomes of sunflower remains unknown and has considerable relevance for the sunflower research community. Recently, a set of linkage group-specific bacterial /binary bacterial artificial chromosome (BAC/BIBAC) clones was identified from two complementary BAC and BIBAC libraries constructed for cultivated sunflower cv. HA89. In the present study, we used these linkage group-specific clones (∼100 kb in size) as probes to in situ hybridize to HA89 mitotic chromosomes at metaphase using the BAC- fluorescence in situ hybridization (FISH) technique. Because a characteristic of the sunflower genome is the abundance of repetitive DNA sequences, a high ratio of blocking DNA to probe DNA was applied to hybridization reactions to minimize the background noise. As a result, all sunflower chromosomes were anchored by one or two BAC/BIBAC clones with specific FISH signals. FISH analysis based on tandem repetitive sequences, such as rRNA genes, has been previously reported; however, the BAC-FISH technique developed here using restriction fragment length polymorphism (RFLP)−derived BAC/BIBAC clones as probes to apply genome-wide analysis is new for sunflower. As chromosome-specific cytogenetic markers, the selected BAC/BIBAC clones that encompass the 17 linkage groups provide a valuable tool for identifying sunflower cytogenetic stocks (such as trisomics) and tracking alien chromosomes in interspecific crosses. This work also demonstrates the potential of using a large-insert DNA library for the development of molecular cytogenetic resources.

sunTILL: a TILLING resource for gene function analysis in sunflower

BACKGROUND

Cultivated sunflower (Helianthus annus L.) is a globally important oilseed crop, subjected to intensive genetic and genomic studies. Although classical mutagenesis has successfully been applied to Helianthus genus in the past, we have developed the first sunflower TILLING resource.

RESULTS

To balance the maximum mutation density with an acceptable plant survival rate, a 'kill curve' analysis was first conducted with different ethylmethanesulfonate (EMS) dosages and different exposure times. According to the germination rate, a treatment with 0.7% EMS for 6 h was chosen. An M2 progeny of 3,651 fertile plants was obtained. Totally, 4.79% of the whole population showed clear aberrant phenotypes. A microsatellite analysis on a representative sample of the original seed stock and mutant lines confirmed the uniformity of the genetic background of plant material. The TILLING procedure was successfully applied to sunflower genome, initially by a CelI-nuclease mismatch cleavage assay coupled with a DNA-pooling level test. To investigate the efficiency of the mutagenic treatment, a pilot screening was carried out on 1,152 M2 lines focusing on four genes, three involved in the fatty acid biosynthetic pathway and one for downy mildew resistance. A total of 9 mutant lines were identified and confirmed by sequencing; thereby, the estimated overall mutation frequency for the pilot assay resulted to be 1/475 kb.

CONCLUSION

A first TILLING population for a high throughput identification of EMS-induced point mutations in sunflower genome has been successfully obtained. This represents a powerful tool to a better understanding of gene function in sunflower.

Cytological characterization of sunflower by in situ hybridization using homologous rDNA sequences and a BAC clone containing highly represented repetitive retrotransposon-like sequences

In the present work we report new tools for the characterization of the complete chromosome complement of sunflower (Helianthus annuus L.), using a bacterial artificial chromosome (BAC) clone containing repetitive sequences with similarity to retrotransposons and a homologous rDNA sequence isolated from the sunflower genome as probes for FISH. The rDNA signal was found in 3 pairs of chromosomes, coinciding with the location of satellites. The BAC clone containing highly represented retroelements hybridized with all the chromosome complement in FISH, and used together with the rDNA probe allowed the discrimination of all chromosome pairs of sunflower. Their distinctive distribution pattern suggests that these probes could be useful for karyotype characterization and for chromosome identification. The karyotype could be subdivided into 3 clear-cut groups of 12 metacentric pairs, 1 submetacentric pair, and 4 subtelocentric pairs, thus resolving previously described karyotype controversies. The use of BAC clones containing single sequences of specific markers and (or) genes associated with important agricultural traits represents an important tool for future locus-specific identification and physical mapping.

Genetic diversity analysis in the section Caulorrhizae (genus Arachis) using microsatellite markers

Diversity in 26 microsatellite loci from section Caulorrhizae germplasm was evaluated by using 33 accessions of A. pintoi Krapov. & W.C. Gregory and ten accessions of Arachis repens Handro. Twenty loci proved to be polymorphic and a total of 196 alleles were detected with an average of 9.8 alleles per locus. The variability found in those loci was greater than the variability found using morphological characters, seed storage proteins and RAPD markers previously used in this germplasm. The high potential of these markers to detect species-specific alleles and discriminate among accessions was demonstrated. The set of microsatellite primer pairs developed by our group for A. pintoi are useful molecular tools for evaluating Section Caulorrhizae germplasm, as well as that of species belonging to other Arachis sections.

Ontology and diversity of transcript-associated microsatellites mined from a globe artichoke EST database

Background

The globe artichoke (Cynara cardunculus var. scolymus L.) is a significant crop in the Mediterranean basin. Despite its commercial importance and its both dietary and pharmaceutical value, knowledge of its genetics and genomics remains scant. Microsatellite markers have become a key tool in genetic and genomic analysis, and we have exploited recently acquired EST (expressed sequence tag) sequence data (Composite Genome Project - CGP) to develop an extensive set of microsatellite markers.

Results

A unigene assembly was created from over 36,000 globe artichoke EST sequences, containing 6,621 contigs and 12,434 singletons. Over 12,000 of these unigenes were functionally assigned on the basis of homology with Arabidopsis thaliana reference proteins. A total of 4,219 perfect repeats, located within 3,308 unigenes was identified and the gene ontology (GO) analysis highlighted some GO term's enrichments among different classes of microsatellites with respect to their position. Sufficient flanking sequence was available to enable the design of primers to amplify 2,311 of these microsatellites, and a set of 300 was tested against a DNA panel derived from 28 C. cardunculus genotypes. Consistent amplification and polymorphism was obtained from 236 of these assays. Their polymorphic information content (PIC) ranged from 0.04 to 0.90 (mean 0.66). Between 176 and 198 of the assays were informative in at least one of the three available mapping populations.

Conclusion

EST-based microsatellites have provided a large set of de novo genetic markers, which show significant amounts of polymorphism both between and within the three taxa of C. cardunculus. They are thus well suited as assays for phylogenetic analysis, the construction of genetic maps, marker-assisted breeding, transcript mapping and other genomic applications in the species.

Genome-wide analysis of microsatellite sequence in seven filamentous fungi

Abundance of microsatellites with repeated unit lengths of 1-6 base pairs in seven fungi: Aspergillus nidulans, Coprinus cinereus, Cryptococcus neoformans (serotype A), Fusarium graminearum, Magnaporthe grisea, Neurospora crassa and Ustilago maydis were investigated on genomic scale. The results showed that each species has its specific profile for different types and different motifs of SSR loci. Ascomycetes fungi M. grisea, N. crassa and basidiomycete fungus U. maydis adopt much more microsatellites than other fungi examined. Total amount of 15,751, 14,788 and 6,854 SSR loci were observed respectively, average density is 406, 389 and 347 per Mbp sequence; overall length of SSR sequence was 0.82%, 0.95% and 0.79% of genomic sequence respectively. While ascomycetes fungus F. graminearum and A. nidulans contains the least SSRs in the genomic DNA, only 4,679 and 4,837 tracts were observed in 36 Mb and 30 Mb genomic sequence respectively. Microsatellite repeats in protein coding regions are investigated in Aspergillus nidulans, Magnaporthe grisea, and Neurospora crassa also, the results show that the difference of different types and motifs among three fungi is very little than that in whole genomic sequence. For trinucleotide repeats, overrepresent (comparing to the total base pair of protein coding region) of AGC, GGC, AGG, ACG and ACC was observed in coding region, frequencies of AAC and AAG were not difference between coding and non-coding region, AAT, AGT and ATG were underrepresent in coding region excepted for A. nidulans, in which ATG was overrepresentative.

Identification, characterization and utilization of unigene derived microsatellite markers in tea (Camellia sinensis L.)

BACKGROUND

Despite great advances in genomic technology observed in several crop species, the availability of molecular tools such as microsatellite markers has been limited in tea (Camellia sinensis L.). The development of microsatellite markers will have a major impact on genetic analysis, gene mapping and marker assisted breeding. Unigene derived microsatellite (UGMS) markers identified from publicly available sequence database have the advantage of assaying variation in the expressed component of the genome with unique identity and position. Therefore, they can serve as efficient and cost effective alternative markers in such species.

RESULTS

Considering the multiple advantages of UGMS markers, 1,223 unigenes were predicted from 2,181 expressed sequence tags (ESTs) of tea (Camellia sinensis L.). A total of 109 (8.9%) unigenes containing 120 SSRs were identified. SSR abundance was one in every 3.55 kb of EST sequences. The microsatellites mainly comprised of di (50.8%), tri (30.8%), tetra (6.6%), penta (7.5%) and few hexa (4.1%) nucleotide repeats. Among the dinucleotide repeats, (GA)n.(TC)n were most abundant (83.6%). Ninety six primer pairs could be designed form 83.5% of SSR containing unigenes. Of these, 61 (63.5%) primer pairs were experimentally validated and used to investigate the genetic diversity among the 34 accessions of different Camellia spp. Fifty one primer pairs (83.6%) were successfully cross transferred to the related species at various levels. Functional annotation of the unigenes containing SSRs was done through gene ontology (GO) characterization. Thirty six (60%) of them revealed significant sequence similarity with the known/putative proteins of Arabidopsis thaliana. Polymorphism information content (PIC) ranged from 0.018 to 0.972 with a mean value of 0.497. The average heterozygosity expected (HE) and observed (Ho) obtained was 0.654 and 0.413 respectively, thereby suggesting highly heterogeneous nature of tea. Further, test for IAM and SMM models for the UGMS loci showed excess heterozygosity and did not show any bottleneck operating in the tea population.

CONCLUSION

UGMS markers identified and characterized in this study provided insight about the abundance and distribution of SSR in the expressed genome of C. sinensis. The identification and validation of 61 new UGMS markers will not only help in intra and inter specific genetic diversity assessment but also be enriching limited microsatellite markers resource in tea. Further, the use of these markers would reduce the cost and facilitate the gene mapping and marker-aided selection in tea. Since, 36 of these UGMS markers correspond to the Arabidopsis protein sequence data with known functions will offer the opportunity to investigate the consequences of SSR polymorphism on gene functions.

Genetic mapping and annotation of genomic microsatellites isolated from globe artichoke

Cynara cardunculus includes three taxa, the globe artichoke (subsp. scolymus L. Hegi), the cultivated cardoon (var. altilis) and their progenitor, the wild cardoon (var. sylvestris). Globe artichoke is an important component of the Mediterranean rural economy, but its improvement through breeding has been rather limited and its genome organization remains largely unexplored. Here, we report the isolation of 61 new microsatellite loci which amplified a total of 208 alleles in a panel of 22 C. cardunculus genotypes. Of these, 51 were informative for linkage analysis and 39 were used to increase marker density in the available globe artichoke genetic maps. Sequence analysis of the 22 loci associated with genes showed that 9 are located within coding sequence, with the repetitive domain probably being involved in DNA binding or in protein-protein interactions. The expression of the genes associated with 9 of the 22 microsatellite loci was demonstrated by RT-PCR.

Isolation of microsatellites from Catharanthus roseus (L.) G. Don using enriched libraries

Catharanthus roseus is an indispensable source of the anticancerous alkaloids-vincristine and vinblastine, even though they are produced in trace amounts in vivo. In order to increase the yield of alkaloids, in vitro tissue culture studies are carried out which result in a large number of lines/cultures. For identification and characterization of the in vitro cultures, microsatellites in the form of STMS (Sequenced Tagged Microsatellite Sites) markers are used for identification of genetic polymorphism. STMS markers are also used for assessment of genetic diversity within natural populations as well as for construction of genetic linkage maps. Isolation of microsatellites and development of STMS markers typically involves library construction and screening, DNA sequencing, polymerase chain reaction (PCR) primer design, and PCR optimization. This chapter details two approaches for the isolation of microsatellite loci. The first approach is based on PCR using microsatellite containing primers which also have degenerate bases at the 5 cent-end that act as anchors preventing the primers from slippage to the 3 cent-end and the subsequent loss of polymorphism. The multi-locus PCR amplified product is cloned and sequenced. Though this method generates a large number of microsatellites, the major drawback is the high redundancy observed in this method. The second approach described in this chapter is based on the construction of a microsatellite enriched library which involves preferential cloning of the microsatellite enriched fraction of genomic DNA. This method therefore necessitates the isolation of microsatellites through hybridization with biotin labeled oligoprobe followed by their capture with streptavidin-coated magnetic beads. In comparison to the first approach, this approach yields less redundant clones with high microsatellite enrichment. Moreover enriched libraries are 40-60 times more efficient than the conventional small insert genomic libraries.

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.

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.

Cytomorphological and molecular diversity in backcross-derived inbred lines of sunflower (Helianthus annuus L.)

A set of 250 distinct, stable, and uniform backcross-derived inbred lines were developed in sunflower through 5 interspecific cross combinations involving 4 wild diploid annual species (Helianthus argophyllus, H. petiolaris, H. annuus, and H. debilis). The presence of the wild-species genome in these inbred lines was confirmed through higher chromosome associations (tri- and quadrivalents) at diakinesis. Maximum structural rearrangements of chromosomes were observed in lines derived from H. petiolaris. Forty morphologically diverse inbred lines along with 2 controls were subjected to measurements of phenotypic and genetic distance using 118 simple sequence repeat (SSR) markers of known map location. A total of 204 alleles were identified and the number of alleles per locus varied between 2 and 5. There were 46 unique alleles and the number of unique alleles was highest in the lines derived from the cross involving H. petiolaris. The polymorphism information content (PIC) values ranged from 0.05 to 0.575. The pair-wise comparison values based on genetic dissimilarity estimates computed using molecular marker data varied between 0.143 and 0.486 among the 42 lines. The results indicate that the sunflower gene pool could benefit from introgression of novel alleles from the latent genetic diversity present in the wild species and particularly through exploitation of the diploid annual H. petiolaris.

Identification of single nucleotide polymorphisms and analysis of linkage disequilibrium in sunflower elite inbred lines using the candidate gene approach

BACKGROUND

Association analysis is a powerful tool to identify gene loci that may contribute to phenotypic variation. This includes the estimation of nucleotide diversity, the assessment of linkage disequilibrium structure (LD) and the evaluation of selection processes. Trait mapping by allele association requires a high-density map, which could be obtained by the addition of Single Nucleotide Polymorphisms (SNPs) and short insertion and/or deletions (indels) to SSR and AFLP genetic maps. Nucleotide diversity analysis of randomly selected candidate regions is a promising approach for the success of association analysis and fine mapping in the sunflower genome. Moreover, knowledge of the distance over which LD persists, in agronomically meaningful sunflower accessions, is important to establish the density of markers and the experimental design for association analysis.

RESULTS

A set of 28 candidate genes related to biotic and abiotic stresses were studied in 19 sunflower inbred lines. A total of 14,348 bp of sequence alignment was analyzed per individual. In average, 1 SNP was found per 69 nucleotides and 38 indels were identified in the complete data set. The mean nucleotide polymorphism was moderate (theta = 0.0056), as expected for inbred materials. The number of haplotypes per region ranged from 1 to 9 (mean = 3.54 +/- 1.88). Model-based population structure analysis allowed detection of admixed individuals within the set of accessions examined. Two putative gene pools were identified (G1 and G2), with a large proportion of the inbred lines being assigned to one of them (G1). Consistent with the absence of population sub-structuring, LD for G1 decayed more rapidly (r2 = 0.48 at 643 bp; trend line, pooled data) than the LD trend line for the entire set of 19 individuals (r2 = 0.64 for the same distance).

CONCLUSION

Knowledge about the patterns of diversity and the genetic relationships between breeding materials could be an invaluable aid in crop improvement strategies. The relatively high frequency of SNPs within the elite inbred lines studied here, along with the predicted extent of LD over distances of 100 kbp (r2 approximately 0.1) suggest that high resolution association mapping in sunflower could be achieved with marker densities lower than those usually reported in the literature.

Identification of single nucleotide polymorphisms and analysis of linkage disequilibrium in sunflower elite inbred lines using the candidate gene approach

BACKGROUND

Association analysis is a powerful tool to identify gene loci that may contribute to phenotypic variation. This includes the estimation of nucleotide diversity, the assessment of linkage disequilibrium structure (LD) and the evaluation of selection processes. Trait mapping by allele association requires a high-density map, which could be obtained by the addition of Single Nucleotide Polymorphisms (SNPs) and short insertion and/or deletions (indels) to SSR and AFLP genetic maps. Nucleotide diversity analysis of randomly selected candidate regions is a promising approach for the success of association analysis and fine mapping in the sunflower genome. Moreover, knowledge of the distance over which LD persists, in agronomically meaningful sunflower accessions, is important to establish the density of markers and the experimental design for association analysis.

RESULTS

A set of 28 candidate genes related to biotic and abiotic stresses were studied in 19 sunflower inbred lines. A total of 14,348 bp of sequence alignment was analyzed per individual. In average, 1 SNP was found per 69 nucleotides and 38 indels were identified in the complete data set. The mean nucleotide polymorphism was moderate (theta = 0.0056), as expected for inbred materials. The number of haplotypes per region ranged from 1 to 9 (mean = 3.54 +/- 1.88). Model-based population structure analysis allowed detection of admixed individuals within the set of accessions examined. Two putative gene pools were identified (G1 and G2), with a large proportion of the inbred lines being assigned to one of them (G1). Consistent with the absence of population sub-structuring, LD for G1 decayed more rapidly (r2 = 0.48 at 643 bp; trend line, pooled data) than the LD trend line for the entire set of 19 individuals (r2 = 0.64 for the same distance).

CONCLUSION

Knowledge about the patterns of diversity and the genetic relationships between breeding materials could be an invaluable aid in crop improvement strategies. The relatively high frequency of SNPs within the elite inbred lines studied here, along with the predicted extent of LD over distances of 100 kbp (r2 approximately 0.1) suggest that high resolution association mapping in sunflower could be achieved with marker densities lower than those usually reported in the literature.

The development of BAC-end sequence-based microsatellite markers and placement in the physical and genetic maps of soybean

The composite map of soybean shared among Soybase, LIS and SoyGD (March 2006) contained 3,073 DNA markers in the "Locus" class. Among the markers were 1,019 class I microsatellite markers with 2-3 bp simple sequence repeats (SSRs) of >10 iterations (BARC-SSR markers). However, there were few class II SSRs (2-5 bp repeats with <10 iterations; mostly SIUC-Satt markers). The aims here were to increase the number of classes I and II SSR markers and to integrate bacterial artificial chromosome (BAC) clones onto the soybean physical map using the markers. Used was 10 Mb of BAC-end sequence (BES) derived from 13,473 reads from 7,050 clones constituting minimum tile path 2 of the soybean physical map ( http://www.soybeangenome.siu.edu ; SoyGD). Identified were 1,053 1-6 bp motif, repeat sequences, 333 from class I (>10 repeats) and 720 from class II (<10 repeats). Potential markers were shown on the MTP_SSR track at Gbrowse. Primers were designed as 20-24 bp oligomers that had Tm of 55 +/- 1 C that would generate 100-500 bp amplicons. About 853 useful primer pairs were established. Motifs were not randomly distributed with biases toward AT rich motifs. Strong biases against the GC motif and all tetra-nucleotide repeats were found. The markers discovered were useful. Among the first 135 targeted for use in genetic map improvement about 60% of class II markers and 75% of class I markers were polymorphic among on the parents of four recombinant inbred line (RIL) populations. Many of the BES-based SSRs were located on the soybean genetic map in regions with few BARC-SSR markers. Therefore, BES-based SSRs represent useful tools for genetic map development in soybean. New members of a consortium to map the markers in additional populations are invited.

Development and characterization of 248 novel microsatellite markers in turbot (Scophthalmus maximus)

The turbot is a flatfish species of great relevance to marine aquaculture in Europe. Only a limited number of microsatellites have been isolated to date in this species. To increase the number of potentially useful mapping markers, we screened simple sequence repeat (SSR) - enriched genomic libraries obtained from several di-, tri-, and tetranucleotide tandem repeat motifs. A total of 248 new polymorphic microsatellites were successfully optimized. The efficiency of the protocol applied (6.4%) was higher than that in other studies of fish that used the same method. Dinucleotide and perfect microsatellites were predominant in this species; the (AC)n motif was the most frequent class of repeat. Polymorphism and structural properties at these loci, together with 30 variable loci previously reported in turbot, were evaluated in 6 wild individuals. The number of alleles per locus ranged from 2 to 10, with an average of 4.046. The microsatellite markers characterized in this study will contribute to the development of the turbot genetic map, which can be used for quantitative trait locus (QTL) identification, marker-assisted selection programs, and other applications to improve its culture.

Mining and characterizing microsatellites from citrus ESTs

Freely available computer programs were arranged in a pipeline to extract microsatellites from public citrus EST sequences, retrieved from the NCBI. In total, 3,278 bi- to hexa-type SSR-containing sequences were identified from 56,199 citrus ESTs. On an average, one SSR was found per 5.2 kb of EST sequence, with the tri-nucleotide motifs as the most abundant. Primer sequences flanking SSR motifs were successfully identified from 2,295 citrus ESTs. Among those, a subset (100 pairs) were synthesized and tested to determine polymorphism and heterozygosity between/within two genera, sweet orange (C. sinensis) and Poncirus (P. trifoliata), which are the parents of the citrus core mapping population selected for an international citrus genomics effort. Eighty-seven pairs of primers gave PCR amplification to the anticipated SSRs, of which 52 and 35 appear to be homozygous and heterozygous, respectively, in sweet orange, and 67 and 20, respectively, in Poncirus. By pairing the loci between the two intergeneric species, it was found that 40 are heterozygous in at least one species with two alleles (9), three alleles (28), or four alleles (3), and the remaining 47 are homozygous in both species with either one allele (31) or two alleles (16). These EST-derived SSRs can be a resource used for understanding of the citrus SSR distribution and frequency, and development of citrus EST-SSR genetic and physical maps. These SSR primer sequences are available upon request.

Molecular analysis of the high stearic acid content in sunflower mutant CAS-14

Increasing the stearic acid content to improve sunflower (Helianthus annuus L.) oil quality is a desirable breeding objective for food-processing applications. CAS-14 is a sunflower mutant line with a high stearic acid content in its seed oil (>35% vs. <6% in currently grown sunflower hybrids), which is controlled by the Es3 gene. However, the expression of the high stearic acid character in CAS-14 is strongly influenced by temperature during seed maturation and it is not uniform along the seed. The objectives of this study were (1) to identify PCR-based molecular markers linked to the Es3 gene from CAS-14, (2) to map this gene on the sunflower genetic map, and (3) to characterize the interaction between CAS-14 and CAS-3, a sunflower high stearic acid (about 26%) mutant line with the Es1 and Es2 genes determining this trait. Two F2 mapping populations were developed from crosses between CAS-14 and P21, a nuclear male sterile line with the Ms11 gene controlling this character, and between CAS-14 and CAS-3. One hundred and thirty-three individuals from P21xCAS-14, and 164 individuals from CAS-3xCAS-14 were phenotyped in F2 and F3 seed generations for fatty acid composition using gas-liquid chromatography, and they were then genotyped with microsatellite [simple sequence repeat (SSR)] and insertion-deletion (INDEL) markers. Bulk segregant analysis in the P21xCAS-14 population identified two markers on LG 8 putatively linked to Es3. A large linkage group was identified using additional markers mapping to LG 8. Es3 mapped to the distal half of LG 8 and was flanked by the SSR markers ORS243 and ORS1161 at genetic distances of 0.5, and 3.9 cM, respectively. The Ms11 gene was also mapped to LG 8 and genetic distance between this gene and Es3 was found to be 7.4 cM. In the CAS-3xCAS-14 population, two QTLs were identified on LG 1 and LG 8, which underlie the Es1 gene from CAS-3 and the Es3 gene from CAS-14, respectively. A significant epistatic interaction between these two QTLs was found. Results from this study provided a basis for determining CAS-14 efficient breeding strategies.

Development, inheritance, and linkage-group assignment of 60 novel microsatellite markers for the gray, short-tailed opossum Monodelphis domestica

Short-tandem-repeat (SSR) or microsatellite polymorphisms are some of the most extensively employed genetic markers in contemporary linkage mapping studies. To date, only a limited number of microsatellites have been isolated in the gray, short-tailed opossum Monodelphis domestica, a South American marsupial widely used for comparative biological and biomedical research. To increase the number of potentially useful mapping markers, we screened 2 microsatellite-enriched genomic libraries containing alternatively (CA)n or (GA)n repeats. A total of 184 clones were sequenced, from which 60 polymorphic microsatellite markers were successfully optimized. The efficiency of this enrichment protocol for M. domestica microsatellite isolation is discussed, and suggestions to improve the outcome are made. All 60 loci showed high allelic diversity, with allele numbers ranging from 2 to 10 in a subset of 33 unrelated animals. Normal Mendelian inheritance was confirmed for all loci by analyzing allelic segregation in 5 two-generation families. One microsatellite appeared to be X linked, and null alleles were found in 5 others. Two-point linkage analyses were implemented using the data on the 5 families, leading to the assignment of 59 of these loci to the existing linkage groups. The 60 novel microsatellites developed in this study will contribute significantly to the M. domestica linkage map, and further QTL mapping studies.Key words: Monodelphis domestica, marsupial, microsatellite, enriched libraries, genetic linkage analysis.

Microsatellite marker polymorphism and mapping in pea (Pisum sativum L.)

This paper aims at providing reliable and cost effective genotyping conditions, level of polymorphism in a range of genotypes and map position of newly developed microsatellite markers in order to promote broad application of these markers as a common set for genetic studies in pea. Optimal PCR conditions were determined for 340 microsatellite markers based on amplification in eight genotypes. Levels of polymorphism were determined for 309 of these markers. Compared to data obtained for other species, levels of polymorphism detected in a panel of eight genotypes were high with a mean number of 3.8 alleles per polymorphic locus and an average PIC value of 0.62, indicating that pea represents a rather polymorphic autogamous species. One of our main objectives was to locate a maximum number of microsatellite markers on the pea genetic map. Data obtained from three different crosses were used to build a composite genetic map of 1,430 cM (Haldane) comprising 239 microsatellite markers. These include 216 anonymous SSRs developed from enriched genomic libraries and 13 SSRs located in genes. The markers are quite evenly distributed throughout the seven linkage groups of the map, with 85% of intervals between the adjacent SSR markers being smaller than 10 cM. There was a good conservation of marker order and linkage group assignment across the three populations. In conclusion, we hope this report will promote wide application of these markers and will allow information obtained by different laboratories worldwide in diverse fields of pea genetics, such as QTL mapping studies and genetic resource surveys, to be easily aligned.

Identification and validation of QTL for Sclerotinia midstalk rot resistance in sunflower by selective genotyping

Midstalk rot, caused by Sclerotinia sclerotiorum (Lib.) de Bary, is an important cause of yield loss in sunflower (Helianthus annuus L.). Objectives of this study were to: (1) estimate the number, genomic positions and genetic effects of quantitative trait loci (QTL) for resistance to midstalk rot in line TUB-5-3234, derived from an interspecific cross; (2) determine congruency of QTL between this line and other sources of resistance; and (3) make inferences about the efficiency of selective genotyping (SG) in detecting QTL conferring midstalk rot resistance in sunflower. Phenotypic data for three resistance (stem lesion, leaf lesion and speed of fungal growth) and two morphological (leaf length and leaf length with petiole) traits were obtained from 434 F3 families from cross CM625 (susceptible) x TUB-5-3234 (resistant) under artificial infection in field experiments across two environments. The SG was applied by choosing the 60 most resistant and the 60 most susceptible F3 families for stem lesion. For genotyping of the respective F2 plants, 78 simple sequence repeat markers were used. Genotypic variances were highly significant for all traits. Heritabilities and genotypic correlations between reMidstalk rot, caused by Sclerotinia sclerotiorum (Lib.) de Bary, is an important cause of yield loss in sunflower (Helianthus annuus L.). Objectives of this study were to: (1) estimate the number, genomic positions and genetic effects of quantitative trait loci (QTL) for resistance to midstalk rot in line TUB-5-3234, derived from an interspecific cross; (2) determine congruency of QTL between this line and other sources of resistance; and (3) make inferences about the efficiency of selective genotyping (SG) in detecting QTL conferring midstalk rot resistance in sunflower. Phenotypic data for three resistance (stem lesion, leaf lesion and speed of fungal growth) and two morphological (leaf length and leaf length with petiole) traits were obtained from 434 F3 families from cross CM625 (susceptible) x TUB-5-3234 (resistant) under artificial infection in field experiments across two environments. The SG was applied by choosing the 60 most resistant and the 60 most susceptible F3 families for stem lesion. For genotyping of the respective F2 plants, 78 simple sequence repeat markers were used. Genotypic variances were highly significant for all traits. Heritabilities and genotypic correlations between resistance traits were moderate to high. Three to four putative QTL were detected for each resistance trait explaining between 40.8% and 72.7% of the genotypic variance (PTS). Two QTL for stem lesion showed large genetic effects and corroborated earlier findings from the cross NDBLOSsel (resistant) x CM625 (susceptible). Our results suggest that SG can be efficiently used for QTL detection and the analysis of congruency for resistance genes across populations.

QTL mapping of Sclerotinia midstalk-rot resistance in sunflower

In many sunflower-growing regions of the world, Sclerotinia sclerotiorum (Lib.) de Bary is the major disease of sunflower (Helianthus annuus L.). In this study, we mapped and characterized quantitative trait loci (QTL) involved in resistance to S. sclerotiorum midstalk rot and two morphological traits. A total of 351 F3 families developed from a cross between a resistant inbred line from the germplasm pool NDBLOS and the susceptible line CM625 were assayed for their parental F2 genotype at 117 codominant simple sequence repeat markers. Disease resistance of the F3 families was screened under artificial infection in field experiments across two sowing times in 1999. For the three resistance traits (leaf lesion, stem lesion, and speed of fungal growth) and the two morphological traits, genotypic variances were highly significant. Heritabilities were moderate to high (h2=0.55-0.89). Genotypic correlations between resistance traits were highly significant (P<0.01) but moderate. QTL were detected for all three resistance traits, but estimated effects at most QTL were small. Simultaneously, they explained between 24.4% and 33.7% of the genotypic variance for resistance against S. sclerotiorum. Five of the 15 genomic regions carrying a QTL for either of the three resistance traits also carried a QTL for one of the two morphological traits. The prospects of marker-assisted selection (MAS) for resistance to S. sclerotiorum are limited due to the complex genetic architecture of the trait. MAS can be superior to classical phenotypic selection only with low marker costs and fast selection cycles.

Quantitative trait loci for broomrape (Orobanche cumana Wallr.) resistance in sunflower

Broomrape (Orobanche cumana Wallr.) is a root parasite of sunflower that is regarded as one of the most important constraints of sunflower production in the Mediterranean region. Breeding for resistance is the most effective method of control. P-96 is a sunflower line which shows dominant resistance to broomrape race E and recessive resistance to the very new race F. The objective of this study was to map and characterize quantitative trait loci (QTL) for resistance to race E and to race F of broomrape in P-96. A population from a cross between P-96 and the susceptible line P-21 was phenotyped for broomrape resistance in four experiments, two for race E and two for race F, by measuring different resistance parameters (resistance or susceptibility, number of broomrape per plant, and proportion of resistant plants per F(3) family). This population was also genotyped with microsatellite and RFLP markers. A linkage map comprising 103 marker loci distributed on 17 linkage groups was developed, and composite interval mapping analyses were performed. In total, five QTL ( or1.1, or3.1, or7.1 or13.1 and or13.2) for resistance to race E and six QTL ( or1.1, or4.1, or5.1, or13.1, or13.2 and or16.1) for resistance to race F of broomrape were detected on 7 of the 17 linkage groups. Phenotypic variance for race E resistance was mainly explained by the major QTL or3.1 associated to the resistance or susceptibility character ( R(2)=59%), while race F resistance was explained by QTL with a small to moderate effect ( R(2) from 15.0% to 38.7%), mainly associated with the number of broomrape per plant. Or3.1 was race E-specific, while or1.1, or13.1 and or13.2 of were non-race specific. Or13.1, and or13.2 were stable across the four experiments. Or3.1, and or7.1 were stable over the two race E experiments and or1.1 and or5.1 over the two race F experiments. The results from this study suggest that resistance to broomrape in sunflower is controlled by a combination of qualitative, race-specific resistance affecting the presence or absence of broomrape and a quantitative non-race specific resistance affecting their number.

Isolation and characterization of 70 novel microsatellite markers from ostrich (Struthiocamelus) genome

Microsatellite markers are widely used in linkage mapping, parentage testing, population genetic studies, and molecular evolution studies in many agricultural species, while only a limited number of ostrich (Struthio camelus) microsatellites have been isolated. Thus, we constructed a random small-insert genomic library and a microsatellite-enriched library containing CA repeats. Fourteen clones containing CA repeats were isolated from 3462 clones in the non-enriched library by radioactive screening and 248 positive clones were isolated from 300 sequenced clones from the enriched library by PCR screening. After the enrichment procedures, the proportion of clones containing CA repeats was raised to 78.8%, compared with 0.4% in the non-enriched libraries, indicating that the enrichment value approaches 200 fold, which decreased the time and cost of cloning. The number of complete simple CA repeats in these positive clones ranged from 5 to 29. The primers for 94 of these microsatellites were developed and used to detect polymorphisms, of which 61 loci exhibited length polymorphisms in 17 unrelated ostrich individuals. The new polymorphic microsatellite markers we have identified and characterized will contribute to the ostrich genetic map, parentage testing, and comparative genomics between avian species.Key words: ostrich, microsatellite markers, enriched library, polymorphism.

Differential representation of sunflower ESTs in enriched organ-specific cDNA libraries in a small scale sequencing project

Background

Subtractive hybridization methods are valuable tools for identifying differentially regulated genes in a given tissue avoiding redundant sequencing of clones representing the same expressed genes, maximizing detection of low abundant transcripts and thus, affecting the efficiency and cost effectiveness of small scale cDNA sequencing projects aimed to the specific identification of useful genes for breeding purposes. The objective of this work is to evaluate alternative strategies to high-throughput sequencing projects for the identification of novel genes differentially expressed in sunflower as a source of organ-specific genetic markers that can be functionally associated to important traits.

Results

Differential organ-specific ESTs were generated from leaf, stem, root and flower bud at two developmental stages (R1 and R4). The use of different sources of RNA as tester and driver cDNA for the construction of differential libraries was evaluated as a tool for detection of rare or low abundant transcripts. Organ-specificity ranged from 75 to 100% of non-redundant sequences in the different cDNA libraries. Sequence redundancy varied according to the target and driver cDNA used in each case. The R4 flower cDNA library was the less redundant library with 62% of unique sequences. Out of a total of 919 sequences that were edited and annotated, 318 were non-redundant sequences. Comparison against sequences in public databases showed that 60% of non-redundant sequences showed significant similarity to known sequences. The number of predicted novel genes varied among the different cDNA libraries, ranging from 56% in the R4 flower to 16 % in the R1 flower bud library. Comparison with sunflower ESTs on public databases showed that 197 of non-redundant sequences (60%) did not exhibit significant similarity to previously reported sunflower ESTs. This approach helped to successfully isolate a significant number of new reported sequences putatively related to responses to important agronomic traits and key regulatory and physiological genes.

Conclusions

The application of suppressed subtracted hybridization technology not only enabled the cost effective isolation of differentially expressed sequences but it also allowed the identification of novel sequences in sunflower from a relative small number of analyzed sequences when compared to major sequencing projects.