Although lily is the second largest flower crop in cutting flower commodity, only six simple sequence repeats SSRs have been reported. Thus, we developed expressed sequence tag derived-SSRs (EST-SSRs) for the Lilium genus. Among 2,235 unique ESTs, 754 ESTs contained SSR motifs, among which 165 ESTs were amenable to primer design. Among these 165 EST-SSRs, 131 EST-SSRs showed amplification in at least one Lilium species, and 76 EST-SSRs showed amplification in at least nine species. Of the 76 EST-SSRs, 47 showed amplification in all Lilium species analyzed. Using 10 breeding lines, we selected 19 EST-SSRs that had the highest number of alleles and polymorphism information content. The polymorphism information content values of these selected EST-SSRs ranged from 0.49 to 0.94 with an average of 0.76, which are higher than other plant species. The phylogenetic dendrogram derived from the amplification profiles of the 19 high polymorphic EST-SSRs was congruent with the genetic background of the 84 selected lily accessions and hybrids, which are available in commerce. Thus, the developed EST-SSRs will be very useful in germplasm management, genetic diversity analysis, cultivar finger printing, and molecular breeding in the lily.

BackgroundPapaya is a major fruit crop in tropical and subtropical regions worldwide. It is trioecious with three sex forms: male, female, and hermaphrodite. Sex determination is controlled by a pair of nascent sex chromosomes with two slightly different Y chromosomes, Y for male and Yh for hermaphrodite. The sex chromosome genotypes are XY (male), XYh (hermaphrodite), and XX (female). The papaya hermaphrodite-specific Yh chromosome region (HSY) is pericentromeric and heterochromatic. Physical mapping of HSY and its X counterpart is essential for sequencing these regions and uncovering the early events of sex chromosome evolution and to identify the sex determination genes for crop improvement.ResultsA reiterate chromosome walking strategy was applied to construct the two physical maps with three bacterial artificial chromosome (BAC) libraries. The HSY physical map consists of 68 overlapped BACs on the minimum tiling path, and covers all four HSY-specific Knobs. One gap remained in the region of Knob 1, the only knob structure shared between HSY and X, due to the lack of HSY-specific sequences. This gap was filled on the physical map of the HSY corresponding region in the X chromosome. The X physical map consists of 44 BACs on the minimum tiling path with one gap remaining in the middle, due to the nature of highly repetitive sequences. This gap was filled on the HSY physical map. The borders of the non-recombining HSY were defined genetically by fine mapping using 1460 F2 individuals. The genetically defined HSY spanned approximately 8.5 Mb, whereas its X counterpart extended about 5.4 Mb including a 900 Kb region containing the Knob 1 shared by the HSY and X. The 8.5 Mb HSY corresponds to 4.5 Mb of its X counterpart, showing 4 Mb (89%) DNA sequence expansion.ConclusionThe 89% increase of DNA sequence in HSY indicates rapid expansion of the Yh chromosome after genetic recombination was suppressed 2–3 million years ago. The genetically defined borders coincide with the common BACs on the minimum tiling paths of HSY and X. The minimum tiling paths of HSY and its X counterpart are being used for sequencing these X and Yh-specific regions.

Lipophilic compounds from Korean perilla ( Perilla frutescens ) seeds were characterized to determine the diversity among their phytochemicals and to analyze relationships between their contents. Twenty-four metabolites consisting of policosanol, phytosterol, tocopherol, and fatty acids were identified. The metabolite profiles were subjected to data mining processes, including principal component analysis (PCA), partial least-squares discriminate analysis (PLS-DA), and Pearson's correlation analysis. PLS-DA could distinguish between all cultivars except between Daesil and Daeyeup cultivars. Linolenic acid contents were positively correlated with β-sitosterol (r = 0.8367, P < 0.0001) and γ-tocopherol contents (r = 0. 7201, P < 0.001) among all perilla grains. The Daesil and Daeyeup cultivars appear to be good candidates for future breeding programs because they have simultaneously high linolenic acid, phytosterol, and tocopherol levels. These results demonstrate the use of metabolite profiling as a tool for assessing the quality of food.

Microsatellites, or simple sequence repeats (SSRs), are highly polymorphic and universally distributed in eukaryotes. SSRs have been used extensively as sequence tagged markers in genetic studies. Recently, the functional and evolutionary importance of SSRs has received considerable attention. Here we report the mining and characterization of the SSRs in papaya genome. We analyzed SSRs from 277.4 Mb of whole genome shotgun (WGS) sequences, 51.2 Mb bacterial artificial chromosome (BAC) end sequences (BES), and 13.4 Mb expressed sequence tag (EST) sequences. The papaya SSR density was one SSR per 0.7 kb of DNA sequence in the WGS, which was higher than that in BES and EST sequences. SSR abundance was dramatically reduced as the repeat length increased. According to SSR motif length, dinucleotide repeats were the most common motif in class I, whereas hexanucleotides were the most copious in class II SSRs. The tri-and hexanucleotide repeats of both classes were greater in EST sequences compared to genomic sequences. In class I SSR, AT and AAT were the most frequent motifs in BES and WGS sequences. By contrast, AG and AAG were the most abundant in EST sequences. For SSR marker development, 9,860 primer pairs were surveyed for amplification and polymorphism. Successful amplification and polymorphic rates were 66.6% and 17.6%, respectively. The highest polymorphic rates were achieved by AT, AG, and ATG motifs. The genome wide analysis of microsatellites revealed their frequency and distribution in papaya genome, which varies among plant genomes. This complete set of SSRs markers throughout the genome will assist diverse genetic studies in papaya and related species.

