Optimization of linkage mapping strategy and construction of a high-density American lotus linkage map

Background

Lotus is a diploid plant with agricultural, medicinal, and ecological significance. Genetic linkage maps are fundamental resources for genome and genetic study, and also provide molecular markers for breeding in agriculturally important species. Genotyping by sequencing revolutionized genetic mapping, the restriction-site associated DNA sequencing (RADseq) allowed rapid discovery of thousands of SNPs markers, and a crucial aspect of the sequence based mapping strategy is the reference sequences used for marker identification.

Results

We assessed the effectiveness of linkage mapping using three types of references for scoring markers: the unmasked genome, repeat masked genome, and gene models. Overall, the repeat masked genome produced the optimal genetic maps. A high-density genetic map of American lotus was constructed using an F₁ population derived from a cross between Nelumbo nucifera ‘China Antique’ and N. lutea ‘AL1’. A total of 4,098 RADseq markers were used to construct the American lotus ‘AL1’ genetic map, and 147 markers were used to construct the Chinese lotus ‘China Antique’ genetic map. The American lotus map has 9 linkage groups, and spans 494.3 cM, with an average distance of 0.7 cM between adjacent markers. The American lotus map was used to anchor scaffold sequences in the N. nucifera ‘China Antique’ draft genome. 3,603 RADseq markers anchored 234 individual scaffold sequences into 9 megascaffolds spanning 67% of the 804 Mb draft genome.

Conclusions

Among the unmasked genome, repeat masked genome and gene models, the optimal reference sequences to call RADseq markers for map construction is repeat masked genome. This high density genetic map is a valuable resource for genomic research and crop improvement in lotus.     

Construction of High Density Sweet Cherry (Prunus avium L.) Linkage Maps Using Microsatellite Markers and SNPs Detected by Genotyping-by-Sequencing (GBS)

Linkage maps are valuable tools in genetic and genomic studies. For sweet cherry, linkage maps have been constructed using mainly microsatellite markers (SSRs) and, recently, using single nucleotide polymorphism markers (SNPs) from a cherry 6K SNP array. Genotyping-by-sequencing (GBS), a new methodology based on high-throughput sequencing, holds great promise for identification of high number of SNPs and construction of high density linkage maps. In this study, GBS was used to identify SNPs from an intra-specific sweet cherry cross. A total of 8,476 high quality SNPs were selected for mapping. The physical position for each SNP was determined using the peach genome, Peach v1.0, as reference, and a homogeneous distribution of markers along the eight peach scaffolds was obtained. On average, 65.6% of the SNPs were present in genic regions and 49.8% were located in exonic regions. In addition to the SNPs, a group of SSRs was also used for construction of linkage maps. Parental and consensus high density maps were constructed by genotyping 166 siblings from a ‘Rainier’ x ‘Rivedel’ (Ra x Ri) cross. Using Ra x Ri population, 462, 489 and 985 markers were mapped into eight linkage groups in ‘Rainier’, ‘Rivedel’ and the Ra x Ri map, respectively, with 80% of mapped SNPs located in genic regions. Obtained maps spanned 549.5, 582.6 and 731.3 cM for ‘Rainier’, ‘Rivedel’ and consensus maps, respectively, with an average distance of 1.2 cM between adjacent markers for both ‘Rainier’ and ‘Rivedel’ maps and of 0.7 cM for Ra x Ri map. High synteny and co-linearity was observed between obtained maps and with Peach v1.0. These new high density linkage maps provide valuable information on the sweet cherry genome, and serve as the basis for identification of QTLs and genes relevant for the breeding of the species.     

Isolation of EST-derived microsatellite markers for genotyping the A and B genomes of wheat

Genetic variation present in 64 durum wheat accessions was investigated by using three sources of microsatellite (SSR) markers: EST-derived SSRs (EST-SSRs) and two sources of SSRs isolated from total genomic DNA. Out of 245 SSR primer pairs screened, 22 EST-SSRs and 20 genomic-derived SSRs were polymorphic and used for genotyping. The EST-SSR primers produced high quality markers, but had the lowest level of polymorphism (25%) compared to the other two sources of genomic SSR markers (53%). The 42 SSR markers detected 189 polymorphic alleles with an average number of 4.5 alleles per locus. The coefficient of similarity ranged from 0.28 to 0.70 and the estimates of similarity varied when different sources of SSR markers were used to genotype the accessions. This study showed that EST-derived SSR markers developed in bread wheat are polymorphic in durum wheat when assaying loci of the A and B genomes. A minumum of ten EST-SSRs generated a very low probability of identity (0.36×10⁻¹²) indicating that these SSRs have a very high discriminatory power. EST-SSR markers directly sample variation in transcribed regions of the genome, which may enhance their value in marker-assisted selection, comparative genetic analysis and for exploiting wheat genetic resources by providing a more-direct estimate of functional diversity. Received: 19 December 2000 / Accepted: 17 April 2001     

