The identification and characterisation of alleles of sucrose phosphate synthase gene family III in sugarcane

Little is known about the extent of allelic diversity of genes in the complex polyploid, sugarcane. Using sucrose phosphate synthase (SPS) Gene (SPS) Family III as an example, we have amplified and sequenced a 400 nt region from this gene from two sugarcane lines that are parents of a mapping population. Ten single nucleotide polymorphisms (SNPs) were identified within the 400 nt region of which seven were present in both lines. In the elite commercial cultivar Q165^A, 10 sequence haplotypes were identified, with four haplotypes recovered at 9% or greater frequency. Based on SNP presence, two clusters of haplotypes were observed. In IJ76-514, a Saccharum officinarum accession, 8 haplotypes were identified with 4 haplotypes recovered at 13% or greater frequency. Again, two clusters of haplotypes were observed. The results suggest that there may be two SPS Gene Family III genes per genome in sugarcane, each with different numbers of different alleles. This suggestion is supported by sequencing results in an elite parental sorghum line, 403463-2-1, in which 4 haplotypes, corresponding to two broad types, were also identified. Primers were designed to the sugarcane SNPs and screened over bulked DNA from high and low Sucrose-containing progeny from a cross between Q165^A and IJ76-514. The SNP frequency did not vary in the two bulked DNA samples, suggesting that these SNPs from this SPS gene family are not associated with variation in sucrose content. Using an ecotilling approach, two of the SPS Gene Family III haplotypes were mapped to two different linkage groups in homology group 1 in Q165^A. Both haplotypes mapped near QTLs for increased sucrose content but were not themselves associated with any sugar-related trait.     

Characterisation of single nucleotide polymorphisms in sugarcane ESTs

Commercial sugarcane cultivars (Saccharum spp. hybrids) are both polyploid and aneuploid with chromosome numbers in excess of 100; these chromosomes can be assigned to 8 homology groups. To determine the utility of single nucleotide polymorphisms (SNPs) as a means of improving our understanding of the complex sugarcane genome, we developed markers to a suite of SNPs identified in a list of sugarcane ESTs. Analysis of 69 EST contigs showed a median of 9 SNPs per EST and an average of 1 SNP per 50 bp of coding sequence. The quantitative presence of each base at 58 SNP loci within 19 contiguous sequence sets was accurately and reliably determined for 9 sugarcane genotypes, including both commercial cultivars and ancestral species, through the use of quantitative light emission technology in pyrophosphate sequencing. Across the 9 genotypes tested, 47 SNP loci were polymorphic and 11 monomorphic. Base frequency at individual SNP loci was found to vary approximately twofold between Australian sugarcane cultivars and more widely between cultivars and wild species. Base quantity was shown to segregate as expected in the IJ76-514 × Q165 sugarcane mapping population, indicating that SNPs that occur on one or two sugarcane chromosomes have the potential to be mapped. The use of SNP base frequencies from five of the developed markers was able to clearly distinguish all genotypes in the population. The use of SNP base frequencies from a further six markers within an EST contig was able to help establish the likely copy number of the locus in two genotypes tested. This is the first instance of a technology that has been able to provide an insight into the copy number of a specific gene locus in hybrid sugarcane. The identification of specific and numerous haplotypes/alleles present in a genotype by pyrophosphate sequencing or alternative techniques ultimately will provide the basis for identifying associations between specific alleles and phenotype and between allele dosage and phenotype in sugarcane.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Construction of a genetic linkage map for Saccharum officinarum incorporating both simplex and duplex markers to increase genome coverage.

