Development and characterization of microsatellite markers for the Cape gooseberry Physalis peruviana

Development of Chloroplast Microsatellite Markers and Evaluation of Genetic Diversity and Population Structure of Cutleaf Groundcherry (Physalis angulata L.) in China

Cutleaf groundcherry (Physalis angulata L.), an annual plant containing a variety of active ingredients, has great medicinal value. However, studies on the genetic diversity and population structure of P. angulata are limited. In this study, we developed chloroplast microsatellite (cpSSR) markers and applied them to evaluate the genetic diversity and population structure of P. angulata. A total of 57 cpSSRs were identified from the chloroplast genome of P. angulata. Among all cpSSR loci, mononucleotide markers were the most abundant (68.24%), followed by tetranucleotide (12.28%), dinucleotide (10.53%), and trinucleotide (8.77%) markers. In total, 30 newly developed cpSSR markers with rich polymorphism and good stability were selected for further genetic diversity and population structure analyses. These cpSSRs amplified a total of 156 alleles, 132 (84.62%) of which were polymorphic. The percentage of polymorphic alleles and the average polymorphic information content (PIC) value of the cpSSRs were 81.29% and 0.830, respectively. Population genetic diversity analysis indicated that the average observed number of alleles (Na), number of effective alleles (He), Nei's gene diversity (h), and Shannon information indices (I) of 16 P. angulata populations were 1.3161, 1.1754, 0.1023, and 0.1538, respectively. Moreover, unweighted group arithmetic mean, neighbor-joining, principal coordinate, and STRUCTURE analyses indicated that 203 P. angulata individuals from 16 populations were grouped into four clusters. A molecular variance analysis (AMOVA) illustrated the considerable genetic variation among populations, while the gene flow (Nm) value (0.2324) indicated a low level of gene flow among populations. Our study not only provided a batch of efficient genetic markers for research on P. angulata but also laid an important foundation for the protection and genetic breeding of P. angulata resources.

Complete Plastome of Physalis angulata var. villosa, Gene Organization, Comparative Genomics and Phylogenetic Relationships among Solanaceae

Physalis angulata var. villosa, rich in withanolides, has been used as a traditional Chinese medicine for many years. To date, few extensive molecular studies of this plant have been conducted. In the present study, the plastome of P. angulata var. villosa was sequenced, characterized and compared with that of other Physalis species, and a phylogenetic analysis was conducted in the family Solanaceae. The plastome of P. angulata var. villosa was 156,898 bp in length with a GC content of 37.52%, and exhibited a quadripartite structure typical of land plants, consisting of a large single-copy (LSC, 87,108 bp) region, a small single-copy (SSC, 18,462 bp) region and a pair of inverted repeats (IR: IRA and IRB, 25,664 bp each). The plastome contained 131 genes, of which 114 were unique and 17 were duplicated in IR regions. The genome consisted of 85 protein-coding genes, eight rRNA genes and 38 tRNA genes. A total of 38 long, repeat sequences of three types were identified in the plastome, of which forward repeats had the highest frequency. Simple sequence repeats (SSRs) analysis revealed a total of 57 SSRs, of which the T mononucleotide constituted the majority, with most of SSRs being located in the intergenic spacer regions. Comparative genomic analysis among nine Physalis species revealed that the single-copy regions were less conserved than the pair of inverted repeats, with most of the variation being found in the intergenic spacer regions rather than in the coding regions. Phylogenetic analysis indicated a close relationship between Physalis and Withania. In addition, Iochroma, Dunalia, Saracha and Eriolarynx were paraphyletic, and clustered together in the phylogenetic tree. Our study published the first sequence and assembly of the plastome of P. angulata var. villosa, reported its basic resources for evolutionary studies and provided an important tool for evaluating the phylogenetic relationship within the family Solanaceae.

Putative Novel Effector Genes Revealed by the Genomic Analysis of the Phytopathogenic Fungus Fusarium oxysporum f. sp. physali (Foph) That Infects Cape Gooseberry Plants

