SNP2CAPS: a SNP and INDEL analysis tool for CAPS marker development

Development and Applications of a High Throughput Genotyping Tool for Polyploid Crops: Single Nucleotide Polymorphism (SNP) Array

Polypoid species play significant roles in agriculture and food production. Many crop species are polyploid, such as potato, wheat, strawberry, and sugarcane. Genotyping has been a daunting task for genetic studies of polyploid crops, which lags far behind the diploid crop species. Single nucleotide polymorphism (SNP) array is considered to be one of, high-throughput, relatively cost-efficient and automated genotyping approaches. However, there are significant challenges for SNP identification in complex, polyploid genomes, which has seriously slowed SNP discovery and array development in polyploid species. Ploidy is a significant factor impacting SNP qualities and validation rates of SNP markers in SNP arrays, which has been proven to be a very important tool for genetic studies and molecular breeding. In this review, we (1) discussed the pros and cons of SNP array in general for high throughput genotyping, (2) presented the challenges of and solutions to SNP calling in polyploid species, (3) summarized the SNP selection criteria and considerations of SNP array design for polyploid species, (4) illustrated SNP array applications in several different polyploid crop species, then (5) discussed challenges, available software, and their accuracy comparisons for genotype calling based on SNP array data in polyploids, and finally (6) provided a series of SNP array design and genotype calling recommendations. This review presents a complete overview of SNP array development and applications in polypoid crops, which will benefit the research in molecular breeding and genetics of crops with complex genomes.

The genome sequence of sweet cherry (Prunus avium) for use in genomics-assisted breeding

