Characterization of the Lycium barbarum fruit transcriptome and development of EST-SSR markers

Combined Transcriptomics and Metabolomics Analysis Reveals the Effect of Selenium Fertilization on Lycium barbarum Fruit

As a beneficial nutrient and essential trace element, selenium plays a significant role in plant growth functions and human protein biosynthesis. Plant selenium enrichment is mainly obtained from both natural soil and exogenous selenium supplementation, while human beings consume selenium-enriched foods for the purposes of selenium supplementation. In this study, different types of selenium fertilizers were sprayed onto Lycium barbarum in Ningxia, and transcriptomics and metabolomics techniques were used to explore the effects of selenium on the fruit differentials and differential genes in Lycium barbarum. Taking the "Ning Qiyi No.1" wolfberry as the research object, sodium selenite, nano-selenium, and organic selenium were sprayed at a concentration of 100 mg·L-1 three times from the first fruiting period to the harvesting period, with a control treatment comprising the spraying of clear water. We determined the major metabolites and differential genes of the amino acids and derivatives, flavonoids, and alkaloids in ripe wolfberries. We found that spraying selenium significantly enhanced the Lycium barbarum metabolic differentiators; the most effective spray was the organic selenium, with 129 major metabolic differentiators and 10 common metabolic pathways screened after spraying. Nano-selenium was the next best fertilizer we screened, with 111 major metabolic differentiators, the same number as organic selenium in terms of differential genes and common metabolite pathways. Sodium selenite was the least effective of the three, with only 59 of its major metabolic differentials screened, but its differential genes and metabolites were enriched for five common pathways.

Analysis of the genetic diversity and population structure of Monochasma savatieri Franch. ex Maxim using novel EST-SSR markers

Background

Monochasma savatieri Franch. ex Maxim is a medicinally valuable herb. However, the collection and protection of the wild germplasm resources of M. savatieri are still insufficient, and their genetic diversity and population structure have been poorly studied.

Results

We collected and examined 46 M. savatieri individuals from Fujian, Hunan, Jiangxi, and Zhejiang provinces for genetic diversity and population structure, using 33 newly developed expressed sequence tag-simple sequence repeat (EST-SSR) markers. Applying these markers, we detected a total of 208 alleles, with an average of 6.303 alleles per locus. The polymorphic information content varied from 0.138 to 0.884 (average: 0.668), indicating a high level of polymorphism. At the population level, there was a low degree of genetic diversity among populations (I = 0.535, He = 0.342), with Zhejiang individuals showing the highest genetic diversity among the four populations (Fst = 0.497), which indicated little gene flow within the M. savatieri populations (Nm = 0.253). Mantel test analysis revealed a significant positive correlation between geographical and genetic distance among populations (R² = 0.3304, p < 0.05), and structure and principal coordinate analyses supported classification of populations into three clusters, which was consistent with the findings of cluster analysis.

Conclusions

As a rare medicinal plants, the protection of M. savatieri does not look optimistic, and accordingly, protective efforts should be beefed up on the natural wild populations. This study provided novel tools and insights for designing effective collection and conservation strategies for M. savatieri.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08832-x.

Basic Characteristics of Flower Transcriptome Data and Derived Novel EST-SSR Markers of Luculia yunnanensis, an Endangered Species Endemic to Yunnan, Southwestern China

