Next Article in Journal
Fatty Acids Composition of Vegetable Oils and Its Contribution to Dietary Energy Intake and Dependence of Cardiovascular Mortality on Dietary Intake of Fatty Acids
Next Article in Special Issue
Genome-Wide Analysis of Sorbitol Dehydrogenase (SDH) Genes and Their Differential Expression in Two Sand Pear (Pyrus pyrifolia) Fruits
Previous Article in Journal
Virulence Factors of Erwinia amylovora: A Review
Previous Article in Special Issue
The Stable Level of Glutamine synthetase 2 Plays an Important Role in Rice Growth and in Carbon-Nitrogen Metabolic Balance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

De Novo Characterization of Flower Bud Transcriptomes and the Development of EST-SSR Markers for the Endangered Tree Tapiscia sinensis

1
School of Life Sciences, Northwest University, 229 Taibai Bei Road, Xi'an 710069, China
2
School of Life Sciences, Luoyang Normal University, 71 Longmen Road, Luoyang 471022, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2015, 16(6), 12855-12870; https://doi.org/10.3390/ijms160612855
Submission received: 15 April 2015 / Revised: 12 May 2015 / Accepted: 27 May 2015 / Published: 5 June 2015
(This article belongs to the Special Issue Plant Molecular Biology)

Abstract

:
Tapiscia sinensis Oliv (Tapisciaceae) is an endangered species native to China famous for its androdioecious breeding system. However, there is a lack of genomic and transcriptome data on this species. In this study, the Tapiscia sinensis transcriptomes from two types of sex flower buds were sequenced. A total of 97,431,176 clean reads were assembled into 52,169 unigenes with an average length of 1116 bp. Through similarity comparison with known protein databases, 36,662 unigenes (70.27%) were annotated. A total of 10,002 (19.17%) unigenes were assigned to 124 pathways using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database. Additionally, 10,371 simple sequence repeats (SSRs) were identified in 8608 unigenes, with 16,317 pairs of primers designed for applications. 150 pairs of primers were chosen for further validation, and the 68 pairs (45.5%) were able to produce clear polymorphic bands. Six polymorphic SSR markers were used to Bayesian clustering analysis of 51 T. sinensis individuals. This is the first report to provide transcriptome information and to develop large-scale SSR molecular markers for T. sinensis. This study provides a valuable resource for conservation genetics and functional genomics research on T. sinensis for future work.

1. Introduction

Tapiscia sinensis (Tapisciaceae) is a rare and ancient tree species native to China. It is sporadically distributed in the Yangtze River valley and the southwestern provinces of China, as well as in the northern regions of Vietnam [1,2]. Tapiscia sinensis is a woody perennial androdioecious species with both male and hermaphrodite individuals in the population [3]. Male and hermaphroditic flowers blossom at similar times and both the male and hermaphroditic flowers produce viable pollen. The flower and fruit grow synchronously on the hermaphroditic individuals of T. sinensis [4,5]. Historically, Tapiscia was usually placed in Staphyleaceae under Sapindales, and was subsequently classified as Tapisciaceae under Huerteales in APG III; however, its system location has not yet been determined [6,7]. T. sinensis is a Quaternary glacial relict plant and could be a good model for inferring the colonization and evolutionary history of post-glacial of plant species in subtropical China [3,8]. In addition, the T. sinensis leaves are rich in flavonoid, which have anti-ulcer, antispasmodic, anti-inflammatory and lipid-lowering biological activities, and has, thus, been of great value in medicine [9]. The T. sinensis tree is beautiful, with a straight trunk and luxuriant foliage, and it can be used as a garden greening tree species. Therefore, T. sinensis is a rare species that has scientific, economic and esthetic value, with excellent germplasm resources. However, due to habitat fragmentation, naturally poor regeneration and the influences of human activity, T. sinensis is currently on the International Union for Conservation of Nature (IUCN) red list of threatened species in the vulnerable A1c category [10].
Understanding the genetic relationship and population structure of the T. sinensis natural populations is important to develop conservation strategies for this rare and interesting species. However, research on T. sinensis at the molecular level remains scarce. At present, only eleven T. sinensis simple sequence repeat (SSR) markers have been developed [11]. There is still a serious lack of genomic data on this species, which has seriously hindered developments in the fields of molecular biology and genetics.
A transcriptome includes the full range of mRNA, tRNA, rRNA, and other noncoding RNA molecules expressed by one or a group of cells, organs or tissues in a specific environment or particular developmental stage. Transcriptome sequencing has significant value in molecular marker development, functional gene discovery and differential gene expression [12,13,14]. Since 2005, next-generation sequencing (NGS) technologies, such as Illumina/Solexa, Roche/454 FLX (Flagship Genome Sequencer) and ABI/SOLiD (Supported Oligo Ligation Detetion), have become available [15]. These technologies have proven to be reliable, affordable, and time-saving methods to study non-model organisms’ genomes and transcriptomes. Furthermore, transcriptome sequencing has allowed for the convenient development of additional SSR molecular markers. Compared to traditional methods, the development of SSR markers using transcriptome sequencing techniques is affordable and efficient.
In this study, we used the Illumina sequencing technology to analyze the flower buds transcriptome of T. sinensis and to develop a set of simple sequence repeat (SSR) markers. This is the first study attempting to characterize the transcriptome of Tapiscia. The transcriptome sequences reported here will provide valuable resources for the development of molecular markers, the molecular genetics research, and novel gene discovery in T. sinensis.