The vascular wilt disease caused by the fungus Fusarium oxysporum f. sp. physali (Foph) is one of the most limiting factors for the production and export of cape gooseberry (Physalis peruviana) in Colombia. A transcriptomic analysis of a highly virulent strain of F. oxysporum in cape gooseberry plants, revealed the presence of secreted in the xylem (SIX) effector genes, known to be involved in the pathogenicity of other formae speciales (ff. spp.) of F. oxysporum. This pathogenic strain was classified as a new f. sp. named Foph, due to its specificity for cape gooseberry hosts. Here, we sequenced and assembled the genome of five strains of F. oxysporum from a fungal collection associated to the cape gooseberry crop (including Foph), focusing on the validation of the presence of SIX homologous and on the identification of putative effectors unique to Foph. By comparative and phylogenomic analyses based on single-copy orthologous, we found that Foph is closely related to F. oxysporum ff. spp., associated with solanaceous hosts. We confirmed the presence of highly identical homologous genomic regions between Foph and Fol that contain effector genes and identified six new putative effector genes, specific to Foph pathogenic strains. We also conducted a molecular characterization using this set of putative novel effectors in a panel of 36 additional stains of F. oxysporum including two of the four sequenced strains, from the fungal collection mentioned above. These results suggest the polyphyletic origin of Foph and the putative independent acquisition of new candidate effectors in different clades of related strains. The novel effector candidates identified in this genomic analysis, represent new sources involved in the interaction between Foph and cape gooseberry, that could be implemented to develop appropriate management strategies of the wilt disease caused by Foph in the cape gooseberry crop.

Optimization of the genotyping-by-sequencing SNP calling for diversity analysis in cape gooseberry (Physalis peruviana L.) and related taxa

A robust Genotyping-By-Sequencing (GBS) pipeline platform was examined to provide accurate discovery of Single Nucleotide Polymorphisms (SNPs) in a cape gooseberry (Physalis peruviana L.) and related taxa germplasm collection. A total of 176 accessions representing, wild, weedy, and commercial cultivars as well as related taxa from the Colombian germplasm bank and other world repositories were screened using GBS. The pipeline parameters mnLCov of 0.5 and a mnScov of 0.7, tomato and potato genomes, and cape gooseberry transcriptome for read alignments, were selected to better assess diversity and population structure in cape gooseberry and related taxa. A total of 7,425 SNPs, derived from P. peruviana common tags (unique 64 bp sequences shared between selected species), were used. Within P. peruviana, five subpopulations with a high genetic diversity and allele fixation (H_E: 0.35 to 0.36 and F_IS: -0.11 to -0.01, respectively) were detected. Conversely, low genetic differentiation (F_ST: 0.01 to 0.05) was also observed, indicating a high gene flow among subpopulations. These results contribute to the establishment of adequate conservation and breeding strategies for Cape gooseberry and closely related Physalis species.

Genome-wide microsatellite characteristics of five human Plasmodium species, focusing on Plasmodium malariae and P. ovale curtisi

Microsatellites can be utilized to explore genotypes, population structure, and other genomic features of eukaryotes. Systematic characterization of microsatellites has not been a focus for several species of Plasmodium, including P. malariae and P. ovale, as the majority of malaria elimination programs are focused on P. falciparum and to a lesser extent P. vivax. Here, five human malaria species (P. falciparum, P. vivax, P. malariae, P. ovale curtisi, and P. knowlesi) were investigated with the aim of conducting in-depth categorization of microsatellites for P. malariae and P. ovale curtisi. Investigation of reference genomes for microsatellites with unit motifs of 1–10 base pairs indicates high diversity among the five Plasmodium species. Plasmodium malariae, with the largest genome size, displays the second highest microsatellite density (1421 No./Mbp; 5% coverage) next to P. falciparum (3634 No./Mbp; 12% coverage). The lowest microsatellite density was observed in P. vivax (773 No./Mbp; 2% coverage). A, AT, and AAT are the most commonly repeated motifs in the Plasmodium species. For P. malariae and P. ovale curtisi, microsatellite-related sequences are observed in approximately 18–29% of coding sequences (CDS). Lysine, asparagine, and glutamic acids are most frequently coded by microsatellite-related CDS. The majority of these CDS could be related to the gene ontology terms “cell parts,” “binding,” “developmental processes,” and “metabolic processes.” The present study provides a comprehensive overview of microsatellite distribution and can assist in the planning and development of potentially useful genetic tools for further investigation of P. malariae and P. ovale curtisi epidemiology.

