High Resolution Melting Analysis Is a More Sensitive and Effective Alternative to Gel-Based Platforms in Analysis of SSR - An Example in Citrus

High resolution melting curve analysis (HRM) has been used as an efficient, accurate and cost-effective tool to detect single nucleotide polymorphisms (SNPs) or insertions or deletions (INDELs). However, its efficiency, accuracy and applicability to discriminate microsatellite polymorphism have not been extensively assessed. The traditional protocols used for SSR genotyping include PCR amplification of the DNA fragment and the separation of the fragments on electrophoresis-based platform. However, post-PCR handling processes are laborious and costly. Furthermore, SNPs present in the sequences flanking repeat motif cannot be detected by polyacrylamide-gel-electrophoresis based methods. In the present study, we compared the discriminating power of HRM with the traditional electrophoresis-based methods and provided a panel of primers for HRM genotyping in Citrus. The results showed that sixteen SSR markers produced distinct polymorphic melting curves among the Citrus spp investigated through HRM analysis. Among those, 10 showed more genotypes by HRM analysis than capillary electrophoresis owing to the presence of SNPs in the amplicons. For the SSR markers without SNPs present in the flanking region, HRM also gave distinct melting curves which detected same genotypes as were shown in capillary electrophoresis (CE) analysis. Moreover, HRM analysis allowed the discrimination of most of the 15 citrus genotypes and the resulting genetic distance analysis clustered them into three main branches. In conclusion, it has been approved that HRM is not only an efficient and cost-effective alternative of electrophoresis-based method for SSR markers, but also a method to uncover more polymorphisms contributed by SNPs present in SSRs. It was therefore suggested that the panel of SSR markers could be used in a variety of applications in the citrus biodiversity and breeding programs using HRM analysis. Furthermore, we speculate that the HRM analysis can be employed to analyse SSR markers in a wide range of applications in all other species.     

Development of EST-SSR markers in castor bean (Ricinus communis L.) and their utilization for genetic purity testing of hybrids

Expressed sequence tag (EST) databases offer opportunity for the rapid development of simple sequence repeat (SSR) markers in crops. Sequence assembly and clustering of 57?895 ESTs of castor bean resulted in the identification of 10?960 unigenes (6459 singletons and 4501 contigs) having 7429 SSRs. On an average, the unigenes contained 1 SSR for every 1.23?kb of unigene sequence. The identified SSRs mostly consisted of dinucleotide (62.4%) and trinucleotide (33.5%) repeats. The AG class was the most common among the dinucleotide motifs (68.9%), whereas the AAG class (25.9%) was predominant among the trinucleotide motifs. A total of 611 primer pairs were designed for the SSRs, having repeat length more than or equal to 20 nucleotides, of which a set of 130 markers were tested and 92 of these yielding robust amplicons were analyzed for their utility in genetic purity assessment of castor bean hybrids. Nine markers were able to detect polymorphism between the parental lines of nine commercial castor bean hybrids (DCH-32, DCH-177, DCH-519, GCH-2, GCH-4, GCH-5, GCH-6, GCH-7, and RHC-1), and their utility in genetic purity testing was demonstrated. These novel EST-SSR markers would be a valuable addition to the growing molecular marker resources that could be used in genetic improvement programmes of castor bean.     

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.     

Construction and characterization of a cDNA library from floral organs and fruitlets of <Emphasis Type="Italic">Citrus reticulata</Emphasis>

To explore and isolate genes related to flowering and fruit development, we constructed a cDNA library from floral organs and fruitlets of Ponkan mandarin (Citrus reticulata Blanco). A total of 661 high-quality expressed sequence tags (ESTs) were generated and submitted to GenBank with the accession numbers from GO343532 to GO344192. All these ESTs were assembled into 43 contigs and 296 singletons (totally 339 unigenes). The BLAST2GO software was employed to annotate the unigenes, among which 77 ones had no significant homology with the sequences in NCBI non-redundant proteins database by BLASTX analysis. Additionally, gene ontology (GO) analysis revealed an overview of sequences distribution, which implied some specially expressed genes related to flower and fruit development. Furthermore, some abundantly expressed unigenes involved in several crucial metabolic pathways related to fruit quality were highlighted and three types of homologues of miraculin-like protein2 were analyzed by both semiquantitative RT-PCR and real-time PCR. The results showed different expression profiles of these genes, which meant that they contribute distinctly to fruit development.     

Development and annotation of perennial Triticeae ESTs and SSR markers

