Identification, characterization and utilization of unigene derived microsatellite markers in tea (Camellia sinensis L.)

Evaluation and characteristic analysis of SSRs from the transcriptomic sequences of Perilla crop (Perilla frutescens L.)

Perilla crop is a self-fertilizing annual plant, cultivated and used mainly in East Asia. Perilla frutescens var. frutescens seeds are rich in unsaturated fatty acids, which have health benefits, and Perilla frutescens var. crispa leaves are rich in anthocyanins. However, genomic analysis such as whole genome sequencing or genetic mapping has not been performed on Perilla crop. This current study confirms the abundance and diversity of 15,991 simple sequence repeats (SSRs) classified in previous studies in the Perilla genome, selects and designs 1,538 SSR primer sets, and confirms which SSR primer sets exhibit high polymorphism. Of the 15,991 SSRs classified, there were 9,910 (62%) dinucleotide repeats, 5,652 (35.3%) trinucleotide repeats, and 429 (2.7%) tetranucleotide repeats. Among these, the most identified was (CT)n with a total of 4,817. The 15,991 SSRs had 4 to 26 repeats. Four repeats were the most frequent with 11,084 (69.3%). A total of 1,538 SSR primers were selected and designed to confirm polymorphism, of which 157 showed persistent and clear polymorphism. Among these 157 SSR primer sets, 98 (62.4%) were dinucleotide repeats, 39 (24.8%) were trinucleotide repeats, and 20 (12.7%) were tetranucleotide repeats. Among 549 SSR primers that showed polymorphism, trinucleotide repeats showed persistent polymorphism at a high rate. Therefore, when developing SSR primer sets for Perilla crop in the future, it is recommended that trinucleotide repeats be selected first. These research results will be helpful in future genomic analysis and development of SSR primers in Perilla crop.

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.

Conservation of an Agrobacterium cT-DNA insert in Camellia section Thea reveals the ancient origin of tea plants from a genetically modified ancestor

Introduction

Many higher plants contain cellular T-DNA (cT-DNA) sequences from Agrobacterium and have been called "natural genetically modified organisms" (nGMOs). Among these natural transformants, the tea plant Camellia sinensis var. sinensis cv. Shuchazao contains a single 5.5 kb T-DNA fragment (CaTA) with three inactive T-DNA genes, with a 1 kb inverted repeat at the ends. Camellia plants are allogamous, so that each individual may contain two different CaTA alleles.

Methods

142 Camellia accessions, belonging to 10 of 11 species of the section Thea, were investigated for the presence of CaTA alleles.

Results discussion

All accessions were found to contain the CaTA insert, showing that section Thea derives from a single transformed ancestor. Allele phasing showed that 82 accessions each contained two different CaTA alleles, 60 others had a unique allele. A phylogenetic tree of these 225 alleles showed two separate groups, A and B, further divided into subgroups. Indel distribution corresponded in most cases with these groups. The alleles of the different Camellia species were distributed over groups A and B, and different species showed very similar CaTA alleles. This indicates that the species boundaries for section Thea may not be precise and require revision. The nucleotide divergence of the indirect CaTA repeats indicates that the cT-DNA insertion took place about 15 Mio years ago, before the emergence of section Thea. The CaTA structure of a C. fangchengensis accession has an exceptional structure. We present a working model for the origin and evolution of nGMO plants derived from allogamous transformants.

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.

Comparative transcriptomic analysis and antioxidant defense mechanisms in clusterbean (Cyamopsis tetragonoloba (L.) Taub.) genotypes with contrasting drought tolerance

To understand drought tolerance mechanism(s) in clusterbean (Cyamopsis tetragonoloba), we conducted physiological, biochemical, and de novo comparative transcriptome analysis of drought-tolerant (RGC-1002) and drought-sensitive (RGC-1066) genotypes subjected to 30 days of drought stress. Relative water content (RWC) was maintained in tolerant genotype but was reduced in sensitive genotype. Leaf pigment concentrations were higher in tolerant genotype. Net photosynthesis was significantly decreased in sensitive genotype but insignificant reduction was found in tolerant genotype. Enzymatic antioxidant (GR, APX, DHAR) activities were enhanced in tolerant genotype, while there were insignificant changes in these enzymes in sensitive genotype. The ratios of antioxidant molecules (ASC/DHA and GSH/GSSG) were higher in tolerant genotype as compared to sensitive genotype. In sensitive genotype, 6625 differentially expressed genes (DEGs) were upregulated and 5365 genes were downregulated. In tolerant genotype, 5206 genes were upregulated and 2793 genes were downregulated. In tolerant genotype, transketolase family protein, phosphoenolpyruvate carboxylase 3, temperature-induced lipocalin, and cytochrome oxidase were highly upregulated. Moreover, according to Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the drought tolerance may be attributed to upregulated starch and sucrose metabolism-related genes in tolerant genotype. Finally, quantitative real-time PCR confirmed the reproducibility of the RNA-seq data.

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 a new set of genic SSR markers in the genus Gentiana: in silico mining, characterization and validation

Gentiana is an important genus of around 360 medicinally important species, majority of which are not well characterized. Despite its importance, very few genomic resources are available for Gentiana L. Till date, the number of informative and robust simple sequence repeat (SSR)-based markers is limited and very few efforts have been made for their development. A set of robust, freely accessible and informative SSR markers for Gentiana is a pre-requisite for any molecular systematic as well as improvement studies in this group of pharmacologically valuable plants. In view of the importance of these plants, Expressed Sequence Tag (EST) sequences of 18 Gentiana species were surveyed for the development of a large set of non-redundant SSR markers. A total of 5808 perfect SSR with an average length of 17 bp and relative abundance of 214 loci/Mb were identified in the analysed 47,487 EST sequences using Krait software. Mapping of the ESTs resulted in gene ontology annotations of 49.14% of the sequences. Based on these perfect SSRs, 2902 primer pairs were designed, and 60 markers were randomly selected and validated on a set of Gentiana kurroo Royle accessions. Among the screened markers, 39 (65%) were found to be cross-species transferable. This is the first report of the largest set of functional, novel genic SSR markers in Gentiana, which will be a valuable resource for future characterization, genotype identification, conservation and genomic studies in the various species of this group of important medicinal plants.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02969-4.

Underpinning the molecular programming attributing heat stress associated thermotolerance in tea (Camellia sinensis (L.) O. Kuntze)

The most daunting issue of global climate change is the deleterious impact of extreme temperatures on tea productivity and quality, which has resulted in a quest among researchers and growers. The current study aims to unravel molecular programming underpinning thermotolerance by characterizing heat tolerance and sensitivity response in 20 tea cultivars. The significantly higher negative influence of heat stress was recorded in a sensitive cultivar with reduced water retention (47%), chlorophyll content (33.79%), oxidation potential (32.48%), and increase in membrane damage (76.4%). Transcriptional profiling of most tolerant and sensitive cultivars identified 78 differentially expressed unigenes with chaperon domains, including low and high molecular weight heat shock protein (HSP) and heat shock transcription factors (HSFs) involved in heat shock response (HSR). Further, predicted transcriptional interactome network revealed their key role in thermotolerance via well-co-ordinated transcriptional regulation of aquaporins, starch metabolism, chlorophyll biosynthesis, calcium, and ethylene mediated plant signaling system. The study identified the key role of HSPs (CsHSP90) in regulating HSR in tea, wherein, structure-based molecular docking revealed the inhibitory role of geldanamycin (GDA) on CsHSP90 by blocking ATP binding site at N-terminal domain of predicted structure. Subsequently, GDA mediated leaf disc inhibitor assay further affirmed enhanced HSR with higher expression of CsHSP17.6, CsHSP70, HSP101, and CsHSFA2 genes in tea. Through the current study, efforts were made to extrapolate a deeper understanding of chaperons mediated regulation of HSR attributing thermotolerance in tea.