We determined the genome sequence of sweet cherry (Prunus avium) using next-generation sequencing technology. The total length of the assembled sequences was 272.4 Mb, consisting of 10,148 scaffold sequences with an N50 length of 219.6 kb. The sequences covered 77.8% of the 352.9 Mb sweet cherry genome, as estimated by k-mer analysis, and included >96.0% of the core eukaryotic genes. We predicted 43,349 complete and partial protein-encoding genes. A high-density consensus map with 2,382 loci was constructed using double-digest restriction site–associated DNA sequencing. Comparing the genetic maps of sweet cherry and peach revealed high synteny between the two genomes; thus the scaffolds were integrated into pseudomolecules using map- and synteny-based strategies. Whole-genome resequencing of six modern cultivars found 1,016,866 SNPs and 162,402 insertions/deletions, out of which 0.7% were deleterious. The sequence variants, as well as simple sequence repeats, can be used as DNA markers. The genomic information helps us to identify agronomically important genes and will accelerate genetic studies and breeding programs for sweet cherries. Further information on the genomic sequences and DNA markers is available in DBcherry (http://cherry.kazusa.or.jp (8 May 2017, date last accessed)).

Candidate genes and molecular markers associated with heat tolerance in colonial Bentgrass

Elevated temperature is a major abiotic stress limiting the growth of cool-season grasses during the summer months. The objectives of this study were to determine the genetic variation in the expression patterns of selected genes involved in several major metabolic pathways regulating heat tolerance for two genotypes contrasting in heat tolerance to confirm their status as potential candidate genes, and to identify PCR-based markers associated with candidate genes related to heat tolerance in a colonial (Agrostis capillaris L.) x creeping bentgrass (Agrostis stolonifera L.) hybrid backcross population. Plants were subjected to heat stress in controlled-environmental growth chambers for phenotypic evaluation and determination of genetic variation in candidate gene expression. Molecular markers were developed for genes involved in protein degradation (cysteine protease), antioxidant defense (catalase and glutathione-S-transferase), energy metabolism (glyceraldehyde-3-phosphate dehydrogenase), cell expansion (expansin), and stress protection (heat shock proteins HSP26, HSP70, and HSP101). Kruskal-Wallis analysis, a commonly used non-parametric test used to compare population individuals with or without the gene marker, found the physiological traits of chlorophyll content, electrolyte leakage, normalized difference vegetative index, and turf quality were associated with all candidate gene markers with the exception of HSP101. Differential gene expression was frequently found for the tested candidate genes. The development of candidate gene markers for important heat tolerance genes may allow for the development of new cultivars with increased abiotic stress tolerance using marker-assisted selection.

PCR-based assays for validation of single nucleotide polymorphism markers in rice and mungbean

BACKGROUND

Single nucleotide polymorphism (SNP) markers are the method of choice for genetic analyses including diversity and quantitative trait loci (QTL) studies. Marker validation is essential for QTL studies, but the cost and workload are considerable when large numbers of markers need to be verified. Marker systems with low development costs would be most suitable for this task.

RESULTS

We have tested allele specific polymerase chain reaction (PCR), tetra markers and a genotyping tool based on the single strand specific nuclease CEL-I to verify randomly selected SNP markers identified previously either with a SNP array or by genotyping by sequencing in rice and mungbean, respectively. The genotyping capacity of allele-specific PCR and tetra markers was affected by the sequence context surrounding the SNP; SNPs located in repeated sequences and in GC-rich stretches could not be correctly identified. In contrast, CEL-I digestion of mixed fragments produced from test and reference DNA reliably pinpointed the correct genotypes, yet scoring of the genotypes became complicated when multiple SNPs were present in the PCR fragments. A cost analysis showed that as long the sample number remains small, CEL-I genotyping is more cost-effective than tetra markers.

CONCLUSIONS

CEL-I genotyping performed better in terms of genotyping accuracy and costs than tetra markers. The method is highly useful for validating SNPs in small to medium size germplasm panels.

Molecular Mapping of Reduced Plant Height Gene Rht24 in Bread Wheat

Height is an important trait related to plant architecture and yield potential in bread wheat (Triticum aestivum L.). We previously identified a major quantitative trait locus QPH.caas-6A flanked by simple sequence repeat markers Xbarc103 and Xwmc256 that reduced height by 8.0-10.4%. Here QPH.caas-6A, designated as Rht24, was confirmed using recombinant inbred lines (RILs) derived from a Jingdong 8/Aikang 58 cross. The target sequences of Xbarc103 and Xwmc256 were used as queries to BLAST against International Wheat Genome Sequence Consortium database and hit a super scaffold of approximately 208 Mb. Based on gene annotation of the scaffold, three gene-specific markers were developed to genotype the RILs, and Rht24 was narrowed to a 1.85 cM interval between TaAP2 and TaFAR. In addition, three single nucleotide polymorphism (SNP) markers linked to Rht24 were identified from SNP chip-based screening in combination with bulked segregant analysis. The allelic efficacy of Rht24 was validated in 242 elite wheat varieties using TaAP2 and TaFAR markers. These showed a significant association between genotypes and plant height. Rht24 reduced plant height by an average of 6.0-7.9 cm across environments and were significantly associated with an increased TGW of 2.0-3.4 g. The findings indicate that Rht24 is a common dwarfing gene in wheat breeding, and TaAP2 and TaFAR can be used for marker-assisted selection.

Single nucleotide polymorphism in sugar pathway and disease resistance genes in sugarcane

Single nucleotide polymorphism in sugar pathway and disease resistance genes showing genetic association with sugar content and red rot resistance would be useful in marker-assisted genetic improvement of sugarcane. Validation and genotyping of potential sequence variants in candidate genes are necessary to understand their functional significance and trait association potential. We discovered, characterized, validated and genotyped SNPs and InDels in sugar pathway and disease resistance genes of Saccharum complex and sugarcane varieties using amplicon sequencing and CAPS assays. The SNPs were abundant in the non-coding 3'UTRs than 5'UTRs and coding sequences depicting a strong bias toward C to T transition substitutions than transversions. Sequencing of cloned amplicons validated 61.6 and 45.2 % SNPs detected in silico in 21 sugar pathway and 16 disease resistance genes, respectively. Sixteen SNPs in four sugar pathway genes and 10 SNPs in nine disease resistance genes were validated through cost-effective CAPS assay. Functional and adaptive significance of SNP and protein haplotypes identified in sugar pathway and disease resistance genes was assessed by correlating their allelic variation with missense amino acid substitutions in the functional domains, alteration in protein structure models and possible modulation of catalytic enzyme activity in contrasting high and low sugar and moderately red rot resistant and highly susceptible sugarcane genotypes. A strong genetic association of five SNPs in the sugar pathway and disease resistance genes, and an InDel marker in the promoter sequence of sucrose synthase-2 gene, with sugar content and red rot resistance, was evident. The functionally relevant SNPs and InDels, detected and validated in sugar pathway and disease resistance genes, and genic CAPS markers designed, would be of immense use in marker-assisted genetic improvement of sugarcane for sugar content and disease resistance.

Disruption of OsSULTR3;3 reduces phytate and phosphorus concentrations and alters the metabolite profile in rice grains

Two low phytic acid (lpa) mutants have been developed previously with the aim to improve the nutritional value of rice (Oryza sativa) grains. In the present study, the impacts of lpa mutations on grain composition and underlying molecular mechanisms were investigated. Comparative compositional analyses and metabolite profiling demonstrated that concentrations of both phytic acid (PA) and total phosphorus (P) were significantly reduced in lpa brown rice, accompanied by changes in other metabolites and increased concentrations of nutritionally relevant compounds. The lpa mutations modified the expression of a number of genes involved in PA metabolism, as well as in sulfate and phosphate homeostasis and metabolism. Map-based cloning and complementation identified the underlying lpa gene to be OsSULTR3;3. The promoter of OsSULTR3;3 is highly active in the vascular bundles of leaves, stems and seeds, and its protein is localized in the endoplasmic reticulum. No activity of OsSULTR3;3 was revealed for the transport of phosphate, sulfate, inositol or inositol 1,4,5 triphosphate by heterologous expression in either yeast or Xenopus oocytes. The findings reveal that OsSULTR3;3 plays an important role in grain metabolism, pointing to a new route to generate value-added grains in rice and other cereal crops.

CisSERS: Customizable In Silico Sequence Evaluation for Restriction Sites

High-throughput sequencing continues to produce an immense volume of information that is processed and assembled into mature sequence data. Data analysis tools are urgently needed that leverage the embedded DNA sequence polymorphisms and consequent changes to restriction sites or sequence motifs in a high-throughput manner to enable biological experimentation. CisSERS was developed as a standalone open source tool to analyze sequence datasets and provide biologists with individual or comparative genome organization information in terms of presence and frequency of patterns or motifs such as restriction enzymes. Predicted agarose gel visualization of the custom analyses results was also integrated to enhance the usefulness of the software. CisSERS offers several novel functionalities, such as handling of large and multiple datasets in parallel, multiple restriction enzyme site detection and custom motif detection features, which are seamlessly integrated with real time agarose gel visualization. Using a simple fasta-formatted file as input, CisSERS utilizes the REBASE enzyme database. Results from CisSERS enable the user to make decisions for designing genotyping by sequencing experiments, reduced representation sequencing, 3’UTR sequencing, and cleaved amplified polymorphic sequence (CAPS) molecular markers for large sample sets. CisSERS is a java based graphical user interface built around a perl backbone. Several of the applications of CisSERS including CAPS molecular marker development were successfully validated using wet-lab experimentation. Here, we present the tool CisSERS and results from in-silico and corresponding wet-lab analyses demonstrating that CisSERS is a technology platform solution that facilitates efficient data utilization in genomics and genetics studies.

Development of cleaved amplified polymorphic sequence markers and a CAPS-based genetic linkage map in watermelon (Citrullus lanatus [Thunb.] Matsum. and Nakai) constructed using whole-genome re-sequencing data

Cleaved amplified polymorphic sequence (CAPS) markers are useful tools for detecting single nucleotide polymorphisms (SNPs). This study detected and converted SNP sites into CAPS markers based on high-throughput re-sequencing data in watermelon, for linkage map construction and quantitative trait locus (QTL) analysis. Two inbred lines, Cream of Saskatchewan (COS) and LSW-177 had been re-sequenced and analyzed by Perl self-compiled script for CAPS marker development. 88.7% and 78.5% of the assembled sequences of the two parental materials could map to the reference watermelon genome, respectively. Comparative assembled genome data analysis provided 225,693 and 19,268 SNPs and indels between the two materials. 532 pairs of CAPS markers were designed with 16 restriction enzymes, among which 271 pairs of primers gave distinct bands of the expected length and polymorphic bands, via PCR and enzyme digestion, with a polymorphic rate of 50.94%. Using the new CAPS markers, an initial CAPS-based genetic linkage map was constructed with the F2 population, spanning 1836.51 cM with 11 linkage groups and 301 markers. 12 QTLs were detected related to fruit flesh color, length, width, shape index, and brix content. These newly CAPS markers will be a valuable resource for breeding programs and genetic studies of watermelon.

Development of genome-wide InDel markers and their integration with SSR, DArT and SNP markers in single barley map

Background

Development of molecular markers such as SSR (simple sequence repeat), DArT (diversity arrays technology) and SNP (single nucleotide polymorphism) is fundamental for linkage map construction and QTL mapping. However, DArT and SNP genotyping require special tools, and detection of SSR polymorphisms requires time-consuming polyacrylamide electrophoresis. Furthermore, many markers have been mapped in different populations such that their genetic positions are inconsistent. Recently, InDel (insertion and deletion) markers have become popular in genetic map construction and map-based cloning.

Results

Aligning genomic DNA sequences in two barley cultivars (Morex and Barke) identified 436,640 InDels. We designed 1140 InDel markers across the barley genome with an average genetic distance of 1 cM, each having a unique location in the barley genome. High-resolution melting (HRM) technology was used to genotype 55 InDel markers; those PCR amplicons with melting temperature differences >0.3 °C were ideal for HRM genotyping. The 1140 InDel markers together with 383 SSRs, 3909 gene-based SNPs and 1544 DArT markers were integrated into single barley genetic map according to their physical map positions.

Conclusions

High-density InDel markers with specific genome locations were developed, with 6976 molecular markers (SSRs, DArTs, SNPs and InDels) integrated into single barley genetic map. HRM genotyping of the InDel markers each with single PCR band will facilitate quick map construction and gene fine-mapping.

Electronic supplementary material

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

PLANT VOLATILES. Biosynthesis of monoterpene scent compounds in roses

The scent of roses (Rosa x hybrida) is composed of hundreds of volatile molecules. Monoterpenes represent up to 70% percent of the scent content in some cultivars, such as the Papa Meilland rose. Monoterpene biosynthesis in plants relies on plastid-localized terpene synthases. Combining transcriptomic and genetic approaches, we show that the Nudix hydrolase RhNUDX1, localized in the cytoplasm, is part of a pathway for the biosynthesis of free monoterpene alcohols that contribute to fragrance in roses. The RhNUDX1 protein shows geranyl diphosphate diphosphohydrolase activity in vitro and supports geraniol biosynthesis in planta.

Fine mapping of shattering locus Br2 reveals a putative chromosomal inversion polymorphism between the two lineages of Aegilops tauschii

This work laid the foundation for cloning of shattering gene Br2 and provided first line of evidence that two major Aegilops tauschii lineages are differentiated by an inversion polymorphism. Chromosome inversions often accompany population differentiation and capture local adaptation during speciation. Aegilops tauschii, the D-genome donor species of hexaploid wheat, consists of two genetically isolated lineages, L1 and L2, but little is known about the genetic mechanisms underlying the population differentiation in this diploid species. During fine mapping of the shattering gene Br2 using a large F2 population derived from a cross between TA1604 (an L1 accession) and AL8/78 (an L2 accession), we found contrasting patterns of crossover distribution in the Br2 interval and neighboring regions despite the high local gene synteny with Brachypodium distachyon and rice. Br2 was localized in a 0.08-cM interval, and 13 marker loci formed a block, where single-crossovers were completely suppressed, but double-crossovers were enriched with a recombination rate of ~11 cM/Mb. In contrast, in a neighboring region no double-crossover was recovered, but single-crossover rate reached 24 cM/Mb, which is much higher than the genome-wide average. This result suggests a putative inversion polymorphism between the parental lines in the Br2 region. Genotyping using the markers from the Br2 region divided a collection of 55 randomly sampled A. tauschii accessions into two major groups, and they are largely genetically isolated. The two groups correspond to the L1 and L2 lineages based on their geographic distribution patterns. This provides first evidence that inversions may underlie the evolution of A. tauschii lineages. The presence of inter-lineage inversions may complicate map-based cloning in A. tauschii and transfer of useful traits to wheat.

Development of STS and CAPS markers for variety identification and genetic diversity analysis of tea germplasm in Taiwan

BACKGROUND

Tea (Camellia sinensis) is an important economic crop in Taiwan. Particularly, two major commercial types of tea (Paochong tea and Oolong tea) which are produced in Taiwan are famous around the world, and they must be manufactured with specific cultivars. Nevertheless, many elite cultivars have been illegally introduced to foreign countries. Because of the lower cost, large amount of "Taiwan-type tea" are produced and imported to Taiwan, causing a dramatic damage in the tea industry. It is very urgent to develop the stable, fast and reliable DNA markers for fingerprinting tea cultivars in Taiwan and protecting intellectual property rights for breeders. Furthermore, genetic diversity and phylogenetic relationship evaluations of tea germplasm in Taiwan are imperative for parental selection in the cross-breeding program and avoidance of genetic vulnerability.

RESULTS

Two STS and 37 CAPS markers derived from cytoplasmic genome and ESTs of tea have been developed in this study providing a useful tool for distinguishing all investigated germplasm. For identifying 12 prevailing tea cultivars in Taiwan, five core markers, including each one of mitochondria and chloroplast, and three nuclear markers, were developed. Based on principal coordinate analysis and cluster analysis, 55 tea germplasm in Taiwan were divided into three groups: sinensis type (C. sinensis var. sinensis), assamica type (C. sinensis var. assamica) and Taiwan wild species (C. formosensis). The result of genetic diversity analysis revealed that both sinensis (0.44) and assamica (0.41) types had higher genetic diversity than wild species (0.25). The close genetic distance between the first (Chin-Shin-Oolong) and the third (Shy-Jih-Chuen) prevailing cultivars was found, and many recently released varieties are the descents of Chin-Shin-Oolong. This implies the potential risk of genetic vulnerability for tea cultivation in Taiwan.

CONCLUSIONS

We have successfully developed a tool for tea germplasm discrimination and genetic diversity analysis, as well as a set of core markers for effective identification of prevailing cultivars in Taiwan. According to the results of phylogenetic analysis on prevailing tea cultivars, it is necessary to broaden genetic diversity from wild species or plant introduction in future breeding programs.

Analysis of bymovirus resistance genes on proximal barley chromosome 4HL provides the basis for precision breeding for BaMMV/BaYMV resistance

Unlocking allelic diversity of the bymovirus resistance gene rym11 located on proximal barley chromosome 4HL and diagnostic markers provides the basis for precision breeding for BaMMV/BaYMV resistance. The recessive resistance gene rym11 on barley chromosome 4HL confers broad-spectrum and complete resistance to all virulent European isolates of Barley mild mosaic virus and Barley yellow mosaic virus (BaMMV/BaYMV). As previously reported, rym11-based resistance is conferred by a series of alleles of naturally occurring deletions in the gene HvPDIL5-1, encoding a protein disulfide isomerase-like protein. Here, a novel resistance-conferring allele of rym11 is reported that, in contrast to previously identified resistance-conferring variants of the gene HvPDIL5-1, carries a single non-synonymous amino acid substitution. Allelism was confirmed by crossing to genotypes carrying previously known rym11 alleles. Crossing rym11 genotypes with a cultivar carrying the recessive resistance gene rym1, which was reported to reside on the same chromosome arm 4HL like rym11, revealed allelism of both loci. This allelic state was confirmed by re-sequencing HvPDIL5-1 in the rym1 genotype, detecting the haplotype of the rym11-d allele. Diagnostic PCR-based markers were established to differentiate all seven resistance-conferring alleles of the rym11 locus providing precise tools for marker-assisted selection (MAS) of rym11 in barley breeding.

Mapping-by-sequencing accelerates forward genetics in barley

Mapping-by-sequencing has emerged as a powerful technique for genetic mapping in several plant and animal species. As this resequencing-based method requires a reference genome, its application to complex plant genomes with incomplete and fragmented sequence resources remains challenging. We perform exome sequencing of phenotypic bulks of a mapping population of barley segregating for a mutant phenotype that increases the rate of leaf initiation. Read depth analysis identifies a candidate gene, which is confirmed by the analysis of independent mutant alleles. Our method illustrates how the genomic resources of barley together with exome resequencing can underpin mapping-by-sequencing.

Mapping-by-sequencing accelerates forward genetics in barley

Mapping-by-sequencing has emerged as a powerful technique for genetic mapping in several plant and animal species. As this resequencing-based method requires a reference genome, its application to complex plant genomes with incomplete and fragmented sequence resources remains challenging. We perform exome sequencing of phenotypic bulks of a mapping population of barley segregating for a mutant phenotype that increases the rate of leaf initiation. Read depth analysis identifies a candidate gene, which is confirmed by the analysis of independent mutant alleles. Our method illustrates how the genomic resources of barley together with exome resequencing can underpin mapping-by-sequencing.

Development of SNP markers for genes of the phenylpropanoid pathway and their association to kernel and malting traits in barley

BACKGROUND

Flavonoids are an important class of secondary compounds in angiosperms. Next to certain biological functions in plants, they play a role in the brewing process and have an effect on taste, color and aroma of beer. The aim of this study was to reveal the haplotype diversity of candidate genes involved in the phenylpropanoid biosynthesis pathway in cultivated barley varieties (Hordeum vulgare L.) and to determine associations to kernel and malting quality parameters.

RESULTS

Five genes encoding phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), chalcone synthase (CHS), flavanone 3-hydroxylase (F3H) and dihydroflavonol reductase (DFR) of the phenylpropanoid biosynthesis pathway were partially resequenced in 16 diverse barley reference genotypes. Their localization in the barley genome, their genetic structure, and their genetic variation e.g. single nucleotide polymorphism (SNP) and Insertion/Deletion (InDel) patterns were revealed. In total, 130 SNPs and seven InDels were detected. Of these, 21 polymorphisms were converted into high-throughput pyrosequencing markers. The resulting SNP and haplotype patterns were used to calculate associations with kernel and malting quality parameters.