Luculia yunnanensis (Rubiaceae), an evergreen shrub or small tree, is endemic to China and confined to Nujiang Prefecture, Yunnan Province. This plant is of high ornamental value owing to its attractive pink flowers, sweet fragrance, and long flowering period. Due to the influence of climate change and human factors, the distribution range of L. yunnanensis has exhibited a significant shrinking trend, and it has become a vulnerable species that is in urgent need of conservation and rational utilization research. In this study, the flower transcriptome sequencing of L. yunnanensis was conducted using an Illumina HiSeq platform. We designed and developed a series of EST-SSR primers based on the flower transcriptome data of L. yunnanensis. The results showed that 98,389 unigenes were obtained from the L. yunnanensis flower transcriptome, all of which were aligned with sequences in public databases. Nr, Nt, Pfam, KOG/COG, Swiss-Prot, KEGG, and GO annotated 31,859, 13,853, 22,684, 10,947, 21,416, 9722, and 23,390 unigenes, respectively. The MISA (Microsatellite) tool was used to identify SSR loci from all unigenes, and a total of 15,384 SSRs were identified. Repeat motifs were given priority with mononucleotides, dinucleotides, and trinucleotides. The 81 primer pairs were synthesized randomly, of which 44 pairs showed effective amplification. A total of 17 primers showed stable amplification, and rich polymorphism was observed in 6 populations. We concluded via genetic diversity analysis that the average effective number of alleles (Ne), Shannon’s information index (I), and polymorphism information content (PIC) were 1.925, 0.837, and 0.403, respectively. In conclusion, 17 EST-SSR primers can be used for subsequent population genetic diversity analysis and molecular-marker-assisted breeding, which is of great significance for formulating resource conservation and utilization strategies for L. yunnanensis.

Discrimination of Geographical Origins of Wolfberry (Lycium barbarum L.) Fruits Using Stable Isotopes, Earth Elements, Free Amino Acids, and Saccharides

To develop sophisticated approaches for distinguishing goji origins, 325 wolfberry fruit samples of a certain cultivar, plant age, drying method, and collection season were gathered from 26 producing areas across Northwest China in 2017 and 2018. We employed 49 indices, including stable isotopes, earth elements, soluble amino acids, and saccharides, to identify the regions of origin of these goji fruits. Analysis of variance (ANOVA) and heritability analysis were used to assess the effects of the environment (producing areas), cultivar, plant age, drying process, and collection season. Samples from the same place can be classified and partially discriminated using principal component analysis (PCA). We were able to distinguish fruits produced in Zhongning County from those produced in the other five producing provinces using orthogonal projection to latent structure-discriminant analysis (OPLS-DA). Calcium (Ca), manganese (Mn), ornithine (Orn), cystine (Cys-Cys), glutamate (Glu), phenylalanine (Phe), phosphoserine (Ps), serine (Ser), lysine (Lys), taurine (Tau), proline (Pro), and tyrosine (Tyr) indices were chosen using S-plots and heritability analysis, and their repeatability was established with samples collected in 2018. The indices selected in this study can distinguish goji berries produced in Zhongning County from fruits originating from five other Provinces with high repeatability, which was validated with various cultivars, drying methods, harvest seasons, and plant ages and with heritability analysis.

Germplasm Resources and Strategy for Genetic Breeding of Lycium Species: A Review

Lycium species (goji), belonging to Solanaceae, are widely spread in the arid to semiarid environments of Eurasia, Africa, North and South America, among which most species have affinal drug and diet functions, resulting in their potential to be a superior healthy food. However, compared with other crop species, scientific research on breeding Lycium species lags behind. This review systematically introduces the present germplasm resources, cytological examination and molecular-assisted breeding progress in Lycium species. Introduction of the distribution of Lycium species around the world could facilitate germplasm collection for breeding. Karyotypes of different species could provide a feasibility analysis of fertility between species. The introduction of mapping technology has discussed strategies for quantitative trait locus (QTL) mapping in Lycium species according to different kinds of traits. Moreover, to extend the number of traits and standardize the protocols of trait detection, we also provide 1,145 potential traits (275 agronomic and 870 metabolic) in different organs based on different reference studies on Lycium, tomato and other Solanaceae species. Finally, perspectives on goji breeding research are discussed and concluded. This review will provide breeders with new insights into breeding Lycium species.

Development and Application of EST-SSR Markers in Cephalotaxus oliveri From Transcriptome Sequences

