Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea [Cajanus cajan (L.) Millspaugh]

The MYB Transcription Factor GmMYB78 Negatively Regulates Phytophthora sojae Resistance in Soybean

Phytophthora root rot is a devastating disease of soybean caused by Phytophthora sojae. However, the resistance mechanism is not yet clear. Our previous studies have shown that GmAP2 enhances sensitivity to P. sojae in soybean, and GmMYB78 is downregulated in the transcriptome analysis of GmAP2-overexpressing transgenic hairy roots. Here, GmMYB78 was significantly induced by P. sojae in susceptible soybean, and the overexpressing of GmMYB78 enhanced sensitivity to the pathogen, while silencing GmMYB78 enhances resistance to P. sojae, indicating that GmMYB78 is a negative regulator of P. sojae. Moreover, the jasmonic acid (JA) content and JA synthesis gene GmAOS1 was highly upregulated in GmMYB78-silencing roots and highly downregulated in overexpressing ones, suggesting that GmMYB78 could respond to P. sojae through the JA signaling pathway. Furthermore, the expression of several pathogenesis-related genes was significantly lower in GmMYB78-overexpressing roots and higher in GmMYB78-silencing ones. Additionally, we screened and identified the upstream regulator GmbHLH122 and downstream target gene GmbZIP25 of GmMYB78. GmbHLH122 was highly induced by P. sojae and could inhibit GmMYB78 expression in resistant soybean, and GmMYB78 was highly expressed to activate downstream target gene GmbZIP25 transcription in susceptible soybean. In conclusion, our data reveal that GmMYB78 triggers soybean sensitivity to P. sojae by inhibiting the JA signaling pathway and the expression of pathogenesis-related genes or through the effects of the GmbHLH122-GmMYB78-GmbZIP25 cascade pathway.

De novo transcriptome assembly of Dalbergia sissoo Roxb. (Fabaceae) under Botryodiplodia theobromae-induced dieback disease

Dalbergia sissoo Roxb. (Shisham) is a timber-producing species of economic, cultural, and medicinal importance in the Indian subcontinent. In the past few decades, Shisham's dieback disease caused by the fungus Botryodiplodia theobromae has become an evolving issue in the subcontinent endangering its survival. To gain insights into this issue, a standard transcriptome assembly was deployed to assess the response of D. sissoo at the transcriptomic level under the stress of B. theobromae infection. For RNA isolation, the control and infected leaf tissue samples were taken from 1-year-old greenhouse-grown D. sissoo plants after 20 days of stem-base spore inoculation. cDNA synthesis was performed from these freshly isolated RNA samples that were then sent for sequencing. About 18.14 Gb (Giga base) of data was generated using the BGISEQ-500 sequencing platform. In terms of Unigenes, 513,821 were identified after a combined assembly of all samples and then filtering the abundance. The total length of Unigenes, their average length, N50, and GC-content were 310,523,693 bp, 604 bp, 1,101 bp, and 39.95% respectively. The Unigenes were annotated using 7 functional databases i.e., 200,355 (NR: 38.99%), 164,973 (NT: 32.11%), 123,733 (Swissprot: 24.08%), 142,580 (KOG: 27.75%), 139,588 (KEGG: 27.17%), 99,752 (GO: 19.41%), and 137,281 (InterPro: 26.72%). Furthermore, the Transdecoder detected 115,762 CDS. In terms of SSR (Simple Sequence Repeat) markers, 62,863 of them were distributed on 51,508 Unigenes and on the predicted 4673 TF (Transcription Factor) coding Unigenes. A total of 16,018 up- and 19,530 down-regulated Differentially Expressed Genes (DEGs) were also identified. Moreover, the Plant Resistance Genes (PRGs) had a count of 9230. We are hopeful that in the future, these identified Unigenes, SSR markers, DEGs and PRGs will provide the prerequisites for managing Shisham dieback disease, its breeding, and in tree improvement programs.

Development of simple sequence repeat markers for sugarcane from data mining of expressed sequence tags

Sugarcane (Saccharum spp. hybrids) is a worldwide acclaimed important agricultural crop used primarily for sugar production and biofuel. Sugarcane's genetic complexity, aneuploidy, and extreme heterozygosity make it a challenging crop in developing improved varieties. The molecular breeding programs promise to develop nutritionally improved varieties for both direct consumption and commercial application. Therefore, to address these challenges, the development of simple sequence repeats (SSRs) has been proven to be a powerful molecular tool in sugarcane. This study involved the collection of 285216 expressed sequence tags (ESTs) from sugarcane, resulting in 23666 unigenes, including 4547 contigs. Our analysis identified 4120 unigenes containing a total of 4960 SSRs, with the most abundant repeat types being monomeric (44.33%), dimeric (13.10%), and trimeric (39.68%). We further chose 173 primers to analyze the banding pattern in 10 sugarcane accessions by PAGE analysis. Additionally, functional annotation analysis showed that 71.07%, 53.6%, and 10.3% unigenes were annotated by Uniport, GO, and KEGG, respectively. GO annotations and KEGG pathways were distributed across three functional categories: molecular (46.46%), cellular (33.94%), and biological pathways (19.6%). The cluster analysis indicated the formation of four distinct clusters among selected sugarcane accessions, with maximum genetic distance observed among the varieties. We believe that these EST-SSR markers will serve as valuable references for future genetic characterization, species identification, and breeding efforts in sugarcane.

De Novo Transcriptome Profiling for the Generation and Validation of Microsatellite Markers, Transcription Factors, and Database Development for Andrographis paniculata

Andrographis paniculata belongs to the family Acanthaceae and is known for its medicinal properties owing to the presence of unique constituents belonging to the lactones, diterpenoids, diterpene glycosides, flavonoids, and flavonoid glycosides groups of chemicals. Andrographolide, a major therapeutic constituent of A. paniculata, is extracted primarily from the leaves of this plant and exhibits antimicrobial and anti-inflammatory activities. Using 454 GS-FLX pyrosequencing, we have generated a whole transcriptome profile of entire leaves of A. paniculata. A total of 22,402 high-quality transcripts were generated, with an average transcript length and N50 of 884 bp and 1007 bp, respectively. Functional annotation revealed that 19,264 (86%) of the total transcripts showed significant similarity with the NCBI-Nr database and were successfully annotated. Out of the 19,264 BLAST hits, 17,623 transcripts were assigned GO terms and distributed into three major functional categories: molecular function (44.62%), biological processes (29.19%), and cellular component (26.18%) based on BLAST2GO. Transcription factor analysis showed 6669 transcripts, belonging to 57 different transcription factor families. Fifteen TF genes that belong to the NAC, MYB, and bHLH TF categories were validated by RT PCR amplification. In silico analysis of gene families involved in the synthesis of biochemical compounds having medicinal values, such as cytochrome p450, protein kinases, heat shock proteins, and transporters, was completed and a total of 102 different transcripts encoding enzymes involved in the biosynthesis of terpenoids were predicted. Out of these, 33 transcripts belonged to terpenoid backbone biosynthesis. This study also identified 4254 EST-SSRs from 3661 transcripts, representing 16.34% of the total transcripts. Fifty-three novel EST-SSR markers generated from our EST dataset were used to assess the genetic diversity among eighteen A. paniculata accessions. The genetic diversity analysis revealed two distinct sub-clusters and all accessions based on the genetic similarity index were distinct from each other. A database based on EST transcripts, EST-SSR markers, and transcription factors has been developed using data generated from the present study combined with available transcriptomic resources from a public database using Meta transcriptome analysis to make genomic resources available in one place to the researchers working on this medicinal plant.

Small cardamom genome: development and utilization of microsatellite markers from a draft genome sequence of Elettaria cardamomum Maton

Small cardamom (Elettaria cardamomum Maton), the queen of spices, is the third most expensive spice in the world after saffron and vanilla, valued highly for its aroma and taste. This perennial herbaceous plant is a native of coastal parts of Southern India and displays a significant amount of morphological diversity. Its genetic potential has not been exploited due to lack of genomic resources limiting our understanding of the genome and important metabolic pathways which give it the economic advantage in the spice industry. Here, we report upon the de novo assembled, draft whole genome sequence of cardamom variety, Njallani Green Gold. We used a hybrid assembly strategy using the reads from the Oxford Nanopore, Illumina and 10x Genomics GemCode sequencing chemistries. The assembled genome length was 1.06 Gb (gigabases) which is close to the estimated genome size of cardamom. More than 75% of the genome was captured in 8000 scaffolds with a N50 of 0.15 Mb. The genome appears to have a high repeat content and 68055 gene models were predicted. The genome is close to Musa species and displays an expansion and contraction in different gene families. The draft assembly was used for in silico mining of simple sequence repeats (SSRs). A total of 2,50,571 SSRs were identified of which 2,18,270 were perfect SSRs and 32,301 were compound SSRs. Among the perfect SSRs, trinucleotides were most abundant (1,25,329) and hexanucleotide repeats appear least (2,380). From the 2,50,571 SSRs mined, 2,27,808 primer pairs were designed based on flanking sequence information. Wet lab validation was performed for 246 SSR loci and based on their amplification profiles, 60 SSR markers were used for diversity analysis of a set of 60 diverse cardamom accessions. The average number of alleles detected per locus were 14.57 with a minimum of 4 and maximum of 30 alleles. Population structure analysis revealed the presence of high degree of admixtures which could primarily be due to cross-pollination prevalent in this species. The SSR markers identified would help in the development of gene or trait-linked markers which can be subsequently used for marker-assisted breeding for crop improvement in cardamom. The information on utilization of the SSR loci for generation of markers has been developed into a public database, 'cardamomSSRdb' that is freely available for use by the cardamom community.

