Floral transcriptome sequencing for SSR marker development and linkage map construction in the tea plant (Camellia sinensis)

Core collection construction and genetic diversity analysis of tea plant (Camellia sinensis [L:] O. Kuntze) accessions in Huangshan city using SSR markers

Assessing genetic diversity and building a core collection is essential to advancing tea plant breeding. In this study, ten SSR markers exhibiting robust amplification and polymorphism were employed to genotype 292 tea accessions sourced from various regions in Huangshan city. The results revealed significant genetic variation, encompassing 180 alleles. Genetic structure was evaluated using neighbor-joining clustering, principal coordinate analysis, and Structure analyses, which categorized the tea accessions into two primary clusters. The genetic diversity within these clusters demonstrated high similarity, likely due to their close geographical proximity. A core collection was established utilizing Core Hunter software, resulting in the selection of 35% of the accessions to effectively represent the genetic diversity of the entire collection. This core collection comprises 102 tea accessions, preserving a high percentage of allele richness and genetic diversity. This research offers valuable insights into genomics research and the sustainable management of tea plant genetic resources in Huangshan city.

Development of SSR markers for genetic diversity analysis and species identification in Polygonatum odoratum (Mill.) Druce based on transcriptome sequences

Polygonatum odoratum (Mill.) Druce is a well-known traditional Chinese herb belonging to the Polygonatum. However, the understanding of the genetic diversity of this species at the molecular level is limited due to the lack of transcriptomic and genomic information. In this study, 37,387 unigenes were assembled based on the transcriptome sequencing of the rhizome of Polygonatum odoratum (Mill.) Druce., and 11,021 single- sequence repeats (SSR) motifs, mainly consisting of single-nucleotide repeats (44.44%), dinucleotides (31.06%), and trinucleotides (22.59%), were identified. Based on these SSR motifs, 9,987 primer pairs of SSR markers were designed and 68 SSR markers were randomly selected for verification, of which 21 SSR markers showed polymorphisms among the 24 Polygonatum odoratum germplasms. Ninety-four alleles were detected: the observed alleles ranged from 2 to 11, the effective alleles varied from 1.086 8 to 4.916 8, the Shannon diversity index was 0.173 2~1.749 7, and the polymorphism information content PIC ranged from 0.076 7 to 0.803 9. Based on our analysis of genetic diversity (SSR genotypes) and population structure, we divided the 24 germplasm resources into two groups, indicating that the germplasm with similar geographical origins can be grouped together. In addition, the primers 'YZ14' and 'YZ47' could effectively distinguished the related species: Polygonatum kingianum Coll.et Hemsl., Polygonatum sibiricum Red., Polygonatum cyrtonema Hua, Polygonatum zanlanscianense Pamp. and Polygonatum odoratum (Mill.) Druce. This is the first study in which a dataset of expressed sequence tag (EST)-SSR markers is constructed for the Polygonatum odoratum (Mill.) Druce, and these newly developed EST-SSR markers provided a very efficient tool for genetic relationship analysis, species identification and marker-assisted selection breeding of Polygonatum odoratum (Mill.) Druce.

High-density genetic map construction and QTL mapping of a zigzag-shaped stem trait in tea plant (Camellia sinensis)

The highly unique zigzag-shaped stem phenotype in tea plants boasts significant ornamental value and is exceptionally rare. To investigate the genetic mechanism behind this trait, we developed BC1 artificial hybrid populations. Our genetic analysis revealed the zigzag-shaped trait as a qualitative trait. Utilizing whole-genome resequencing, we constructed a high-density genetic map from the BC1 population, incorporating 5,250 SNP markers across 15 linkage groups, covering 3,328.51 cM with an average marker interval distance of 0.68 cM. A quantitative trait locus (QTL) for the zigzag-shaped trait was identified on chromosome 4, within a 61.2 to 97.2 Mb range, accounting for a phenotypic variation explained (PVE) value of 13.62%. Within this QTL, six candidate genes were pinpointed. To better understand their roles, we analyzed gene expression in various tissues and individuals with erect and zigzag-shaped stems. The results implicated CsXTH (CSS0035625) and CsCIPK14 (CSS0044366) as potential key contributors to the zigzag-shaped stem formation. These discoveries lay a robust foundation for future functional genetic mapping and tea plant genetic enhancement.

The Potential Role of Genic-SSRs in Driving Ecological Adaptation Diversity in Caragana Plants

Caragana, a xerophytic shrub genus widely distributed in northern China, exhibits distinctive geographical substitution patterns and ecological adaptation diversity. This study employed transcriptome sequencing technology to investigate 12 Caragana species, aiming to explore genic-SSR variations in the Caragana transcriptome and identify their role as a driving force for environmental adaptation within the genus. A total of 3666 polymorphic genic-SSRs were identified across different species. The impact of these variations on the expression of related genes was analyzed, revealing a significant linear correlation (p < 0.05) between the length variation of 264 polymorphic genic-SSRs and the expression of associated genes. Additionally, 2424 polymorphic genic-SSRs were located in differentially expressed genes among Caragana species. Through weighted gene co-expression network analysis, the expressions of these genes were correlated with 19 climatic factors and 16 plant functional traits in various habitats. This approach facilitated the identification of biological processes associated with habitat adaptations in the studied Caragana species. Fifty-five core genes related to functional traits and climatic factors were identified, including various transcription factors such as MYB, TCP, ARF, and structural proteins like HSP90, elongation factor TS, and HECT. The roles of these genes in the ecological adaptation diversity of Caragana were discussed. Our study identified specific genomic components and genes in Caragana plants responsive to heterogeneous habitats. The results contribute to advancements in the molecular understanding of their ecological adaptation, lay a foundation for the conservation and development of Caragana germplasm resources, and provide a scientific basis for plant adaptation to global climate change.

The MADS-box transcription factor CsAGL9 plays essential roles in seed setting in Camellia sinensis

The number of seed setting (NSS) is an important biological trait that affects tea propagation and yield. In this study, the NSS of an F1 tea population (n = 324) generated via a cross between 'Longjing 43' and 'Baihaozao' was investigated at two locations in two consecutive years. Quantitative trait locus (QTL) mapping of the NSS was performed, and 10 major QTLs were identified. In total, 318 genes were found in these 10 QTLs intervals, and 11 key candidate genes were preliminarily identified. Among them, the MADS-box transcription factor AGAMOUS LIKE 9 (CsAGL9, CSS0037962) located in the most stable QTL (qNSS2) was identified as a key gene affecting the NSS. CsAGL9 overexpression in Arabidopsis promoted early flowering and significantly decreased the length and number of pods and number of seeds per pod. Transcriptome analysis demonstrated that the auxin pathway, a key hormone pathway regulating plant reproduction, was highly affected in the transgenic lines. The auxin pathway was likewise the most prominent in the gene co-expression network study of CsAGL9 in tea plants. In summary, we identified CsAGL9 is essential for seed setting using QTL mapping integrated with RNA-seq, which shed a new light on the mechanism NSS of seed setting in tea plants.

SSR-based molecular diagnosis for Taiwan tea cultivars and its application in identifying cultivar composition of the processed tea

