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

Advances in Physalis molecular research: applications in authentication, genetic diversity, phylogenetics, functional genes, and omics

The plants of the genus Physalis L. have been extensively utilized in traditional and indigenous Chinese medicinal practices for treating a variety of ailments, including dermatitis, malaria, asthma, hepatitis, and liver disorders. The present review aims to achieve a comprehensive and up-to-date investigation of the genus Physalis, a new model crop, to understand plant diversity and fruit development. Several chloroplast DNA-, nuclear ribosomal DNA-, and genomic DNA-based markers, such as psbA-trnH, internal-transcribed spacer (ITS), simple sequence repeat (SSR), random amplified microsatellites (RAMS), sequence-characterized amplified region (SCAR), and single nucleotide polymorphism (SNP), were developed for molecular identification, genetic diversity, and phylogenetic studies of Physalis species. A large number of functional genes involved in inflated calyx syndrome development (AP2-L, MPF2, MPF3, and MAGO), organ growth (AG1, AG2, POS1, and CNR1), and active ingredient metabolism (24ISO, DHCRT, P450-CPL, SR, DUF538, TAS14, and 3β-HSB) were identified contributing to the breeding of novel Physalis varieties. Various omic studies revealed and functionally identified a series of reproductive organ development-related factors, environmental stress-responsive genes, and active component biosynthesis-related enzymes. The chromosome-level genomes of Physalis floridana Rydb., Physalis grisea (Waterf.) M. Martínez, and Physalis pruinosa L. have been recently published providing a valuable resource for genome editing in Physalis crops. Our review summarizes the recent progress in genetic diversity, molecular identification, phylogenetics, functional genes, and the application of omics in the genus Physalis and accelerates efficient utilization of this traditional herb.

Comparative transcriptomic analysis reveals key components controlling spathe color in Anthurium andraeanum (Hort.)

Anthurium andraeanum (Hort.) is an important ornamental in the tropical cut-flower industry. However, there is currently insufficient information to establish a clear connection between the genetic model(s) proposed and the putative genes involved in the differentiation between colors. In this study, 18 cDNA libraries related to the spathe color and developmental stages of A. andraeanum were characterized by transcriptome sequencing (RNA-seq). For the de novo transcriptome, a total of 114,334,082 primary sequence reads were obtained from the Illumina sequencer and were assembled into 151,652 unigenes. Approximately 58,476 transcripts were generated and used for comparative transcriptome analysis between three cultivars that differ in spathe color (‘Sasha’ (white), ‘Honduras’ (red), and ‘Rapido’ (purple)). A large number of differentially expressed genes (8,324), potentially involved in multiple biological and metabolic pathways, were identified, including genes in the flavonoid and anthocyanin biosynthetic pathways. Our results showed that the chalcone isomerase (CHI) gene presented the strongest evidence for an association with differences in color and the highest correlation with other key genes (flavanone 3-hydroxylase (F3H), flavonoid 3’5’ hydroxylase (F3’5’H)/ flavonoid 3’-hydroxylase (F3’H), and leucoanthocyanidin dioxygenase (LDOX)) in the anthocyanin pathway. We also identified a differentially expressed cytochrome P450 gene in the late developmental stage of the purple spathe that appeared to determine the difference between the red- and purple-colored spathes. Furthermore, transcription factors related to putative MYB-domain protein that may control anthocyanin pathway were identified through a weighted gene co-expression network analysis (WGCNA). The results provided basic sequence information for future research on spathe color, which have important implications for this ornamental breeding strategies.

Integrating Omics and Gene Editing Tools for Rapid Improvement of Traditional Food Plants for Diversified and Sustainable Food Security

Indigenous communities across the globe, especially in rural areas, consume locally available plants known as Traditional Food Plants (TFPs) for their nutritional and health-related needs. Recent research shows that many TFPs are highly nutritious as they contain health beneficial metabolites, vitamins, mineral elements and other nutrients. Excessive reliance on the mainstream staple crops has its own disadvantages. Traditional food plants are nowadays considered important crops of the future and can act as supplementary foods for the burgeoning global population. They can also act as emergency foods in situations such as COVID-19 and in times of other pandemics. The current situation necessitates locally available alternative nutritious TFPs for sustainable food production. To increase the cultivation or improve the traits in TFPs, it is essential to understand the molecular basis of the genes that regulate some important traits such as nutritional components and resilience to biotic and abiotic stresses. The integrated use of modern omics and gene editing technologies provide great opportunities to better understand the genetic and molecular basis of superior nutrient content, climate-resilient traits and adaptation to local agroclimatic zones. Recently, realizing the importance and benefits of TFPs, scientists have shown interest in the prospection and sequencing of TFPs for their improvements, cultivation and mainstreaming. Integrated omics such as genomics, transcriptomics, proteomics, metabolomics and ionomics are successfully used in plants and have provided a comprehensive understanding of gene-protein-metabolite networks. Combined use of omics and editing tools has led to successful editing of beneficial traits in several TFPs. This suggests that there is ample scope for improvement of TFPs for sustainable food production. In this article, we highlight the importance, scope and progress towards improvement of TFPs for valuable traits by integrated use of omics and gene editing techniques.