Cephalotaxus oliveri is an endemic conifer of China, which has medicinal and ornamental value. However, the limited molecular markers and genetic information are insufficient for further genetic studies of this species. In this study, we characterized and developed the EST-SSRs from transcriptome sequences for the first time. The results showed that a total of 5089 SSRs were identified from 36446 unigenes with a density of one SSR per 11.1 kb. The most common type was trinucleotide repeats, excluding mononucleotide repeats, followed by dinucleotide repeats. AAG/CTT and AT/AT exhibited the highest frequency in the trinucleotide and dinucleotide repeats, respectively. Of the identified SSRs, 671, 1125, and 1958 SSRs were located in CDS, 3′UTR, and 5′UTR, respectively. Functional annotation showed that the SSR-containing unigenes were involved in growth and development with various biological functions. Among successfully designed primer pairs, 238 primer pairs were randomly selected for amplification and validation of EST-SSR markers and 47 primer pairs were identified as polymorphic. Finally, 28 high-polymorphic primers were used for genetic analysis and revealed a moderate level of genetic diversity. Seven natural C. oliveri sampling sites were divided into two genetic groups. Furthermore, the 28 EST-SSRs had 96.43, 71.43, and 78.57% of transferability rate in Cephalotaxus fortune, Ametotaxus argotaenia, and Pseudotaxus chienii, respectively. These markers developed in this study lay the foundation for further genetic and adaptive evolution studies in C. oliveri and related species.

Wolfberry genomes and the evolution of Lycium (Solanaceae)

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

Yeast Synthetic Biology for the Production of Lycium barbarum Polysaccharides

The fruit of Lycium barbarum L. (goji berry) is used as traditional Chinese medicine, and has the functions of immune regulation, anti-tumor, neuroprotection, anti-diabetes, and anti-fatigue. One of the main bioactive components is L. barbarum polysaccharide (LBP). Nowadays, LBP is widely used in the health market, and it is extracted from the fruit of L. barbarum. The planting of L. barbarum needs large amounts of fields, and it takes one year to harvest the goji berry. The efficiency of natural LBP production is low, and the LBP quality is not the same at different places. Goji berry-derived LBP cannot satisfy the growing market demands. Engineered Saccharomyces cerevisiae has been used for the biosynthesis of some plant natural products. Recovery of LBP biosynthetic pathway in L. barbarum and expression of them in engineered S. cerevisiae might lead to the yeast LBP production. However, information on LBP biosynthetic pathways and the related key enzymes of L. barbarum is still limited. In this review, we summarized current studies about LBP biosynthetic pathway and proposed the strategies to recover key enzymes for LBP biosynthesis. Moreover, the potential application of synthetic biology strategies to produce LBP using engineered S. cerevisiae was discussed.

LbNR-Derived Nitric Oxide Delays Lycium Fruit Coloration by Transcriptionally Modifying Flavonoid Biosynthetic Pathway

Anthocyanin-derived fleshy fruit pigmentation has become an excellent system for studying the regulatory network underlying fruit ripening and quality. The transcriptional control of anthocyanin biosynthesis by MYB-bHLH-WDR complexes has been well established, but the intermediate signals through which the environmental or developmental cues regulate these transcription factors remain poorly understood. Here we found that nitric oxide (NO) production during Lycium fruit ripening decreased progressively presenting a negative relationship with anthocyanins. After cloning of the nitric reductase (NR) gene from Lycium barbarum (LbNR) plants, we demonstrated that LbNR-derived NO partially inhibited anthocyanin biosynthesis but enhanced proanthocyanidin (PA) accumulation, and delayed fruit coloration. Application of the NO donor, sodium nitroprusside (SNP), produced a similar effect. The endogenous or exogenous NO downregulated the transcripts both of the regulatory genes and the structural genes that related to anthocyanin biosynthesis, while upregulated both of those genes that related to PA biosynthesis. Given there is a significant negative relationship between the levels of anthocyanins and PAs during Lycium fruit ripening, NO not only inhibited anthocyanin de novo biosynthesis but redirected the flavonoid biosynthetic pathway from anthocyanins to PA production. Two types of LrMYB transcription factors of opposite nature, namely anthocyanin-specific and PA-specific, which belong to the R2R3-MYB subfamily and 1R-MYB subfamily, respectively, were identified from L. ruthenicum fruits. It was further found that NO acts by antagonizing the ABA signaling, a phytohormone we have previously shown playing a positive role in Lycium fruit coloration. Our results provided particularly novel information about NO-ABA-anthocyanin interplay during Lycium fruit development and ripening, which may fill a gap between the developmental cues and the transcriptional regulation of anthocyanin biosynthesis.