CONCLUSIONS

SNP patterns were found to be highly variable for the investigated genes. The developed high-throughput markers are applicable for assessing the genetic variability and for the determination of haplotype patterns in a set of barley accessions. The candidate genes PAL, C4H and F3H were shown to be associated to several malting properties like glassiness (PAL), viscosity (C4H) or to final attenuation (F3H).

BSTA: a targeted approach combines bulked segregant analysis with next- generation sequencing and de novo transcriptome assembly for SNP discovery in sunflower

Background

Sunflower belongs to the largest plant family on earth, the genomically poorly explored Compositae. Downy mildew Plasmopara halstedii (Farlow) Berlese & de Toni is one of the major diseases of cultivated sunflower (Helianthus annuus L.). In the search for new sources of downy mildew resistance, the locus Pl_ARG on linkage group 1 (LG1) originating from H. argophyllus is promising since it confers resistance against all known races of the pathogen. However, the mapping resolution in the Pl_ARG region is hampered by significantly suppressed recombination and by limited availability of polymorphic markers. Here we examined a strategy developed for the enrichment of molecular markers linked to this specific genomic region. We combined bulked segregant analysis (BSA) with next-generation sequencing (NGS) and de novo assembly of the sunflower transcriptome for single nucleotide polymorphism (SNP) discovery in a sequence resource combining reads originating from two sunflower species, H. annuus and H. argophyllus.