cgSSR marker-based genome-wide association study identified genomic regions for panicle characters and yield in rice (Oryza sativa L.)

BACKGROUND

To improve production efficiency, positive alleles corresponding to yield-related attributes must be accumulated in a single elite background. We designed and used cgSSR markers, which are superior to random SSR markers in genome-wide association study, to identify genomic regions that contribute to panicle characters and grain yield in this study.

RESULTS

As evidenced by the high polymorphic information content value and gene diversity coefficient, the new cgSSR markers were determined to be highly informative. These cgSSR markers were employed to generate genotype data for an association panel evaluated for four panicle characters and grain yield over three seasons. For five traits, 17 significant marker-trait associations on six chromosomes were discovered. The percentage of phenotypic variance that could be explained ranged from 4% to 13%. Unrelated gene-derived markers had a strong association with target traits as well.

CONCLUSION

Trait-associated cgSSR markers derived from corresponding or related genes ensure their utility in direct allele selection, while other linked markers aid in allele selection indirectly by altering the phenotype of interest. Through a marker-assisted breeding approach, these marker-trait associations can be leveraged to accumulate favourable alleles for yield enhancement in rice. © 2022 Society of Chemical Industry.

Occurrence of simple sequence repeats in cDNA sequences of safflower (Carthamus tinctorius) reveals the importance of SSR-containing genes for cell biology and dynamic response to environmental cues

Safflower (Carthamus tinctorius) is a diploid crop plant belonging to the family Asteraceae and is well known as one of important oilseed crops due to edible oil containing unsaturated fatty acids. In recent years it is gaining increased attention for food, pharmaceutical and industrial uses, and hence the updating its breeding methods is necessary. Genic simple sequence repeats (SSRs) in addition of being desire molecular markers, are supposed to influence gene function and the respective phenotype. This study aimed to identify SSRs in cDNA sequences and further analysis of the functional features of the SSR-containing genes to elucidate their role in biological and cellular processes. We identified 1,841 SSR regions in 1,667 cDNA sequences. Among all types of repeats, trinucleotide repeats were the most abundant (35.7%), followed by hexanucleotide (29.6%) and dinucleotide repeats (22.0%). Thirty five SSR primer pairs were validated by PCR reaction, detected a high rate of polymorphism (>57%) among safflower accessions, physically mapped on safflower genome and could clearly discriminate the cultivated accessions from wild relatives. The cDNA-derived SSR markers are suitable for evaluation of genetic diversity, linkage and association mapping studies and genome-based breeding programmes. Occurrence of SSR repeats in biologically-important classes of proteins such as kinases, transferases and transcription factors was inferred from functional analyses, which along with variability of their repeat copies, can endow the cell and whole organism the flexibility of facing with continuously changing environment, and indicate a structure-based evolution mechanism of the genome which acts as an up-to-dating tool for the cell and whole origanism, which is realized in GO terms such as involvement of most SSR-containing genes in biological, cellular and metabolic processes, especially in response to stimulus, response to stress, interaction to other organisms and defense responses.

Association mapping for abiotic stress tolerance using heat- and drought-related syntenic markers in okra

BACKGROUND

Considerable production losses are caused by heat and drought stress in okra. Germplasm evaluation at genetic level is essential for the selection of promising genotypes. Lack of genomic information of okra limits the use of genetic markers. However, syntenic markers of some related family could be used for molecular characterization of major economic traits.

METHODS AND RESULTS

Herein, 56 okra genotypes were evaluated for drought and heat tolerance. Sixty-one expressed sequence tags (ESTs) identified for heat and drought tolerance in cotton were searched from literature surveys and databases. The identified ESTs were BLAST searched into okra unigene database. Primers of selected okra unigenes were synthesized and amplified in all genotypes using standard polymerase chain reaction (PCR) protocol. Marker trait association (MTA) of the syntenic unigenes were identified between genotypic and phenotypic data on the basis of linkage disequilibrium Functional syntenic analysis revealed that out of these 61 cotton ESTs 55 had functional homology with okra unigenes. These 55 unigenes were used as markers for further analysis (amplification). Okra genotypes showed significance variations for all the physo-morphological parameters under heat and drought stress. Genotypes Perbhani Karanti, IQRA-III, Selection Super Green, Anmol and Line Bourd performed better under drought stress whereas genotypes Perbhani Karanti, IQRA-III, Green Gold, OK-1501 and Selection Super Green showed heat tolerance. Fifty markers showed amplification in okra. Fifty-six okra genotypes were clustered into three distinct populations. LD analysis has shown most significant linkage between markers Unigene43786 and Unigene3662. MTAs using MLM and GLM models revealed that 23 markers have significant associations (p < 0.05) with different traits under control and stressed conditions. Relative water content is associated with four markers (Unigene10673, Unigene99547, Unigene152901, and Unigene129684) under drought conditions. Whereas, Electrolyte leakage was associated with 3 markers (Unigene109922, Unigene28667 and Unigene146907) under heat stress.

CONCLUSION

These identified unigenes may be helpful in the development of drought and heat tolerant genotypes in okra.

De novo Transcriptome Analysis of Drought-Adapted Cluster Bean (Cultivar RGC-1025) Reveals the Wax Regulatory Genes Involved in Drought Resistance

Cluster bean (Cyamopsis tetragonoloba L.) is one of the multipurpose underexplored crops grown as green vegetable and for gum production in dryland areas. Cluster bean is known as relatively tolerant to drought and salinity stress. To elucidate the molecular mechanisms involved in the drought tolerance of cluster bean cultivar RGC-1025, RNA sequencing (RNA-seq) of the drought-stressed and control samples was performed. De novo assembly of the reads resulted in 66,838 transcripts involving 203 pathways. Among these transcripts, differentially expressed gene (DEG) analysis resulted in some of the drought-responsive genes expressing alpha dioxygenase 2, low temperature-induced 65 kDa protein (LDI65), putative vacuolar amino acid transporter, and late embryogenesis abundant protein (LEA 3). The analysis also reported drought-responsive transcription factors (TFs), such as NAC, WRKY, GRAS, and MYB families. The relative expression of genes by qRT-PCR revealed consistency with the DEG analysis. Key genes involved in the wax biosynthesis pathway were mapped using the DEG data analysis. These results were positively correlated with epicuticular wax content and the wax depositions on the leaf surfaces, as evidenced by scanning electron microscope (SEM) image analysis. Further, these findings support the fact that enhanced wax deposits on the leaf surface had played a crucial role in combating the drought stress in cluster beans under drought stress conditions. In addition, this study provided a set of unknown genes and TFs that could be a source of engineering tolerance against drought stress in cluster beans.

Single trait versus principal component based association analysis for flowering related traits in pigeonpea

Pigeonpea, a tropical photosensitive crop, harbors significant diversity for days to flowering, but little is known about the genes that govern these differences. Our goal in the current study was to use genome wide association strategy to discover the loci that regulate days to flowering in pigeonpea. A single trait as well as a principal component based association study was conducted on a diverse collection of 142 pigeonpea lines for days to first and fifty percent of flowering over 3 years, besides plant height and number of seeds per pod. The analysis used seven association mapping models (GLM, MLM, MLMM, CMLM, EMLM, FarmCPU and SUPER) and further comparison revealed that FarmCPU is more robust in controlling both false positives and negatives as it incorporates multiple markers as covariates to eliminate confounding between testing marker and kinship. Cumulatively, a set of 22 SNPs were found to be associated with either days to first flowering (DOF), days to fifty percent flowering (DFF) or both, of which 15 were unique to trait based, 4 to PC based GWAS while 3 were shared by both. Because PC1 represents DOF, DFF and plant height (PH), four SNPs found associated to PC1 can be inferred as pleiotropic. A window of ± 2 kb of associated SNPs was aligned with available transcriptome data generated for transition from vegetative to reproductive phase in pigeonpea. Annotation analysis of these regions revealed presence of genes which might be involved in floral induction like Cytochrome p450 like Tata box binding protein, Auxin response factors, Pin like genes, F box protein, U box domain protein, chromatin remodelling complex protein, RNA methyltransferase. In summary, it appears that auxin responsive genes could be involved in regulating DOF and DFF as majority of the associated loci contained genes which are component of auxin signaling pathways in their vicinity. Overall, our findings indicates that the use of principal component analysis in GWAS is statistically more robust in terms of identifying genes and FarmCPU is a better choice compared to the other aforementioned models in dealing with both false positive and negative associations and thus can be used for traits with complex inheritance.

RNA-Seq Analysis of Developing Grains of Wheat to Intrigue Into the Complex Molecular Mechanism of the Heat Stress Response