Saccharum officinarum L. is an octoploid with 80 chromosomes and a basic chromosome number of x = 10. It has high stem sucrose and contributes 80% of the chromosomes to the interspecific sugarcane cultivars that are grown commercially for sucrose. A genetic linkage map was developed for S. officinarum (clone IJ76-514) using a segregating population generated from a cross between Q165 (a commercial sugarcane cultivar) and IJ76-514. In total, 40 AFLP and 72 SSR primer pairs were screened across the population, revealing 595 polymorphic bands inherited from IJ76-514. These 595 markers displayed a frequency distribution different from all other sugarcane genetic maps produced, with only 40% being simplex markers (segregated 1:1). Of these 240 simplex markers, 178 were distributed on 47 linkage groups (LGs) and 62 remained unlinked. With the addition of 234 duplex markers and 80 biparental simplex markers (segregating 3:1), 534 markers formed 123 LGs. Using the multi-allelic SSR markers, repulsion phase linkage, and alignment with the Q165 linkage map, 105 of the 123 LGs could be grouped into 10 homology groups (HGs). These 10 HGs were further assigned to the 8 HGs observed in cultivated sugarcane and S. spontaneum. Analysis of repulsion phase linkage indicated that IJ76-514 is neither a complete autopolyploid nor an allopolyploid. Detection of 28 repulsion linkages that occurred between 6 pairs of LGs located in 4 HGs suggested the occurrence of limited preferential chromosome pairing in this species.     

A tailed PCR procedure for cost-effective, two-order multiplex sequencing of candidate genes in polyploid plants

Complex polyploid crop genomes can be recalcitrant towards conventional DNA sequencing approaches for allele mining in candidate genes for valuable traits. In the past, this has greatly complicated the transfer of knowledge on promising candidate genes from model plants to even closely related polyploid crops. Next-generation sequencing offers diverse solutions to overcome such difficulties. Here, we present a method for multiplexed 454 sequencing in gene-specific PCR amplicons that can simultaneously address multiple homologues of given target genes. We devised a simple two-step PCR procedure employing a set of barcoded M13/T7 universal fusion primers that enable a cost-effective and efficient amplification of large numbers of target gene amplicons. Sequencing-ready amplicons are generated that can be simultaneously sequenced in pools comprising multiple amplicons from multiple genotypes. High-depth sequencing allows resolution of the resulting sequence reads into contigs representing multiple homologous loci, with only insignificant off-target capture of paralogues or PCR artefacts. In a case study, the procedure was tested in the complex polyploid genome of Brassica napus for a set of nine genes identified in Arabidopsis as candidates for regulation of seed development and oil content. Up to six copies of these genes were expected in B.?napus. SNP discovery was performed by pooled multiplex sequencing of 30 amplicons in 20 diverse B.?napus accessions with interesting trait variation for oil content, providing a basis for comparative mapping to relevant quantitative trait loci and for subsequent marker-assisted breeding.     

The Complex Genetic Structure of Sugarcane Limits Identification of Additional SNP-Defined Simplex Alleles in Microsatellite Loci

Cultivated sugarcane possesses a large number of chromosomes (typically >80) which can be organised into eight homology groups, each containing approximately 12 homo(eo)logous chromosomes. Currently, microsatellite (SSR) markers are the most easily used markers for marker-assisted selection and other genetic applications. However, only SSR alleles that segregate as simplex and duplex markers can be incorporated into sugarcane maps using populations of ～300 progeny. Consequently only a subset of possible alleles have been mapped for a given SSR locus and are available for subsequent QTL analyses. Three sugarcane SSR loci, mSSCIR8, mSSCIR17 and mSSCIR18, were amplified, cloned and sequenced from the parents of an Australian sugarcane mapping population, IJ76-514 and . The sequences were examined to identify nucleotide sequence polymorphisms in the flanking regions that could provide additional simplex SSR allele markers. Alignment of the sequences revealed SNP, indel and repeat length variation and the pattern of sequence variation suggested multiple alleles for each SSR, including all of the alleles previously scored by fragment length. While the flanking regions of the SSR loci contained numerous SNPs and indels, none defined new simplex alleles within multiplex fragments and no new SSR simplex alleles were mapped. Furthermore, many of the sequence-defined alleles appear to be spurious and may have arisen from PCR-mediated recombination. These results confirm the difficulties associated with characterising allelic diversity in a polyploid species and the complexity of mapping in sugarcane.     

