Genome-wide SNP discovery in tetraploid alfalfa using 454 sequencing and high resolution melting analysis
- Yuanhong Han†1,
- Yun Kang†2,
- Ivone Torres-Jerez2,
- Foo Cheung3, 4,
- Christopher D Town4,
- Patrick X Zhao2,
- Michael K Udvardi2 and
- Maria J Monteros1Email author
© Han et al; licensee BioMed Central Ltd. 2011
Received: 22 March 2011
Accepted: 6 July 2011
Published: 6 July 2011
Single nucleotide polymorphisms (SNPs) are the most common type of sequence variation among plants and are often functionally important. We describe the use of 454 technology and high resolution melting analysis (HRM) for high throughput SNP discovery in tetraploid alfalfa (Medicago sativa L.), a species with high economic value but limited genomic resources.
The alfalfa genotypes selected from M. sativa subsp. sativa var. 'Chilean' and M. sativa subsp. falcata var. 'Wisfal', which differ in water stress sensitivity, were used to prepare cDNA from tissue of clonally-propagated plants grown under either well-watered or water-stressed conditions, and then pooled for 454 sequencing. Based on 125.2 Mb of raw sequence, a total of 54,216 unique sequences were obtained including 24,144 tentative consensus (TCs) sequences and 30,072 singletons, ranging from 100 bp to 6,662 bp in length, with an average length of 541 bp. We identified 40,661 candidate SNPs distributed throughout the genome. A sample of candidate SNPs were evaluated and validated using high resolution melting (HRM) analysis. A total of 3,491 TCs harboring 20,270 candidate SNPs were located on the M. truncatula (MT 3.5.1) chromosomes. Gene Ontology assignments indicate that sequences obtained cover a broad range of GO categories.
We describe an efficient method to identify thousands of SNPs distributed throughout the alfalfa genome covering a broad range of GO categories. Validated SNPs represent valuable molecular marker resources that can be used to enhance marker density in linkage maps, identify potential factors involved in heterosis and genetic variation, and as tools for association mapping and genomic selection in alfalfa.
DNA pyrosequencing using 454 sequencing technology enables sequencing millions of high-quality DNA bases per sequencing run [1, 2]. This method has been used successfully for transcriptome sequencing and identification of single nucleotide polymorphisms (SNP) in many plant species including maize (Zea mays L.) , sugarcane (Saccarum spp.) , eucalyptus (Eucalyptus grandis) , and the model legume and close relative of alfalfa, M. truncatula. The M. truncatula genome is one of three legume genomes including Lotus japonicus and soybean (Glycine max) that have been assembled and annotated [7, 8]. SNPs are the most common sequence variation among plants and are often functionally important. SNPs can be converted into genetic markers that can be inexpensively assayed using high-throughput approaches . Traditionally, molecular markers have been used to determine genetic relatedness between plant materials, to assist in the identification of novel sources of genetic variation, to study evolutionary relationships, to confirm the pedigree and identity of new varieties, in population structure analysis and association genetics, to locate quantitative trait loci (QTLs) and genes of interest, and for marker-assisted breeding . The value and uses of DNA markers have been shaped in large part by innovations in marker technologies that increase throughput and reduce costs per data point . A large number of molecular markers can be used in high-throughput genotyping platforms in association mapping studies to dissect complex traits and in molecular breeding approaches at the whole genome level .
Alfalfa (Medicago sativa L.) is one of the most important forage legume species worldwide and the third most valuable crop in the USA ($8 billion per annum). Alfalfa is a high yielding perennial species that requires little or no nitrogen fertilizer because of its ability to carry out symbiotic nitrogen fixation and can be harvested multiple times during the growing season. Therefore, it has been an important component of sustainable agricultural systems for many years and has recently been promoted as a potential bioenergy crop . Cultivated alfalfa is tetraploid (2n = 4x = 32) and displays tetrasomic inheritance . Alfalfa is partially self-incompatible and populations are extremely polymorphic due to their high degree of outcrossing. Inbreeding severely depresses plant vigor and fertility in tetraploid alfalfa due to the loss of complementary gene interactions [15, 16], preventing the development of inbred lines. A significant level of sequence conservation was reported between alfalfa and M. truncatula allowing estimates of marker colinearity between the two species . Simple sequence repeat (SSR) markers developed from M. truncatula and some alfalfa sequences are currently available [19–21]. Alfalfa like many other crop species, lacks validated SNP markers which are required in large numbers for map-based gene isolation, association genetics and genomic selection approaches .
High-resolution melting (HRM) curve analysis has proven to be a highly sensitive method for mutation discovery and SNP genotyping . Nucleic acid melting is tracked by monitoring the fluorescence of the samples across a defined temperature range generating high-resolution melting profiles that are used to identify the presence of sequence variation within the amplicon . Base-pair mismatches shift the stability of a duplex by varying amounts depending on the particular mismatch, the mismatch position, and neighboring base pairs . HRM was successfully implemented to assay SNP variation in diploid and tetraploid alfalfa .
We describe the generation of thousands of alfalfa ESTs (expressed sequence tags) obtained using 454 sequencing technology and the first iteration of genome-wide SNP identification in tetraploid alfalfa, a non-model plant species that currently lacks a sequenced genome. A sample of candidate SNPs were experimentally validated using a HRM platform.
Sequence assembly and annotation
Alfalfa sequence length distribution pre-and post-assembly
Total base count (bp)
Minimum length (bp)
Maximum length (bp)
Average length (bp)
Total base count (bp)
> 2000 bp
Minimum length (bp)
Maximum length (bp)
Average length (bp)
Alignment of alfalfa sequences to the M. truncatula genome
Summary of candidate SNP containing assemblies and their location on the MT3.5 chromosomes
TC with SNPs
Number of SNPs
Mapped to Mt. 3.5
When combined with the MtGI results, 46,086 (85%) alfalfa sequences had a blastn or blastx hit in MtGI and/or MT3.5 databases. Of the remaining 8,130 alfalfa sequences, 992 of them had a blastx hit (cut off E value 1e-5) with 763 proteins in the soybean Glyma1.0 high confidence gene protein database  and 732 alfalfa sequences had a blastx hit (cut off E value 1e-5) with 572 L. japonicus proteins  further supporting the notion of shared sequences among legumes. Only 12 of the remaining 7,057 alfalfa sequences without a match in any of the aforementioned databases had a hit in the Arabidopsis TAIR 9.0 protein database (cut off E value 1e-5). The remaining 7,045 sequences represent a total length of 2.3 Mb (8% of the total sequences post-assembly), did not have a significant match when blasted against the NCBI non-redundant nucleotide database suggesting they may be miss-assembled or novel sequences.
SNP discovery and validation
SNP markers have a broad range of applications and a large number of markers are needed to implement association genetics approaches and genomic selection in plants with low levels of LD [12, 29]. Large scale SNP discovery efforts in a wide range of crop plants are currently limited to a few species. Previous SNP discovery approaches in other species include sequence comparisons between the Arabidopsis ecotypes Landsberg and Columbia  and subspecies of rice . Although alfalfa is the most widely cultivated forage legume, it lacks genomics resource development compared to other crop species, partly because it is an autotetraploid, perennial outcrossing species. Thus, a whole genome SNP discovery effort is among the most critical endeavors to develop the necessary tools and integrate molecular breeding approaches in alfalfa. In this study, we utilized 454 sequencing for genome-wide SNP discovery between and within two genotypes from M. sativa subsp. sativa (Chilean) and M. sativa subsp. falcata (Wisfal), which have contrasting phenotypes for water stress tolerance. Sequencing and SNP discovery based on transcriptome sequencing of water stressed plants increased the likelihood of having sequence coverage in genes relevant to water-stress and identification of allelic differences in these genes between the two alfalfa genotypes evaluated in this study.
Sequence assembly and alignment to M. truncatula
The close phylogenetic relationship between M. truncatula and alfalfa, same basic chromosome number (x = 8), high degree of sequence similarity and marker co-linearity [17, 18] facilitated using the M. truncatula genome sequence as a scaffold to align alfalfa transcriptome sequences. Despite the size and complexity of the alfalfa genome, alignments of alfalfa 454 sequences throughout the M. truncatula genome were apparent (Figure 2). The M. truncatula genome sequence project targeted the euchromatic region and available data indicate that 66% of the genes captured in the current pseudomolecules and 24% of the genes in M. truncatula captured in the fragmented Illumina assemblies are found within the euchromatin. Therefore, the M. truncatula genome provides a valuable resource for gene annotation of alfalfa sequences. The identification of some alfalfa sequences without alignment to the M. truncatula genome is not unexpected considering that M. truncatula is a selfing, annual diploid species and alfalfa is an outcrossing, perennial tetraploid species with a considerably larger genome. Additionally, non-alignment of alfalfa sequences to the M. truncatula genome could occur in regions where there have been authentic genome rearrangements between the two species. We demonstrated the feasibility of using genomic resources from model species to quickly assemble transcriptome sequences from related species. These alfalfa resources can also facilitate opportunities to further understand legume genome structure and evolution.
SNP discovery and validation
454 sequencing is a high-throughput approach to identify sequence variation and was utilized successfully here to identify SNP variation in tetraploid alfalfa. Alfalfa is genetically very diverse and the existence of heterosis between M. sativa subsp. sativa and M. sativa subsp. falcata was previously reported [32, 33]. Because Chilean and Wisfal are highly heterozygous outcrossing lines, transcriptome sequencing of alfalfa genotypes enabled the identification of gene-associated SNPs within and between genotypes covering a range of gene functions and location in the genome. Reports that 454 sequencing technology may be susceptible to indel-type errors  prompted us to consider only base substitutions (i.e. SNPs) with stringent SNP prediction criteria.
The HRM platform proved to be a suitable method for validation and genotyping of a large number of candidate SNPs. The HRM melting curve profiles enables identification of homozygous and heterozygous SNP genotypes . We report a 62.5% rate of validation of single candidate SNPs (120 out of the 192) indicating that transcriptome sequencing of a crop legume with economic value and assembly of those transcript sequences using the genome of a closely related model species is a viable strategy for genome-wide SNP discovery efforts. Studies in eucalyptus reported a 83% rate of validation from the candidate SNPs . Additional SNP in alfalfa were detected using HRM that were not predicted based on the bioinformatics analysis from the 454 sequences. These findings were partly due to our efforts to minimize the false discovery rate in sequences with low depth of sequence coverage. Therefore, the challenge for SNP discovery using 454 sequences in a highly heterozygous species such as alfalfa resides in limitations from the depth of sequence coverage rather than sequencing errors. In contrast to a diploid species, tetraploid species have four homologous chromosomes and therefore deeper sequence coverage is needed to accurately ascertain the SNP variation within and among genotypes. Although an increase in the rate of SNP prediction can be obtained using less stringent SNP parameters, it also increases the likelihood of false positive SNP calls. The discrepancy between SNP detection and validation could also be the result of preferential allele expression or allele silencing in the pooled cDNA sample that was used because the initial SNP identification was based on cDNA sequencing whereas validation was made using PCR fragments amplified from genomic DNA. The identified SNPs complement our current work identifying genomic regions associated with water stress tolerance in backcross populations developed from the same two genotypes used for 454 sequencing.
Previous efforts to identify sequence polymorphism in alfalfa include array hybridizations that identify single feature polymorphisms (SFP) using M. truncatula Affy Chips . The approach described here enables the discovery of SNPs within genes that may not be present on an array developed for the model M. truncatula. Additionally, SNP discovery efforts based on transcriptome sequencing in a crop species may result in the identification of SNPs with functional relevance to agriculturally important traits such as persistence, likely not present in an annual model system.
The large number of candidate SNPs with potential functional relevance identified in this study provides a valuable tool for comparative genomic studies between M. truncatula and alfalfa, and potentially to other forage legume species including white and red clover. Evaluations of SNP variability, frequency, and distribution, as well as estimates of linkage disequilibrium in alfalfa are now attainable research targets. Applications of crop-related SNP discovery efforts include genetic diversity analysis, increasing the marker density of genetic maps to support high-resolution association and linkage mapping of target traits, facilitation of future genome assembly of this highly complex genome, integration of genetic and physical maps, and increasing the feasibility of integrating marker-dense molecular breeding approaches such as genomic selection in this important forage legume.
In order to enhance crop productivity, we need to improve our understanding of the genes involved in plant growth, development, and adaptation to a changing environment including suboptimal soil conditions and biotic and abiotic stress conditions. The 454 EST-based SNP discovery and validation pipeline described here focused on two tetraploid alfalfa subspecies with contrasting water stress tolerance provides a platform for the discovery of agriculturally-relevant genes and corresponding SNPs. The availability of a large number of validated SNPs is a valuable molecular marker resource for genetics, high-resolution linkage and association mapping, and molecular breeding approaches using genomic selection, for practical and agriculturally-relevant improvements in alfalfa which are not feasible using model systems.
Two tetraploid alfalfa genotypes (sampled from M. sativa subsp. sativa var. 'Chilean' and M. sativa subsp. falcata var. 'Wisfal'  were clonally propagated and grown in metromix and sand mixture (4:1 v:v). Plant tissue was sampled at seven time points during a progressively developing water deficit including: well-watered, three, seven, nine, 11 and 14 days after withholding water, and one-day re-watering on day 15. Roots and shoots were sampled separately. Two to four clonal propagules of each genotype were sampled at each time point. For each genotype, cDNA was prepared from shoot and root tissue of each clonally propagated sample grown under either well-watered or water-stressed conditions, and then equal amounts of cDNA from each sampling point were pooled for sequencing using 454 sequencing technology.
cDNA preparation and sequencing
Total RNA was extracted using TRIzol® (Invitrogen, Carlsbad, CA) following the manufacturer's recommendations, digested with DNAseI (Ambion, Austin, TX), column purified with RNeasy MinElute Cleanup kit (Qiagen, Valencia, CA) and quantified using a NanoDrop Spectrophotometer ND-100 (NanoDrop Technologies, Willington, DE). Total RNA integrity was assessed using a Bioanalyzer 2100 (Agilent, Santa Clara, CA). Messenger RNA (mRNA) isolation was performed using Poly(A)Purist(tm) MAG (Magnetic mRNA Purification Kit, Ambion, Austin, TX). For the first strand cDNA synthesis, three microgram of mRNA was used with Oligo (dT)12-18 and Superscript II Reverse Transcriptase followed by a second cDNA synthesis (dscDNA) using Superscript First and Double Stranded cDNA synthesis Kit (Invitrogen, Carlsbad, CA). Double stranded cDNA from roots and shoots from the same genotype was pooled and deep sequenced using Titanium 454 technology at the J. Craig Venter Institute (JCVI) with each genotype occupying half of a 454 sequencing plate.
The raw sequences were assembled for each genotype separately and then combined (Table 1) using TGICL clustering utility  with The Paracel Transcript Assembler instead of CAP3 and a modified version of the TIGR Gene Index pipeline. The blast searches were performed using TCs and singletons from the assembly with the two genotypes combined (Table 1). All TCs and singletons were blasted against the MtGI 10.0 database using Blastn with e-10 as a cut-off expected value. The TCs and singletons were assigned to the GO categories based on their top hit.
SNP discovery and validation
Candidate SNP calls were predicted using Perl scripts to parse the transcript assembly files and requiring a minimum coverage of at least two 454 reads for each form of the SNP at the variant position. A confidence score was generated based on the percent of nucleotides matching the surrounding regions of the location of the SNP. Low confidence regions with less than 90% identity where a SNP was called were not included in the final SNP predictions. SNPs were quality controlled by remapping the 454 reads back onto the assembly and checked against the ACE alignment files. Indels were ignored and the process was biased towards getting good quality SNPs vs. getting all the SNPs at the expense of false SNP calls. Primer3 was used to design primers targeting candidate SNPs. Criteria for primer design include a predicted annealing temperature (Tm) of 59°C to 61°C, primer length ranging between 18-24 bp and PCR amplicon lengths of 40 to 200 bp. All PCR reactions were performed in 384-well plates using a 9700 Thermal Cycler (Applied Biosystems, Foster City, CA, USA) using a total volume of 5 mL per well. The PCR reaction mixture consisted of 5 ng of genomic DNA, 0.25 μM of forward and reverse primer, 1X LightScanner High Sensitivity Master Mix (Idaho Technologies, Salt Lake, UT, USA) and 10 mL mineral oil. After an initial denaturation step of 2 min at 95°C, 45 PCR cycles were performed with 30 s of denaturation at 94°C and 30 s at the target annealing temperature, followed by a final hold at 4°C. Samples were then transferred to a LightScanner 384-well system (Idaho Technologies, Salt Lake, UT) and a melting cycle was performed by increasing the temperature at 0.1°C s-1 from 56 to 95°C. Melting data was analyzed and visualized with the LightScanner Software with CALL-IT 2.0 (Idaho Technologies, Salt Lake, UT) using the small amplicon module.
Re-sequencing PCR amplicons
The PCR amplicons from twenty candidate SNP primers were cleaned using Exo/SAP-IT PCR product cleanup kit following the manufacturer's instructions (Affymetrix, Inc., Cleveland, OH) and directly sequenced using a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA) using both forward and reverse primers. The amplicons were analyzed using an ABI 3730 automated sequencer (PE-ABI, Foster City, CA). The sequence data was analyzed and aligned for SNP discovery with Sequencher 4.8 (Gene Codes).
This study was supported by The Samuel Roberts Noble Foundation, grants from the Oklahoma Bioenergy Center to M.J. Monteros and M.K. Udvardi, and a grant from Forage Genetics International to M.J. Monteros. We would like to thank Shelby Bidwell for uploading the 454 sequences to Genebank and assisting in sequence assembly, Cook Ly for assistance with preparing figures, and Mark Sorrells for providing comments on the manuscript.
- Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer MLI, Jarvie TP, Jirage KB, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, Lohman KL, Lu H, Makhijani VB, et al: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.PubMed CentralPubMedGoogle Scholar
- Barbazuk WB, Schnable PS: SNP discovery by transcriptome pyrosequencing. Methods Mol Biol. 2011, 729: 225-246. 10.1007/978-1-61779-065-2_15.PubMedView ArticleGoogle Scholar
- Barbazuk WB, Scott JE, Hsin DC, Li L, Patrick SS: SNP discovery via 454 transcriptome sequencing. Plant J. 2007, 51: 910-918. 10.1111/j.1365-313X.2007.03193.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Bundock PC, Eliott FG, Ablett G, Benson AD, Casu RE, Aitken KS, Henry RJ: Targeted single nucleotide polymorphism (SNP) discovery in a highly polyploid plant species using 454 sequencing. Plant Biotechnol J. 2009, 7: 347-354. 10.1111/j.1467-7652.2009.00401.x.PubMedView ArticleGoogle Scholar
- Novaes E, Drost D, Farmerie W, Pappas G, Grattapaglia D, Sederoff R, Kirst M: High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome. BMC Genomics. 2008, 9: 312-10.1186/1471-2164-9-312.PubMed CentralPubMedView ArticleGoogle Scholar
- Cheung F, Haas B, Goldberg S, May G, Xiao Y, Town C: Sequencing Medicago truncatula expressed sequenced tags using 454 Life Sciences technology. BMC Genomics. 2006, 7: 272-10.1186/1471-2164-7-272.PubMed CentralPubMedView ArticleGoogle Scholar
- Young ND, Udvardi M: Translating Medicago truncatula genomics to crop legumes. Curr Opin Plant Biol. 2009, 12: 193-201. 10.1016/j.pbi.2008.11.005.PubMedView ArticleGoogle Scholar
- Cannon SB, May GD, Jackson SA: Three sequenced legume genomes and many crop species: rich opportunities for translational genomics. Plant Physiol. 2009, 151: 970-977. 10.1104/pp.109.144659.PubMed CentralPubMedView ArticleGoogle Scholar
- Kwok PY: Methods for genotyping single nucleotide polymorphisms. Annu Rev Genomics Hum Genet. 2001, 2: 235-258. 10.1146/annurev.genom.2.1.235.PubMedView ArticleGoogle Scholar
- Ganal MW, Altmann T, Röder MS: SNP identification in crop plants. Curr Opin Plant Biol. 2009, 12: 211-217. 10.1016/j.pbi.2008.12.009.PubMedView ArticleGoogle Scholar
- Walker DR, Monteros MJ, Yates JL: Chapter 6: Molecular breeding. Genetics, Genomics and Breeding of Crop Plants (Soybean). Edited by: Bilyeu K, Ratnaparkhe M, Kole C. 2010, Science Publishers IncGoogle Scholar
- Heffner EL, Sorrells MK, Jannink JL: Genomic selection for crop improvement. Crop Sci. 2009, 49: 1-12. 10.2135/cropsci2008.08.0512.View ArticleGoogle Scholar
- Ebert J: Alfalfa's bioenergy appeal. Ethanol Producer Magazine. 2007, BBI International, 88-94.Google Scholar
- McCoy TJ, Bingham ET: Cytology and cytogenetics of alfalfa. Alfalfa and alfalfa improvement. Edited by: Hanson AA, Barnes, DK, Hill, RR. 1988, Madison, Wisconsin: American Society of Agronomy, 29: 737-776.Google Scholar
- Bingham ET, Groose RW, Woodfield DR, Kidwell KK: Complementary gene interactions in alfalfa are greater in autotetraploids than diploids. Crop Sci. 1994, 34: 823-829. 10.2135/cropsci1994.0011183X003400040001x.View ArticleGoogle Scholar
- Osborn TC, Brouwer D, McCoy TJ: Molecular marker analysis in alfalfa. Biotechnology and the improvement of forage legumes. Edited by: McKersie BD, Brown DCW. 1997, Guelph, Canada: CAB International, 91-109.Google Scholar
- Zhu H, Choi HK, Cook DR, Shoemaker RC: Bridging model and crop legumes through comparative genomics. Plant Physiol. 2005, 137: 1189-1196. 10.1104/pp.104.058891.PubMed CentralPubMedView ArticleGoogle Scholar
- Choi HK, Kim D, Uhm T, Limpens E, Lim H, Mun JH, Kalo P, Penmetsa RV, Seres A, Kulikova O, Roe BA, Bisseling T, Kiss GB, Cook DR: A sequence-based genetic map of Medicago truncatula and comparison of marker colinearity with M. sativa. Genetics. 2004, 166: 1463-1502. 10.1534/genetics.166.3.1463.PubMed CentralPubMedView ArticleGoogle Scholar
- Diwan N, Bouton JH, Kochert G, Cregan PB: Mapping of simple sequence repeat (SSR) DNA markers in diploid and tetraploid alfalfa. Theor Appl Genet. 2000, 101: 165-172. 10.1007/s001220051465.View ArticleGoogle Scholar
- Eujayl I, Sledge M, Wang L, May G, Chekhovskiy K, Zwonitzer J, Mian M: Medicago truncatula EST-SSRs reveal cross-species genetic markers for Medicago spp. Theor Appl Genet. 2004, 108: 414-422. 10.1007/s00122-003-1450-6.PubMedView ArticleGoogle Scholar
- Flajoulot S, Ronfort J, Baudouin P, Barre P, Huguet T, Huyghe C, Julier B: Genetic diversity among alfalfa (Medicago sativa) cultivars coming from a breeding program, using SSR markers. Theor Appl Genet. 2005, 111: 1420-1429. 10.1007/s00122-005-0074-4.PubMedView ArticleGoogle Scholar
- Cho MH, Ciulla D, Klanderman BJ, Raby BA, Silverman EK: High-resolution melting curve analysis of genomic and whole-genome amplified DNA. Clin Chem. 2008, 54: 2055-2058. 10.1373/clinchem.2008.109744.PubMed CentralPubMedView ArticleGoogle Scholar
- deSilva D, Blackett J: Assay: High-resolution melting and unlabeled probes. Gen Eng & Biotech News. 2007, 27:Google Scholar
- SantaLucia J, Allawi HT, Seneviratne PA: Improved nearest-neighbor parameters for predicting DNA duplex stability. Biochemistry. 1996, 35: 3555-3562. 10.1021/bi951907q.PubMedView ArticleGoogle Scholar
- Han Y, Khu D, Monteros M: High-resolution melting analysis for SNP genotyping in diploid and tetraploid alfalfa (Medicago sativa L.). Mol Breeding. 2011, 10.1007/s11032-011-9566-x.Google Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.PubMedView ArticleGoogle Scholar
- Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L: Genome sequence of the palaeopolyploid soybean. Nature. 463: 178-183.
- Sato S, Nakamura Y, Kaneko T, Asamizu E, Kato T, Nakao M, Sasamoto S, Watanabe A, Ono A, Kawashima K, Fujishiro T, Katoh M, Kohara M, Kishida Y, Minami C, Nakayama S, Nakazaki N, Shimizu Y, Shinpo S, Takahashi C, Wada T, Yamada M, Ohmido N, Hayashi M, Fukui K, Baba T, Nakamichi T, Mori H, Tabata S: Genome structure of the legume, Lotus japonicus. DNA Res. 2008, 1-13.Google Scholar
- Yu J, Buckler ES: Genetic association mapping and genome organization of maize. Curr Opin Biotechnol. 2006, 17: 155-160. 10.1016/j.copbio.2006.02.003.PubMedView ArticleGoogle Scholar
- Jander G, Norris SR, Rounsley SD, Bush DF, Levin IM, Last RL: Arabidopsis map-based cloning in the post-genome era. Plant Physiol. 2002, 129: 440-450. 10.1104/pp.003533.PubMed CentralPubMedView ArticleGoogle Scholar
- Feltus AF, Wan J, Schulze SR, Estill JC, Jiang N, Paterson AH: An SNP resource for rice genetics and breeding based on subspecies indica and japonica genome alignments. Genome Res. 2004, 14: 1812-1819. 10.1101/gr.2479404.PubMed CentralPubMedView ArticleGoogle Scholar
- Riday H, Brummer EC: Forage yield heterosis in alfalfa. Crop Sci. 2002, 42: 716-723. 10.2135/cropsci2002.0716.View ArticleGoogle Scholar
- Robins JG, Luth D, Campbell TA, Bauchan GR, He C, Viands DR, Hansen JL, Brummer EC: Genetic mapping of biomass production in tetraploid alfalfa. Crop Sci. 2007, 47: 1-10. 10.2135/cropsci2005.11.0401.View ArticleGoogle Scholar
- Montgomery J, Wittwer CT, Palais R, Zhou L: Simultaneous mutation scanning and genotyping by high-resolution DNA melting analysis. Nat Protocols. 2007, 2: 59-66.View ArticleGoogle Scholar
- Yang SS, Xu WW, Tesfaye M, Lamb JFS, Jung HJG, Samac DA, Vance CP, Gronwald JW: Single-feature polymorphism discovery in the transcriptome of tetraploid alfalfa. Plant Genome. 2009, 2: 224-232. 10.3835/plantgenome2009.03.0014.View ArticleGoogle Scholar
- Sledge M, Ray I, Jiang G: An expressed sequence tag SSR map of tetraploid alfalfa (Medicago sativa L.). Theor Appl Genet. 2005, 111: 980-992. 10.1007/s00122-005-0038-8.PubMedView ArticleGoogle Scholar