Taiwan specialty teas are produced with distinct manufacturing processes from specific cultivars of tea plants in Camellia. Due to the widespread transplantation of Taiwan tea cultivars and active international trading of tea materials, an accurate and reliable method to identify tea cultivars at the border is vital to protect the image of premium Taiwan specialty teas. In this study, we introduced the Taiwan Tea Variety Identification (TTVID) kit, a capillary electrophoresis-based multiplex PCR assay consisting of 12 simple sequence repeat (SSR) markers. A database composing these 12 SSR loci genotypes in 144 cultivars was established for marker assessment and molecular diagnosis. The power of discrimination on a locus ranged from 0.7894 to 0.966 and the combined match probability of 12 SSR loci was 5.34e-14. Cultivar pairwise comparison among 144 accessions showed that over 90.6% of the pairs had differential genotypes on at least 10 of 12 SSR loci. Further assessment showed that the TTVID kit could unambiguously recognize the cultivars mixed in the loose-leaf teas processed with various degrees of fermentation and roasting. Our results suggested that this TTVID kit effectively identified cultivar composition in loose-leaf tea and is helpful for border control in preventing adulteration and fraud in the Taiwan tea market.

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

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

Transcriptome Analysis of the Late-Acting Self-Incompatibility Associated with RNase T2 Family in Camellia oleifera

The Camellia oil tree (Camellia oleifera Abel.) is an important nonwood forest species in China, and the majority of its cultivars are late-acting self-incompatibility (LSI) types. Although several studies have examined the mechanism of LSI, the process is quite complicated and unclear. In this study, pollen tube growth and fruit setting of two Camellia oil tree cultivars Huashuo (HS) and Huajin (HJ) were investigated after non and self-pollination, and transcriptomic analysis of the ovaries was performed 48 h after self-pollination to identify the potential genes implicated in the LSI of Camellia oil trees. The results showed that the fruit set of HS was significantly higher than that of HJ after self-pollination. Transcriptomic analysis revealed that plant hormone signal transduction, the phosphatidylinositol signaling system, ATP-binding cassette (ABC) transporters, reactive oxygen species (ROS) metabolism, and Ca2+ signaling were mainly contributed in the LSI of reaction of Camellia oil tree. Moreover, nine RNase T2 genes were identified from the transcriptome analysis, which also showed that CoRNase7 participated in the self-incompatibility reaction in HS. Based on phylogenetic analysis, CoRNase6 was closely related to S-RNase from coffee, and CoRNase7 and CoRNase8 were closely related to S-RNase from Camellia sinensis. The 9 RNase T2 genes successfully produced proteins in prokaryotes. Subcellular localization indicated that CoRNase1 and CoRNase5 were cytoplasmic proteins, while CoRNase7 was a plasma membrane protein. These results screened the main metabolic pathways closely related to LSI in Camellia oil tree, and SI signal transduction might be regulated by a large molecular regulatory network. The discovery of T2 RNases provided evidence that Camellia oil tree might be under RNase-based gametophytic self-incompatibility.

Genome-wide characterization and development of SSR markers for genetic diversity analysis in northwestern Himalayas Walnut (Juglans regia L.)

In the present study, we designed and validated genome-wide polymorphic SSR markers (110 SSRs) by mining the walnut genome. A total of 198,924 SSR loci were identified. Among these, successful primers were designed for 162,594 (81.73%) SSR loci. Dinucleotides were the most predominant accounting for 88.40% (175,075) of total SSRs. The SSR frequency was 377.312 SSR/Mb and it showed a decreasing trend from dinucleotide to octanucleotide motifs. We identified 20 highly polymorphic SSR markers and used them to genotype 72 walnut accessions. Over all, we obtained 118 alleles that ranged from 2 to 12 with an average value of 5.9. The higher SSR PIC values indicate their robustness in discriminating walnut genotypes. Heat map, PCA, and population structure categorized 72 walnut genotypes into 2 distinct clusters. The genetic variation within population was higher than among population as inferred by analysis of molecular variance (AMOVA). For walnut improvement, it is necessary to have a large repository of SSRs with high discriminative power. The present study reports 150,000 SSRs, which is the largest SSR repository for this important nut crop. Scientific communities may use this repository for walnut improvement such as QTL mapping, genetic studies, linkage map construction, and marker-assisted selection.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03563-6.

Molecular markers in tea plant (Camellia sinensis): Applications to evolution, genetic identification, and molecular breeding

Tea plants have a long cultivation history in the world, and the beverage (tea) made from its leaves is well known in the world. Due to the characteristics of self-incompatibility, long-term natural and artificial hybridization, tea plants have a very complex genetic background, which make the classification of tea plants unclear. Molecular marker, one type of genetic markers, has the advantages of stable inheritance, large amount of information, and high reliability. The development of molecular marker has facilitated the understanding of complex tea germplasm resources. So far, molecular markers had played important roles in the study of the origin and evolution, the preservation and identification of tea germplasms, and the excellent cultivars breeding of tea plants. However, the information is scattered, making it difficult to understand the advance of molecular markers in tea plants. In this paper, we summarized the development process and types of molecular markers in tea plants. In addition, the application advance of these molecular markers in tea plants was reviewed. Perspectives of molecular markers in tea plants were also systematically provided and discussed. The elaboration of molecular markers in this paper should help us to renew understanding of its application in tea plants.

Genome Survey and SSR Analysis of Camellia nitidissima Chi (Theaceae)

Camellia nitidissima Chi (CNC), a species of golden Camellia, is well known as "the queen of camellias." It is an ornamental, medicinal, and edible plant grown in China. In this study, we conducted a genome survey sequencing analysis and simple sequence repeat (SSR) identification of CNC using the Illumina sequencing platform. The 21-mer analysis predicted its genome size to be 2,778.82 Mb, with heterozygosity and repetition rates of 1.42% and 65.27%, respectively. The CNC genome sequences were assembled into 9,399,197 scaffolds, covering ∼2,910 Mb and an N50 of 869 base pair. Its genomic characteristics were found to be similar to those of Camellia oleifera. In addition, 1,940,616 SSRs were identified from the genome data, including mono-(61.85%), di-(28.71%), tri-(6.51%), tetra-(1.85%), penta-(0.57%), and hexanucleotide motifs (0.51%). We believe these data will provide a useful foundation for the development of novel molecular markers for CNC as well as for further whole-genome sequencing of CNC.

Self-Incompatibility of Camellia weiningensis Y.K. Li

This study compared the pollen tube growth, fruit setting, and seed setting characteristics of Camellia weiningensis Y.K. Li. under self- and cross-pollination to identify its self-incompatibility characteristics and types. C. weiningensis pollen tube growth was observed by fluorescence and scanning electron microscopy, and a field experiment with manual pollination verified fruit and seed setting characteristics. Both self- and cross-pollinated pollen germinated from the stigma. At 72 h after cross-pollination, the pollen tube reached the style base, with tube growth showing a slow-fast-slow pattern. The tube growth speed was maximal, 343.36 μm·h−1, at 12–24 h after pollination. For self-pollination, the pollen did not germinate on the stigma 4 h before pollination. At 12–24 h after pollination, the growth rate was maximal at 263.36 μm·h−1. At 96 h, a small amount of pollen reached the style base and stagnated. The pollen tube end showed callose reactions, such as abnormal swelling, distortion, and brightness. In the field experiment, the fruit setting rate under cross-pollination was 68.5%, while that under self-pollination was 15.3%. When the fruit grew to maturity, the growth dynamics of the transverse and longitudinal diameters showed a “slow-fast-slow”, S-shaped curve. The number of aborted selfed and outcrossed seeds was 13.9 and 4.7, respectively. Thus, C. weiningensis showed self-incompatibility. The self-incompatibility reaction occurred at the style base and represented prezygotic self-incompatibility. The self-incompatibility of C. weiningensis is one of the main reasons for its low seed setting rate, which should be fully considered in cross breeding.