Genome-Wide Analysis of Heat Shock Transcription Factors in Ziziphus jujuba Identifies Potential Candidates for Crop Improvement Under Abiotic Stress

Plant heat shock transcription factors (Hsfs) play a significant role in adoption under abiotic stress conditions by modulating the expression of several stress-responsive genes. Analysis of the Hsf gene family will serve to understand the molecular mechanism which is involved in response to abiotic stress. The Ziziphus species grows in warm and dry regions and is inherently tolerant to abiotic stress conditions; thus, Ziziphus is a highly enriched source of genes conferring abiotic stress tolerance. Therefore, the present study provides a comprehensive genome-wide analysis of the Hsf gene family in Z. jujuba. Identified 21 non-redundant Hsf genes were grouped into three major classes (classes A, B, and C) based on the phylogenetic analysis. Promoter and gene ontology analysis suggested that ZjHsfs perform diverse functions in response to abiotic stress conditions. Two paralogous pairs resulting from tandem gene duplication events were identified. Also, physio-chemical properties of chromosomal locations, gene structure, motifs, and protein domain organization of Hsfs were analyzed. Real-time PCR expression analyses revealed that most of the Z. jujuba Hsf genes are differentially expressed in response to heat stress. The analysis suggested ZjHsf-2, ZjHsf-3, ZjHsf-5, ZjHsf-7, ZjHsf-8, ZjHsf-10, ZjHsf-12, ZjHsf-17, and ZjHsf-18 were the outstanding candidate genes for imparting heat stress tolerance and for future functional analysis. The present analysis laid the foundation for understanding the molecular mechanism of the Hsf gene family regulating Z. jujuba development and tolerance to abiotic stress conditions.

Transcriptome wide SSR discovery cross-taxa transferability and development of marker database for studying genetic diversity population structure of Lilium species

Lily belongs to family liliaceae, which mainly propagates vegetatively. Therefore, sufficient number of polymorphic, informative, and functional molecular markers are essential for studying a wide range of genetic parameters in Lilium species. We attempted to develop, characterize and design SSR (simple sequence repeat) markers using online genetic resources for analyzing genetic diversity and population structure of Lilium species. We found di-nucleotide repeat motif were more frequent (4684) within 0.14 gb (giga bases) transcriptome than other repeats, of which was two times higher than tetra-repeat motifs. Frequency of di-(AG/CT), tri-(AGG/CTT), tetra-(AAAT), penta-(AGAGG), and hexa-(AGAGGG) repeats was 34.9%, 7.0%, 0.4%, 0.3%, and 0.2%, respectively. A total of 3607 non-redundant SSR primer pairs was designed based on the sequences of CDS, 5'-UTR and 3'-UTR region covering 34%, 14%, 23%, respectively. Among them, a sub set of primers (245 SSR) was validated using polymerase chain reaction (PCR) amplification, of which 167 primers gave expected PCR amplicon and 101 primers showed polymorphism. Each locus contained 2 to 12 alleles on average 0.82 PIC (polymorphic information content) value. A total of 87 lily accessions was subjected to genetic diversity analysis using polymorphic SSRs and found to separate into seven groups with 0.73 to 0.79 heterozygosity. Our data on large scale SSR based genetic diversity and population structure analysis may help to accelerate the breeding programs of lily through utilizing different genomes, understanding genetics and characterizing germplasm with efficient manner.

Development of EST-SSR markers and association mapping with floral traits in Syringa oblata

BACKGROUND

Lilac (Syringa oblata) is an important woody plant with high ornamental value. However, very limited genetic marker resources are currently available, and little is known about the genetic architecture of important ornamental traits for S. oblata, which is hindering its genetic studies. Therefore, it is of great significance to develop effective molecular markers and understand the genetic architecture of complex floral traits for the genetic research of S. oblata.

RESULTS

In this study, a total of 10,988 SSRs were obtained from 9864 unigene sequences with an average of one SSR per 8.13 kb, of which di-nucleotide repeats were the dominant type (32.86%, 3611). A set of 2042 primer pairs were validated, out of which 932 (45.7%) exhibited successful amplifications, and 248 (12.1%) were polymorphic in eight S. oblata individuals. In addition, 30 polymorphic EST-SSR markers were further used to assess the genetic diversity and the population structure of 192 cultivated S. oblata individuals. Two hundred thirty-four alleles were detected, and the PIC values ranged from 0.23 to 0.88 with an average of 0.51, indicating a high level of genetic diversity within this cultivated population. The analysis of population structure showed two major subgroups in the association population. Finally, 20 significant associations were identified involving 17 markers with nine floral traits using the mixed linear model. Moreover, marker SO104, SO695 and SO790 had significant relationship with more than one trait.

CONCLUSION

The results showed newly developed markers were valuable resource and provided powerful tools for genetic breeding of lilac. Beyond that, our study could serve an efficient foundation for further facilitate genetic improvement of floral traits for lilac.

Development of novel g-SSR markers in guava (Psidium guajava L.) cv. Allahabad Safeda and their application in genetic diversity, population structure and cross species transferability studies

Dearth of genomic resources particularly, microsatellite markers in nutritionally and commercially important fruit crop, guava necessitate the development of the novel genomic SSR markers through the library enrichment techniques. Three types of 3' -biotinylated oligonucleotide probes [(CT)14, (GT)12, and (AAC)8] were used to develop microsatellite enriched libraries. A total of 153 transformed colonies were screened of which 111 positive colonies were subjected for Sanger sequencing. The clones having more than five motif repeats were selected for primer designing and a total of 38 novel genomic simple sequence repeats could be identified. The g-SSRs had the motif groups ranging from monomer to pentamer out of which dimer group occurred the most (89.47%). Out of 38 g-SSRs markers developed, 26 were found polymorphic, which showed substantial genetic diversity among the guava genotypes including wild species. The average number of alleles per locus, major allele frequency, gene diversity, expected heterozygosity and polymorphic information content of 26 SSRs were 3.46, 0.56, 0.53, 0.29 and 0.46, respectively. The rate of cross-species transferability of the developed g-SSR loci varied from 38.46 to 80.77% among the studied wild Psidium species. Generation of N-J tree based on 26 SSRs grouped the 40 guava genotypes into six clades with two out-groups, the wild guava species showed genetic distinctness from cultivated genotypes. Furthermore, population structure analysis grouped the guava genotypes into three genetic groups, which were partly supported by PCoA and N-J tree. Further, AMOVA and PCoA deciphered high genetic diversity among the present set of guava genotypes including wild species. Thus, the developed novel g-SSRs were found efficient and informative for diversity and population structure analyses of the guava genotypes. These developed novel g-SSR loci would add to the new genomic resource in guava, which may be utilized in genomic-assisted guava breeding.

Cross-transferability of SSR markers developed in Rhododendron species of Himalaya

Rhododendron is a genus of evergreen woody ornamental plants of northern hemisphere with strong cold resistance, attractive flowers and high altitude adaptation capacity. The genus originated and diversified from Sino-Himalayan region and spread across the world, and has high species diversity in Northeast India. To assess cross-species amplification, we tested 32 microsatellites markers in fifteen taxa of the genus Rhododendron of North-eastern Himalaya, of which fourteen microsatellites were newly developed from Rhododendron simsii, and eighteen microsatellites were previously developed from Rhododendron catawbiense and Rhododendron mucronatum var. ripense. Nine pairs of primers were amplified successfully in all species, however, none of them was failed for amplification in any of the species. The average observed heterozygosity, expected heterozygosity and PIC value were recorded as 0.310, 0.433 and 0.379 respectively. Clustering based on neighbour-joining analysis revealed the potential of these markers to segregate species according to their subgenus level, however, subspecies exhibited closeness with each other. Cross-application of these microsatellite loci will provide a potentially useful tool to investigate the genetic diversity, population structure, gene flow, phylogenetics and evolutionary relationships in species of genus Rhododendron.