Combining transcriptome analysis and GWAS for identification and validation of marker genes in the Physalis peruviana-Fusarium oxysporum pathosystem

Vascular wilt, caused by the pathogen Fusarium oxysporum f. sp. physali (Foph), is a major disease of cape gooseberry (Physalis peruviana L.) in Andean countries. Despite the economic losses caused by this disease, there are few studies related to molecular mechanisms in the P. peruviana—Foph pathosystem as a useful tool for crop improvement. This study evaluates eight candidate genes associated with this pathosystem, using real-time quantitative PCR (RT-qPCR). The genes were identified and selected from 1,653 differentially expressed genes (DEGs) derived from RNA-Seq analysis and from a previous genome-wide association study (GWAS) of this plant-pathogen interaction. Based on the RT-qPCR analysis, the tubuline (TUB) reference gene was selected for its highly stable expression in cape gooseberry. The RT-qPCR validation of the candidate genes revealed the biological variation in their expression according to their known biological function. Three genes related to the first line of resistance/defense responses were highly expressed earlier during infection in a susceptible genotype, while three others were overexpressed later, mostly in the tolerant genotype. These genes are mainly involved in signaling pathways after pathogen recognition, mediated by hormones such as ethylene and salicylic acid. This study provided the first insight to uncover the molecular mechanism from the P. peruviana—Foph pathosystem. The genes validated here have important implications in the disease progress and allow a better understanding of the defense response in cape gooseberry at the molecular level. Derived molecular markers from these genes could facilitate the identification of tolerant/susceptible genotypes for use in breeding schemes.

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.

Transcriptomic variation of the flower-fruit transition in Physalis and Solanum

Gene expression variations in response to fertilization between Physalis and Solanum might play essential roles in species divergence and fruit evolution. Fertilization triggers variation in fruit development and morphology. The Chinese lantern, a morphological novelty derived from the calyx, is formed upon fertilization in Physalis but is not observed in Solanum. The underlying genetic variations are largely unknown. Here, we documented the developmental and morphological differences in the flower and fruit between Physalis floridana and Solanum pimpinellifolium and then evaluated both the transcript sequence variation and gene expression at the transcriptomic level at fertilization between the two species. In Physalis transcriptomic analysis, 468 unigenes were identified as differentially expressed genes (DEGs) that were strongly regulated by fertilization across 3 years. In comparison with tomato, 14,536 strict single-copy orthologous gene pairs were identified between P. floridana and S. pimpinellifolium in the flower-fruit transcriptome. Nine types of gene variations with specific GO-enriched patterns were identified, covering 58.82% orthologous gene pairs that were DEGs in either trend or dosage at the flower-fruit transition between the two species, which could adequately distinguish Solanum and Physalis, implying that differential gene expression at fertilization might play essential roles during the divergence and fruit evolution of Solanum-Physalis. Virus-induced gene silencing analyses revealed the developmental roles of some transcription factor genes in fertility, Chinese lantern development, and fruit weight control in Physalis. This study presents the first floral transcriptomic resource of Physalis, and reveals some candidate genetic variations accounting for the early fruit developmental evolution in Physalis in comparison to Solanum.

Bioactive compounds induced in Physalis angulata L. by methyl-jasmonate: an investigation of compound accumulation patterns and biosynthesis-related candidate genes

We employed both metabolomic and transcriptomic approaches to explore the accumulation patterns of physalins, flavonoids and chlorogenic acid in Physalis angulata and revealed the genes associated with the biosynthesis of bioactive compounds under methyl-jasmonate (MeJA) treatment. Physalis angulata L. is an annual Solanaceae plant with a number of medicinally active compounds. Despite the potential pharmacological benefits of P. angulata, the scarce genomic information regarding this plant has limited the studies on the mechanisms of bioactive compound biosynthesis. To facilitate the basic understanding of the main chemical constituent biosynthesis pathways, we performed both metabolomic and transcriptomic approaches to reveal the genes associated with the biosynthesis of bioactive compounds under methyl-jasmonate (MeJA) treatment. Untargeted metabolome analysis showed that most physalins, flavonoids and chlorogenic acid were significantly upregulated. Targeted HPLC-MS/MS analysis confirmed variations in the contents of two important representative steroid derivatives (physalins B and G), total flavonoids, neochlorogenic acid, and chlorogenic acid between MeJA-treated plants and controls. Transcript levels of a few steroid biosynthesis-, flavonoid biosynthesis-, and chlorogenic acid biosynthesis-related genes were upregulated, providing a potential explanation for MeJA-induced active ingredient synthesis in P. angulata. Systematic correlation analysis identified a number of novel candidate genes associated with bioactive compound biosynthesis. These results may help to elucidate the regulatory mechanism underlying MeJA-induced active compound accumulation and provide several valuable candidate genes for further functional study.