Genome-wide identification and characterization of functionally relevant microsatellite markers from transcription factor genes of Tea (Camellia sinensis (L.) O. Kuntze)

Tea, being one of the most popular beverages requires large set of molecular markers for genetic improvement of quality, yield and stress tolerance. Identification of functionally relevant microsatellite or simple sequence repeat (SSR) marker resources from regulatory “Transcription factor (TF) genes” can be potential targets to expedite molecular breeding efforts. In current study, 2776 transcripts encoding TFs harbouring 3687 SSR loci yielding 1843 flanking markers were identified from traits specific transcriptome resource of 20 popular tea cultivars. Of these, 689 functionally relevant SSR markers were successfully validated and assigned to 15 chromosomes (Chr) of CSS genome. Interestingly, 589 polymorphic markers including 403 core-set of TF-SSR markers amplified 2864 alleles in key TF families (bHLH, WRKY, MYB-related, C2H2, ERF, C3H, NAC, FAR1, MYB and G2-like). Their significant network interactions with key genes corresponding to aroma, quality and stress tolerance suggests their potential implications in traits dissection. Furthermore, single amino acid repeat reiteration in CDS revealed presence of favoured and hydrophobic amino acids. Successful deployment of markers for genetic diversity characterization of 135 popular tea cultivars and segregation in bi-parental population suggests their wider utility in high-throughput genotyping studies in tea.

Development of SSR Databases Available for Both NGS and Capillary Electrophoresis in Apple, Pear and Tea

Developing new varieties in fruit and tea breeding programs is very costly and labor-intensive. Thus, establishing a variety discrimination system is important for protecting breeders’ rights and producers’ profits. Simple sequence repeat (SSR) databases that can be utilized for both next-generation sequencing (SSR-GBS) and polymerase chain reaction–capillary electrophoresis (PCR-CE) would be very useful in variety discrimination. In the present study, SSRs with tri-, tetra- and pentanucleotide repeats were examined in apple, pear and tea. Out of 37 SSRs that showed clear results in PCR-CE, 27 were suitable for SSR-GBS. Among the remaining markers, there was allele dropout for some markers that caused differences between the results of PCR-CE and SSR-GBS. For the selected 27 markers, the alleles detected by SSR-GBS were comparable to those detected by PCR-CE. Furthermore, we developed a computational pipeline for automated genotyping using SSR-GBS by setting a value “α” for each marker, a criterion whether a genotype is homozygous or heterozygous based on allele frequency. The set of 27 markers contains 10, 8 and 9 SSRs for apple, pear and tea, respectively, that are useful for both PCR-CE and SSR-GBS and suitable for automation. The databases help researchers discriminate varieties in various ways depending on sample size, markers and methods.

Development of genic SSR marker resources from RNA-seq data in Camellia japonica and their application in the genus Camellia

Camellia is a genus of flowering plants in the family Theaceae, and several species in this genus have economic importance. Although a great deal of molecular makers has been developed for molecular assisted breeding in genus Camellia in the past decade, the number of simple sequence repeats (SSRs) publicly available for plants in this genus is insufficient. In this study, a total of 28,854 potential SSRs were identified with a frequency of 4.63 kb. A total of 172 primer pairs were synthesized and preliminarily screened in 10 C. japonica accessions, and of these primer pairs, 111 were found to be polymorphic. Fifty-one polymorphic SSR markers were randomly selected to perform further analysis of the genetic relationships of 89 accessions across the genus Camellia. Cluster analysis revealed major clusters corresponding to those based on taxonomic classification and geographic origin. Furthermore, all the genotypes of C. japonica separated and consistently grouped well in the genetic structure analysis. The results of the present study provide high-quality SSR resources for molecular genetic breeding studies in camellia plants.

QTL Mapping for Leaf Area of Tea Plants (Camellia sinensis) Based on a High-Quality Genetic Map Constructed by Whole Genome Resequencing

High-quality genetic maps play important roles in QTL mapping and molecular marker-assisted breeding. Tea leaves are not only important vegetative organs but are also the organ for harvest with important economic value. However, the key genes and genetic mechanism of regulating leaf area have not been clarified. In this study, we performed whole-genome resequencing on "Jinxuan," "Yuncha 1" and their 96 F1 hybrid offspring. From the 1.84 Tb of original sequencing data, abundant genetic variation loci were identified, including 28,144,625 SNPs and 2,780,380 indels. By integrating the markers of a previously reported genetic map, a high-density genetic map consisting of 15 linkage groups including 8,956 high-quality SNPs was constructed. The total length of the genetic map is 1,490.81 cM, which shows good collinearity with the genome. A total of 25 representative markers (potential QTLs) related to leaf area were identified, and there were genes differentially expressed in large and small leaf samples near these markers. GWAS analysis further verified the reliability of QTL mapping. Thirty-one pairs of newly developed indel markers located near these potential QTLs showed high polymorphism and had good discrimination between large and small leaf tea plant samples. Our research will provide necessary support and new insights for tea plant genetic breeding, quantitative trait mapping and yield improvement.

Inheritance and quantitative trait loci analyses of the anthocyanins and catechins of Camellia sinensis cultivar 'Ziyan' with dark-purple leaves

Owing to the potential health benefits, anthocyanin-rich teas (Camellia sinensis) have attracted interest over the past decade. Previously, we developed the cultivar 'Ziyan,' which has dark-purple leaves because of the accumulation of a high amount of anthocyanins. In this study, we performed a genetic analysis of this anthocyanin-rich tea cultivar and 176 of its naturally pollinated offspring. For two consecutive years, we quantified the anthocyanins and catechins of 'Ziyan' and the offspring population. While >60% of the offspring accumulated less than half of the amount of anthocyanins of 'Ziyan,' 17 (2018) and 15 (2019) individuals exceeded 'Ziyan' in anthocyanin content. A negative correlation between anthocyanin and total catechin content (r = -0.59, P < 0.001) was observed. The population was genotyped with 131 SSR markers spanning all linkage groups of the C. sinensis genome. Kruskal-Wallis tests identified 10 markers significantly associated with anthocyanins, catechins and their ratios in both years. Quantitative trait locus (QTL) analyses using the interval mapping method detected 13 QTLs, suggesting the dark-purple trait of 'Ziyan' is because of the pyramiding of anthocyanin-promoting alleles on at least five linkage groups. Two genetic loci reversely related to anthocyanin and total catechin contents were identified. This study provides valuable information for genetic improvement of purple tea cultivars and for fine-mapping related genes.

De novo Transcriptome Assembly, Gene Annotation and SSR Marker Development in the Moon Seed Genus Menispermum (Menispermaceae)