Polymorphic markers for identification of parasite population in Plasmodium malariae

Background

Molecular genotyping in Plasmodium serves many aims including providing tools for studying parasite population genetics and distinguishing recrudescence from reinfection. Microsatellite typing, insertion-deletion (INDEL) and single nucleotide polymorphisms is used for genotyping, but only limited information is available for Plasmodium malariae, an important human malaria species. This study aimed to provide a set of genetic markers to facilitate the study of P. malariae population genetics.

Methods

Markers for microsatellite genotyping and pmmsp1 gene polymorphisms were developed and validated in symptomatic P. malariae field isolates from Myanmar (N = 37). Fragment analysis was used to determine allele sizes at each locus to calculate multiplicity of infections (MOI), linkage disequilibrium, heterozygosity and construct dendrograms. Nucleotide diversity (π), number of haplotypes, and genetic diversity (H_d) were assessed and a phylogenetic tree was constructed. Genome-wide microsatellite maps with annotated regions of newly identified markers were constructed.

Results

Six microsatellite markers were developed and tested in 37 P. malariae isolates which showed sufficient heterozygosity (0.530–0.922), and absence of linkage disequilibrium (I_A^S=0.03, p value > 0.05) (N = 37). In addition, a tandem repeat (VNTR)-based pmmsp1 INDEL polymorphisms marker was developed and assessed in 27 P. malariae isolates showing a nucleotide diversity of 0.0976, haplotype gene diversity of 0.698 and identified 14 unique variants. The size of VNTR consensus repeat unit adopted as allele was 27 base pairs. The markers Pm12_426 and pmmsp1 showed greatest diversity with heterozygosity scores of 0.920 and 0.835, respectively. Using six microsatellites markers, the likelihood that any two parasite strains would have the same microsatellite genotypes was 8.46 × 10⁻⁴ and was further reduced to 1.66 × 10⁻⁴ when pmmsp1 polymorphisms were included.

Conclusions

Six novel microsatellites genotyping markers and a set of pmmsp1 VNTR-based INDEL polymorphisms markers for P. malariae were developed and validated. Each marker could be independently or in combination employed to access genotyping of the parasite. The newly developed markers may serve as a useful tool for investigating parasite diversity, population genetics, molecular epidemiology and for distinguishing recrudescence from reinfection in drug efficacy studies.

A Transcript and Metabolite Atlas of Blackcurrant Fruit Development Highlights Hormonal Regulation and Reveals the Role of Key Transcription Factors

Blackcurrant fruit collected at six stages of development were assessed for changes in gene expression using custom whole transcriptome microarrays and for variation in metabolite content using a combination of liquid chromatography-mass spectrometry and gas chromatography-mass spectrometry. Principal components analysis demonstrated that fruit development could be clearly defined according to their transcript or metabolite profiles. During early developmental stages, metabolite profiles were dominated by amino acids and tannins, whilst transcript profiles were enriched in functions associated with cell division, anatomical structure morphogenesis and cell wall metabolism. During mid fruit development, fatty acids accumulated and transcript profiles were consistent with seed and embryo development. At the later stages, sugars and anthocyanins accumulated consistent with transcript profiles that were associated with secondary metabolism. Transcript data also indicated active signaling during later stages of fruit development. A targeted analysis of signaling networks revealed a dynamic activation and repression of almost 60 different transcripts encoding transcription factors across the course of fruit development, many of which have been demonstrated as pivotal to controlling such processes in other species. Transcripts associated with cytokinin and gibberellin were highly abundant at early fruit development, whilst those associated with ABA and ethylene tended to be more abundant at later stages. The data presented here provides an insight into fruit development in blackcurrant and provides a foundation for further work in the elucidation of the genetic basis of fruit quality.