Fruit ripening in Lycium barbarum and Lycium ruthenicum is associated with distinct gene expression patterns

Goji berries have been used as food and medicine for millennia. Due to their high morphological similarity, fruits of two distinct species belonging to the family Solanaceae, Lycium barbarum (LB) and Lycium chinense (Chinese boxthorn), are usually marketed together as goji berries, but nearly 90% of all commercially available goji berries belong to the former species. A third closely related species, a wild perennial thorny shrub native to north‐western China, Lycium ruthenicum (LR; known as Russian box thorn, and its fruit as black wolfberry), has become a popular choice for combating soil desertification and for alleviating soil salinity/alkalinity due to its high resistance to the harsh environment of saline deserts. Despite the phylogenetic closeness of LB and LR, their fruits are very different. To identify the genes involved in these distinct phenotypes, here we studied expression patterns of 22 transcriptional regulators that may be crucial drivers of these differences during five developmental stages. BAM1 may contribute to higher sugar content in LB. High expression of BFRUCT in ripe LR is likely to be an evolutionary adaptation to fruit ripening in an arid environment. Two arogenate dehydratase paralogues, CHS and LDOX, are probably crucial elements of the mechanism by which LR accumulates much higher levels of anthocyanin. DXS2 (carotenoid accumulation in LB) and CCD4 (carotenoid degradation in ripe LR fruit) may be crucial drivers behind the much higher content of carotenoids in LB. EIL3 and ERF5 are two transcription factors that may contribute to the higher abiotic stress resilience of LR. GATA22‐like appears to have more important roles in growth than ripening in LB fruit and vice versa in LR. HAT5‐like exhibited opposite temporal patterns in two fruits: high in the 1st stage in LB and high in the 5th stage in LR. PED1 was expressed at a much lower level in LR. Finally, we hypothesise that the poorly functionally characterised SCL32 gene may play a part in the increased resistance to environmental stress of LR. We suggest that BAM1, BFRUCT, EIL3, ERF5, ADT paralogues (for functional redundancy), PED1, GATA22‐like, HAT5‐like and SCL32 warrant further functional studies.

The effects of Lycium barbarum water extract and fish collagen on milk proteolysis and in vitro angiotensin I-converting enzyme inhibitory activity of yogurt

Plain and Lycium barbarum yogurt were made in the presence and absence of fish collagen. Yogurt samples were analyzed for acidification, milk protein proteolysis, angiotensin I-converting enzyme (ACE) inhibitory activity, and sensory evaluation during refrigerated storage for up to 21 days. The o-phthaldialdehyde peptides amount of L. barbarum yogurt both in the presence and absence of fish collagen were significantly increased during 14 days of storage. SDS-PAGE showed improvement in whey proteins degradation of L. barbarum yogurt with/without fish collagen after 3 weeks of storage. L. barbarum yogurt in absence of fish collagen was acting as a great ACE inhibitor reached up to 85% on day 7 of storage. The incorporation of L. barbarum and/or fish collagen affected to a small extent the overall sensory characteristics of yogurt. Yogurt supplemented with L. barbarum and/or fish collagen may lead to the improvement in the production and formulation of yogurt differing in their anti-ACE activity.

Comparative transcriptome analysis of two contrasting wolfberry genotypes during fruit development and ripening and characterization of the LrMYB1 transcription factor that regulates flavonoid biosynthesis

Background

Lycium barbarum and L. ruthenicum have been used as traditional medicinal plants in China and other Asian counties for centuries. However, the molecular mechanisms underlying fruit development and ripening, as well as the associated production of medicinal and nutritional components, have been little explored in these two species.

Results