Results

A computational pipeline developed for SNP calling and pattern detection identified 219 candidate genes. For a proof of concept, 42 resistance gene-like sequences were subjected to experimental SNP validation. Using a high-resolution mapping population, 12 SNP markers were mapped to LG1. We successfully verified candidate sequences either co-segregating with or closely flanking Pl_ARG.

Conclusions

This study is the first successful example to improve bulked segregant analysis with de novo transcriptome assembly using next generation sequencing. The BSTA pipeline we developed provides a useful guide for similar studies in other non-model organisms. Our results demonstrate this method is an efficient way to enrich molecular markers and to identify candidate genes in a specific mapping interval.

Genome-wide association studies using single nucleotide polymorphism markers developed by re-sequencing of the genomes of cultivated tomato

With the aim of understanding relationship between genetic and phenotypic variations in cultivated tomato, single nucleotide polymorphism (SNP) markers covering the whole genome of cultivated tomato were developed and genome-wide association studies (GWAS) were performed. The whole genomes of six tomato lines were sequenced with the ABI-5500xl SOLiD sequencer. Sequence reads covering ∼13.7× of the genome for each line were obtained, and mapped onto tomato reference genomes (SL2.40) to detect ∼1.5 million SNP candidates. Of the identified SNPs, 1.5% were considered to confer gene functions. In the subsequent Illumina GoldenGate assay for 1536 SNPs, 1293 SNPs were successfully genotyped, and 1248 showed polymorphisms among 663 tomato accessions. The whole-genome linkage disequilibrium (LD) analysis detected highly biased LD decays between euchromatic (58 kb) and heterochromatic regions (13.8 Mb). Subsequent GWAS identified SNPs that were significantly associated with agronomical traits, with SNP loci located near genes that were previously reported as candidates for these traits. This study demonstrates that attractive loci can be identified by performing GWAS with a large number of SNPs obtained from re-sequencing analysis.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

Conserved synteny-based anchoring of the barley genome physical map

Gene order is largely collinear in the small-grained cereals, a feature which has proved helpful in both marker development and positional cloning. The accuracy of a virtual gene order map ("genome zipper") for barley (Hordeum vulgare), developed by combining a genetic map of this species with a large number of gene locations obtained from the maps constructed in other grass species, was evaluated here both at the genome-wide level and at the fine scale in a representative segment of the genome. Comparing the whole genome "genome zipper" maps with a genetic map developed by using transcript-derived markers, yielded an accuracy of >94 %. The fine-scale comparison involved a 14 cM segment of chromosome arm 2HL. One hundred twenty-eight genes of the "genome zipper" interval were analysed. Over 95 % (45/47) of the polymorphic markers were genetically mapped and allocated to the expected region of 2HL, following the predicted order. A further 80 of the 128 genes were assigned to the correct chromosome arm 2HL by analysis of wheat-barley addition lines. All 128 gene-based markers developed were used to probe a barley bacterial artificial chromosome (BAC) library, delivering 26 BAC contigs from which all except two were anchored to the targeted zipper interval. The results demonstrate that the gene order predicted by the "genome zipper" is remarkably accurate and that the "genome zipper" represents a highly efficient informational resource for the systematic identification of gene-based markers and subsequent physical map anchoring of the barley genome.

Using the candidate gene approach for detecting genes underlying seed oil concentration and yield in soybean