Heat stress is one of the significant constraints affecting wheat production worldwide. To ensure food security for ever-increasing world population, improving wheat for heat stress tolerance is needed in the presently drifting climatic conditions. At the molecular level, heat stress tolerance in wheat is governed by a complex interplay of various heat stress-associated genes. We used a comparative transcriptome sequencing approach to study the effect of heat stress (5°C above ambient threshold temperature of 20°C) during grain filling stages in wheat genotype K7903 (Halna). At 7 DPA (days post-anthesis), heat stress treatment was given at four stages: 0, 24, 48, and 120 h. In total, 115,656 wheat genes were identified, including 309 differentially expressed genes (DEGs) involved in many critical processes, such as signal transduction, starch synthetic pathway, antioxidant pathway, and heat stress-responsive conserved and uncharacterized putative genes that play an essential role in maintaining the grain filling rate at the high temperature. A total of 98,412 Simple Sequences Repeats (SSR) were identified from de novo transcriptome assembly of wheat and validated. The miRNA target prediction from differential expressed genes was performed by psRNATarget server against 119 mature miRNA. Further, 107,107 variants including 80,936 Single nucleotide polymorphism (SNPs) and 26,171 insertion/deletion (Indels) were also identified in de novo transcriptome assembly of wheat and wheat genome Ensembl version 31. The present study enriches our understanding of known heat response mechanisms during the grain filling stage supported by discovery of novel transcripts, microsatellite markers, putative miRNA targets, and genetic variant. This enhances gene functions and regulators, paving the way for improved heat tolerance in wheat varieties, making them more suitable for production in the current climate change scenario.

Genetic Diversity Study on Geographical Populations of the Multipurpose Species Elsholtzia stauntonii Using Transferable Microsatellite Markers

Elsholtzia stauntonii Benth. (Lamiaceae) is an economically important ornamental, medicinal and aromatic plant species. To meet the increasing market demand for E. stauntonii, it is necessary to assess genetic diversity within the species to accelerate the process of genetic improvement. Analysis of the transferability of simple sequence repeat (SSR) markers from related species or genera is a fast and economical method to evaluate diversity, and can ensure the availability of molecular markers in crops with limited genomic resources. In this study, the cross-genera transferability of 497 SSR markers selected from other members of the Lamiaceae (Salvia L., Perilla L., Mentha L., Hyptis Jacq., Leonurus L., Pogostemon Desf., Rosmarinus L., and Scutella L.) to E. stauntonii was 9.05% (45 primers). Among the 45 transferable markers, 10 markers revealed relatively high polymorphism in E. stauntonii. The genetic variation among 825 individuals from 18 natural populations of E. stauntonii in Hebei Province of China was analyzed using the 10 polymorphic SSR markers. On the basis of the SSR data, the average number of alleles (N _A), expected heterozygosity (H _E), and Shannon's information index (I) of the 10 primers pairs were 7.000, 0.478, and 0.688, respectively. Lower gene flow (N _m = 1.252) and high genetic differentiation (F _st = 0.181) were detected in the populations. Analysis of molecular variance (AMOVA) revealed that most of the variation (81.47%) was within the populations. Integrating the results of STRUCTURE, UPGMA (Unweighted Pair Group Method with Arithmetic Mean) clustering, and principal coordinate analysis, the 825 samples were grouped into two clusters associated with geographical provenance (southwestern and northeastern regions), which was consistent with the results of a Mantel test (r = 0.56, p < 0.001). Overall, SSR markers developed in related genera were effective to study the genetic structure and genetic diversity in geographical populations of E. stauntonii. The results provide a theoretical basis for conservation of genetic resources, genetic improvement, and construction of a core collection for E. stauntonii.

Genetic analysis of sucrose concentration in soybean seeds using a historical soybean genomic panel

KEY MESSAGE

Significant QTL for sucrose concentration have been identified using a historical soybean genomic panel, which could aid in the development of food-grade soybean cultivars. Soybean (Glycine max (L.) Merr) is a crop of global importance for both human and animal consumption, which was domesticated in China more than 6000 years ago. A concern about losing genetic diversity as a result of decades of breeding has been expressed by soybean researchers. In order to develop new cultivars, it is critical for breeders to understand the genetic variability present for traits of interest in their program germplasm. Sucrose concentration is becoming an increasingly important trait for the production of soy-food products. The objective of this study was to use a genome-wide association study (GWAS) to identify putative QTL for sucrose concentration in soybean seed. A GWAS panel consisting of 266 historic and current soybean accessions was genotyped with 76 k genotype-by-sequencing (GBS) SNP data and phenotyped in four field locations in Ontario (Canada) from 2015 to 2017. Seven putative QTL were identified on chromosomes 1, 6, 8, 9, 10, 13 and 14. A key gene related to sucrose synthase (Glyma.06g182700) was found to be associated with the QTL located on chromosome 6. This information will facilitate efforts to increase the available genetic variability for sucrose concentration in soybean breeding programs and develop new and improved high-sucrose soybean cultivars suitable for the soy-food industry.

Full-Length SMRT Transcriptome Sequencing and SSR Analysis of Bactrocera dorsalis (Hendel)

Simple Summary

In this study, a full-length transcriptome was analyzed with single-molecule real-time (SMRT) sequencing, which was first used to discover simple sequence repeat (SSR) genetic markers from B. dorsalis. Moreover, SSR markers from isoforms were screened for the identification of species diversity. These results could provide molecular biology methods for further population research.

Abstract

Bactrocera dorsalis (Hendel), as one of the most notorious and destructive invasive agricultural pests in the world, causes damage to over 250 different types of fruits and vegetables throughout tropical and subtropical areas. PacBio single-molecule real-time (SMRT) sequencing was used to generate the full-length transcriptome data of B. dorsalis. A total of 40,319,890 subreads (76.6 Gb, clean reads) were generated, including 535,241 circular consensus sequences (CCSs) and 386,916 full-length non-concatemer reads (FLNCs). Transcript cluster analysis of the FLNC reads revealed 22,780 high-quality reads (HQs). In total, 12,274 transcripts were functionally annotated based on four different databases. A total of 1978 SSR loci were distributed throughout 1714 HQ transcripts, of which 1926 were complete SSRs and 52 were complex SSRs. Among the total SSR loci, 2–3 nucleotide repeats were dominant, occupying 83.62%, of which di- and tri- nucleotide repeats were 39.38% and 44.24%, respectively. We detected 105 repeat motifs, of which AT/AT (50.19%), AC/GT (39.15%), CAA/TTG (32.46%), and ACA/TGT (10.86%) were the most common in di- and tri-nucleotide repeats. The repeat SSR motifs were 12–190 bp in length, and 1638 (88.02%) were shorter than 20 bp. According to the randomly selected microsatellite sequence, 80 pairs of primers were designed, and 174 individuals were randomly amplified by PCR using primers. The number of primers that had amplification products with clear bands and showed good polymorphism came to 41, indicating that this was a feasible way to explore SSR markers from the transcriptomic data of B. dorsalis. These results lay a foundation for developing highly polymorphic microsatellites for researching the functional genomics, population genetic structure, and genetic diversity of B. dorsalis.

Development, characterization, functional annotation and validation of genomic and genic-SSR markers using de novo next generation sequencing in Melia dubia Cav

Melia dubia Cav. (Meliaceae), a fast-growing tropical tree finds use in plywood, pulp and high-value solid wood products. To increase its productivity, we must essentially capture genetic diversity and identify genotypes with superior wood properties. This study aimed to develop novel microsatellite markers from genomic data and validate the markers in M. dubia. Direct Seq-to-SSR approach was adopted and using an in-house Perl script, 426,390 SSR markers identified. For validation, selected 151 markers, of which 50 were genomic markers chosen randomly, and 101 were genic markers identified through BLAST2GO. Amplification was observed in all loci, and 81.4% generated high-quality, reproducible amplicons of the expected size. Out of 50 genomic markers, we used ten highly polymorphic markers to assess genetic diversity among 75 genotypes from three populations. One hundred fourteen alleles were recorded, with a moderate level of diversity and a positive fixation index. Twenty-nine genic markers representing 13 enzymes showing polymorphism for wood stiffness were selected for diversity assessment of 24 genotypes (12 genotypes each with high and low-stress wave velocity). The product size ranged from 87 to 279, covering the majority of the genome. Cluster and structure analysis segregated ~ 80% of the genotypes based on the trait. This is the first report of the development of genic markers from a genomic survey and has proved efficient in differentiating genotypes based on the trait. The markers developed in this study will be useful for genetic mapping, diversity estimation, marker-assisted selection for desired traits and breeding for wood traits in M. dubia.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02858-w.

Full-length SMRT transcriptome sequencing and microsatellite characterization in Paulownia catalpifolia