2. Results and Discussion

2.1. Sequencing and de novo Assembly

Next-generation sequencing technology made the transcriptome sequencing lower cost and higher throughput than traditional Sanger sequencing. As one representative of NGS, the Illumina HiSeq™ 2000 sequencing platform has been widely used for transcriptome sequencing in many non-model organisms, including Liriodendron chinense, Cucurbita moschata and Houttuynia cordata [16,17,18].
In this project, we obtained 97,431,176 clean reads, with a mean length of 100 base pairs (bp). The percentage of Q20, N, and GC was 98.31%, 0.00%, and 45.59%, respectively. The clean reads were de novo assembled using Trinity into 55,711 contigs, with a mean length of 1133 bp and an N50 of 1768 bp. The total number of unigenes was 52,169, with a mean length of 1116 bp and an N50 of 1768 bp. Of all of the unigenes, the majority (30,046 unigenes) were in the range of 201 to 1000 bp, which accounted for 57.59%. 13,723 of the unigenes (26.3%) ranged from 1001 to 2000 bp in length, while 8400 of the unigenes (16.1%) had lengths of more than 2000 bp (Figure 1).
The average 1116 bp length of the unigenes was longer than that reported for pumpkin transcriptomes (765 bp) using the Illumina HiSeq™ 2000 sequencing platform [17]. It was also greater than that of the Chinese jujube (473.4 bp) and Siberian apricot (651.6 bp) transcriptomes, using the 454 GS FLX Titanium genomic sequencer platform [19,20].
Figure 1. Length distributions of the unigenes in the assembled transcriptome.
Figure 1. Length distributions of the unigenes in the assembled transcriptome.
Ijms 16 12855 g001

2.2. Sequence Annotation and Classification

We searched the 52,169 unigenes against the protein databases. 36,662 unigenes (70.27%) were annotated, 15,507 unigenes cannot be aligned to any database remained without annotation. In the annotated unigenes, 36,608 (70.17%), 27,353 (52.43%), 14,415 (27.63%) and 10,002 (19.17%) unigenes showed homology with sequences in the Nr (Non redundant), Swiss-Prot, COG (Clusters of Orthologous Groups of proteins), and KEGG databases, respectively (Figure 2). In the annotated unigenes in the Nr database, most samples (13,776 unigenes) matched the “theobroma cacao” proteins, followed by “vitis vinifera” (8474 unigenes), “fragaria vesca subsp. Vesca” (2058 unigenes), “cucumis sativus” (1887 unigenes) and “oryza sativa Japonica Group” (1317 unigenes). In the 15,507 un-annotated unigenes, 835 unigenes were established nucleotide sequence direction (5′-3′) and amino sequence of the predicted coding region by using ESTScan [21].
Figure 2. The unigenes showed homology with sequences in the Nr, Swiss-Prot, COG (Clusters of Orthologous Groups of proteins) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases.
Figure 2. The unigenes showed homology with sequences in the Nr, Swiss-Prot, COG (Clusters of Orthologous Groups of proteins) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases.
Ijms 16 12855 g002
Based on Nr annotation, 15,690 unigenes were assigned gene ontology (GO) with 83,842 functional terms. In total, 36,497 terms were assigned to the biological processes group, 30,925 terms were assigned to the cellular components group and 16,420 terms were assigned to the molecular functions group (Figure 3). Of the biological process category, metabolic process (8526, 23.36%) and cellular process (8406, 23.03%) were the most frequent terms used. Genes related to response to stimulus (3476 unigenes), anatomical structure formation (270 unigenes), and rhythmic process (105 unigenes) were also identified in this group. In the class of “response to stimulus”, unigene0008775, unigene0005239, and unigene0049385, that were respectively annotated as “Pedicel, carpel, stamen, petal differentiation and expansion stage, group 2-like protein”, “Phytochrome and flowering time regulatory protein isoform 2”, and “Protein early flowering 3”, which putatively related to response to environmental stimulus were discovered. In the class of “rhythmic process”, unigene0044321 and unigene0017817, which were respectively annotated as “circadian clock-associated FKF1” and “Protein EARLY FLOWERING 4 OS”, were found. In the class of “anatomical structure formation”, unigene0011847 was annotated as “Sex determination protein tasselseed-2-like” was found. These unigenes may involve in flower development, differentiation and flowering time determination in T. sinensis. Regarding the cellular components category, the cell and cell parts were the most prominent terms. In the molecular functions category, catalytic activity (7861 unigenes) and binding (6905 unigenes) were the most highly represented terms.
Figure 3. Unigenes were assigned to Gene Ontology.
Figure 3. Unigenes were assigned to Gene Ontology.
Ijms 16 12855 g003
Furthermore, all unigenes were subjected to a search against the Clusters of Orthologous Groups of proteins (COG) database for functional prediction and classification, and 27,672 sequences were assigned to COG classifications (Figure 4). Of the 25 COG categories, the cluster for “general function prediction only” (4902, 17.71%) was the largest group, followed by “transcription” (2704, 9.77%), “replication, recombination and repair” (2451, 8.86%), “signal transduction mechanisms” (1991, 7.19%) and “posttranslational modification, protein turnover, chaperones” (1955, 7.06%). Additionally, 1605 unigenes were assigned to “carbohydrate transport and metabolism” group, which include unigene0027247 that was annotated as “Glycogen synthase”, unigene0027034 annotated as “6-phosphofructokinase”, and unigene0026072 annotated as “Glyceraldehyde-3-phosphate dehydrogenase/erythrose-4-phosphate dehydrogenase”; 743 unigenes were assigned to “energy production and conversion”, which include unigene0032350 annotated as “Fumarase”, unigene0032531 annotated as “Glycerol-3-phosphate dehydrogenase”, unigene0038751 annotated as “Nicotinamide-adenine dinucleotid dehydrogenase, Flavin adenine dinucleotide containing subunit”, and unigene0039174 annotated as “Phosphoenolpyruvate carboxykinase (ATP)”. These key genes laid the foundation for the study of T. sinensis.
Figure 4. Unigenes were assigned to COG classifications.
Figure 4. Unigenes were assigned to COG classifications.
Ijms 16 12855 g004
To identify the biological pathways in T. sinensis, the annotated unigenes were mapped to the reference of typical pathways in the KEGG database. In this study, 10,002 (19.17%) unigenes were assigned to 124 pathways, which include pathways about “Starch and sucrose metabolism” (224 unigenes), “Terpenoid backbone biosynthesis” (72 unigenes), and “Carbon fixation in photosynthetic organisms” (136 unigenes). The top five pathways in KEGG analysis consist of “Metabolic pathways” (ko01100), “Biosynthesis of secondary metabolites” (ko01110), “Ribosomes” (ko03010), “Spliceosome” (ko03040), and “RNA transport” (ko03013) (Table S1). It is worth noting that 260 (2.6%) unigenes in the ‘‘Plant hormone signal transduction’’ pathway were found in the pathways of auxin, cytokinine, gibberellins, abscisic acid, ethylene, brassinosteroid, jasmonic acid and salicylic acid, they may be involved in the regulation of “Cell enlargement”, “Plant growth”, “Shoot initiation”, and “Senescence” in T. sinensis flowers (Figure S1).