Novel intron length polymorphic (ILP) markers from starch biosynthesis genes reveal genetic relationships in Indian wheat varieties and related species

Introns experience lesser selection pressure, thus are liable for higher polymorphism. Intron Length Polymorphic (ILP) markers designed from exon-flanking introns exploits this polymorphic potential and have been proved to be a robust co-dominant marker in eukaryotes. Wheat is among the most consumed cereal crop by majority of the word population. It is a rich source of calories in the form of stored starch. In the current study, starch biosynthesis genes were mined for development of ILP markers and their subsequent utilization for genetic characterization of popular Indian wheat varieties and transferability to wild relatives. Sixty-one markers generated 122 alleles and showed 77-88.5% transferability (mean PIC: 0.36) to the related species. A subset of markers showed clear genetic distinctions (Avg. genetic dissimilarity = 0.42) among Indian wheat varieties, signifying the importance of novel ILPs. 'Kenphad25' showed maximum genetic dissimilarity with 'K 8962' (0.82), while maximum genetic similarity was observed between 'Safed Lerma' and 'RAJ 4037' (0.1). This is the first report of ILP markers in wheat and will be a useful genomic resource for future germplasm conservation and molecular breeding studies.

Combating Climate Change in the Kenyan Tea Industry

Climate change triggered by global warming poses a major threat to agricultural systems globally. This phenomenon is characterized by emergence of pests and diseases, extreme weather events, such as prolonged drought, high intensity rains, hailstones and frosts, which are becoming more frequent ultimately impacting negatively to agricultural production including rain-fed tea cultivation. Kenya is predominantly an agricultural based economy, with the tea sector generating about 26% of the total export earnings and about 4% gross domestic product (GDP). In the recent years, however, the country has witnessed unstable trends in tea production associated with climate driven stresses. Toward mitigation and adaptation of climate change, multiple approaches for impact assessment, intensity prediction and adaptation have been advanced in the Kenyan tea sub-sector. Further, pressure on tea breeders to release improved climate-compatible cultivars for the rapidly deteriorating environment has resulted in the adoption of a multi-targeted approach seeking to understand the complex molecular regulatory networks associated with biotic and abiotic stresses adaptation and tolerance in tea. Genetic modeling, a powerful tool that assists in breeding process, has also been adopted for selection of tea cultivars for optimal performance under varying climatic conditions. A range of physiological and biochemical responses known to counteract the effects of environmental stresses in most plants that include lowering the rates of cellular growth and net photosynthesis, stomatal closure, and the accumulation of organic solutes such as sugar alcohols, or osmolytes have been used to support breeding programs through screening of new tea cultivars suitable for changing environment. This review describes simulation models combined with high resolution climate change scenarios required to quantify the relative importance of climate change on tea production. In addition, both biodiversity and ecosystem based approaches are described as a part of an overall adaptation strategy to mitigate adverse effects of climate change on tea in Kenya and gaps highlighted for urgent investigations.

Development of polymorphic EST-SSR markers and their applicability in genetic diversity evaluation in Rhododendron arboreum

The genus Rhododendron, known for large impressive flowers is widely distributed throughout the world. Rhododendrons have limited genetic information, despite of comprising high species diversity, morphological overlap and weak genetic barrier. In present study, expressed sequence tag (EST) data from Rhododendron catawbiense Michx (Subgenus Hymenanthes, Section Ponticum) and Rhododendron mucronatum var. ripense (Makino) E.H. Wilson (Subgenus Tsutsusi, Section Tsutsusi) were utilized for mining and identification of the SSRs for genetic diversity analysis of R. arboreum Smith (Subgenus Tsutsusi, Section Tsutsusi). A total of 249 SSRs were developed from 1767 contigs. Di-nucleotide was found to be most abundant repeat followed by tri- and tetra-nucleotide repeats. The motif AG/CT was most common di-nucleotide motif (31.73%), whereas, AAC/GTT (8.43%), ACG/CGT (8.03%), AAG/CTT (7.23%) and AGG/CCT (6.43%) were most abundant tri-nucleotide repeat motif. Among these SSRs, 168 sequences were only fit into the criteria to design flanking primer pairs. A total of 30 randomly selected primer pairs were utilized for validation and genetic diversity study in 36 genotypes of R. arboreum collected from western Himalayan region. In aggregate, 26 SSR markers (86.66%) produced good and repeatable amplifications. Expected heterozygosity (H_E) ranged from 0.322 to 0.841 and observed heterozygosity (H_O) ranged from 0.327 to 1.000 and PIC value ranged from 0.008 to 0.786. These primers were able to distinguish the geographic differences of occurrence based on cluster analysis. These developed EST-SSRs can be useful in future population genetics analysis and micro-evolutionary studies in Rhododendron species.

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

Development of polymorphic simple sequence repeat markers in Huperzia serrata (Lycopodiaceae)

PREMISE

The natural population size of Huperzia serrata (Lycopodiaceae) has dramatically decreased and the species has become endangered due to overexploitation. Here, we developed simple sequence repeat (SSR) markers for H. serrata to survey both its genetic diversity and population structure.

METHODS AND RESULTS

Based on 177 individuals, 120 SSR primer pairs were developed and optimized from five regions of the H. serrata transcriptomic data. Of these primer pairs, 20 were successfully amplified and 10 showed obvious polymorphism. These polymorphic loci were investigated to study the genetic diversity of H. serrata. Two to 11 alleles per locus were identified, the level of observed heterozygosity ranged from 0.00 to 1.00, and the level of expected heterozygosity ranged from 0.19 to 0.79. All loci were successfully amplified in H. crispata, H. sutchueniana, and H. selago.

CONCLUSIONS

The 10 polymorphic primer pairs developed here will be valuable for studies of the endangered H. serrata and other related species.

Circular RNA architecture and differentiation during leaf bud to young leaf development in tea (Camellia sinensis)

Circular RNA (circRNA) discovery, expression patterns and experimental validation in developing tea leaves indicates its correlation with circRNA-parental genes and potential roles in ceRNA interaction network. Circular RNAs (circRNAs) have recently emerged as a novel class of abundant endogenous stable RNAs produced by circularization with regulatory potential. However, identification of circRNAs in plants, especially in non-model plants with large genomes, is challenging. In this study, we undertook a systematic identification of circRNAs from different stage tissues of tea plant (Camellia sinensis) leaf development using rRNA-depleted circular RNA-seq. By combining two state-of-the-art detecting tools, we characterized 3174 circRNAs, of which 342 were shared by each approach, and thus considered high-confidence circRNAs. A few predicted circRNAs were randomly chosen, and 20 out of 24 were experimental confirmed by PCR and Sanger sequencing. Similar in other plants, tissue-specific expression was also observed for many C. sinensis circRNAs. In addition, we found that circRNA abundances were positively correlated with the mRNA transcript abundances of their parental genes. qRT-PCR validated the differential expression patterns of circRNAs between leaf bud and young leaf, which also indicated the low expression abundance of circRNAs compared to the standard mRNAs from the parental genes. We predicted the circRNA-microRNA interaction networks, and 54 of the differentially expressed circRNAs were found to have potential tea plant miRNA binding sites. The gene sets encoding circRNAs were significantly enriched in chloroplasts related GO terms and photosynthesis/metabolites biosynthesis related KEGG pathways, suggesting the candidate roles of circRNAs in photosynthetic machinery and metabolites biosynthesis during leaf development.