Genetic control of yield related stalk traits in sugarcane

A major focus of sugarcane variety improvement programs is to increase sugar yield, which can be accomplished by either increasing the sugar content of the cane or by increasing cane yield, as the correlation between these traits is low. We used a cross between an Australian sugarcane variety Q165, and a Saccharum officinarum accession, IJ76-514, to dissect the inheritance of yield-related traits in the complex polyploid sugarcane. A population of 227 individuals was grown in a replicated field trial and evaluated over 3 years for stalk weight, stalk diameter, stalk number, stalk length and total biomass. Over 1,000 AFLP and SSR markers were scored across the population and used to identify quantitative trait loci (QTL). In total, 27 regions were found that were significant at the 5% threshold using permutation tests with at least one trait; individually, they explained from 4 to 10% of the phenotypic variation and up to 46% were consistent across years. With the inclusion of digeneic interactions, from 28 to 60% of the variation was explained for these traits. The 27 genomic regions were located on 22 linkage groups (LGs) in six of the eight homology groups (HGs) indicating that a number of alleles or quantitative trait alleles (QTA) at each QTL contribute to the trait; from one to three alleles had an effect on the traits for each QTL identified. Alleles of a candidate gene, TEOSINTE BRANCHED 1 (TB1), the major gene controlling branching in maize, were mapped in this population using either an SSR or SNP markers. Two alleles showed some association with stalk number, but unlike maize, TB1 is not a major gene controlling branching in sugarcane but only has a minor and variable effect.     

Conversion of AFLP markers to high-throughput markers in a complex polyploid,sugarcane

The application of DNA markers linked to traits of commercial value in sugarcane may increase the efficiency of sugarcane breeding. The majority of markers generated for quantitative trait locus mapping in sugarcane have been single sequence repeats or AFLPs (amplified fragment length polymorphisms). Since AFLP markers are not adapted for large-scale implementation in plant breeding, our objective was to assess the feasibility of converting AFLP markers to fast, cheap and reliable PCR-based assays in a complex polyploid, sugarcane. Three AFLP markers were selected on the basis of an association to resistance to the fungal pathogen Ustilago scitaminea, the causal agent of smut in sugarcane. We developed an approach which enabled the identification of polymorphisms in these AFLP markers. Towards this goal, we employed GenomeWalking and 454 sequencing to isolate sequences adjacent to the linked AFLP markers and identify SNP (single nucleotide polymorphisms) haplotypes present in the homo(eo)logous chromosomes of sugarcane. One AFLP marker was converted to a cleavage amplified polymorphic sequence marker, another to a SCAR (sequence characteristered amplified region) marker and the final AFLP marker to a SNP PCR-based assay. However, validation of each of the markers in 240 genotypes resulted in 99, 90 and 60% correspondence with the original AFLP marker. These experiments indicate that even in a complex polyploid such as sugarcane, polymorphisms identified by AFLP can be converted to high-throughput marker systems, but due to the complexity this would only be carried out for high-value markers. In some cases, the polymorphisms identified are not transferable to more sequence-specific PCR applications.     

A Survey Sequence Comparison of Saccharum Genotypes Reveals Allelic Diversity Differences