Endogenous Hormones Inhibit Differentiation of Young Ears in Maize (Zea mays L.) Under Heat Stress

Global warming frequently leads to extreme temperatures, which pose a serious threat to the growth, development, and yield formation of crops such as maize. This study aimed to deeply explore the molecular mechanisms of young ear development under heat stress. We selected the heat-tolerant maize variety Zhengdan 958 (T) and heat-sensitive maize variety Xianyu 335 (S), and subjected them to heat stress in the V9 (9th leaf), V12 (12th leaf), and VT (tasseling) growth stages. We combined analysis of the maize phenotype with omics technology and physiological indicators to compare the differences in young ear morphology, total number of florets, floret fertilization rate, grain abortion rate, number of grains, and main metabolic pathways between plants subjected to heat stress and those left to develop normally. The results showed that after heat stress, the length and diameter of young ears, total number of florets, floret fertilization rate, and number of grains all decreased significantly, whereas the length of the undeveloped part at the top of the ear and grain abortion rate increased significantly. In addition, the differentially expressed genes (DEGs) in young ears were significantly enriched in the hormone signaling pathways. The endogenous hormone content in young ears exhibited different changes: zeatin (ZT) and zeatin riboside (ZR) decreased significantly, but gibberellin acid₃ (GA₃), gibberellin acid₄ (GA₄), and abscisic acid (ABA) increased significantly, in ears subjected to heat stress. In the heat-tolerant maize variety, the salicylic acid (SA), and jasmonic acid (JA) content in the vegetative growth stage also increased in ears subjected to heat stress, whereas the opposite effect was observed for the heat-sensitive variety. The changes in endogenous hormone content of young ears that were subjected to heat stress significantly affected ear development, resulting in a reduction in the number of differentiated florets, fertilized florets and grains, which ultimately reduced the maize yield.

Genome-wide transcriptome profiling of transgenic hop (Humulus lupulus L.) constitutively overexpressing HlWRKY1 and HlWDR1 transcription factors

Background

The hop plant (Humulus lupulus L.) is a valuable source of several secondary metabolites, such as flavonoids, bitter acids, and essential oils. These compounds are widely implicated in the beer brewing industry and are having potential biomedical applications. Several independent breeding programs around the world have been initiated to develop new cultivars with enriched lupulin and secondary metabolite contents but met with limited success due to several constraints. In the present work, a pioneering attempt has been made to overexpress master regulator binary transcription factor complex formed by HlWRKY1 and HlWDR1 using a plant expression vector to enhance the level of prenylflavonoid and bitter acid content in the hop. Subsequently, we performed transcriptional profiling using high-throughput RNA-Seq technology in leaves of resultant transformants and wild-type hop to gain in-depth information about the genome-wide functional changes induced by HlWRKY1 and HlWDR1 overexpression.

Results

The transgenic WW-lines exhibited an elevated expression of structural and regulatory genes involved in prenylflavonoid and bitter acid biosynthesis pathways. In addition, the comparative transcriptome analysis revealed a total of 522 transcripts involved in 30 pathways, including lipids and amino acids biosynthesis, primary carbon metabolism, phytohormone signaling and stress responses were differentially expressed in WW-transformants. It was apparent from the whole transcriptome sequencing that modulation of primary carbon metabolism and other pathways by HlWRKY1 and HlWDR1 overexpression resulted in enhanced substrate flux towards secondary metabolites pathway. The detailed analyses suggested that none of the pathways or genes, which have a detrimental effect on physiology, growth and development processes, were induced on a genome-wide scale in WW-transgenic lines.

Conclusions

Taken together, our results suggest that HlWRKY1 and HlWDR1 simultaneous overexpression positively regulates the prenylflavonoid and bitter acid biosynthesis pathways in the hop and thus these transgenes are presented as prospective candidates for achieving enhanced secondary metabolite content in the hop.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5125-8) contains supplementary material, which is available to authorized users.

Rapid improvement of domestication traits in an orphan crop by genome editing