2.3. Genes Putatively Related to Flower Development

In the natural population of T. sinensis, male and hermaphrodite individuals coexist. The flowers of male individuals have normal stamens and degenerated pistil, and the hermaphrodite flowers have normal pistil and stamens (Figure 5). In this study we identified several unigenes that are putative homologs of Arabidopsis thaliana, Populus tomentosa, Theobroma cacao and Vitis vinifera genes that are involved in floral development. These include genes that flowering time determination, such as the MADS box genes SOC1 and Flowering Locus C, and perceive or respond to environmental signals, such as FKF1 and ELF4, and Sex determination protein tasselseed-2 (Table S2).
Figure 5. The male flower and hermaphrodite flower of T. sinensis, (A) and (C) are flowers from male individual; (B) and (D) are bisexual flowers from hermaphrodite individual. Anatomical images (C,D) are the difference between male and hermaphrodite flowers about pistils. Picture (C) shows degenerated pistil, picture (D) shows well-developed pistil.
Figure 5. The male flower and hermaphrodite flower of T. sinensis, (A) and (C) are flowers from male individual; (B) and (D) are bisexual flowers from hermaphrodite individual. Anatomical images (C,D) are the difference between male and hermaphrodite flowers about pistils. Picture (C) shows degenerated pistil, picture (D) shows well-developed pistil.
Ijms 16 12855 g005
Flowering is an important developmental process for the reproduction of the plant. It is regulated by environmental and genetic factors. The MADS box genes SOC1 and Flowering Locus C (FLC) of the plant Arabidopsis thaliana have been shown to have an important role in the integration of molecular flowering time pathways [22,23]. These genes are essential for the correct timing of flowering, and help to ensure that fertilization occurs at the time of maximal reproductive potential. The transition of vegetative to flowering in plants is coordinated by external cues such as day length (photoperiod) or temperature. Day-length sensing involves an interaction between the relative length of day and night, and endogenous rhythms that are controlled by the plant circadian clock. The genes “EARLY FLOWERING 4 (ELF4)” and “FLAVIN-BINDING, KELCH REPEAT, F-BOX (FKF1)” are involved in circadian regulation and photoperiod perception [24,25]. The unigenes were annotated as sex determination protein tasselseed-2 (Ts-2) has been shown to be involved in the arrest of developing pistils in male flowers of the grass family (Poaceae), Zea mays and Tripsacum dactyloides [26,27]. Maize through programmed abortion of preformed organ primordia produces separate unisexual flowers.
By comparing the unigenes that first assembled separately in male and hermaphrodite samples, we found that the expression of the above-mentioned unigenes had no significant difference in the two types of flowers, except unigene0052119 and unigene0050634. The unigene0052119 was annotated as “ELF4-like protein” over expression in flowers of hermaphrodites, which may involve in circadian regulation and photoperiod perception in T. sinensis. The unigene0050634 was annotated as “sex determination protein tasselseed-2-like” over expression in flowers of males, which may involve in the arrest of developing pistils in male flowers. These unigenes are important resources for the study of floral development and flower organ formation of T. sinensis in the future.