Transcriptome sequence analysis and mining of SSRs in Jhar Ber (Ziziphus nummularia (Burm.f.) Wight & Arn) under drought stress

Ziziphus nummularia (Burm.f.) Wight & Arn., a perennial shrub that thrives in the arid regions, is naturally tolerant to drought. However, there are limited studies on the genomics of drought tolerance in Ziziphus sp. In this study, RNA-sequencing of one month old seedlings treated with PEG 6000 was performed using Roche GS-FLX454 Titanium pyrosequencing. A total of 367,176 raw sequence reads were generated, and upon adapter trimming and quality filtration 351,872 reads were assembled de novo into 32,739 unigenes. Further characterization of the unigenes indicated that 73.25% had significant hits in the protein database. Kyoto encyclopedia of genes and genomes database (KEGG) identified 113 metabolic pathways from the obtained unigenes. A large number of drought-responsive genes were obtained and among them differential gene expression of 16 highly induced genes was validated by qRT-PCR analysis. To develop genic-markers, 3,425 simple sequence repeats (SSRs) were identified in 2,813 unigene sequences. The data generated shall serve as an important reservoir for the identification and characterization of drought stress responsive genes for development of drought tolerant crops.

Cross-species transferability of EST-SSR markers developed from the transcriptome of Melilotus and their application to population genetics research

Melilotus is one of the most important legume forages, but the lack of molecular markers has limited the development and utilization of Melilotus germplasm resources. In the present study, 151 M clean reads were generated from various genotypes of Melilotus albus using Illumina sequencing. A total of 19,263 potential EST-SSRs were identified from 104,358 unigene sequences. Moreover, 18,182 primer pairs were successfully designed, and 550 primer pairs were selected using criteria of base repeat type, fragment length and annealing temperature. In addition, 550 primer pairs were screened by using PCR amplification products and used to assess polymorphisms in 15 M. albus accessions. A total of 114 primer pairs were detected as being highly polymorphic, and the average polymorphism information content (PIC) value was 0.79. Furthermore, those 114 polymorphic primer pairs were used to evaluate the transferability to 18 species of the genus Melilotus, and 70 EST-SSR markers were found to be transferable among the 18 Melilotus species. According to the UPGMA dendrogram and STRUCTURE analysis, the 18 Melilotus species were classified into three clusters. This study offers a valuable resource for the genetic diversity and molecular assisted breeding of germplasm resources in the genus Melilotus.

SSR-Based DNA Fingerprinting and Diversity Assessment Among Indian Germplasm of Euryale ferox: an Aquatic Underutilized and Neglected Food Crop

Euryale ferox is native to Southeast Asia and China, and it is one of the important aquatic food crops propagated mostly in eastern part of India. The aim of the present study was to characterize and evaluate the genetic diversity of ex situ collections of E. ferox germplasm from different geographical states of India using microsatellite (simple sequence repeats (SSRs)) markers. Ten SSR markers were analyzed to assess DNA fingerprinting and genetic diversity of 16 cultivated germplasm of E. ferox. Total 37 polymorphic alleles were recorded with an average of 3.7 allele frequency per primer. The polymorphic information content value varied from 0.204 to 0.735 with mean of 0.448. A high range of heterozygosity (Ho 0.228; He 0.512) was detected in the present study. The neighbor-joining (N-J) tree and the principle coordinate analysis showed that the germplasm divided in to three main clusters. The results of the present investigation comply that SSR markers are effective for computing genetic assessment of genetic diversity and similarity with classifying cultivated varieties of E. ferox. Evaluation of genetic diversity among Indian E. ferox germplasm could provide useful information for genetic improvement.

Proteomic analysis of tea plants (Camellia sinensis) with purple young shoots during leaf development

Tea products made from purple leaves are highly preferred by consumers due to the health benefits. This study developed a proteome reference map related to color changes during leaf growth in tea (Camellia sinensis) plant with purple young shoots using two-dimensional electrophoresis (2-DE). Forty-six differentially expressed proteins were detected in the gel and successfully identified by using MALDI-TOF/TOF-MS. The pronounced changes in the proteomic profile between tender purple leaves (TPL) and mature green leaves (MGL) included: 1) the lower activity of proteins associated with CO2 assimilation, energy metabolism and photo flux efficiency and higher content of anthocyanins in TPL than those in MGL may protect tender leaves against photo-damage; 2) the higher abundance of chalcone synthase (CHS), chalcone isomerase (CHI) and flavonol synthase (FLS) likely contributes to the synthesis of anthocyanins, catechins and flavonols in TPL tissues; 3) higher abundance of stress response proteins, such as glutathione S-transferases (GST) and phospholipid hydroperoxide glutathione peroxidase (PHGPx), could enhance the tolerance of TPL tissues to adverse condition in; and 4) the increased abundance of proteins related to protein synthesis, nucleic acids and cell wall proteins should be beneficial for the proliferation and expansion of leaf cell in TPL tissues. qPCR analysis showed that the expression of differentially abundant proteins was regulated at the transcriptional level. Therefore, the results indicated that higher abundance of CHI and CHS may account for the production of the purple-shoot phenotype in Wuyiqizhong 18 and thereby, enhancing the anthocyanin biosynthesis. The higher abundance of glutamine synthetase (GS) proteins related to the theanine biosynthesis may improve the flavor of tea products from TPL materials. Thus, this work should help to understand the molecular mechanisms underlying the changes in leaf color alteration.

Development of genomic tools in a widespread tropical tree, Symphonia globulifera L.f.: a new low-coverage draft genome, SNP and SSR markers

Population genetic studies in tropical plants are often challenging because of limited information on taxonomy, phylogenetic relationships and distribution ranges, scarce genomic information and logistic challenges in sampling. We describe a strategy to develop robust and widely applicable genetic markers based on a modest development of genomic resources in the ancient tropical tree species Symphonia globulifera L.f. (Clusiaceae), a keystone species in African and Neotropical rainforests. We provide the first low-coverage (11X) fragmented draft genome sequenced on an individual from Cameroon, covering 1.027 Gbp or 67.5% of the estimated genome size. Annotation of 565 scaffolds (7.57 Mbp) resulted in the prediction of 1046 putative genes (231 of them containing a complete open reading frame) and 1523 exact simple sequence repeats (SSRs, microsatellites). Aligning a published transcriptome of a French Guiana population against this draft genome produced 923 high-quality single nucleotide polymorphisms. We also preselected genic SSRs in silico that were conserved and polymorphic across a wide geographical range, thus reducing marker development tests on rare DNA samples. Of 23 SSRs tested, 19 amplified and 18 were successfully genotyped in four S. globulifera populations from South America (Brazil and French Guiana) and Africa (Cameroon and São Tomé island, F_ST  = 0.34). Most loci showed only population-specific deviations from Hardy-Weinberg proportions, pointing to local population effects (e.g. null alleles). The described genomic resources are valuable for evolutionary studies in Symphonia and for comparative studies in plants. The methods are especially interesting for widespread tropical or endangered taxa with limited DNA availability.

Transcriptome Analysis Reveals Candidate Genes involved in Blister Blight defense in Tea (Camellia sinensis (L) Kuntze)