High-density Integrated Linkage Map Based on SSR Markers in Soybean

A well-saturated molecular linkage map is a prerequisite for modern plant breeding. Several genetic maps have been developed for soybean with various types of molecular markers. Simple sequence repeats (SSRs) are single-locus markers with high allelic variation and are widely applicable to different genotypes. We have now mapped 1810 SSR or sequence-tagged site markers in one or more of three recombinant inbred populations of soybean (the US cultivar ‘Jack’ × the Japanese cultivar ‘Fukuyutaka’, the Chinese cultivar ‘Peking’ × the Japanese cultivar ‘Akita’, and the Japanese cultivar ‘Misuzudaizu’ × the Chinese breeding line ‘Moshidou Gong 503’) and have aligned these markers with the 20 consensus linkage groups (LGs). The total length of the integrated linkage map was 2442.9 cM, and the average number of molecular markers was 90.5 (range of 70–114) for the 20 LGs. We examined allelic diversity for 1238 of the SSR markers among 23 soybean cultivars or lines and a wild accession. The number of alleles per locus ranged from 2 to 7, with an average of 2.8. Our high-density linkage map should facilitate ongoing and future genomic research such as analysis of quantitative trait loci and positional cloning in addition to marker-assisted selection in soybean breeding.Key words: EST-derived SSR marker, integrated linkage map, microsatellite marker, polymorphism information content     

Comparative analysis of genetic diversity in sacred lotus (<Emphasis Type="Italic">Nelumbo nucifera</Emphasis> Gaertn.) using AFLP and SSR markers

The sacred lotus (Nelumbo nucifera Gaertn.) is an aquatic plant of economic and ornamental importance in China. In this study, we developed twenty novel sacred lotus SSR markers, and used AFLP and SSR markers to investigate the genetic diversity and genetic relationships among 58 accessions of N. nucifera including 15 seed lotus, 12 rhizome lotus, 24 flower lotus and 7 wild lotus. Our results showed that sacred lotus exhibited a low level of genetic diversity, which may attribute to asexual reproduction and long-term artificial selection. A dendrogram based on both AFLP and SSR clustering data showed that: (1) the seed lotus accessions and rhizome lotus accessions were distinctly clustered into different groups, which indicated the significant genetic differentiation between them. This may be attributed to the two modes of reproduction and lack of genetic exchange; (2) the accessions of Thailand wild lotus were separated from other wild lotus accessions. This implied that the Thailand lotus might be genetically differentiated from other wild lotuses. In addition, Mantel test conducted gave highly significant correlation between AFLP-SSR data and each of the AFLP and SSR ones, with the values of r = 0.941 and r = 0.879, respectively, indicating the higher efficiency of the combination of these techniques (AFLP and SSR) in estimation and validation of the genetic diversity among the accession of sacred lotus. This knowledge of the genetic diversity and genetic relatedness of N. nucifera is potentially useful to improve the current strategies in breeding and germplasm conservation to enhance the ornamental and economic value of sacred lotus.     

Genome-Wide SNP Identification and Characterization in Two Soybean Cultivars with Contrasting Mungbean Yellow Mosaic India Virus Disease Resistance Traits

Mungbean yellow mosaic India virus (MYMIV) is a bipartite Geminivirus, which causes severe yield loss in soybean (Glycine max). Considering this, the present study was conducted to develop large-scale genome-wide single nucleotide polymorphism (SNP) markers and identify potential markers linked with known disease resistance loci for their effective use in genomics-assisted breeding to impart durable MYMIV tolerance. The whole-genome re-sequencing of MYMIV resistant cultivar ‘UPSM-534’ and susceptible Indian cultivar ‘JS-335’ was performed to identify high-quality SNPs and InDels (insertion and deletions). Approximately 234 and 255 million of 100-bp paired-end reads were generated from UPSM-534 and JS-335, respectively, which provided ~98% coverage of reference soybean genome. A total of 3083987 SNPs (1559556 in UPSM-534 and 1524431 in JS-335) and 562858 InDels (281958 in UPSM-534 and 280900 in JS-335) were identified. Of these, 1514 SNPs were found to be present in 564 candidate disease resistance genes. Among these, 829 non-synonymous and 671 synonymous SNPs were detected in 266 and 286 defence-related genes, respectively. Noteworthy, a non-synonymous SNP (in chromosome 18, named 18-1861613) at the 149^th base-pair of LEUCINE-RICH REPEAT RECEPTOR-LIKE PROTEIN KINASE gene responsible for a G/C transversion [proline (CCC) to alanine(GCC)] was identified and validated in a set of 12 soybean cultivars. Taken together, the present study generated a large-scale genomic resource such as, SNPs and InDels at a genome-wide scale that will facilitate the dissection of various complex traits through construction of high-density linkage maps and fine mapping. In the present scenario, these markers can be effectively used to design high-density SNP arrays for their large-scale validation and high-throughput genotyping in diverse natural and mapping populations, which could accelerate genomics-assisted MYMIV disease resistance breeding in soybean.     

