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

Genetics, genomics, and breeding of black gram [Vigna mungo (L.) Hepper]

Black gram [Vigna mungo (L.) Hepper] is a highly nutritious grain legume crop, mainly grown in South and Southeast Asia, with the largest area in India, where the crop is challenged by several biotic and abiotic stresses leading to significant yield losses. Improving genetic gains to increase on-farm yields is the primary goal of black gram breeding programs. This could be achieved by developing varieties resistant to major diseases like mungbean yellow mosaic disease, urdbean leaf crinkle virus, Cercospora leaf spot, anthracnose, powdery mildew, and insect pests such as whitefly, cowpea aphids, thrips, stem flies, and bruchids. Along with increasing on-farm yields, incorporating market-preferred traits ensures the adoption of improved varieties. Black gram breeding programs rely upon a limited number of parental lines, leading to a narrow genetic base of the developed varieties. For accelerating genetic gain, there is an urgent need to include more diverse genetic material for improving traits for better adaptability and stress resistance in breeding populations. The present review summarizes the importance of black gram, the major biotic and abiotic stresses, available genetic and genomic resources, major traits for potential crop improvement, their inheritance, and the breeding approaches being used in black gram for the development of new varieties.

Genome-wide identification, characterization and primer designing of simple sequence repeats across Leguminosae family

Legumes are important clade of commercially important family Leguminosae that mainly include medicinal, flowering and edible plants. Although the genomic sequence of legumes is accessible, only the limited number of effective simple sequence repeat markers has been identified by prior research. Additional polymorphic simple sequence repeats marker discovery will aid in the genetics and breeding of legumes. In this study, 13 complete genome sequences were screened for the identification of chromosome-wise simple sequence repeats (SSRs) and 1,866,861 SSRs were identified. Based on the study, it was observed that the number of SSRs in non-coding region was more as compared to coding region and frequency of mononucleotides was highest followed by di-nucleotides while penta- and hexa-nucleotide repeats were least frequent one. The identified genome-wide SSRs and newly developed SSR markers, primers and their mapping will provide a powerful means for genetic researches across Leguminosae plants, including genetic diversity and evolutionary origin analysis, fingerprinting, QTL mapping and marker-assisted selection for breeding as well as comparative genomic analysis studies.

Unique Salt-Tolerance-Related QTLs, Evolved in Vigna riukiuensis (Na+ Includer) and V. nakashimae (Na+ Excluder), Shed Light on the Development of Super-Salt-Tolerant Azuki Bean (V. angularis) Cultivars

Wild relatives of crops have the potential to improve food crops, especially in terms of improving abiotic stress tolerance. Two closely related wild species of the traditional East Asian legume crops, Azuki bean (Vigna angularis), V. riukiuensis "Tojinbaka" and V. nakashimae "Ukushima" were shown to have much higher levels of salt tolerance than azuki beans. To identify the genomic regions responsible for salt tolerance in "Tojinbaka" and "Ukushima", three interspecific hybrids were developed: (A) azuki bean cultivar "Kyoto Dainagon" × "Tojinbaka", (B) "Kyoto Dainagon" × "Ukushima" and (C) "Ukushima" × "Tojinbaka". Linkage maps were developed using SSR or restriction-site-associated DNA markers. There were three QTLs for "percentage of wilt leaves" in populations A, B and C, while populations A and B had three QTLs and population C had two QTLs for "days to wilt". In population C, four QTLs were detected for Na⁺ concentration in the primary leaf. Among the F₂ individuals in population C, 24% showed higher salt tolerance than both wild parents, suggesting that the salt tolerance of azuki beans can be further improved by combining the QTL alleles of the two wild relatives. The marker information would facilitate the transfer of salt tolerance alleles from "Tojinbaka" and "Ukushima" to azuki beans.

De Novo Transcriptome Assembly and EST-SSR Marker Development and Application in Chrysosplenium macrophyllum