2.4. SSR Mining from the T. sinensis Transcriptome

SSRs are useful as molecular markers for genetics and biology researches. All 52,169 unigenes generated from the transcriptome were used to develop potential SSR loci using MISA (MIcroSAtellite identification) software. A total of 10,371 SSRs were identified in 8608 unigenes, with 5439 sequences suit to designed primers, and 16,317 pairs of primers designed for application. Of all of the SSR containing sequences, 1436 sequences contained more than one SSR, and 733 SSRs were present in compound forms (Table 1).
Table 1. Summary of the EST-SSRs identified in the transcriptome sequences.
Table 1. Summary of the EST-SSRs identified in the transcriptome sequences.
ItemNumber
Total number of sequences examined52,169
Total size of examined sequences (bp)58,253,075
Total number of identified SSRs10,371
Number of SSR containing sequences8608
Number of sequences containing more than 1 SSR1436
Number of SSRs present in compound formation733
Of the 10,371 identified SSRs, the most abundant repeat motif types were di-nucleotide repeat motifs (6888, 66.42%), followed by tri-nucleotide (2486, 23.97%), tetra-nucleotide (522, 5.03%), hexa-nucleotide (296, 2.85%) and penta-nucleotide (179, 1.73%) repeat motifs. The frequencies of the SSRs with different numbers of tandem repeats are summarized in Table 2. SSRs with six tandem repeats were the most common (2797, 26.97%), followed by seven tandem repeats (1831, 17.66%), five tandem repeats (1447, 13.95%), and nine tandem repeats (1256, 12.11%). The most common type of SSR motif was AG/CT which accounted for 54.5%, followed by AT/AT and AAG/CTT (6.5%), AC/GT (5.2%), ACC/GGT (4.1%) (Figure 6). This result is consistent with previous research showing di-nucleotide repeats to be the most abundant type, followed by were tri-nucleotide repeats [17,28,29]. The most abundant di-nucleotide motifs AG/CT and tri-nucleotide motifs AAG/CTT were also coincident with earlier research [17].
In the 8608 SSR-containing ESTs, the most of them are functional. These ESTs are useful for the SSR identification of related to flower types and some other important phenotypes, such as stress resistance. The detail and functional description of 8608 SSR-containing EST accessions were shown in Tables S3 and S4.
Figure 6. Frequency distribution of SSRs based on the motif types.
Figure 6. Frequency distribution of SSRs based on the motif types.
Ijms 16 12855 g006
Table 2. The distribution of the identified EST-SSRs in sequences using the MISA (MIcroSAtellite identification) software.
Table 2. The distribution of the identified EST-SSRs in sequences using the MISA (MIcroSAtellite identification) software.
Repeat NumbersMotifTotalRepeat Number (%)
Di-Tri-Tetra-Penta-Hexa
4003681582667927.64
5013191011512144713.95
6202670352511279726.97
71416411004183117.66
8116449112121711.73
912560000125612.11
1082310008247.95
1119510001961.89
˃1182001110.11
Total6888248652217929610,371100
Motif (%)66.4223.975.031.732.85

2.5. SSR Marker Assessment and Application

A suit of 150 pairs of primers was synthetized and assessed in the 24 T. sinensis individuals. The result shows that, 117 pairs (78%) successfully amplified target fragments for SSR-PCR (Simple Sequence Repeat Polymerase Chain Reaction), including 68 pairs that were able to produce clear polymorphic bands (45.5%) and 49 pairs of primers that were able to produce single bands without polymorphism. The results of the SSR-PCR amplification are shown as examples in Figure 7. Of 68 polymorphic SSR markers (N = 24), the detected number of alleles per locus ranged from 2 to 6 with an average of 3.24 alleles per locus, and the observed and expected heterozygosity varied from 0.000 to 1.000 and from 0.042 to 0.782, respectively. Polymorphic information content varied 0.040 to 0.729. Details of the 68 polymorphic SSR markers from T. sinensis are shown in Table S5. In addition, the Bayesian clustering analysis of the 51 individuals was constructed based on six polymorphic SSR markers, with the characteristics shown in Table 3. The analysis indicated that we developed SSR markers can differentiate population genetic differences of different geographic locations, dividing the 51 individuals into two main groups of 29 (fromWuling Mountains) and 22 (from Qinling Mountains) (Figure S2).
In recent years, EST-SSR has been widely used in the studies of plant at the DNA level. In 2014, Guo et al. used forty effective EST-SSR to test the genetic diversity of 15 hexaploid oat germplasm resources, and clustering them into three groups according to species. The forty primers showed higher availability on the study of 31 chromosome ploidy unknown oat germplasm resources, and some diploids and tetraploids of oat germplasm resources were found in them [30]. Wang et al. used 28 EST-SSR derived from alfalfa and soybean to study the genetic relationship of 26 accessions representing four species of Crotalaria, and re determined the phylogenetic classification of the 26 accessions [31]. Recently, nineteen EST-SSR analyses provide insights about genetic relatedness, population structure and gene flow in grass pea of 176 accessions [32]. In this study, we developed EST-SSR could be used for genetic diversity analysis, identification of Germplasm Resources and the molecular phylogeny and evolution in Tapisciaceae.
Figure 7. The results of the SSR-PCR that produced clear bands. M refers to DNA markers pBR322 DNA/MspI. Numbers 1 to 24 represent twenty-four individuals of T. sinensis. The products of primer TS061 are single bands without polymorphism, while those of TS062 are clear polymorphic bands.
Figure 7. The results of the SSR-PCR that produced clear bands. M refers to DNA markers pBR322 DNA/MspI. Numbers 1 to 24 represent twenty-four individuals of T. sinensis. The products of primer TS061 are single bands without polymorphism, while those of TS062 are clear polymorphic bands.
Ijms 16 12855 g007
Table 3. Characteristics of the six primer pairs used for the phylogenetic analysis of the 51 individuals.
Table 3. Characteristics of the six primer pairs used for the phylogenetic analysis of the 51 individuals.
PrimerMotifKHOHEPICHW
TS053(TCA)690.7060.8350.804ND
TS126(TAC)760.7060.7680.723NS
TS103(AAGA)550.5490.6310.573NS
TS149(TTG)1160.6080.7690.725NS
TS060(TGT)660.5290.5010.464NS
TS013(AGA)850.510.6530.592NS
K, No. (Number) of alleles; HO, Observed heterozygosity; HE, Expected heterozygosity; PIC, Polymorphism Information Content; HW, Hardy-Weinberg equilibrium test; NS, Not significant; ND, Not done.
For comparison with existing databases, functional annotation and SSR identification using PCR amplification, the results indicated that the quality of the T. sinensis transcriptome sequencing and assembly was very high. To date, this is the first research on T. sinensis transcriptome sequencing and analysis. These data will lay the foundations for the future research on functional gene identification, molecular genetics, physiology and biochemistry of T. sinensis.