A competitive transcriptome analysis was performed to identify the regulators and pathways involved in the fruit ripening of red wolfberry (L. barbarum) and black wolfberry (L. ruthenicum) using an Illumina sequencing platform. In total, 155,606 genes and 194,385 genes were detected in red wolfberry (RF) and black wolfberry (BF), respectively. Of them, 20,335, 24,469, and 21,056 genes were differentially expressed at three different developmental stages in BF and RF. Functional categorization of the differentially expressed genes revealed that phenylpropanoid biosynthesis, flavonoid biosynthesis, anthocyanin biosynthesis, and sugar metabolism were the most differentially regulated processes during fruit development and ripening in the RF and BF. Furthermore, we also identified 38 MYB transcription factor-encoding genes that were differentially expressed during black wolfberry fruit development. Overexpression of LrMYB1 resulted in the activation of structural genes for flavonoid biosynthesis and led to an increase in flavonoid content, suggesting that the candidate genes identified in this RNA-seq analysis are credible and might offer important utility.

Conclusion

This study provides novel insights into the molecular mechanism of Lycium fruit development and ripening and will be of value to novel gene discovery and functional genomic studies.

Transcriptomic and metabolomic analyses of Lycium ruthenicum and Lycium barbarum fruits during ripening

Red wolfberry (or goji berry, Lycium barbarum; LB) is an important agricultural product with a high content of pharmacologically important secondary metabolites such as phenylpropanoids. A close relative, black wolfberry (L. ruthenicum; LR), endemic to the salinized deserts of northwestern China, is used only locally. The two fruits exhibit many morphological and phytochemical differences, but genetic mechanisms underlying them remain poorly explored. In order to identify the genes of interest for further studies, we studied transcriptomic (Illumina HiSeq) and metabolomic (LC-MS) profiles of the two fruits during five developmental stages (young to ripe). As expected, we identified much higher numbers of significantly differentially regulated genes (DEGs) than metabolites. The highest numbers were identified in pairwise comparisons including the first stage for both species, but total numbers were consistently somewhat lower for the LR. The number of differentially regulated metabolites in pairwise comparisons of developmental stages varied from 66 (stages 3 vs 4) to 133 (stages 2 vs 5) in both species. We identified a number of genes (e.g. AAT1, metE, pip) and metabolites (e.g. rutin, raffinose, galactinol, trehalose, citrulline and DL-arginine) that may be of interest to future functional studies of stress adaptation in plants. As LB is also highly suitable for combating soil desertification and alleviating soil salinity/alkalinity/pollution, its potential for human use may be much wider than its current, highly localized, relevance.

SSR marker development in Clerodendrum trichotomum using transcriptome sequencing

Clerodendrum trichotomum, a member of the Lamiaceae (Verbenaceae) family, is an ornamental plant widely distributed in South Asia. Previous studies have focused primarily on its growth characteristics, stress resistance, and pharmacological applications; however, molecular investigations remain limited. Considering germplasm conservation and the extensive applications of this plant, it is necessary to explore transcriptome resources and SSR makers for C. trichotomum. In the present study, RNA sequencing was used to determine the transcriptome of C. trichotomum. Subsequently, unigene annotations and classifications were obtained, and SSRs were mined with MIcroSAtellite. Finally, primer pairs designed with Oligo 6.0 were selected for polymorphism validation. In total, 127,325,666 high-quality reads were obtained, and 58,345 non-redundant unigenes were generated, of which 36,900 (63.24%) were annotated. Among the annotated unigenes, 35,980 (97.51%) had significant similarity to 607 species in Nr databases. In addition, a total of 6,444 SSRs were identified in 5,530 unigenes, and 200 random primer pairs were designed for polymorphism validation. Furthermore, after primary polymorphism identification, 30 polymorphic primer pairs were selected for the further polymorphism screening, and 200 alleles were identified, 197 of which showed polymorphism. In this work, a large number of unigenes were generated, and numerous SSRs were detected. These findings should be beneficial for further investigations into germplasm conservation and various applications of C. trichotomum. These results should also provide a solid foundation for future molecular biology studies in C. trichotomum.

Impact of Nitrogen Fertilizer Levels on Metabolite Profiling of the Lycium barbarum L. Fruit