To unravel the molecular mechanism of defense against blister blight (BB) disease caused by an obligate biotrophic fungus, Exobasidium vexans, transcriptome of BB interaction with resistance and susceptible tea genotypes was analysed through RNA-seq using Illumina GAIIx at four different stages during ~20-day disease cycle. Approximately 69 million high quality reads were assembled de novo, yielding 37,790 unique transcripts with more than 55% being functionally annotated. Differentially expressed, 149 defense related transcripts/genes, namely defense related enzymes, resistance genes, multidrug resistant transporters, transcription factors, retrotransposons, metacaspases and chaperons were observed in RG, suggesting their role in defending against BB. Being present in the major hub, putative master regulators among these candidates were identified from predetermined protein-protein interaction network of Arabidopsis thaliana. Further, confirmation of abundant expression of well-known RPM1, RPS2 and RPP13 in quantitative Real Time PCR indicates salicylic acid and jasmonic acid, possibly induce synthesis of antimicrobial compounds, required to overcome the virulence of E. vexans. Compendiously, the current study provides a comprehensive gene expression and insights into the molecular mechanism of tea defense against BB to serve as a resource for unravelling the possible regulatory mechanism of immunity against various biotic stresses in tea and other crops.

Assessment of Functional EST-SSR Markers (Sugarcane) in Cross-Species Transferability, Genetic Diversity among Poaceae Plants, and Bulk Segregation Analysis

Expressed sequence tags (ESTs) are important resource for gene discovery, gene expression and its regulation, molecular marker development, and comparative genomics. We procured 10000 ESTs and analyzed 267 EST-SSRs markers through computational approach. The average density was one SSR/10.45 kb or 6.4% frequency, wherein trinucleotide repeats (66.74%) were the most abundant followed by di- (26.10%), tetra- (4.67%), penta- (1.5%), and hexanucleotide (1.2%) repeats. Functional annotations were done and after-effect newly developed 63 EST-SSRs were used for cross transferability, genetic diversity, and bulk segregation analysis (BSA). Out of 63 EST-SSRs, 42 markers were identified owing to their expansion genetics across 20 different plants which amplified 519 alleles at 180 loci with an average of 2.88 alleles/locus and the polymorphic information content (PIC) ranged from 0.51 to 0.93 with an average of 0.83. The cross transferability ranged from 25% for wheat to 97.22% for Schlerostachya, with an average of 55.86%, and genetic relationships were established based on diversification among them. Moreover, 10 EST-SSRs were recognized as important markers between bulks of pooled DNA of sugarcane cultivars through BSA. This study highlights the employability of the markers in transferability, genetic diversity in grass species, and distinguished sugarcane bulks.

Genetic Divergence between Camellia sinensis and Its Wild Relatives Revealed via Genome-Wide SNPs from RAD Sequencing

Tea is one of the most popular beverages across the world and is made exclusively from cultivars of Camellia sinensis. Many wild relatives of the genus Camellia that are closely related to C. sinensis are native to Southwest China. In this study, we first identified the distinct genetic divergence between C. sinensis and its wild relatives and provided a glimpse into the artificial selection of tea plants at a genome-wide level by analyzing 15,444 genomic SNPs that were identified from 18 cultivated and wild tea accessions using a high-throughput genome-wide restriction site-associated DNA sequencing (RAD-Seq) approach. Six distinct clusters were detected by phylogeny inferrence and principal component and genetic structural analyses, and these clusters corresponded to six Camellia species/varieties. Genetic divergence apparently indicated that C. taliensis var. bangwei is a semi-wild or transient landrace occupying a phylogenetic position between those wild and cultivated tea plants. Cultivated accessions exhibited greater heterozygosity than wild accessions, with the exception of C. taliensis var. bangwei. Thirteen genes with non-synonymous SNPs exhibited strong selective signals that were suggestive of putative artificial selective footprints for tea plants during domestication. The genome-wide SNPs provide a fundamental data resource for assessing genetic relationships, characterizing complex traits, comparing heterozygosity and analyzing putatitve artificial selection in tea plants.

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.

Genome-Wide Survey and Analysis of Microsatellite Sequences in Bovid Species

Microsatellites or simple sequence repeats (SSRs) have become the most popular source of genetic markers, which are ubiquitously distributed in many eukaryotic and prokaryotic genomes. This is the first study examining and comparing SSRs in completely sequenced genomes of the Bovidae. We analyzed and compared the number of SSRs, relative abundance, relative density, guanine-cytosine (GC) content and proportion of SSRs in six taxonomically different bovid species: Bos taurus, Bubalus bubalis, Bos mutus, Ovis aries, Capra hircus, and Pantholops hodgsonii. Our analysis revealed that, based on our search criteria, the total number of perfect SSRs found ranged from 663,079 to 806,907 and covered from 0.44% to 0.48% of the bovid genomes. Relative abundance and density of SSRs in these Bovinae genomes were non-significantly correlated with genome size (Pearson, r < 0.420, p > 0.05). Perfect mononucleotide SSRs were the most abundant, followed by the pattern: perfect di- > tri- > penta- > tetra- > hexanucleotide SSRs. Generally, the number of SSRs, relative abundance, and relative density of SSRs decreased as the motif repeat length increased in each species of Bovidae. The most GC-content was in trinucleotide SSRs and the least was in the mononucleotide SSRs in the six bovid genomes. The GC-contents of tri- and pentanucleotide SSRs showed a great deal of similarity among different chromosomes of B. taurus, O. aries, and C. hircus. SSR number of all chromosomes in the B. taurus, O.aries, and C. hircus is closely positively correlated with chromosome sequence size (Pearson, r > 0.980, p < 0.01) and significantly negatively correlated with GC-content (Pearson, r < -0.638, p < 0.01). Relative abundance and density of SSRs in all chromosomes of the three species were significantly negatively correlated with GC-content (Pearson, r < -0.333, P < 0.05) but not significantly correlated with chromosome sequence size (Pearson, r < -0.185, P > 0.05). Relative abundances of the same nucleotide SSR type showed great similarity among different chromosomes of B. taurus, O. aries, and C. hircus.

Genome-Wide Computational Analysis of Musa Microsatellites: Classification, Cross-Taxon Transferability, Functional Annotation, Association with Transposons & miRNAs, and Genetic Marker Potential

The development of organized, informative, robust, user-friendly, and freely accessible molecular markers is imperative to the Musa marker assisted breeding program. Although several hundred SSR markers have already been developed, the number of informative, robust, and freely accessible Musa markers remains inadequate for some breeding applications. In view of this issue, we surveyed SSRs in four different data sets, developed large-scale non-redundant highly informative therapeutic SSR markers, and classified them according to their attributes, as well as analyzed their cross-taxon transferability and utility for the genetic study of Musa and its relatives. A high SSR frequency (177 per Mbp) was found in the Musa genome. AT-rich dinucleotide repeats are predominant, and trinucleotide repeats are the most abundant in transcribed regions. A significant number of Musa SSRs are associated with pre-miRNAs, and 83% of these SSRs are promising candidates for the development of therapeutic SSR markers. Overall, 74% of the SSR markers were polymorphic, and 94% were transferable to at least one Musa spp. Two hundred forty-three markers generated a total of 1047 alleles, with 2-8 alleles each and an average of 4.38 alleles per locus. The PIC values ranged from 0.31 to 0.89 and averaged 0.71. We report the largest set of non-redundant, polymorphic, new SSR markers to be developed in Musa. These additional markers could be a valuable resource for marker-assisted breeding, genetic diversity and genomic studies of Musa and related species.

Development and Characterization of Transcription Factor Gene-Derived Microsatellite (TFGM) Markers in Medicago truncatula and Their Transferability in Leguminous and Non-Leguminous Species