Genome editing holds great promise for increasing crop productivity, and there is particular interest in advancing breeding in orphan crops, which are often burdened by undesirable characteristics resembling wild relatives. We developed genomic resources and efficient transformation in the orphan Solanaceae crop 'groundcherry' (Physalis pruinosa) and used clustered regularly interspaced short palindromic repeats (CRISPR)-CRISPR-associated protein-9 nuclease (Cas9) (CRISPR-Cas9) to mutate orthologues of tomato domestication and improvement genes that control plant architecture, flower production and fruit size, thereby improving these major productivity traits. Thus, translating knowledge from model crops enables rapid creation of targeted allelic diversity and novel breeding germplasm in distantly related orphan crops.

Comparative Metabolomic and Proteomic Analyses Reveal the Regulation Mechanism Underlying MeJA-Induced Bioactive Compound Accumulation in Cutleaf Groundcherry ( Physalis angulata L.) Hairy Roots

Cutleaf groundcherry ( Physalis angulata L.) is an annual plant with a number of medicinal ingredients. However, studies about the secondary metabolism of P. angulata are very limited. An integrated metabolome and proteome approach was used to reveal the variations in the metabolism associated with bioactive compounds under methyl-jasmonate (MeJA) treatment. Application of MeJA to the hairy roots could significantly increase the accumulation of most active ingredients. A targeted approach confirmed the variations in physalins D and H between MeJA treatment and the controls. Increases in the levels of a number of terpenoid backbone biosynthesis and steroid biosynthesis related enzymes, cytochrome P450 monooxygenases and 3β-hydroxysterioid dehydrogenase might provide a potential explanation for the MeJA-induced active ingredient synthesis. Our results may contribute to a deeper understanding of the regulation mechanism underlying the MeJA-induced active compound accumulation in P. angulata.

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.

Association Study Reveals Novel Genes Related to Yield and Quality of Fruit in Cape Gooseberry (Physalis peruviana L.)

Association mapping has been proposed as an efficient approach to assist plant breeding programs to investigate the genetic basis of agronomic traits. In this study, we evaluated 18 traits related to yield, (FWP, NF, FWI, and FWII), fruit size-shape (FP, FA, MW, WMH, MH, HMW, DI, FSI, FSII, OVO, OBO), and fruit quality (FIR, CF, and SST), in a diverse collection of 100 accessions of Physalis peruviana including wild, landraces, and anther culture derived lines. We identified seven accessions with suitable traits: fruit weight per plant (FWP) > 7,000 g/plant and cracked fruits (CF) < 4%, to be used as parents in cape gooseberry breeding program. In addition, the accessions were also characterized using Genotyping By Sequencing (GBS). We discovered 27,982 and 36,142 informative SNP markers based on the alignment against the two cape gooseberry references transcriptomes. Besides, 30,344 SNPs were identified based on alignment to the tomato reference genome. Genetic structure analysis showed that the population could be divided into two or three sub-groups, corresponding to landraces-anther culture and wild accessions for K = 2 and wild, landraces, and anther culture plants for K = 3. Association analysis was carried out using a Mixed Linear Model (MLM) and 34 SNP markers were significantly associated. These results reveal the basis of the genetic control of important agronomic traits and may facilitate marker-based breeding in P. peruviana.

Eocene lantern fruits from Gondwanan Patagonia and the early origins of Solanaceae

The nightshade family Solanaceae holds exceptional economic and cultural importance. The early diversification of Solanaceae is thought to have occurred in South America during its separation from Gondwana, but the family's sparse fossil record provides few insights. We report 52.2-million-year-old lantern fruits from terminal-Gondwanan Patagonia, featuring highly inflated, five-lobed calyces, as a newly identified species of the derived, diverse New World genus Physalis (e.g., groundcherries and tomatillos). The fossils are considerably older than corresponding molecular divergence dates and demonstrate an ancient history for the inflated calyx syndrome. The derived position of these early Eocene fossils shows that Solanaceae were well diversified long before final Gondwanan breakup.

Analysis of gene expression for microminipig liver transcriptomes using parallel long-read technology and short-read sequencing

