Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: an integration of the BMc series

Mapping of adult plant recessive resistance to anthracnose in Indian common bean landrace Baspa/KRC 8

BACKGROUND

The common bean (Phaseolus vulgaris) has become the food of choice owing to its wealthy nutritional profile, leading to a considerable increase in its cultivation worldwide. However, anthracnose has been a major impediment to production and productivity, as elite bean cultivars are vulnerable to this disease. To overcome barriers in crop production, scientists worldwide are working towards enhancing the genetic diversity of crops. One way to achieve this is by introducing novel genes from related crops, including landraces like KRC 8. This particular landrace, found in the North Western Himalayan region, has shown adult plant resistance against anthracnose and also possesses a recessive resistance gene.

METHODS AND RESULTS

In this study, a population of 179 F2:9 RIL individuals (Jawala × KRC 8) was evaluated at both phenotypic and genotypic levels using over 830 diverse molecular markers to map the resistance gene present in KRC 8. We have successfully mapped a resistance gene to chromosome Pv01 using four SSR markers, namely IAC 238, IAC 235, IAC 259, and BM 146. The marker IAC 238 is closely linked to the gene with a distance of 0.29 cM, while the other markers flank the recessive resistance gene at 10.87 cM (IAC 259), 17.80 cM (BM 146), and 25.22 cM (IAC 235). Previously, a single recessive anthracnose resistance gene (co-8) has been reported in the common bean accession AB 136. However, when we performed PCR amplification with our tightly linked marker IAC 238, we got different amplicons in AB 136 and KRC 8. Interestingly, the susceptible cultivar Jawala produced the same amplicon as AB 136. This observation indicated that the recessive gene present in KRC 8 is different from co-8. As the gene is located far away from the Co-1 locus, we suggest naming the recessive gene co-Indb/co-19. Fine mapping of co-Indb in KRC 8 may provide new insights into the cloning and characterization of this recessive gene so that it can be incorporated into future bean improvement programs. Further, the tightly linked marker IAC 238 can be utilized in marker assisted introgression in future bean breeding programs.

CONCLUSION

The novel co-Indb gene present in Himalayan landrace KRC 8, showing adult plant resistance against common bean anthracnose, is independent from all the resistance genes previously located on chromosome Pv01.

Genome-Wide Association Study and Genomic Prediction of Fusarium Wilt Resistance in Common Bean Core Collection

The common bean (Phaseolus vulgaris L.) is a globally cultivated leguminous crop. Fusarium wilt (FW), caused by Fusarium oxysporum f. sp. phaseoli (Fop), is a significant disease leading to substantial yield loss in common beans. Disease-resistant cultivars are recommended to counteract this. The objective of this investigation was to identify single nucleotide polymorphism (SNP) markers associated with FW resistance and to pinpoint potential resistant common bean accessions within a core collection, utilizing a panel of 157 accessions through the Genome-wide association study (GWAS) approach with TASSEL 5 and GAPIT 3. Phenotypes for Fop race 1 and race 4 were matched with genotypic data from 4740 SNPs of BARCBean6K_3 Infinium Bea Chips. After ranking the 157-accession panel and revealing 21 Fusarium wilt-resistant accessions, the GWAS pinpointed 16 SNPs on chromosomes Pv04, Pv05, Pv07, Pv8, and Pv09 linked to Fop race 1 resistance, 23 SNPs on chromosomes Pv03, Pv04, Pv05, Pv07, Pv09, Pv10, and Pv11 associated with Fop race 4 resistance, and 7 SNPs on chromosomes Pv04 and Pv09 correlated with both Fop race 1 and race 4 resistances. Furthermore, within a 30 kb flanking region of these associated SNPs, a total of 17 candidate genes were identified. Some of these genes were annotated as classical disease resistance protein/enzymes, including NB-ARC domain proteins, Leucine-rich repeat protein kinase family proteins, zinc finger family proteins, P-loopcontaining nucleoside triphosphate hydrolase superfamily, etc. Genomic prediction (GP) accuracy for Fop race resistances ranged from 0.26 to 0.55. This study advanced common bean genetic enhancement through marker-assisted selection (MAS) and genomic selection (GS) strategies, paving the way for improved Fop resistance.

Genomics Associated Interventions for Heat Stress Tolerance in Cool Season Adapted Grain Legumes

Cool season grain legumes occupy an important place among the agricultural crops and essentially provide multiple benefits including food supply, nutrition security, soil fertility improvement and revenue for farmers all over the world. However, owing to climate change, the average temperature is steadily rising, which negatively affects crop performance and limits their yield. Terminal heat stress that mainly occurred during grain development phases severely harms grain quality and weight in legumes adapted to the cool season, such as lentils, faba beans, chickpeas, field peas, etc. Although, traditional breeding approaches with advanced screening procedures have been employed to identify heat tolerant legume cultivars. Unfortunately, traditional breeding pipelines alone are no longer enough to meet global demands. Genomics-assisted interventions including new-generation sequencing technologies and genotyping platforms have facilitated the development of high-resolution molecular maps, QTL/gene discovery and marker-assisted introgression, thereby improving the efficiency in legumes breeding to develop stress-resilient varieties. Based on the current scenario, we attempted to review the intervention of genomics to decipher different components of tolerance to heat stress and future possibilities of using newly developed genomics-based interventions in cool season adapted grain legumes.

Screening of mungbean for drought tolerance and transcriptome profiling between drought-tolerant and susceptible genotype in response to drought stress

Mungbean, is a widely cultivated pulse crop in India, experiences severe drought stress during the cultivation period. The mechanism of drought tolerance in mungbean is not well understood. In this presents study we screened 7 widely cultivated mungbean genotypes towards their drought sensitivity at seedling stage and transcriptome sequencing of drought-tolerant and susceptible genotype to understand the drought tolerance mechanism. Our physiological data such as increase in root length, shoot length, fresh weight, dry weight, relative water content (RWC), proline content, MDA content and molecular data in terms of quantitative expression of drought stress responsive genes under 3-d drought stress in mungbean suggests that, K851 seems to be most drought tolerant and PDM-139 as drought susceptible genotype. The transcriptomic study between K-851 and PDM-139 revealed 22,882 differentially expressed genes (DEGs) which were identified under drought stress, and they were mainly mapped to phytohormone signal transduction, carbon metabolism and flavonoid biosynthesis. Out of these, 10,235 genes were up-regulated and 12,647 genes were down-regulated. Furthermore, we found that, the DEGs related to, phytohormone signal transduction, carbon metabolism and flavonoid biosynthesis and they were more induced in K-851. Our data suggested that, the drought tolerant genotype K-851, scavenges the damage of drought stress by producing more amount of osmolytes, ROS scavengers and sugar biosynthesis.