Sugarcane (Saccharum spp.) is a crop of substantial international significance for both food and fuel, however its highly polyploid nature challenges investigation of its genetic composition. Efforts to generate the full sugarcane genome sequence are underway, however in the meantime crop improvement efforts are somewhat limited by the lack of genome sequence resources available for physiological characterization. Low-coverage survey sequence data was generated and assembled for six sugarcane genotypes representing a range of significant S. spontaneum, S. officinarum, and S. hybrid cultivar accessions from around the world. These data were explored to investigate the composition of repetitive sequences and variations in chloroplast genome sequence, as well as assembled into a conglomerate monoploid genome sequence for polymorphism comparison between the genotypes. Almost half (47 %) of the inter-genomic polymorphisms analysed in these data represented poly-allelic variations which cannot be applied in traditional present/absent marker analysis, suggesting that new approaches are required to better understand and access genetic diversity within the Saccharum genus. These results support previous assertions that S. spontaneum is both less repetitive (62 % repetitive k-mers in Mandalay vs. 65 % in IJ76-514) and more highly polymorphic (17 % poly-alleles in Mandalay vs. 10 % poly-alleles in IJ76-514) than S. officinarum, with S. hybrids being intermediate between the two. However, contrary to previous analysis the monoploid genome size of S. spontaneum does not appear to differ significantly from that of S. officinarum as had been expected. This genomic survey assembly will be a very useful resource for sugarcane genomics in the absence of a monoploid or polyploid genome sequence, and will be made available upon request.     

Enrichment of genomic DNA for polymorphism detection in a non-model highly polyploid crop plant

Large polyploid genomes of non-model species remain challenging targets for DNA polymorphism discovery despite the increasing throughput and continued reductions in cost of sequencing with new technologies. For these species especially, there remains a requirement to enrich genomic DNA to discover polymorphisms in regions of interest because of large genome size and to provide the sequence depth to enable estimation of copy number. Various methods of enriching DNA have been utilised, but some recent methods enable the efficient sampling of large regions (e.g. the exome). We have utilised one of these methods, solution-based hybridization (Agilent SureSelect), to capture regions of the genome of two sugarcane genotypes (one Saccharum officinarum and one Saccharum hybrid) based mainly on gene sequences from the close relative Sorghum bicolor. The capture probes span approximately 5.8?megabases (Mb). The enrichment over whole-genome shotgun sequencing was 10-11-fold for the two genotypes tested. This level of enrichment has important consequences for detecting single nucleotide polymorphisms (SNPs) from a single lane of Illumina (Genome Analyzer) sequence reads. The detection of polymorphisms was enabled by the depth of sequence at or near probe sites and enabled the detection of 270?000-280?000 SNPs within each genotype from a single lane of sequence using stringent detection parameters. The SNPs were present in 13?000-16?000 targeted genes, which would enable mapping of a large number of these chosen genes. SNP validation from 454 sequencing and between-genotype confirmations gave an 87%-91% validation rate.     

High resolution melting analysis of almond SNPs derived from ESTs

High resolution melting curve (HRM) is a recent advance for the detection of SNPs. The technique measures temperature induced strand separation of short PCR amplicons, and is able to detect variation as small as one base difference between samples. It has been applied to the analysis and scan of mutations in the genes causing human diseases. In plant species, the use of this approach is limited. We applied HRM analysis to almond SNP discovery and genotyping based on the predicted SNP information derived from the almond and peach EST database. Putative SNPs were screened from almond and peach EST contigs by HRM analysis against 25 almond cultivars. All 4 classes of SNPs, INDELs and microsatellites were discriminated, and the HRM profiles of 17 amplicons were established. The PCR amplicons containing single, double and multiple SNPs produced distinctive HRM profiles. Additionally, different genotypes of INDEL and microsatellite variations were also characterised by HRM analysis. By sequencing the PCR products, 100 SNPs were validated/revealed in the HRM amplicons and their flanking regions. The results showed that the average frequency of SNPs was 1:114 bp in the genic regions, and transition to transversion ratio was 1.16:1. Rare allele frequencies of the SNPs varied from 0.02 to 0.5, and the polymorphic information contents of the SNPs were from 0.04 to 0.53 at an average of 0.31. HRM has been demonstrated to be a fast, low cost, and efficient approach for SNP discovery and genotyping, in particular, for species without much genomic information such as almond.     

Expressed sequence tags for the chicken genome from a normalized 10-day-old White Leghorn whole embryo cDNA library: 1. DNA sequence characterization and linkage analysis