Increasing the oil concentration in soybean seeds has been given more attention in recent years because of demand for both edible oil and biodiesel production. Oil concentration in soybean is a complex quantitative trait regulated by many genes as well as environmental conditions. To identify genes governing seed oil concentration in soybean, 16 putative candidate genes of three important gene families (GPAT: acyl-CoA:sn-glycerol-3-phosphate acyltransferase, DGAT: acyl-CoA:diacylglycerol acyltransferase, and PDAT: phospholipid:diacylglycerol acyltransferase) involved in triacylglycerol (TAG) biosynthesis pathways were selected and their sequences retrieved from the soybean database ( http://www.phytozome.net/soybean ). Three sequence mutations were discovered in either coding or noncoding regions of three DGAT soybean isoforms when comparing the parents of a 203 recombinant inbreed line (RIL) population; OAC Wallace and OAC Glencoe. The RIL population was used to study the effects of these mutations on seed oil concentration and other important agronomic and seed composition traits, including seed yield and protein concentration across three field locations in Ontario, Canada, in 2009 and 2010. An insertion/deletion (indel) mutation in the GmDGAT2B gene in OAC Wallace was significantly associated with reduced seed oil concentration across three environments and reduced seed yield at Woodstock in 2010. A mutation in the 3' untranslated (3'UTR) region of GmDGAT2C was associated with seed yield at Woodstock in 2009. A mutation in the intronic region of GmDGAR1B was associated with seed yield and protein concentration at Ottawa in 2010. The genes identified in this study had minor effects on either seed yield or oil concentration, which was in agreement with the quantitative nature of the traits. However, the novel gene-specific markers designed in the present study can be used in soybean breeding for marker-assisted selection aimed at increasing seed yield and oil concentration with no significant impact on seed protein concentration.

Genetic analysis of tolerance to boron toxicity in the legume Medicago truncatula

BACKGROUND

Medicago truncatula Gaertn. (barrel medic) is cultivated as a pasture legume for its high protein content and ability to improve soils through nitrogen fixation. Toxic concentrations of the micronutrient Boron (B) in agricultural soils hamper the production of cereal and leguminous crops. In cereals, the genetic analysis of B tolerance has led to the development of molecular selection tools to introgress and maintain the B tolerance trait in breeding lines. There is a comparable need for selection tools in legumes that grow on these toxic soils, often in rotation with cereals.

RESULTS

Genetic variation for B tolerance in Medicago truncatula was utilised to generate two F2 populations from crosses between tolerant and intolerant parents. Phenotyping under B stress revealed a close correlation between B tolerance and biomass production and a segregation ratio explained by a single dominant locus. M. truncatula homologues of the Arabidopsis major intrinsic protein (MIP) gene AtNIP5;1 and the efflux-type transporter gene AtBOR1, both known for B transport, were identified and nearby molecular markers screened across F2 lines to verify linkage with the B-tolerant phenotype. Most (95%) of the phenotypic variation could be explained by the SSR markers h2_6e22a and h2_21b19a, which flank a cluster of five predicted MIP genes on chromosome 4. Three CAPS markers (MtBtol-1,-2,-3) were developed to dissect the region further. Expression analysis of the five predicted MIPs indicated that only MtNIP3 was expressed when leaf tissue and roots were assessed. MtNIP3 showed low and equal expression in the roots of tolerant and intolerant lines but a 4-fold higher expression level in the leaves of B-tolerant cultivars. The expression profile correlates closely with the B concentration measured in the leaves and roots of tolerant and intolerant plants. Whereas no significant difference in B concentration exists between roots of tolerant and intolerant plants, the B concentration in the leaves of tolerant plants is less than half that of intolerant plants, which further supports MtNIP3 as the best candidate for the tolerance trait-defining gene in Medicago truncatula.

CONCLUSION

The close linkage of the MtNIP3 locus to B toxicity tolerance provides a source of molecular selection tools to pasture breeding programs. The economical importance of the locus warrants further investigation of the individual members of the MIP gene cluster in other pasture and in grain legumes.

RAD-Seq derived markers flank the shell colour and banding loci of the Cepaea nemoralis supergene

Studies on the classic shell colour and banding polymorphism of the land snail Cepaea played a crucial role in establishing the importance of natural selection in maintaining morphological variation. Cepaea is also a pre-eminent model for ecological genetics because the outward colour and banding phenotype is entirely genetically determined, primarily by a 'supergene' of at least five loci. Unfortunately, progress in understanding the evolution and maintenance of the Cepaea polymorphism stalled, partly because of a lack of genetic markers. With a view to re-establish Cepaea as a prominent model of molecular ecology, we made six laboratory crosses of Cepaea nemoralis, five of which segregated for shell ground colour (C) and the presence or absence of bands (B). First, scoring of colour and banding in 323 individuals found no recombination between the C and B loci of the supergene. Second, using restriction site-associated DNA sequencing (RAD-Seq) of two parents and 22 offspring, we identified 44 anonymous markers putatively linked to the colour (C) and banding (B) loci. The genotype of eleven of the most promising RAD-Seq markers was independently validated in the same 22 offspring, then up to a further 146 offspring were genotyped. The closest RAD-Seq markers scored are within ~0.6 centimorgan (cM) of the C-B supergene linkage group, with the combined loci together forming a 35.8 cM linkage map of markers that flank both sides of the Cepaea C-B supergene.

Asymmetric functional divergence of young, dispersed gene duplicates in Arabidopsis thaliana

One prediction of the classic Ohno model of gene duplication predicts that new genes form from the asymmetric functional divergence of a newly arisen, redundant duplicate locus. In order to understand the mechanisms which give rise to functional divergence of newly formed dispersed duplicates, we assessed the expression and molecular evolutionary divergence of a suite of 19 highly similar dispersed duplicates in Arabidopsis thaliana. These duplicates have a K sil equal to or less than 5 % and are specific to the A. thaliana lineage; thus, they predictably represent some of the youngest duplicates in the A. thaliana genome. We found that the majority of young duplicate loci exhibit asymmetric expression patterns, with the daughter locus exhibiting reduced expression across all tissues analyzed relative to the progenitor locus or simply not expressed. Furthermore, daughter loci, on the whole, have significantly more nonsynonymous substitutions than the progenitor loci. We also identified four pairs of loci which exhibit significant (P < 0.05) evolutionary rate asymmetry, three of which exhibit elevated dN/dS in the duplicate copy. We suggest, based on these data, that functional diversification initially takes the form of asymmetric regulatory divergence that can be a direct consequence of the mode of duplication. The reduced and/or absence of expression in the daughter copy relaxes functional constraint on its protein coding sequence leading to the asymmetric accumulation of nonsynonymous mutations. Thus, our data both affirm Ohno's prediction while explaining the mechanism by which functional divergence initially occurs following duplication for dispersed gene duplicates.

A toolkit for bulk PCR-based marker design from next-generation sequence data: application for development of a framework linkage map in bulb onion (Allium cepa L.)

Background

Although modern sequencing technologies permit the ready detection of numerous DNA sequence variants in any organisms, converting such information to PCR-based genetic markers is hampered by a lack of simple, scalable tools. Onion is an example of an under-researched crop with a complex, heterozygous genome where genome-based research has previously been hindered by limited sequence resources and genetic markers.

Results

We report the development of generic tools for large-scale web-based PCR-based marker design in the Galaxy bioinformatics framework, and their application for development of next-generation genetics resources in a wide cross of bulb onion (Allium cepa L.). Transcriptome sequence resources were developed for the homozygous doubled-haploid bulb onion line ‘CUDH2150’ and the genetically distant Indian landrace ‘Nasik Red’, using 454™ sequencing of normalised cDNA libraries of leaf and shoot. Read mapping of ‘Nasik Red’ reads onto ‘CUDH2150’ assemblies revealed 16836 indel and SNP polymorphisms that were mined for portable PCR-based marker development. Tools for detection of restriction polymorphisms and primer set design were developed in BioPython and adapted for use in the Galaxy workflow environment, enabling large-scale and targeted assay design. Using PCR-based markers designed with these tools, a framework genetic linkage map of over 800cM spanning all chromosomes was developed in a subset of 93 F₂ progeny from a very large F₂ family developed from the ‘Nasik Red’ x ‘CUDH2150’ inter-cross. The utility of tools and genetic resources developed was tested by designing markers to transcription factor-like polymorphic sequences. Bin mapping these markers using a subset of 10 progeny confirmed the ability to place markers within 10 cM bins, enabling increased efficiency in marker assignment and targeted map refinement. The major genetic loci conditioning red bulb colour (R) and fructan content (Frc) were located on this map by QTL analysis.

Conclusions

The generic tools developed for the Galaxy environment enable rapid development of sets of PCR assays targeting sequence variants identified from Illumina and 454 sequence data. They enable non-specialist users to validate and exploit large volumes of next-generation sequence data using basic equipment.

Characterization, fine mapping and expression profiling of Ragged leaves1 in maize

The Ragged leaves1 (Rg1) maize mutant frequently develops lesions on leaves, leaf sheaths, and ear bracts. Lesion formation is independent of biotic stress. High-level accumulation of H(2)O(2) revealed by staining Rg1 leaves, with 3',3'-diaminobenzidine and trypan blue, suggested that lesion formation appeared to be due to cell death. Rg1 was initially mapped to an interval around 70.5 Mb in bin 3.04 on the short arm of chromosome 3. Utilizing 15 newly developed markers, Rg1 was delimitated to an interval around 17 kb using 16,356 individuals of a BC1 segregating population. There was only one gene, rp3, predicted in this region according to the B73 genome. Analysis of transcriptome data revealed that 441 genes significantly up-regulated in Rg1 leaves were functionally over-represented. Among those genes, several were involved in the production of reactive oxygen species (ROS). Our results suggested that lesions of Rg1 maize arose probably due to an aberrant rust resistance allele of Rp3, which elicited the accumulation of ROS independent of biotic stress.

An improved allele-specific PCR primer design method for SNP marker analysis and its application

BACKGROUND

Although Single Nucleotide Polymorphism (SNP) marker is an invaluable tool for positional cloning, association study and evolutionary analysis, low SNP detection efficiency by Allele-Specific PCR (AS-PCR) still restricts its application as molecular marker like other markers such as Simple Sequence Repeat (SSR). To overcome this problem, primers with a single nucleotide artificial mismatch introduced within the three bases closest to the 3'end (SNP site) have been used in AS-PCR. However, for one SNP site, nine possible mismatches can be generated among the three bases and how to select the right one to increase primer specificity is still a challenge.

RESULTS

In this study, different from the previous reports which used a limited quantity of primers randomly (several or dozen pairs), we systematically investigated the effects of mismatch base pairs, mismatch sites and SNP types on primer specificity with 2071 primer pairs, which were designed based on SNPs from Brassica oleracea 01-88 and 02-12. According to the statistical results, we (1) found that the primers designed with SNP (A/T), in which the mismatch (CA) in the 3rd nucleotide from the 3' end, had the highest allele-specificity (81.9%). This information could be used when designing primers from a large quantity of SNP sites; (2) performed the primer design principle which forms the one and only best primer for every SNP type. This is never reported in previous studies. Additionally, we further identified its availability in rapeseed (Brassica napus L.) and sesame (Sesamum indicum). High polymorphism percent (75%) of the designed primers indicated it is a general method and can be applied in other species.

CONCLUSION

The method provided in this study can generate primers more effectively for every SNP site compared to other AS-PCR primer design methods. The high allele-specific efficiency of the SNP primer allows the feasibility for low- to moderate- throughput SNP analyses and is much suitable for gene mapping, map-based cloning, and marker-assisted selection in crops.

High resolution mapping of Dense spike-ar (dsp.ar) to the genetic centromere of barley chromosome 7H

Spike density in barley is under the control of several major genes, as documented previously by genetic analysis of a number of morphological mutants. One such class of mutants affects the rachis internode length leading to dense or compact spikes and the underlying genes were designated dense spike (dsp). We previously delimited two introgressed genomic segments on chromosome 3H (21 SNP loci, 35.5 cM) and 7H (17 SNP loci, 20.34 cM) in BW265, a BC(7)F(3) nearly isogenic line (NIL) of cv. Bowman as potentially containing the dense spike mutant locus dsp.ar, by genotyping 1,536 single nucleotide polymorphism (SNP) markers in both BW265 and its recurrent parent. Here, the gene was allocated by high-resolution bi-parental mapping to a 0.37 cM interval between markers SC57808 (Hv_SPL14)-CAPSK06413 residing on the short and long arm at the genetic centromere of chromosome 7H, respectively. This region putatively contains more than 800 genes as deduced by comparison with the collinear regions of barley, rice, sorghum and Brachypodium, Classical map-based isolation of the gene dsp.ar thus will be complicated due to the infavorable relationship of genetic to physical distances at the target locus.

RAD tag sequencing as a source of SNP markers in Cynara cardunculus L

BACKGROUND

The globe artichoke (Cynara cardunculus L. var. scolymus) genome is relatively poorly explored, especially compared to those of the other major Asteraceae crops sunflower and lettuce. No SNP markers are in the public domain. We have combined the recently developed restriction-site associated DNA (RAD) approach with the Illumina DNA sequencing platform to effect the rapid and mass discovery of SNP markers for C. cardunculus.

RESULTS

RAD tags were sequenced from the genomic DNA of three C. cardunculus mapping population parents, generating 9.7 million reads, corresponding to ~1 Gbp of sequence. An assembly based on paired ends produced ~6.0 Mbp of genomic sequence, separated into ~19,000 contigs (mean length 312 bp), of which ~21% were fragments of putative coding sequence. The shared sequences allowed for the discovery of ~34,000 SNPs and nearly 800 indels, equivalent to a SNP frequency of 5.6 per 1,000 nt, and an indel frequency of 0.2 per 1,000 nt. A sample of heterozygous SNP loci was mapped by CAPS assays and this exercise provided validation of our mining criteria. The repetitive fraction of the genome had a high representation of retrotransposon sequence, followed by simple repeats, AT-low complexity regions and mobile DNA elements. The genomic k-mers distribution and CpG rate of C. cardunculus, compared with data derived from three whole genome-sequenced dicots species, provided a further evidence of the random representation of the C. cardunculus genome generated by RAD sampling.

CONCLUSION

The RAD tag sequencing approach is a cost-effective and rapid method to develop SNP markers in a highly heterozygous species. Our approach permitted to generate a large and robust SNP datasets by the adoption of optimized filtering criteria.

High resolution mapping of Dense spike-ar (dsp.ar) to the genetic centromere of barley chromosome 7H

Spike density in barley is under the control of several major genes, as documented previously by genetic analysis of a number of morphological mutants. One such class of mutants affects the rachis internode length leading to dense or compact spikes and the underlying genes were designated dense spike (dsp). We previously delimited two introgressed genomic segments on chromosome 3H (21 SNP loci, 35.5 cM) and 7H (17 SNP loci, 20.34 cM) in BW265, a BC(7)F(3) nearly isogenic line (NIL) of cv. Bowman as potentially containing the dense spike mutant locus dsp.ar, by genotyping 1,536 single nucleotide polymorphism (SNP) markers in both BW265 and its recurrent parent. Here, the gene was allocated by high-resolution bi-parental mapping to a 0.37 cM interval between markers SC57808 (Hv_SPL14)-CAPSK06413 residing on the short and long arm at the genetic centromere of chromosome 7H, respectively. This region putatively contains more than 800 genes as deduced by comparison with the collinear regions of barley, rice, sorghum and Brachypodium, Classical map-based isolation of the gene dsp.ar thus will be complicated due to the infavorable relationship of genetic to physical distances at the target locus.

Development and use of genic molecular markers (GMMs) for construction of a transcript map of chickpea (Cicer arietinum L.)

A transcript map has been constructed by the development and integration of genic molecular markers (GMMs) including single nucleotide polymorphism (SNP), genic microsatellite or simple sequence repeat (SSR) and intron spanning region (ISR)-based markers, on an inter-specific mapping population of chickpea, the third food legume crop of the world and the first food legume crop of India. For SNP discovery through allele re-sequencing, primer pairs were designed for 688 genes/expressed sequence tags (ESTs) of chickpea and 657 genes/ESTs of closely related species of chickpea. High-quality sequence data obtained for 220 candidate genic regions on 2-20 genotypes representing 9 Cicer species provided 1,893 SNPs with an average frequency of 1/35.83 bp and 0.34 PIC (polymorphism information content) value. On an average 2.9 haplotypes were present in 220 candidate genic regions with an average haplotype diversity of 0.6326. SNP2CAPS analysis of 220 sequence alignments, as mentioned above, provided a total of 192 CAPS candidates. Experimental analysis of these 192 CAPS candidates together with 87 CAPS candidates identified earlier through in silico mining of ESTs provided scorable amplification in 173 (62.01%) cases of which predicted assays were validated in 143 (82.66%) cases (CGMM). Alignments of chickpea unigenes with Medicago truncatula genome were used to develop 121 intron spanning region (CISR) markers of which 87 yielded scorable products. In addition, optimization of 77 EST-derived SSR (ICCeM) markers provided 51 scorable markers. Screening of easily assayable 281 markers including 143 CGMMs, 87 CISRs and 51 ICCeMs on 5 parental genotypes of three mapping populations identified 104 polymorphic markers including 90 markers on the inter-specific mapping population. Sixty-two of these GMMs together with 218 earlier published markers (including 64 GMM loci) and 20 other unpublished markers could be integrated into this genetic map. A genetic map developed here, therefore, has a total of 300 loci including 126 GMM loci and spans 766.56 cM, with an average inter-marker distance of 2.55 cM. In summary, this is the first report on the development of large-scale genic markers including development of easily assayable markers and a transcript map of chickpea. These resources should be useful not only for genome analysis and genetics and breeding applications of chickpea, but also for comparative legume genomics.

An integrated pipeline of open source software adapted for multi-CPU architectures: use in the large-scale identification of single nucleotide polymorphisms

The large amounts of EST sequence data available from a single species of an organism as well as for several species within a genus provide an easy source of identification of intra- and interspecies single nucleotide polymorphisms (SNPs). In the case of model organisms, the data available are numerous, given the degree of redundancy in the deposited EST data. There are several available bioinformatics tools that can be used to mine this data; however, using them requires a certain level of expertise: the tools have to be used sequentially with accompanying format conversion and steps like clustering and assembly of sequences become time-intensive jobs even for moderately sized datasets. We report here a pipeline of open source software extended to run on multiple CPU architectures that can be used to mine large EST datasets for SNPs and identify restriction sites for assaying the SNPs so that cost-effective CAPS assays can be developed for SNP genotyping in genetics and breeding applications. At the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), the pipeline has been implemented to run on a Paracel high-performance system consisting of four dual AMD Opteron processors running Linux with MPICH. The pipeline can be accessed through user-friendly web interfaces at http://hpc.icrisat.cgiar.org/PBSWeb and is available on request for academic use. We have validated the developed pipeline by mining chickpea ESTs for interspecies SNPs, development of CAPS assays for SNP genotyping, and confirmation of restriction digestion pattern at the sequence level.