3. Materials and Methods

3.1. Plant Material and RNA Extraction

The inflorescences of three male and three hermaphrodite individuals were collected separately at Northwest University, Shanxi, China (34°15′N, 108°55′E, H 400 m) in April 2014. The samples were between 0.2–1.0 cm long and inserted a lot of small flower buds. The total RNA of above listed six samples were isolated separately using the Trizol Kit (Promega, Beijing, China) following by the manufacturer’s instructions, then were treated with RNase-free DNase I (Takara Biotech Incorporation, Otsu, Japan) for 30 min at 37 °C to remove residual DNA. RNA quality was verified using a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and were also checked by RNase free agarose gel electrophoresis. The equal amounts of RNA from the three same sex individuals were mixed together. Then, the two total RNA samples, one from three male individuals and another from three hermaphroditic individuals, were prepared to following experiment.
The subsequent steps, which include two cDNA libraries construction, Illumina sequencing, and first assembled of the two samples were carried out separately.

3.2. cDNA Library Construction and Illumina Sequencing

In brief, Poly (A) mRNA was enriched by Oligo (dT) magnetic beads. All mRNA was cut into short fragments by adding fragmentation buffer. Taking these short fragments as templates, random hexamer primers were used to synthesize first-strand cDNA. Second-strand cDNA was synthesized using the first-strand cDNA as templates in the presence of buffer, dNTPs, RNase H and DNA polymerase I. The short fragments were purified using the QIAquick PCR extraction kit (Qiagen, Shanghai, China) and resolved with EB buffer for end reparation and A tailing. The short fragments were then connected with sequencing adapters. Following agarose gel electrophoresis and extraction of cDNA from gels, the cDNA fragments were purified and enriched by PCR to construct the final cDNA library for sequencing. Illumina sequencing was performed at Gene Denovo Co. Guangzhou, China using the HiSeq™ 2000 platform. The sequencing data is available from the the National Center for Biotechnology Information (NCBI) Short Read Archive (SRA, Bethesda, Maryland, MD, USA. Accession numbers: PRJNA284864).

3.3. De Novo Assembly and Annotation

The raw sequencing reads were first filtered by removing invalid reads, including reads with adaptor contamination, with ambiguous “N” bases at a ratio greater than 5% and reads with more than 50% base with quality lower than 20 in one sequence. The short clean reads were assembled using Trinity software to construct contigs [33]. These contigs were used for the further process of sequence clustering with the software TGICL (TIGR Gene Indices clustering tools) to form unigenes [34]. The reads from the male flower buds and hermaphrodite flower buds were first assembled separately, then assembled together to acquire non-redundant unigenes that were as long as possible.
All of the unigenes were aligned to sequences in the four protein databases, such as Nr (NCBI non-redundant database), Swiss-Prot, COG and KEGG (the Kyoto Encyclopedia of Genes and Genomes), using blastx with an E-value threshold <1 × 10−5. The sequence direction of the unigenes was determined according to the best alignment results. When the results were conflicted among databases, the direction was determined consecutively by Nr, Swiss-Prot, KEGG and COG. When a unigene would not align to any database, ESTScan was used to predict coding regions and determine sequence direction [21]. GO annotations for the unigenes were achieved using the Blast2GO software (Valencia, Spain) [35]. Functional classification of the unigenes was performed using WEGO (web tool for plotting GO annotations) software (Hangzhou, China) [36].

3.4. SSR Detection and Primer Design