Paulownia catalpifolia is an important, fast-growing timber species known for its high density, color and texture. However, few transcriptomic and genetic studies have been conducted in P. catalpifolia. In this study, single-molecule real-time sequencing technology was applied to obtain the full-length transcriptome of P. catalpifolia leaves treated with varying degrees of drought stress. The sequencing data were then used to search for microsatellites, or simple sequence repeats (SSRs). A total of 28.83 Gb data were generated, 25,969 high-quality (HQ) transcripts with an average length of 1624 bp were acquired after removing the redundant reads, and 25,602 HQ transcripts (98.59%) were annotated using public databases. Among the HQ transcripts, 16,722 intact coding sequences, 149 long non-coding RNAs and 179 alternative splicing events were predicted, respectively. A total of 7367 SSR loci were distributed throughout 6293 HQ transcripts, of which 763 complex SSRs and 6604 complete SSRs. The SSR appearance frequency was 28.37%, and the average distribution distance was 5.59 kb. Among the 6604 complete SSR loci, 1-3 nucleotide repeats were dominant, occupying 97.85% of the total SSR loci, of which mono-, di- and tri-nucleotide repeats were 44.68%, 33.86% and 19.31%, respectively. We detected 112 repeat motifs, of which A/T (42.64%), AG/CT (12.22%), GA/TC (9.63%), GAA/TTC (1.57%) and CCA/TGG (1.54%) were most common in mono-, di- and tri-nucleotide repeats, respectively. The length of the repeat SSR motifs was 10-88 bp, and 4997 (75.67%) were ≤ 20 bp. This study provides a novel full-length transcriptome reference for P. catalpifolia and will facilitate the identification of germplasm resources and breeding of new drought-resistant P. catalpifolia varieties.

Identification of Transcription Factor Genes and Functional Characterization of PlMYB1 From Pueraria lobata

Kudzu, Pueraria lobata, is a traditional Chinese food and medicinal herb that has been commonly used since ancient times. Kudzu roots are rich sources of isoflavonoids, e.g., puerarin, with beneficial effects on human health. To gain global information on the isoflavonoid biosynthetic regulation network in kudzu, de novo transcriptome sequencings were performed using two genotypes of kudzu with and without puerarin accumulation in roots. RNAseq data showed that the genes of the isoflavonoid biosynthetic pathway were significantly represented in the upregulated genes in the kudzu with puerarin. To discover regulatory genes, 105, 112, and 143 genes encoding MYB, bHLH, and WD40 transcription regulators were identified and classified, respectively. Among them, three MYB, four bHLHs, and one WD40 gene were found to be highly identical to their orthologs involved in flavonoid biosynthesis in other plants. Notably, the expression profiles of PlMYB1, PlHLH3-4, and PlWD40-1 genes were closely correlated with isoflavonoid accumulation profiles in different tissues and cell cultures of kudzu. Over-expression of PlMYB1 in Arabidopsis thaliana significantly increased the accumulation of anthocyanins in leaves and proanthocyanidins in seeds, by activating AtDFR, AtANR, and AtANS genes. Our study provided valuable comparative transcriptome information for further identification of regulatory or structural genes involved in the isoflavonoid pathway in P. lobata, as well as for bioengineering of bioactive isoflavonoid compounds.

Developments and Prospects in Imperative Underexploited Vegetable Legumes Breeding: A Review

Vegetable legumes are an essential source of carbohydrates, vitamins, and minerals, along with health-promoting bioactive chemicals. The demand for the use of either fresh or processed vegetable legumes is continually expanding on account of the growing consumer awareness about their well-balanced diet. Therefore, sustaining optimum yields of vegetable legumes is extremely important. Here we seek to present d etails of prospects of underexploited vegetable legumes for food availability, accessibility, and improved livelihood utilization. So far research attention was mainly focused on pulse legumes' performance as compared to vegetable legumes. Wild and cultivated vegetable legumes vary morphologically across diverse habitats. This could make them less known, underutilized, and underexploited, and make them a promising potential nutritional source in developing nations where malnutrition still exists. Research efforts are required to promote underexploited vegetable legumes, for improving their use to feed the ever-increasing population in the future. In view of all the above points, here we have discussed underexploited vegetable legumes with tremendous potential; namely, vegetable pigeon pea (Cajanus cajan), cluster bean (Cyamopsis tetragonoloba), winged bean (Psophocarpus tetragonolobus), dolichos bean (Lablab purpureus), and cowpea (Vigna unguiculata), thereby covering the progress related to various aspects such as pre-breeding, molecular markers, quantitative trait locus (QTLs), genomics, and genetic engineering. Overall, this review has summarized the information related to advancements in the breeding of vegetable legumes which will ultimately help in ensuring food and nutritional security in developing nations.

Development of EST-SSR markers for Pongamia pinnata by transcriptome database mining: cross-species amplification and genetic diversity

EST-SSR markers were developed from Pongamia pinnata transcriptome libraries. We have successfully utilised EST-SSRs to study the genetic diversity of Indian P. pinnata germplasms and transferability study on legume plants. P. pinnata is a non-edible oil, seed-bearing leguminous tree well known for its multipurpose benefits and acts as a potential source for medicine and biodiesel preparation. Moreover, the plant is not grazable by animal and wildly grown in different agro climatic condition of India. Recently, it is much used in reforestation and rehabilitation of marginal and coal mined land in different part of India. Due to increasing demand for cultivation, understanding of the genetic diversity is important parameter for further breeding and cultivation program. In this investigation, an attempt has been undertaken to develop novel EST-SSR markers by analyzing the assembled transcriptome from previously published Illumina libraries of P. pinnata, which is cross transferrable to legume plants. Twenty EST-SSR markers were developed from oil yielding and secondary metabolite biosynthesis genes. To our knowledge, this is the first EST-SSR marker based genetic diversity study on Indian P. pinnata germplasms. The genetic diversity parameter analysis of P. pinnata showed that the Gangetic plain and Eastern India are highly diverse compared to the Central Deccan and Western germplasms. The lowest genetic diversity in the Western region may be due to the pressure of lower precipitation, high-temperature stress and reduced groundwater availability. Nevertheless, the highest genetic diversity of Gangetic plain and Eastern India may be due to the higher groundwater availability, high precipitation, higher temperature fluctuations and growing by the side of glacier-fed river water. Thus, our study shows the evidence of natural selection on the genetic diversity of P. pinnata germplasms of the Indian subcontinent.

Characterization and Development of Microsatellite Markers in Pseudotaxus chienii (Taxaceae) Based on Transcriptome Sequencing

Pseudotaxus chienii (Taxaceae) is an endangered conifer species endemic to China. However, a lack of suitable molecular markers hinders the genomic and genetic studies on this species. Here, we characterized and developed the microsatellite markers from a newly sequenced P. chienii transcriptome. A total of 21,835 microsatellite loci were detected from 161,131 non-redundant unigene sequences, and the frequency of SSRs was 13.55%, with an average of one SSR loci per 9.18 kb. Mono-nucleotide, di-nucleotide, and tri-nucleotide were the dominant repeat types, accounting for 50.06, 13.49, and 29.39% of the total SSRs, respectively. In terms of distribution location, the coding regions (CDS) with few microsatellites and mainly consisted of tri-nucleotides. There were significant differences in the length of microsatellite among genic regions and motif types. Functional annotation showed that the unigenes containing microsatellites had a wide range of biological functions, most of which were related to basic metabolism, and a few might be involved in expression regulation of gene and response to environmental stress. In addition, 375 primer pairs were randomly selected and synthesized for the amplification and validation of microsatellite markers. Seventy-seven primer pairs were successfully amplified and 40 primer pairs were found to be polymorphic. Finally, 20 pairs of primers with high polymorphism were selected to assess the genetic diversity in four P. chienii populations. In addition, the newly developed microsatellite markers exhibited high transferability (70%) in Amentotaxus argotaenia. Our study could enable further genetic diversity analysis and functional gene mining on Taxaceae.

SSREnricher: a computational approach for large-scale identification of polymorphic microsatellites based on comparative transcriptome analysis