Transcriptome sequencing, de novo assembly, characterisation of wild accession of blackgram (Vigna mungo var. silvestris) as a rich resource for development of molecular markers and validation of SNPs by high resolution melting (HRM) analysis

BACKGROUND

Blackgram [Vigna mungo (L.) Hepper], is an important legume crop of Asia with limited genomic resources. We report a comprehensive set of genic simple sequence repeat (SSR) and single nucleotide polymorphism (SNPs) markers using Illumina MiSeq sequencing of transcriptome and its application in genetic variation analysis and mapping.

RESULTS

Transcriptome sequencing of immature seeds of wild blackgram, V. mungo var. silvestris by Illumina MiSeq technology generated 1.9 × 10⁷ reads, which were assembled into 40,178 transcripts (TCS) with an average length of 446 bp covering 2.97 GB of the genome. A total of 38,753 CDS (Coding sequences) were predicted from 40,178 TCS and 28,984 CDS were annotated through BLASTX and mapped to GO and KEGG database resulting in 140 unique pathways. The tri-nucleotides were most abundant (39.9%) followed by di-nucleotide (30.2%). About 60.3 and 37.6% of SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Among SNPs, the most abundant substitution type were transitions (Ts) (61%) followed by transversions (Tv) type (39%), with a Ts/Tv ratio of 1.58. A total of 2306 DEGs were identified by RNA Seq between wild and cultivar and validation was done by quantitative reverse transcription polymerase chain reaction. In this study, we genotyped SNPs with a validation rate of 78.87% by High Resolution Melting (HRM) Assay.

CONCLUSION

In the present study, 1621genic-SSR and 1844 SNP markers were developed from immature seed transcriptome sequence of blackgram and 31 genic-SSR markers were used to study genetic variations among different blackgram accessions. Above developed markers contribute towards enriching available genomic resources for blackgram and aid in breeding programmes.

Milestones achieved in response to drought stress through reverse genetic approaches

Drought stress is the most important abiotic stress that constrains crop production and reduces yield drastically. The germplasm of most of the cultivated crops possesses numerous unknown drought stress tolerant genes. Moreover, there are many reports suggesting that the wild species of most of the modern cultivars have abiotic stress tolerant genes. Due to climate change and population booms, food security has become a global issue. To develop drought tolerant crop varieties knowledge of various genes involved in drought stress is required. Different reverse genetic approaches such as virus-induced gene silencing (VIGS), clustered regularly interspace short palindromic repeat (CRISPR), targeting induced local lesions in genomes (TILLING) and expressed sequence tags (ESTs) have been used extensively to study the functionality of different genes involved in response to drought stress. In this review, we described the contributions of different techniques of functional genomics in the study of drought tolerant genes.

Bta-miR-130a regulates the biosynthesis of bovine milk fat by targeting peroxisome proliferator-activated receptor gamma

Milk fat determines the quality of milk and is also a main targeted trait in dairy cow breeding. Recent studies have revealed important regulatory roles of microRNAs (miRNA) in milk fat synthesis in the mammary gland. However, the role of miRNA in bovine mammary epithelial cells (BMEC) remains largely unknown. In this study, we found that the overexpression of miR-130a significantly decreased cellular triacylglycerol (TAG) levels and suppressed lipid droplet formation, whereas the inhibition of miR-130a resulted in greater lipid droplet formation and TAG accumulation in BMEC. MiR-130a also significantly affected mRNA expression related to milk fat metabolism. Specifically, the overexpression of miR-130a reduced the mRNA expression of , , , and , whereas the downregulation of miR-130a increased the mRNA expression of , , , , , and . Furthermore, western blot analysis revealed the protein level of PPARG in miR-130a mimic and inhibitor transfection groups to be consistent with the mRNA expression response. Finally, luciferase reporter assays verified that PPARG was the direct target of miR-130a. This study provides the first experimental evidence that miR-130a directly affects TAG synthesis in BMEC by targeting PPARG, suggesting that miR-130a potentially could be used to improve beneficial milk components in dairy cows.

Dissection of Resistance Genes to Pseudomonas syringae pv. phaseolicola in UI3 Common Bean Cultivar

Few quantitative trait loci have been mapped for resistance to Pseudomonas syringae pv. phaseolicola in common bean. Two F₂ populations were developed from the host differential UI3 cultivar. The objective of this study was to further characterize the resistance to races 1, 5, 7 and 9 of Psp included in UI3. Using a QTL mapping approach, 16 and 11 main-effect QTLs for pod and primary leaf resistance were located on LG10, explaining up to 90% and 26% of the phenotypic variation, respectively. The homologous genomic region corresponding to primary leaf resistance QTLs detected tested positive for the presence of resistance-associated gene cluster encoding nucleotide-binding and leucine-rich repeat (NL), Natural Resistance Associated Macrophage (NRAMP) and Pentatricopeptide Repeat family (PPR) proteins. It is worth noting that the main effect QTLs for resistance in pod were located inside a 3.5 Mb genomic region that included the Phvul.010G021200 gene, which encodes a protein that has the highest sequence similarity to the RIN4 gene of Arabidopsis, and can be considered an important candidate gene for the organ-specific QTLs identified here. These results support that resistance to Psp from UI3 might result from the immune response activated by combinations of R proteins, and suggest the guard model as an important mechanism in pod resistance to halo blight. The candidate genes identified here warrant functional studies that will help in characterizing the actual defense gene(s) in UI3 genotype.

Mapping QTLs for drought tolerance in a SEA 5 x AND 277 common bean cross with SSRs and SNP markers

The common bean is characterized by high sensitivity to drought and low productivity. Breeding for drought resistance in this species involves genes of different genetic groups. In this work, we used a SEA 5 x AND 277 cross to map quantitative trait loci associated with drought tolerance in order to assess the factors that determine the magnitude of drought response in common beans. A total of 438 polymorphic markers were used to genotype the F8 mapping population. Phenotyping was done in two greenhouses, one used to simulate drought and the other to simulate irrigated conditions. Fourteen traits associated with drought tolerance were measured to identify the quantitative trait loci (QTLs). The map was constructed with 331 markers that covered all 11 chromosomes and had a total length of 1515 cM. Twenty-two QTLs were discovered for chlorophyll, leaf and stem fresh biomass, leaf biomass dry weight, leaf temperature, number of pods per plant, number of seeds per plant, seed weight, days to flowering, dry pod weight and total yield under well-watered and drought (stress) conditions. All the QTLs detected under drought conditions showed positive effects of the SEA 5 allele. This study provides a better understanding of the genetic inheritance of drought tolerance in common bean.

Comparative Transcriptome Analysis Identified Candidate Genes Related to Bailinggu Mushroom Formation and Genetic Markers for Genetic Analyses and Breeding