The 52,169 unigenes were searched for Simple Sequence Repeats (SSRs) using the MISA program (MIcroSAtellite identification tool). The parameters were based on the minimum number of repeats, which were set as follows: six repeat units for di-nucleotides, five for tri-nucleotides, and four for tetra-, penta- and hexa-nucelotides. The maximum interruption distance between two SSRs was specified as 100 bases. Primer pairs of each SSR loci were designed using the Primer3 program in the flanking regions of SSRs [37]. Primers were designed based on the following criteria: (1) GC content ranged between 40% and 60%; (2) Primer length between 18 and 27 bp; (3) Melting temperature between 55 and 63 °C; and (4) The expected PCR product sizes ranged from 100 to 280 bp.

3.5. SSR Validation and Assessment

To validate the EST-SSR markers developed from the transcriptome, fifty-one T. sinensis adult individuals, of which 22 from the Qinling Mountains (QL) (33°31′N, 108°35′E, H 1348 m) and 29 from the Wuling Mountains (WL) of China (28°47′N, 110°15′E, H 490 m), were collected.
A total of 150 primer pairs were synthesized and twenty-four of the fifty-one samples were selected randomly to validate the polymorphism of the SSR markers. After detecting Hardy-Weinberg equilibrium and linkage disequilibrium between microsatellites by using Cervus 3.0 and Genepop [38,39], six polymorphic SSR markers were selected for Bayesian clustering analysis on the 51 T. sinensis individuals using Structure v2.3 software (Stanford, CA, USA). [40]. The analysis was ran with the parameters as: 30,000 Burn-in Period, 100,000 MCMC Reps (Markov Chain Monte Carlo Repetitions), three runs were conducted for each K (1–8). STRUCTURE HARVESTER were used to calculate ΔK values and identified the suitable value of K = 2 for the population structure [41].

4. Conclusions

In this study, the transcriptome of T. sinensis was analyzed using high-throughput Illumina sequencing. Through de novo assembly and sequence annotation, we obtained 52,169 sequences and identified some important characteristics of the T. sinensis transcriptome. In addition, we identified and characterized 10,371 potential SSR markers. A set of 150 SSR markers were selected to detect polymorphisms among the 24 T. sinensis individuals. A total of 117 pairs (78%) had successfully amplified target fragments, with 68 pairs showing polymorphic loci. This study is the first application of Illumina sequencing technology used to investigate the transcriptome and assemble RNA-seq reads without a reference genome for T. sinensis. The SSRs and sequences will provide useful resources for conservation genetics and functional genomics research on T. sinensis in the future.

Supplementary Files

Supplementary File 1

Acknowledgments

This work was supported by the National Science Foundation of China (31270428), National Foundation for Talents Cultivation in Fundamental Science of China (J1210063), Training Project of Innovative Talents for Graduate Students of Northwest University (YZZ13064) and the Youth Research Foundation of Luoyang Normal University (2013-QNJJ-004).

Author Contributions