Microsatellite (SSR) markers are the most popular markers for genetic analyses and molecular selective breeding in plants and animals. However, the currently available methods to develop SSRs are relatively time-consuming and expensive. One of the most factors is low frequency of polymorphic SSRs. In this study, we developed a software, SSREnricher, which composes of six core analysis procedures, including SSR mining, sequence clustering, sequence modification, enrichment containing polymorphic SSR sequences, false-positive removal and results output and multiple sequence alignment. After running of transcriptome sequences on this software, a mass of polymorphic SSRs can be identified. The validation experiments showed almost all markers (>90%) that were identified by the SSREnricher as putative polymorphic markers were indeed polymorphic. The frequency of polymorphic SSRs identified by SSREnricher was significantly higher (P < 0.05) than that of traditional and HTS approaches. The software package is publicly accessible on GitHub (https://github.com/byemaxx/SSREnricher).

Genomics-assisted breeding for pigeonpea improvement

The review outlines advances in pigeonpea genomics, breeding and seed delivery systems to achieve yield gains at farmers' field. Pigeonpea is a nutritious and stress-tolerant grain legume crop of tropical and subtropical regions. Decades of breeding efforts in pigeonpea have resulted in development of a number of high-yielding cultivars. Of late, the development of CMS-based hybrid technology has allowed the exploitation of heterosis for yield enhancement in this crop. Despite these positive developments, the actual on-farm yield of pigeonpea is still well below its potential productivity. Growing needs for high and sustainable pigeonpea yields motivate scientists to improve the breeding efficiency to deliver a steady stream of cultivars that will provide yield benefits under both ideal and stressed environments. To achieve this objective in the shortest possible time, it is imperative that various crop breeding activities are integrated with appropriate new genomics technologies. In this context, the last decade has seen a remarkable rise in the generation of important genomic resources such as genome-wide markers, high-throughput genotyping assays, saturated genome maps, marker/gene-trait associations, whole-genome sequence and germplasm resequencing data. In some cases, marker/gene-trait associations are being employed in pigeonpea breeding programs to improve the valuable yield and market-preferred traits. Embracing new breeding tools like genomic selection and speed breeding is likely to improve genetic gains. Breeding high-yielding pigeonpea cultivars with key adaptation traits also calls for a renewed focus on systematic selection and utilization of targeted genetic resources. Of equal importance is to overcome the difficulties being faced by seed industry to take the new cultivars to the doorstep of farmers.

Single-molecule long-read sequencing of the full-length transcriptome of Rhododendron lapponicum L

Rhododendron lapponicum L. is a familiar ornamental plant worldwide with important ornamental and economic value. However, a full-length R. lapponicum transcriptome is still lacking. In the present study, we used the Pacific Biosciences single-molecule real-time sequencing technology to generate the R. lapponicum transcriptome. A total of 346,270 full-length non-chimeric reads were generated, from which we obtained 75,002 high-quality full-length transcripts. We identified 55,255 complete open reading frames, 7,140 alternative splicing events and 2,011 long non-coding RNAs. In gene annotation analyses, 71,155, 33,653, 30,359 and 31,749 transcripts were assigned to the Nr, GO, COG and KEGG databases, respectively. Additionally, 3,150 transcription factors were detected. KEGG pathway analysis showed that 96 transcripts were identified coding for the enzymes associated with anthocyanin synthesis. Furthermore, we identified 64,327 simple sequence repeats from 45,319 sequences, and 150 pairs of primers were randomly selected to develop SSR markers. This study provides a large number of full-length transcripts, which will facilitate the further study of the genetics of R. lapponicum.

A 62K genic-SNP chip array for genetic studies and breeding applications in pigeonpea (Cajanus cajan L. Millsp.)

Pigeonpea is the second most important pulse legume crop for food and nutritional security of South Asia that requires accelerated breeding using high throughput genomic tools. Single nucleotide polymorphisms (SNPs) are highly suitable markers for this purpose because of their bi-allelic nature, reproducibility and high abundance in the genome. Here we report on development and use of a pigeonpea 62 K SNP chip array ‘CcSNPnks’ for Affymetrix GeneTitan^® platform. The array was designed after filtering 645,662 genic-SNPs identified by re-sequencing of 45 diverse genotypes and has 62,053 SNPs from 9629 genes belonging to five different categories, including 4314 single-copy genes unique to pigeonpea, 4328 single-copy genes conserved between soybean and pigeonpea, 156 homologs of agronomically important cloned genes, 746 disease resistance and defense response genes and 85 multi-copy genes of pigeonpea. This fully genic chip has 28.94% exonic, 33.04% intronic, 27.56% 5′UTR and 10.46% 3′UTR SNPs and incorporates multiple SNPs per gene allowing gene haplotype network analysis. It was used successfully for the analysis of genetic diversity and population structure of 95 pigeonpea varieties and high resolution mapping of 11 yield related QTLs for number of branches, pod bearing length and number of seeds per pod in a biparental RIL population. As an accurate high-density genotyping tool, ‘CcSNPnks’ chip array will be useful for high resolution fingerprinting, QTL mapping and genome wide as well as gene-based association studies in pigeonpea.

Genetic Divergence and Relationship Among Opisthopappus Species Identified by Development of EST-SSR Markers

Opisthopappus Shih is an endemic and endangered genus restricted to the Taihang Mountains that has important ornamental and economic value. According to the Flora Reipublicae Popularis Sinicae (FRPS, Chinese version), this genus contains two species (Opisthopappus longilobus and Opisthopappus taihangensis), whereas in the Flora of China (English version) only one species O. taihangensis is present. The interspecific phylogenetic relationship remains unclear and undefined, which might primarily be due to the lack of specific molecular markers for phylogenetic analysis. For this study, 2644 expressed sequence tag-simple sequence repeats (EST-SSRs) from 33,974 unigenes using a de novo transcript assembly of Opisthopappus were identified with a distribution frequency of 7.78% total unigenes. Thereinto, mononucleotides (1200, 45.39%) were the dominant repeat motif, followed by trinucleotides (992, 37.52%), and dinucleotides (410, 15.51%). The most dominant trinucleotide repeat motif was ACC/GGT (207, 20.87%). Based on the identified EST-SSRs, 245 among 1444 designed EST-SSR primers were selected for the development of potential molecular markers. Among these markers, 63 pairs of primers (25.71%) generated clear and reproducible bands with expected sizes. Eventually, 11 primer pairs successfully amplified all individuals from the studied populations. Through the EST-SSR markers, a high level of genetic diversity was detected between Opisthopappus populations. A significant genetic differentiation between the O. longilobus and O. taihangensis populations was found. All studied populations were divided into two clusters by UPGMA, NJ, STRUCTURE, and PCoA. These results fully supported the view of the FRPS, namely, that O. longilobus and O. taihangensis should be regarded as two distinct species. Our study demonstrated that transcriptome sequences, as a valuable tool for the quick and cost-effective development of molecular markers, was helpful toward obtaining comprehensive EST-SSR markers that could contribute to an in-depth assessment of the genetic and phylogenetic relationships between Opisthopappus species.

Characteristics of Microsatellites Mined from Transcriptome Data and the Development of Novel Markers in Paeonia lactiflora

The insufficient number of available simple sequence repeats (SSRs) inhibits genetic research on and molecular breeding of Paeonia lactiflora, a flowering crop with great economic value. The objective of this study was to develop SSRs for P. lactiflora with Illumina RNA sequencing and assess the role of SSRs in gene regulation. The results showed that dinucleotides with AG/CT repeats were the most abundant type of repeat motif in P. lactiflora and were preferentially distributed in untranslated regions. Significant differences in SSR size were observed among motif types and locations. A large number of unigenes containing SSRs participated in catalytic activity, metabolic processes and cellular processes, and 28.16% of all transcription factors and 21.74% of hub genes for inflorescence stem straightness were found to contain SSRs. Successful amplification was achieved with 89.05% of 960 pairs of SSR primers, 55.83% of which were polymorphic, and most of the 46 tested primers had a high level of transferability to the genus Paeonia. Principal component and cluster dendrogram analyses produced results consistent with known genealogical relationships. This study provides a set of SSRs with abundant information for future accession identification, marker-trait association and molecular assisted breeding in P. lactiflora.

J, A. S, Bhat KV, Gaikwad AB. Genome wide identification and characterization of microsatellite markers in black pepper (Piper nigrum): A valuable resource for boosting genomics applications

Black pepper is one of the most valued and widely used spices in the world and dominates multi-billion dollar global spices trade. India is amongst the major producers, consumers and exporters of black pepper. In spite of its commercial and cultural importance, black pepper has received meagre attention in terms of generation of genomic resources. Availability of markers distributed throughout the genome would facilitate and accelerate genetic studies, QTL identification, genetic enhancement and crop improvement in black pepper. In this perspective, the sequence information from the recently sequenced black pepper (Piper nigrum) genome has been used for identification and characterisation of Simple Sequence Repeats (SSRs). Total 69,126 SSRs were identified from assembled genomic sequence of P. nigrum. The SSR frequency was 158 per MB making it, one SSR for every 6.3 kb in the assembled genome. Among the different types of microsatellite repeat motifs, dinucleotides were the most abundant (48.6%), followed by trinucleotide (23.7%) and compound repeats (20.62%). A set of 85 SSRs were used for validation, of which 74 produced amplification products of expected size. Genetic diversity of 30 black pepper accessions using 50 SSRs revealed four distinct clusters. Further, the cross species transferability of the SSRs was checked in nine other Piper species. Out of 50 SSRs used, 19 and 31 SSRs were amplified in nine and seven species, respectively. Thus the identified SSRs may have application in other species of the genus Piper where genome sequence is not available yet. Present study reports the first NGS based genomic SSRs in black pepper and thus constitute a valuable resource for a whole fleet of applications in genetics and plant breeding studies such as genetic map construction, QTL identification, map-based gene cloning, marker-assisted selection and evolutionary studies in Piper nigrum and related species.

A Novel Software and Method for the Efficient Development of Polymorphic SSR Loci Based on Transcriptome Data

Traditional methods for developing polymorphic microsatellite loci without reference sequences are time-consuming and labor-intensive, and the polymorphisms of simple sequence repeat (SSR) loci developed from expressed sequence tag (EST) databases are generally poor. To address this issue, in this study, we developed a new software (PSSRdt) and established an effective method for directly obtaining polymorphism details of SSR loci by analyzing diverse transcriptome data. The new method includes three steps, raw data processing, PSSRdt application, and loci extraction and verification. To test the practicality of the method, we successfully obtained 1940 potential polymorphic SSRs from the transcript dataset combined with 44 pea aphid transcriptomes. Fifty-two SSR loci obtained by the new method were selected for validating the polymorphic characteristics by genotyping in pea aphid individuals. The results showed that over 92% of SSR loci were polymorphic and 73.1% of loci were highly polymorphic. Our new software and method provide an innovative approach to microsatellite development based on RNA-seq data, and open a new path for the rapid mining of numerous loci with polymorphism to add to the body of research on microsatellites.

Comprehensive Stress-Based De Novo Transcriptome Assembly and Annotation of Guar (Cyamopsis tetragonoloba (L.) Taub.): An Important Industrial and Forage Crop

The forage crop Guar (Cyamopsis tetragonoloba (L.) Taub.) has the ability to endure heat, drought, and mild salinity. A complete image on its genic architecture will promote our understanding about gene expression networks and different tolerance mechanisms at the molecular level. Therefore, whole mRNA sequence approach on the Guar plant was conducted to provide a snapshot of the mRNA information in the cell under salinity, heat, and drought stresses to be integrated with previous transcriptomic studies. RNA-Seq technology was employed to perform a 2 × 100 paired-end sequencing using an Illumina HiSeq 2500 platform for the transcriptome of leaves of C. tetragonoloba under normal, heat, drought, and salinity conditions. Trinity was used to achieve a de novo assembly followed by gene annotation, functional classification, metabolic pathway analysis, and identification of SSR markers. A total of 218.2 million paired-end raw reads (~44 Gbp) were generated. Of those, 193.5M paired-end reads of high quality were used to reconstruct a total of 161,058 transcripts (~266 Mbp) with N50 of 2552 bp and 61,508 putative genes. There were 6463 proteins having >90% full-length coverage against the Swiss-Prot database and 94% complete orthologs against Embryophyta. Approximately, 62.87% of transcripts were blasted, 50.46% mapped, and 43.50% annotated. A total of 4715 InterProScan families, 3441 domains, 74 repeats, and 490 sites were detected. Biological processes, molecular functions, and cellular components comprised 64.12%, 25.42%, and 10.4%, respectively. The transcriptome was associated with 985 enzymes and 156 KEGG pathways. A total of 27,066 SSRs were gained with an average frequency of one SSR/9.825 kb in the assembled transcripts. This resulting data will be helpful for the advanced analysis of Guar to multi-stress tolerance.

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.

An "Axiom Cajanus SNP Array" based high density genetic map and QTL mapping for high-selfing flower and seed quality traits in pigeonpea

BACKGROUND

Pigeonpea has considerable extent of insect-aided natural out-crossing that impedes genetic purity of seeds. Pre-anthesis cleistogamy in pigeonpea promotes self-pollination which helps in maintaining genetic purity. The cleistogamous flowers are linked with shriveled seeds, an undesirable trait from variety adoption point of view, and breeding using genomics tools can help in overcoming this constraint. Therefore, in order to identify genomic regions governing these target traits, one recombinant inbred line (RIL) population was developed using contrasting parents (ICPL 99010 and ICP 5529) for flower shape and shriveled seeds. The RILs were phenotyped for two years and genotyped using the Axiom Cajanus SNP Array.

RESULTS

Out of the 56,512 unique sequence variations on the array, the mapping population showed 8634 single nucleotide polymorphism (SNPs) segregating across the genome. These data facilitated generation of a high density genetic map covering 6818 SNPs in 974 cM with an average inter-marker distance of 0.1 cM, which is the lowest amongst all pigeonpea genetic maps reported. Quantitative trait loci (QTL) analysis using this genetic map and phenotyping data identified 5 QTLs associated with cleistogamous flower, 3 QTLs for shriveled seed and 1 QTL for seed size. The phenotypic variance explained by these QTLs ranged from 9.1 to 50.6%. A consistent QTL "qCl3.2" was identified for cleistogamous flower on CcLG03 covering a span of 42 kb in the pigeonpea genome. Epistatic QTLs were also identified for cleistogamous flower and shriveled seed traits.

CONCLUSION

Identified QTLs and genomic interactions for cleistogamous flower, shriveled seed and seed size will help in incorporating the required floral architecture in pigeonpea varieties/lines. Besides, it will also be useful in understanding the molecular mechanisms, and map-based gene cloning for the target traits.

Elymus nutans genes for seed shattering and candidate gene-derived EST-SSR markers for germplasm evaluation

BACKGROUND

Elymus nutans and E. sibiricus are two important forage grasses of the genus Elymus. But they are difficult to grow for commercial seed production due to serious seed shattering. We conducted a comparative transcriptome analysis of abscission zone to find possible transcription changes associated with seed shattering, explore candidate genes involved in seed shattering and identify candidate gene-based EST-SSR markers for germplasm evaluation.

RESULTS

cDNA libraries from abscission zone (AZ) and non-abscission zone (NAZ) tissues of E. nutans were constructed and sequenced. A total of 111,667 unigenes were annotated and 7644 differentially expressed transcripts (DETs) were predicted, corresponding to 6936 up-regulated in AZ and 708 down-regulated in NAZ. We identified 489 candidate genes related to transcription factor, cell wall hydrolysis or modification, hydrolase activity, phytohormone signaling and response, lignin biosynthesis, and signal transduction or protein turnover. Eleven similar candidate genes involved in polygalacturonase activity, hydrolase activity, and mitogen-activated protein kinase were up-regulated in the abscission zone of the two Elymus species, suggesting these genes may have specific function for abscission zone development and seed shattering. A total of 67 polymorphic EST-SSR markers were developed and characterized based on the sequences of these candidate genes. Fourteen polymorphic EST-SSR primers were finally used to study genetic diversity in 48 E. nutans genotypes with contrasting seed shattering habit. The dendrogram based on molecular data showed that most accessions with similar seed shattering degree tended to group together.

CONCLUSIONS

The expression data generated from this study provides an important resource for future molecular biological research. Many DETs were associated with abscission zone development, and EST-SSR loci related to candidate genes may have potential application in identifying trait-associated markers in E. nutans in the future.

Characterization of the Tibet plateau Jerusalem artichoke (Helianthus tuberosus L.) transcriptome by de novo assembly to discover genes associated with fructan synthesis and SSR analysis

BACKGROUND

Jerusalem artichoke (Helianthus tuberosus L.) is a characteristic crop in the Qinghai-Tibet Plateau which has rapidly developed and gained socioeconomic importance in recent years. Fructans are abundant in tubers and represent the foundation for their formation, processing and utilization of yield; and are also widely used in new sugar-based materials, bioenergy processing, ecological management, and functional feed. To identify key genes in the metabolic pathway of fructans in Jerusalem artichoke, high-throughput sequencing was performed using Illumina Hi Seq™ 2500 equipment to construct a transcriptome library.

RESULTS

Qinghai-Tibet Plateau Jerusalem artichoke "Qingyu No.1" was used as the material; roots, stems, leaves, flowers and tubers of Jerusalem artichoke in its flowering stage were mixed into a mosaic of the Jerusalem artichoke transcriptome library, obtaining 63,089 unigenes with an average length of 713.6 bp. Gene annotation through the Nr, Swiss Prot, GO, KOG and KEGG databases revealed 34.95 and 46.91% of these unigenes had similar sequences in the Nr and Swiss Prot databases. The GO classification showed the Jerusalem artichoke unigenes were divided into three ontologies, with a total of 49 functional groups encompassing biological processes, cellular components, and molecular functions. Among them, there were more unigenes involved in the functional groups for cellular processes, metabolic processes, and single-organism processes. 38,999 unigenes were annotated by KOG and divided into 25 categories according to their functions; the most common annotation being general function prediction. A total of 13,878 unigenes (22%) were annotated in the KEGG database, with the largest proportion corresponding to pathways related to carbohydrate metabolism. A total of 12 unigenes were involved in the synthesis and degradation of fructan. Cluster analysis revealed the candidate 12 unigene proteins were dispersed in the 5 major families of proteins involved in fructan synthesis and degradation. The synergistic effect of INV gene is necessary during fructose synthesis and degradation in Jerusalem artichoke tuber development. The sequencing data from the transcriptome of this species can provide a reliable data basis for the identification and assessment of the expression of the members of the INV gene family.A simple sequence repeat (SSR) loci search was performed on the transcriptome data of Jerusalem artichoke, identifying 6635 eligible SSR loci with a large proportion of dinucleotide and trinucleotide repeats, and the most different motifs were repeated 5 times and 6 times. Dinucleotide and trinucleotide repeat motifs were the most frequent, with AG/CT and ACC/GGT repeat motifs accounting for the highest proportion.

CONCLUSIONS

In this study, a database search of the transcriptome of the Jerusalem artichoke from the Qinghai Tibet Plateau was conducted by high throughput sequencing technology to obtain important transcriptional and SSR loci information. This allowed characterization of the overall expression features of the Jerusalem artichoke transcriptome, identifying the key genes involved in metabolism in this species. In turn, this offers a foundation for further research on the regulatory mechanisms of fructan metabolism in Jerusalem artichoke.

Development and Characterization of EST-SSR Markers From RNA-Seq Data in Phyllostachys violascens

Bamboo are woody grass species containing important economic and ecological values. Lei bamboo (Phyllostachys violascens) is a kind of shoot-producing bamboo species with the highest economic yield per unit area. However, identifying different varieties of Lei bamboo based on morphological characteristics is difficult. Microsatellites play an important role in plant identification and genetic diversity analysis and are superior to other molecular markers. In this study, we identified 18,356 expressed sequence tag-simple sequence repeat (EST-SSR) loci in Lei bamboo transcriptome data. A total of 11,264 primer pairs were successfully designed from unigenes of all EST-SSR loci, and 96 primer pairs were randomly selected and synthesized. A total of 54 primer pairs were used for classifying 16 Lei bamboo varieties and 10 different Phyllostachys species. The number of polymorphism alleles among the 54 primer pairs ranged from 3 to 12 for P. violascens varieties and 3 to 20 for Phyllostachys. The phylogenetic tree based on polymorphism alleles successfully distinguished 16 P. violascens varieties and 10 Phyllostachys species. Our study provides abundant EST-SSR resources that are useful for genetic diversity analysis and molecular verification of bamboo and suggests that SSR markers developed from Lei bamboo are more efficient and reliable than ISSR, SRAP or AFLP markers.

Comparative analysis of microsatellite, SNP, and InDel markers in four Rhododendron species based on RNA-seq

Rhododendron possesses valuable horticultural and medicinal properties. However, the genetic studies have been hindered due to the lack of genetic markers. Based on RNA-seq, large-scale molecular markers were developed from four Rhododendron species endemic to Dabie Mountains (central China): R. fortunei (5.25 Gb; SSRs, 12,756, one/2.37 kb, 147 types; SNPs, 38,313; InDels, 3,174), R. simsii (5.80 Gb; SSRs, 13,294, one/2.58 kb, 167 types; SNPs, 136,590; InDels, 6,258), R. mariesii (6.53 Gb; SSRs, 15,724, one/2.51 kb, 170 types; SNPs, 44,942; InDels, 4,126), and R. molle (4.35 Gb; SSRs, 10,214, one/2.49 kb, 110 types; SNPs, 77,829; InDels, 3,416). Di-nucleotide repeats were the main type (59.126%-64.314%), and AG/CT repeat (55.18%-61.22%) was the most. In particular, 89 species-specific types had been found. Furthermore, C:G→T:A mutation was the main SNP type (30.475%-34.99%). However, C:G→G:C mutation was the least type in R. fortunei, while T:A→G:C mutation was the least in the other three species. InDels with length of 3 nt was most in R. fortunei, but 1 nt InDels were the main type in the other three species. Twelve microsatellite markers developed from R. simsii reveled high genetic diversity in the four populations, and heterozygote excess was observed. This research would benefit the genetic study, molecular marker-assisted selection, and breeding studies in Rhododendron species.

Development of novel EST-SSR markers in the macaúba palm (Acrocomia aculeata) using transcriptome sequencing and cross-species transferability in Arecaceae species

BACKGROUND

The macaúba palm is a novel feedstock for oil production suitable for multiple uses, including as biodiesel and in the food and cosmetic industries. As an efficient alternative, the macaúba palm has limited genomic resources, particularly expressed sequence tag (EST) markers. We report a comprehensive set of validated EST-simple sequence repeat (SSR) markers by using transcriptome sequencing, its application in genetic diversity analysis and cross transferability in other palm trees with environmental and economic importance.

RESULTS

In this study, a total of 418 EST-SSRs were identified to be unique for one transcript and region; 232 EST-SSRs were selected, with trinucleotide repeats being the most frequent motif, representing 380 (90.9%), followed by composited (4.5%), di- (3.6%), and hexanucleotides (3.6%). A total of 145 EST-SSRs (62.5%) were validated for consistent amplification in seventeen macaúba palm samples, and 100 were determined to be polymorphic with PIC values ranging from 0.25 to 0.77. Genetic diversity analysis was performed with the 20 most informative EST-SSR markers showing a distinct separation of the different groups of macaúba palm. Additionally, these 145 markers were transferred in six other palm species resulting in transferability rates of 99% (144) in Acrocomia intumescens, 98% (143) in Acrocomia totai, 80.7% (117 EST-EST) in African oil palm (Elaeis guineensis) and peach palm (Bactris gasipaes) samples, 70% (102) in the juçara palm (Euterpe edulis) and 71.7% (104) in the hat palm (Sabal causiarum). Analysis of genetic distance showed a high separation in accordance with geographic location, establishing distinct groups by genera.

CONCLUSIONS

The EST markers identified in our study are a valuable resource and provide a genomic tool for genetic mapping and further genetic studies, as well as evaluation of co-location between QTLs and functionally associated markers.

Augmentation of crop productivity through interventions of omics technologies in India: challenges and opportunities

With the continuous increase in the population of developing countries and decline of natural resources, there is an urgent need to qualitatively and quantitatively augment crop productivity by using new tools and technologies for improvement of agriculturally important traits. The new scientific and technological omics-based approaches have enabled us to deal with several issues and challenges faced by modern agricultural system and provided us novel opportunities for ensuring food and nutritional security. Recent developments in sequencing techniques have made available huge amount of genomic and transcriptomic data on model and cultivated crop plants including Arabidopsis thaliana, Oryza sativa, Triticum aestivum etc. The sequencing data along with other data generated through several omics platforms have significantly influenced the disciplines of crop sciences. Gene discovery and expression profiling-based technologies are offering enormous opportunities to the scientific community which can now apply marker-assisted selection technology to assess and enhance diversity in their collected germplasm, introgress essential traits from new sources and investigate genes that control key traits of crop plants. Utilization of omics science and technologies for crop productivity, protection and management has recently been receiving a lot of attention; the majority of the efforts have been put into signifying the possible applications of various omics technologies in crop plant sciences. This article highlights the background of challenges and opportunities for augmentation of crop productivity through interventions of omics technologies in India.

Identification, analysis and development of salt responsive candidate gene based SSR markers in wheat

BACKGROUND

Salinity severely limits wheat production in many parts of the world. Development of salt tolerant varieties represents the most practical option for enhancing wheat production from these areas. Application of marker assisted selection may assist in fast tracking development of salt tolerant wheat varieties. However, SSR markers available in the public domain are not specifically targeted to functional regions of wheat genome, therefore large numbers of these need to be analysed for identification of markers associated with traits of interest. With the availability of a fully annotated wheat genome assembly, it is possible to develop SSR markers specifically targeted to genic regions. We performed extensive analysis to identify candidate gene based SSRs and assessed their utility in characterizing molecular diversity in a panel of wheat genotypes.

RESULTS

Our analysis revealed, 161 SSR motifs in 94 salt tolerance candidate genes of wheat. These SSR motifs were nearly equally distributed on the three wheat sub-genomes; 29.8% in A, 35.7% in B and 34.4% in D sub-genome. The maximum number of SSR motifs was present in exons (31.1%) followed by promoters (29.8%), 5'UTRs (21.1%), introns (14.3%) and 3'UTRs (3.7%). Out of the 65 candidate gene based SSR markers selected for validation, 30 were found polymorphic based on initial screening and employed for characterizing genetic diversity in a panel of wheat genotypes including salt tolerant and susceptible lines. These markers generated an average of 2.83 alleles/locus. Phylogenetic analysis revealed four clusters. Salt susceptible genotypes were mainly represented in clusters I and III, whereas high and moderate salt tolerant genotypes were distributed in the remaining two clusters. Population structure analysis revealed two sub-populations, sub-population 1 contained the majority of salt tolerant whereas sub-population 2 contained majority of susceptible genotypes. Moreover, we observed reasonably higher transferability of SSR markers to related wheat species.

CONCLUSION

We have developed salt responsive gene based SSRs in wheat for the first time. These were highly useful in unravelling functional diversity among wheat genotypes with varying responses to salt stress. The identified gene based SSR markers will be valuable genomic resources for genetic/association mapping of salinity tolerance in wheat.

Analysis of evolution and genetic diversity of sweetpotato and its related different polyploidy wild species I. trifida using RAD-seq

BACKGROUND

Sweetpotato (Ipomoea batatas (L.) Lam.) is one of the most important crops from the family of Convolvulaceae. It is widely reported that cultivated sweetpotato was originated from Ipomoea trifida. However, diploid, tetraploid and hexaploid I. trifida were found in nature. The relationship, between them, and among them and sweetpotato, is remaining unclear.

RESULTS

In the present study, we detected the genome diversity and relationship of sweetpotato and different polyploidy types I. trifida using Restriction-site Associated DNA Sequencing (RAD-seq). A total of 38,605 RAD-tags containing 832,204 SNPs had been identified. These tags were annotated using five public databases, about 11,519 tags were aligned to functional genes in various pathways. Based on SNP genotype, phylogenetic relation analysis results confirmed that cultivated sweetpotato has a closer relationship with I. trifida 6× than with I. trifida 4X and I. trifida 2×. Besides, 5042 SSRs were detected in I. trifida 6×, and 3202 pairs of high-quality SSR primers were developed. A total of 68 primers were randomly selected and synthesized, of which 61 were successfully amplified.

CONCLUSION

These results provided new evidence that cultivated sweetpotato originated from I. trifida 6×, and that I. trifida 6× evolved from I. trifida 4X and I. trifida 2×. Therefore, using I. trifida 6× as the model plant of sweetpotato research should be more practical than using I. trifida 2× in the future. Meanwhile, sequence information and markers from the present study will be helpful for sweetpotato and I. trifida studies in the future.

Comparative transcriptome analysis reveals the gene expression profiling in bighead carp (Aristichthys nobilis) in response to acute nitrite toxicity

OBJECTIVE

Nitrite exposure induces growth inhibition, metabolic disturbance, oxidative stress, organic damage, and infection-mediated mortality of aquatic organism. This study aimed to investigate the mechanism in responses to acute nitrite toxicity in bighead carp (Aristichthys nobilis, A. nobilis) by RNA-seq analysis.

METHODS

Bighead carps were exposed to water with high nitrite content (48.63 mg/L) for 72 h, and fish livers and gills were separated for RNA-seq analysis. De novo assembly was performed, and differentially expressed genes (DEGs) between control and nitrite-exposed fishes were identified. Furthermore, enrichment analysis was performed for DEGs to annotate the molecular functions.

RESULTS

A total of 406,135 transcripts and 352,730 unigenes were tagged after de novo assembly. Accordingly, 4108 and 928 DEGs were respectively identified in gill and liver in responses to nitrite exposure. Most of these DEGs were up-regulated DEGs. Enrichment analysis showed these DEGs were mainly associated with immune responses and nitrogen metabolism.

CONCLUSIONS

We suggested that the nitrite toxicity-induced DEGs were probably related to dysregulation of nitrogen metabolism and immune responses in A. nobilis, particularly in gill.

Revealing the complex genetic structure of cultivated amaryllis (Hippeastrum hybridum) using transcriptome-derived microsatellite markers

Although amaryllis (Hippeastrum hybridum) plants are commonly used in physiological and ecological research, the extent of their genomic and genetic resources remains limited. The development of molecular markers is therefore of great importance to accelerate genetic improvements in Hippeastrum species. In this study, a total of 269 unique genes were defined that might regulate the flower spathe development of amaryllis. In addition, 2000 simple sequence repeats (SSRs) were detected based on 171,462 de novo assembled unigenes from transcriptome data, and 66,4091 single nucleotide polymorphisms (SNPs) were also detected as putative molecular markers. Twenty-one SSR markers were screened to evaluate the genetic diversity and population structure of 104 amaryllis accessions. A total of 98 SSR loci were amplified for all accessions. The results reveal that Nei’s gene diversity (H) values of these markers ranged between 0.055 and 0.394, whereas the average values of Shannon’s Information index (I) ranged between 0.172 and 0.567. Genetic tree analysis further demonstrates that all accessions can be grouped into three main clusters, which can be further divided into two subgroups. STRUCTURE-based analysis revealed that the highest ΔK values were observed when K = 5, K = 6, K = 7 and K = 8. The results of this study enable large-scale transcriptomics and classification of Hippeastrum genetic polymorphisms and will be useful in the future for resource conservation and production.

Genome-wide characterization of simple sequence repeats in Pyrus bretschneideri and their application in an analysis of genetic diversity in pear

Background

Pear (Pyrus spp.) is an economically important temperate fruit tree worldwide. In the past decade, significant progress has been made in pear molecular genetics based on DNA research, but the number of molecular markers is still quite limited, which hardly satisfies the increasing needs of geneticists and breeders.

Results

In this study, a total of 156,396 simple sequence repeat (SSR) loci were identified from a genome sequence of Pyrus bretschneideri ‘Dangshansuli’. A total of 101,694 pairs of SSR primers were designed from the SSR loci, and 80,415 of the SSR loci were successfully located on 17 linkage groups (LGs). A total of 534 primer pairs were synthesized and preliminarily screened in four pear cultivars, and of these, 332 primer pairs were selected as clear, stable, and polymorphic SSR markers. Eighteen polymorphic SSR markers were randomly selected from the 332 polymorphic SSR markers in order to perform a further analysis of the genetic diversity among 44 pear cultivars. The 14 European pears and their hybrid materials were clustered into one group (European pear group); 29 Asian pear cultivars were clustered into one group (Asian pear group); and the Zangli pear cultivar ‘Deqinli’ from Yunnan Province, China, was grouped in an independent group, which suggested that the cultivar ‘Deqinli’ is a distinct and valuable germplasm resource. The population structure analysis partitioned the 44 cultivars into two populations, Pop 1 and Pop 2. Pop 2 was further divided into two subpopulations. Results from the population structure analysis were generally consistent with the results from the UPGMA cluster analysis.

Conclusions

The results of the present study showed that the use of next-generating sequencing to develop SSR markers is fast and effective, and the developed SSR markers can be utilized by researchers and breeders for future pear improvement.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4822-7) contains supplementary material, which is available to authorized users.