The microminipig is one of the smallest minipigs that has emerged as a possible experimental animal model, because it shares many anatomical and/or physiological similarities with humans, including the coronary artery distribution in the heart, the digestive physiology, the kidney size and its structure, and so on. However, information on gene expression profiles, including those on drug-metabolizing phase I and II enzymes, in the microminipig is limited. Therefore, the aim of the present study was to identify transcripts in microminipig livers and to determine gene expression profiles. De novo assembly and expression analyses of microminipig transcripts were conducted with liver samples from three male and three female microminipigs using parallel long-read and short-read sequencing technologies. After unique sequences had been automatically aligned by assembling software, the mean contig length of 50843 transcripts was 707 bp. The expression profiles of cytochrome P450 (P450) 1A2, 2C, 2E1 and 3A genes in livers in microminipigs were similar to those in humans. Liver carboxylesterase (CES) precursor, liver CES-like, UDP-glucuronosyltransferase (UGT) 2C1-like, amine sulfotransferase (SULT)-like, N-acetyltransferases (NAT8) and glutathione S-transferase (GST) A2 genes, which are relatively unknown genes in pigs and/or humans, were expressed strongly. Furthermore, no significant gender differences were observed in the gene expression profiles of phase I enzymes, whereas UGT2B17, SULT1E1, SULT2A1, amine SULT-like, NAT8 and GSTT4 genes were different between males and females among phase II enzyme genes under the present sample conditions. These results provide a foundation for mechanistic studies and the use of microminipigs as model animals for drug development in the future. Copyright © 2016 John Wiley & Sons, Ltd.

Workflow and web application for annotating NCBI BioProject transcriptome data

Abstract

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

Combination of transcriptomic and proteomic approaches helps to unravel the protein composition of Chelidonium majus L. milky sap

A novel annotated Chelidonium majus L. transcriptome database composed of 23,004 unique coding sequences allowed to significantly improve the sensitivity of proteomic C. majus assessments, which showed novel defense-related proteins characteristic to its latex. To date, the composition of Chelidonium majus L. milky sap and biosynthesis of its components are poorly characterized. We, therefore, performed de novo sequencing and assembly of C. majus transcriptome using Illumina technology. Approximately, 119 Mb of raw sequence data was obtained. Assembly resulted in 107,088 contigs, with N50 of 1913 bp and N90 of 450 bp. Among 34,965 unique coding sequences (CDS), 23,004 obtained CDS database served as a basis for further proteomic analyses. The database was then used for the identification of proteins from C. majus milky sap, and whole plant extracts analyzed using liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) approach. Of about 334 different putative proteins were identified in C. majus milky sap and 1155 in C. majus whole plant extract. The quantitative comparative analysis confirmed that C. majus latex contains proteins connected with response to stress conditions and generation of precursor metabolites and energy. Notable proteins characteristic to latex include major latex protein (MLP, presumably belonging to Bet v1-like superfamily), polyphenol oxidase (PPO, which could be responsible for browning of the sap after exposure to air), and enzymes responsible for anthocyanidin, phenylpropanoid, and alkaloid biosynthesis.

De Novo Assembly and Comparative Transcriptome Analysis Provide Insight into Lysine Biosynthesis in Toona sinensis Roem

Toona sinensis Roem is a popular leafy vegetable in Chinese cuisine and is also used as a traditional Chinese medicine. In this study, leaf samples were collected from the same plant on two development stages and then used for high-throughput Illumina RNA-sequencing (RNA-Seq). 125,884 transcripts and 54,628 unigenes were obtained through de novo assembly. A total of 25,570 could be annotated with known biological functions, which indicated that the T. sinensis leaves and shoots were undergoing multiple developmental processes especially for active metabolic processes. Analysis of differentially expressed unigenes between the two libraries showed that the lysine biosynthesis was an enriched KEGG pathway, and candidate genes involved in the lysine biosynthesis pathway in T. sinensis leaves and shoots were identified. Our results provide a primary analysis of the gene expression files of T. sinensis leaf and shoot on different development stages and afford a valuable resource for genetic and genomic research on plant lysine biosynthesis.

Genome Sequence of the Banana Plant Growth-Promoting Rhizobacterium Pseudomonas fluorescens PS006

Pseudomonas fluorescens is a well-known plant growth-promoting rhizobacterium (PGPR). We report here the first whole-genome sequence of PGPR P. fluorescens evaluated in Colombian banana plants. The genome sequences contains genes involved in plant growth and defense, including bacteriocins, 1-aminocyclopropane-1-carboxylic acid (ACC) deaminase, and genes that provide resistance to toxic compounds.

The first de novo transcriptome of pepino (Solanum muricatum): assembly, comprehensive analysis and comparison with the closely related species S. caripense, potato and tomato

BACKGROUND

Solanum sect. Basarthrum is phylogenetically very close to potatoes (Solanum sect. Petota) and tomatoes (Solanum sect. Lycopersicon), two groups with great economic importance, and for which Solanum sect. Basarthrum represents a tertiary gene pool for breeding. This section includes the important regional cultigen, the pepino (Solanum muricatum), and several wild species. Among the wild species, S. caripense is prominent due to its major involvement in the origin of pepino and its wide geographical distribution. Despite the value of the pepino as an emerging crop, and the potential for gene transfer from both the pepino and S. caripense to potatoes and tomatoes, there has been virtually no genomic study of these species.