NMR-based metabolic study of fruits of Physalis peruviana L. grown in eight different Peruvian ecosystems

The berry of Physalis peruviana L. (Solanaceae) represents an important socio-economical commodity for Latin America. The absence of a clear phenotype renders it difficult to trace its place of origin. In this study, Cape gooseberries from eight different regions within the Peruvian Andes were profiled for their metabolism implementing a NMR platform. Twenty-four compounds could be unequivocally identified and sixteen quantified. One-way ANOVA and post-hoc Tukey test revealed that all of the quantified metabolites changed significantly among regions: Bambamarca I showed the most accumulated significant differences. The coefficient of variation demonstrated high phenotypic plasticity for amino acids, while sugars displayed low phenotypic plasticity. Correlation analysis highlighted the closely coordinated behavior of the amino acid profile. Finally, PLS-DA revealed a clear separation among the regions based on their metabolic profiles, accentuating the discriminatory capacity of NMR in establishing significant phytochemical differences between producing regions of the fruit of P. peruviana L.

Workflow and web application for annotating NCBI BioProject transcriptome data

The volume of transcriptome data is growing exponentially due to rapid improvement of experimental technologies. In response, large central resources such as those of the National Center for Biotechnology Information (NCBI) are continually adapting their computational infrastructure to accommodate this large influx of data. New and specialized databases, such as Transcriptome Shotgun Assembly Sequence Database (TSA) and Sequence Read Archive (SRA), have been created to aid the development and expansion of centralized repositories. Although the central resource databases are under continual development, they do not include automatic pipelines to increase annotation of newly deposited data. Therefore, third-party applications are required to achieve that aim. Here, we present an automatic workflow and web application for the annotation of transcriptome data. The workflow creates secondary data such as sequencing reads and BLAST alignments, which are available through the web application. They are based on freely available bioinformatics tools and scripts developed in-house. The interactive web application provides a search engine and several browser utilities. Graphical views of transcript alignments are available through SeqViewer, an embedded tool developed by NCBI for viewing biological sequence data. The web application is tightly integrated with other NCBI web applications and tools to extend the functionality of data processing and interconnectivity. We present a case study for the species Physalis peruviana with data generated from BioProject ID 67621.

Database URL: http://www.ncbi.nlm.nih.gov/projects/physalis/

Characterization of the Kenaf (Hibiscus cannabinus) Global Transcriptome Using Illumina Paired-End Sequencing and Development of EST-SSR Markers

Kenaf (Hibiscus cannabinus L.) is an economically important natural fiber crop grown worldwide. However, only 20 expressed tag sequences (ESTs) for kenaf are available in public databases. The aim of this study was to develop large-scale simple sequence repeat (SSR) markers to lay a solid foundation for the construction of genetic linkage maps and marker-assisted breeding in kenaf. We used Illumina paired-end sequencing technology to generate new EST-simple sequences and MISA software to mine SSR markers. We identified 71,318 unigenes with an average length of 1143 nt and annotated these unigenes using four different protein databases. Overall, 9324 complementary pairs were designated as EST-SSR markers, and their quality was validated using 100 randomly selected SSR markers. In total, 72 primer pairs reproducibly amplified target amplicons, and 61 of these primer pairs detected significant polymorphism among 28 kenaf accessions. Thus, in this study, we have developed large-scale SSR markers for kenaf, and this new resource will facilitate construction of genetic linkage maps, investigation of fiber growth and development in kenaf, and also be of value to novel gene discovery and functional genomic studies.

Characterization of the Kenaf (Hibiscus cannabinus) Global Transcriptome Using Illumina Paired-End Sequencing and Development of EST-SSR Markers