Papaya is a major fruit crop in the tropics and has recently evolved sex chromosomes. Towards sequencing the papaya sex chromosomes, two bacterial artificial chromosome (BAC) libraries were constructed from papaya male and female genomic DNA. The female BAC library was constructed using restriction enzyme BstY I and consists of 36,864 clones with an average insert size of 104 kb, providing 10.3x genome equivalents. The male BAC library was constructed using restriction enzyme EcoR I and consists of 55,296 clones with an average insert size of 101 kb, providing 15.0x genome equivalents. The male BAC library was used in constructing the physical map of the male-specific region of the male Y chromosome (MSY) and in filling gaps and extending the physical map of the hermaphrodite-specific region of the Yh chromosome (HSY) and the X chromosome physical map. The female BAC library was used to extend the X physical map gap. The MSY, HSY, and X physical maps offer a unique opportunity to study chromosomal rearrangements, Y chromosome degeneration, and dosage compensation of the papaya nascent sex chromosomes.

To broaden and delve into the genomic information of , an important medicinal plant in many Asian countries, RNA sequencing (RNA-seq) analysis was performed and a total of 16,638 non-redundant unigenes (≥ 300 bp) with an average length of 755 bp were generated by de novo assembly from 17,580,456 trimmed clear reads. The functional categorization of the identified unigenes by a gene ontology (GO) term resulted in 2305 genes in the cellular component, 5577 in the biological processes, and 8056 in the molecular functions, respectively. The top sub-category in biological processes was the metabolic process with 4374 genes. Among annotated genes, 3006 were mapped to 123 metabolic pathways by the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathway analysis tool. The search for simple sequence repeats (SSRs) resulted in 845 SSRs from 749 SSR-containing unigenes and the most abundant SSR motifs was AAG/CTT with 179 occurrences. Twelve SSR markers were tested for cross transferability among five species; eight of them exhibited polymorphism. Taken together, these data provide valuable resources for genomic or genetic studies of species and other relative studies. The transcriptome shotgun assembly data have been deposited at DDBJ/EMBL/GenBank under the accession GGEM00000000.

Based on sequence similarity to Arabidopsis inositol polyphosphate 5-phosphatases (5PTases) involved in abiotic stress responses and development, four tomato cDNAs (Le5PT1-4) encoding putative 5PTase proteins were identified. The predicted protein sequences of the Le5PTs include conserved catalytic domains required for 5PTase enzyme activity. Le5PT1, 2, and 3 showed high amino acid sequence identity with At5PTase2, At5PTase1 and At5PTase3, and At5PTase5 and At5PTase6, respectively. The expression of Le5PT1 was downregulated soon after initiation of dehydration and salt stress as well as exposure to polyethylene glycol (PEG) and NaCl, but not by exogenous ABA treatment. On the other hand, the expression of Le5PT2 gradually increased with time in all treatments. Transgenic tobacco plants overexpressing Le5PT1 exhibited reduced growth in height, leaf area, and dry weight compared to wild type plants. Transgenic plants also had lower water use efficiency (WUE) than wild type and the downregulation of the drought-responsive gene, NtERD10B. Together these results suggest that Le5PT1 may have a negative role in response to water deficit through the repression of drought-inducible genes that in turn affects plant growth and development.

Pistacia chinensis Bunge has not only been used as a medicinal plant to treat various illnesses but its young shoots and leaves have also been used as vegetables. In addition, P. chinensis is used as a rootstock for Pistacia vera (pistachio). Here, the transcriptome of P. chinensis was sequenced to enrich genetic resources and identify secondary metabolite biosynthetic pathways using Illumina RNA-seq methods. De novo assembly resulted in 18,524 unigenes with an average length of 873 bp from 19 million RNA-seq reads. A Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation tool assigned KO (KEGG orthology) numbers to 6,553 (36.2%) unigenes, among which 4,061 unigenes were mapped into 391 different metabolic pathways. For terpenoid backbone and carotenoid biosynthesis pathways, 44 and 22 unigenes encode enzymes corresponding to 30 and 16 entries, respectively. Twentytwo unigenes encode proteins for 16 entries of the carotenoid biosynthesis pathway. As for the phenylpropanoid and flavonoid biosynthesis pathways, 63 and 24 unigenes were homologous to 17 and 14 entry proteins, respectively. Mining of simple sequence repeat identified 2,599 simple sequence repeats from P. chinensis unigenes. The results of the present study provide a valuable resource for in-depth studies on comparative and functional genomics to unravel the underlying mechanisms of the medicinal properties of Pistacia L.