The yield and quality of goji (Lycium barbarum L.) fruit are heavily dependent on fertilizer, especially the availability of nitrogen, phosphorus, and potassium (N, P, and K, respectively). In this study, we performed a metabolomic analysis of the response of goji berry to nitrogen fertilizer levels using an Ultra Performance Liquid Chromatography-Electrospray Ionization-Tandem Mass Spectrometry (UPLC-ESI-MS/MS) method. There was no significant difference in the fruit yield or the commodity grade between N0 (42.5 g/plant), N1 (85 g/plant), and N2 (127.5 g/plant). The primary nutrients of the goji berry changed with an increasing nitrogen fertilization. Comparative metabolomic profiling of three nitrogen levels resulted in the identification of 612 metabolites, including amino acids, flavonoids, carbohydrates, organic acids, and lipids/alcohols, among others, of which 53 metabolites (lipids, fatty acids, organic acids, and phenolamides) demonstrated significant changes. These results provide new insights into the molecular mechanisms of the relationship between yield and quality of goji berry and nitrogen fertilizer.

De novo characterization of the Goji berry (Lycium barbarium L.) fruit transcriptome and analysis of candidate genes involved in sugar metabolism under different CO2 concentrations

Goji berry (Lycium barbarum L.) is one of the important economic crops due to its exceptional nutritional value and medicinal benefits. Although reduced sugar levels in goji berry exposed to long-term elevated carbon dioxide (CO2) have been documented, the underlying molecular mechanisms remain unknown. The objective of this study was to explore the transcriptome of goji berry fruit under ambient and elevated CO2 concentrations and further to screen the differentially expressed genes (DEGs) for functions related to sugar metabolism. Fruit samples from goji berry exposed to ambient (400 μmol mol-1) and elevated (700 μmol mol-1) levels of CO2 for 120 days were analyzed for total sugar, carotenoid and flavone analysis. In this study, a reduction in total sugar and carotenoid levels in the fruits grown under elevated CO2 levels were observed. Fruit samples were also used to construct cDNA libraries using a HiSeqTM2500 platform. Consequently, 81,100 unigenes were assembled, of which 35,111 (43.3%) were annotated using various databases. Through DEGs analysis, it was found that 55 genes were upregulated and 18 were down-regulated in response to elevated CO2 treatment. Genes involved in the sugar metabolism and the related pathways were identified by Gene Ontology and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analysis. Furthermore, three genes, LBGAE (Lycium barbarum UDP-glucuronate 4-epimerase), LBGALA (Lycium barbarum alpha-galactosidase) and LBMS (Lycium barbarum malate synthase), associated with sugar metabolism were identified and discussed with respect to the reduction in the total sugar levels along with the enzymes acid invertase (AI), sucrose synthase (SS) and sucrose phosphate synthase (SPS) of the sucrose metabolism. This study can provide gene sources for elucidating the molecular mechanisms of sugar metabolism in the fruit of goji berry under elevated CO2.

Transcriptome Sequencing of Different Avocado Ecotypes: de novo Transcriptome Assembly, Annotation, Identification and Validation of EST-SSR Markers

Avocado (Persea americana Mill.) could be considered as an important tropical and subtropical woody oil crop with high economic and nutritional value. Despite the importance of this species, genomic information is currently unavailable for avocado and closely related congeners. In this study, we generated more than 216 million clean reads from different avocado ecotypes using Illumina HiSeq high-throughput sequencing technology. The high-quality reads were assembled into 154,310 unigenes with an average length of 922 bp. A total of 55,558 simple sequence repeat (SSR) loci detected among the 43,270 SSR-containing unigene sequences were used to develop 74,580 expressed sequence tag (EST)-SSR markers. From these markers, a subset of 100 EST-SSR markers was randomly chosen to identify polymorphic EST-SSR markers in 28 avocado accessions. Sixteen EST-SSR markers with moderate to high polymorphism levels were detected, with polymorphism information contents ranging from 0.33 to 0.84 and averaging 0.63. These 16 polymorphic EST-SSRs could clearly and effectively distinguish the 28 avocado accessions. In summary, our study is the first presentation of transcriptome data of different avocado ecotypes and comprehensive study on the development and analysis of a set of EST-SSR markers in avocado. The application of next-generation sequencing techniques for SSR development is a potentially powerful tool for genetic studies.

Transcriptomic analysis of Lycium ruthenicum Murr. during fruit ripening provides insight into structural and regulatory genes in the anthocyanin biosynthetic pathway