Kenaf (Hibiscus cannabinus L.) is an economically important natural fiber crop grown worldwide. However, only 20 expressed tag sequences (ESTs) for kenaf are available in public databases. The aim of this study was to develop large-scale simple sequence repeat (SSR) markers to lay a solid foundation for the construction of genetic linkage maps and marker-assisted breeding in kenaf. We used Illumina paired-end sequencing technology to generate new EST-simple sequences and MISA software to mine SSR markers. We identified 71,318 unigenes with an average length of 1143 nt and annotated these unigenes using four different protein databases. Overall, 9324 complementary pairs were designated as EST-SSR markers, and their quality was validated using 100 randomly selected SSR markers. In total, 72 primer pairs reproducibly amplified target amplicons, and 61 of these primer pairs detected significant polymorphism among 28 kenaf accessions. Thus, in this study, we have developed large-scale SSR markers for kenaf, and this new resource will facilitate construction of genetic linkage maps, investigation of fiber growth and development in kenaf, and also be of value to novel gene discovery and functional genomic studies.

Genetic diversity and population structure in Physalis peruviana and related taxa based on InDels and SNPs derived from COSII and IRG markers

The genus Physalis is common in the Americas and includes several economically important species, among them Physalis peruviana that produces appetizing edible fruits. We studied the genetic diversity and population structure of P. peruviana and characterized 47 accessions of this species along with 13 accessions of related taxa consisting of 222 individuals from the Colombian Corporation of Agricultural Research (CORPOICA) germplasm collection, using Conserved Orthologous Sequences (COSII) and Immunity Related Genes (IRGs). In addition, 642 Single Nucleotide Polymorphism (SNPs) markers were identified and used for the genetic diversity analysis. A total of 121 alleles were detected in 24 InDels loci ranging from 2 to 9 alleles per locus, with an average of 5.04 alleles per locus. The average number of alleles in the SNP markers was two. The observed heterozygosity for P. peruviana with InDel and SNP markers was higher (0.48 and 0.59) than the expected heterozygosity (0.30 and 0.41). Interestingly, the observed heterozygosity in related taxa (0.4 and 0.12) was lower than the expected heterozygosity (0.59 and 0.25). The coefficient of population differentiation F_ST was 0.143 (InDels) and 0.038 (SNPs), showing a relatively low level of genetic differentiation among P. peruviana and related taxa. Higher levels of genetic variation were instead observed within populations based on the AMOVA analysis. Population structure analysis supported the presence of two main groups and PCA analysis based on SNP markers revealed two distinct clusters in the P. peruviana accessions corresponding to their state of cultivation. In this study, we identified molecular markers useful to detect genetic variation in Physalis germplasm for assisting conservation and crossbreeding strategies.

Cape Gooseberry [Physalis peruviana L.] Calyces Ameliorate TNBS Acid-induced Colitis in Rats

BACKGROUND AND AIMS

Physalis peruviana [cape gooseberry] is highly appreciated for its commercial value. The Colombian ecotype is in great demand in the international market, particularly for the unique morphological characteristics of the calyx, which has extended use as a traditional herbal remedy in Colombia because of its anti-inflammatory properties. In this work, the anti-inflammatory activity of the total ethereal extract of Physalis peruviana calyces was evaluated in preventive and therapeutic protocols in a TNBS acid-induced colitis rat model.

METHODS

Colitis was induced by intrarectal administration of TNBS. An evaluation of macroscopic and histopathological parameters in colonic tissue was performed, along with the determination of myeloperoxidase enzyme activity, cytokine levels and gene expression. Additionally, effects on nitric oxide release by lipopolysaccharide [LPS]-stimulated RAW264.7 macrophages and the scavenging activity of DPPH and ABTS free radicals were determined.

RESULTS

The treatment with the Physalis peruviana extract produced a significant improvement in the colonic tissue at both macroscopic and histological levels. IL-1β and TNF-α production was reduced by the extract in both experimental approaches. The groups treated with Physalis peruviana showed a tendency to MUC2 up-regulation and down-regulation of COX-2, iNOS, NLRP3, IL-1β, IL-6 and IL-10 expression. Nitric oxide release in RAW264.7 macrophages was significantly inhibited.