Bailinggu (Pleurotus tuoliensis) is a major, commercially cultivated mushroom and widely used for nutritional, medicinal, and industrial applications. Yet, the mushroom’s genetic architecture and the molecular mechanisms underlying its formation are largely unknown. Here we performed comparative transcriptomic analysis during Bailinggu’s mycelia, primordia, and fruiting body stages to identify genes regulating fruiting body development and develop EST-SSR markers assessing the genetic value of breeding materials. The stage-specific and differentially expressed unigenes (DEGs) involved in morphogenesis, primary carbohydrate metabolism, cold stimulation and blue-light response were identified using GO and KEGG databases. These unigenes might help Bailinggu adapt to genetic and environmental factors that influence fructification. The most pronounced change in gene expression occurred during the vegetative-to-reproductive transition, suggesting that is most active and key for Bailinggu development. We then developed 26 polymorphic and informative EST-SSR markers to assess the genetic diversity in 82 strains of Bailinggu breeding materials. These EST-SSRs exhibited high transferability in closely related species P. eryngii var. ferulae and var. eryngii. Genetic population structure analysis indicated that China’s Bailinggu has low introgression with these two varieties and likely evolved independently. These findings provide new genes, SSR markers, and germplasm to enhance the breeding of commercially cultivated Bailinggu.

De novo transcriptomic analysis of cowpea (Vigna unguiculata L. Walp.) for genic SSR marker development

Background

Cowpea [Vigna unguiculata (L.) Walp.] is one of the most important legumes in tropical and semi-arid regions. However, there is relatively little genomic information available for genetic research on and breeding of cowpea. The objectives of this study were to analyse the cowpea transcriptome and develop genic molecular markers for future genetic studies of this genus.

Results

Approximately 54 million high-quality cDNA sequence reads were obtained from cowpea based on Illumina paired-end sequencing technology and were de novo assembled to generate 47,899 unigenes with an N50 length of 1534 bp. Sequence similarity analysis revealed 36,289 unigenes (75.8%) with significant similarity to known proteins in the non-redundant (Nr) protein database, 23,471 unigenes (49.0%) with BLAST hits in the Swiss-Prot database, and 20,654 unigenes (43.1%) with high similarity in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Further analysis identified 5560 simple sequence repeats (SSRs) as potential genic molecular markers. Validating a random set of 500 SSR markers yielded 54 polymorphic markers among 32 cowpea accessions.

Conclusions

This transcriptomic analysis of cowpea provided a valuable set of genomic data for characterizing genes with important agronomic traits in Vigna unguiculata and a new set of genic SSR markers for further genetic studies and breeding in cowpea and related Vigna species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0531-5) contains supplementary material, which is available to authorized users.

Development of Gene-Based SSR Markers in Winged Bean (Psophocarpus tetragonolobus (L.) DC.) for Diversity Assessment

Winged bean (Psophocarpus tetragonolobus) is an herbaceous multipurpose legume grown in hot and humid countries as a pulse, vegetable (leaves and pods), or root tuber crop depending on local consumption preferences. In addition to its different nutrient-rich edible parts which could contribute to food and nutritional security, it is an efficient nitrogen fixer as a component of sustainable agricultural systems. Generating genetic resources and improved lines would help to accelerate the breeding improvement of this crop, as the lack of improved cultivars adapted to specific environments has been one of the limitations preventing wider use. A transcriptomic de novo assembly was constructed from four tissues: leaf, root, pod, and reproductive tissues from Malaysian accessions, comprising of 198,554 contigs with a N50 of 1462 bp. Of these, 138,958 (70.0%) could be annotated. Among 9682 genic simple sequence repeat (SSR) motifs identified (excluding monomer repeats), trinucleotide-repeats were the most abundant (4855), followed by di-nucleotide (4500) repeats. A total of 18 SSR markers targeting di- and tri-nucleotide repeats have been validated as polymorphic markers based on an initial assessment of nine genotypes originated from five countries. A cluster analysis revealed provisional clusters among this limited, yet diverse selection of germplasm. The developed assembly and validated genic SSRs in this study provide a foundation for a better understanding of the plant breeding system for the genetic improvement of winged bean.

miR-27a controls triacylglycerol synthesis in bovine mammary epithelial cells by targeting peroxisome proliferator-activated receptor gamma

Growing evidence has revealed that microRNA are central elements in milk fat synthesis in mammary epithelial cells. A negative regulator of adipocyte fat synthesis, miR-27a has been reported to be involved in the regulation of milk fat synthesis in goat mammary epithelial cells; however, the regulatory role of miR-27a in bovine milk fat synthesis remains unclear. In the present study, primary bovine mammary epithelial cells (BMEC) were harvested from mid-lactation cows and cultured in Dulbecco's modified Eagle's medium/F-12 medium with 10% fetal bovine serum, 5 μg/mL of insulin, 1 μg/mL of hydrocortisone, 2 μg/mL of prolactin, 1 μg/mL of progesterone, 100 U/mL of penicillin, and 100 μg/mL of streptomycin. We found that the overexpression of miR-27a significantly suppressed lipid droplet formation and decreased the cellular triacylglycerol (TAG) levels, whereas inhibition of miR-27a resulted in a greater lipid droplet formation and TAG accumulation in BMEC. Meanwhile, overexpression of miR-27a inhibited mRNA expression of peroxisome proliferator-activated receptor gamma (PPARG), CCAAT/enhancer-binding protein beta (C/EBPβ), perilipin 2 (PLIN2), and fatty acid binding protein 3 (FABP3), whereas miR-27a downregulation increased PPARG, C/EBPβ, FABP3, and CCAAT enhancer binding protein alpha (C/EBPα) mRNA expression. Furthermore, Western blot analysis revealed the protein level of PPARG in miR-27a mimic and inhibitor transfection groups to be consistent with the mRNA expression response. Moreover, luciferase reporter assays verified that PPARG was the direct target of miR-27a. In summary, these results indicate that miR-27a has the ability to control TAG synthesis in BMEC via targeting PPARG, suggesting that miR-27a could potentially be used to improve beneficial milk components in dairy cows.

Linkage and mapping of quantitative trait loci associated with angular leaf spot and powdery mildew resistance in common beans

Angular leaf spot (ALS) and powdery mildew (PWM) are two important fungi diseases causing significant yield losses in common beans. In this study, a new genetic linkage map was constructed using single sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs), in a segregating population derived from the AND 277 x SEA 5 cross, with 105 recombinant inbred lines. Phenotypic evaluations were performed in the greenhouse to identify quantitative trait loci (QTLs) associated with resistance by means of the composite interval mapping analysis. Four QTLs were identified for ALS resistance. The QTL ALS11^AS, linked on the SNP BAR 5054, mapped on chromosome Pv11, showed the greatest effect (R² = 26.5%) on ALS phenotypic variance. For PWM resistance, two QTLs were detected, PWM2^AS and PWM11^AS, on Pv2 and Pv11, explaining 7% and 66% of the phenotypic variation, respectively. Both QTLs on Pv11 were mapped on the same genomic region, suggesting that it is a pleiotropic region. The present study resulted in the identification of new markers closely linked to ALS and PWM QTLs, which can be used for marker-assisted selection, fine mapping and positional cloning.