RESULTS

Using Illumina HiSeq 2000, RNA-Seq was performed with a pool of three tissues (young leaf, flowers in pre-anthesis and mature fruits) from S. muricatum and S. caripense, generating almost 111,000,000 reads among the two species. A high quality de novo transcriptome was assembled from S. muricatum clean reads resulting in 75,832 unigenes with an average length of 704 bp. These unigenes were functionally annotated based on similarity of public databases. We used Blast2GO, to conduct an exhaustive study of the gene ontology, including GO terms, EC numbers and KEGG pathways. Pepino unigenes were compared to both potato and tomato genomes in order to determine their estimated relative position, and to infer gene prediction models. Candidate genes related to traits of interest in other Solanaceae were evaluated by presence or absence and compared with S. caripense transcripts. In addition, by studying five genes, the phylogeny of pepino and five other members of the family, Solanaceae, were studied. The comparison of S. caripense reads against S. muricatum assembled transcripts resulted in thousands of intra- and interspecific nucleotide-level variants. In addition, more than 1000 SSRs were identified in the pepino transcriptome.

CONCLUSIONS

This study represents the first genomic resource for the pepino. We suggest that the data will be useful not only for improvement of the pepino, but also for potato and tomato breeding and gene transfer. The high quality of the transcriptome presented here also facilitates comparative studies in the genus Solanum. The accurate transcript annotation will enable us to figure out the gene function of particular traits of interest. The high number of markers (SSR and nucleotide-level variants) obtained will be useful for breeding programs, as well as studies of synteny, diversity evolution, and phylogeny.

Comparative Characterization of the Leaf Tissue of Physalis alkekengi and Physalis peruviana Using RNA-seq and Metabolite Profiling

The genus Physalis in the Solanaceae family contains several species of benefit to humans. Examples include P. alkekengi (Chinese-lantern plant, hôzuki in Japanese) used for medicinal and for decorative purposes, and P. peruviana, also known as Cape gooseberry, which bears an edible, vitamin-rich fruit. Members of the Physalis genus are a valuable resource for phytochemicals needed for the development of medicines and functional foods. To fully utilize the potential of these phytochemicals we need to understand their biosynthesis, and for this we need genomic data, especially comprehensive transcriptome datasets for gene discovery. We report the de novo assembly of the transcriptome from leaves of P. alkekengi and P. peruviana using Illumina RNA-seq technologies. We identified 75,221 unigenes in P. alkekengi and 54,513 in P. peruviana. All unigenes were annotated with gene ontology (GO), Enzyme Commission (EC) numbers, and pathway information from the Kyoto Encyclopedia of Genes and Genomes (KEGG). We classified unigenes encoding enzyme candidates putatively involved in the secondary metabolism and identified more than one unigenes for each step in terpenoid backbone- and steroid biosynthesis in P. alkekengi and P. peruviana. To measure the variability of the withanolides including physalins and provide insights into their chemical diversity in Physalis, we also analyzed the metabolite content in leaves of P. alkekengi and P. peruviana at five different developmental stages by liquid chromatography-mass spectrometry. We discuss that comprehensive transcriptome approaches within a family can yield a clue for gene discovery in Physalis and provide insights into their complex chemical diversity. The transcriptome information we submit here will serve as an important public resource for further studies of the specialized metabolism of Physalis species.

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.

Genome Sequence of the Banana Plant Growth-Promoting Rhizobacterium Bacillus amyloliquefaciens BS006

Bacillus amyloliquefaciens is an important plant growth-promoting rhizobacterium (PGPR). We report the first whole-genome sequence of PGPR Bacillus amyloliquefaciens evaluated in Colombian banana plants. The genome sequences encode genes involved in plant growth and defense, including bacteriocins, ribosomally synthesized antibacterial peptides, in addition to genes that provide resistance to toxic compounds.

Evolutionarily Dynamic, but Robust, Targeting of Resistance Genes by the miR482/2118 Gene Family in the Solanaceae