CONCLUSIONS

The Physalis peruviana extract showed intestinal anti-inflammatory activity in the TNBS-induced colitis model, placing this species' calyx, a natural derivative, as a promising source of metabolites that could be used in treatment for inflammatory bowel disease.

Searching for repeats, as an example of using the generalised Ruzzo-Tompa algorithm to find optimal subsequences with gaps

Some biological sequences contain subsequences of unusual composition; e.g. some proteins contain DNA binding domains, transmembrane regions and charged regions, and some DNA sequences contain repeats. The linear-time Ruzzo-Tompa (RT) algorithm finds subsequences of unusual composition, using a sequence of scores as input and the corresponding 'maximal segments' as output. In principle, permitting gaps in the output subsequences could improve sensitivity. Here, the input of the RT algorithm is generalised to a finite, totally ordered, weighted graph, so the algorithm locates paths of maximal weight through increasing but not necessarily adjacent vertices. By permitting the penalised deletion of unfavourable letters, the generalisation therefore includes gaps. The program RepWords, which finds inexact simple repeats in DNA, exemplifies the general concepts by out-performing a similar extant, ad hoc tool. With minimal programming effort, the generalised Ruzzo-Tompa algorithm could improve the performance of many programs for finding biological subsequences of unusual composition.

Identification of immunity related genes to study the Physalis peruviana--Fusarium oxysporum pathosystem

The Cape gooseberry (Physalisperuviana L) is an Andean exotic fruit with high nutritional value and appealing medicinal properties. However, its cultivation faces important phytosanitary problems mainly due to pathogens like Fusarium oxysporum, Cercosporaphysalidis and Alternaria spp. Here we used the Cape gooseberry foliar transcriptome to search for proteins that encode conserved domains related to plant immunity including: NBS (Nucleotide Binding Site), CC (Coiled-Coil), TIR (Toll/Interleukin-1 Receptor). We identified 74 immunity related gene candidates in P. peruviana which have the typical resistance gene (R-gene) architecture, 17 Receptor like kinase (RLKs) candidates related to PAMP-Triggered Immunity (PTI), eight (TIR-NBS-LRR, or TNL) and nine (CC–NBS-LRR, or CNL) candidates related to Effector-Triggered Immunity (ETI) genes among others. These candidate genes were categorized by molecular function (98%), biological process (85%) and cellular component (79%) using gene ontology. Some of the most interesting predicted roles were those associated with binding and transferase activity. We designed 94 primers pairs from the 74 immunity-related genes (IRGs) to amplify the corresponding genomic regions on six genotypes that included resistant and susceptible materials. From these, we selected 17 single band amplicons and sequenced them in 14 F. oxysporum resistant and susceptible genotypes. Sequence polymorphisms were analyzed through preliminary candidate gene association, which allowed the detection of one SNP at the PpIRG-63 marker revealing a nonsynonymous mutation in the predicted LRR domain suggesting functional roles for resistance.

Development and characterization of 1,827 expressed sequence tag-derived simple sequence repeat markers for ramie (Boehmeria nivea L. Gaud)