Development and Application of Transcriptome-Derived Microsatellites in Actinidia eriantha (Actinidiaceae)

Actinidia eriantha Benth. is a diploid perennial woody vine native to China and is recognized as a valuable species for commercial kiwifruit improvement with high levels of ascorbic acid as well as having been used in traditional Chinese medicine. Due to the lack of genomic resources for the species, microsatellite markers for population genetics studies are scarce. In this study, RNASeq was conducted on fruit tissue of A. eriantha, yielding 5,678,129 reads with a total output of 3.41 Gb. De novo assembly yielded 69,783 non-redundant unigenes (41.3 Mb), of which 21,730 were annotated using protein databases. A total of 8,658 EST-SSR loci were identified in 7,495 unigene sequences, for which primer pairs were successfully designed for 3,842 loci (44.4%). Among these, 183 primer pairs were assayed for PCR amplification, yielding 69 with detectable polymorphism in A. eriantha. Additionally, 61 of the 69 polymorphic loci could be successfully amplified in at least one other Actinidia species. Of these, 14 polymorphic loci (mean N_A = 6.07 ± 2.30) were randomly selected for assessing levels of genetic diversity and population structure within A. eriantha. Finally, a neighbor-joining tree and Bayesian clustering analysis showed distinct clustering into two groups (K = 2), agreeing with the geographical distributions of these populations. Overall, our results will facilitate further studies of genetic diversity within A. eriantha and will aid in discriminating outlier loci involved in local adaptation.

Dynamics of a Novel Highly Repetitive CACTA Family in Common Bean (Phaseolus vulgaris)

Transposons are ubiquitous genomic components that play pivotal roles in plant gene and genome evolution. We analyzed two genome sequences of common bean (Phaseolus vulgaris) and identified a new CACTA transposon family named pvCACTA1. The family is extremely abundant, as more than 12,000 pvCACTA1 elements were found. To our knowledge, this is the most abundant CACTA family reported thus far. The computational and fluorescence in situ hybridization (FISH) analyses indicated that the pvCACTA1 elements were concentrated in terminal regions of chromosomes and frequently generated AT-rich 3 bp target site duplications (TSD, WWW, W is A or T). Comparative analysis of the common bean genomes from two domesticated genetic pools revealed that new insertions or excisions of pvCACTA1 elements occurred after the divergence of the two common beans, and some of the polymorphic elements likely resulted in variation in gene sequences. pvCACTA1 elements were detected in related species but not outside the Phaseolus genus. We calculated the molecular evolutionary rate of pvCACTA1 transposons using orthologous elements that indicated that most transposition events likely occurred before the divergence of the two gene pools. These results reveal unique features and evolution of this new transposon family in the common bean genome.

Development of Gene-Based SSR Markers in Rice Bean (Vigna umbellata L.) Based on Transcriptome Data

Rice bean (Vigna umbellata (Thunb.) Ohwi & Ohashi) is a warm season annual legume mainly grown in East Asia. Only scarce genomic resources are currently available for this legume crop species and no simple sequence repeat (SSR) markers have been specifically developed for rice bean yet. In this study, approximately 26 million high quality cDNA sequence reads were obtained from rice bean using Illumina paired-end sequencing technology and assembled into 71,929 unigenes with an average length of 986 bp. Of these unigenes, 38,840 (33.2%) showed significant similarity to proteins in the NCBI non-redundant protein and nucleotide sequence databases. Furthermore, 30,170 (76.3%) could be classified into gene ontology categories, 25,451 (64.4%) into Swiss-Prot categories and 21,982 (55.6%) into KOG database categories (E-value < 1.0E-5). A total of 9,301 (23.5%) were mapped onto 118 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database. A total of 3,011 genic SSRs were identified as potential molecular markers. AG/CT (30.3%), AAG/CTT (8.1%) and AGAA/TTCT (20.0%) are the three main repeat motifs. A total of 300 SSR loci were randomly selected for validation by using PCR amplification. Of these loci, 23 primer pairs were polymorphic among 32 rice bean accessions. A UPGMA dendrogram revealed three major clusters among 32 rice bean accessions. The large number of SSR-containing sequences and genic SSRs in this study will be valuable for the construction of high-resolution genetic linkage maps, association or comparative mapping and genetic analyses of various Vigna species.

Population Genetics of the Rubber-Producing Russian Dandelion (Taraxacum kok-saghyz)

The Russian dandelion, Taraxacum kok-saghyz (TKS), is a perennial species native to Central Asia that produces high quality, natural rubber. Despite its potential to help maintain a stable worldwide rubber supply, little is known about genetic variation in this species. To facilitate future germplasm improvement efforts, we developed simple-sequence repeat (SSR) markers from available expressed-sequence tag (EST) data and used them to investigate patterns of population genetic diversity in this nascent crop species. We identified numerous SSRs (1,510 total) in 1,248 unigenes from a larger set of 6,960 unigenes (derived from 16,441 ESTs) and designed PCR primers targeting 767 of these loci. Screening of a subset of 192 of these primer pairs resulted in the identification of 48 pairs that appeared to produce single-locus polymorphisms. We then used the most reliable 17 of these primer pairs to genotype 176 individuals from 17 natural TKS populations. We observed an average of 4.8 alleles per locus with population-level expected heterozygosities ranging from 0.28 to 0.50. An average pairwise F_ST of 0.11 indicated moderate but statistically significant levels of genetic differentiation, though there was no clear geographic patterning to this differentiation. We also tested these 17 primer pairs in the widespread common dandelion, T. officinale, and a majority successfully produced apparently single-locus amplicons. This result demonstrates the potential utility of these markers for genetic analyses in other species in the genus.

Major Contribution of Flowering Time and Vegetative Growth to Plant Production in Common Bean As Deduced from a Comparative Genetic Mapping