Fruit development in Lycium ruthenicum Murr. involves a succession of physiological and biochemical changes reflecting the transcriptional modulation of thousands of genes. Although recent studies have investigated the dynamic transcriptomic responses during fruit ripening in L. ruthenicum, most have been limited in scope, and thus systematic data representing the structural genes and transcription factors involved in anthocyanin biosynthesis are lacking. In this study, the transcriptomes of three ripening stages associated with anthocyanin accumulation, including S1 (green ripeness stage), S2 (skin color change) and S3 (complete ripeness stage) in L. ruthenicum were investigated using Illumina sequencing. Of a total of 43,573 assembled unigenes, 12,734 were differentially expressed during fruit ripening in L. ruthenicum. Twenty-five significantly differentially expressed structural genes (including PAL, C4H, 4CL, CHS, CHI, F3H, F3'H, F3'5'H, DFR, ANS and UFGT) were identified that might be associated with anthocyanin biosynthesis. Additionally, several transcription factors, including MYB, bHLH, WD40, NAC, WRKY, bZIP and MADS, were correlated with the structural genes, implying their important interaction with anthocyanin biosynthesis-related genes. Our findings provide insight into anthocyanin biosynthesis and regulation patterns in L. ruthenicum and offer a systematic basis for elucidating the molecular mechanisms governing anthocyanin biosynthesis in L. ruthenicum.

Development and application of EST-SSRs markers for analysis of genetic diversity in erect milkvetch (Astragalus adsurgens Pall.)

Erect milkvetch (Astragalus adsurgens Pall.) is a major legume forage plant widely grown in Northern China. However, the lack of molecular markers has limited its research into its genetic diversity and work on germplasm improvement. In this study, a total of 39,163 EST-SSR loci were identified from 30,262 unigene sequences in the erect milkvetch transcriptome using Illumina sequencing. Moreover, 22,367 EST-SSR primer pairs (PPs) were successfully designed. In addition, 100 PPs were synthesized and preliminarily screened in two accessions; of these, 90 were determined to be clear and stable EST-SSR markers. Fifty-one PPs were randomly selected in order to assess the genetic diversity of 27 erect milkvetch accessions. The average polymorphism information content of the 51 PPs was 0.682. Greater genetic diversity was detected in accessions from Inner Mongolia and in the group of landrace and wild erect milkvetch accessions. This study provides an important resource for germplasm improvement and genetic diversity analysis in erect milkvetch.

De Novo assembly, characterization and development of EST-SSRs from Bletilla striata transcriptomes profiled throughout the whole growing period

Bletilla striata is an endangered orchid that has been used for millennia as a medicinal herb, in cosmetics and as a horticultural plant. To construct the first nucleotide database for this species and to develop abundant EST-SSR markers for facilitating further studies, various tissues and organs of plants in the main developmental stages were harvested for mRNA isolation and subsequent RNA sequencing. A total of 106,054,784 clean reads were generated by using Illumina paired-end sequencing technology. The reads were assembled into 127,261 unigenes by the Trinity package; the unigenes had an average length of 612 bp and an N50 of 957 bp. Of these unigenes, 67,494 (51.86%) were annotated in a series of databases. Of these annotated unigenes, 41,818 and 24,615 were assigned to gene ontology categories and clusters of orthologous groups, respectively. Additionally, 20,764 (15.96%) unigenes were mapped onto 275 pathways using the KEGG database. In addition, 25,935 high-quality EST-SSR primer pairs were developed from the 15,433 unigenes by MISA mining. To validate the accuracy of the newly designed markers, 87 of 100 randomly selected primers were effectively amplified; 63 of those yielded PCR products of the expected size, and 25 yielded products with significant amounts of polymorphism among the 4 landraces. Furthermore, the transferability test of the 25 polymorphic markers was performed in 6 individuals of two closely related genus Phalaenopsis and dendrobium. Which results showed a total of 5 markers can successfully amplified among these populations. This research provides a comprehensive nucleotide database and lays a solid foundation for functional gene mining and genomic research in B. striata. The developed EST-SSR primers could facilitate phylogenetic studies and breeding.