Chrysosplenium macrophyllum Oliv., belonging to the family Saxifragaceae, is a traditional and unique Chinese herbal medicine. However, the lack of adequate molecular markers has hampered the progress regarding population genetics and evolution within this species. In this research, we used the DNBSEQ-T7 Sequencer (MGI) sequencing assay to analyze the transcriptome profiles of C. macrophyllum. SSR markers were developed on the basis of transcriptomic sequences and further validated on C. macrophyllum and other Chrysosplenium species. The genetic diversity and structure of the 12 populations were analyzed by using polymorphic expressed sequence tag simple sequence repeat (EST-SSR) markers. A potential pool of 3127 non-redundant EST-SSR markers were identified for C. macrophyllum in this study. The developed EST-SSR markers had high amplification rates and cross-species transferability in Chrysosplenium. Our results also showed that the natural populations of C. macrophyllum had a high level of genetic diversity. Genetic distance, principal component analysis, and popular structure analysis revealed that all 60 samples clustered into two major groups that were consistent with their geographical origins. This study provided a batch of highly polymorphic EST-SSR molecular markers that were developed via transcriptome sequencing. These markers will be of great significance for the study of the genetic diversity and evolutionary history of C. macrophyllum and other Chrysosplenium species.

RNA-Seq and genetic diversity analysis of faba bean (Vicia faba L.) varieties in China

Background

Faba bean (Vicia faba L) is one of the most important legumes in the world. However, there is relatively little genomic information available for this species owing to its large genome. The lack of data impedes the discovery of molecular markers and subsequent genetic research in faba bean. The objective of this study was to analyze the faba bean transcriptome, and to develop simple sequence repeat (SSR) markers to determine the genetic diversity of 226 faba bean varieties derived from different regions in China.

Methods

Faba bean varieties with different phenotype were used in transcriptome analysis. The functions of the unigenes were analyzed using various database. SSR markers were developed and the polymorphic markers were selected to conduct genetic diversity analysis.

Results

A total of 92.43 Gb of sequencing data was obtained in this study, and 133,487 unigene sequences with a total length of 178,152,541 bp were assembled. A total of 5,200 SSR markers were developed on the basis of RNA-Seq analysis. Then, 200 SSR markers were used to evaluate polymorphisms. In total, 103 (51.5%) SSR markers showed significant and repeatable bands between different faba bean varieties. Clustering analysis revealed that 226 faba bean materials were divided into five groups. Genetic diversity analysis revealed that the relationship between different faba beans in China was related, especially in the same region. These results provided a valuable data resource for annotating genes to different categories and developing SSR markers.

Progress of Genomics-Driven Approaches for Sustaining Underutilized Legume Crops in the Post-Genomic Era

Legume crops, belonging to the Fabaceae family, are of immense importance for sustaining global food security. Many legumes are profitable crops for smallholder farmers due to their unique ability to fix atmospheric nitrogen and their intrinsic ability to thrive on marginal land with minimum inputs and low cultivation costs. Recent progress in genomics shows promise for future genetic gains in major grain legumes. Still it remains limited in minor legumes/underutilized legumes, including adzuki bean, cluster bean, horse gram, lathyrus, red clover, urd bean, and winged bean. In the last decade, unprecedented progress in completing genome assemblies of various legume crops and resequencing efforts of large germplasm collections has helped to identify the underlying gene(s) for various traits of breeding importance for enhancing genetic gain and contributing to developing climate-resilient cultivars. This review discusses the progress of genomic resource development, including genome-wide molecular markers, key breakthroughs in genome sequencing, genetic linkage maps, and trait mapping for facilitating yield improvement in underutilized legumes. We focus on 1) the progress in genomic-assisted breeding, 2) the role of whole-genome resequencing, pangenomes for underpinning the novel genomic variants underlying trait gene(s), 3) how adaptive traits of wild underutilized legumes could be harnessed to develop climate-resilient cultivars, 4) the progress and status of functional genomics resources, deciphering the underlying trait candidate genes with putative function in underutilized legumes 5) and prospects of novel breeding technologies, such as speed breeding, genomic selection, and genome editing. We conclude the review by discussing the scope for genomic resources developed in underutilized legumes to enhance their production and play a critical role in achieving the "zero hunger" sustainable development goal by 2030 set by the United Nations.

Transcriptome Analysis Reveals Key Pathways and Candidate Genes Controlling Seed Development and Size in Ricebean (Vigna umbellata)