Determinacy growth habit and accelerated flowering traits were selected during or after domestication in common bean. Both processes affect several presumed adaptive traits such as the rate of plant production. There is a close association between flowering initiation and vegetative growth; however, interactions among these two crucial developmental processes and their genetic bases remain unexplored. In this study, with the aim to establish the genetic relationships between these complex processes, a multi-environment quantitative trait locus (QTL) mapping approach was performed in two recombinant inbred line populations derived from inter-gene pool crosses between determinate and indeterminate genotypes. Additive and epistatic QTLs were found to regulate flowering time, vegetative growth, and rate of plant production. Moreover, the pleiotropic patterns of the identified QTLs evidenced that regions controlling time to flowering traits, directly or indirectly, are also involved in the regulation of plant production traits. Further QTL analysis highlighted one QTL, on the lower arm of the linkage group Pv01, harboring the Phvul.001G189200 gene, homologous to the Arabidopsis thaliana TERMINAL FLOWER1 (TFL1) gene, which explained up to 32% of phenotypic variation for time to flowering, 66% for vegetative growth, and 19% for rate of plant production. This finding was consistent with previous results, which have also suggested Phvul.001G189200 (PvTFL1y) as a candidate gene for determinacy locus. The information here reported can also be applied in breeding programs seeking to optimize key agronomic traits, such as time to flowering, plant height and an improved reproductive biomass, pods, and seed size, as well as yield.

SNP Assay Development for Linkage Map Construction, Anchoring Whole-Genome Sequence, and Other Genetic and Genomic Applications in Common Bean

A total of 992,682 single-nucleotide polymorphisms (SNPs) was identified as ideal for Illumina Infinium II BeadChip design after sequencing a diverse set of 17 common bean (Phaseolus vulgaris L) varieties with the aid of next-generation sequencing technology. From these, two BeadChips each with >5000 SNPs were designed. The BARCBean6K_1 BeadChip was selected for the purpose of optimizing polymorphism among market classes and, when possible, SNPs were targeted to sequence scaffolds in the Phaseolus vulgaris 14× genome assembly with sequence lengths >10 kb. The BARCBean6K_2 BeadChip was designed with the objective of anchoring additional scaffolds and to facilitate orientation of large scaffolds. Analysis of 267 F2 plants from a cross of varieties Stampede × Red Hawk with the two BeadChips resulted in linkage maps with a total of 7040 markers including 7015 SNPs. With the linkage map, a total of 432.3 Mb of sequence from 2766 scaffolds was anchored to create the Phaseolus vulgaris v1.0 assembly, which accounted for approximately 89% of the 487 Mb of available sequence scaffolds of the Phaseolus vulgaris v0.9 assembly. A core set of 6000 SNPs (BARCBean6K_3 BeadChip) with high genotyping quality and polymorphism was selected based on the genotyping of 365 dry bean and 134 snap bean accessions with the BARCBean6K_1 and BARCBean6K_2 BeadChips. The BARCBean6K_3 BeadChip is a useful tool for genetics and genomics research and it is widely used by breeders and geneticists in the United States and abroad.

Development and Validation of EST-SSR Markers from the Transcriptome of Adzuki Bean (Vigna angularis)

The adzuki bean (Vigna angularis (Ohwi) Ohwi and Ohashi) is an important grain legume of Asia. It is cultivated mainly in China, Japan and Korea. Despite its importance, few genomic resources are available for molecular genetic research of adzuki bean. In this study, we developed EST-SSR markers for the adzuki bean through next-generation sequencing. More than 112 million high-quality cDNA sequence reads were obtained from adzuki bean using Illumina paired-end sequencing technology, and the sequences were de novo assembled into 65,950 unigenes. The average length of the unigenes was 1,213 bp. Among the unigenes, 14,547 sequences contained a unique simple sequence repeat (SSR) and 3,350 sequences contained more than one SSR. A total of 7,947 EST-SSRs were identified as potential molecular markers, with mono-nucleotide A/T repeats (99.0%) as the most abundant motif class, followed by AG/CT (68.4%), AAG/CTT (30.0%), AAAG/CTTT (26.2%), AAAAG/CTTTT (16.1%), and AACGGG/CCCGTT (6.0%). A total of 500 SSR markers were randomly selected for validation, of which 296 markers produced reproducible amplicons with 38 polymorphic markers among the 32 adzuki bean genotypes selected from diverse geographical locations across China. The large number of SSR-containing sequences and EST-SSR markers will be valuable for genetic analysis of the adzuki bean and related Vigna species.

De novo Assembly, Characterization of Immature Seed Transcriptome and Development of Genic-SSR Markers in Black Gram [Vigna mungo (L.) Hepper]

Black gram [V. mungo (L.) Hepper] is an important legume crop extensively grown in south and south-east Asia, where it is a major source of dietary protein for its predominantly vegetarian population. However, lack of genomic information and markers has become a limitation for genetic improvement of this crop. Here, we report the transcriptome sequencing of the immature seeds of black gram cv. TU94-2, by Illumina paired end sequencing technology to generate transcriptome sequences for gene discovery and genic-SSR marker development. A total of 17.2 million paired-end reads were generated and 48,291 transcript contigs (TCS) were assembled with an average length of 443 bp. Based on sequence similarity search, 33,766 TCS showed significant similarity to known proteins. Among these, only 29,564 TCS were annotated with gene ontology (GO) functional categories. A total number of 138 unique KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways were identified, of which majority of TCS are grouped into purine metabolism (678) followed by pyrimidine metabolism (263). A total of 48,291 TCS were searched for SSRs and 1,840 SSRs were identified in 1,572 TCS with an average frequency of one SSR per 11.9 kb. The tri-nucleotide repeats were most abundant (35%) followed by di-nucleotide repeats (32%). PCR primer pairs were successfully designed for 933 SSR loci. Sequences analyses indicate that about 64.4% and 35.6% of the SSR motifs were present in the coding sequences (CDS) and untranslated regions (UTRs) respectively. Tri-nucleotide repeats (57.3%) were preferentially present in the CDS. The rate of successful amplification and polymorphism were investigated using selected primers among 18 black gram accessions. Genic-SSR markers developed from the Illumina paired end sequencing of black gram immature seed transcriptome will provide a valuable resource for genetic diversity, evolution, linkage mapping, comparative genomics and marker-assisted selection in black gram.

Identification and analysis of novel salt responsive candidate gene based SSRs (cgSSRs) from rice (Oryza sativa L.)

BACKGROUND

Majority of the Asian people depend on rice for nutritional energy. Rice cultivation and yield are severely affected by soil salinity stress worldwide. Marker assisted breeding is a rapid and efficient way to develop improved variety for salinity stress tolerance. Genomic microsatellite markers are an elite group of markers, but there is possible uncertainty of linkage with the important genes. In contrast, there are better possibilities of linkage detection with important genes if SSRs are developed from candidate genes. To the best of our knowledge, there is no such report on SSR markers development from candidate gene sequences in rice. So the present study was aimed to identify and analyse SSRs from salt responsive candidate genes of rice.

RESULTS