Plants are exposed to pathogens around the clock. A common resistance response in plants upon pathogen detection is localized cell death. Given the irreversible nature of this response, multiple layers of negative regulation are present to prevent the untimely or misexpression of resistance genes. One layer of negative regulation is provided by a recently discovered microRNA (miRNA) gene family, miR482/2118. This family targets the transcripts of resistance genes in plants. We investigated the evolutionary history and specificity of this miRNA gene family within the Solanaceae. This plant family includes many important crop species, providing a set of well-defined resistance gene repertoires. Across 14 species from the Solanaceae, we identified eight distinct miR482/2118 gene family members. Our studies show conservation of miRNA type and number in the group of wild tomatoes and, to a lesser extent, throughout the Solanaceae. The eight orthologous miRNA gene clusters evolved under different evolutionary constraints, allowing for individual subfunctionalization of the miRNAs. Despite differences in the predicted targeting behavior of each miRNA, the miRNA-R-gene network is robust due to its high degree of interconnectivity and redundant targeting. Our data suggest that the miR482/2118 gene family acts as an evolutionary buffer for R-gene sequence diversity.

Next Generation Sequencing and Transcriptome Analysis Predicts Biosynthetic Pathway of Sennosides from Senna (Cassia angustifolia Vahl.), a Non-Model Plant with Potent Laxative Properties

Senna (Cassia angustifolia Vahl.) is a world's natural laxative medicinal plant. Laxative properties are due to sennosides (anthraquinone glycosides) natural products. However, little genetic information is available for this species, especially concerning the biosynthetic pathways of sennosides. We present here the transcriptome sequencing of young and mature leaf tissue of Cassia angustifolia using Illumina MiSeq platform that resulted in a total of 6.34 Gb of raw nucleotide sequence. The sequence assembly resulted in 42230 and 37174 transcripts with an average length of 1119 bp and 1467 bp for young and mature leaf, respectively. The transcripts were annotated using NCBI BLAST with 'green plant database (txid 33090)', Swiss Prot, Kyoto Encylcopedia of Genes & Genomes (KEGG), Cluster of Orthologous Gene (COG) and Gene Ontology (GO). Out of the total transcripts, 40138 (95.0%) and 36349 (97.7%) from young and mature leaf, respectively, were annotated by BLASTX against green plant database of NCBI. We used InterProscan to see protein similarity at domain level, a total of 34031 (young leaf) and 32077 (mature leaf) transcripts were annotated against the Pfam domains. All transcripts from young and mature leaf were assigned to 191 KEGG pathways. There were 166 and 159 CDS, respectively, from young and mature leaf involved in metabolism of terpenoids and polyketides. Many CDS encoding enzymes leading to biosynthesis of sennosides were identified. A total of 10,763 CDS differentially expressing in both young and mature leaf libraries of which 2,343 (21.7%) CDS were up-regulated in young compared to mature leaf. Several differentially expressed genes found functionally associated with sennoside biosynthesis. CDS encoding for many CYPs and TF families were identified having probable roles in metabolism of primary as well as secondary metabolites. We developed SSR markers for molecular breeding of senna. We have identified a set of putative genes involved in various secondary metabolite pathways, especially those related to the synthesis of sennosides which will serve as an important platform for public information about gene expression, genomics, and functional genomics in senna.

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.

Qualitative de novo analysis of full length cDNA and quantitative analysis of gene expression for common marmoset (Callithrix jacchus) transcriptomes using parallel long-read technology and short-read sequencing

The common marmoset (Callithrix jacchus) is a non-human primate that could prove useful as human pharmacokinetic and biomedical research models. The cytochromes P450 (P450s) are a superfamily of enzymes that have critical roles in drug metabolism and disposition via monooxygenation of a broad range of xenobiotics; however, information on some marmoset P450s is currently limited. Therefore, identification and quantitative analysis of tissue-specific mRNA transcripts, including those of P450s and flavin-containing monooxygenases (FMO, another monooxygenase family), need to be carried out in detail before the marmoset can be used as an animal model in drug development. De novo assembly and expression analysis of marmoset transcripts were conducted with pooled liver, intestine, kidney, and brain samples from three male and three female marmosets. After unique sequences were automatically aligned by assembling software, the mean contig length was 718 bp (with a standard deviation of 457 bp) among a total of 47,883 transcripts. Approximately 30% of the total transcripts were matched to known marmoset sequences. Gene expression in 18 marmoset P450- and 4 FMO-like genes displayed some tissue-specific patterns. Of these, the three most highly expressed in marmoset liver were P450 2D-, 2E-, and 3A-like genes. In extrahepatic tissues, including brain, gene expressions of these monooxygenases were lower than those in liver, although P450 3A4 (previously P450 3A21) in intestine and P450 4A11- and FMO1-like genes in kidney were relatively highly expressed. By means of massive parallel long-read sequencing and short-read technology applied to marmoset liver, intestine, kidney, and brain, the combined next-generation sequencing analyses reported here were able to identify novel marmoset drug-metabolizing P450 transcripts that have until now been little reported. These results provide a foundation for mechanistic studies and pave the way for the use of marmosets as model animals for drug development in the future.