Transcription factors (TFs) are critical adaptor molecules that regulate many plant processes by controlling gene expression. The recent increase in the availability of TF data has made TFs a valuable resource for genic functional microsatellite marker development. In the present study, we developed TF gene-derived microsatellite (TFGM) markers for Medicago truncatula and assessed their cross-species transferability. A total of 203 SSRs were identified from 1467 M. truncatula TF coding sequences, 87.68% of which were trinucleotide repeats, followed by mono- (4.93%) and hexanucleotide repeats (1.48%). Further, 142 TFGM markers showed a high level of transferability to the leguminous (55.63%-85.21%) and non-leguminous (28.17%-50.00%) species. Polymorphisms of 27 TFGM markers were evaluated in 44 alfalfa accessions. The allele number per marker ranged from two to eight with an average of 4.41, and the PIC values ranged from 0.08 to 0.84 with an average of 0.60. Considering the high polymorphism, these TFGM markers developed in our study will be valuable for genetic relationship assessments, marker-assisted selection and comparative genomic studies in leguminous and non-leguminous species.

Identification, validation and cross-species transferability of novel Lavandula EST-SSRs

We identified and characterized EST-SSRs with strong discrimination power against Lavandula angustifolia and Lavandula x intermedia . The markers also showed considerable cross-species transferability rate into six related Lavandula species. Lavenders (Lavandula) are important economical crops grown around the globe for essential oil production. In an attempt to develop genetic markers for these plants, we analyzed over 13,000 unigenes developed from L. angustifolia and L. x intermedia EST databases, and identified 3,459 simple sequence repeats (SSR), which were dominated by trinucleotides (41.2 %) and dinucleotides (31.45 %). Approximately, 19 % of the unigenes contained at least one SSR marker, over 60 % of which were localized in the UTRs. Only 252 EST-SSRs were 18 bp or longer from which 31 loci were validated, and 24 amplified discrete fragments with 85 % polymorphism in L. x intermedia and L. angustifolia. The average number of alleles in L. x intermedia and L. angustifolia were 3.42 and 3.71 per marker with average PIC values of 0.47 and 0.52, respectively. These values suggest a moderate to strong level of informativeness for the markers, with some loci producing unique fingerprints. The cross-species transferability rate of the markers ranges 50-100 % across eight species. The utility of these markers was assessed in eight Lavandula species and 15 L. angustifolia and L. x intermedia cultivars, and the dendrogram deduced from their similarity indexes successfully delineated the species into their respective sections and the cultivars into their respective species. These markers have potential for application in fingerprinting, diversity studies and marker-assisted breeding of Lavandula.

Identification and characterization of microsatellites in expressed sequence tags and their cross transferability in different plants

Expressed sequence tags (EST) are potential source for the development of genic microsatellite markers, gene discovery, comparative genomics, and other genomic studies. In the present study, 7630 ESTs were examined from NCBI for SSR identification and characterization. A total of 263 SSRs were identified with an average density of one SSR/4.2 kb (3.4% frequency). Analysis revealed that trinucleotide repeats (47.52%) were most abundant followed by tetranucleotide (19.77%), dinucleotide (19.01%), pentanucleotide (9.12%), and hexanucleotide repeats (4.56%). Functional annotation was done through homology search and gene ontology, and 35 EST-SSRs were selected. Primer pairs were designed for evaluation of cross transferability and polymorphism among 11 plants belonging to five different families. Total 402 alleles were generated at 155 loci with an average of 2.6 alleles/locus and the polymorphic information content (PIC) ranged from 0.15 to 0.92 with an average of 0.75. The cross transferability ranged from 34.84% to 98.06% in different plants, with an average of 67.86%. Thus, the validation study of annotated 35 EST-SSR markers which correspond to particular metabolic activity revealed polymorphism and evolutionary nature in different families of Angiospermic plants.

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.

Transcriptome analysis of the oil-rich tea plant, Camellia oleifera, reveals candidate genes related to lipid metabolism

Background

Rapidly driven by the need for developing sustainable sources of nutritionally important fatty acids and the rising concerns about environmental impacts after using fossil oil, oil-plants have received increasing awareness nowadays. As an important oil-rich plant in China, Camellia oleifera has played a vital role in providing nutritional applications, biofuel productions and chemical feedstocks. However, the lack of C. oleifera genome sequences and little genetic information have largely hampered the urgent needs for efficient utilization of the abundant germplasms towards modern breeding efforts of this woody oil-plant.

Results

Here, using the 454 GS-FLX sequencing platform, we generated approximately 600,000 RNA-Seq reads from four tissues of C. oleifera. These reads were trimmed and assembled into 104,842 non-redundant putative transcripts with a total length of ∼38.9 Mb, representing more than 218-fold of all the C. oleifera sequences currently deposited in the GenBank (as of March 2014). Based on the BLAST similarity searches, nearly 42.6% transcripts could be annotated with known genes, conserved domains, or Gene Ontology (GO) terms. Comparisons with the cultivated tea tree, C. sinensis, identified 3,022 pairs of orthologs, of which 211 exhibited the evidence under positive selection. Pathway analysis detected the majority of genes potentially related to lipid metabolism. Evolutionary analysis of omega-6 fatty acid desaturase (FAD2) genes among 20 oil-plants unexpectedly suggests that a parallel evolution may occur between C. oleifera and Olea oleifera. Additionally, more than 2,300 simple sequence repeats (SSRs) and 20,200 single-nucleotide polymorphisms (SNPs) were detected in the C. oleifera transcriptome.

Conclusions

The generated transcriptome represents a considerable increase in the number of sequences deposited in the public databases, providing an unprecedented opportunity to discover all related-genes associated with lipid metabolic pathway in C. oleifera. It will greatly enhance the generation of new varieties of C. oleifera with increased yields and high quality.

Comparison of hematopoietic cancer stem cells with normal stem cells leads to discovery of novel differentially expressed SSRs

Tandem repeat expansion in the transcriptomics level has been considered as one of the underlying causes of different cancers. Cancer stem cells are a small portion of cancer cells within the main neoplasm and can remain alive during chemotherapy and re-induce tumor growth. The EST-SSR background of cancer stem cells and possible roles of expressed SSRs in altering normal stem cells to cancer ones have not been investigated yet. Here, SSR distributions in hematopoietic normal and cancer stem cells were compared based on the expressed EST-SSR. One hundred eighty nine and 223 EST-SSRs were identified in cancer and normal stem cells, respectively. The EST-SSR expression pattern was significantly different between normal and cancer stem cells. The frequencies of AC/GT and TA/TA EST-SSRs were about 10% higher in cancer than normal stem cells. Remarkably, the number of triplets in cancer stem cells was 1.5 times higher than that in normal stem cells. GAT EST-SSR was frequent in cancer stem cells, but, conversely, normal stem cells did not express GAT EST-SSR. We suggest this EST-SSR as a novel triplet in cancer stem cell induction. Translating EST-SSRs to amino acids demonstrated that Asp and Ile were more abundant in cancer stem cells compared to normal stem cells. Finally, Gene Ontology (GO) enrichment analysis was carried out on genes containing triplet SSRs and showed that SSRs intentionally visit some specific GO classes. Interestingly, a NF-kappa (nuclear factor-kB) binding transcription factor was significantly hit by SSR instability which is a hallmark for leukemia stem cells. NF-kappa is an over represented transcription factor during cancer progression. It seems that there is a crosstalk between the NF-kB transcription factor and expressed GAT tandem repeat which negatively regulate apoptosis. In addition to better understanding of tumorigenesis, the findings of this study offer new DNA markers for diagnostic purposes and identifying at risk populations. In addition, a new approach for gene discovery in cancer by target analysis of differentially expressed EST-SSRs between cancer and normal stem cells is presented here.