In the present study, based on the comprehensive literature survey, we selected 220 different salt responsive genes of rice. Out of them, 106 genes were found to contain 180 microsatellite loci with, tri-nucleotide motifs (56%) being most abundant, followed by di-(41%) and tetra nucleotide (2.8%) motifs. Maximum loci were found in the coding sequences (37.2%), followed by in 5'UTR (26%), intron (21.6%) and 3'UTR (15%). For validation, 19 primer sets were evaluated to detect polymorphism in diversity analysis among the two panels consisting of 17 salt tolerant and 17 susceptible rice genotypes. Except one, all primer sets exhibited polymorphic nature with an average of 21.8 alleles/primer and with a mean PIC value of 0.28. Calculated genetic similarity among genotypes was ranged from 19%-89%. The generated dendrogram showed 3 clusters of which one contained entire 17 susceptible genotypes and another two clusters contained all tolerant genotypes.

CONCLUSION

The present study represents the potential of salt responsive candidate gene based SSR (cgSSR) markers to be utilized as novel and remarkable candidate for diversity analysis among rice genotypes differing in salinity response.

Assessment of genetic diversity and population structure of mung bean (Vigna radiata) germplasm using EST-based and genomic SSR markers

Mung bean is an important legume crop in tropical and subtropical countries of Asia and has high nutritional and economic value. However the genetic diversity of mung bean is poorly characterized. In this study, our goal was to develop and use microsatellite simple sequence repeat (SSR) markers for germplasm evaluation. In total, 500 novel expression sequence tag EST-based SSRs (eSSRs) and genomic SSRs (gSSRs) were developed from mung bean transcriptome and genome sequences. Of these, only 58 were useful for diversity evaluation in a panel of 157 cultivated and wild mung bean accessions from different collection sites in East Asia. A total of 2.66 alleles were detected on average per locus which shows that polymorphism is generally low for the species. The average polymorphic information content (PIC) of gSSRs was higher than eSSRs and most of the polymorphic gSSRs were composed of di- and tri-nucleotide repeats (52.4% and 38.1% of all loci, respectively). The genotypes were differentiated into nine subgroups by cluster analysis, and the wild mung bean accessions separated well from the cultivated accessions. Analysis of molecular variance indicated that 22% of variance was observed among populations and 78% was due to differences within populations. Clustering, population structure analyses showed that non-Chinese cultivated and wild mung bean accessions were separated from Chinese accessions, but no geographical distinctions existed between genotypes collected in China. Interestingly, the average PIC value of cultivated mung bean (0.36) was higher than that of wild mung bean (0.25) showing that further collecting and wide crosses are necessary for mung bean improvement.

Transcriptome sequencing of mung bean (Vigna radiate L.) genes and the identification of EST-SSR markers

Mung bean (Vigna radiate (L.) Wilczek) is an important traditional food legume crop, with high economic and nutritional value. It is widely grown in China and other Asian countries. Despite its importance, genomic information is currently unavailable for this crop plant species or some of its close relatives in the Vigna genus. In this study, more than 103 million high quality cDNA sequence reads were obtained from mung bean using Illumina paired-end sequencing technology. The processed reads were assembled into 48,693 unigenes with an average length of 874 bp. Of these unigenes, 25,820 (53.0%) and 23,235 (47.7%) showed significant similarity to proteins in the NCBI non-redundant protein and nucleotide sequence databases, respectively. Furthermore, 19,242 (39.5%) could be classified into gene ontology categories, 18,316 (37.6%) into Swiss-Prot categories and 10,918 (22.4%) into KOG database categories (E-value < 1.0E-5). A total of 6,585 (8.3%) were mapped onto 244 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database. Among the unigenes, 10,053 sequences contained a unique simple sequence repeat (SSR), and 2,303 sequences contained more than one SSR together in the same expressed sequence tag (EST). A total of 13,134 EST-SSRs were identified as potential molecular markers, with mono-nucleotide A/T repeats being the most abundant motif class and G/C repeats being rare. In this SSR analysis, we found five main repeat motifs: AG/CT (30.8%), GAA/TTC (12.6%), AAAT/ATTT (6.8%), AAAAT/ATTTT (6.2%) and AAAAAT/ATTTTT (1.9%). A total of 200 SSR loci were randomly selected for validation by PCR amplification as EST-SSR markers. Of these, 66 marker primer pairs produced reproducible amplicons that were polymorphic among 31 mung bean accessions selected from diverse geographical locations. The large number of SSR-containing sequences found in this study will be valuable for the construction of a high-resolution genetic linkage maps, association or comparative mapping and genetic analyses of various Vigna species.

Transcriptome sequencing of different narrow-leafed lupin tissue types provides a comprehensive uni-gene assembly and extensive gene-based molecular markers

Narrow-leafed lupin (NLL; Lupinus angustifolius L.) is an important grain legume crop that is valuable for sustainable farming and is becoming recognized as a human health food. NLL breeding is directed at improving grain production, disease resistance, drought tolerance and health benefits. However, genetic and genomic studies have been hindered by a lack of extensive genomic resources for the species. Here, the generation, de novo assembly and annotation of transcriptome datasets derived from five different NLL tissue types of the reference accession cv. Tanjil are described. The Tanjil transcriptome was compared to transcriptomes of an early domesticated cv. Unicrop, a wild accession P27255, as well as accession 83A:476, together being the founding parents of two recombinant inbred line (RIL) populations. In silico predictions for transcriptome-derived gene-based length and SNP polymorphic markers were conducted and corroborated using a survey assembly sequence for NLL cv. Tanjil. This yielded extensive indel and SNP polymorphic markers for the two RIL populations. A total of 335 transcriptome-derived markers and 66 BAC-end sequence-derived markers were evaluated, and 275 polymorphic markers were selected to genotype the reference NLL 83A:476 × P27255 RIL population. This significantly improved the completeness, marker density and quality of the reference NLL genetic map.

Analysis of BAC-end sequences in common bean (Phaseolus vulgaris L.) towards the development and characterization of long motifs SSRs

The increasing volume of genomic data on the Phaseolus vulgaris species have contributed to its importance as a model genetic species and positively affected the investigation of other legumes of scientific and economic value. To expand and gain a more in-depth knowledge of the common bean genome, the ends of a number of bacterial artificial chromosome (BAC) were sequenced, annotated and the presence of repetitive sequences was determined. In total, 52,270 BESs (BAC-end sequences), equivalent to 32 Mbp (~6 %) of the genome, were processed. In total, 3,789 BES-SSRs were identified, with a distribution of one SSR (simple sequence repeat) per 8.36 kbp and 2,000 were suitable for the development of SSRs, of which 194 were evaluated in low-resolution screening. From 40 BES-SSRs based on long motifs SSRs (≥ trinucleotides) analyzed in high-resolution genotyping, 34 showed an equally good amplification for the Andean and for the Mesoamerican genepools, exhibiting an average gene diversity (H E) of 0.490 and 5.59 alleles/locus, of which six classified as Class I showed a H E ≥ 0.7. The PCoA and structure analysis allowed to discriminate the gene pools (K = 2, FST = 0.733). From the 52,270 BESs, 2 % corresponded to transcription factors and 3 % to transposable elements. Putative functions for 24,321 BESs were identified and for 19,363 were assigned functional categories (gene ontology). This study identified highly polymorphic BES-SSRs containing tri- to hexanucleotides motifs and bringing together relevant genetic characteristics useful for breeding programs. Additionally, the BESs were incorporated into the international genome-sequencing project for the common bean.