The moonseed genus Menispermum L. (Menispermaceae) is disjunctly distributed in East Asia and eastern North America. Although Menispermum has important medicinal value, genetic and genomic information is scarce, with very few available molecular markers. In the current study, we used Illumina transcriptome sequencing and de novo assembly of the two Menispermum species to obtain in-depth genetic knowledge. From de novo assembly, 53,712 and 78,921 unigenes were generated for M. canadense and M. dauricum, with 37,527 (69.87%) and 55,211 (69.96%) showing significant similarities against the six functional databases, respectively. Moreover, 521 polymorphic EST-SSRs were identified. Of them, 23 polymorphic EST-SSR markers were selected to investigate the population genetic diversity within the genus. The newly developed EST-SSR markers also revealed high transferability among the three examined Menispermaceae species. Overall, we provide the very first transcriptomic analyses of this important medicinal genus. In addition, the novel microsatellite markers developed here will aid future studies on the population genetics and phylogeographic patterns of Menispermum at the intercontinental geographical scale.

TeaMiD: a comprehensive database of simple sequence repeat markers of tea

Tea is a highly cross-pollinated, woody, perennial tree. High heterozygosity combined with a long gestational period makes conventional breeding a cumbersome process. Therefore, marker-assisted breeding is a better alternative approach when compared with conventional breeding. Considering the large genome size of tea (~3 Gb), information about simple sequence repeat (SSR) is scanty. Thus, we have taken advantage of the recently published tea genomes to identify large numbers of SSR markers in the tea. Besides the genomic sequences, we identified SSRs from the other publicly available sequences such as RNA-seq, GSS, ESTs and organelle genomes (chloroplasts and mitochondrial) and also searched published literature to catalog validated set of tea SSR markers. The complete exercise yielded a total of 935 547 SSRs. Out of the total, 82 SSRs were selected for validation among a diverse set of tea genotypes. Six primers (each with four to six alleles, an average of five alleles per locus) out of the total 27 polymorphic primers were used for a diversity analysis in 36 tea genotypes with mean polymorphic information content of 0.61-0.76. Finally, using all the information generated in this study, we have developed a user-friendly database (TeaMiD; http://indianteagenome.in:8080/teamid/) that hosts SSR from all the six resources including three nuclear genomes of tea and transcriptome sequences of 17 Camellia wild species. Database URL: http://indianteagenome.in:8080/teamid/.

De novo assembly, transcriptome characterization, and simple sequence repeat marker development in duckweed Lemna gibba

Lemna gibba is a species of duckweed showing great potential in bioenergy production and wastewater treatment. However, the relevant transcriptomic and genomic resources are very limited for this species, which dramatically hinders its genetic diversity and genome mapping researches. In this work, ~ 233.5 million clean reads were generated from L. gibba by Illumina paired-end sequencing, and subsequently they were de novo assembled into 131,870 unigenes, of which 61,622 were annotated and 43,319 were expressed with Fragments Per Kilobase of transcript per Million fragments mapped (FPKM) > 5. In total, 19,297 simple sequence repeats (SSRs) were identified from 15,261 SSR-containing unigenes. Dinucleotide (78.4%) were the most abundant SSRs, followed by tri- (14.9%), tetra- (4.1%), and penta-nucleotides (1.5%). The top three motifs were AG/CT (69.9%), AC/GT (6.5%), and ATC/ATG (4.9%). Further analysis revealed that the presence of SSR motif was independent of the expression level for a given gene. Based on the sequence of these SSR-containing unigenes, a total of 10,292 SSR markers were developed, of which only 2671 were further retained after removing those derived from unannotated or extra-low expressed (e.g., FPKM ≤ 5) unigenes. Finally, a subset of 70 SSR markers was randomly selected and examined in nine diverse L. gibba genotypes for the PCR amplification and polymorphism, as well as in other duckweed species for the inter-specifically amplifiability. This work is the first report on the transcriptome-based large-scale SSR markers development and analysis in L. gibba. The transcriptome generated and the SSR markers developed in this work will provide a valuable resource for genetic diversity assessment in L. gibba and also for species relationship investigation in Lemnaceae family.

Optimized sequencing depth and de novo assembler for deeply reconstructing the transcriptome of the tea plant, an economically important plant species

Background

Tea is the oldest and among the world’s most popular non-alcoholic beverages, which has important economic, health and cultural values. Tea is commonly produced from the leaves of tea plants (Camellia sinensis), which belong to the genus Camellia of family Theaceae. In the last decade, many studies have generated the transcriptomes of tea plants at different developmental stages or under abiotic and/or biotic stresses to investigate the genetic basis of secondary metabolites that determine tea quality. However, these results exhibited large differences, particularly in the total number of reconstructed transcripts and the quality of the assembled transcriptomes. These differences largely result from limited knowledge regarding the optimized sequencing depth and assembler for transcriptome assembly of structurally complex plant species genomes.

Results

We employed different amounts of RNA-sequencing data, ranging from 4 to 84 Gb, to assemble the tea plant transcriptome using five well-known and representative transcript assemblers. Although the total number of assembled transcripts increased with increasing sequencing data, the proportion of unassembled transcripts became saturated as revealed by plant BUSCO datasets. Among the five representative assemblers, the Bridger package shows the best performance in both assembly completeness and accuracy as evaluated by the BUSCO datasets and genome alignment. In addition, we showed that Bridger and BinPacker harbored the shortest runtimes followed by SOAPdenovo and Trans-ABySS.

Conclusions

The present study compares the performance of five representative transcript assemblers and investigates the key factors that affect the assembly quality of the transcriptome of the tea plants. This study will be of significance in helping the tea research community obtain better sequencing and assembly of tea plant transcriptomes under conditions of interest and may thus help to answer major biological questions currently facing the tea industry.

Development of polymorphic EST-SSR markers and characterization of the autotetraploid genome of sainfoin (Onobrychis viciifolia)

Background

Sainfoin (Onobrychis viciifolia) is a highly nutritious, tannin-containing, and tetraploid forage legume. Due to the lack of detailed transcriptomic and genomic information on this species, genetic and breeding projects for sainfoin improvement have been significantly hindered.

Methods

In this study, a total of 24,630,711 clean reads were generated from 14 different sainfoin tissues using Illumina paired-end sequencing technology and deposited in the NCBI SRA database (SRX3763386). From these clean reads, 77,764 unigene sequences were obtained and 6,752 EST-SSRs were identified using de novo assembly. A total of 2,469 primer pairs were designed, and 200 primer pairs were randomly selected to analyze the polymorphism in five sainfoin wild accessions.

Results

Further analysis of 40 sainfoin individuals from the five wild populations using 61 EST-SSR loci showed that the number of alleles per locus ranged from 4 to 15, and the expected heterozygosity varied from 0.55 to 0.91. Additionally, by counting the EST-SSR band number and sequencing the three or four bands in one sainfoin individual, sainfoin was confirmed to be autotetraploid. This finding provides a high level of information about this plant.

Discussion

Through this study, 61 EST-SSR markers were successfully developed and shown to be useful for genetic studies and investigations of population genetic structures and variabilities among different sainfoin accessions.

A coupled role for CsMYB75 and CsGSTF1 in anthocyanin hyperaccumulation in purple tea

Cultivars of purple tea (Camellia sinensis) that accumulate anthocyanins in place of catechins are currently attracting global interest in their use as functional health beverages. RNA-seq of normal (LJ43) and purple Zijuan (ZJ) cultivars identified the transcription factor CsMYB75 and phi (F) class glutathione transferase CsGSTF1 as being associated with anthocyanin hyperaccumulation. Both genes mapped as a quantitative trait locus (QTL) to the purple bud leaf color (BLC) trait in F1 populations, with CsMYB75 promoting the expression of CsGSTF1 in transgenic tobacco (Nicotiana tabacum). Although CsMYB75 elevates the biosynthesis of both catechins and anthocyanins, only anthocyanins accumulate in purple tea, indicating selective downstream regulation. As glutathione transferases in other plants are known to act as transporters (ligandins) of flavonoids, directing them for vacuolar deposition, the role of CsGSTF1 in selective anthocyanin accumulation was investigated. In tea, anthocyanins accumulate in multiple vesicles, with the expression of CsGSTF1 correlated with BLC, but not with catechin content, in diverse germplasm. Complementation of the Arabidopsis tt19-8 mutant, which is unable to express the orthologous ligandin AtGSTF12, restored anthocyanin accumulation, but did not rescue the transparent testa phenotype, confirming that CsGSTF1 did not function in catechin accumulation. Consistent with a ligandin function, transient expression of CsGSTF1 in Nicotiana occurred in the nucleus, cytoplasm and membrane. Furthermore, RNA-Seq of the complemented mutants exposed to 2% sucrose as a stress treatment showed unexpected roles for anthocyanin accumulation in affecting the expression of genes involved in redox responses, phosphate homeostasis and the biogenesis of photosynthetic components, as compared with non-complemented plants.

Understanding Camellia sinensis using Omics Technologies along with Endophytic Bacteria and Environmental Roles on Metabolism: A Review

Camellia sinensis is the most consumed beverage worldwide. It contains a wide variety of secondary metabolites, such as alkaloids, saponins, tannins, catechins, and polyphenols, generated through a condensation reaction of cinnamic acid with three malonyl-CoA groups. In addition to the metabolic processes occurring within this plant, there are also some plant-associated bacterial endophytes. These bacteria reside in the living tissues of the host plants without causing any harmful effect to them, thereby stimulating secondary metabolite production with a diverse range of biological effects. Omics technologies reveal understanding of the biological phenomena of transcriptomics, proteomics, and metabolomics. In this sense, the present review aims to provide a comprehensive review of various methods used to identify distinct plant compounds, namely transcriptomic, proteomic, and metabolomic analysis. The role of endophytic bacteria in C. sinensis metabolism, and C. sinensis antioxidant and antimicrobial effects, are also carefully highlighted.

High-density SNP linkage map construction and QTL mapping for flavonoid-related traits in a tea plant (Camellia sinensis) using 2b-RAD sequencing

BACKGROUND

Flavonoids are important components that confer upon tea plants a unique flavour and health functions. However, the traditional breeding method for selecting a cultivar with a high or unique flavonoid content is time consuming and labour intensive. High-density genetic map construction associated with quantitative trait locus (QTL) mapping provides an effective way to facilitate trait improvement in plant breeding. In this study, an F1 population (LJ43×BHZ) was genotyped using 2b-restriction site-associated DNA (2b-RAD) sequencing to obtain massive single nucleotide polymorphism (SNP) markers to construct a high-density genetic map for a tea plant. Furthermore, QTLs related to flavonoids were identified using our new genetic map.

RESULTS

A total of 13,446 polymorphic SNP markers were developed using 2b-RAD sequencing, and 4,463 of these markers were available for constructing the genetic linkage map. A 1,678.52-cM high-density map at an average interval of 0.40 cM with 4,217 markers, including 427 frameset simple sequence repeats (SSRs) and 3,800 novel SNPs, mapped into 15 linkage groups was successfully constructed. After QTL analysis, a total of 27 QTLs related to flavonoids or caffeine content (CAF) were mapped to 8 different linkage groups, LG01, LG03, LG06, LG08, LG10, LG11, LG12, and LG13, with an LOD from 3.14 to 39.54, constituting 7.5% to 42.8% of the phenotypic variation.

CONCLUSIONS

To our knowledge, the highest density genetic map ever reported was constructed since the largest mapping population of tea plants was adopted in present study. Moreover, novel QTLs related to flavonoids and CAF were identified based on the new high-density genetic map. In addition, two markers were located in candidate genes that may be involved in flavonoid metabolism. The present study provides valuable information for gene discovery, marker-assisted selection breeding and map-based cloning for functional genes that are related to flavonoid content in tea plants.

Unraveling the Roles of Regulatory Genes during Domestication of Cultivated Camellia: Evidence and Insights from Comparative and Evolutionary Genomics

With the increasing power of DNA sequencing, the genomics-based approach is becoming a promising resolution to dissect the molecular mechanism of domestication of complex traits in trees. Genus Camellia possesses rich resources with a substantial value for producing beverage, ornaments, edible oil and more. Currently, a vast number of genetic and genomic research studies in Camellia plants have emerged and provided an unprecedented opportunity to expedite the molecular breeding program. In this paper, we summarize the recent advances of gene expression and genomic resources in Camellia species and focus on identifying genes related to key economic traits such as flower and fruit development and stress tolerances. We investigate the genetic alterations and genomic impacts under different selection programs in closely related species. We discuss future directions of integrating large-scale population and quantitative genetics and multiple omics to identify key candidates to accelerate the breeding process. We propose that future work of exploiting the genomic data can provide insights related to the targets of domestication during breeding and the evolution of natural trait adaptations in genus Camellia.

Insight into Catechins Metabolic Pathways of Camellia sinensis Based on Genome and Transcriptome Analysis

Tea is an important economic crop with a 3.02 Gb genome. It accumulates various bioactive compounds, especially catechins, which are closely associated with tea flavor and quality. Catechins are biosynthesized through the phenylpropanoid and flavonoid pathways, with 12 structural genes being involved in their synthesis. However, we found that in Camellia sinensis the understanding of the basic profile of catechins biosynthesis is still unclear. The gene structure, locus, transcript number, transcriptional variation, and function of multigene families have not yet been clarified. Our previous studies demonstrated that the accumulation of flavonoids in tea is species, tissue, and induction specific, which indicates that gene coexpression patterns may be involved in tea catechins and flavonoids biosynthesis. In this paper, we screened candidate genes of multigene families involved in the phenylpropanoid and flavonoid pathways based on an analysis of genome and transcriptome sequence data. The authenticity of candidate genes was verified by PCR cloning, and their function was validated by reverse genetic methods. In the present study, 36 genes from 12 gene families were identified and were accessed in the NCBI database. During this process, some intron retention events of the CsCHI and CsDFR genes were found. Furthermore, the transcriptome sequencing of various tea tissues and subcellular location assays revealed coexpression and colocalization patterns. The correlation analysis showed that CsCHIc, CsF3'H, and CsANRb expression levels are associated significantly with the concentration of soluble PA as well as the expression levels of CsPALc and CsPALf with the concentration of insoluble PA. This work provides insights into catechins metabolism in tea and provides a foundation for future studies.

Mining and Development of Novel SSR Markers Using Next Generation Sequencing (NGS) Data in Plants

Microsatellites, or simple sequence repeats (SSRs), are one of the most informative and multi-purpose genetic markers exploited in plant functional genomics. However, the discovery of SSRs and development using traditional methods are laborious, time-consuming, and costly. Recently, the availability of high-throughput sequencing technologies has enabled researchers to identify a substantial number of microsatellites at less cost and effort than traditional approaches. Illumina is a noteworthy transcriptome sequencing technology that is currently used in SSR marker development. Although 454 pyrosequencing datasets can be used for SSR development, this type of sequencing is no longer supported. This review aims to present an overview of the next generation sequencing, with a focus on the efficient use of de novo transcriptome sequencing (RNA-Seq) and related tools for mining and development of microsatellites in plants.

De Novo Assembly of Transcriptome and Development of Novel EST-SSR Markers in Rhododendron rex Lévl. through Illumina Sequencing

Transcriptome sequences generated by next-generation sequencing (NGS) technologies can be utilized to rapidly detect and characterize a large number of gene-based microsatellites from different plants. Rhododendron rex Lévl. is a perennial woody species from the family Ericaceae and an endangered plant with high ornamental value endemic to Southwestern China. Nevertheless, the genetic and genomic information of R. rex remain unknown. In this study, we performed transcriptome sequencing for R. rex leaf samples, and generated large transcript sequences for functional characterization and development gene-associated SSR markers. A total of 164,242 unigenes were assembled and 115,089 (70.07%) unigenes were successfully annotated in public databases. In addition, a total of 15,314 potential EST-SSRs were identified, and the frequency of SSRs in the R. rex unigenes was 9.32%, with an average of one EST-SSR per 5.65 kb. The most abundant type was repeated di-nucleotide (54.63%), followed by mono- (26.03%) and tri-nucleotide (18.51%) repeats. Based on the SSR-containing sequence, 100 primer pairs were randomly selected and synthesized and used for assessment of the polymorphism. Thirty-six primer pairs were polymorphic and revealed polymorphism among 20 individuals from four R. rex populations. A total of 197 alleles were identified, with an average of 5.472 alleles per locus. The Polymorphism Information Content ranged from 0.154 to 0.870, with a mean of 0.482. The newly developed EST-SSR markers exhibited high transferability (58.33-83.33%) among the six subgenera. Thus, these novel EST-SSR markers developed would provide valuable sequence resources for population structure, genetic diversity analysis, and genetic resource assessments of R. rex and its related species.

Natural products against Alzheimer's disease: Pharmaco-therapeutics and biotechnological interventions

Alzheimer's disease (AD) is a severe, chronic and progressive neurodegenerative disease associated with memory and cognition impairment ultimately leading to death. It is the commonest reason of dementia in elderly populations mostly affecting beyond the age of 65. The pathogenesis is indicated by accumulation of the amyloid-beta (Aβ) plaques and neurofibrillary tangles (NFT) in brain tissues and hyperphosphorylation of tau protein in neurons. The main cause is considered to be the formation of reactive oxygen species (ROS) due to oxidative stress. The current treatment provides only symptomatic relief by offering temporary palliative therapy which declines the rate of cognitive impairment associated with AD. Inhibition of the enzyme acetylcholinesterase (AChE) is considered as one of the major therapeutic strategies offering only symptomatic relief and moderate disease-modifying effect. Other non-cholinergic therapeutic approaches include antioxidant and vitamin therapy, stem cell therapy, hormonal therapy, use of antihypertensive or lipid-lowering medications and selective phosphodiesterase (PDE) inhibitors, inhibition of β-secretase and γ-secretase and Aβ aggregation, inhibition of tau hyperphosphorylation and intracellular NFT, use of nonsteroidal anti-inflammatory drugs (NSAIDs), transition metal chelators, insulin resistance drugs, etanercept, brain-derived neurotrophic factor (BDNF) etc. Medicinal plants have been reported for possible anti-AD activity in a number of preclinical and clinical trials. Ethnobotany, being popular in China and in the Far East and possibly less emphasized in Europe, plays a substantial role in the discovery of anti-AD agents from botanicals. Chinese Material Medica (CMM) involving Chinese medicinal plants has been used traditionally in China in the treatment of AD. Ayurveda has already provided numerous lead compounds in drug discovery and many of these are also undergoing clinical investigations. A number of medicinal plants either in their crude forms or as isolated compounds have exhibited to reduce the pathological features associated with AD. In this present review, an attempt has been made to elucidate the molecular mode of action of various plant extracts, phytochemicals and traditional herbal formulations investigated against AD as reported in various preclinical and clinical tests. Herbal synergism often found in polyherbal formulations were found effective to combat disease heterogeneity as found in complex pathogenesis of AD. Finally a note has been added to describe biotechnological improvement, genetic and genomic resources and mathematical and statistical techniques for empirical model building associated with anti-AD plant secondary metabolites and their source botanicals.

Gene expression analysis of bud and leaf color in tea

Purple shoot tea attributing to the high anthocyanin accumulation is of great interest for its wide health benefits. To better understand potential mechanisms involved in purple buds and leaves formation in tea plants, we performed transcriptome analysis of six green or purple shoot tea individuals from a F1 population using the Illumina sequencing method. Totally 292 million RNA-Seq reads were obtained and assembled into 112,233 unigenes, with an average length of 759 bp and an N50 of 1081 bp. Moreover, totally 2193 unigenes showed significant differences in expression levels between green and purple tea samples, with 1143 up- and 1050 down-regulated in the purple teas. Further real time PCR analysis confirmed RNA-Seq results. Our study identified 28 differentially expressed transcriptional factors and A CsMYB gene was found to be highly similar to AtPAP1 in Arabidopsis. Further analysis of differentially expressed genes involved in anthocyanin biosynthesis and transportation showed that the late biosynthetic genes and genes involved in anthocyanin transportation were largely affected but the early biosynthetic genes were less or none affected. Overall, the identification of a large number of differentially expressed genes offers a global view of the potential mechanisms associated with purple buds and leaves formation, which will facilitate molecular breeding in tea plants.

An efficient identification strategy of clonal tea cultivars using long-core motif SSR markers

Microsatellites, or simple sequence repeats (SSRs), especially those with long-core motifs (tri-, tetra-, penta-, and hexa-nucleotide) represent an excellent tool for DNA fingerprinting. SSRs with long-core motifs are preferred since neighbor alleles are more easily separated and identified from each other, which render the interpretation of electropherograms and the true alleles more reliable. In the present work, with the purpose of characterizing a set of core SSR markers with long-core motifs for well fingerprinting clonal cultivars of tea (Camellia sinensis), we analyzed 66 elite clonal tea cultivars in China with 33 initially-chosen long-core motif SSR markers covering all the 15 linkage groups of tea plant genome. A set of 6 SSR markers were conclusively selected as core SSR markers after further selection. The polymorphic information content (PIC) of the core SSR markers was >0.5, with ≤5 alleles in each marker containing 10 or fewer genotypes. Phylogenetic analysis revealed that the core SSR markers were not strongly correlated with the trait 'cultivar processing-property'. The combined probability of identity (PID) between two random cultivars for the whole set of 6 SSR markers was estimated to be 2.22 × 10(-5), which was quite low, confirmed the usefulness of the proposed SSR markers for fingerprinting analyses in Camellia sinensis. Moreover, for the sake of quickly discriminating the clonal tea cultivars, a cultivar identification diagram (CID) was subsequently established using these core markers, which fully reflected the identification process and provided the immediate information about which SSR markers were needed to identify a cultivar chosen among the tested ones. The results suggested that long-core motif SSR markers used in the investigation contributed to the accurate and efficient identification of the clonal tea cultivars and enabled the protection of intellectual property.

Renaissance in phytomedicines: promising implications of NGS technologies

Medicinal plant research is growing significantly in faith to discover new and more biologically compatible phytomedicines. Deposition of huge genome/trancriptome sequence data assisted by NGS technologies has revealed the new possibilities for producing upgraded bioactive molecules in medicinal plants. Growing interest of investors and consumers in the herbal drugs raises the need for extensive research to open the facts and details of every inch of life canvas of medicinal plants to produce improved quality of phytomedicines. As in agriculture crops, knowledge emergence from medicinal plant's genome/transcriptome, can be used to assure their amended quality and these improved varieties are then transported to the fields for cultivation. Genome studies generate huge sequence data which can be exploited further for obtaining information regarding genes/gene clusters involved in biosynthesis as well as regulation. This can be achieved rapidly at a very large scale with NGS platforms. Identification of new RNA molecules has become possible, which can lead to the discovery of novel compounds. Sequence information can be combined with advanced phytochemical and bioinformatics tools to discover functional herbal drugs. Qualitative and quantitative analysis of small RNA species put a light on the regulatory aspect of biosynthetic pathways for phytomedicines. Inter or intra genomic as well as transcriptomic interactive processes for biosynthetic pathways can be elucidated in depth. Quality management of herbal material will also become rapid and high throughput. Enrichment of sequence information will be used to engineer the plants to get more efficient phytopharmaceuticals. The present review comprises of role of NGS technologies to boost genomic studies of pharmaceutically important plants and further, applications of sequence information aiming to produce enriched phytomedicines. Emerging knowledge from the medicinal plants genome/transcriptome can give birth to deep understanding of the processes responsible for biosynthesis of medicinally important compounds.

Transcriptome analysis reveals self-incompatibility in the tea plant (Camellia sinensis) might be under gametophytic control

BACKGROUND

Self-incompatibility (SI) is under genetic control and prevents inbreeding depression in angiosperms. SI mechanisms are quite complicated and still poorly understood in many plants. Tea (Camellia sinensis L.) belonging to the family of Theaceae, exhibits high levels of SI and high heterozygosity. Uncovering the molecular basis of SI of the tea plant may enhance breeding and simplify genomics research for the whole family.

RESULTS

The growth of pollen tubes following selfing and crossing was observed using fluorescence microscopy. Self-pollen tubes grew slower than cross treatments from 24 h to 72 h after pollination. RNA-seq was employed to explore the molecular mechanisms of SI and to identify SI-related genes in C. sinensis. Self and cross-pollinated styles were collected at 24 h, 48 h and 72 h after pollination. Six RNA-seq libraries (SP24, SP48, SP72, CP24 CP48 and CP72; SP = self-pollinated, CP = cross-pollinated) were constructed and separately sequenced. In total, 299.327 million raw reads were generated. Following assembly, 63,762 unigenes were identified, and 27,264 (42.76 %) unigenes were annotated in five public databases: NR, KOG, KEGG, Swiss-Port and GO. To identify SI-related genes, the fragments per kb per million mapped reads (FPKM) values of each unigene were evaluated. Comparisons of CP24 vs. SP24, CP48 vs. SP48 and CP72 vs. SP72 revealed differential expression of 3,182, 3,575 and 3,709 genes, respectively. Consequently, several ubiquitin-mediated proteolysis, Ca(2+) signaling, apoptosis and defense-associated genes were obtained. The temporal expression pattern of genes following CP and SP was analyzed; 6 peroxidase, 1 polyphenol oxidase and 7 salicylic acid biosynthetic process-related genes were identified. The RNA-seq data were validated by qRT-PCR of 15 unigenes. Finally, a unigene (CL25983Contig1) with strong homology to the S-RNase was analyzed. It was mainly expressed in styles, with dramatically higher expression in self-pollinated versus cross-pollinated tissues at 24 h post-pollination.

CONCLUSIONS

The present study reports the transcriptome of styles after cross- and self-pollination in tea and offers novel insights into the molecular mechanism behind SI in C. sinensis. We believe that this RNA-seq dataset will be useful for improvement in C. sinensis as well as other plants in the Theaceae family.

Development and cross-species transferability of EST-SSR markers in Siberian wildrye (Elymus sibiricus L.) using Illumina sequencing

Siberian wildrye (Elymus sibiricus L.) is a perennial, self-fertilizing grass that plays an important role in animal husbandry and environmental sustenance. However, the transcriptomic and genomic information on this species is very limited, which hinders genetic and breeding studies. In the present study, 76,686,804 clean reads were generated from 11 different tissue samples of E. sibiricus by Illumina paired-end sequencing, and the reads were deposited into the NCBI SRA database (SRX574376). A total of 8,769 EST-SSRs were identified from 94,458 unigene sequences, which were obtained by de novo assembly. Moreover, 1,078 primer pairs were successfully designed, and 500 pairs were randomly selected to assess polymorphisms in 15 E. sibiricus accessions. A total of 112 primer pairs were polymorphic, and the polymorphism information content (PIC) values ranged from 0.39 to 0.81, indicating a high level of informativeness. Furthermore, these 112 polymorphic primer pairs were used to evaluate the transferability to 13 other related species, and 55 EST-SSR markers were found to be polymorphic among these 13 Elymus species. This study collected the global sequence data for E. sibiricus, and the newly developed markers will prove valuable in facilitating genetic diversity in E. sibiricus and related Elymus species.

Biotechnological advances in tea (Camellia sinensis [L.] O. Kuntze): a review

This article presents a comprehensive review on the success and limitations of biotechnological approaches aimed at genetic improvement of tea with a purpose to explore possibilities to address challenging areas. Tea is a woody perennial tree with a life span of more than 100 years. Conventional breeding of tea is slow and limited primarily to selection which leads to narrowing down of its genetic base. Harnessing the benefits of wild relatives has been negligible due to low cross-compatibility, genetic drag and undesirable alleles for low yield. Additionally, being a recalcitrant species, in vitro propagation of tea is constrained too. Nevertheless, maneuvering with tissue/cell culture techniques, a considerable success has been achieved in the area of micropropagation, somatic embryogenesis as well as genetic transformation. Besides, use of molecular markers, "expressomics" (transcriptomics, proteomics, metabolomics), map-based cloning towards construction of physical maps, generation of expressed sequenced tags (ESTs) have facilitated the identification of QTLs and discovery of genes associated with abiotic or biotic stress tolerance and agronomic traits. Furthermore, the complete genome (or at least gene space) sequence of tea is expected to be accessible in the near future which will strengthen combinational approaches for improvement of tea. This review presents a comprehensive account of the success and limitations of the biotechnological tools and techniques hitherto applied to tea and its wild relatives. Expectedly, this will form a basis for making further advances aimed at genetic improvement of tea in particular and of economically important woody perennials in general.

Transcriptome Analysis Identifies Candidate Genes Related to Triacylglycerol and Pigment Biosynthesis and Photoperiodic Flowering in the Ornamental and Oil-Producing Plant, Camellia reticulata (Theaceae)

Camellia reticulata, which is native to Southwest China, is famous for its ornamental flowers and high-quality seed oil. However, the lack of genomic information for this species has largely hampered our understanding of its key pathways related to oil production, photoperiodic flowering process and pigment biosynthesis. Here, we first sequenced and characterized the transcriptome of a diploid C. reticulata in an attempt to identify genes potentially involved in triacylglycerol biosynthesis (TAGBS), photoperiodic flowering, flavonoid biosynthesis (FlaBS), carotenoid biosynthesis (CrtBS) pathways. De novo assembly of the transcriptome provided a catalog of 141,460 unigenes with a total length of ~96.1 million nucleotides (Mnt) and an N50 of 1080 nt. Of them, 22,229 unigenes were defined as differentially expressed genes (DEGs) across five sequenced tissues. A large number of annotated genes in C. reticulata were found to have been duplicated, and differential expression patterns of these duplicated genes were commonly observed across tissues, such as the differential expression of SOC1_a, SOC1_b, and SOC1_c in the photoperiodic flowering pathway. Up-regulation of SAD_a and FATA genes and down-regulation of FAD2_a gene in the TAGBS pathway in seeds may be relevant to the ratio of monounsaturated fatty acid (MUFAs) to polyunsaturated fatty acid (PUFAs) in seed oil. MYBF1, a transcription regulator gene of the FlaBS pathway, was found with great sequence variation and alteration of expression patterns, probably resulting in functionally evolutionary differentiation in C. reticulata. MYBA1_a and some anthocyanin-specific biosynthetic genes in the FlaBS pathway were highly expressed in both flower buds and flowers, suggesting important roles of anthocyanin biosynthesis in flower development. Besides, a total of 40,823 expressed sequence tag simple sequence repeats (EST-SSRs) were identified in the C. reticulata transcriptome, providing valuable marker resources for further basic and applied researches on this economically important Camellia plant.

Large-scale developing of simple sequence repeat markers and probing its correlation with ramie (Boehmeria nivea L.) fiber quality

Marker-assisted selection is an important component of the discipline of molecular breeding. Using DNA markers to assist in plant breeding, the efficiency and precision could be greatly increased. However, the scarcity number of identified DNA markers has hindered the research and the breeding process of ramie (Boehmeria nivea L.) in many aspects, especially fiber quality, one of the top-priority breeding objectives of ramie. In this study, 4230 SSR loci were identified in 3969 unigenes (6.80 % of 58,369), which were de novo assembled from the transcriptome involving different ramie fiber developmental stages. Among these SSRs, the dinucleotides (1599, 37.80 %) and trinucleotides (772, 18.25 %) were most abundant; the motifs AG/CT (1140, 26.94 %), AT/AT (407, 9.62 %) and AGA/TCT (246, 8.31 %) comprised the three most abundant repeats. A total of 2431 primer pairs were designed flanking the SSRs and 1050 of them were employed in PCR amplification for their usefulness using three ramie cultivars. The results showed that 88.10 % of these primers could generate positive PCR bands in any of the three cultivars. Further phylogenetic analysis that conducted from the PCR amplification of 52 specifically sifted SSR primers within 17 cultivars approved that the possible correlation may exist between the primers and ramie fiber quality. These developed SSR markers could be applied in downstream studies, like genetic and physical maps, quantitative trait loci mapping, genetic diversity studies and cultivar fingerprinting, and breeding processes of ramie with better fiber quality under further confirmation of the correlation with ramie fiber quality.

Large-Scale SNP Discovery and Genotyping for Constructing a High-Density Genetic Map of Tea Plant Using Specific-Locus Amplified Fragment Sequencing (SLAF-seq)

Genetic maps are important tools in plant genomics and breeding. The present study reports the large-scale discovery of single nucleotide polymorphisms (SNPs) for genetic map construction in tea plant. We developed a total of 6,042 valid SNP markers using specific-locus amplified fragment sequencing (SLAF-seq), and subsequently mapped them into the previous framework map. The final map contained 6,448 molecular markers, distributing on fifteen linkage groups corresponding to the number of tea plant chromosomes. The total map length was 3,965 cM, with an average inter-locus distance of 1.0 cM. This map is the first SNP-based reference map of tea plant, as well as the most saturated one developed to date. The SNP markers and map resources generated in this study provide a wealth of genetic information that can serve as a foundation for downstream genetic analyses, such as the fine mapping of quantitative trait loci (QTL), map-based cloning, marker-assisted selection, and anchoring of scaffolds to facilitate the process of whole genome sequencing projects for tea plant.

Optimal assembly strategies of transcriptome related to ploidies of eukaryotic organisms

BACKGROUND

Several de novo transcriptome assemblers have been developed recently to assemble the short reads generated from the next-generation sequencing platforms and different strategies were employed for assembling transcriptomes of various eukaryotes without genome sequences. Though there are some comparisons among these de novo assembly tools for assembling transcriptomes of different eukaryotic organisms, there is no report about the relationship between assembly strategies and ploidies of the organisms.

RESULTS

When we de novo assembled transcriptomes of sweet potato (hexaploid), Trametes gallica (a diploid fungus), Oryza meyeriana (a diploid wild rice), five assemblers, including Edena, Oases, Soaptrans, IDBA-tran and Trinity, were used in different strategies (Single-Assembler Single-Parameter, SASP; Single-Assembler Multiple-Parameters, SAMP; Combined De novo Transcriptome Assembly, CDTA, that is multiple assembler multiple parameter). It was found that CDTA strategy has the best performance compared with other two strategies for assembling transcriptome of the hexaploid sweet potato, whereas SAMP strategy with assembler Oases is better than other strategies for assembling transcriptomes of diploid fungus and the wild rice transcriptomes.

CONCLUSION

Based on the results from ours and others, it is suggested that CDTA strategy is better used for transcriptome assembly of polyploidy organisms and SAMP strategy of Oases is outperformed for those diploid organisms without genome sequences.

De novo assembly and characterization of the fruit transcriptome of Chinese jujube (Ziziphus jujuba Mill.) Using 454 pyrosequencing and the development of novel tri-nucleotide SSR markers

Chinese jujube (Ziziphus jujuba Mill.) is an economically important deciduous tree that has high therapeutic value and health benefits. However, a lack of sequence data and molecular markers have constrained genetic and breeding studies for better fruit quality and other traits in Chinese jujube. In this study, two combined cDNA libraries of ‘Dongzao’ fruit representing the early and late stages of fruit development were constructed and sequenced on the 454 GS FLX Titanium platform. In total, 1,124,197 reads were generated and then de novo assembled into 97,479 unigenes. A total of 52,938 unigenes were homologous to genes in the NCBI non-redundant sequence database. A total of 33,123 unigenes were assigned to one or more Gene Ontology terms, and 16,693 unigenes were classified into 319 Kyoto Encyclopedia of Genes and Genomes pathways. The results showed that the Smirnoff-Wheeler pathway was the main pathway for the biosynthesis of ascorbic acid in Chinese jujube. The number of differentially expressed genes between the two stages of fruit development was 1,764, among which 974 and 790 genes were up-regulated and down-regulated, respectively. Furthermore, 9,893 sequences were identified containing SSRs. 93 primer pairs designed from the sequences with a tri-nucleotide repeat showed successful PCR amplification and could be validated in Chinese jujube accessions and Z. mauritiana Lam and Z. acidojujuba as well, of which 71 primer pairs were polymorphic. The obtained transcriptome provides a most comprehensive resource currently available for gene discovery and the development of functional markers in Z. jujuba. The newly developed microsatellite markers could be used in applications such as genetic linkage analysis and association studies, diversity analysis, and marker-assisted selection in Chinese jujube and related species.