Ricebean (Vigna umbellata) is a lesser known pulse with well-recognized potential. Recently, it has emerged as a legume with endowed nutritional potential because of high concentration of quality protein and other vital nutrients in its seeds. However, the genes and pathways involved in regulating seed development and size are not understood in this crop. In our study, we analyzed the transcriptome of two genotypes with contrasting grain size (IC426787: large seeded and IC552985: small seeded) at two different time points, namely, 5 and 10 days post-anthesis (DPA). The bold seeded genotype across the time points (B5_B10) revealed 6,928 differentially expressed genes (DEGs), whereas the small seeded genotype across the time point (S5_S10) contributed to 14,544 DEGs. We have also identified several candidate genes for seed development–related traits like seed size and 100-seed weight. On the basis of similarity search and domain analysis, some candidate genes (PHO1, cytokinin dehydrogenase, A-type cytokinin, and ARR response negative regulator) related to 100-seed weight and seed size showed downregulation in the small seeded genotype. The MapMan and KEGG analysis confirmed that auxin and cytokinin pathways varied in both the contrasting genotypes and can therefore be the regulators of the seed size and other seed development–related traits in ricebeans. A total of 51 genes encoding SCF^TIR1/AFB, Aux/IAA, ARFs, E3 ubiquitin transferase enzyme, and 26S proteasome showing distinct expression dynamics in bold and small genotypes were also identified. We have also validated randomly selected SSR markers in eight accessions of the Vigna species (V. umbellata: 6; Vigna radiata: 1; and Vigna mungo: 1). Cross-species transferability pattern of ricebean–derived SSR markers was higher in V. radiata (73.08%) than V. mungo (50%). To the best of our knowledge, this is the first transcriptomic study conducted in this crop to understand the molecular basis of any trait. It would provide us a comprehensive understanding of the complex transcriptome dynamics during the seed development and gene regulatory mechanism of the seed size determination in ricebeans.

Draft genome sequence of the pulse crop blackgram [Vigna mungo (L.) Hepper] reveals potential R-genes

Blackgram [Vigna mungo (L.) Hepper] (2n = 2x = 22), an important Asiatic legume crop, is a major source of dietary protein for the predominantly vegetarian population. Here we construct a draft genome sequence of blackgram, for the first time, by employing hybrid genome assembly with Illumina reads and third generation Oxford Nanopore sequencing technology. The final de novo whole genome of blackgram is ~ 475 Mb (82% of the genome) and has maximum scaffold length of 6.3 Mb with scaffold N50 of 1.42 Mb. Genome analysis identified 42,115 genes with mean coding sequence length of 1131 bp. Around 80.6% of predicted genes were annotated. Nearly half of the assembled sequence is composed of repetitive elements with retrotransposons as major (47.3% of genome) transposable elements, whereas, DNA transposons made up only 2.29% of the genome. A total of 166,014 SSRs, including 65,180 compound SSRs, were identified and primer pairs for 34,816 SSRs were designed. Out of the 33,959 proteins, 1659 proteins showed presence of R-gene related domains. KIN class was found in majority of the proteins (905) followed by RLK (239) and RLP (188). The genome sequence of blackgram will facilitate identification of agronomically important genes and accelerate the genetic improvement of blackgram.

A chromosome-scale assembly of the black gram (Vigna mungo) genome

Black gram (Vigna mungo) is an important short duration grain legume crop. Black gram seeds provide an inexpensive source of dietary protein. Here, we applied the 10X Genomics linked-read technology to obtain a de novo whole genome assembly of V. mungo cultivated variety Chai Nat 80 (CN80). The preliminary assembly contained 12,228 contigs and had an N50 length of 5.2 Mb. Subsequent scaffolding using the long-range Chicago and HiC techniques yielded the first high-quality, chromosome-level assembly of 499 Mb comprising 11 pseudomolecules. Comparative genomics analyses based on sequence information from single-copy orthologous genes revealed that black gram and mungbean (Vigna radiata) diverged about 2.7 million years ago . The transversion rate (4DTv) analysis in V. mungo revealed no evidence supporting a recent genome-wide duplication event observed in the tetraploid créole bean (Vigna reflexo-pilosa). The proportion of repetitive elements in the black gram genome is slightly lower than the numbers reported for related Vigna species. The majority of long terminal repeat retrotransposons appeared to integrate into the genome within the last five million years. We also examined alternative splicing events in V. mungo using full-length transcript sequences. While intron retention was the most prevalent mode of alternative splicing in several plant species, alternative 3' acceptor site selection represented the majority of events in black gram. Our high-quality genome assembly along with the genomic variation information from the germplasm provides valuable resources for accelerating the development of elite varieties through marker-assisted breeding and for future comparative genomics and phylogenetic studies in legume species.