Wen-Zhe Liu and Xiao-Jun Zhou conceived and designed the experiments. Xiao-Jun Zhou wrote the paper. Yue-Yue Wang, Ya-Nan Xu and Rong-Shan Yan performed the experiments. Peng Zhao analyzed the data. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bremer, B.; Bremer, K.; Chase, M.; Fay, M.; Reveal, J.; Soltis, D.; Soltis, P.; Stevens, P. An update of the angiosperm phylogeny group classification for the orders and families of flowering plants: APG III. Bot. J. Linn. Soc. 2009, 161, 105–121. [Google Scholar]
  2. Di, W.Z.; Yu, Z. The First Countries to Protect Rare and Endangered Plants in Shanxi Province; Northwest University Press: Xi’an, China, 1989; pp. 200–203. [Google Scholar]
  3. Li, D.Z.; Cai, J.; Wen, J. Flora of China; Science Press: Beijing, China, 1996; p. 496. [Google Scholar]
  4. Liu, W.Z.; Kang, H.Q.; Zheng, H.C.; Feng, Y.Z. An investigation on the sexual reproductive cycle in Tapiscia sinensis. J. Syst. Evol. 2008, 46, 175–182. [Google Scholar]
  5. Lü, W.; Liu, W.Z. Pollination biology in androdioecious species Tapiscia sinensis. Chin. Bull. Bot. 2010, 45, 713–722. (In Chinese) [Google Scholar]
  6. Krause, J. Staphyleaceae. In Die natürlichen Pflanzenfamilien, 2nd ed.; Duncker & Humblot 20b: Berlin, Germany, 1942; pp. 255–321. [Google Scholar]
  7. Matthews, M.L.; Endress, P.K. Comparative floral structure and systematics in Crossosomatales (Crossosomataceae, Stachyuraceae, Staphyleaceae, Aphloiaceae, Geissolomataceae, Ixerbaceae, Strasburgeriaceae). Bot. J. Linn. Soc. 2005, 147, 1–46. [Google Scholar] [CrossRef]
  8. Tang, Y.; Li, L. On historical elements and tertiary precursors of eastern Asiatic angiospermous flora-based on phytogeography of staphyleaceae, morinaceae and caprifoliaceae. Acta. Phytotaxon. Sin. 1996, 34, 453–478. (In Chinese) [Google Scholar]
  9. Xie, C.P. A review of research advances in rare and endangered plant Tapiscia sinensis. Subtrop. Plant Sci. 2006, 35, 71–74. (In Chinese) [Google Scholar]
  10. Tapiscia sinensis. The IUCN red list of threatened species. Available Online: http://www.iucnredlist.org. (accessed on 14 April 2015).
  11. Zhang, P.F.; Wang, Z.G.; Liu, Y.L.; Tian, H.; Yao, X.H.; Zhang, J.J. Development and characterization of 11 polymorphic microsatellite markers in Tapiscia sinensis (Staphyleaceae). Appl. Plant Sci. 2014, 1, 1300051. [Google Scholar]
  12. Zhang, J.; Liang, S.; Duan, J.; Wang, J.; Chen, S.; Cheng, Z.; Zhang, Q.; Liang, X.; Li, Y. De novo assembly and characterisation of the transcriptome during seed development, and generation of genic-SSR markers in peanut (Arachis hypogaea L.). BMC Genomics 2012, 13, 90. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, Z.; Fang, B.; Chen, J.; Zhang, X.; Luo, Z.; Huang, L.; Chen, X.; Li, Y. De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweetpotato (Ipomoea batatas). BMC Genomics 2010, 11, 726. [Google Scholar] [CrossRef] [PubMed]
  14. Lu, T.; Lu, G.; Fan, D.; Zhu, C.; Li, W.; Zhao, Q.; Feng, Q.; Zhao, Y.; Guo, Y.; Li, W.; et al. Function annotation of the rice transcriptome at singlenucleotide resolution by RNA-seq. Genome Res. 2010, 20, 1238–1249. [Google Scholar] [CrossRef] [PubMed]
  15. Fu, N.; Wang, Q.; Shen, H.L. De novo assembly, gene annotation and marker development using Illumina paired-end transcriptome sequences in celery (Apium graveolens L.). PLoS ONE 2013, 8, e57686. [Google Scholar] [CrossRef] [PubMed]
  16. Yang, Y.; Xu, M.; Luo, Q.; Wang, J.; Li, H. De novo transcriptome analysis of Liriodendron chinense petals and leaves by Illumina sequencing. Gene 2014, 534, 155–162. [Google Scholar] [CrossRef] [PubMed]
  17. Wu, T.; Luo, S.; Wang, R.; Zhong, Y.; Xu, X.; Lin, Y.E.; He, X.; Sun, B.; Huang, H. The first Illumina-based de novo transcriptome sequencing and analysis of pumpkin (Cucurbita moschata Duch.) and SSR marker development. Mol. Breed. 2014, 34, 1437–1447. [Google Scholar] [CrossRef]
  18. Wei, L.; Li, S.; Liu, S.; He, A.; Wang, D.; Wang, J.; Tang, Y.; Wu, X. Transcriptome analysis of Houttuynia cordata Thunb. by illumina paired-end RNA sequencing and SSR marker discovery. PLoS ONE 2014, 9, e84105. [Google Scholar] [CrossRef] [PubMed]
  19. Li, Y.; Xu, C.; Lin, X.; Cui, B.; Wu, R.; Pang, X. 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. PLoS ONE 2014, 9, e106438. [Google Scholar] [CrossRef] [PubMed]
  20. Dong, S.; Liu, Y.; Niu, J.; Ning, Y.; Lin, S.; Zhang, Z. De novo transcriptome analysis of the Siberian apricot (Prunus sibirica L.) and search for potential SSR markers by 454 pyrosequencing. Gene 2014, 544, 220–227. [Google Scholar] [CrossRef] [PubMed]
  21. Iseli, C.; Jongeneel, C.V.; Bucher, P. ESTScan: A program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. ISMB 1999, 99, 138–148. [Google Scholar]
  22. Onouchi, H.; Igeño, M.I.; Périlleux, C.; Graves, K.; Coupland, G. Mutagenesis of plants overexpressing CONSTANS demonstrates novel interactions among arabidopsis flowering-time genes. Plant Cell 2000, 12, 885–900. [Google Scholar] [CrossRef] [PubMed]
  23. Michaels, S.D.; Amasino, R.M. FLOWERING LOCUS C encodes a novel MADS domain protein that acts as a repressor of flowering. Plant Cell 1999, 11, 949–956. [Google Scholar] [CrossRef] [PubMed]
  24. Doyle, M.R.; Davis, S.J.; Bastow, R.M.; McWatters, H.G.; Kozma-Bognár, L.; Nagy, F.; Millar, A.J.; Amasino, R.M. The ELF4 gene controls circadian rhythms and flowering time in Arabidopsis thaliana. Nature 2002, 419, 74–77. [Google Scholar] [CrossRef] [PubMed]
  25. Li, F.; Zhang, X.; Hu, R.; Wu, F.; Ma, J.; Meng, Y.; Fu, Y. Identification and molecular characterization of FKF1 and GI homologous genes in soybean. PLoS ONE 2013, 8, e79036. [Google Scholar] [CrossRef] [PubMed]
  26. Kinney, M.S.; Columbus, J.T.; Friar, E.A. Molecular evolution of the maize sex-determining gene TASSELSEED2 in Bouteloua (Poaceae). Mol. Phylogenet. Evol. 2003, 29, 519–528. [Google Scholar] [CrossRef] [PubMed]
  27. DeLong, A.; Calderon-Urrea, A.; Dellaporta, S.L. Sex determination gene TASSELSEED2 of maize encodes a short-chain alcohol dehydrogenase required for stage-specific floral organ abortion. Cell 1993, 74, 757–768. [Google Scholar] [CrossRef]
  28. Li, D.; Deng, Z.; Qin, B.; Liu, X.; Men, Z. De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree (Hevea brasiliensis Muell. Arg.). BMC Genomics 2012, 13, 192. [Google Scholar] [CrossRef] [PubMed]
  29. Triwitayakorn, K.; Chatkulkawin, P.; Kanjanawattanawong, S.; Sraphet, S.; Yoocha, T.; Sangsrakru, D.; Chanprasert, J.; Ngamphiw, C.; Jomchai, N.; Therawattanasuk, K.; et al. Transcriptome sequencing of Hevea brasiliensis for development of microsatellite markers and construction of a genetic linkage map. DNA Res. 2011, 18, 471–482. [Google Scholar] [CrossRef] [PubMed]
  30. Guo, H.Y.; Jia, J.Q.; Zhang, X.H.; La, H.X.; Yang, W.D. Development and application of EST-SSR markers in oat. Plant Sci. J. 2014, 32, 240–250. [Google Scholar]
  31. Wang, M.L.; Mosjidis, J.A.; Morris, J.B.; Dean, R.E.; Jenkins, T.M.; Pederson, G.A. Genetic diversity of Crotalaria germplasm assessed through phylogenetic analysis of EST-SSR markers. Genome 2006, 49, 707–715. [Google Scholar] [CrossRef] [PubMed]
  32. Soren, K.R.; Yadav, A.; Pandey, G.; Gangwar, P.; Parihar, A.K.; Bohra, A.; Dixit, G.P.; Datta, S.; Singh, N.P. EST-SSR analysis provides insights about genetic relatedness, population structure and gene flow in grass pea (Lathyrus sativus). Plant Breed. 2015. [Google Scholar] [CrossRef]
  33. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  34. Pertea, G.; Huang, X.; Liang, F.; Antonescu, V.; Sultana, R.; Karamycheva, S.; Lee, Y.; White, J.; Cheung, F.; Parvizi, B.; et al. TIGR Gene Indices clustering tools (TGICL): A software system for fast clustering of large EST datasets. Bioinformatics 2003, 19, 651–652. [Google Scholar] [CrossRef] [PubMed]
  35. Conesa, A.; Götz, S.; García-Gómez, J.M.; Terol, J.; Talón, M.; Robles, M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 2005, 21, 3674–3676. [Google Scholar] [CrossRef] [PubMed]
  36. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Zhang, Z.; Wang, J.; Li, S.; Li, R.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, 293–297. [Google Scholar] [CrossRef] [PubMed]
  37. Rozen, S.; Skaletsky, H. Primer3 on the WWW for general users and for biologist programmers. In Bioinformatics Methods and Protocols; Humana Press: New Jersey, NJ, USA, 1999; pp. 365–386. [Google Scholar]
  38. Kalinowski, S.T.; Taper, M.L.; Marshall, T.C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol. Ecol. 2007, 16, 1099–1106. [Google Scholar] [CrossRef] [PubMed]
  39. Rousset, F. Genepop’ 007: A complete re-implementation of the Genepop software for Windows and Linux. Mol. Ecol. Resour. 2008, 8, 103–106. [Google Scholar] [CrossRef] [PubMed]
  40. Documentation for Structure Software. Available online: http://pritchardlab.stanford.edu/structure.html (accessed on 2 February 2010).
  41. Earl, D.A.; von Holdt, B.M. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Resour. 2012, 4, 359–361. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Zhou, X.-J.; Wang, Y.-Y.; Xu, Y.-N.; Yan, R.-S.; Zhao, P.; Liu, W.-Z. De Novo Characterization of Flower Bud Transcriptomes and the Development of EST-SSR Markers for the Endangered Tree Tapiscia sinensis. Int. J. Mol. Sci. 2015, 16, 12855-12870. https://doi.org/10.3390/ijms160612855