De novo assembly of the common bean transcriptome using short reads for the discovery of drought-responsive genes

The common bean (Phaseolus vulgaris L.) is one of the most important food legumes, far ahead of other legumes. The average grain yield of the common bean worldwide is much lower than its potential yields, primarily due to drought in the field. However, the gene network that mediates plant responses to drought stress remains largely unknown in this species. The major goals of our study are to identify a large scale of genes involved in drought stress using RNA-seq. First, we assembled 270 million high-quality trimmed reads into a non-redundant set of 62,828 unigenes, representing approximately 49 Mb of unique transcriptome sequences. Of these unigenes, 26,501 (42.2%) common bean unigenes had significant similarity with unigenes/predicted proteins from other legumes or sequenced plants. All unigenes were functionally annotated within the GO, COG and KEGG pathways. The strategy for de novo assembly of transcriptome data generated here will be useful in other legume plant transcriptome studies. Second, we identified 10,482 SSRs and 4,099 SNPs in transcripts. The large number of genetic markers provides a resource for gene discovery and development of functional molecular markers. Finally, we found differential expression genes (DEGs) between terminal drought and optimal irrigation treatments and between the two different genotypes Long 22-0579 (drought tolerant) and Naihua (drought sensitive). DEGs were confirmed by quantitative real-time PCR assays, which indicated that these genes are functionally associated with the drought-stress response. These resources will be helpful for basic and applied research for genome analysis and crop drought resistance improvement in the common bean.

BAC-end microsatellites from intra and inter-genic regions of the common bean genome and their correlation with cytogenetic features

Highly polymorphic markers such as simple sequence repeats (SSRs) or microsatellites are very useful for genetic mapping. In this study novel SSRs were identified in BAC-end sequences (BES) from non-contigged, non-overlapping bacterial artificial clones (BACs) in common bean (Phaseolus vulgaris L.). These so called "singleton" BACs were from the G19833 Andean gene pool physical map and the new BES-SSR markers were used for the saturation of the inter-gene pool, DOR364×G19833 genetic map. A total of 899 SSR loci were found among the singleton BES, but only 346 loci corresponded to the single di- or tri-nucleotide motifs that were likely to be polymorphic (ATT or AG motifs, principally) and useful for primer design and individual marker mapping. When these novel SSR markers were evaluated in the DOR364×G19833 population parents, 136 markers revealed polymorphism and 106 were mapped. Genetic mapping resulted in a map length of 2291 cM with an average distance between markers of 5.2 cM. The new genetic map was compared to the most recent cytogenetic analysis of common bean chromosomes. We found that the new singleton BES-SSR were helpful in filling peri-centromeric spaces on the cytogenetic map. Short genetic distances between some new singleton-derived BES-SSR markers was common showing suppressed recombination in these regions compared to other parts of the genome. The correlation of singleton-derived SSR marker distribution with other cytogenetic features of the bean genome is discussed.

Annotation and sequence diversity of transposable elements in common bean (Phaseolus vulgaris)

Common bean (Phaseolus vulgaris) is an important legume crop grown and consumed worldwide. With the availability of the common bean genome sequence, the next challenge is to annotate the genome and characterize functional DNA elements. Transposable elements (TEs) are the most abundant component of plant genomes and can dramatically affect genome evolution and genetic variation. Thus, it is pivotal to identify TEs in the common bean genome. In this study, we performed a genome-wide transposon annotation in common bean using a combination of homology and sequence structure-based methods. We developed a 2.12-Mb transposon database which includes 791 representative transposon sequences and is available upon request or from www.phytozome.org. Of note, nearly all transposons in the database are previously unrecognized TEs. More than 5,000 transposon-related expressed sequence tags (ESTs) were detected which indicates that some transposons may be transcriptionally active. Two Ty1-copia retrotransposon families were found to encode the envelope-like protein which has rarely been identified in plant genomes. Also, we identified an extra open reading frame (ORF) termed ORF2 from 15 Ty3-gypsy families that was located between the ORF encoding the retrotransposase and the 3'LTR. The ORF2 was in opposite transcriptional orientation to retrotransposase. Sequence homology searches and phylogenetic analysis suggested that the ORF2 may have an ancient origin, but its function is not clear. These transposon data provide a useful resource for understanding the genome organization and evolution and may be used to identify active TEs for developing transposon-tagging system in common bean and other related genomes.

Simple sequence repeats amplification

The technique of SSR amplification is a prerequisite to generate the molecular profiles of various alleles of an individual or genotype. Amplification is the multifold duplication and accumulation of a targeted region which is achieved by polymerase chain reaction. It needs ingredients such as buffer, MgCl2, dNTPs, primers, and DNA polymerase enzyme. The utilization of these essential PCR components in optimal concentrations determines the success of amplification. Thus SSRs, as primers, play an important role in enhancing the amplification and thereby generating the genotype profile. With the advent of technology, fluorophore-labeled primers along with automated capillary electrophoresis system have enhanced the efficiency of detection.

Isolation and characterization of 13 new polymorphic microsatellite markers in the Phaseolus vulgaris L. (Common Bean) genome

In this study, 13 polymorphic microsatellite markers were isolated from the Phaseolus vulgaris L. (common bean) by using the Fast Isolation by AFLP of Sequence COntaining Repeats (FIASCO) protocol. These markers revealed two to seven alleles, with an average of 3.64 alleles per locus. The polymorphic information content (PIC) values ranged from 0.055 to 0.721 over 13 loci, with a mean value of 0.492, and 7 loci having PIC greater than 0.5. The expected heterozygosity (H_E) and observed heterozygosity (H_O) levels ranged from 0.057 to 0.814 and from 0.026 to 0.531, respectively. Cross-species amplification of the 13 prime pairs was performed in its related specie of Vigna unguiculata L. Seven out of all these markers showed cross-species transferability. These markers will be useful for future genetic diversity and population genetics studies for this agricultural specie and its related species.

Use of the advanced backcross-QTL method to transfer seed mineral accumulation nutrition traits from wild to Andean cultivated common beans