SSR marker development and intraspecific genetic divergence exploration of Chrysanthemum indicum based on transcriptome analysis

BACKGROUND

Chrysanthemum indicum L., an important ancestral species of the flowering plant chrysanthemum, can be used as medicine and for functional food development. Due to the lack of hereditary information for this species and the difficulty of germplasm identification, we herein provide new genetic insight from the perspective of intraspecific transcriptome comparison and present single sequence repeat (SSR) molecular marker recognition technology.

RESULTS

Through the study of a diploid germplasm (DIWNT) and a tetraploid germplasm (DIWT), the following outcome were obtained. (1) A significant difference in Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations for specific homologous genes was observed using the OrthoMCL method for the identification of homologous gene families between the two cytotypes. Ka/Ks analysis of common, single-copy homologous family members also revealed a greater difference among genes that experienced positive selection than among those experiencing positive selection. (2) Of more practical value, 2575 SSR markers were predicted and partly verified. We used TaxonGap as a visual tool to inspect genotype uniqueness and screen for high-performance molecular loci; we recommend four primers of 65 randomly selected primers with a combined identification success rate of 88.6% as priorities for further development of DNA fingerprinting of C. indicum germplasm.

CONCLUSIONS

The SSR technology based on next-generation sequencing was proved to be successful in the identification of C. indicum germplasms. And the information on the intraspecfic genetic divergence generated by transcriptome comparison deepened the understanding of this complex species' nature.