Expressed sequence tags (ESTs) provide a rapid and reliable method for gene discovery as well as a resource for the large-scale analysis of gene expression of known and unknown genes. Here we describe a normalized cDNA library developed from a 10-day-old White Leghorn chicken whole embryo. The utility of the library was evaluated by partial sequencing of 99 randomly selected insert-containing clones and the analysis of EST-targeted genomic regions for single nucleotide polymorphisms (SNPs) in the East Lansing chicken reference DNA mapping panel. Using stringent match criteria of percent identity of 80 or higher across a length of 50 or more bases, 46 ESTs matched database sequences including previously reported Gallus gallus genes. Thirty-seven of the 50 primer pairs developed from 50 unique ESTs amplified a single fragment. The size of the 37 amplicons ranged from 276 to 693 bp for a total of 17,508 and an average of 473. About 70% of the SNPs detected were either G-->A or C-->T transition. The number of SNPs detected within the amplicons from EST-targeted genomic regions ranged from 0 to 4 for a total of 65 and a frequency of about 1 every 470 bases. About 35% of the amplicons contained only 1 SNP, while 19% had 4 SNPs. Using the SNPs that were informative in the East Lansing reference panel, 17 ESTs were mapped on the East Lansing chicken genetic map. The ESTs described, as well as the nucleotide variants identified within the EST-targeted genomic regions, represent significant resources for genome analysis in the chicken.     

Single nucleotide polymorphism (SNP) discovery in porcine expressed genes

High-throughput genotyping of swine populations is a potentially efficient method for establishing animal lineage and identification of loci important to animal health and efficient pork production. Markers were developed based upon single nucleotide polymorphisms (SNPs), which are abundant and amenable to automated genotyping platforms. The focus of this research was SNP discovery in expressed porcine genes providing markers to develop the porcine/human comparative map. Locus specific amplification (LSA) and comparative sequencing were used to generate PCR products and allelic information from parents of a swine reference family. Discovery of 1650 SNPs in 403 amplicons and strategies for optimizing LSA-based SNP discovery using alternative methods of PCR primer design, data analysis, and germplasm selection that are applicable to other populations and species are described. These data were the first large-scale assessment of frequency and distribution of porcine SNPs.     

Multiplex restriction amplicon sequencing: a novel next‐generation sequencing‐based marker platform for high‐throughput genotyping

To enable rapid selection of traits in marker‐assisted breeding, markers must be technically simple, low‐cost, high‐throughput and randomly distributed in a genome. We developed such a technology, designated as Multiplex Restriction Amplicon Sequencing (MRASeq), which reduces genome complexity by polymerase chain reaction (PCR) amplification of amplicons flanked by restriction sites. The first PCR primers contain restriction site sequences at 3’‐ends, preceded by 6‐10 bases of specific or degenerate nucleotide sequences and then by a unique M13‐tail sequence which serves as a binding site for a second PCR that adds sequencing primers and barcodes to allow sample multiplexing for sequencing. The sequences of restriction sites and adjacent nucleotides can be altered to suit different species. Physical mapping of MRASeq SNPs from a biparental population of allohexaploid wheat (Triticum aestivum L.) showed a random distribution of SNPs across the genome. MRASeq generated thousands of SNPs from a wheat biparental population and natural populations of wheat and barley (Hordeum vulgare L.). This novel, next‐generation sequencing‐based genotyping platform can be used for linkage mapping to screen quantitative trait loci (QTL), background selection in breeding and many other genetics and breeding applications of various species.     

Use of bovine EST data and human genomic sequences to map 100 gene-specific bovine markers