SNP-RFLPing 2: an updated and integrated PCR-RFLP tool for SNP genotyping

Background

PCR-restriction fragment length polymorphism (RFLP) assay is a cost-effective method for SNP genotyping and mutation detection, but the manual mining for restriction enzyme sites is challenging and cumbersome. Three years after we constructed SNP-RFLPing, a freely accessible database and analysis tool for restriction enzyme mining of SNPs, significant improvements over the 2006 version have been made and incorporated into the latest version, SNP-RFLPing 2.

Results

The primary aim of SNP-RFLPing 2 is to provide comprehensive PCR-RFLP information with multiple functionality about SNPs, such as SNP retrieval to multiple species, different polymorphism types (bi-allelic, tri-allelic, tetra-allelic or indels), gene-centric searching, HapMap tagSNPs, gene ontology-based searching, miRNAs, and SNP500Cancer. The RFLP restriction enzymes and the corresponding PCR primers for the natural and mutagenic types of each SNP are simultaneously analyzed. All the RFLP restriction enzyme prices are also provided to aid selection. Furthermore, the previously encountered updating problems for most SNP related databases are resolved by an on-line retrieval system.

Conclusions

The user interfaces for functional SNP analyses have been substantially improved and integrated. SNP-RFLPing 2 offers a new and user-friendly interface for RFLP genotyping that can be used in association studies and is freely available at http://bio.kuas.edu.tw/snp-rflping2.