AMA Style

Zhou X-J, Wang Y-Y, Xu Y-N, Yan R-S, Zhao P, Liu W-Z. De Novo Characterization of Flower Bud Transcriptomes and the Development of EST-SSR Markers for the Endangered Tree Tapiscia sinensis. International Journal of Molecular Sciences. 2015; 16(6):12855-12870. https://doi.org/10.3390/ijms160612855

Chicago/Turabian Style

Zhou, Xiao-Jun, Yue-Yue Wang, Ya-Nan Xu, Rong-Shan Yan, Peng Zhao, and Wen-Zhe Liu. 2015. "De Novo Characterization of Flower Bud Transcriptomes and the Development of EST-SSR Markers for the Endangered Tree Tapiscia sinensis" International Journal of Molecular Sciences 16, no. 6: 12855-12870. https://doi.org/10.3390/ijms160612855

APA Style

Zhou, X. -J., Wang, Y. -Y., Xu, Y. -N., Yan, R. -S., Zhao, P., & Liu, W. -Z. (2015). De Novo Characterization of Flower Bud Transcriptomes and the Development of EST-SSR Markers for the Endangered Tree Tapiscia sinensis. International Journal of Molecular Sciences, 16(6), 12855-12870. https://doi.org/10.3390/ijms160612855

Article Metrics

Back to TopTop