Development and Characterization of Polymorphic EST-SSR and Genomic SSR Markers for Tibetan Annual Wild Barley

Tibetan annual wild barley is rich in genetic variation. This study was aimed at the exploitation of new SSRs for the genetic diversity and phylogenetic analysis of wild barley by data mining. We developed 49 novel EST-SSRs and confirmed 20 genomic SSRs for 80 Tibetan annual wild barley and 16 cultivated barley accessions. A total of 213 alleles were generated from 69 loci with an average of 3.14 alleles per locus. The trimeric repeats were the most abundant motifs (40.82%) among the EST-SSRs, while the majority of the genomic SSRs were di-nuleotide repeats. The polymorphic information content (PIC) ranged from 0.08 to 0.75 with a mean of 0.46. Besides this, the expected heterozygosity (He) ranged from 0.0854 to 0.7842 with an average of 0.5279. Overall, the polymorphism of genomic SSRs was higher than that of EST-SSRs. Furthermore, the number of alleles and the PIC of wild barley were both higher than that of cultivated barley, being 3.12 vs 2.59 and 0.44 vs 0.37. Indicating more polymorphism existed in the Tibetan wild barley than in cultivated barley. The 96 accessions were divided into eight subpopulations based on 69 SSR markers, and the cultivated genotypes can be clearly separated from wild barleys. A total of 47 SSR-containing EST unigenes showed significant similarities to the known genes. These EST-SSR markers have potential for application in germplasm appraisal, genetic diversity and population structure analysis, facilitating marker-assisted breeding and crop improvement in barley.     

Development of EST-SSR markers and their utility in revealing cryptic diversity in safflower (Carthamus tinctorius L.)

With the aim of developing additional genomic resources in safflower, a set of 41,011 ESTs of safflower were mined for the presence of SSRs. 18,773 SSR containing ESTs (SSR-ESTs) were identified and were analyzed to remove redundant sequences leading to identification of 8,810 non-redundant SSR-ESTs (categorized into 6104 singletons and 2,706 contigs) having 13,085 non-redundant SSRs. The average number of non-redundant SSRs per EST was 0.32 and they predominantly consisted of dinucleotide (57.7 %), and trinucleotide (37.7 %) repeat motifs. 500 primer pairs were designed for the non-redundant EST-SSRs of which, 151 were tested. 60 markers which gave robust amplicons, were validated in a set of 19 Carthamus lines. A subset of EST-SSR markers, having average polymorphism information content (PIC) ≥0.4 could precisely elucidate the pedigree relatedness among these lines. Further, these markers exhibited high cross-species transferability to five other wild species of Carthamus. The markers reported here would be a valuable addition to existing safflower marker resources aiding in hastening its improvement.     

A database of simple sequence repeats from cereal and legume expressed sequence tags mined in silico: survey and evaluation