Triticeae contains hundreds of species of both annual and perennial types. Although substantial genomic tools are available for annual Triticeae cereals such as wheat and barley, the perennial Triticeae lack sufficient genomic resources for genetic mapping or diversity research. To increase the amount of sequence information available in the perennial Triticeae, three expressed sequence tag (EST) libraries were developed and annotated for Pseudoroegneria spicata, a mixture of both Elymus wawawaiensis and E. lanceolatus, and a Leymus cinereus x L. triticoides interspecific hybrid. The ESTs were combined into unigene sets of 8 780 unigenes for P. spicata, 11 281 unigenes for Leymus, and 7 212 unigenes for Elymus. Unigenes were annotated based on putative orthology to genes from rice, wheat, barley, other Poaceae, Arabidopsis, and the non-redundant database of the NCBI. Simple sequence repeat (SSR) markers were developed, tested for amplification and polymorphism, and aligned to the rice genome. Leymus EST markers homologous to rice chromosome 2 genes were syntenous on Leymus homeologous groups 6a and 6b (previously 1b), demonstrating promise for in silico comparative mapping. All ESTs and SSR markers are available on an EST information management and annotation database (http://titan.biotec.uiuc.edu/triticeae/).     

柑橘EST-SSR分子标记分析

对来源于甜橙（Citrus sinensis Osbeck）、枳壳（Poncirus trifoliata Raf．）和其他柑橘非冗余EST数据库的38124条单-基因（Unigene）序列进行了简单重复序列SSRs（Simple Sequence Repeat）搜索,所分析的柑橘非冗余核酸序列总长23．29Mb,从中获得了8218条SSR,其中包括单碱基重复4913条（59．8％）,2碱基重复1419条（17．3％）,3碱基重复1709条（20．8％）,4碱基重复114条（1．39％）,5碱基重复23条（0．28％）,6碱基重复40条（0．49％）。大约每2．8kb长度的单-基因序列中即存在1个SSR,即平均4．6个单-基因中存在1个SSR。随碱基重复单元（motif）的不同,SSR的最大长度在40-105之间,全部重复序列的平均长度为20．9bp。各种SSR（1-,2-,3-,4-,5-,6-核苷酸重复）的发生频率在甜橙和枳壳间非常接近。其中单碱基重复序列是最丰富的重复单元,其次为3碱基重复。在所得的SSR的重复单元中,富含A碱基的重复单元的分布占据优势地位,出现的频率与密度均较高,而富含CG碱基的重复单元出现频率和密度较低。用25对EST-SSR引物对6个柑橘品种的多样性进行了PCR检测,结果表明,所有25对引物在6个柑橘品种间均扩增到多样性条带,证实通过柑橘EST数据库的发掘能够高效地筛选到基因水平的SSR标记。     

Analysis of Microsatellites in Citrus Unigenes

Simple sequence repeats (SSRs) were investigated in the unigene sequences from expressed sequence tags (EST) of sweet orange (Citrus sinensis osbeck), trifoliate orange (Poncirus trifoliata Raf.) and other citrus species and cultivars. A total of 37 802 citrus unigene sequences corresponding to 23.29 Mb were searched, resulting in the identification of 8 218 SSRs. Among them there were 4 913 (59.8%) mono-, 1 419 (17.3%) di-, 1 709 (20.8%) tri-, 114 (1.39%) tetra-, 23 (0.28%) penta- and 40 (0.49%) hexa-nucleotide SSRs. The estimated frequency of SSRs was approximately 1/2.8 kb, which could be extrapolated to 1 SSR-containing unigene in 4.6 unigenes. The maximum length of the SSR ranged from 40 to 105 bp depending on the repeating numbers of the motif in the SSR. The overall average length of SSRs was 20.9 bp. The frequencies of different SSR types (di-, tri-, tetra-, and penta-nucleotide repeats) were very similar between sweet orange and trifoliate orange. The mononucelotide repeats appeared to be the most abundant SSRs within sweet orange and trifoliate orange, followed by trimeric repeats. The adenine rich repeats such as A/T, AG, AT, AAG, AAAT, AAAG, AAAT, AAAAG, AAAAT etc. were predominant in each type of SSRs (mono-, di-, tri-, tetra-, and penta-), whereas the C/G, CG, CCG repeats were less abundant. Twenty-five primer pairs flanking EST-SSR loci were designed to detect the possible polymorphism of six citrus cultivars including sweet orange and trifoliate orange. The PCR result with all these 25 primer pairs revealed the existence of polymorphism within six citrus cultivars confirming that citrus EST database could be efficiently exploited for the development of gene-derived SSR markers.     

EST-SNP genotyping of citrus species using high-resolution melting curve analysis

Citrus taxonomy is very complex mainly due to specific aspects of its reproductive biology. A number of studies have been performed using various molecular markers in order to evaluate the level of genetic variability in Citrus. SNP markers have been used for genetic diversity assessment using a variety of different methods. Recently, the availability of EST database and whole genome sequences has made it possible to develop more markers such as SNPs. In the present study, the high-resolution melting curve analysis (HRM) was used to detect SNPs or INDELs in Citrus genus for the first time. We aimed to develop a panel of SNPs to differentiate Citrus genotypes which can also be applied to Citrus biodiversity studies. The results showed that 21 SNP containing markers produced distinct polymorphic melting curves among the Citrus spp. investigated through HRM analysis. It was proved that HRM is an efficient, cost-effective, and accurate method for discriminating citrus SNPs as well as a method to analyze more polymorphisms in a single PCR amplicon, representing a useful tool for genetic, biodiversity, and breeding studies. SNPs developed based on Citrus sinensis EST database showed a good transferability within the Citrus genus. Moreover, HRM analysis allowed the discrimination of citrus genotypes at specific level and the resulting genetic distance analysis clustered these genotypes into three main branches. The results suggested that the panel of SNP markers could be used in a variety of applications in citrus biodiversity assessment and breeding programs using HRM analysis.     

Analysis of expressed sequence tags from grapevine flower and fruit and development of simple sequence repeat markers

A total of 6,230 EST sequences were produced from 7,561 clones in a cDNA library generated from grapevine (Vitis vinifera cv. ‘Summer Black’) flower and fruit tissues in this study. After cluster and assembly analysis of the datasets, 3,582 unigenes (GenBank accession numbers GW836604–GW840185) were established, among which 381 were new grapevine EST sequences. Out of the 381 new ESTs, 289 could be mapped on the 19 grapevine chromosomes. 913 unique ESTs with known or putative functions were assigned to 11 putative cellular roles. 540 potentially workable grapevine EST-SSRs were developed from 3,582 unigenes and about 42.6% of these unigenes were identified as true-to-type SSR loci and could amplify polymorphic bands from 22 individual plants of V. vinifera L, indicating that grapevine EST datasets are a valuable source for the development of functional simple sequence repeat (SSR) markers.     

Transferability of the EST-SSRs developed on Nules clementine (<Emphasis Type="Italic">Citrus clementina</Emphasis> Hort ex Tan) to other <Emphasis Type="Italic">Citrus</Emphasis> species and their effectiveness for genetic mapping

Background

During the last decade, numerous microsatellite markers were developed for genotyping and to identify closely related plant genotypes. In citrus, previously developed microsatellite markers were arisen from genomic libraries and more often located in non coding DNA sequences. To optimize the use of these EST-SSRs as genetic markers in genome mapping programs and citrus systematic analysis, we have investigated their polymorphism related to the type (di or trinucleotide) or their position in the coding sequences.

Results

Among 11000 unigenes from a Clementine EST library, we have found at least one microsatellite sequence (repeated units size ranged from 2 to 6 nucleotides) in 1500 unigenes (13.6%). More than 95% of these SSRs were di or trinucleotides. If trinucleotide microsatellites were encountered trough all part of EST sequences, dinucleotide microsatellites were preferentially (50%) concentrated in the 5' 100th nucleotides. We assessed the polymorphism of 41 EST-SSR, by PCR amplification droved with flanking primers among ten Citrus species plus 3 from other genera. More than 90% of EST-SSR markers were polymorphic. Furthermore, dinucleotide microsatellite markers were more polymorphic than trinucleotide ones, probably related to their distribution that was more often located in the 5' UnTranslated Region (UTR). We obtained a good agreement of diversity relationships between the citrus species and relatives assessed with EST-SSR markers with the established taxonomy and phylogeny. To end, the heterozygosity of each genotype and all dual combinations were studied to evaluate the percentage of mappable markers. Higher values (> 45%) were observed for putative Citrus inter-specific hybrids (lime lemon, or sour orange) than for Citrus basic true species (mandarin, pummelo and citron) (<30%). Most favorable combinations for genome mapping were observed in those involving interspecific hybrid genotypes. Those gave higher levels of mappable markers (>70%) with a significant proportion suitable for synteny analysis.

Conclusion

Fourty one new EST-SSR markers were produced and were available for citrus genetic studies. Whatever the position of the SSR in the ESTs the EST-SSR markers we developed are powerful to investigate genetic diversity and genome mapping in citrus.

     

A high-throughput apple SNP genotyping platform using the GoldenGate™ assay

EST data generated from 14 apple genotypes were downloaded from NCBI and mapped against a reference EST assembly to identify Single Nucleotide Polymorphisms (SNPs). Mapping of these SNPs was undertaken using 90% of sequence similarity and minimum coverage of four reads at each SNP position. In total, 37,807 SNPs were identified with an average of one SNP every 187 bp from a total of 6888 unique EST contigs. Identified SNPs were checked for flanking sequences of ≥ 60 bp along both sides of SNP alleles for reliable design of a custom high-throughput genotyping assay. A total of 12,299 SNPs, representing 6525 contigs, fit the selected criterion of ≥ 60 bp sequences flanking a SNP position. Of these, 1411 SNPs were validated using four apple genotypes. Based on genotyping assays, it was estimated that 60% of SNPs were valid SNPs, while 26% of SNPs might be derived from paralogous regions.