Isolation and characterization of microsatellite markers in Garcinia gummi-gutta by next-generation sequencing and cross-species amplification

Garcinia gummi-gutta (L.) Roxb. (Clusiaceae) is an endemic, semidomesticated, fruit-yielding tree species distributed in the Western Ghats of India and Sri Lanka. Various bioactive phytochemicals, such as garcinol, benzophenones and xanthones are isolated from G. gummi-gutta and have shown antibacterial, antiviral and antioxidant activities. We sequenced the total genomic DNA using Illumina Hiseq 2000 platform and examined 241,141,804 bp high quality data, assembled into 773,889 contigs. In these contigs, 27,313 simple-sequence repeats (SSRs) were identified, among which mononucleotide repeats were predominant (44.98%) followed by dinucleotide and trinucleotide repeats. Primers were designed for 9964 microsatellites among which 32 randomly selected SSR primer pairs were standardized for amplification. Polymerase chain reaction (PCR) amplification of genomic DNA in 30 G. gummi-gutta genotypes revealed polymorphic information content (PIC) across all 32 loci ranging from 0.867 to 0.951, with a mean value of 0.917. The observed and expected heterozygosity ranged from 0.00 to 0.63 and 0.896 to 0.974, respectively. Alleles per locus ranged from 12 to 27. This is the first report on the development of genomic SSR markers in G. gummi-gutta using next-generation sequencing technology. The genomic SSR markers developed in this study will be useful in identification, mapping, diversity and breeding studies.

Genetic Variation and Association Analysis of the SSR Markers Linked to the Major Drought-Yield QTLs of Rice

Drought is one of the major abiotic stresses, which hampers the production of rice worldwide. Informative molecular markers are valuable tools for improving the drought tolerance in various varieties of rice. The present study was conducted to evaluate the informative simple sequence repeat (SSR) markers in a diverse set of rice genotypes. The genetic diversity analyses of the 83 studied rice genotypes were performed using 34 SSR markers closely linked to the major quantitative trait loci (QTLs) of grain yield under drought stress (qDTYs). In general, our results indicated high levels of polymorphism. In addition, we screened these rice genotypes at the reproductive stage under both drought stress and nonstressful conditions. The results of the regression analysis demonstrated a significant relationship between 11 SSR marker alleles and the plant paddy weight under stressful conditions. Under the nonstressful conditions, 16 SSR marker alleles showed a significant correlation with the plant paddy weight. Finally, four markers (RM279, RM231, RM166, and RM231) demonstrated a significant association with the plant paddy weight under both stressful and nonstressful conditions. These informative-associated alleles may be useful for improving the crop yield under both drought stress and nonstressful conditions in breeding programs.