Construction of a SSR-based genetic map and identification of QTLs for catechins content in tea plant (Camellia sinensis)

Catechins are the most important bioactive compounds in tea, and have been demonstrated to possess a wide variety of pharmacological activities. To characterize quantitative trait loci (QTLs) for catechins content in the tender shoots of tea plant, we constructed a moderately saturated genetic map using 406 simple sequence repeat (SSR) markers, based on a pseudo-testcross population of 183 individuals derived from an intraspecific cross of two Camellia sinensis varieties with diverse catechins composition. The map consisted of fifteen linkage groups (LGs), corresponding to the haploid chromosome number of tea plant (2n = 2x = 30). The total map length was 1,143.5 cM, with an average locus spacing of 2.9 cM. A total of 25 QTLs associated with catechins content were identified over two measurement years. Of these, nine stable QTLs were validated across years, and clustered into four main chromosome regions on LG03, LG11, LG12 and LG15. The population variability explained by each QTL was predominantly at moderate-to-high levels and ranged from 2.4% to 71.0%, with an average of 17.7%. The total number of QTL for each trait varied from four to eight, while the total population variability explained by all QTLs for a trait ranged between 38.4% and 79.7%. This is the first report on the identification of QTL for catechins content in tea plant. The results of this study provide a foundation for further cloning and functional characterization of catechin QTLs for utilization in improvement of tea plant.

Development and characterization of simple sequence repeat (SSR) markers based on RNA-sequencing of Medicago sativa and in silico mapping onto the M. truncatula genome

Sufficient codominant genetic markers are needed for various genetic investigations in alfalfa since the species is an outcrossing autotetraploid. With the newly developed next generation sequencing technology, a large amount of transcribed sequences of alfalfa have been generated and are available for identifying SSR markers by data mining. A total of 54,278 alfalfa non-redundant unigenes were assembled through the Illumina HiSeqTM 2000 sequencing technology. Based on 3,903 unigene sequences, 4,493 SSRs were identified. Tri-nucleotide repeats (56.71%) were the most abundant motif class while AG/CT (21.7%), AGG/CCT (19.8%), AAC/GTT (10.3%), ATC/ATG (8.8%), and ACC/GGT (6.3%) were the subsequent top five nucleotide repeat motifs. Eight hundred and thirty- seven EST-SSR primer pairs were successfully designed. Of these, 527 (63%) primer pairs yielded clear and scored PCR products and 372 (70.6%) exhibited polymorphisms. High transferability was observed for ssp falcata at 99.2% (523) and 71.7% (378) in M. truncatula. In addition, 313 of 527 SSR marker sequences were in silico mapped onto the eight M. truncatula chromosomes. Thirty-six polymorphic SSR primer pairs were used in the genetic relatedness analysis of 30 Chinese alfalfa cultivated accessions generating a total of 199 scored alleles. The mean observed heterozygosity and polymorphic information content were 0.767 and 0.635, respectively. The codominant markers not only enriched the current resources of molecular markers in alfalfa, but also would facilitate targeted investigations in marker-trait association, QTL mapping, and genetic diversity analysis in alfalfa.

Genetic diversity and domestication origin of tea plant Camellia taliensis (Theaceae) as revealed by microsatellite markers

BACKGROUND

Tea is one of the most popular beverages in the world. Many species in the Thea section of the Camellia genus can be processed for drinking and have been domesticated. However, few investigations have focused on the genetic consequence of domestication and geographic origin of landraces on tea plants using credible wild and planted populations of a single species. Here, C. taliensis provides us with a unique opportunity to explore these issues.

RESULTS

Fourteen nuclear microsatellite loci were employed to determine the genetic diversity and domestication origin of C. taliensis, which were represented by 587 individuals from 25 wild, planted and recently domesticated populations. C. taliensis showed a moderate high level of overall genetic diversity. The greater reduction of genetic diversity and stronger genetic drift were detected in the wild group than in the recently domesticated group, indicating the loss of genetic diversity of wild populations due to overexploitation and habitat fragmentation. Instead of the endangered wild trees, recently domesticated individuals were used to compare with the planted trees for detecting the genetic consequence of domestication. A little and non-significant reduction in genetic diversity was found during domestication. The long life cycle, selection for leaf traits and gene flow between populations will delay the emergence of bottleneck in planted trees. Both phylogenetic and assignment analyses suggested that planted trees may have been domesticated from the adjacent central forest of western Yunnan and dispersed artificially to distant places.

CONCLUSIONS

This study contributes to the knowledge about levels and distribution of genetic diversity of C. taliensis and provides new insights into genetic consequence of domestication and geographic origin of planted trees of this species. As an endemic tea source plant, wild, planted and recently domesticated C. taliensis trees should all be protected for their unique genetic characteristics, which are valuable for tea breeding.

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

Despite the worldwide consumption and high economic importance of tea, the plant (Camellia sinensis) is not well studied in molecular biology. Under the few circumstances in which the plant is studied, C. sinensis flowers, which are important for reproduction and cross-breeding, receive less emphasis than investigation of its leaves or roots. Using high-throughput Illumina RNA sequencing, we analyzed a C. sinensis floral transcriptome, and 26.9 million clean reads were assembled into 75,531 unigenes averaging 402 bp. Among them, 50,792 (67.2%) unigenes were annotated with a BLAST search against the NCBI Non-Redundant (NR) database and 10,290 (16.67%) were detected that contained one or more simple sequence repeats (SSRs). From these SSR-containing sequences, 2,439 candidate SSR markers were developed and 720 were experimentally tested, validating 431 (59.9%) novel polymorphic SSR markers for C. sinensis. Then, a consensus SSR-based linkage map was constructed that covered 1,156.9 cM with 237 SSR markers distributed in 15 linkage groups. Both transcriptome information and the genetic map of C. sinensis presented here offer a valuable foundation for molecular biology investigations such as functional gene isolation, quantitative trait loci mapping, and marker-assisted selection breeding in this important species.

Revisiting an important component of plant genomes: microsatellites

Microsatellites are some of the most highly variable repetitive DNA tracts in genomes. Few studies focus on whether the characteristic instability of microsatellites is linked to phenotypic effects in plants. We summarise recent data to investigate how microsatellite variations affect gene expression and hence phenotype. We discuss how the basic characteristics of microsatellites may contribute to phenotypic effects. In summary, microsatellites in plants are universal and highly mutable, they coexist and coevolve with transposable elements, and are under selective pressure. The number of motif nucleotides, the type of motif and transposon activity all contribute to the nonrandom generation and decay of microsatellites, and to conservation and distribution biases. Although microsatellites are generated by accident, they mature through responses to environmental change before final decay. This process is mediated by organism adjustment mechanisms, which maintain a balance between birth versus death and growth versus decay in microsatellites. Close relationships also exist between the physical structure, variation and functionality of microsatellites: in most plant species, sequences containing microsatellites are associated with catalytic activity and binding functions, are expressed in the membrane and organelles, and participate in the developmental and metabolic processes. Microsatellites contribute to genome structure and functional plasticity, and may be considered to promote species evolution in plants in response to environmental changes. In conclusion, the generation, loss, functionality and evolution of microsatellites can be related to plant gene expression and functional alterations. The effect of microsatellites on phenotypic variation may be as significant in plants as it is in animals.

Development, cross-species/genera transferability of novel EST-SSR markers and their utility in revealing population structure and genetic diversity in sugarcane