Iron deficiency anemia and zinc deficiency are major health concerns across the world and can be addressed by biofortification breeding of higher mineral concentration in staple crops, such as common bean. Wild common beans have for the most part had higher average seed mineral concentration than cultivars of this species but have small un-commercial seeds. A logical approach for the transfer of the seed mineral trait from wild beans to cultivated beans is through the advanced backcross breeding approach. The goal of this study was to analyze a population of 138 BC(2)F(3:5) introgression lines derived from the very high iron wild genotype G10022 backcrossed into the genetic background of the commercial-type variety 'Cerinza', a large-red seeded bush bean cultivar of the Andean genepool. In addition to measuring seed mineral accumulation traits and the quantitative trait loci (QTL) controlling these traits we were interested in simultaneously testing the adaptation of the introgression lines in two replicated yield trials. We found the cross to have high polymorphism and constructed an anchored microsatellite map for the population that was 1,554-cM long and covered all 11 linkage groups of the common bean genome. Through composite interval mapping (CIM) and single point analysis (SPA), we identified associations of markers and mineral traits on b01, b06, b07, b08, b10 and b11 for seed iron concentration, and markers on b01, b04 and b10 for seed zinc concentration. The b07 and b08 QTL aligned with previous QTL for iron concentration. A large number of QTL were found for seed weight (9 with CIM and 36 with SPA analysis) and correlations between seed size and mineral content affected the identification of iron and zinc contents' QTL on many linkage groups. Segregation distortion around domestication genes made some areas difficult to introgress. However, in conclusion, the advanced backcross program produced some introgression lines with high mineral accumulation traits using a wild donor parent.

Marker-based linkage map of Andean common bean (Phaseolus vulgaris L.) and mapping of QTLs underlying popping ability traits

BACKGROUND

Nuña bean is a type of ancient common bean (Phaseolus vulgaris L.) native to the Andean region of South America, whose seeds possess the unusual property of popping. The nutritional features of popped seeds make them a healthy low fat and high protein snack. However, flowering of nuña bean only takes place under short-day photoperiod conditions, which means a difficulty to extend production to areas where such conditions do not prevail. Therefore, breeding programs of adaptation traits will facilitate the diversification of the bean crops and the development of new varieties with enhanced healthy properties. Although the popping trait has been profusely studied in maize (popcorn), little is known about the biology and genetic basis of the popping ability in common bean. To obtain insights into the genetics of popping ability related traits of nuña bean, a comprehensive quantitative trait loci (QTL) analysis was performed to detect single-locus and epistatic QTLs responsible for the phenotypic variance observed in these traits.

RESULTS

A mapping population of 185 recombinant inbred lines (RILs) derived from a cross between two Andean common bean genotypes was evaluated for three popping related traits, popping dimension index (PDI), expansion coefficient (EC), and percentage of unpopped seeds (PUS), in five different environmental conditions. The genetic map constructed included 193 loci across 12 linkage groups (LGs), covering a genetic distance of 822.1 cM, with an average of 4.3 cM per marker. Individual and multi-environment QTL analyses detected a total of nineteen single-locus QTLs, highlighting among them the co-localized QTLs for the three popping ability traits placed on LGs 3, 5, 6, and 7, which together explained 24.9, 14.5, and 25.3% of the phenotypic variance for PDI, EC, and PUS, respectively. Interestingly, epistatic interactions among QTLs have been detected, which could have a key role in the genetic control of popping.

CONCLUSIONS

The QTLs here reported constitute useful tools for marker assisted selection breeding programs aimed at improving nuña bean cultivars, as well as for extending our knowledge of the genetic determinants and genotype x environment interaction involved in the popping ability traits of this bean crop.

Characterization and compilation of polymorphic simple sequence repeat (SSR) markers of peanut from public database

Background

There are several reports describing thousands of SSR markers in the peanut (Arachis hypogaea L.) genome. There is a need to integrate various research reports of peanut DNA polymorphism into a single platform. Further, because of lack of uniformity in the labeling of these markers across the publications, there is some confusion on the identities of many markers. We describe below an effort to develop a central comprehensive database of polymorphic SSR markers in peanut.

Findings

We compiled 1,343 SSR markers as detecting polymorphism (14.5%) within a total of 9,274 markers. Amongst all polymorphic SSRs examined, we found that AG motif (36.5%) was the most abundant followed by AAG (12.1%), AAT (10.9%), and AT (10.3%).The mean length of SSR repeats in dinucleotide SSRs was significantly longer than that in trinucleotide SSRs. Dinucleotide SSRs showed higher polymorphism frequency for genomic SSRs when compared to trinucleotide SSRs, while for EST-SSRs, the frequency of polymorphic SSRs was higher in trinucleotide SSRs than in dinucleotide SSRs. The correlation of the length of SSR and the frequency of polymorphism revealed that the frequency of polymorphism was decreased as motif repeat number increased.

Conclusions

The assembled polymorphic SSRs would enhance the density of the existing genetic maps of peanut, which could also be a useful source of DNA markers suitable for high-throughput QTL mapping and marker-assisted selection in peanut improvement and thus would be of value to breeders.

Identification and analysis of common bean (Phaseolus vulgaris L.) transcriptomes by massively parallel pyrosequencing

BACKGROUND

Common bean (Phaseolus vulgaris) is the most important food legume in the world. Although this crop is very important to both the developed and developing world as a means of dietary protein supply, resources available in common bean are limited. Global transcriptome analysis is important to better understand gene expression, genetic variation, and gene structure annotation in addition to other important features. However, the number and description of common bean sequences are very limited, which greatly inhibits genome and transcriptome research. Here we used 454 pyrosequencing to obtain a substantial transcriptome dataset for common bean.

RESULTS

We obtained 1,692,972 reads with an average read length of 207 nucleotides (nt). These reads were assembled into 59,295 unigenes including 39,572 contigs and 19,723 singletons, in addition to 35,328 singletons less than 100 bp. Comparing the unigenes to common bean ESTs deposited in GenBank, we found that 53.40% or 31,664 of these unigenes had no matches to this dataset and can be considered as new common bean transcripts. Functional annotation of the unigenes carried out by Gene Ontology assignments from hits to Arabidopsis and soybean indicated coverage of a broad range of GO categories. The common bean unigenes were also compared to the bean bacterial artificial chromosome (BAC) end sequences, and a total of 21% of the unigenes (12,724) including 9,199 contigs and 3,256 singletons match to the 8,823 BAC-end sequences. In addition, a large number of simple sequence repeats (SSRs) and transcription factors were also identified in this study.

CONCLUSIONS

This work provides the first large scale identification of the common bean transcriptome derived by 454 pyrosequencing. This research has resulted in a 150% increase in the number of Phaseolus vulgaris ESTs. The dataset obtained through this analysis will provide a platform for functional genomics in common bean and related legumes and will aid in the development of molecular markers that can be used for tagging genes of interest. Additionally, these sequences will provide a means for better annotation of the on-going common bean whole genome sequencing.