A comprehensive resource of drought- and salinity- responsive ESTs for gene discovery and marker development in chickpea (Cicer arietinum L.)

BACKGROUND

Chickpea (Cicer arietinum L.), an important grain legume crop of the world is seriously challenged by terminal drought and salinity stresses. However, very limited number of molecular markers and candidate genes are available for undertaking molecular breeding in chickpea to tackle these stresses. This study reports generation and analysis of comprehensive resource of drought- and salinity-responsive expressed sequence tags (ESTs) and gene-based markers.

RESULTS

A total of 20,162 (18,435 high quality) drought- and salinity- responsive ESTs were generated from ten different root tissue cDNA libraries of chickpea. Sequence editing, clustering and assembly analysis resulted in 6,404 unigenes (1,590 contigs and 4,814 singletons). Functional annotation of unigenes based on BLASTX analysis showed that 46.3% (2,965) had significant similarity (< or =1E-05) to sequences in the non-redundant UniProt database. BLASTN analysis of unique sequences with ESTs of four legume species (Medicago, Lotus, soybean and groundnut) and three model plant species (rice, Arabidopsis and poplar) provided insights on conserved genes across legumes as well as novel transcripts for chickpea. Of 2,965 (46.3%) significant unigenes, only 2,071 (32.3%) unigenes could be functionally categorised according to Gene Ontology (GO) descriptions. A total of 2,029 sequences containing 3,728 simple sequence repeats (SSRs) were identified and 177 new EST-SSR markers were developed. Experimental validation of a set of 77 SSR markers on 24 genotypes revealed 230 alleles with an average of 4.6 alleles per marker and average polymorphism information content (PIC) value of 0.43. Besides SSR markers, 21,405 high confidence single nucleotide polymorphisms (SNPs) in 742 contigs (with > or = 5 ESTs) were also identified. Recognition sites for restriction enzymes were identified for 7,884 SNPs in 240 contigs. Hierarchical clustering of 105 selected contigs provided clues about stress- responsive candidate genes and their expression profile showed predominance in specific stress-challenged libraries.

CONCLUSION

Generated set of chickpea ESTs serves as a resource of high quality transcripts for gene discovery and development of functional markers associated with abiotic stress tolerance that will be helpful to facilitate chickpea breeding. Mapping of gene-based markers in chickpea will also add more anchoring points to align genomes of chickpea and other legume species.

Fine mapping, physical mapping and development of diagnostic markers for the Rrs2 scald resistance gene in barley

The Rrs2 gene confers resistance to the fungal pathogen Rhynchosporium secalis which causes leaf scald, a major barley disease. The Rrs2 gene was fine mapped to an interval of 0.08 cM between markers 693M6_6 and P1D23R on the distal end of barley chromosome 7HS using an Atlas (resistant) x Steffi (susceptible) mapping population of 9,179 F(2)-plants. The establishment of a physical map of the Rrs2 locus led to the discovery that Rrs2 is located in an area of suppressed recombination within this mapping population. The analysis of 58 barley genotypes revealed a large linkage block at the Rrs2 locus extending over several hundred kb which is present only in Rrs2 carrying cultivars. Due to the lack of recombination in the mapping population and the presence of a Rrs2-specific linkage block, we assume a local chromosomal rearrangement (alien introgression or inversion) in Rrs2 carrying varieties. The variety analysis led to the discovery of eight SNPs which were diagnostic for the Rrs2 phenotype. Based on these SNPs diagnostic molecular markers (CAPS and pyrosequencing markers) were developed which are highly useful for marker-assisted selection in resistance gene pyramiding programmes for Rhynchosporium secalis resistance in barley.

Demonstration of genetic exchange during cyclical development of Leishmania in the sand fly vector

Genetic exchange has not been shown to be a mechanism underlying the extensive diversity of Leishmania parasites. We report here evidence that the invertebrate stages of Leishmania are capable of having a sexual cycle consistent with a meiotic process like that described for African trypanosomes. Hybrid progeny were generated that bore full genomic complements from both parents, but kinetoplast DNA maxicircles from one parent. Mating occurred only in the sand fly vector, and hybrids were transmitted to the mammalian host by sand fly bite. Genetic exchange likely contributes to phenotypic diversity in natural populations, and analysis of hybrid progeny will be useful for positional cloning of the genes controlling traits such as virulence, tissue tropism, and drug resistance.

Linkage mapping of putative regulator genes of barley grain development characterized by expression profiling

BACKGROUND

Barley (Hordeum vulgare L.) seed development is a highly regulated process with fine-tuned interaction of various tissues controlling distinct physiological events during prestorage, storage and dessication phase. As potential regulators involved within this process we studied 172 transcription factors and 204 kinases for their expression behaviour and anchored a subset of them to the barley linkage map to promote marker-assisted studies on barley grains.

RESULTS

By a hierachical clustering of the expression profiles of 376 potential regulatory genes expressed in 37 different tissues, we found 50 regulators preferentially expressed in one of the three grain tissue fractions pericarp, endosperm and embryo during seed development. In addition, 27 regulators found to be expressed during both seed development and germination and 32 additional regulators are characteristically expressed in multiple tissues undergoing cell differentiation events during barley plant ontogeny. Another 96 regulators were, beside in the developing seed, ubiquitously expressed among all tissues of germinating seedlings as well as in reproductive tissues. SNP-marker development for those regulators resulted in anchoring 61 markers on the genetic linkage map of barley and the chromosomal assignment of another 12 loci by using wheat-barley addition lines. The SNP frequency ranged from 0.5 to 1.0 SNP/kb in the parents of the various mapping populations and was 2.3 SNP/kb over all eight lines tested. Exploration of macrosynteny to rice revealed that the chromosomal orders of the mapped putative regulatory factors were predominantly conserved during evolution.

CONCLUSION

We identified expression patterns of major transcription factors and signaling related genes expressed during barley ontogeny and further assigned possible functions based on likely orthologs functionally well characterized in model plant species. The combined linkage map and reference expression map of regulators defined in the present study offers the possibility of further directed research of the functional role of regulators during seed development in barley.

Isolation and initial characterization of GW5, a major QTL associated with rice grain width and weight

Grain weight is a major determinant of crop grain yield and is controlled by naturally occurring quantitative trait loci (QTLs). We earlier identified a major QTL that controls rice grain width and weight, GW5, which was mapped to a recombination hotspot on rice chromosome 5. To gain a better understanding of how GW5 controls rice grain width, we conducted fine mapping of this locus and uncovered a 1 212-bp deletion associated with the increased grain width in the rice cultivar Asominori, in comparison with the slender grain rice IR24. In addition, genotyping analyses of 46 rice cultivars revealed that this deletion is highly correlated with the grain-width phenotype, suggesting that the GW5 deletion might have been selected during rice domestication. GW5 encodes a novel nuclear protein of 144 amino acids that is localized to the nucleus. Furthermore, we show that GW5 physically interacts with polyubiquitin in a yeast two-hybrid assay. Together, our results suggest that GW5 represents a major QTL underlying rice width and weight, and that it likely acts in the ubiquitin-proteasome pathway to regulate cell division during seed development. This study provides novel insights into the molecular mechanisms controlling rice grain development and suggests that GW5 could serve as a potential tool for high-yield breeding of crops.Cell Research (2008) 18:1199-1209. doi: 10.1038/cr.2008.307; published online 18 November 2008.

Molecular marker systems for Oenothera genetics