De novo assembly and characterization of leaf transcriptome for the development of functional molecular markers of the extremophile multipurpose tree species Prosopis alba

Background

Prosopis alba (Fabaceae) is an important native tree adapted to arid and semiarid regions of north-western Argentina which is of great value as multipurpose species. Despite its importance, the genomic resources currently available for the entire Prosopis genus are still limited. Here we describe the development of a leaf transcriptome and the identification of new molecular markers that could support functional genetic studies in natural and domesticated populations of this genus.

Results

Next generation DNA pyrosequencing technology applied to P. alba transcripts produced a total of 1,103,231 raw reads with an average length of 421 bp. De novo assembling generated a set of 15,814 isotigs and 71,101 non-assembled sequences (singletons) with an average of 991 bp and 288 bp respectively. A total of 39,000 unique singletons were identified after clustering natural and artificial duplicates from pyrosequencing reads.

Regarding the non-redundant sequences or unigenes, 22,095 out of 54,814 were successfully annotated with Gene Ontology terms. Moreover, simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were searched, resulting in 5,992 and 6,236 markers, respectively, throughout the genome. For the validation of the the predicted SSR markers, a subset of 87 SSRs selected through functional annotation evidence was successfully amplified from six DNA samples of seedlings. From this analysis, 11 of these 87 SSRs were identified as polymorphic. Additionally, another set of 123 nuclear polymorphic SSRs were determined in silico, of which 50% have the probability of being effectively polymorphic.

Conclusions

This study generated a successful global analysis of the P. alba leaf transcriptome after bioinformatic and wet laboratory validations of RNA-Seq data.

The limited set of molecular markers currently available will be significantly increased with the thousands of new markers that were identified in this study. This information will strongly contribute to genomics resources for P. alba functional analysis and genetics. Finally, it will also potentially contribute to the development of population-based genome studies in the genera.

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.

Bioinformatics challenges in de novo transcriptome assembly using short read sequences in the absence of a reference genome sequence

Plant natural product research can be facilitated through genome and transcriptome sequencing approaches that generate informative sequence and expression datasets that enable characterization of biochemical pathways of interest. As the overwhelming majority of plant-derived natural products are derived from species with little, if any, sequence and/or genomic resources, the ability to perform whole genome shotgun sequencing and assembly has been and will continue to be transformative as access to a genome sequence provides molecular resources and a context for discovery and characterization of biosynthetic pathways. Due to the reduced size and complexity of the transcriptome relative to the genome, transcriptome sequencing provides a rapid, inexpensive approach to access gene sequences, gene expression abundances, and gene expression patterns in any species, including those that lack a reference genome sequence. To date, successful applications of RNA sequencing in conjunction with de novo transcriptome assembly has enabled identification of new genes in an array of biochemical pathways in plants. While sequencing technologies are well developed, challenges remain in the handling and analysis of transcriptome sequences. In this Highlight article, we provide an overview of the bioinformatics challenges associated with transcriptome analyses using short read sequences and how to address these issues in plant species that lack a reference genome.

Transcriptome characterisation of Pinus tabuliformis and evolution of genes in the Pinus phylogeny

BACKGROUND

The Chinese pine (Pinus tabuliformis) is an indigenous conifer species in northern China but is relatively underdeveloped as a genomic resource; thus, limiting gene discovery and breeding. Large-scale transcriptome data were obtained using a next-generation sequencing platform to compensate for the lack of P. tabuliformis genomic information.

RESULTS

The increasing amount of transcriptome data on Pinus provides an excellent resource for multi-gene phylogenetic analysis and studies on how conserved genes and functions are maintained in the face of species divergence. The first P. tabuliformis transcriptome from a normalised cDNA library of multiple tissues and individuals was sequenced in a full 454 GS-FLX run, producing 911,302 sequencing reads. The high quality overlapping expressed sequence tags (ESTs) were assembled into 46,584 putative transcripts, and more than 700 SSRs and 92,000 SNPs/InDels were characterised. Comparative analysis of the transcriptome of six conifer species yielded 191 orthologues, from which we inferred a phylogenetic tree, evolutionary patterns and calculated rates of gene diversion. We also identified 938 fast evolving sequences that may be useful for identifying genes that perhaps evolved in response to positive selection and might be responsible for speciation in the Pinus lineage.

CONCLUSIONS

A large collection of high-quality ESTs was obtained, de novo assembled and characterised, which represents a dramatic expansion of the current transcript catalogues of P. tabuliformis and which will gradually be applied in breeding programs of P. tabuliformis. Furthermore, these data will facilitate future studies of the comparative genomics of P. tabuliformis and other related species.

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.