Ramie (Boehmeria nivea L. Gaud) is one of the most important natural fiber crops, and improvement of fiber yield and quality is the main goal in efforts to breed superior cultivars. However, efforts aimed at enhancing the understanding of ramie genetics and developing more effective breeding strategies have been hampered by the shortage of simple sequence repeat (SSR) markers. In our previous study, we had assembled de novo 43,990 expressed sequence tags (ESTs). In the present study, we searched these previously assembled ESTs for SSRs and identified 1,685 ESTs (3.83%) containing 1,878 SSRs. Next, we designed 1,827 primer pairs complementary to regions flanking these SSRs, and these regions were designated as SSR markers. Among these markers, dinucleotide and trinucleotide repeat motifs were the most abundant types (36.4% and 36.3%, respectively), whereas tetranucleotide, pentanucleotide, and hexanucleotide motifs represented <10% of the markers. The motif AG/CT was the most abundant, accounting for 28.74% of the markers. One hundred EST-SSR markers (97 SSRs located in genes encoding transcription factors and 3 SSRs in genes encoding cellulose synthases) were amplified using polymerase chain reaction for detecting 24 ramie varieties. Of these 100 markers, 98 markers were successfully amplified and 81 markers were polymorphic, with 2–6 alleles among the 24 varieties. Analysis of the genetic diversity of all 24 varieties revealed similarity coefficients that ranged from 0.51 to 0.80. The EST-SSRs developed in this study represent the first large-scale development of SSR markers for ramie. These SSR markers could be used for development of genetic and physical maps, quantitative trait loci mapping, genetic diversity studies, association mapping, and cultivar fingerprinting.

Identification of single-copy orthologous genes between Physalis and Solanum lycopersicum and analysis of genetic diversity in Physalis using molecular markers

The genus Physalis includes a number of commercially important edible and ornamental species. Its high nutritional value and potential medicinal properties leads to the increased commercial interest in the products of this genus worldwide. However, lack of molecular markers prevents the detailed study of genetics and phylogeny in Physalis, which limits the progress of breeding. In the present study, we compared the DNA sequences between Physalis and tomato, and attempted to analyze genetic diversity in Physalis using tomato markers. Blasting 23180 DNA sequences derived from Physalis against the International Tomato Annotation Group (ITAG) Release2.3 Predicted CDS (SL2.40) discovered 3356 single-copy orthologous genes between them. A total of 38 accessions from at least six species of Physalis were subjected to genetic diversity analysis using 97 tomato markers and 25 SSR markers derived from P. peruviana. Majority (73.2%) of tomato markers could amplify DNA fragments from at least one accession of Physalis. Diversity in Physalis at molecular level was also detected. The average Nei's genetic distance between accessions was 0.3806 with a range of 0.2865 to 0.7091. These results indicated Physalis and tomato had similarity at both molecular marker and DNA sequence levels. Therefore, the molecular markers developed in tomato can be used in genetic study in Physalis.

The Physalis peruviana leaf transcriptome: assembly, annotation and gene model prediction

Background

Physalis peruviana commonly known as Cape gooseberry is a member of the Solanaceae family that has an increasing popularity due to its nutritional and medicinal values. A broad range of genomic tools is available for other Solanaceae, including tomato and potato. However, limited genomic resources are currently available for Cape gooseberry.

Results

We report the generation of a total of 652,614 P. peruviana Expressed Sequence Tags (ESTs), using 454 GS FLX Titanium technology. ESTs, with an average length of 371 bp, were obtained from a normalized leaf cDNA library prepared using a Colombian commercial variety. De novo assembling was performed to generate a collection of 24,014 isotigs and 110,921 singletons, with an average length of 1,638 bp and 354 bp, respectively. Functional annotation was performed using NCBI’s BLAST tools and Blast2GO, which identified putative functions for 21,191 assembled sequences, including gene families involved in all the major biological processes and molecular functions as well as defense response and amino acid metabolism pathways. Gene model predictions in P. peruviana were obtained by using the genomes of Solanum lycopersicum (tomato) and Solanum tuberosum (potato). We predict 9,436 P. peruviana sequences with multiple-exon models and conserved intron positions with respect to the potato and tomato genomes. Additionally, to study species diversity we developed 5,971 SSR markers from assembled ESTs.

Conclusions

We present the first comprehensive analysis of the Physalis peruviana leaf transcriptome, which will provide valuable resources for development of genetic tools in the species. Assembled transcripts with gene models could serve as potential candidates for marker discovery with a variety of applications including: functional diversity, conservation and improvement to increase productivity and fruit quality. P. peruviana was estimated to be phylogenetically branched out before the divergence of five other Solanaceae family members, S. lycopersicum, S. tuberosum, Capsicum spp, S. melongena and Petunia spp.