The genus Oenothera has an outstanding scientific tradition. It has been a model for studying aspects of chromosome evolution and speciation, including the impact of plastid nuclear co-evolution. A large collection of strains analyzed during a century of experimental work and unique genetic possibilities allow the exchange of genetically definable plastids, individual or multiple chromosomes, and/or entire haploid genomes (Renner complexes) between species. However, molecular genetic approaches for the genus are largely lacking. In this study, we describe the development of efficient PCR-based marker systems for both the nuclear genome and the plastome. They allow distinguishing individual chromosomes, Renner complexes, plastomes, and subplastomes. We demonstrate their application by monitoring interspecific exchanges of genomes, chromosome pairs, and/or plastids during crossing programs, e.g., to produce plastome-genome incompatible hybrids. Using an appropriate partial permanent translocation heterozygous hybrid, linkage group 7 of the molecular map could be assigned to chromosome 9.8 of the classical Oenothera map. Finally, we provide the first direct molecular evidence that homologous recombination and free segregation of chromosomes in permanent translocation heterozygous strains is suppressed.

Expression profiling and mapping of defence response genes associated with the barley-Pyrenophora teres incompatible interaction

Barley net- and spot-form of net blotch disease are caused by two formae of the hemibiotrophic fungus Pyrenophora teres (P. t. f. teres and P. t. f. maculata). In the present study, suppression subtractive hybridization (SSH) was used in combination with quantitative real-time reverse transcriptase PCR to identify and profile the expression of defence response (DR) genes in the early stages of both barley-P. teres incompatible and compatible interactions. From a pool of 307 unique gene transcripts identified by SSH, 45 candidate DR genes were selected for temporal expression profiling in infected leaf epidermis. Differential expression profiles were observed for 28 of the selected candidates, which were grouped into clusters depending on their expression profiles within the first 48 h after inoculation. The expression profiles characteristic of each gene cluster were very similar in both barley-P. t. f. teres and barley-P. t. f. maculata interactions, indicating that resistance to both pathogens could be mediated by induction of the same group of DR genes. Chromosomal map locations for 21 DR genes were identified using four doubled-haploid mapping populations. The mapped DR genes were distributed across all seven barley chromosomes, with at least one gene mapping to within 15 cM of another on chromosomes 1H, 2H, 5H and 7H. Additionally, some DR genes appeared to co-localize with loci harbouring known resistance genes or quantitative trait loci for net blotch resistance on chromosomes 6H and 7H, as well as loci associated with resistance to other barley diseases. The DR genes are discussed with respect to their map locations and potential functional role in contributing to net blotch disease resistance.

A new diagnostic SSR marker for selection of the Rym4/Rym5 locus in barley breeding

Genomic sequence AY661558, representing a part of the BAC contig of the Rym4/Rym5 locus conferring resistance to the barley yellow mosaic virus complex (BaMMV/BaYMV), was exploited in order to develop SSR markers for practical barley breeding. Out of 57 SSR motifs found within this sequence, primers were designed and tested for the 5 SSRs with the highest repeat length. The polymorphic SSR marker QLB1 co-segregated with rym4 and rym5 phenotypes in respective high-resolution mapping populations developed for the construction of the original BAC contig. The primers targeted 2 sites located 756 bp and 5173 bp downstream of the translation initiation factor 4E (Hv-eIF4E). Physical linkage of the QLB1 marker to the Rym4/Rym5 locus was confirmed experimentally on Morex BAC 519J14, a seed BAC of Hv-eIF4E, and BAC 801A11, which is located proximally to Hv-eIF4E. QLB1 revealed 7 alleles in a set of 100 winter barley lines and cultivars. Five alleles were found within 673 advanced breeding lines derived from applied Polish winter barley breeding programmes, which corresponds to a PIC value of 0.684. No recombinants between Rym4/5 and QLB1 were detected, suggesting that QLB1 can be used efficiently in marker-assisted selection of the Hv-eIF4E-mediated bymovirus resistance.

Mapping EST-derived SSRs and ESTs involved in resistance to bacterial blight in Manihot esculenta

Cassava (Manihot esculenta Crantz) is a major root crop widely grown in the tropics. Cassava bacterial blight, caused by Xanthomonas axonopodis pv. manihotis (Xam), is an important disease in Latin America and Africa resulting in significant losses. The preferred control method is the use of resistant genotypes. Mapping expressed sequence tags (ESTs) and determining their co-localization with quantitative trait loci (QTLs) may give additional evidence of the role of the corresponding genes in resistance or defense. Twenty-one EST-derived simple sequence repeats (SSRs) were mapped in 16 linkage groups. ESTs showing similarities with candidate resistance genes or defense genes were also mapped using strategies such as restriction fragment length polymorphisms, cleaved amplified polymorphic sequences, and allele-specific primers. In total, 10 defense-related genes and 2 bacterial artificial chromosomes (BACs) containing resistance gene candidates (RGCs) were mapped in 11 linkage groups. Two new QTLs associated with resistance to Xam strains CIO121 and CIO151 were detected in linkage groups A and U, respectively. The QTL in linkage group U explained 61.6% of the phenotypic variance and was associated with an RGC-containing BAC. No correlation was found between the new EST-derived SSRs or other mapped ESTs and the new or previously reported QTLs.

EST-derived single nucleotide polymorphism markers for assembling genetic and physical maps of the barley genome

In a panel of seven genotypes, 437 expressed sequence tag (EST)-derived DNA fragments were sequenced. Single nucleotide polymorphisms (SNPs) that were polymorphic between the parents of three mapping populations were mapped by heteroduplex analysis and a genome-wide consensus map comprising 216 EST-derived SNPs and 4 InDel (insertion/deletion) markers was constructed. The average frequency of SNPs amounted to 1/130 bp and 1/107.8 bp for a set of randomly selected and a set of mapped ESTs, respectively. The calculated nucleotide diversities (pi) ranged from 0 to 40.0 x 10(-3) (average 3.1 x 10(-3)) and 0.52 x 10(-3) to 39.51 x 10(-3) (average 4.37 x 10(-3)) for random and mapped ESTs, respectively. The polymorphism information content value for mapped SNPs ranged from 0.24 to 0.50 with an average of 0.34. As expected, combination of SNPs present in an amplicon (haplotype) exhibited a higher information content ranging from 0.24 to 0.85 with an average of 0.50. Cleaved amplified polymorphic sequence assays (including InDels) were designed for a total of 87 (39.5%) SNP markers. The high abundance of SNPs in the barley genome provides avenues for the systematic development of saturated genetic maps and their integration with physical maps.

A chlorophyll-deficient rice mutant with impaired chlorophyllide esterification in chlorophyll biosynthesis

Chlorophyll (Chl) synthase catalyzes esterification of chlorophyllide to complete the last step of Chl biosynthesis. Although the Chl synthases and the corresponding genes from various organisms have been well characterized, Chl synthase mutants have not yet been reported in higher plants. In this study, a rice (Oryza Sativa) Chl-deficient mutant, yellow-green leaf1 (ygl1), was isolated, which showed yellow-green leaves in young plants with decreased Chl synthesis, increased level of tetrapyrrole intermediates, and delayed chloroplast development. Genetic analysis demonstrated that the phenotype of ygl1 was caused by a recessive mutation in a nuclear gene. The ygl1 locus was mapped to chromosome 5 and isolated by map-based cloning. Sequence analysis revealed that it encodes the Chl synthase and its identity was verified by transgenic complementation. A missense mutation was found in a highly conserved residue of YGL1 in the ygl1 mutant, resulting in reduction of the enzymatic activity. YGL1 is constitutively expressed in all tissues, and its expression is not significantly affected in the ygl1 mutant. Interestingly, the mRNA expression of the cab1R gene encoding the Chl a/b-binding protein was severely suppressed in the ygl1 mutant. Moreover, the expression of some nuclear genes associated with Chl biosynthesis or chloroplast development was also affected in ygl1 seedlings. These results indicate that the expression of nuclear genes encoding various chloroplast proteins might be feedback regulated by the level of Chl or Chl precursors.