A system to use bovine EST data in conjunction with human genomic sequence to improve the bovine linkage map over the entire genome or on specific chromosomes was evaluated. Bovine EST sequence was used to provide primer sequences corresponding to bovine genes, while human genomic sequence directed primer design to flank introns and produce amplicons of appropriate size for efficient direct sequencing. The sequence tagged sites (STS) produced in this way from the four sires of the MARC reference families were examined for single nucleotide polymorphisms (SNPs) that could be used to map the corresponding genes. With this approach, along with a primer/extension mass spectrometry SNP genotyping assay, 100 ESTs were placed on the bovine genetic linkage map. The first 70 were chosen at random from bovine EST–human genomic comparisons. An additional 30 ESTs were successfully mapped to bovine Chromosome 19 (BTA19), and comparison of the resulting BTA19 map to the position of the corresponding human orthologs on the HSA17 draft sequences revealed differences in the spacing and order of genes. Over 80% of successful amplicons contained SNPs, indicating that this is an efficient approach to generating EST-associated genetic markers. We have demonstrated the feasibility of constructing a linkage map based on SNPs associated with ESTs and the plausibility of utilizing EST, comparative mapping information, and human sequence data to target regions of the bovine genome for SNP marker development.     

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

Key message

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.

Abstract

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.

     

Ordered catenation of sequence-tagged sites and multiplexed SNP genotyping by sequencing

We describe a method for the efficient genotyping of SNPs, involving sequencing of ordered and catenated sequence-tagged sites (OCS). In OCS, short genomic segments, each containing an SNP, are amplified by PCR using primers that carry specially designed extra nucleotides at their 5′-ends. Amplification products are then combined and converted to a concatamer in a defined order by a second round of thermal cycling. The concatenation takes place because the 5′-ends of each amplicon are designed to be complementary to the ends of the presumptive neighboring amplicons. The primer sequences for OCS are chosen using newly developed dedicated software, OCS Optimizer. Using sets of SNPs, we show that at least 10 STSs can be concatenated in a predefined order and all SNPs in the STSs are accurately genotyped by one two-way sequencing reaction.     

A comprehensive genetic map of sugarcane that provides enhanced map coverage and integrates high-throughput Diversity Array Technology (DArT) markers

Background

Sugarcane genetic mapping has lagged behind other crops due to its complex autopolyploid genome structure. Modern sugarcane cultivars have from 110-120 chromosomes and are in general interspecific hybrids between two species with different basic chromosome numbers: Saccharum officinarum (2n = 80) with a basic chromosome number of 10 and S. spontaneum (2n = 40-128) with a basic chromosome number of 8. The first maps that were constructed utilised the single dose (SD) markers generated using RFLP, more recent maps generated using AFLP and SSRs provided at most 60% genome coverage. Diversity Array Technology (DArT) markers are high throughput allowing greater numbers of markers to be generated.

Results

Progeny from a cross between a sugarcane variety Q165 and a S. officinarum accession IJ76-514 were used to generate 2467 SD markers. A genetic map of Q165 was generated containing 2267 markers, These markers formed 160 linkage groups (LGs) of which 147 could be placed using allelic information into the eight basic homology groups (HGs) of sugarcane. The HGs contained from 13 to 23 LGs and from 204 to 475 markers with a total map length of 9774.4 cM and an average density of one marker every 4.3 cM. Each homology group contained on average 280 markers of which 43% were DArT markers 31% AFLP, 16% SSRs and 6% SNP markers. The multi-allelic SSR and SNP markers were used to place the LGs into HGs.

Conclusions

The DArT array has allowed us to generate and map a larger number of markers than ever before and consequently to map a larger portion of the sugarcane genome. This larger number of markers has enabled 92% of the LGs to be placed into the 8 HGs that represent the basic chromosome number of the ancestral species, S. spontaneum. There were two HGs (HG2 and 8) that contained larger numbers of LGs verifying the alignment of two sets of S. officinarum chromosomes with one set of S. spontaneum chromosomes and explaining the difference in basic chromosome number between the two ancestral species. There was also evidence of more complex structural differences between the two ancestral species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-152) contains supplementary material, which is available to authorized users.