Sugarcane (Saccharum spp. hybrid) with complex polyploid genome requires a large number of informative DNA markers for various applications in genetics and breeding. Despite the great advances in genomic technology, it is observed in several crop species, especially in sugarcane, the availability of molecular tools such as microsatellite markers are limited. Now-a-days EST-SSR markers are preferred to genomic SSR (gSSR) as they represent only the functional part of the genome, which can be easily associated with desired trait. The present study was taken up with a new set of 351 EST-SSRs developed from the 4085 non redundant EST sequences of two Indian sugarcane cultivars. Among these EST-SSRs, TNR containing motifs were predominant with a frequency of 51.6%. Thirty percent EST-SSRs showed homology with annotated protein. A high frequency of SSRs was found in the 5'UTR and in the ORF (about 27%) and a low frequency was observed in the 3'UTR (about 8%). Two hundred twenty-seven EST-SSRs were evaluated, in sugarcane, allied genera of sugarcane and cereals, and 134 of these have revealed polymorphism with a range of PIC value 0.12 to 0.99. The cross transferability rate ranged from 87.0% to 93.4% in Saccharum complex, 80.0% to 87.0% in allied genera, and 76.0% to 80.0% in cereals. Cloning and sequencing of EST-SSR size variant amplicons revealed that the variation in the number of repeat-units was the main source of EST-SSR fragment polymorphism. When 124 sugarcane accessions were analyzed for population structure using model-based approach, seven genetically distinct groups or admixtures thereof were observed in sugarcane. Results of principal coordinate analysis or UPGMA to evaluate genetic relationships delineated also the 124 accessions into seven groups. Thus, a high level of polymorphism adequate genetic diversity and population structure assayed with the EST-SSR markers not only suggested their utility in various applications in genetics and genomics in sugarcane but also enriched the microsatellite marker resources in sugarcane.

Three EST-SSR markers associated with QTL for the growth of the clam Meretrix meretrix revealed by selective genotyping

The clam Meretrix meretrix is a member of widely cultured, commercially important clams. A marker-trait association analysis was performed using expressed sequence tag (EST) simple sequence repeat (SSR) markers for marker-assisted selection in M. meretrix. Three markers, MM1272, MM2034, and MM7721, were found to be significantly associated with quantitative trait loci (QTLs) controlling shell length (P < 0.0001) in clams of a fast-growing population (JSF) and a control population (JSC). The 144-bp allele of MM1272, the 154-bp allele of MM2034, and the 152- and 165-bp alleles of MM7721 showed a significantly higher frequency in the JSF population (17.65, 36.41, 28.67, and 29.33 %) than in the JSC population (4.65, 8.33, 3.47, and 5.56 %). The three markers showed lower values for the number of alleles and observed heterozygosity as well as a higher proportion of homozygotes in JSF than in JSC population. The three markers have been further confirmed in the high and low tails of another population (09G₃SPSB); similarly, lower values for the number of alleles and observed heterozygosity as well as a higher proportion of homozygotes were found in 09G₃SPSB(H). The putative functions of the three gene fragments containing MM1272, MM2034, and MM7721 also suggested that the three SSR-containing genes might be involved in growth of M. meretrix. All the results suggest that the three EST-SSR markers associated with growth QTLs would be useful for marker-assisted selection in M. meretrix breeding.

Development of simple sequence repeat markers and diversity analysis in alfalfa (Medicago sativa L.)

Efficient and robust molecular markers are essential for molecular breeding in plant. Compared to dominant and bi-allelic markers, multiple alleles of simple sequence repeat (SSR) markers are particularly informative and superior in genetic linkage map and QTL mapping in autotetraploid species like alfalfa. The objective of this study was to enrich SSR markers directly from alfalfa expressed sequence tags (ESTs). A total of 12,371 alfalfa ESTs were retrieved from the National Center for Biotechnology Information. Total 774 SSR-containing ESTs were identified from 716 ESTs. On average, one SSR was found per 7.7 kb of EST sequences. Tri-nucleotide repeats (48.8 %) was the most abundant motif type, followed by di-(26.1 %), tetra-(11.5 %), penta-(9.7 %), and hexanucleotide (3.9 %). One hundred EST-SSR primer pairs were successfully designed and 29 exhibited polymorphism among 28 alfalfa accessions. The allele number per marker ranged from two to 21 with an average of 6.8. The PIC values ranged from 0.195 to 0.896 with an average of 0.608, indicating a high level of polymorphism of the EST-SSR markers. Based on the 29 EST-SSR markers, assessment of genetic diversity was conducted and found that Medicago sativa ssp. sativa was clearly different from the other subspecies. The high transferability of those EST-SSR markers was also found for relative species.

Comparative mapping in the Fagaceae and beyond with EST-SSRs

BACKGROUND

Genetic markers and linkage mapping are basic prerequisites for comparative genetic analyses, QTL detection and map-based cloning. A large number of mapping populations have been developed for oak, but few gene-based markers are available for constructing integrated genetic linkage maps and comparing gene order and QTL location across related species.

RESULTS

We developed a set of 573 expressed sequence tag-derived simple sequence repeats (EST-SSRs) and located 397 markers (EST-SSRs and genomic SSRs) on the 12 oak chromosomes (2n = 2x = 24) on the basis of Mendelian segregation patterns in 5 full-sib mapping pedigrees of two species: Quercus robur (pedunculate oak) and Quercus petraea (sessile oak). Consensus maps for the two species were constructed and aligned. They showed a high degree of macrosynteny between these two sympatric European oaks. We assessed the transferability of EST-SSRs to other Fagaceae genera and a subset of these markers was mapped in Castanea sativa, the European chestnut. Reasonably high levels of macrosynteny were observed between oak and chestnut. We also obtained diversity statistics for a subset of EST-SSRs, to support further population genetic analyses with gene-based markers. Finally, based on the orthologous relationships between the oak, Arabidopsis, grape, poplar, Medicago, and soybean genomes and the paralogous relationships between the 12 oak chromosomes, we propose an evolutionary scenario of the 12 oak chromosomes from the eudicot ancestral karyotype.

CONCLUSIONS

This study provides map locations for a large set of EST-SSRs in two oak species of recognized biological importance in natural ecosystems. This first step toward the construction of a gene-based linkage map will facilitate the assignment of future genome scaffolds to pseudo-chromosomes. This study also provides an indication of the potential utility of new gene-based markers for population genetics and comparative mapping within and beyond the Fagaceae.

Expressed sequence tags from organ-specific cDNA libraries of tea (Camellia sinensis) and polymorphisms and transferability of EST-SSRs across Camellia species

Tea is one of the most popular beverages in the world and the tea plant, Camellia sinensis (L.) O. Kuntze, is an important crop in many countries. To increase the amount of genomic information available for C. sinensis, we constructed seven cDNA libraries from various organs and used these to generate expressed sequence tags (ESTs). A total of 17,458 ESTs were generated and assembled into 5,262 unigenes. About 50% of the unigenes were assigned annotations by Gene Ontology. Some were homologous to genes involved in important biological processes, such as nitrogen assimilation, aluminum response, and biosynthesis of caffeine and catechins. Digital northern analysis showed that 67 unigenes were expressed differentially among the seven organs. Simple sequence repeat (SSR) motif searches among the unigenes identified 1,835 unigenes (34.9%) harboring SSR motifs of more than six repeat units. A subset of 100 EST-SSR primer sets was tested for amplification and polymorphism in 16 tea accessions. Seventy-one primer sets successfully amplified EST-SSRs and 70 EST-SSR loci were polymorphic. Furthermore, these 70 EST-SSR markers were transferable to 14 other Camellia species. The ESTs and EST-SSR markers will enhance the study of important traits and the molecular genetics of tea plants and other Camellia species.