De novo assembly and characterization of transcriptome towards understanding molecular mechanism associated with MYMIV-resistance in Vigna mungo - A computational study.. [DOI: 10.1101/844639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The fast climate change affects yield in Vigna mungo via enhancing both biotic and abiotic stresses. Out of all factors, the yellow mosaic disease has the most damaging effect. However, due to lack of reference genome of Vigna mungo, the complete mechanism associated with MYMIV (Mungbean Yellow Mosaic Indian Virus) resistance in Vigna mungo remain elusive to date. Considering this, the authors made an attempt to release new transcriptome and its annotation by employing computational approaches. Quality assessment of the generated transcriptomes reveals that it successfully aligned with 99.03% of the raw reads and hence can be employed for future research. Functional annotation of the transcriptome reveals that 31% and ∼14% of the total transcripts encode lncRNAs and protein-coding sequences, respectively. Further, analysis reveals that, out of total transcripts, only 4536 and 78808 are significantly down and up-regulated during MYMIV infection in Vigna mungo, respectively. These significant transcripts are mainly associated with ribosome, spliceosome, glycolysis /gluconeogenesis, RNA transport, oxidative phosphorylation, protein processing in the endoplasmic reticulum, MAPK signaling pathway - plant, methionine and cysteine metabolism, purine metabolism and RNA degradation. Unlike the previous study, this is for the first time, the present study identified these pathways may play key role in MYMIV resistance in Vigna mungo. Thus, information and transcriptomes data available in the present study make a significant contribution to understanding the genomic structure of Vigna mungo, enabling future analyses as well as downstream applications of gene expression, sequence evolution, and genome annotation.

De novo assembly and characterization of transcriptome towards understanding molecular mechanism associated with MYMIV-resistance in Vigna mungo - A computational study.. [DOI: 10.1101/844639]

The fast climate change affects yield in Vigna mungo via enhancing both biotic and abiotic stresses. Out of all factors, the yellow mosaic disease has the most damaging effect. However, due to lack of reference genome of Vigna mungo, the complete mechanism associated with MYMIV (Mungbean Yellow Mosaic Indian Virus) resistance in Vigna mungo remain elusive to date. Considering this, the authors made an attempt to release new transcriptome and its annotation by employing computational approaches. Quality assessment of the generated transcriptomes reveals that it successfully aligned with 99.03% of the raw reads and hence can be employed for future research. Functional annotation of the transcriptome reveals that 31% and ∼14% of the total transcripts encode lncRNAs and protein-coding sequences, respectively. Further, analysis reveals that, out of total transcripts, only 4536 and 78808 are significantly down and up-regulated during MYMIV infection in Vigna mungo, respectively. These significant transcripts are mainly associated with ribosome, spliceosome, glycolysis /gluconeogenesis, RNA transport, oxidative phosphorylation, protein processing in the endoplasmic reticulum, MAPK signaling pathway - plant, methionine and cysteine metabolism, purine metabolism and RNA degradation. Unlike the previous study, this is for the first time, the present study identified these pathways may play key role in MYMIV resistance in Vigna mungo. Thus, information and transcriptomes data available in the present study make a significant contribution to understanding the genomic structure of Vigna mungo, enabling future analyses as well as downstream applications of gene expression, sequence evolution, and genome annotation.

De novo assembly and characterization of transcriptome towards understanding molecular mechanism associated with MYMIV-resistance in Vigna mungo - A computational study.. [DOI: 10.1101/844639] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The fast climate change affects yield in Vigna mungo via enhancing both biotic and abiotic stresses. Out of all factors, the yellow mosaic disease has the most damaging effect. However, due to lack of reference genome of Vigna mungo, the complete mechanism associated with MYMIV (Mungbean Yellow Mosaic Indian Virus) resistance in Vigna mungo remain elusive to date. Considering this, the authors made an attempt to release new transcriptome and its annotation by employing computational approaches. Quality assessment of the generated transcriptomes reveals that it successfully aligned with 99.03% of the raw reads and hence can be employed for future research. Functional annotation of the transcriptome reveals that 31% and ∼14% of the total transcripts encode lncRNAs and protein-coding sequences, respectively. Further, analysis reveals that, out of total transcripts, only 4536 and 78808 are significantly down and up-regulated during MYMIV infection in Vigna mungo, respectively. These significant transcripts are mainly associated with ribosome, spliceosome, glycolysis /gluconeogenesis, RNA transport, oxidative phosphorylation, protein processing in the endoplasmic reticulum, MAPK signaling pathway - plant, methionine and cysteine metabolism, purine metabolism and RNA degradation. Unlike the previous study, this is for the first time, the present study identified these pathways may play key role in MYMIV resistance in Vigna mungo. Thus, information and transcriptomes data available in the present study make a significant contribution to understanding the genomic structure of Vigna mungo, enabling future analyses as well as downstream applications of gene expression, sequence evolution, and genome annotation.

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.

Microsatellite Borders and Micro-sequence Conservation in Juglans

Walnuts (Juglans spp.) are economically important nut and timber species with a worldwide distribution. Using the published Persian walnut genome as a reference for the assembly of short reads from six Juglans species and several interspecific hybrids, we identified simple sequence repeats in 12 Juglans nuclear and organellar genomes. The genome-wide distribution and polymorphisms of nuclear and organellar microsatellites (SSRs) for most Juglans genomes have not been previously studied. We compared the frequency of nuclear SSR motifs and their lengths across Juglans, and identified section-specific chloroplast SSR motifs. Primer pairs were designed for more than 60,000 SSR-containing sequences based on alignment against assembled scaffold sequences. Of the >60,000 loci, 39,000 were validated by e-PCR using unique primer pairs. We identified primers containing 100% sequence identity in multiple species. Across species, sequence identity in the SSR-flanking regions was generally low. Although SSRs are common and highly dispersed in the genome, their flanking sequences are conserved at about 90 to 95% identity within Juglans and within species. In a few rare cases, flanking sequences are identical across species of Juglans. This comprehensive report of nuclear and organellar SSRs in Juglans and the generation of validated SSR primers will be a useful resource for future genetic analyses, walnut breeding programs, high-level taxonomic evaluations, and genomic studies in Juglandaceae.

Green gram and black gram: prospects of cultivation and breeding in Russian Federation

Diversifcation of crop production in the Russian Federation could be partly achieved by the introduction and production of minor and underutilized crops. Green gram or mung bean (Vigna radiata(L.) R. Wilczek) and black gram or urd (V. mungo(L.) Hepper) are grain legume crops cultivated in limited areas in the Russian Federation. Meanwhile, green gram occupies about 8.5 % of the world production area under pulses (without soybean). It is cultivated mainly in countries of Southeast Asia, but production is expanding to the entire subtropical belt of the globe. In our country these crops can be successfully grown on irrigation in a number of regions in the southern area of the European part and the Russian Far East, where the temperatures during their vegetation are about 28–30 °C and always above 15 °C. The purpose of this paper is to summarize the world’s experience in breeding improvement of mung bean and urd as crops with promise for cultivation in certain soil and climatic zones of the Russian Federation. The world production, use of these high-protein crops, history and peculiarities of breeding, including in the USSR, are covered. To expand the production of both crops in the Russian Federation, their popularization and development of breeding are required. Basic requirements for modern varieties include resistance to biotic and abiotic stressors which can be introgressed from wild relatives. The great importance of both crops in the Asian countries led to the rapid development of molecular researches there. The genome of black gram has been fully sequenced, the genome of green gram has been partly sequenced. Some genes and QTL of adaptability traits have been marked and mapped in a number of wild species of the genusVigna. The role of wild relatives in the breeding of crops concerned is discussed. In the world genebanks, signifcant genetic resources of mung bean and urd have been accumulated. All this creates prerequisites for the development of marker-assistant and genomic breeding.

VigSatDB: genome-wide microsatellite DNA marker database of three species of Vigna for germplasm characterization and improvement

Genus Vigna represented by more than 100 species is a source of nutritious edible seeds and sprouts that are rich sources of protein and dietary supplements. It is further valuable because of therapeutic attributes due to its antioxidant and anti-diabetic properties. A highly diverse and an extremely ecological niche of different species can be valuable genomic resources for productivity enhancement. It is one of the most underutilized crops for food security and animal feeds. In spite of huge species diversity, only three species of Vigna have been sequenced; thus, there is a need for molecular markers for the remaining species. Computational approach of microsatellite marker discovery along with evaluation of polymorphism utilizing available genomic data of different genotypes can be a quick and an economical approach for genomic resource development. Cross-species transferability by e-PCR over available genomes can further prioritize the potential SSR markers, which could be used for genetic diversity and population differentiation of the remaining species saving cost and time. We present VigSatDB—the world’s first comprehensive microsatellite database of genus Vigna, containing >875 K putative microsatellite markers with 772 354 simple and 103 865 compound markers mined from six genome assemblies of three Vigna species, namely, Vigna radiata (Mung bean), Vigna angularis (Adzuki bean) and Vigna unguiculata (Cowpea). It also contains 1976 validated published markers. Markers can be selected on the basis of chromosomes/location specificity, and primers can be generated using Primer3core tool integrated at backend. Efficacy of VigSatDB for microsatellite loci genotyping has been evaluated by 15 markers over a panel of 10 diverse genotype of V. radiata. Our web genomic resources can be used in diversity analysis, population and varietal differentiation, discovery of quantitative trait loci/genes, marker-assisted varietal improvement in endeavor of Vigna crop productivity and management.

Characterization of new microsatellite markers based on the transcriptome sequencing of Clematis finetiana

Background

Clematis is the biggest genus in the family Ranunculaceae with about 300 species. Clematis is also a globally important commercial group of flowers, especially in the United States and European countries. Their petals with different colors and shapes make the genus the “Queen of the Vines”. However, the genomic information and phylogeny of Clematis based on existing molecular studies are limited. In this paper, new microsatellites (SSR) markers were identified from the transcriptome data of C. finetiana obtained using the Illumina paired-end sequencing technology.

Results

Sequences on a total of 71,900 high-quality unigenes with the mean length of 865 bp were produced in this study. There were 6192unigenes annotated and classified into 49 functional sub-groups in three main ontology categories in GO (Gen Ontology) database,14,022 unigenes mapped to COGs (Clusters of Orthologous Groups) database and classified into 25 functional categories, and 21,494 unigenes obtained and divided into 128 pathways of KEGG (Kyoto Encyclopedia of Genes) Database. A total of 7532 SSRs were discovered from 6337 unigenes. We randomly tested 210 primer pairs, of which 52 primer pairs were able to generate specific products, and 19 possessed polymorphism in the 13 wild populations of six species from Clematis, which were used as a test material.

Conclusions

The dataset of C. finetiana transcriptome and the identified new SSR markers will promote genetic research and breeding effort in Clematis.

Electronic supplementary material

The online version of this article (10.1186/s41065-018-0060-x) contains supplementary material, which is available to authorized users.

Characterizing the transcriptome and microsatellite markers for almond (Amygdalus communis L.) using the Illumina sequencing platform

BACKGROUND

The almond tree (Prunus amygdalus Batsch) is an important nut tree grown in subtropical regions that produces nutrient-rich nuts. However, a paucity of genomic information and DNA markers has restricted the development of modern breeding technologies for almond trees.

RESULTS

In this study, almonds were sequenced with Illumina paired-end sequencing technology to obtain transcriptome data and develop simple sequence repeats (SSR) markers. We generated approximately 64 million clean reads from the various tissues of mixed almonds, and a total of 42,135 unigenes with an average length of 988 bp were obtained in the present study. A total of 27,586 unigenes (57.7% of all unigenes generated) were annotated using several databases. A total of 112,812 unigenes were annotated with the Gene Ontology (GO) database and assigned to 82 functional sub-groups, and 29,075 unigenes were assigned to the KOG database and classified into 25 function classifications. There were 9470 unigenes assigned to 129 Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways from five categories in the KEGG pathway database. We further identified 8641 SSR markers from 48,012 unigenes. A total of 100 SSR markers were randomly selected to validate quality, and 82 markers could amplify the specific products of A. communis L., whereas 70 markers were successfully transferable to five species (A. ledebouriana, A. mongolica, A. pedunculata, A. tangutica, and A. triloba).

CONCLUSIONS

Our study was the first to produce public transcriptome data from almonds. The development of SSR markers will promote genetics research and breeding programmes for almonds.

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 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.

Transcriptome sequencing and marker development in winged bean (Psophocarpus tetragonolobus; Leguminosae)

Winged bean, Psophocarpus tetragonolobus (L.) DC., is similar to soybean in yield and nutritional value but more viable in tropical conditions. Here, we strengthen genetic resources for this orphan crop by producing a de novo transcriptome assembly and annotation of two Sri Lankan accessions (denoted herein as CPP34 [PI 491423] and CPP37 [PI 639033]), developing simple sequence repeat (SSR) markers, and identifying single nucleotide polymorphisms (SNPs) between geographically separated genotypes. A combined assembly based on 804,757 reads from two accessions produced 16,115 contigs with an N50 of 889 bp, over 90% of which has significant sequence similarity to other legumes. Combining contigs with singletons produced 97,241 transcripts. We identified 12,956 SSRs, including 2,594 repeats for which primers were designed and 5,190 high-confidence SNPs between Sri Lankan and Nigerian genotypes. The transcriptomic data sets generated here provide new resources for gene discovery and marker development in this orphan crop, and will be vital for future plant breeding efforts. We also analyzed the soybean trypsin inhibitor (STI) gene family, important plant defense genes, in the context of related legumes and found evidence for radiation of the Kunitz trypsin inhibitor (KTI) gene family within winged bean.