Simple sequence repeats (SSRs) or microsatellites are an important class of molecular markers for genome analysis and plant breeding applications. In this paper, the SSR distributions within ESTs from the legumes soybean (Glycine max, representing 135.86 Mb), medicago (Medicago truncatula, 121.1 Mb) and lotus (Lotus japonicus, 45.4 Mb) have been studied relative to the distributions in cereals such as sorghum (Sorghum bicolor, 98.9 Mb), rice (Oryza sativa, 143.9 Mb) and maize (Zea mays, 183.7 Mb). The relative abundance, density, composition and putative annotations of di-, tri-, tetra- and penta-nucleotide repeats have been compared and SSR containing ESTs (SSR-ESTs) have been clustered to give a non-redundant set of EST-SSRs, available in a database. Further, a subset of such candidate EST-SSRs from sorghum have been tested for their ability to detect polymorphism between Striga-susceptible, stay-green drought tolerant mapping population parent 'E 36-1' and its Striga-resistant, non-stay-green counterpart 'N13'. Primer sets for 64% of the EST-SSRs tested produced a clear and specific PCR product band and 34% of these detected scorable polymorphism between the N13 and E 36-1 parental lines. Over half of these markers have been genotyped on 94 RILs from the (N13 x E 36-1)-based mapping population, with 42 markers mapping onto the ten sorghum linkage groups. This establishes the value of this database as a resource of molecular markers for practical applications in cereal and legume genetics and breeding. The primer pairs for non-redundant EST-SSRs have been designed and are freely available through the database (http://intranet.icrisat.org/gt1/ssr/ssrdatabase.html).     

Utility of EST-derived SSRs as population genetics markers in a beetle

Microsatellite, or simple sequence repeat (SSR), loci can be identified by mining expressed sequence tag (EST) databases, and where these are available, marker development time and expense can be decreased considerably over conventional strategies of probing the entire genome. However, it is unclear whether they provide information on population structure similar to that generated by anonymous genomic SSRs. We performed comparative population genetic analyses between EST-derived SSRs (EST-SSRs) and anonymous SSRs developed from genomic DNA for the same set of populations of the insect Diabrotica virgifera, a beetle in the family Chrysomelidae. Compared with noncoding, nontranscribed regions, EST-SSRs were generally less polymorphic but had reduced occurrence of null alleles and greater cross-species amplification. Neutrality tests suggested the loci were not under positive selection. Across all populations and all loci, the genomic and EST-SSRs performed similarly in estimating genetic diversity, F(IS), F(ST), population assignment and exclusion tests, and detection of distinct populations. These findings, therefore, indicate that the EST-SSRs examined can be used with confidence in future genetic studies of Diabrotica populations and suggest that EST libraries can be added as a valuable source of markers for population genetics studies in insects and other animals.     

A microsatellite map of white clover

The white clover (Trifolium repens) nuclear genome (n=2x=16) is an important yet under-characterised genetic environment. We have developed simple sequence repeat (SSR) genetic markers for the white clover genome by mining an expressed sequence tag (EST) database and by isolation from enriched genomic libraries. A total of 2,086 EST-derived SSRs (EST-SSRs) were identified among 26,480 database accessions. Evaluation of 792 EST-SSR primer pairs resulted in 566 usable EST-SSRs. Of these, 335 polymorphic EST-SSRs, used in concert with 30 genomic SSRs, detected 493 loci in the white clover genome using 92 F₁ progeny from a pair cross between two highly heterozygous genotypes—364/7 and 6525/5. Map length, as estimated using the joinmap algorithm, was 1,144 cM and spanned all 16 homologues. The R (red leaf) locus was mapped to linkage group B1 and is tightly linked to the microsatellite locus prs318c. The eight homoeologous pairs of linkage groups within the white clover genome were identified using 96 homoeologous loci. Segregation distortion was detected in four areas (groups A1, D1, D2 and H2). Marker locus density varied among and within linkage groups. This is the first time EST-SSRs have been used to build a whole-genome functional map and to describe subgenome organisation in an allopolyploid species, and T. repens is the only Trifolieae species to date to be mapped exclusively with SSRs. This gene-based microsatellite map will enable the resolution of quantitative traits into Mendelian characters, the characterisation of syntenic relationships with other genomes and acceleration of white clover improvement programmes.     

Developing and characterising Ricinus communis SSR markers by data mining of whole-genome sequences

Ricinus communis is a versatile industrial oil crop that is cultivated worldwide. Genetic improvement and marker-assisted breeding of castor bean have been slowed owing to the lack of abundant and efficient molecular markers. As co-dominant markers, simple sequence repeats (SSRs) are useful for genetic evaluation and molecular breeding. The recently released whole-genome sequence of castor bean provides useful genomic resources for developing markers on a genome-wide scale. In the present study, the distribution and frequency of microsatellites in the castor bean genome were characterised and numerous SSR markers were developed using genomic data mining. In total, 18,647 SSR loci at a density of one SSR per 18.89 Kb in the castor bean genome sequence (representing approximately 352.27 Mb) were identified. Dinucleotide repeats were the most frequently observed microsatellites, although the AAT repeat motif was also prevalent. Using six cultivars as screening samples, 670 polymorphic SSR markers from 1,435 primer pairs (46.7 %) were developed. Trinucleotide motif loci contained a higher proportion of polymorphisms (48.5 %) than dinucleotide motif loci (39.2 %). The polymorphism level in the SSR loci was positively correlated with the increasing number of repeat units in the microsatellites. The phylogenetic relationship among 32 varieties was evaluated using the developed SSR markers. Cultivars developed at the same institute clustered together, suggesting that these cultivars have a narrow genetic background. The large number of SSR markers developed in this study will be useful for genetic mapping and for breeding improved castor-oil plants. These markers will also facilitate genetic and genomic studies of Euphorbiaceae.     

Assessment of genetic diversity and population structure of Chinese wild almond, Amygdalus nana, using EST- and genomic SSRs

Amygdalus nana L., commonly known as wild almond, is an endangered wild relative of cultivated almond, which has great potential in almond crop breeding. In this study, we used microsatellite (SSR) loci derived from both expressed sequence tag (EST) and anonymous genomic sequence to explore the genetic diversity and population structure of A. nana in Xinjiang of China. Seven natural populations were collected across the whole distribution of A. nana in China, including populations from both inside (four populations) and outside (three populations) the established protected areas. A total of 22 and 19 alleles were detected from the seven pairs of EST and genomic SSR loci, respectively. Generally, the genomic SSRs showed lower levels of variation than EST-SSRs, which may partially due to the higher cross-species transferability in EST-SSRs than in genomic SSRs. The population-level genetic diversity (A = 1.84, P = 50.00%, H_o = 0.3491, H_E = 0.2271) was lower than cultivated almond and several wild fruit species with similar breeding system. Most of the genetic variation (82.16%) was partitioned within populations. In particular, the population collected from Tacheng County (outside the protected areas) had the highest levels of genetic diversity and had significantly different genetic constitution from other populations.     

EST-derived genic molecular markers: development and utilization for generating an advanced transcript map of chickpea

Well-saturated linkage maps especially those based on expressed sequence tag (EST)-derived genic molecular markers (GMMs) are a pre-requisite for molecular breeding. This is especially true in important legumes such as chickpea where few simple sequence repeats (SSR) and even fewer GMM-based maps have been developed. Therefore, in this study, 2,496 ESTs were generated from chickpea seeds and utilized for the development of 487 novel EST-derived functional markers which included 125 EST-SSRs, 151 intron targeted primers (ITPs), 109 expressed sequence tag polymorphisms (ESTPs), and 102 single nucleotide polymorphisms (SNPs). Whereas EST-SSRs, ITPs, and ESTPs were developed by in silico analysis of the developed EST sequences, SNPs were identified by allele resequencing and their genotyping was performed using the Illumina GoldenGate Assay. Parental polymorphism was analyzed between C. arietinum ICC4958 and C. reticulatum PI489777, parents of the reference chickpea mapping population, using a total of 872 markers: 487 new gene-based markers developed in this study along with 385 previously published markers, of which 318 (36.5%) were found to be polymorphic and were used for genotyping. The genotypic data were integrated with the previously published data of 108 markers and an advanced linkage map was generated that contained 406 loci distributed on eight linkage groups that spanned 1,497.7 cM. The average marker density was 3.68 cM and the average number of markers per LG was 50.8. Among the mapped markers, 303 new genomic locations were defined that included 177 gene-based and 126 gSSRs (genomic SSRs) thereby producing the most advanced gene-rich map of chickpea solely based on co-dominant markers.     

First Microsatellite Markers Developed from Cupuassu ESTs: Application in Diversity Analysis and Cross-Species Transferability to Cacao

The cupuassu tree (Theobroma grandiflorum) (Willd. ex Spreng.) Schum. is a fruitful species from the Amazon with great economical potential, due to the multiple uses of its fruit´s pulp and seeds in the food and cosmetic industries, including the production of cupulate, an alternative to chocolate. In order to support the cupuassu breeding program and to select plants presenting both pulp/seed quality and fungal disease resistance, SSRs from Next Generation Sequencing ESTs were obtained and used in diversity analysis. From 8,330 ESTs, 1,517 contained one or more SSRs (1,899 SSRs identified). The most abundant motifs identified in the EST-SSRs were hepta- and trinucleotides, and they were found with a minimum and maximum of 2 and 19 repeats, respectively. From the 1,517 ESTs containing SSRs, 70 ESTs were selected based on their functional annotation, focusing on pulp and seed quality, as well as resistance to pathogens. The 70 ESTs selected contained 77 SSRs, and among which, 11 were polymorphic in cupuassu genotypes. These EST-SSRs were able to discriminate the cupuassu genotype in relation to resistance/susceptibility to witches’ broom disease, as well as to pulp quality (SST/ATT values). Finally, we showed that these markers were transferable to cacao genotypes, and that genome availability might be used as a predictive tool for polymorphism detection and primer design useful for both Theobroma species. To our knowledge, this is the first report involving EST-SSRs from cupuassu and is also a pioneer in the analysis of marker transferability from cupuassu to cacao. Moreover, these markers might contribute to develop or saturate the cupuassu and cacao genetic maps, respectively.